Modelization and optimized speckle detection scheme in photorefractive self-referenced acousto-optic imaging

Fabrice Devaux; Jean-Pierre Huignard; François Ramaz

doi:10.1364/OE.22.010682

1. Introduction

Optical imaging through thick biological tissues of several centimeters requires indirect methods in order to retrieve a millimeter resolution over the entire volume. The major drawback comes from scattering that generates multiple optical paths, so that the origin of the photons that propagate through the medium is lost. Many attempts have been made over the past few decades based on flux techniques such as diffuse optical tomography [1, 2], but the extraction of the signal is not straightforward, since it requires solving an inverse scattering problem. A more recent concept couples light with an ultrasonic beam within the medium [3]. Dealing with the acousto-optic (AO) effect, the ultrasound (US) locally modulates the optical path, and thus the phase of light crossing the US field. This phase modulation mainly creates two weak optical sidebands (< 1%), shifted from the career frequency ω_L of the US frequency +/ − ω_US. The magnitude of these so-called tagged photons depends on the optical properties of the medium, such as absorption and scattering. A 3D image of the optical properties of the medium is thus performed with a scan of the US beam. Mainly two approaches are implemented to detect the optical sidebands. The first one consists in a highly spectrally resolved flux measurement with a confocal Fabry-Perot interferometer [4], or in spectral holeburning filtering at liquid helium temperature [5]. A second approach uses a coherent detection of the tagged photons, initially performed by digital off-axis holography with camera detection (CCD, CMOS) [6, 7]. This coherent detection requires accurate sampling of the speckle and the use of a fast camera, since the speckle decorrelates in living tissues on a millisecond time scale or less [8]. For these reasons, wavefront adaptive setups have been implemented with photorefractive (PR) crystals and large-area single detectors, operating much faster than a camera [9]. In dynamic PR holography, a refractive index grating is recorded within the volume of the crystal, resulting from the interference between the tagged photons and a plane wave off-axis reference beam. Consequently, the reference diffracts a spatial replica of the tagged photons’ wavefront, which is time averaged according to the relative response time τ_PR of the crystal (1–20ms), depending on light irradiance (a few 100mW/cm²) [10]. The use of a reference beam to record the hologram generally brings an important coherent background noise (beam fanning), which limits the signal-to-noise ratio (SNR) of detection, since the AO signal is less than two orders of magnitude compared to the DC-flux components when dynamic axial resolution along the ultrasonic beam is considered. In this paper, we will show that AO imaging with PR detection can also be performed without applying a reference beam onto the crystal. This type of interaction was demonstrated by Kamshilin et al. [11] for application of PR correlation filtering to the detection of vibrating rough surfaces or vibrating multimode fibers. The proof of principle for AO imaging is illustrated by Benoità la Guillaume et al. [12] with a dynamic 2D image of ink inclusions embedded within 40mm of a multiply scattering phantom. The PR crystal used is the sillenite type Bi₁₂SiO₂₀ (BSO), known for its high photosensitivity in conventional holographic experiments in the blue-green spectral region. We outline that this configuration is much simpler to implement after having optimized the light-collection efficiency and the speckle size incident on the crystal. The experimental setup detailed in Ref. [12] is shown in Fig. 1. We apply an excitation burst (carrier at 2.23 MHz) on the US transducer to get a good axial resolution of 1.35 mm while the lateral resolution, given by the tranverse profile of the US, is 2 mm. The used gel phantom has a reduced scattering coefficient μ′_s of about 10 cm⁻¹ with two included china ink cylinders to simulate local absorbers. The detection of the AO signal is performed with a large area photodiode of 3.6×3.6mm² placed behind the optical system (L₃L₄) that improves the collection of the speckle light. As analyzed in [12], the self referenced interaction can be interpreted as follows: in absence of US the incident speckled wavefront creates in the PR BSO a photoinduced space charge field which modulates the refractive index through the diffusion of carriers at no applied electric field. This complex phase volume hologram results from the interference of the incident speckle field with self-diffracted beam fanning noise fields. The speckle pattern photo-induces micrometer size random refractive index dots within the BSO volume. The operating conditions of the self-referenced PR holography for acousto-optic imaging are explained as follows: when the US transducer is off, the incident speckle field is perfectly phase matched with the recorded phase volume hologram. It thus results an optimum coherent diffraction on the detector. When the US transducer is now operating over short pulse duration the incident speckle structure is changed and it diffracts a wavefront which does not perfectly match with the complex phase hologram recorded in the volume of the BSO. In such conditions the diffracted intensity on the detector is reduced or enhanced as shown in [12] and permits an acousto-optic imaging of a phantom embedded in the scattering medium. It was also shown experimentally in [12] that in the case of self referenced holographic configuration an optimum signal detection is obtained when the two following conditions are fulfilled: first, the speckle grain size recorded in the PR crystal must be of the order of the Debye screening length. Second, the position of the detector in the transverse plane must be optimized due to the spatial frequency response of the photorefractive effect when speckle is involved. The objective of the paper is to develop hereafter an analytic model and numerical simulations which confirm these optimum conditions for signal detection in the conditions of self referenced speckle beam holography with a PR crystal operation at zero electric field.

Fig. 1 Experimental setup.

Download Full Size | PDF

2. Theoretical model and numerical simulations

According to the experiment described in [12], we first define the incident speckle field (i.e. the reference speckle field) and its intensity inside the BSO crystal respectively as: S_r(r) and I_r(r)= |S_r(r)|², where r is the position vector inside the PR medium. This random distribution of coherent light intensity I_r photoinduces an internal space charge electric field in the whole volume of the PR crystal noted E_sc(r). When the classical approximations are made in the PR model where electrons are considered as the single specy of charge carriers [13], the Ohm’s law governing the current density and the rate equations of the PR charge carriers lead to the following differential equation [14, 15] :

E_{sc} (r) - \frac{1}{k_{0}^{2}} \nabla (\nabla \cdot E_{sc} (r)) = - \frac{k_{B} T}{e} \nabla ln [I_{r} (r) + I_{d}],

where ∇ is the nabla operator, k_B is the Boltzmann constant, T is the temperature, e is the elementary charge, k₀ is the Debye screening wave number (typically for BSO k₀ ≃ 4.65μm⁻¹ [16]) and I_d is a uniform background illumination responsible of dark current photoconductivity. Using the Fourier transform properties, Eq. (1) becomes:

ℱ 𝒯 {E_{sc} (r)} = - \frac{k_{B} T}{e} \frac{i k}{1 + \frac{{| k |}^{2}}{k_{0}^{2}}} ℱ 𝒯 {ln [I_{r} (r + I_{d})]},

where ℱ𝒯 denotes the 3D Fourier transform and k is the wavevector. Assuming that the variation of the PR variables along the propagation axis can be neglected and a two-dimensional illumination where I_r(x, y) ≪ I_d, Eq. (2) may be rewritten as:

{\begin{array}{l} {\tilde{E}}_{{sc}_{x}} (k_{t}) = & - \frac{i α k_{x}}{1 + \frac{{| k_{t} |}^{2}}{k_{0}^{2}}} {\tilde{I}}_{r} (k_{t}) \\ {\tilde{E}}_{{sc}_{y}} (k_{t}) = & - \frac{i α k_{y}}{1 + \frac{{| k_{t} |}^{2}}{k_{0}^{2}}} {\tilde{I}}_{r} (k_{t}) \end{array}

where Ẽ_{sc_x}(k_t), Ẽ_{sc_y}(k_t) and Ĩ_r(k_t) are respectively the 2D Fourier transforms, calculated in the (x, y) plane, of the transverse components of the space charge field and of the speckle intensity. k_t = k_x + k_y is the transverse component of the wavevector and

α = \frac{k_{B} T}{I_{d} e}

is a constant. These equations show that each internal space charge field component corresponds to the product of an odd function, presenting extreme values for ±k₀, with Ĩ_r(k_t) that corresponds to the autocorrelation function of the speckle intensity.

Following the analytical expression of the internal space charge field given by Eq. (3), numerical calculations are performed to illustrate the electric field components. Figures 2(a) and 2(b) represent respectively the intensity of a typical speckle pattern in the (x, y) spatial domain and its spectrum in the (k_x, k_y) spatial frequency domain. Figure 2(c) depicts the modulus of the Fourier transform of the speckle intensity |Ĩ_r(k_t)|. The calculations are performed at the wavelength λ = 0.532μm and with the sampling steps $Δ x = Δ y = \frac{λ}{4}$ in the (x, y) plane. The speckles are simulated using a complex field amplitude with a real part and an imaginary part that obey to a gaussian random distribution. The $\frac{1}{e^{2}}$ half width σ_s of the speckle grains is adjusted by filtering, in the Fourier domain, the spectrum of the speckle with a gaussian filter. We verified that the statistical property of the intensity fluctuations are in agreement with the thermal statistic of speckles [17] where the probability density function for intensity verifies the distribution law $p (I) = \frac{1}{〈 I 〉} e^{- \frac{I}{〈 I 〉}}$ with 〈I〉 the mean intensity. We also verified that standard deviation of the spatial fluctuations of intensity is equal to the mean intensity : σ_I = 〈I〉. Speckles, where σ_s is set in the range [0.1μm – 3μm], are considered here. In Fig. 2(a), the mean size of the speckle grains is σ_s = 0.21μm that corresponds to the Debye screening length such that $σ_{s} = l_{D} = \frac{1}{k_{0}}$ . In the (k_x, k_y) spatial frequency domain, the axes are normalized by k₀, the Debye screening wave number. In Fig. 2(c), as the Fourier transform of the speckle intensity exhibits a strong Dirac peak at the zero spatial frequency, the dynamic of the color bar is adjusted to make visible the pedestal of |Ĩ_r|. Figure 3 shows the shape of the two transverse components of the internal space charge field in the (k_x, k_y) and in the (x, y) domains. According to Eq. (3), we assume that the pedestal of the autocorrelation function of the speckle intensity has a gaussian shape with a $\frac{1}{e^{2}}$ half width σ_k related to the size σ_s of the speckle grains. We can easily demonstrate that the amplitude of the space charge field components exibit maxima for the spatial frequencies k_x = ±σ_k and k_x = ±k₀ when the size of the speckle grains is respectively greater and smaller (or equal) than the Debye screening length.

Fig. 2 (a) Intensity of a typical speckle pattern in the (x, y) spatial domain where the $\frac{1}{e^{2}}$ half width of the speckle grains is σ_s = 0.21μm. (b) Intensity of the spatial frequency spectrum of the speckle field. (c) Modulus of the Fourier transform of the speckle intensity. In the spatial frequency domain, the axes are normalized by k₀, the Debye screening wave number.

Download Full Size | PDF

Fig. 3 (a) and (c) Modulus of the transverse components of the space charge electric field in the (k_x, k_y) spatial frequency domain. (b) and (d) Modulus the transverse components of the space charge electric field in the (x, y) spatial domain.

Download Full Size | PDF

The crystallographic orientation of the BSO crystal used in [12] is shown in Fig. 1. The light is incident on the plane (x, y) where the x axis is perpendicular to the plane (0,0,1) and the y axis is perpendicular to the plane (1,−1,0). In [12], measurements of the angular dependency of the PR signal were performed with a photodiode placed in the far field of the PR crystal and that scanned the PR signal along the k_x spatial frequencies where k_y ≈ 0. For this spatial frequency range, the internal space charge field has only a component along the x axis (i.e. the Ẽ_{sc_x}(k_t) component) and two birefringence axes appear, one parallel to the laboratory x axis, and one parallel to the y axis. The refractive indices found by the light polarized along these axes and propagating in the directions such that k_y = 0, are n₀ + Δn_x and n₀ + Δn_y where n₀ is the linear refractive index of the BSO crystal and [13]

{\begin{array}{l} Δ n_{x} (k_{x}) = & 0 \\ Δ n_{y} (k_{x}) = & - \frac{1}{2} r_{41} n_{0}^{3} {\tilde{E}}_{{sc}_{x}} (k_{x}) . \end{array}

As scattered light is strongly depolarized, we emphasize that only the component of the speckle field polarized along the y axis reads the modulation of the refractive index Δn_y induced by the component Ẽ_{sc_x} of the internal space charge field. These results suggest that the x component of the speckle field provides no contribution to the PR signal. Then, in the experimental setup, the use of a polarizer oriented along the y axis at the output of the BSO crystal could improve the signal to noise ratio of the PR signal. For the y axis polarization considered in this paper, BSO exhibits both photoinduced birefringence and large optical activity at the green wavelength. In such conditions wave progation of an incident depolarized speckle beam is complex and its detailed analysis may be out of the scope of this paper. When taking account of the optical activity for linearly polarized plane waves, it leads to a reduction of the BSO holographic efficiency [16].

Because the modulation of the refractive index is weak (10⁻⁵ to 10⁻⁶), the PR phase hologram induced by the speckle and read by the y polarized light can be expressed by :

H_{y} (x, y) = exp (- i \frac{2 π Δ n_{y} L}{λ}) ≃ 1 - i \frac{2 π Δ n_{y} L}{λ}

where L is the thickness of the hologram. Then, the y component of the output speckle field is given by:

S_{r y}^{out} (x, y) = S_{r y}^{in} (x, y) H_{y} (x, y) .

Using the Fourier transform expression of the space charge field component given by the Eq. (3) and the expression of the phase function of the hologram given by the Eq. (5), the (k_x, k_y = 0) spatial frequency spectrum of the output speckle is given by:

\begin{array}{l} {\tilde{S}}_{r y}^{out} (k_{x}) & = {\tilde{S}}_{r y}^{in} (k_{x}) \otimes (δ (k_{x}) + \frac{β k_{x}}{1 + \frac{k_{x}^{2}}{k_{0}^{2}}} {\tilde{I}}_{r} (k_{x})) \\ = {\tilde{S}}_{r y}^{in} (k_{x}) + {\tilde{S}}_{r y}^{in} (k_{x}) \otimes \frac{β k_{x}}{1 + \frac{k_{x}^{2}}{k_{0}^{2}}} {\tilde{I}}_{r} (k_{x}), \end{array}

where ⊗ denotes the convolution product and δ(k_x) is the Dirac function centered at the zero spatial frequency and β a new constant including all the constants of the previous equations. We can see that the amplitude of the output spectrum of the speckle diffracted by the photorefractive hologram is made of the incident speckle field (i.e. the zero order of diffraction), and a second term corresponding to the diffraction of the incident speckle by its own hologram. From an algebric point of view, this term is of the same sign as k_x and explains the origin of the reduction or the enhancement of the PR signal with respect to the position of the detector along the k_x axis. We stress out that this term is null in the mean direction of the incident light (i.e. the zero spatial frequency) so that no PR signal can be recorded in this direction. This property can be physically explained if considering the speckle field within the crystal as a superposition of plane waves whose incidences are included within the aperture of the lenses that shine the cristal : a multi two-wave mixing process is performed along the mean propagation axis z, leading by symmetry to a photorefractive gain which can be positive or negative, depending on the sign of the angle of the plane wave with respect to z. As a consequence, the superposition of all these contributions cancels the effect along z.

Because of the small phase modulation induced by the hologram, the diffracted speckle is very similar to the reference speckle and the contribution of the PR signal in the output speckle is very weak. Then, the PR signal is obtained by subtracting the background level of light such that:

I_{PR} (k_{x}) = {| {\tilde{S}}_{r y}^{out} (k_{x}) |}^{2} - {| {\tilde{S}}_{r y}^{in} (k_{x}) |}^{2} .

In order to study the diffraction efficiency of the PR hologram, numerical simulations are performed with different sizes of the speckle grains varying in the range [0.1μm – 3μm] and speckles with a mean intensity kept constant. Figures 4(a) and 4(b) show, for different sizes of the speckle grains, the variation of the modulus of the x component of the space charge field and of the PR signal with respect to the normalized spatial frequency $\frac{k_{x}}{k_{0}}$ . In agreement with Eq. (3), the amplitude and the position of the maximum of the space charge field depend on the size of speckle grains. The best efficiency is obtained at the spatial frequencies k_x = ±k₀ when the size of the speckle grains corresponds approximately to the Debye screening length ( $σ_{s} \approx \frac{1}{k_{0}} = 0.21 μ m$ ). The experimental results reported in Fig. 2 of [12] show the angular dependency of the PR signal measured with a large area photodiode of 3.6×3.6 mm² placed at a distance of 25mm of the PR crystal (in this case the optical system L₃L₄ in the experimental setup is removed). In order to fit with these experimental conditions, the PR signal curves presented in Fig. 4(b) are the result of the convolution product between the PR signal calculated with Eq. (8) and a square function, with a width Δk_D = 1.7 μm⁻¹, that simulates the 0.14 rad angular aperture of the photodiode. Moreover, to smooth out the intensity fluctuations of the PR signal due to the small number of speckle grains in the area of the simulated photodiode, we average the calculated PR signal over more than 10 numerical realizations to obtain these curves. Finally, the PR signal is plotted simultaneously with respect to the normalized spatial frequency $\frac{k_{x}}{k_{0}}$ and with $θ = asin (\frac{k_{x} λ}{2 π})$ , the external angle (outside the crystal). In Fig. 4(b), the shape of the PR signal is in a good agreement with the results presented in [12] where amplitude of the PR signal was maximized for a position of the detector corresponding to an external angle θ = ±17° when the size of speckle grains was on the order of l_D. These results also suggest that experimental parameters used in [12] were not optimized. Indeed, Fig. 4(b) shows that the PR signal is maximized for larger speckle grains (0.4μm ≤ σ_s ≤ 0.9μm) while the amplitude of the space charge field is maximized when the size of the speckle grains fits with the Debye screening length. Secondly, the efficiency of the PR signal is optimized when the detection is performed in a direction corresponding to a spatial frequency smaller than k₀. Physically, this result differs from the two-wave mixing experiments where the photorefractive gain is optimized when the grating wave number K equals k₀ [16]. Indeed, in this case the terms related to the speckle in Eq. (7) will be replaced by Dirac functions δ (k_x ± K) that leads to a diffraction efficiency which is governed only by the impulse response of the PR crystal $\frac{β k_{x}}{1 + \frac{k_{x}^{2}}{k_{0}^{2}}}$ and maximized when K = k₀. With speckles, convolution product in Eq. (7) implies that diffraction efficiency depends on the impulse response of the PR crystal as well as the spatial frequency width of the speckle $σ_{k} = \frac{1}{σ_{s}}$ . Finally, finite size of the detector is also a factor to be taken into account to optimize measurement of the PR signal. Figure 4c shows the amplitude of the PR signal convolved with a square function with a width Δk_D and calculated as a function of the normalized spatial frequency $\frac{k_{x}}{k_{0}}$ and the normalized half width of the gaussian envelop of the speckle spectrum $\frac{σ_{k}}{k_{0}}$ . The black dotted curve depicts the position of the detector that gives the higher PR signal as a fonction of σ_k. Vertical and horizontal black dotted lines make visible the couple of parameters (σ_k = 0.20k₀, k_x = 0.35k₀) optimizing the PR signal. We can stress out that this optimized couple of parameters depends on the angular aperture of the detector.

Fig. 4 Variation of the modulus of the x component of the space charge field (a) and of the PR signal (b) as a function of the normalized spatial frequency $\frac{k_{x}}{k_{0}}$ (k_y = 0) and of the mean size σ_s of the speckle grains. (c) Normalized PR signal as a function of the normalized spatial frequency $\frac{k_{x}}{k_{0}}$ and the normalized half width of the gaussian envelop of the speckle spectrum $\frac{σ_{k}}{k_{0}}$ . The black dotted curve depicts the optimized position of the detector that gives the higher PR signal as a fonction of σ_k.

Download Full Size | PDF

Now, in order to modelize the ultrasound modulated optical tomography by self-referenced photorefractive holography reported in [12], let us consider the case where an incident speckle tagged by a short US pulse, noted $S_{t}^{in}$ , is diffracted by the hologram induced by the reference speckle S_r. Because the interaction between the short US pulse with the light propagating in the scattering medium leads to changes of the speckle pattern, the tagged speckle is transiently not perfectly matched with the complex phase hologram recorded in the whole volume of the BSO crystal. As a consequence, the level of the PR signal varies during the time of the US pulse.

In our simulations, the tagged speckle is modelised by the coherent superposition of the reference speckle S_r with a ”noise” noted S_n (i.e. a fully decorrelated speckle) so that the amplitude of the tagged speckle may be written in the (x, y) domain as:

S_{t}^{in} = \sqrt{γ} S_{n} + \sqrt{1 - γ} S_{r},

where γ represents the ratio of the incident light interacting with the US pulse in the scatering medium. Then, the intensity of the tagged speckle is given by:

I_{t}^{in} = γ I_{n} + (1 - γ) I_{r} + 2 \sqrt{γ (1 - γ)} ℜ 𝔢 {S_{n} S_{r}^{*}} .

Because the noise speckle S_n is fully decorrelated with the reference speckle S_r, the mean value of the interference term between the two speckles is null. Figure 5 shows the typical result of the interference between two speckles fully decorrelated. Therefore, this term can be neglected when γ ≪ 1 so that the intensity of the tagged speckle is assumed to be the incoherent sum of the intensities of the two speckles. According to this assumption, the detected intensity of the tagged speckle that reads the PR hologram after US interaction can be written, in the spatial frequency domain, as :

{| {\tilde{S}}_{t}^{out} |}^{2} ≃ γ {| {\tilde{S}}_{n}^{out} |}^{2} + (1 - γ) {| {\tilde{S}}_{r}^{out} |}^{2} .

Experimentally, the position of the photodiode is adjusted along the k_x dimension to maximize the PR signal. Then, the measured AO signal is directly the intensity of the spectrum of the tagged speckle at this spatial frequency and the DC background corresponds to the intensity of the untagged speckle measured at the same spatial frequency. Hence, we will define the relative AO signal as :

R_{AO} = \frac{\iint_{Δ k_{D}} {| {\tilde{S}}_{t}^{out} |}^{2} d^{2} k - \iint_{Δ k_{D}} {| {\tilde{S}}_{r}^{out} |}^{2} d^{2} k}{\iint_{Δ k_{D}} {| {\tilde{S}}_{r}^{out} |}^{2} d^{2} k}

where

\iint_{Δ k_{D}}

denotes the integration of the intensity of the spectrum of the diffracted speckle field over the spatial frequency bandwidth Δk_D of the photodiode. For numerical calculation, the mean radius of the speckle grains in the (x, y) is adjusted to 0.21 μm (i.e. the Debye screening length). Figure 6(a) represents the typical shape of the intensity of the spatial frequency spectrum of a tagged speckle and the two yellow squares symbolise the optimized position and the area of the detector centered on the spatial frequencies

k_{x} = \pm \frac{k_{0}}{2}

(see Fig. 4(b)). For these two positions of the photodiode and using Eq. (12), we calculate the variation of the relative AO signals as a function of the level γ of the ”noise” in the tagged speckle. Assuming that the scattered light totally filled the volume of the scattering medium, the maximum value taken by γ is determined very approximately to 10⁻³ by calculating the ratio between the focusing volume of the US pulse (∼ 5mm³) and the total volume of the scattering medium (< 40³ mm³). Figure 6(b) gives the variation of the relative AO signals, respectively noted R_AO+ and R_AO− when the detector is centered on

+ \frac{k_{0}}{2}

or

- \frac{k_{0}}{2}

, as a function of the level of noise in the tagged speckle. Because of the small number of speckle grains in the area of the simulated detector, the curves correspond to the AO signals averaged over few tens of numerical realizations. In agreement with the shape of the PR signal given by Fig. 4(b) and Eq. (11), the relative AO signal measured around

+ \frac{k_{0}}{2}

and

- \frac{k_{0}}{2}

respectively decreases and increases linearly, when the level of noise increases.

Fig. 5 (a) and (b) intensities in the (x, y) domain of the reference speckle and the fully decorrelated noise speckle. (c) Real part of the interference term between the reference and the noise speckles.

Download Full Size | PDF

Fig. 6 (a) Intensity of the spatial frequency spectrum of a tagged speckle. The two yellow squares represent the detector area in the (k_x, k_y) domain when it is centered respectively on the spatial frequencies $k_{x} = \pm \frac{k_{0}}{2}$ . (b) Amplitude of the AO signals as a function of the level γ of noise in the tagged speckle for the two positions of the detector.

Download Full Size | PDF

All these results give a satisfactory quantitative interpretation of the experimental results reported in [12]. However, the developed model mainly outlines the optimal values of the speckle grain size and of the detector position and, for the moment, does not permit to quantify the holographic efficiency and the AO signal dynamic range. In future work, several experimental parameters, such as fraction of tagged photons, beams polarization, light collection efficiency on the detector need to be carefully considered in order to refine amplitude consistency between the model and experiment.

3. Conclusion

The self referenced wavefront acousto-optic imaging technique has been analysed in the condition where a single speckle field is projected on a PR crystal. In these experiments, no reference beam is applied and the incident speckle creates a complex space charge field inducing a phase volume hologram in the BSO. The model leads to determine the parameters controling the efficiency of the hologram. We show that the PR signal efficiency is maximized for a spatial frequency that depends on the crystal Debye screening wave number, the width of the spatial frequency spectrum of the speckle but also on the angular aperture of the detector. It results that there is an optimum position of the photodetector and an optimum size of the speckle grains to detect the speckle pattern in presence of the tagged photons in an acousto-optic imaging experiment. The numerical simulations agree well with the experimental results previously reported. The method is robust and easy to implement in comparison with conventional holography. It can be applied with other types of nonlinear media in particular operating in the near IR region for acousto-optic imaging in highly scattering biological tissues and in turbid media and vibrometry.

Acknowledgments

This work was supported by the LABEX WIFI (Laboratory of Excellence within the French Program Investment for the Future) under references ANR-10-LABX-24 and ANR-10-IDEX-0001-02 PSL*, the ANR project under the name ICLM (Coherent Imaging of Multi-scattered Light, ANR-2011-BS0417) and the Conseil Régional de Franche-Comté. Computations have been performed on the supercomputer facilities of the Mésocentre de calcul de Franche–Comté.

References and links

1. S. R. Arridge, “Optical tomography in medical imaging,” Inverse Probl. 15, R41–R93 (1999). [CrossRef]

2. A. Bratchenia, R. Molenaar, T. G. van Leeuwen, and R. P. H. Kooyman, “Acousto-optic-assisted diffuse optical tomography,” Opt. Lett. 36, 1539–1541 (2011). [CrossRef] [PubMed]

3. L. Wang, S. L. Jacques, and X. Zhao, “Continuous-wave ultrasonic modulation of scattered laser light to image objects in turbid media,” Opt. Lett. 20, 629–631 (1995). [CrossRef] [PubMed]

4. G. Rousseau, A. Blouin, and J.-P. Monchalin, “Ultrasound modulated optical imaging using a high-power pulsed laser and a double-pass confocal Fabry-Perot interferometer,” Opt. Lett. 34, 3445–3447 (2009). [CrossRef] [PubMed]

5. X. Xu, S.-R. Kothapalli, H. Liu, and L. V. Wang, “Spectral hole burning for ultrasound-modulated optical tomography of thick tissue,” J. Biomed. Opt. 15, 066018 (2010). [CrossRef]

6. M. Gross, P. Goy, and M. Al-Koussa, “Shot-noise detection of ultrasound-tagged photons in ultrasound-modulated optical imaging,” Opt. Lett. 28, 2482–2484 (2003). [CrossRef] [PubMed]

7. E. Benoit à la Guillaume, S. Farahi, E. Bossy, M. Gross, and F. Ramaz, “Acousto-optical coherence tomography with a digital holographic detection scheme,” Opt. Lett. 37, 3216–3218 (2012). [CrossRef] [PubMed]

8. A. Lev and B. Sfez, “In vivo demonstration of the ultrasound-modulated light technique,” J. Opt. Soc. Am. A 20, 2347–2354 (2003). [CrossRef]

9. S. Farahi, E. Benoit à la Guillaume, A. A. Grabar, J.-P. Huignard, and F. Ramaz, “Time resolved three-dimensional acousto-optic imaging of thick scattering media,” Opt. Lett. 37, 2754–2756 (2012). [CrossRef] [PubMed]

10. M. Lesaffre, F. Jean, F. Ramaz, A. C. Boccara, P. Delaye, and G. Roosen, “In situ monitoring of the photorefractive response time in a self-adaptive holography setup developed for acousto-optic imaging,” Opt. Express 15, 1030–1042 (2007). [CrossRef] [PubMed]

11. A. A. Kamshilin, T. Jaaskelainen, and Y. N. Kulchin, “Adaptive correlation filter for stabilization of interference fiber optic sensors,” Appl. Phys. Lett. 73, 705–707 (1998). [CrossRef]

12. E. Benoit à la Guillaume, U. Bortolozzo, J. P. Huignard, S. Residori, and F. Ramaz, “Dynamic ultrasound modulated optical tomography by self-referenced photorefractive holography,” Opt. Lett. 38, 287–289 (2013). [CrossRef] [PubMed]

13. P. Yeh, Introduction to Photorefractive Nonlinear Optics (John Wiley, 1993).

14. O. Daniel, J.-M. C. Jonathan, and G. Roosen, “Photorefractrive effect in the Fourier plane,” Opt. Mater. 4, 294–298 (1995). [CrossRef]

15. O. Daniel, A. Stelmach, J.-M. C. Jonathan, and G. Roosen, “Whole-beam method analysis of photorefractive effect in correlators,” Opt. Commun. 113, 559–567 (1995). [CrossRef]

16. J. A. Gómez, H. Lorduy Gómez, and Á. Salazar, “Novel procedure for the simultaneous determination of the Debye length and electro-optic coefficient for an optically active photorefractive Bi₁₂SiO₂₀ crystal,” Opt. Commun. 284, 460–466 (2011). [CrossRef]

17. J. W. Goodman, Laser Speckle and Related Phenomena (Springer, 1975).

Modelization and optimized speckle detection scheme in photorefractive self-referenced acousto-optic imaging

Abstract

1. Introduction

2. Theoretical model and numerical simulations

3. Conclusion

Acknowledgments

References and links

Cited By

Figures (6)

Equations (12)

Optics Express