Anisotropic nonlinear refractive index measurement of a photorefractive crystal via spatial self-phase modulation

TL;DR

This paper introduces a simple method to measure and control the anisotropic nonlinear refractive index in a photorefractive crystal using spatial self-phase modulation and far-field diffraction analysis.

Contribution

The study presents an accessible technique for accurately measuring and controlling the anisotropic nonlinear refractive index in photorefractive crystals through far-field diffraction patterns.

Findings

Measured saturation intensity and maximum refractive index modification.

Characterized the anisotropic nonlinear response of the crystal.

Demonstrated control of refractive index modification via external electric field.

Abstract

We show that the refractive index modification photoinduced in a biased nonlinear photorefractive crystal can be accurately measured and controlled by means of a background incoherent illumination and an external electric field. The proposed easy-to-implement method is based on the far-field measurement of the diffraction patterns of a laser beam propagating through a self-defocusing medium undergoing spatial self-phase modulation. For various experimental conditions, both saturation intensity and maximum refractive index modification have been measured. We also clearly evidence and characterise the anisotropic nonlinear response of the crystal in the stationary regime.

Tables

Table 1. Table 1: Summary of the fitted parameters extracted using an isotropic saturable model on the data presented on Fig.  4 (a) and Fig.  4 (b).

	Symbols	$\left|\Delta n_{\text{max}}\right|(\times {10}^{-4})$	$I_{\text{sat}}(\mathrm{mW}.{\mathrm{cm}}^2)$
Same $E_{\text{ext}}$	$\circ$	$2.99\pm 0.02$	$9.0\pm 0.2$
Fig. 4(a)	$△$	$3.03\pm 0.08$	$25.7\pm 0.9$
	$◆$	$2.64\pm 0.06$	$74.7\pm 3.3$
Same $I_{\text{sat}}$	$△$	$3.03\pm 0.08$	$25.7\pm 0.9$
Fig. 4(b)	$\circ$	$1.86\pm 0.05$	$23.9\pm 0.9$
	$◆$	$1.09\pm 0.06$	$24.9\pm 2.2$

Equations

$$\Delta n(I)=-\frac{1}{2}{n}_0^3{r}_{33}{E}_{\mathrm{sc}}[I(x,y)],$$

$${\nabla }_{\perp }^2{\varphi }_0+{\nabla }_{\perp }\ln (1+\tilde{I})\cdot {\nabla }_{\perp }{\varphi }_0=|{\mathbf{E}}_{\mathrm{ext}}|{\partial }_x\ln (1+\tilde{I}),$$

$$\Delta n(I)=-\frac{1}{2}{n}_0^3{r}_{33}{E}_{\mathrm{ext}}\frac{1}{1+\tilde{I}}.$$

$$i{\partial }_{z}E=-\frac{1}{2k_0{n}_0}{\nabla }_{\perp }^{2}E-i\alpha E+k_0\Delta n(I)E,$$

$$i{\partial }_{z}A=-\frac{1}{2k_0{n}_0}{\nabla }_{\perp }^{2}A-i\alpha A+k_0\Delta n_{\max }\left(\frac{\tilde{I}}{1+\tilde{I}}\right)A.$$

$$\Delta n(I)=\frac{1}{2}{n}_0^3{r}_{33}{E}_{\mathrm{ext}}\frac{\tilde{I}}{1+\tilde{I}}=\Delta n_{\max }\frac{\tilde{I}}{1+\tilde{I}},$$

$$\Delta {\Psi }_{\text{NL}}(x,y)=k_0{\int }_0^{L_z}\Delta n(r,z)dz,$$

$$\Delta {\Psi }_{\text{NL}}(x,y)=\Delta {\Psi }_0\exp \left(-r^2/w_p^2\right),$$

$$N_{\mathrm{rings}}\simeq \frac{\Delta {\Psi }_0}{2\pi }.$$

Full text

I Introduction

In the past few decades, photorefractive (PR) crystals Yeh (1993) have been extensively used as a suitable medium for nonlinear optics applications Günter and Huignard (2007). With a controllable bias voltage and under incoherent illumination, the photoinduced optical properties of PR crystals can indeed be efficiently engineered. Numerous nonlinear optics experiments have been carried out, mainly focused on the study of the linear and nonlinear light transport properties. On the one hand, non-diffracting optical beams allow to fabricate non-permanent photonic lattices, i.e. arrays of coupled waveguides. It has been widely used to observe various families of optical discrete solitons Fleischer et al. (2003, 2005); Lederer et al. (2008), to study wave properties of quasicrystals Boguslawski et al. (2016) and graphene-like structures Peleg et al. (2007); Song et al. (2019) as well as to evidence weak or strong transverse localisation in disordered lattices Schwartz et al. (2007); Boguslawski et al. (2017). On the other hand, in homogeneous configuration, nonlinear patterns formation dynamics Denz et al. (1998); Marsal et al. (2010); Caullet et al. (2012) as well as quantum hydrodynamical features of light Wan et al. (2007); Sun et al. (2012); Michel et al. (2018) have been investigated.

In order to get better insights on such light transport properties and allow quantitative comparison with theory and numerics, a key experimental parameter to monitor is the amplitude of the photoinduced refractive index modification. In principle, this value can be accurately controlled through an external voltage and an incoherent illumination but, surprisingly, is roughly estimated in most experiments (typically of the order of $10^{-4}$). Moreover, when measured, the calibration method employed and the associated uncertainties are rarely addressed. Recently, Armijo et al. proposed an absolute calibration method to measure the refractive index of photoinduced 1D and 2D photonic lattices Armijo et al. (2014). These measurements are very consistent with theory but specific to lattice configurations in a focusing, highly saturated regime. In a previous work, we also addressed a method based on the displacement of a weak optical defect in a strong optical gaussian environment to extract the maximum refractive index and the saturation intensity Michel et al. (2018). The measurements were in quite good accordance with the expected theoretical values but not easy to implement and time consuming.

When an ultrashort pulse travels in a nonlinear medium new frequencies appear within its spectrum as the input intensity increases, via the self-phase modulation effect Agrawal (2012); Finot et al. (2018). Each new frequency is associated to a 2$\pi$-phase shift and thus to a given value of the nonlinear refractive index. Similarly, when a cw spatially limited wavepacket propagates in a nonlinear medium, new spatial frequencies are generated as the input intensity increases. In this case, it simply consists of rings appearing in the far-field, as illustrated in Fig. 1. This approach already allowed to measure cw nonlinear effects in liquid crystals Durbin et al. (1981); di Pietro et al. (2018), graphene suspensions Shan et al. (2019) as well as hot atomic vapours Santic et al. (2018); Fontaine et al. (2018) to cite a few.

Here, we exploit spatial self-phase-modulation (sSPM) to measure the nonlinear photorefractive response of a Strontium Barium Niobate (SBN) crystal. Both saturation intensity and maximum refractive index modification have been measured with respect to the input intensity, external voltage and incoherent illumination amplitude. We also clearly evidenced and measured the anisotropic nonlinear response of the crystal, in good agreement with reported measurement in lattices Armijo et al. (2014).

The paper is organised as follow. We first briefly describe the nonlinear response of the PR crystal, and the nonlinear propagation of a laser beam polarised along its crystallographic axis. Then we present the sSPM effect, which describes the nonlinear interference pattern formed in the far-field at the output of the nonlinear medium. This effect results from the nonlinear phase accumulated by the laser beam along its propagation. Then we show how this approach allows to measure the absolute nonlinear refractive optical index of the crystal. Experimental and numerical results are compared.

II Photorefractive model and nonlinear CW laser propagation

In this section, we shortly describe the nonlinear photorefractive effect Denz et al. (2003) occuring when a photorefractive crystal is illuminated by a laser beam. In the regions of high intensity, impurities charge carriers go from the valence to the conduction band of the crystal and become mobile. They therefore move to different locations in the crystal, and two transport mechanisms are considered. First, charges can be diffused toward a lower charge concentration place. Second, they can be drifted whether by the internal electric field generated by the photovoltaic effect, and/or by an external electric field. The driven motion of the charge carriers occurs in a preferential direction, which is naturally determined by the direction of the polar axis ($c$-axis) of the crystal. An exeternal electric field applied along the $c$-axis reinforces the directed motion of the charges. Finally, the electrons (or holes) recombine with free sites of the crystal lattice, creating a new distribution of charges, resulting in a spatially distributed space charge electric field. This field modulates the optical refractive index of the crystal, through the linear electro-optic Pockels effect. Depending on the orientation of the external electric field with respect to the $c$-axis, the refractive index change can be either positive (focusing effect) or negative (defocusing effect). For a linear polarisation in the transverse $(x,y)$ plane, parallel to the crystallographic $c$-axis, $z$ being the direction of propagation, the variation of the refractive index in the transverse plane is then given by:

[TABLE]

where $n_{0}=2.36$ in the linear refractive index of the bulk crystal along the $c$-axis. In this case, the electro-optic coefficient is $r_{33}=250$ pm/V and $E_{\rm sc}=|\mathbf{E}_{\rm sc}(I)|$ is the amplitude of the space charge electric field which depends on the intensity $I(x,y)$ of the optical beam.

II.1 Anisotropic description

In a biased photorefractive crystal with an external electric field $\mathbf{E}_{\rm ext}$ applied along the $c$-axis, and in the framework of the band-transport model proposed by Kukhtarev et al. Kukhtarev et al. (1978), the space charge electric field can be expressed in terms of an electrostatic potential $\phi$. It consists in a light induced potential $\phi_{0}$ and an external bias term $-|\mathbf{E}_{\rm ext}|x$, considering the $c$-axis aligned with the $x$-axis. This results in an anisotropic description of the photorefractive effect, leading to $\mathbf{E}_{\rm sc}=\mathbf{E}_{\rm ext}-\boldsymbol{\nabla}\phi_{0}$ Zozulya and Anderson (1995). The quantity $\mathbf{E}_{\rm sc}-\mathbf{E}_{\rm ext}=-\boldsymbol{\nabla}\phi_{0}$ is often referred as a screening field within the crystal, generated by the light induced potential Diebel et al. (2014). Considering that the system is in a steady state, that the photovoltaic effect is negligible, and that the drift effect overcomes diffusion for the charge carriers migration, the potential equation for the electrostatic field $\phi_{0}$ is given by

[TABLE]

where $\boldsymbol{\nabla}_{\perp}$ is the gradient in the transverse plane in cartesian coordinates. $\tilde{I}=I/I_{\rm sat}$ is the laser beam intensity normalised to the saturation intensity (often associated to a dark intensity), which accounts for the ratio between the thermal and the photoinduced excitations. In experiments, an incoherent white light illumination is usually used as a background to artificially create a contribution to the thermal excitation in a controllable way.

II.2 Isotropic approximation

If the spatial dimension of the optical beam is much larger than the spatial scale of the electric field, the system adopts a (1+1) dimensional geometry, with a quasi-homogeneous illumination of the crystal Denz et al. (2003). In this case, the description of the photorefractive effect is slightly different than the previous case and leads to an isotropic 1D description of the phenomenon Christodoulides and Carvalho (1995). The space charge electric field thus comes to a simpler expression, $E_{\rm sc}=E_{\rm ext}/(1+\tilde{I})$. In this model, the refractive index variation reads:

[TABLE]

In the latter expression, the light intensity induces a screening of the refractive index variation. Indeed, when the laser intensity is zero, the absolute refractive index variation $|\Delta n|$ is maximised. On the contrary, this variation tends to zero when $\tilde{I}$ becomes much greater than unity. The screening effect saturates for high enough laser intensities, and the level of saturation is controlled by the parameter $I_{\rm sat}$. Although not strictly valid in the (2+1)D configuration, this isotropic description will be considered for the analysis of the experimental results.

II.3 Nonlinear cw laser propagation

We model the propagation of a laser beam of amplitude $E(x,y,z)$ and intensity $I=|E|^{2}$, polarised along the $c$-axis of the crystal with the nonlinear Schrödinger equation (NLSE):

[TABLE]

where $z$ is the propagation direction, $k_{0}=2\pi/\lambda_{0}$ with $\lambda_{0}$ the vacuum wavelength of the laser beam, $n_{0}$ the base index of refraction of the bulk crystal, $\Delta n(I)$ the nonlinear variation of the refractive index described by equations (1) to (3), and $\alpha$ the linear absorption coefficient of the bulk material. The absorption coefficient has been measured experimentally and is estimated at $\alpha=1.319$ cm*-1* at 532 nm.

Applying a change of variable, that is $E~{}=~{}A\cdot\mathrm{e}^{ik_{0}\Delta n_{\rm max}z}$ with $\Delta n_{\rm max}=\frac{1}{2}n_{0}^{3}r_{33}E_{\rm ext}$ leads to a new propagation equation for $A$ :

[TABLE]

This corresponds to a change of reference for the bulk refractive index, and defines a new expression for the nonlinear variation of the refractive index that reads :

[TABLE]

in which the saturation of $\Delta n$ with respect to the laser intensity appears clearly. The coefficient $\Delta n_{\rm max}$ thus corresponds to the maximum variation of the nonlinear refractive index so that $\Delta n_{\rm max}=\frac{1}{2}n_{0}^{3}r_{33}E_{\rm ext}$.

III Spatial self-phase modulation and nonlinear refractive index measurement

In order to fully characterise the nonlinear index variation induced by a laser beam, we take advantage of the spatial self-phase modulation (sSPM) to measure the accumulated nonlinear phase shift of a Gaussian beam propagating through the SBN crystal. sSPM is a nonlinear effect which refers to the modulation of the phase in the transverse plane of a laser beam propagating in a nonlinear medium. This phase modulation is proportionnal to the transverse intensity profile of the laser beam. Given the relation (6), the variation of the intensity in the transverse plane as well as along the propagation results in a nonlinear phase shift $\Delta\Psi_{\text{NL}}(x,y)$ that can be written as:

[TABLE]

with $L_{z}$ the length of the photorefractive crystal. For a laser beam with a transverse Gaussian profile of radius at $1/{\rm e}^{2}$ denoted $w_{p}$, one can show that, in a non-saturated regime,

[TABLE]

$r$ being the radial coordinate such as $\mathbf{r}=(x,y)$. Consequently, for each point $r_{1}$ on the curve of $\Delta\Psi_{\text{NL}}$, it is always possible to find a second point, $r_{2}$, having the same slope. Hence, fields coming from regions close to $r_{1}$ and $r_{2}$ present the same wavevector and two-wave interferences are expected. If $\Delta\Psi_{0}\gg 2\pi$, the diffraction pattern the far-field is composed of a set of concentric rings. Bright or dark rings, that is to say maximum constructive (resp. destructive) interference, occur when $\Delta\Psi_{\text{NL}}(r_{1})-\Delta\Psi_{\text{NL}}(r_{2})=m\pi$, with $m$ an even (resp. odd) integer. Finally, the number of rings $N_{\mathrm{rings}}$ can be estimated thanks to the relation Durbin et al. (1981):

[TABLE]

In the following, we present the use of sSPM to access experimentally the nonlinear refractive index of SBN by counting the rings appearing in the diffraction pattern in the far-field.

IV Experimental setup

The experimental set-up is depicted in Fig. 1. The nonlinear medium is a $5\times 5\times 20\,\mathrm{mm}^{3}$ SBN:61 crystal, doped with cerium (0.002 %) to enhance its photoconductivity. The crystal is shone by a laser beam (wavelength $\lambda_{0}=532\,\mathrm{nm}$) with a Gaussian transverse profile whose radius at $1/\mathrm{e}^{2}$ is $220\,\mu\mathrm{m}$ at the entrance of the PR crystal, propagating along the $z$ axis. The position of the crystal along $z$ is fixed in order to keep constant the incident radius of curvature of the wavefronts estimated to be $21\,\mathrm{cm}$ on the entrance facet, and considered constant along the propagation of the light through the crystal. The polarisation of the light is maintained horizontal (i.e. along $x$), and parallel to the $c$-axis to maximise the electrooptic coefficient $r_{33}$ of SBN. Two electrodes allow to create an external electric field $E_{\rm ext}$ along the $c$-axis. Changing $E_{\rm ext}$ results in a modification of the maximal nonlinear refractive index $\Delta n_{\mathrm{max}}$, as indicated in Eq.(3). The dynamics of the saturation of the photorefractive effect is controlled by illuminating the crystal with an incoherent white light, providing a control over $I_{\rm sat}$ and so the dark intensity. As explained previously, the white light enhances thermal effects that increase the charge carriers’ mobility, leading to an easier recombination process. The presence of white light thus increases $I_{\rm sat}$, hence one needs a higher laser intensity to saturate the refractive index variation. The illumination remains unmoved over all experiments, only its power is varied, in order to ensure the best possible control over the variation of $I_{\rm sat}$. The maximum accessible white light power is estimated to be $3\,\mathrm{W}$, provided by a halogen lamp of nominal power $60\,\mathrm{W}$. Practically, the nonlinearity is thus controlled by changing the intensity $I$ of the laser, the electric field $E_{\rm ext}$ and the saturation intensity $I_{\mathrm{sat}}$. Finally, the Fourier transform of the outgoing light is imaged on a sCMOS camera placed one meter away from the output facet. This distance is large enough to considerer being in the far-field.

V Results and discussion

Figure 2(a) presents the intensity spatial distribution of the spatial spectrum measured in the far-field ($k_{x}$, $k_{y}$), for various input intensities $I$. A ring structure appears as $I$ increases. The shape is clearly not symmetric. As will be discussed later, this is due to the anisotropic response of the crystal.

First, let us focus on the behaviour along the $k_{x}$ axis which corresponds to the $c$-axis conjugate. We plot the $k_{x}$-profiles on the Fig. 3(a). To improve the signal-to-noise ratio, we integrate the signal over 40 pixels along the $k_{y}$ axis. With increasing optical intensities, more and more rings appear, indicating an increase of the nonlinear phase shift. Typical profiles at $I=9$ mW.cm*-2* and $I=49$ mW.cm*-2* are shown in the Fig. 3(b-c). The gray points count twice the number of rings $N_{\rm ring}$. The measurement of the nonlinear refractive index based on sSPM requires a precise analysis of the curvature of the central spot. Indeed, as detailed in Shen (1984); Deng et al. (2005); Ramirez et al. (2010), the nonlinear phase shift $\Delta\Psi_{0}$ is $2\pi{}N_{\mathrm{rings}}$ for a bright (i.e. concave) central spot and $(2N_{\mathrm{rings}}-1)\pi$ for the dark (i.e. convex) case. Note that in both cases, the contribution to the nonlinear phase shift of the appearance of a new ring is $2\pi$ Durbin et al. (1981), if the curvature of the central spot remains unchanged. For example, in Fig. 3(b), we have $N_{\rm ring}=4$ and a convex curve at $k_{x}=0$ (see inset) which gives $\Delta\Psi_{0}=7\pi$. In the Fig. 3(c), we have $N_{\rm ring}=1$, a concave curve at $k_{x}=0$ and thus $\Delta\Psi_{0}=2\pi$. Finally, one obtains the nonlinear index of refraction $\Delta n$ as a function of the nonlinear phase, $\Delta n=\Delta\Psi_{0}/(k_{0}L_{z})$, with $L_{z}=20\,\mathrm{mm}$.

Moreover, the diameter of the rings saturates after a certain intensity. As discussed above, such a saturation effect is expected for a photorefractive nonlinearity (in both isotropic and anisotropic descriptions). Then, at very high optical intensity (above 300 mW.cm*-2*), the power inside the central spot increases, reducing the contrast of the outer rings. This observation can be explained considering Eq. (6). In the regime where $I\gg I_{\rm sat}$, the nonlinear refractive index is independent of $I$, implying that the medium becomes homogeneous. Therefore, no interference patterns can be observed in the far-field and one gets a clear Gaussian diffracted pattern. This behavior has been also found in our simulations. For a Kerr nonlinearity, the numerics did not reveal such a behavior. This indicates that it is only related to the saturable nature of the PR nonlinearity.

From the rings counting procedure and the curvature analysis at $k_{x}=0$, the absolute value of $\Delta{}n$ is plotted in the Fig. 4 versus the optical intensity $I$ of the incoming beam for different sets of external electric field and white light intensity. In both configurations, we observe an increase of $|\Delta{}n|$ for small $I$ and a clear saturation at high intensity. Let us first consider the case of Fig. 4(a) where $E_{\text{ext}}=1.5\times 10^{5}\,\mathrm{V.m}^{-1}$ is fixed. From left (blue-filled circles) to right (red diamonds), the white light power is increased from $0.36\,\mathrm{W}$ to $1.92\,\mathrm{W}$. The values reached on the plateau, $\Delta n_{\text{max}}$, by the three curves are very close to each other. Moreover when the white light intensity gets higher, the optical intensities needed to reach the plateau gets higher as well. The data are fitted taking into account the isotropic model defined by eq. (6) along the $c$-axis of the crystal. The absolute variation of the nonlinear refractive index is then $|\Delta n(I)|=|\Delta n_{\text{max}}|\times\tilde{I}/(1+\tilde{I})$. The three fitted curves are shown on Fig. 4(a) and the fitted parameters are presented in Tab. 1. The values extracted from the experimental data are in good agreement with the theoretical value, $|\Delta n_{\text{max}}|=2.32\times 10^{-4}$ given by Eq. (6).

We present on Fig. 4(b) the case where the white light is set at $P_{\rm WL}=1.9\,\mathrm{W}$ such as $I_{\rm sat}$ is fixed to $25.7\,\mathrm{mW.cm}^{2}$ and different external electric fields $[0.5,\,1,\,1.5]\times 10^{5}\,\mathrm{V.m}^{-1}$, from left to right. Here, it is clear that the plateau is reached for the same incident optical intensity, which is expected as the saturation intensity is fixed. On the other hand, the maximum value obtained for $\Delta n_{\rm max}$ gets smaller as $E_{\rm ext}$ decreases. The extracted free parameters of the fits, $|\Delta n_{\rm max}|$ and $I_{\rm sat}$, are summarised in Tab. 1.

As a summary, results presented in Figs. 4 (a) and (b) and in Tab. 1 show that we have a control over the dynamics and asymptotic value of the nonlinear refractive index of the crystal with respect to the laser beam intensity. Indeed, we are able to measure an absolute value for the maximal nonlinear index of refraction in good agreement with the theory. Moreover, the estimation of the saturation intensity $I_{\text{sat}}$ seems to be accurate and grows, as predicted, with the white light intensity. Those measurements thus suggest that the spatial self-phase modulation is a suitable tool to predict the behaviour of the photorefractive nonlinearity.

As mentioned previously and shown on Fig. 2(a), the recorded images present a shape radically different from the expected circular rings. To get better understanding on the observed anisotropy, we solved numerically the coupled Eqs. (4) and (2). The results shown in Fig. 2(b) for $I=100$ mW.cm*-2* is, at least qualitatively, in very good accordance with the corresponding experimental image. To compare with, we present in Fig. 2(c), the numerics obtained at the same input intensity in the isotropic approximation, solving directly Eq. (5). This allows us to clearly attribute the non-symmetric shape of the spectrum to the anisotropic response of the crystal, described by Eqs. (1) and (2).

To be more quantitative, we apply the same experimental procedure on profiles obtained along $k_{y}$ to extract the nonlinear refractive index for each input intensity. Although not shown here, a saturation is also observed. We present on Fig. 5 the ratios of $|\Delta n_{\text{max}}|$ along $k_{x}$ over along $k_{y}$ extracted from the fits for all the different sets of parameters $E_{\text{ext}}$ and $I_{\text{sat}}$ (blue-filled circles). Those values are compared to the aspect ratios extracted from the images (Fig. 2(a)), defined as the ratio between the diameter of the outer ring along the $k_{x}$-direction and along the perpendicular direction and averaged over the 10 images having the highest intensity (red-filled diamonds). We observe a very good agreement between the two methods used to estimate the anisotropy of the nonlinear crystal, and the value of 1.6 is coherent with the literature Armijo et al. (2014). Moreover, our simulations indicate a slightly higher anisotropy of 2.0. This measurement suggests that one can take benefit of the spatial self-phase modulation to access and characterise properly the anisotropy of a photorefractive crystal such as SBN. If we want to describe entirely the diffracted pattern in the far-field, the anisotropic description of the PR effect is required. However, if we are only interested in measuring the nonlinear refractive index of the medium along the $c$-axis ($k_{x}$ conjugated), the isotropic approximation is quantitatively pertinent.

VI Conclusion

To conclude, we have presented an absolute measurement of the nonlinear refractive index of a SBN photorefractive crystal by studying diffraction patterns in the far-field taking benefit of spatial self-phase modulation. When a gaussian wavepacket is initially injected in a medium with isotropic nonlinear response, the expected signature of the sSPM consists in the appearance of new spatial frequencies in the far-field with the observations of concentric rings. We have shown experimentally that for the photorefractive crystal, the spectral intensity distribution is not symmetric. By means of numerical simulations, we confirmed that such signal is due to the anisotropic response of the crystal. We analysed the spatial spectrum and measured the nonlinear refractive index versus input optical intensity for various experimental conditions (external voltage and white light illumination). Our measurements are in very good agreement with the isotropic model of the photorefractive effect when we consider the spectrum along $k_{x}$ which is the conjugated of the $c$-axis. Finally, the anisotropy we measured is 1.6.

The proposed experimental technique is quick, easy to implement and might be applicable to any nonlinear crystals. This is of great importance in the context of linear and nonlinear light transport investigations where precise calibration of the photoinduced nonlinear refractive index is required.

It’s worth mentioning that we didn’t discuss the transient regime of the photorefractive response which, in the present case, is of the order of few tens of seconds. During the transient, the refractive index goes from zero to its maximum value. This regime opens new perspectives in probing the nonlinear light dynamics and will be the subject of future studies.

Acknowledgments

The authors thank G. Labeyrie for fruitful discussions. This work has been supported by the the Region Sud, the French government, through the UCAJEDI Investments in the Future project managed by the National Research Agency (ANR) with the reference number ANR-15-IDEX-01 and the European Union’s Horizon 2020 research and innovation programme under grant agreement No 820392 in the PhoQuS project framework.