## Coherent control of Optical limiting in atomic systems

Mohsen Ghaderi Goran Abad , Mahsa Mahdieh , Mohadeseh Veisi , Hamid Nadjari , and 1 1 1 1 Mohammad Mahmoudi 1,*

1 Department of Physics, University of Zanjan, University Blvd., 45371-38791, Zanjan, Iran * mahmoudi@znu.ac.ir

## ABSTRACT

Generation and control of the reverse saturable absorption (RSA) and optical limiting (OL) are investigated in a four-level Y-type quantum system. It is demonstrated that the applied laser fields induce the RSA and it can be coherently controlled by either intensity or frequency of the applied laser fields. The effect of the static magnetic field on the induced RSA is studied and we obtain that it has a constructive role in determining the intensity range in which the OL is established in the system. In addition, we find that the OL threshold can be decreased either by increasing the length of the medium or by getting the atomic system denser. Finally, Z-scan technique is presented to confirm our theoretical results. The proposed scheme can be used in designing the coherent optical limiters with controllable threshold and intensity range of OL.

Atomic coherence offers a systematic basis for the fundamental concepts of the coherent control. Optical devices are generally based on manipulation of the phase, polarization, propagation direction and also intensity of the optical light which can be done by atomic coherence due to the applied laser fields. Intensity manipulation of the light is a major basis of designing all-optical switches 1, 2 and optical limiters 2,3 . Optical limiters have been receiving much attention for the increasing demand in protection of optical components in optically based devices. All optical sensors involved in the optical devices may be vulnerable for their sensitivity to the light intensity. In fact, human eye, sensors and other optical sensitive elements have intensity threshold above which laser-induced damage happens. In optical limiters the transmission of the light reduces or even becomes constant for the input intensities higher than the threshold intensity. Thus, the presence of optical limiters to restrict the intensity of the incident laser beam is completely requisite prior to the sensors and direct viewing devices. Optical limiters protect the sensors from the damages due to the higher intensity laser pulse by extension their intensity range to operate under rougher conditions.

Various techniques, i.e., two photon absorption 4-6 , nonlinear scattering 7, 8 have been reported for generating the OL. The basic mechanism to establish the OL is the RSA 9-12 in which, unlike the saturable absorption (SA), the absorption increases by increasing the incident intensity. In general, the RSA can occur when the absorption of the excited state is large compared to the absorption of the ground state. On the contrary, SA is the dominant phenomenon in the system. Note that in SA materials, the absorption reduces by increasing the intensity due to depletion of the ground state, leading the materials to be more transparent. The RSA and OL have been observed in numerous compounds such as organic materials 3,13 , C 60 solution 14 , nano materials 15-18 and semiconductors 19 . In addition, most studies on the OL have been carried out considering molecular systems with the aid of the rate equation approach 20-23 .

In this paper, unlike the previous reported works, we introduce a coherently controllable optical limiter using the atomic systems in a four-level Y-type configuration. It is demonstrated that the RSA is coherently induced by the laser fields and the conditions are provided for preparing the OL in the system. It is shown that all characteristics of the induced OL such as the intensity range and the threshold intensity can be controlled by either intensity or frequency of the laser fields. Moreover, the effect of the static magnetic field on the OL is studied and it is illustrated that the RSA and the OL regions are extended by increasing the magnitude of the static magnetic field. In addition, we show that the OL threshold decreases either by increasing the length of the medium or by getting the atomic medium denser. Finally, Z-scan technique is presented to confirm our theoretical results.

## Model and Equations

The proposed realistic atomic system is a four-level Y-type quantum system which can be established in 5 S 1 / 2 , 5 P 1 / 2 and 5 D 3 / 2 lines of 87 Rb atoms as shown in Fig. 1. Two lower states are | 1 〉 = | 5 S 1 / 2 , ( F = 1 , mF = ) 0 〉 and | 2 〉 = | 5 P 1 / 2 , ( F = 2 , mF = ) 0 〉 , separated by 377 THz. Two degenerate states | 3 〉 = | 5 D 3 / 2 , ( F ′ = 2 , mF ′ = -1 ) 〉 and | 4 〉 = | 5 D 3 / 2 , ( F ′ = 2 , mF ′ = ) 1 〉 are

2

chosen as the excited states. Here, F and F ′ are the quantum numbers of the total angular momentum and mF F ( ′ ) denotes magnetic quantum number of the corresponding states. A weak linearly probe field, /vector Ep ( z t , ) = /vector ε p ( ) z e -i ( ω pt -kpz ) + c c . ., with wave vector kp and polarization in ˆ direction, drives the transition x | 1 〉 ↔ | 2 〉 with the Rabi frequency Ω p = /vector µ 21 ./vector ε p / ¯. h The transition | 2 〉 ↔ | 3 〉 ( 2 | 〉 ↔ | 4 ) is coupled by the strong left (right) circularly polarized coupling field with the Rabi frequency 〉 Ω s = /vector µ 32 . ˆ ε -Es / ¯ h ( Ω c = /vector µ 42 . ˆ ε + Ec / ¯). h ε p and Ei ( i = s c , ) are the amplitude of the probe and coupling fields, respectively. ε ± stands for the polarization unit vectors of the coupling fields. A static magnetic field is also employed to remove the degeneracy of the states | 3 〉 and | 4 〉 by 2¯ h ∆ B = 2 mFgs µ BB where µ B is Bohr magneton and gs is Land´ factor. e ☎

F

,

Figure 1. A schematic of a four-level Y-type quantum system driven by a weak probe field and two coupling fields with Rabi frequencies Ω p , Ω s and Ω c respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01186v1-figure-7524a4d4-aede-4ec2-bcee-9458111bf78b.jpg)

The Hamiltonian of the considered system in the dipole and rotating wave approximations can be written as

$$V _ { I } = - \hbar { ( } \Omega _ { p } ^ { * } e ^ { - i \Delta _ { p } ^ { * } } | 2 \rangle \langle 1 | + \Omega _ { s } ^ { * } e ^ { - i ( \Delta _ { s } + \Delta _ { p } ) t } | 3 \rangle \langle 2 | + \Omega _ { c } ^ { * } e ^ { - i ( \Delta _ { c } - \Delta _ { p } ) t } | 4 \rangle \langle 2 | ) + h. c.,$$

where ∆ p = ω p -ω 21, ∆ s = ω s -ω 32 and ∆ c = ω c -ω 42 are the detunings of the applied fields frequencies and the central frequencies of the corresponding transitions. The density matrix equations of motion can be written as follows

- ˙ ρ 11 = i Ω ρ p 21 -i Ω ρ ∗ p 12 + γ 21 ρ 22 + γ 31 ρ 33 + γ 41 ρ 44 ,
- ˙ ρ 44 = i Ω ρ ∗ c 24 -i Ω ρ c 42 -Γ ρ 2 44 ,
- ˙ ρ 33 = i Ω ρ ∗ s 23 -i Ω ρ s 32 -Γ ρ 1 33 ,
- ˙ ρ 21 = i Ω ∗ p ( ρ 11 -ρ 22 ) + i Ω ρ s 31 + i Ω ρ c 41 -[ γ 21 2 -i ∆ p ] ρ 21 ,

$$\dot { \rho } _ { 3 1 } \ = \ i \Omega _ { s } ^ { * } \rho _ { 2 1 } - i \Omega _ { p } ^ { * } \rho _ { 3 2 } - [ \frac { \Gamma _ { 1 } } { 2 } - i ( \Delta _ { s } + \Delta _ { B } + \Delta _ { p } ) ] \rho _ { 3 1 }, \\ \Gamma \ + \..$$

$$\dot { \rho } _ { 3 2 } \ = \ i \Omega _ { s } ^ { * } ( \rho _ { 2 2 } - \rho _ { 3 3 } ) - i \Omega _ { p } \rho _ { 3 1 } - i \Omega _ { c } ^ { * } \rho _ { 3 4 } - [ \frac { \Gamma _ { 1 } + \gamma _ { 2 1 } } { 2 } - i ( \Delta _ { s } + \Delta _ { B } ) ] \rho _ { 3 2 },$$

$$\dot { \rho } _ { 4 1 } \ = \ i \Omega _ { c } ^ { * } \rho _ { 2 1 } - i \Omega _ { p } ^ { * } \rho _ { 4 2 } - [ \frac { \Gamma _ { 2 } } { 2 } - i ( \Delta _ { c } + \Delta _ { p } - \Delta _ { B } ) ] \rho _ { 4 1 }, \\ \Gamma \ + \..$$

- ˙ ρ 42 = i Ω ∗ c ( ρ 22 -ρ 44 ) -i Ω ρ p 41 -i Ω ρ ∗ s 43 -[ Γ 2 + γ 21 2 -i ( ∆ c -∆ B )] ρ 42 ,

˙

ρ

Γ

1

+

Γ

2

43

=

i

Ω ρ

∗

c

23

-

i

Ω ρ

s

42

-

[

2

+ (

i

∆

s

-

∆

c

+

2

∆

B

)]

ρ

43

,

- ˙ ρ 22 = -( ˙ ρ 11 + ˙ ρ 33 + ˙ ρ 44 ) ,

(2)

where Γ 1 = γ 31 + γ 32 and Γ 2 = γ 41 + γ 42. The parameter γ i 1 ( γ i 2 )( i = 3 4 , ) denotes the spontaneous decay rate from the excited state | i 〉 to the lower states | 1 〉 ( 2 | 〉 ). The polarization vector in the atomic medium is given by

$$\vec { P } ( z, t ) = \chi _ { p } \vec { \varepsilon } _ { p } e ^ { - i ( \theta _ { p } t - k _ { p } z ) } + c. c.$$

Here, χ p is the susceptibility representing the response of the medium to the probe field.

Let us now solve the wave equation for the probe field which can be written as

$$\nabla ^ { 2 } \vec { E } _ { p } - \mu _ { 0 } \varepsilon _ { 0 } \frac { \partial ^ { 2 } \vec { E } _ { p } } { \partial t ^ { 2 } } - \mu _ { 0 } \frac { \partial ^ { 2 } \vec { P } } { \partial t ^ { 2 } } = 0.$$

Inserting equation (3) into equation (4) and using slowly varying approximation, we simplify the equation (4) as

$$\frac { \partial \varepsilon _ { p } } { \partial z } = i 2 \pi \omega _ { p } ( \mu _ { 0 } \varepsilon _ { 0 } ) ^ { 1 / 2 } \varepsilon _ { p } \chi _ { p }.$$

Thus the solution of the equation (4) substituting c = 1 / √ µ ε 0 0 and kp = ω p / c leads the output probe field amplitude to become as

$$\varepsilon _ { p } ( z = l ) = \varepsilon _ { p } ( 0 ) e ^ { \int 2 \pi k _ { p } l x _ { p } }.$$

χ p can be related to the probe transition coherence ρ 21 defined as

$$\chi _ { p } = \frac { n \mu _ { 2 1 } ^ { 2 } \rho _ { 2 1 } } { \hbar { \Omega } _ { p } },$$

where n is density of atoms and ρ 21 is calculated from equation (2). Therefore, the equation (6) reduces to

$$\varepsilon _ { p } ( z = l ) = \varepsilon _ { p } ( 0 ) e ^ { l \frac { \alpha l p _ { 2 1 } y } { 2 \alpha _ { p } } },$$

where α l = 4 π n µ 2 21 kl / ¯ h γ is the resonant absorption. By introducing the normalized susceptibility Sp = ρ 21 γ / Ω p , the output probe field amplitude takes the form

$$\varepsilon _ { p } ( z = l ) = \varepsilon _ { p } ( 0 ) e ^ { i \frac { \alpha } { 2 } s _ { p } }.$$

Finally, the normalized transmission of the probe field is given by

$$T & = \frac { | \varepsilon _ { p } ( z = l ) | ^ { 2 } } { | \varepsilon _ { p } ( 0 ) | ^ { 2 } } = e ^ { - i a l l m | S _ { p } | }. \\ \int _ { \Delta } \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{$$

The normalized susceptibility Sp is clearly a complex quantity that its imaginary part stands for the absorption of the probe field. The intensity region in which the imaginary part of Sp increases with increase the input intensity denotes the RSA region. Equation (10) displays the transmission behavior of the light which is going to be based on the study of the OL properties of the quantum system.

The Z-scan technique is widely used to study the nonlinear refractive index 24 as well as the OL properties of various materials 25 . In experiment, a Z-scan setup includes a laser field with a transverse Gaussian profile focused by using a lens. The sample is then moved along the propagation direction of the focused Gaussian field. It is clear that the sample experiences maximum intensity at the focal point ( z = ) 0 , which gradually decreases in either direction from the focus. The Z-scan technique shows the transmission based on the scanning of the sample position relative to the focal plane of the lens. The incident probe field is a Gaussian laser field with the Rabi frequency Ω p

$$\Omega _ { p } ( z, r, t = 0 ) = \Omega _ { \rho 0 } \frac { w _ { 0 } ^ { 2 } } { w ^ { 2 } ( z ) } e x p [ - ( \frac { 2 r ^ { 2 } } { w ^ { 2 } ( z ) } ) ],$$

where Ω p 0 is the probe Rabi frequency at the focal point (beam waist), w 0 = 0 1 . mm is the beam waist radius at focus, w z ( ) = w 0 [ 1 +( z / z 0 ) 2 1 ] / 2 is the beam radius at z (the distance of the sample from the focal point) and z 0 = π w 2 0 / λ is the diffraction length of the beam. It should be noted that the Z-scan measurements in our work are carried out for the probe field at 800 nm wavelength corresponding to the transition 5 S 1 / 2 ↔ 5 P 1 / 2 .

With calculations of equation (11) numerically for a pulsed Gaussian beam, normalized transmission as a function of position can be obtained as

$$T ( z ) = \frac { 4 } { w ^ { 2 } } \int _ { 0 } ^ { * } t T e x p [ - ( \frac { 2 ^ { 2 } } { w ^ { 2 } ( z ) } ) ] d r.$$

## Results and Discussion

Here, we are going to present our numerical results describing the absorption behavior of the system. We are interested in investigating the SA and RSA regions in the system to provide the appropriate conditions for inducing the OL. All the parameters are scaled by γ which is 2 π × 5 75 MHz for the transition 5 . S 1 / 2 ↔ 5 P 1 / 2 of 87 Rb atoms. Figure 2 shows the absorption of the probe field versus the incident intensity of the probe field for different values of Ω s . The dotted line is for Ω s = 0, the dot-dashed line for Ω s = 0 5 . γ , the dashed line for Ω s = γ and the solid line for Ω s = 2 . The other used parameters γ are Ω p = 0 01 . γ , Ω c = 65 , γ ∆ c = 100 , γ ∆ p = 1 5 . γ , ∆ s = 0 and ∆ B = 0. It is seen that when the coupling field Ω s is switched off, the absorption of the probe field decreases by increasing the intensity of the input probe field and the SA is dominant in the absence of the Ω s . By switching on the Ω s , the RSA region is induced and a peak is generated in the absorption of the probe field, leading to separate the RSA and SA regions. Moreover, the absorption peak enhances by increasing the Ω s , so it makes us possible to control the RSA phenomenon. Generally, in the RSA region (left side of the peak), the absorption increases by growing the intensity of the input laser field. An investigation on Fig. 2 shows that the RSA can switch to the SA for the intense input laser field.

Figure 2. Absorption of the probe field versus intensity of the input probe field for different values of Ω s . The used parameters are Ω p = 0 01 . γ , Ω c = 65 , γ ∆ c = 100 , γ ∆ p = 1 5 . γ , ∆ s = 0, ∆ B = 0, Ω s = 0 (dotted), Ω s = 0 5 . γ (dot-dashed), Ω s = γ (dashed) and Ω s = 2 γ (solid).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01186v1-figure-2928266f-aaf1-45a8-9030-116e127f403d.jpg)

Figure 3. Absorption behavior of the probe field versus intensity of the input probe field for different values of ∆ B . The taken parameters are Ω p = 0 01 . γ , Ω c = 65 , γ ∆ c = 100 , γ ∆ p = 1 5 . γ , ∆ s = 0, Ω s = 2 , γ ∆ B = 0 (dotted), ∆ B = γ (dot-dashed), ∆ B = 2 γ (dashed) and ∆ B = 3 γ (solid).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01186v1-figure-88ab956c-b51b-4bda-8859-f6ef305d18af.jpg)

The extension of the RSA with respect to the SA region is another scenario that we can do by applying the static magnetic field. After making the RSA in the system, it is important that the system maintains the RSA behavior in wider range of the intensity of the input field. The constructive role of the static magnetic field in the RSA is presented in Fig. 3 for ∆ B = 0

Figure 4. Transmission of the probe field versus intensity of the input probe field for Ω s = 0 (dotted), Ω s = 0 5 . γ (dot-dashed), Ω s = γ (dashed) and Ω s = 2 γ (solid). The other used parameters are the same used in Fig. 2 accompanied by α l = 800 . γ

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01186v1-figure-0f3c25e0-7312-45bc-a274-e1f4df112518.jpg)

(dotted), ∆ B = γ (dot-dashed), ∆ B = 2 γ (dashed) and ∆ B = 3 γ (solid). The other taken parameters are Ω p = 0 01 . γ , Ω c = 65 , γ ∆ c = 100 , γ ∆ p = 1 5 . γ , ∆ s = 0 and Ω s = 2 . We result that the RSA region is extended by increasing the magnitude of the static γ magnetic field and the RSA is established in a larger range of the input probe field intensity, which promises the extension of the OL range.

In the following, we investigate the effect of the coupling field on the transmission of the probe field. In Fig. 4, the transmission of the probe field versus the intensity of the input probe field is shown for α l = 800 , γ Ω s = 0 (dotted), Ω s = 0 5 . γ (dot-dashed), Ω s = γ (dashed) and Ω s = 2 γ (solid). The other used parameters are those taken in Fig. 2. Figure 4 shows that in the absence of the coupling field, the transmission of the probe field grows with increase of the input probe field intensity going the system toward transparency. Thus, the system cannot be used as an optical limiter. As proved in Fig. 2, the RSA was induced and intensified by increasing the Ω s . In the RSA domain, the transmission of the probe field keeps constant or even reduces by increasing the coupling field intensity. Hence, it is demonstrated that the OL is coherently induced and controlled in the system. In addition, a bird's eye view of Fig. 4 reveals that increase of the Ω s can lead to decrease the threshold of the OL induced in the system.

Figure 5. Effect of the static magnetic field as ∆ B = 0 (dotted), ∆ B = γ (dot-dashed), ∆ B = 2 γ (dashed) and ∆ B = 3 γ (solid) on the transmission of the probe field versus input intensity of the probe field. The other parameters used are the same as those in Fig. 3 accompanied by α l = 800 . γ

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01186v1-figure-6fc53046-9308-4e13-b009-d616d79bafb2.jpg)

One of the important features that distinguish an optical limiter from the rest is the intensity range in which the optical limiter operates. Here, we show that the applying the static magnetic field can extend the OL range. Figure 5 depicts the transmission of the probe field versus input intensity for different values of the static magnetic field. The taken parameters are the same used in Fig. 3. Figure 5 shows that the OL range can be controlled by the static magnetic field. It is seen that in the

absence of the static magnetic field, the OL is even established in small range of the input intensity. Thus, using the suggested atomic optical limiter, applying the static magnetic field makes the optical devices and sensors safe from damages in a larger range of input intensity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01186v1-figure-c0503698-2f15-4e5b-b094-a64b8b4053d9.jpg)

in

Figure 6. Transmission of the probe field versus input intensity of the probe field for different values of the resonance absorption α l = 200 γ (dotted), α l = 400 γ (dot-dashed), α l = 600 γ (dashed) and α l = 800 γ (solid). The other taken parameters are Ω p = 0 01 . γ , Ω c = 65 , γ ∆ c = 100 , γ ∆ p = 1 5 . γ , ∆ s = 0, Ω s = 2 γ and ∆ B = 3 . γ

Control of the OL threshold is another advantageous of the suggested optical limiter. In Fig. 6, the effect of the resonant absorption, α l , is studied on the transmission of the probe field plotted versus the intensity of the input probe field. Resonance absorption is directly related to the length of the medium and density of atoms. It is observed that by increasing α l , the OL threshold decreases. Decrease of the OL threshold makes it possible that the presented optical limiter can be set to use in optical devices which needs the optical limiters with lower OL threshold.

Figure 7. Slope of the transmission of the probe field versus the static magnetic field and intensity of the input probe field. The used parameters are Ω p = 0 01 . γ , Ω c = 65 , γ ∆ c = 100 , γ ∆ p = 1 5 . γ , ∆ s = 0, Ω s = 2 γ and α l = 800 . γ

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01186v1-figure-9197dc95-4c5d-428d-9a7e-4317bbdf807f.jpg)

In order to gain a deeper insight, the slope of transmission of the probe field, D(T), is displayed in Fig. 7 as a function of the intensity of the input probe field and the static magnetic field. The used parameters are Ω p = 0 01 . γ , Ω c = 65 , γ ∆ c = 100 , γ ∆ p = 1 5 . γ , ∆ s = 0, Ω s = 2 γ and α l = 800 . γ It is worth to note that for the OL range, the slope of the transmission is zero and even negative. Otherwise, the slope of the transmission is positive for the SA region. Figure 7 delineates the behavior of the RSA and corresponding the OL induced in the system as well as the SA for all values of the static magnetic field and intensity of the input field. This figure helps us to determine the OL range needed for different optical devices by selecting the appropriate parameters. The OL line, shown in the Fig. 7, presents the zero slop of transmission in the end of the RSA region. The left side of the OL line specifies the OL range, while the right side determines the SA region.

## Z-scan technique

Finally, we employ Z-scan technique to confirm the validity of the obtained theoretical results. The dip in Z-scan transmission curve corresponds to the OL, while the peak stands for the SA effect. In Fig. 8, Z-scan technique measurements is displayed to investigate the z-dependent transmission for different values of input intensity. As seen in Fig. 8, Z-scan technique demonstrates that the OL happens for the input intensities in which the RSA takes place shown in Fig. 3. In contrast, for the intensities through the SA region, Z-scan profiles doesn't show the OL around the focal point. The theoretical Z-scan results are in good agreement with the results mentioned in Figs. 2-6.

Figure 8. Z-scan measurements of the obtained the OL using input gaussian probe field at 800 nm wavelength for different values of input intensity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01186v1-figure-3d422339-83f9-4367-8467-3d8261f03fe0.jpg)

## Conclusion

In summary, coherent generation and control of the RSA and OL are reported in a four-level Y-type quantum system. It was shown that the RSA is coherently induced by applied laser fields. We showed that, consequently, the OL is coherently induced through the RSA region so that all characteristics of the induced OL such as the intensity range and the threshold intensity can be controlled by either intensity or frequency of the laser fields. In addition, we proved that the static magnetic field has a constructive role in extending the RSA region and the OL range. Besides, it was demonstrated that the OL threshold can be controlled by either increasing the length of the medium or getting the atomic medium denser. Finally the obtained theoretical results was confirmed by Z-scan technique. Our presented scheme can be used in designing the optical limiters with controllable intensity range and OL threshold.

## References

- 1. Li, C., Zhang, L., Wang, R., Song, Y. &amp; Wang, Y. Dynamics of reverse saturable absorption and all-optical switching in C 60. J. Opt. Soc. Am. B 11 , 1356-1360 (1994).
- 2. Tutt, L. W., &amp; Boggess, T. F. A review of optical limiting mechanisms and devices using organics, fullerenes, semiconductors and other materials. Prog. Quantum Electron. 17 , 299-338 (1393).
- 3. Sun, Y. P., &amp; Riggs, J. E. Organic and inorganic optical limiting materials. From fullerenes to nanoparticles. Int. Rev. Phys. Chem. 18 , 43-90 (1999).
- 4. Li, Q., Liu, C., Liu, Z., &amp; Gong, Q. Broadband optical limiting and two-photon absorption properties of colloidal GaAs nanocrystals. Opt. Express 13 , 1833-1838 (2005).
- 5. Ehrlich, J. E., Wu, X. L., Lee, I. Y., Hu, Z. Y., R¨ckel, H., Marder, S. R., &amp; Perry, J. W. Two-photon absorption and o broadband optical limiting with bis-donor stilbenes. Opt. Lett. 22 , 1843-1845 (1997).
- 6. Rao, S. Optical limiting in the presence of simultaneous one and two photon absorption. Optik 157 , 900-905 (2018).
- 7. Venkatram, N., Rao, D. N., &amp; Akundi, M. A. Nonlinear absorption, scattering and optical limiting studies of CdS nanoparticles. Opt. Express 13 , 867-872 (2005).

- 8. Venkatram, N., Kumar, R. S. S., &amp; Narayana Rao, D. Nonlinear absorption and scattering properties of cadmium sulphide nanocrystals with its application as a potential optical limiter. J. Appl. Phys. 100 , 074309 (2006).
- 9. Kobyakov, A., Hagan, D. J., &amp; Van Stryland, E. W. Analytical approach to dynamics of reverse saturable absorbers. J. Opt. Soc. Am. B 17 , 1884-1893 (2000).
- 10. Dong, P., &amp; Tang, S. H. Reverse-saturated absorption and narrow line in the absorption spectrum of a V-type three-level medium. Phys. Lett. A 290 , 255-260 (2001).
- 11. Lepkowicz, R., Kobyakov, A., Hagan, D. J., &amp; Van Stryland, E. W. Picosecond optical limiting in reverse saturable absorbers: a theoretical and experimental study. J. Opt. Soc. Am. B 19 , 94-101 (2002).
- 12. Azzam, S. I., &amp; Kildishev, A. V. Time-domain dynamics of reverse saturable absorbers with application to plasmonenhanced optical limiters. Nanophotonics 8 , 145-151 (2018).
- 13. Sun, Z., Tong, M., Zeng, H., Ding, L., Wang, Z., Dai, J., Bian, G., &amp; Xu, Z. Nanosecond reverse saturable absorption and optical limiting in (Me 4 N) 2 [Cd (dmit)(Sph) 2]. J. Opt. Soc. Am. B 18 , 1464-1468 (2001).
- 14. Mishra, S. R., Rawat, H. S., &amp; Mehendale, S. C. Reverse saturable absorption and optical limiting in C 60 solution in the near-infrared. Appl. Phys. Lett. 71 , 46-48 (1997).
- 15. Pan, H., Chen, W., Feng, Y. P., Ji, W., &amp; Lin, J. Optical limiting properties of metal nanowires. Appl. Phys. Lett. 88 , 223106 (2006).
- 16. Wang, J., &amp; Blau, W. J. Inorganic and hybrid nanostructures for optical limiting. J. Opt. 11 , 024001 (2009).
- 17. Liu, Z., Zhang, B., &amp; Chen, Y. Recent Progress in Two-Dimensional Nanomaterials for Laser Protection. Chemistry 1 , 17-43 (2019).
- 18. Muller, O., Pichot, V ., Merlat, L., &amp; Spitzer, D. (2019). Optical limiting properties of surface functionalized nanodiamonds probed by the Z-scan method. Sci. Rep. 9 , 519 (2019).
- 19. Tintu, R., Nampoori, V. P. N., Radhakrishnan, P., &amp; Thomas, S. Reverse saturable absorption in nano colloidal Ge28Sb12Se60 chalcogenide glass. J. Non-Cryst. Solids 375 , 2888-2891 (2011).
- 20. Kiran, P. P., Reddy, D. R., Maiya, B. G., Dharmadhikari, A. K., Kumar, G. R., &amp; Desai, N. R. Enhanced optical limiting and nonlinear absorption properties of azoarene-appended phosphorus (V) tetratolylporphyrins. App. Opt. 41 , 7631-7636 (2002).
- 21. Li, F., Zheng, Q., Yang, G., Dai, N., &amp; Lu, P. Optical limiting properties of two phthalocyanines using 1064-nm laser in solution. Mater. Lett. 62 , 3059-3062 (2008).
- 22. Allam, S. R. Theoretical study on nonlinear properties of four level systems under nano-second illumination. Laser Phys. , 25 , 055701 (2015).
- 23. Abdullah, M., Bakhtiar, H., Krishnan, G., Aziz, M. S. A., Danial, W. H., &amp; Islam, S. Transition from saturable absorption to reverse saturable absorption of carmoisine dye under low-powered continuous wave laser excitation. OPT LASER TECHNOL , 115 , 97-103 (2019).
- 24. Sheik-Bahae, M., Said, A. A., Wei, T. H., Hagan, D. J., &amp; Van Stryland, E. W. Sensitive measurement of optical nonlinearities using a single beam. IEEE journal of quantum electronics, IEEE J. Quantum Electron. 26 , 760-769 (1990).
- 25. Liu, Z. B., Liu, Y . L., Zhang, B., Zhou, W. Y ., Tian, J. G., Zang, W. P., &amp; Zhang, C. P. Nonlinear absorption and optical limiting properties of carbon disulfide in a short-wavelength region. J. Opt. Soc. Am. B 24 , 1101-1104 (2007).

## Author contributions

M. Mahmoudi conceived the idea of the research and directed the project. All authors developed the research conceptions, analysed, and discussed the obtained results. M. Ghaderi Goran Abad, M. Mahdieh and M. Veisi performed the calculations. M. Ghaderi Goran Abad wrote the paper with major input from M. Mahmoudi. H. Nadjari has a major role in description of the results.

## Additional information

Competing financial interests: The authors declare no competing financial and non-financial interests.
Mohsen Ghaderi Goran Abad Mahsa Mahdieh Mohadeseh Veisi Hamid Nadjari Mohammad Mahmoudi
Coherent control of Optical limiting in atomic systems
Generation and control of the reverse saturable absorption (RSA) and optical
limiting (OL) are investigated in a four-level Y-type quantum system. It is
demonstrated that the applied laser fields induce the RSA and it can be
coherently controlled by either intensity or frequency of the applied laser
fields. The effect of the static magnetic field on the induced RSA is studied
and we obtain that it has a constructive role in determining the intensity
range in which the OL is established in the system. In addition, we find that
the OL threshold can be decreased either by increasing the length of the medium
or by getting the atomic system denser. Finally, Z-scan technique is presented
to confirm our theoretical results. The proposed scheme can be used in
designing the coherent optical limiters with controllable threshold and
intensity range of OL.