Acousto-optic modulation in ion implanted semiconductor plasmas having SDDC

The acousto-optic modulation is the convenient and widely used means of controlling intensity and/or phase of propagating radiation [1, 2]. The conventional way of controlling phonons takes advantage of size confinement. These are used to tailor the propagation properties of acoustic and optical phonons, as well as their interactions with optical fields. The concept of transverse modulational instability originates from a space-time analogy that exists when the dispersion is replaced by diffraction and instability of a plane wave in self focusing Kerr medium.

A large number of attempts have been made on the investigation of modulational interactions. It has found that growth rate of instability in strain dependent dielectric (SDDC) semiconductor crystals such as BaTiO₃ is large enough [3, 4, 5, 6, 7, 8]. These modulations have number of applications including the impression of information onto optical pulses, mode-locking, and optical beam deflection [9, 11]. An important field of study in nonlinear acoustics is the amplification/attenuation and frequency mixing of waves in III-V semiconductors because of its immediate relevance to problems of optical communication systems. In most of the cases of nonlinear optical interactions, the SDDC materials are normally ignored. Looking at the potential of the semiconductors in modern optoelectronic technology, the analytical investigations of some basic nonlinear processes in such crystals are of considerable significance due to its vast technological potentialities.

Ion implantation a process in which use of ions is made to dope and modify semiconductor materials. The colloids that act as third species or foreign particles are the result of the implantation of any metal ion inside the medium. Colloidal plasmas are a new and fascinating field of plasma physics. These colloids acquire a negative charge through the sticking of high mobility free electrons on them. The negatively charged colloidal grains (CGs) are assumed to be of uniform size and smaller than both the wavelength under study and the carrier Debye length [12, 13]. The high mobility charge carrier makes diffusion effects even more relevant in semiconductor technology as they (charge carriers) travel significant distances before recombining. Therefore inclusion of carrier diffusion in theoretical studies of nonlinear wave-wave interactions seems to be very important from the fundamental as well as application view points and thus attracted many researchers in the last decades [14].

Present analysis is based on coupled mode theory for investigating the effects of negatively charged CGs on the acousto-optic modulation due to parametric four wave mixing process in ion implanted semiconductors having SDDC. The intense laser beam electrostrively generates an acoustic wave within the semiconductor medium that induces an interaction between the free carriers (through electron plasma wave) and the acoustic phonons (through material vibration). This interaction induces a strong space charge field that modulates the pump beam. Thus the optical and acoustic waves present in an acousto-optic modulator can be strongly amplified through nonlinear optical pumping. The presence of charged CGs in semiconductor plasma medium add new dimensions to the analysis presented in n-doped semiconductors with strain dependent dielectric constants (SDDC). It is found that the presence of colloidal grains (CGs) plays an effective role in changing the threshold intensity and effective gain constant.

The well-known hydrodynamic description of semiconductor plasmas has been considered to study the acousto-optic (AO) modulation in ion-implanted n-type semiconductor plasma having SDDC (for which k_al ≪ 1, k_a and l being the acoustic wave number and mean free path of an electron, respectively). In order to study the crystal, a uniform (|k₀| ≈ 0 and homogeneous optical pump field E0→=x^E0expiω0t$\begin{array}{} \displaystyle \overrightarrow {{E_0}} = \hat x{E_0}\exp \left( {i{\omega _0}t} \right) \end{array}$ that propagates with parametrically generated acoustic wave within the medium is considered. Since these acoustic grating results in a proportional refractive index variation through the medium photo- elastic response, the incident optical field will be diffracted by this grating to generate an additional field within the medium. The diffracted beam is either frequency up-shifted (anti –stokes mode) or down-shifted (stokes mode) depending on the orientation of the incident wave. In presence of strain, stokes and anti stokes mode can be coupled over a long interaction path. This coupled wave propagates as a solitary wave form in the non dispersive regime of acoustic wave and can be amplified under appropriate phase matching condition.

We proceed with the following basic equations describing acousto-optic modulation under one dimensional configuration (along x-axis):

∂v0∂t+υv0=−eE0.m$$\begin{array}{} \displaystyle \frac{{\partial {v_0}}}{{\partial t}} + \upsilon {v_0} = - \frac{{e{E_{0.}}}}{m} \end{array}$$(1)

∂v1∂t+υv1+v0∂∂xv1=−eE1m−kBTmn0.∇n1$$\begin{array}{} \displaystyle \frac{{\partial {v_1}}}{{\partial t}} + \upsilon {v_1} + \left( {{v_0}\frac{\partial }{{\partial x}}} \right){v_1} = - \frac{{e{E_1}}}{m} - \frac{{{k_B}{T_{}}}}{{m{n_0}}}.\nabla {n_1} \end{array}$$(2)

∂n1∂t+δdn0∂v1∂x+v0∂n1∂x+D∂2n1∂x2=0$$\begin{array}{} \displaystyle \frac{{\partial {n_1}}}{{\partial t}} + {\delta _d}{n_0}\frac{{\partial {v_1}}}{{\partial x}} + {v_0}\frac{{\partial {n_1}}}{{\partial x}} + D\frac{{{\partial ^2}{n_1}}}{{\partial {x^2}}} = 0 \end{array}$$(3)

in which diffusion coefficient

D=kBTeμ$$\begin{array}{} \displaystyle D = \frac{{{k_B}T}}{e}\mu \end{array}$$(4)

Equations (1) and (2) represent the zeroth and first order oscillatory fluid velocities under influence of the respective fields. m and υ represents the effective mass and momentum transfer collision frequency of electrons. Equation (3) is the continuity equation, where n₀ and n₁ are the equilibrium and perturbed carrier densities. In III-V semiconductor plasmas (i.e. n_0e ≈ n_0h ≈ n₀) embedded with CGs, the electrons and holes from all directions colloidal with CGs and get stick onto them. However, greater number of electrons stick onto the CGs as compared to holes in given time interval. The charge imbalance parameter is define as δd=Zden0dn0h.$\begin{array}{} \displaystyle \delta _{d} = \frac{{{Z_d}e{n_{0d}}}}{{{n_{0h}}}}. \end{array}$ The plasma is quasi neutral and the conservation of particle number density must always holds. Thus charged should be hold:n_0ee = n_0he – Z_den_0d

where n_0d is concentration of CGs and Z_d is charge state of colloids [15, 16, 17]. In equation (4), μ=e/mv$\begin{array}{} \mu\left(={^{e_/}}\!_{mv}\right) \end{array}$ is the electron mobility, k_B is the Boltzmann’s constant and T the electron temperature in ⁰K.

∂2u∂t2=cρ∂2u∂x2−(ϵgE0)ρ⋅∂E1∗∂x$$\begin{array}{} \displaystyle \frac{\partial^2 u}{\partial t^2} = \frac{c}{\rho} \frac{\partial^2 u }{\partial x^2} - \frac{(\epsilon g E_0)}{\rho} \cdot \frac{\partial E_1^*}{\partial x} \end{array}$$(5)

∂E1∂x=−n1eε+gE0.∂2u∗∂x2$$\begin{array}{} \displaystyle \frac{{\partial {E_1}}}{{\partial x}} = - \frac{{{n_1}e}}{\varepsilon } + g{E_{0.}}\frac{{{\partial ^2}{u^ * }}}{{\partial {x^2}}} \end{array}$$(6)

Pao=−εgE0.Δu∗$$\begin{array}{} \displaystyle {P_{ao}} = - \varepsilon g{E_{0.}}\Delta {u^ * } \end{array}$$(7)

The equation (5) describes the lattice displacement in an ion implanted semiconductor plasma in which ρ, u and c being mass density of the crystal, lattice displacement and crystal elastic stiffness, respectively. The migration of charge that leads to strong space charge field E⃗₁ is determined by the Poisson’s equation (6) in which D⃗ = εE⃗ (1 + gS) and ε = ε₀.ε_s where D⃗ is the electric displacement, ε_s is dielectric constant in absence of any strain S(= du/dx) and g(= ε_s/3) is coupling constant due to SDDC. The electrostatic force thus produced is the origin of acousto-optic strain within the medium given by equation (7), the nonlinear polarization P_ao.

Physically in acousto-optic modulation process a carrier density perturbation is created in the medium under the influence of a strong pump beam, which is associated with phonon mode and varies as the acoustic frequency. The equation for density fluctuation of the coupled electron plasma wave in ion implanted semiconductor is obtained from equations (1) to (7) using linearized perturbation theory as

∂2n1∂t2+υ∂n1∂t+υD∂2n1∂x2+ω¯p2n1+ikPaoeδdωp2=ik±n1E¯b$$\begin{array}{} \displaystyle \frac{{{\partial ^2}{n_1}}}{{\partial {t^2}}} + \upsilon \frac{{\partial {n_1}}}{{\partial t}} + \upsilon D\frac{{{\partial ^2}{n_1}}}{{\partial {x^2}}} + \bar \omega _p^2{n_1} + \frac{{ik{P_{ao}}}}{e}{\delta _d}\omega _p^2 = i{k_ \pm }{n_1}{\bar E_b} \end{array}$$(8)

In which ω¯p2=δdωp2+k±2kBTm$\begin{array}{} \displaystyle \bar \omega _p^2 = {\delta _d}\left( {\omega _p^2 + k_ \pm ^2\left( {\frac{{{k_B}T}}{m}} \right)} \right) \end{array}$ being modified plasma frequency, ω<!-- ? -->p2=n0e2mε<!-- e -->$\begin{array}{} \omega _p^2 = \left( {\frac{{{n_0}{e^2}}}{{m\varepsilon }}} \right) \end{array}$ is electron plasma frequency and E¯b=eE0m.$\begin{array}{} \displaystyle {\bar E_b} = \frac{{e{E_0}}}{m}. \end{array}$

The pump beam is thus phase modulated by the density perturbations to produce enforced disturbances at the upper (ω_a + ω₀) and lower (ω_a – ω₀) sideband frequencies. The higher order components are filtered out by assuming a long interaction path (by considering the crystal to be of infinite extent). The modulation process under consideration must also fulfill the phase matching conditions k_± = k_a ± k_a and ω_± = ω_a ± ω_a known as momentum and energy conservation relation, respectively under spatially uniform laser irradiations |k₀| ≈ 0 such that |k_±| = |k_a ± k₀| ≈ |k_a| = k (say). The density modulation oscillating at the upper(ω_a + ω₀) and lower(ω_a – ω₀) sideband frequencies can be represented after simplification as

n1ω±,k±=−ikPaoeδdωp2ω¯p2−ω±2−υk±2D−iυω±−ik±E¯b−1$$\begin{array}{} \displaystyle {n_1}\left( {{\omega _ \pm },{k_ \pm }} \right) = - \frac{{ik{P_{ao}}}}{e}{\delta _d}\omega _p^2{\left[ {\bar \omega _p^2 - \omega _ \pm ^2 - \upsilon k_ \pm ^2D - i\upsilon {\omega _ \pm } - i{k_ \pm }{{\bar E}_b}} \right]^{ - 1}} \end{array}$$(9)

The density perturbations oscillating at the forced frequency in equation (9) are obtained under quasi-static approximation and by neglecting the Doppler shift under the assumption that ω₀ ≫ υ ≫ k₀v₀. We also neglected the contribution of transition dipole momentum in the analysis of modulational instability to study the effect of nonlinear current density.

The induced nonlinear current densities for upper and lower sidebands may be expressed as

J+ω+,k+=−eD∂n1∂xω+,k+$$\begin{array}{} \displaystyle {J_ + }\left( {{\omega _ + },{k_ + }} \right) = - eD\frac{{\partial {n_1}}}{{\partial x}}\left( {{\omega _ + },{k_ + }} \right) \end{array}$$(10a)

and

J−ω−,k−=−eD∂n1∂xω−,k−$$\begin{array}{} \displaystyle {J_ - }\left( {{\omega _ - },{k_ - }} \right) = - eD\frac{{\partial {n_1}}}{{\partial x}}\left( {{\omega _ - },{k_ - }} \right) \end{array}$$(10b)

In centrosymmetric system, the four-wave parametric interaction involving the incident pump, the upper and lower side band signals and induced acousto-optical idler wave characterized by the cubic nonlinear susceptibility tensor effectively results in the modulational instability of the pump. Hence the induced cubic nonlinear optical polarization at the modulated frequencies

P_NL(ω_±, k_±) as the time integral of the nonlinear current density J_NL(ω_±, k_±) as we have

PNLω±,k±=∫JNLω±,k±dt$$\begin{array}{} \displaystyle {P_{NL}}\left( {{\omega _ \pm },{k_ \pm }} \right) = \int {{J_{NL}}\left( {{\omega _ \pm },{k_ \pm }} \right)dt} \end{array}$$(11)

The effective nonlinear third order polarization has contribution from both individual side band and can be represented as

Peffω±,k±=PNLω+,k++PNLω−,k−$$\begin{array}{} \displaystyle {P_{eff}}\left( {{\omega _ \pm },{k_ \pm }} \right) = {P_{NL}}\left( {{\omega _ + },{k_ + }} \right) + {P_{NL}}\left( {{\omega _ - },{k_ - }} \right) \end{array}$$(12)

The effective nonlinear third order polarization after algebraic simplification as

Peff=−2νϵ2g2k4Dδdωp2|E0|2E1ρ(ωa2−ka2νa2)(Δ2−k2E¯b2+ω02ν2)+(2iE¯bΔ)(Δ2−k2E¯b2+ω02ν2)2+4k2E¯b2Δ2$$\begin{array}{} \displaystyle P_{eff} = - \frac{2 \nu \epsilon^2 g^2 k^4 D \delta_d \omega_p^2 |E_0|^2 E_1}{\rho (\omega_a^2 - k_a^2 \nu_a^2)} \left [ \frac{(\Delta^2 - k^2 \bar{E}_b^2 + \omega_0^2 \nu^2) + (2i\bar{E}_b \Delta)}{(\Delta^2 - k^2 \bar{E}_b^2 + \omega_0^2 \nu^2)^2 + 4k^2 \bar{E}_b^2 \Delta^2}\right] \end{array}$$(13)

Thus the effective acousto-optic nonlinear susceptibility of medium including the contribution due to the drift and diffusion of charge carriers can be obtained as

ξeff=−2νϵ2g2k4Dδdωp2ρ(ωa2−ka2νa2)(Δ2−k2E¯b2+ω02ν2)+(2iE¯bΔ)(Δ2−k2E¯b2+ω02ν2)2+4k2E¯b2Δ2$$\begin{array}{} \displaystyle \xi_{eff} = - \frac{2 \nu \epsilon^2 g^2 k^4 D \delta_d \omega_p^2 }{\rho (\omega_a^2 - k_a^2 \nu_a^2)} \left [ \frac{(\Delta^2 - k^2 \bar{E}_b^2 + \omega_0^2 \nu^2) + (2i\bar{E}_b \Delta)}{(\Delta^2 - k^2 \bar{E}_b^2 + \omega_0^2 \nu^2)^2 + 4k^2 \bar{E}_b^2 \Delta^2}\right] \end{array}$$(14)

Equation (14) may be separated into real and imaginary parts (ξ_eff.) = (ξ_eff. + i(ξ_eff.)) as,

[ξeff]real=−2νϵ2g2k4Dδdωp2ρ(ωa2−ka2νa2)(Δ2−k2E¯b2+ω02ν2)(Δ2−k2E¯b2+ω02ν2)2+4k2E¯b2Δ2$$\begin{array}{} \displaystyle [\xi_{eff}]_{real} = - \frac{2 \nu \epsilon^2 g^2 k^4 D \delta_d \omega_p^2 }{\rho (\omega_a^2 - k_a^2 \nu_a^2)} \left [ \frac{(\Delta^2 - k^2 \bar{E}_b^2 + \omega_0^2 \nu^2)}{(\Delta^2 - k^2 \bar{E}_b^2 + \omega_0^2 \nu^2)^2 + 4k^2 \bar{E}_b^2 \Delta^2}\right] \end{array}$$(15)

and

[ξeff]imag.=−2νϵ2g2k4Dδdωp2ρ(ωa2−ka2νa2)2E¯bΔ(Δ2−k2E¯b2+ω02ν2)2+4k2E¯b2Δ2$$\begin{array}{} \displaystyle [\xi_{eff}]_{imag.} = - \frac{2 \nu \epsilon^2 g^2 k^4 D \delta_d \omega_p^2 }{\rho (\omega_a^2 - k_a^2 \nu_a^2)} \left [ \frac{2\bar{E}_b \Delta}{(\Delta^2 - k^2 \bar{E}_b^2 + \omega_0^2 \nu^2)^2 + 4k^2 \bar{E}_b^2 \Delta^2}\right] \end{array}$$(16)

In order to explore the possibility of the modulational amplification, the nonlinear steady state growth rate of the modulated waveform of the pump exceeding a threshold value is obtained through the relation

[geff]ao=−k2ϵ[ξeff.]imga.|E0|2$$\begin{array}{} \displaystyle [g_{eff}]_{ao} = - \frac{k}{2 \epsilon} [\xi_{eff.}]_{imga.}|E_0|^2 \end{array}$$(17)

Thus from equations (16) and (17), in order to attain a growth of the modulated signal, it may be infer that [ξ_eff]_imag. should be negative. Thus condition for achieving a positive growth rate is as follows:

k2E¯b2>(Δ2+ω02ν2),$$\begin{array}{} \displaystyle k^2 \bar{E}_b^2 > (\Delta^2 + \omega_0^2 \nu^2), \end{array}$$(18)

It is evident from above discussions that not only the presence of particle diffusion is an necessity to induce instability but also the value of applied pump intensity must be well above the threshold defined by equation (18). The threshold value of the pump amplitude require for the onset of the modulational amplification is obtained as

Eth=mek⋅Δ2+ω02ν21/2$$\begin{array}{} \displaystyle E_{th} = \frac{m}{ek} \cdot \left( \Delta^2 + \omega_0^2 \nu^2 \right)^{1/2} \end{array}$$(19)

The corresponding pump intensity can be obtained by using the relation

Ith=12ηϵ0c|E0th|2.$$\begin{array}{} \displaystyle I_{th} = \frac{1}{2} \eta \epsilon_0 c |E_{0th}|^2. \end{array}$$(20)

The numerical calculations are performed for the n-type semiconductor sample (BaTiO₃) at 300 K duly irradiated by 10.6 μm CO₂ laser.

The following material parameters have been considered as follows: m = 0.0145m₀ (m₀ being the free electron rest mass), m_d = 1.67 × 10^–27kg · ϵ_s = 2000, ρ = 4 × 10³kg · m^–3, η = 2.4n₀ = 10²⁵m^–3, ν_e = 5 × 10¹¹s^–1, ω₀ = 1.78 × 10¹³s^–1, ω_a = 1.6 × 10¹³s^–1, and ν_a = 3 × 10³m · s^–1.

The threshold characteristics are illustrated in Figures 1. Figure 1 shows the variation of threshold intensity I_th with wave number (equation 20). The I_th decreases abruptly as the k increases in ion implanted semiconductors with chosen values of charge imbalance parameter δ_d. It can be seen that at k ≈ 1 × 10⁶m^–1 the intensity is I_th ≈ 2.451 × 10¹³W · m^–2. The I_th decreases abruptly with wave number k up to k ≫ 4 × 10⁷m^–1 and then afterwards decreases slowly. It is evident from the figure that as the fraction of negative charge stick on to the CGs δ_d, decreases (i.e. δ_d ≈ 1 > 0.90 > 0.80), the I_th gets lowered.

Fig. 1

The nature of dispersion arising due to third order real effective susceptibility [χ_eff.]_real for ion implanted semiconductor plasmas has been displayed in Figure 2 with respect carrier concentration[ in terms of ωp/ω0$\begin{array}{} \displaystyle ^{\omega_{p}}/_{\omega_{0}} \end{array}$ ] whenE₀ = 5 × 10⁷Vm^–1 (equation 15). It can be observed that there exists a distinct anomalous parametric dispersion regime that varies with n₀ and δ_d both. The [χ_eff.]_real, can be both positive and negative under the anomalous regime. For ωp/ω0$\begin{array}{} \displaystyle ^{\omega_{p}}/_{\omega_{0}} \end{array}$ < 8.136, [χ_eff.]_real is a positive quantity which initially increases and then achieves a maximum value at about ωp/ω0$\begin{array}{} \displaystyle ^{\omega_{p}}/_{\omega_{0}} \end{array}$ ≈ 5.88. A slight turning in ωp/ω0$\begin{array}{} \displaystyle ^{\omega_{p}}/_{\omega_{0}} \end{array}$ beyond this point causes a sharp fall in [χ_eff.]_real, which first vanishes at near ωp/ω0$\begin{array}{} \displaystyle ^{\omega_{p}}/_{\omega_{0}} \end{array}$ ≈ 10. After this resonance condition[χ_eff.]_real, decreases very sharply and achieves a minimum value.

Fig. 2

The variation of effective gain [g_eff.]_ao with pump electric field E₀ is depicted in Figure 3 when k₀ = 5 × 10⁷m^–1, D = 0.896, n₀ = 10²⁶m^–3 (equation 17). It may be inferred for this set of data that at initially the [g_eff.]_ao increases abruptly with increase in E₀ and attains maximum value at E₀ < 2.5 × 10⁷Vm^–1. Beyond this point E₀ > 4 × 10⁷Vm^–1 the [g_eff.]_ao decreases slowly and almost constant in further increase. On changing the fraction δ_d, the [g_eff.]_ao shows anomalous behavior as shown in figure. Figure 4 shows the effective gain [g_eff.]_ao curve with wave number k₀. The [g_eff.]_ao increase gradually as increase in k₀. On decrease in value of δ_d, the [g_eff.]_ao increases slightly

Fig. 3

Fig. 4

In present analysis, we have analytically studied the influence of CGs on the threshold intensity, effective susceptibility and effective gain constant. It is observed from the study that significant change in threshold and gain characteristics when the charge imbalance parameter is slightly changed. The presence of CGs plays an effective role in changing the threshold intensity and effective gain constant.