Noninvasive diagnostic imaging technologies to detect cancer and other diseases come in a wide variety. Most of these technologies like MRI and x-ray lack optical contrast that can be beneficial for detection of small lesions without the use of contrast agents or ionizing radiation.1
Optical techniques such as optical coherence tomography (OCT) and diffuse optical tomography (DOT) show good contrast, but the measuring depth and resolution, respectively, are limited by the strong scattering of tissue. Measuring with high resolution at penetration depths larger than 1 mm is challenging due to the high optical scattering that induces a strong attenuation of ballistic photons. These optical techniques are noninvasive and use nonionizing radiation in contrast to x-ray. The absorption coefficient of light is wavelength dependent, making spectroscopy possible, which can be used to determine blood oxygenation levels.
Ultrasound (US) in the few MHz ranges has a scattering coefficient two to three orders of magnitude2 less than light, which comparatively allows for superior penetration depth while retaining spatial resolution. There is, however, a drawback in which US lacks the benefits of optical contrast.
Two hybrid forms of imaging techniques that use the interplay of sound and light are emerging: acousto-optics (AO) and photoacoustics (PA). These techniques combine the high resolution of US with the strong contrast found in purely optical techniques while remaining noninvasive in nature.
The principle of PA relies on a light pulse irradiating an object of interest. Some parts of the object will absorb this light and will expand thermoelastically. The resulting expansion causes an US wave that is detected by a broadband US probe.3 The resolution is limited by the light pulse duration, the acoustic tissue properties, and the frequency transfer function of the US detection. The higher the detected frequencies the higher the resolution will be. The longer the light pulse duration, the lower the frequencies of the PA signal and the lower the resolution. The acoustic attenuation typically increases with the acoustic frequency, therefore the resolution deteriorates with increasing depth. Typical resolutions are micrometers for PA microscopy with a measuring depth of millimeters4 and up to millimeters for PA mammography5 on a measuring depth of centimeters. While the resolution is achieved by the ultrasound, the contrast is provided by optical absorption; the amount of detected pressure depends on the optical absorption inside the object under investigation.
A second technique that combines ultrasound and light is AO. The possibility of tagging light was investigated by Marks et al.6 Leutz and Maret,7 Wang et al.8 and by a patent from Dolfi and Micheron.9 The sample is illuminated with coherent light, and hence a speckle pattern is formed due to complete randomization of the phases of the electric field of the light escaping from the medium.
To measure only photons that traveled through a volume of interest, ultrasound (US) is employed to encapsulate this bounded region from the characteristic focal zone of the transducer. The photons passing through the US focus will be frequency shifted or, in other words, phase modulated and thus tagged by the US. A variety of tagging strategies have been developed, along with different detection schemes to find the tagged photons in the sea of untagged photons. The tagged photons detected after exiting the medium provide information on the optical properties present in the vicinity of the volume of interest. When the tagging zone is placed in a tissue volume with high optical absorption compared with the rest of the object, less tagged photons will be detected. This in contrast with when the ultrasound focus is placed in the neighboring tissue. In this review we describe the progress made in AO imaging since the previous review by Wang10 in 2003: theoretical models, the detection methods to detect the tagged photons, and applications of AO imaging. For further reading, are view by Elson et al.11 is also focusing on the detection methods, scanning, and applications. They give, however, no equations from the models in the modeling section.
In recent years, considerable improvements in AO have been realized. Some groups developed a method to use AO in a reflection mode setup.1213.14.–15 Others16,17 investigated the decay of the modulation of the signal as function of the imaging depth and enhancement of AO with micro bubbles, which can be important in the development of new applications.
Acousto-Optic Tagging Mechanisms
Light can be modulated in a scattering medium using ultrasound by sixdifferent mechanisms. The relative influence of each of these mechanisms depends on the optical properties of the sample and light source and the acoustic properties of sample and transducer.
The first two mechanisms (see Fig. 1) require light to have a sufficient coherence length. The first mechanism relies on the modulation of the optical path length. The scatterers inside the medium are assumed to oscillate with the US frequency at the US focus due to the pressure wave. The optical path length varies with the distance between the scatterers, and therefore the speckle pattern will vary with the applied US frequency.7 This mechanism is thus equivalent to Doppler scattering and relies on coherent light with a long coherence length to obtain a clear speckle pattern. The mean optical path length is, depending on the optical properties, in the order of 10 to 30 cm for an object with a thickness of 3 cm.18,19 The coherence length must be at least in the same order of magnitude. Leutz and Maret7 have developed a theory that partially modeled this mechanism. The theory is only valid when the mean free path is much greater than the light wavelength, and the particle displacement must be much smaller than the light wavelength. Unfortunately, the second mechanism is dominant in this case.
The second mechanism is based on refractive index variations due to the acoustic pressure wave. As a result the optical path lengths, and therefore the phases, are modulated. Consequently, the intensity of the resulting speckle pattern is also modulated.10
The relative influence of particle vibration and refractive index changes on the induced phase change depends on the optic and acoustic properties of the medium and the used waves.2021.–22 Properties such as scattering coefficient and the wavelengths of the ultrasound and probing light change the relative influence of the first and second mechanisms. The strengths of both mechanisms are comparable up to a critical point from where the acoustic wave vector becomes larger than a critical fraction of the mean free path of the photons and the second mechanism becomes dominant.20
The otherfourmechanisms (see Fig. 2) are incoherent phenomena. As a result of US associated local pressure variations, the medium is continuously locally compressed and decompressed. Local optical properties such as absorption (mechanism 3), position, and scattering cross-section of scattering particles (4 and 5) and the refractive index6 of the medium will oscillate with the US frequency due to variations in density.10 The four mechanisms each result in variation of the fluence distribution inside the object. In many situations the associated variation in signal strength is too low to be measured.10 However, Krishnan23 proposed a theory that fluorescent light is modulated by the gradient of the refractive index by an acoustic lens effect (mechanism 6). Kobayashi et al.24 showed the feasibility of detecting this modulated intensity in fluorescence.
To understand AO modulation, and to explore the relative importance of the above mechanisms for modulating the light, different theoretical models have been developed. Most models use the more classical approach of phase changes induced by ultrasound, as described by Secs. 3.2 and 3.3. A different approach used by Mahan et al.25 is more from a quantum mechanics perspective. Here the fraction of tagged and untagged photons is calculated (see Sec. 3.1).
Fraction of Tagged Light by Brillouin Scattering
Mahan et al.25 used the mechanism of refractive index changes due to the ultrasound to model the tagged light. The periodic refractive index changes causes Brillouin scattering of the photon stream. A small fraction of these photons will scatter inelastic and are frequency shifted with the US frequency. The intensity of the tagged photon stream is a fraction of the untagged photon stream. They start by deriving an expression for the signal strength of the untagged photons. Step two is calculating the fraction for different geometries. For a slab with thickness
Theoretical Predictions of Observable Optical Quantities During US Modulation
The ultrasound modulation of light has been modeled for its consequence for two observable aspects of speckle patterns: temporal fluctuations and reduction of speckle contrast, which are coupled to one another. The modulation of the speckle pattern formed at the detector is often calculated by using the temporal autocorrelation function of the electric field. The intensity of a speckle consists of a large DC component and a smaller component, the one in which we are interested, that varies in time with the US frequency. Leutz et al.7 and others10,20,21,26,27 expressed this autocorrelation function as:10 Point-like scatterers are undergoing a collective motion due to the ultrasound and Brownian motion. The autocorrelation function length can be split in one part describing the exponential decay due to Brownian motion and a second part that describes the autocorrelation as a function of the periodic ultrasound wave. When both effects are assumed independent from the other then the properties of both factors can be explored separately.
The scatterers in the medium may tend to diffuse due to the Brownian motion. The single particle relaxation time is given by where is the particle diffusion constant and the wave vector of the light. The autocorrelation of the electric field with a single path length will decay exponentially in time according to:7
This decay time limits the time in which an accurate AO measurement can be performed in most situations. The Brownian motion can be neglected for time scales in the order of the cycle time of the ultrasound.28 Wang20,21 investigated the contribution of phase modulation by refractive index changes and scatterer displacement induced by the ultrasound. He concluded that both effects have a similar contribution up to a critical point. When the mean free paths of the photons between two scattering events become small compared with the acoustic wavelength, the second mechanism becomes increasingly more important. This was verified with a Monte Carlo simulation.20
Sakadžić et al.26,2930.–31 and Wang20 developed an analytical solution for this autocorrelation function Eq. (2), mostly focusing on the first two mechanisms. For a single pathlength the contribution to the autocorrelation function is given by:2010,26 and the assumption of weak modulation, the modulation depth, defined as ratio of the intensity at the signal frequency and the unmodulated intensity,10 is found to be proportional to the acoustic amplitude squared.10 In general the maximum variation of the autocorrelation function will increase with increasing acoustic power and decrease with increasing US frequency while keeping the US power constant. Increasing the optical absorption coefficients will lead to a smaller variation according to:22,26
The intensity and variation of the speckle contrast, which is related to the autocorrelation function, are derived by Zemp et al.27 The statistical properties of this speckle pattern are derived from the autocorrelation function. The first order statistics describe the time average intensity and is approximately equal to . The second order statistics describe the variation in speckle contrast. The variation is approximatelyRef. 27. All these theories assume homogeneous acoustic waves, except for Sakadžić et al.2930.–31 in which localized acoustic waves were used.
Monte Carlo Simulations
Besides the analytical solutions Monte Carlo (MC) simulations are used to gain more insight into the principles of AO. These models simulate photon distributions inside the medium and include scattering and absorption of the light. One or more of the tagging mechanisms can be implemented, and the refractive index and scatterer motion can be modulated. The results of measurements and analytical models are often compared and tested with MC simulation results.3233.–34 Wang20 and others33 developed and modified MC models for AO applications. Especially the model of Wang et al. is widely used in Biomed.Op. i.e.,35,36 and is supported by an analytical model.21,26,29,30
Leung et al.37 compared the speed of this model implemented on both central processing unit (CPU) and graphics processing unit (GPU, Nvidia GeForce 9800) with Compute Unified Device Architecture (CUDA). They found an increase in speed, depending on simulation settings and hardware, of a factor 72 for the CUDA implementation.
Detection of Tagged Photons
A detection technique must be able to distinguish the tagged or modulated photons from the untagged or unmodulated photons to obtain information from the US focus. Longer integration times will lead to higher signal-to-noise ratios. However, speckle decorrelation alters the speckle pattern and destroys the correlation of this pattern between the start and end of the detection time. Sources of speckle decorrelation are small changes in the optical paths inside the object, i.e., Brownian motion of scatterers and blood flow or muscular activity in living tissues. This decorrelation diminishes the detectability of tagged photons. Developed detection techniques can be divided in time- and frequency domain techniques.
Wang et al.8 measured light modulated at a US frequency of 1 MHz with a photomultiplier. The intensity of a single speckle will oscillate with the US frequency of 1 MHz, and this frequency is observed by recording the light flux with the photomultiplier tube. They used a 1.1-mm aperture at 5-cm distance from the object and a photomultiplier tube at a distance of 10 cm from this aperture. The advantage of a photomultiplier tube is the speed and sensitivity, so the modulation of a speckle can be detected real time at the US frequency and well within the speckle decorrelation time. The disadvantage is the small detection surface, which in turn reduces the number of tagged photons that can be detected, leading to small signal-to-noise ratio due to the influence of shot noise. In a nutshell: increasing the number of speckle grains on the same single detector increases the random modulated signal standard deviation (signal of interest) as and the average power impinging the detector by . If the detection is shot noise limited (that is rather unlikely to happen) the noise is proportional to , and the signal-to-noise ratio is not increased by increasing the number of speckle grains on the detector. If the noise is linked to the laser fluctuations it is proportional to and the signal-to-noise ratio will decrease when increasing the number of speckle grains on the detector.
A parallel speckle detection technique was developed by Leveque et al.38 Instead of a single detector element, a CCD camera is used in which the illumination pulses and the ultrasound are timed with a fixed phase delay. The system was designed to make speckle size comparable to the pixel size. The SNR is improved by averaging the signal from all CCD pixels, making it a larger detector surface. The major disadvantage of a CCD camera as a detector is the reduced frame rate, hence the intensity of a speckle cannot be followed in time at the US frequency. However, by recording the light intensity for different phases of the ultrasound, the modulated intensity and the unmodulated intensity can be estimated. This can be achieved by varying the delay between US and the start of the exposure of the camera. In the case of four-phase measurement the exposure time of the camera is up to one quarter of the US period time.
The signal for each phase from Fig. 3 is given by:39 showed similar relations for two and three phases.
Parallel speckle detection is an efficient technique for detecting modulated light. However, this detection method is sensitive tospeckle decorrelation, which lowers the SNR by lowering the signal and increasing the noise.18 Li et al.41 showed a laser speckle contrast detection scheme that is less sensitive for speckle decorrelation when short integration times are used. The speckle contrast is defined as the standard deviation of the intensity in the pattern divided by the average intensity of this pattern.42 Kothapalli et al.43 tested their Monte Carlo implementation and measured with a speckle contrast setup; they further investigated the linear relation between speckle contrast and the local scattering coefficient at the US focus and found a good agreement between simulation and experiment.
Instead of following the intensity of the speckles in time domain, it is possible to filter out the modulated light with interferometers. The Fabry-Perot interferometer and Mach-Zehnder based techniques are used to detect, or select, only the modulated light at the US frequency. By optically removing noninteresting frequencies, the SNR should be increased.
Leutz et al.7 used a Fabry-Perot interferometer in combination with a photomultiplier tube to detect photons with a frequency shift of 2.17 and 27.3 MHz from a laser source with a wavelength of 514.5 nm. To resolve this relative small frequency difference, they used a FP etalon (mirror reflectivity: 99.3%) with a mirror separation distance of 15 cm, obtaining a resolution of 12 MHz. The disadvantage of this technique is the great loss of signal photons due to the small etendue (the product of the solid angle and area of an aperture) when pinholes are used to select a parallel light beam with a certain frequency. However long-cavity44 confocal Fabry-Perot interferometers are used, which have a high etendue. Double-pass confocal45 Fabry-Perot interferometers are also used to separate the faint spectral line of the tagged photons from the strong untagged spectral line.
The heterodyne parallel speckle detection method adds a reference beam or local oscillator (LO) beam to the standard parallel speckle detection scheme,4647.–48 making it a Mach-Zehnder interferometer, as depicted in Fig 4. This technique is shot-noise limited when using a large heterodyne reference intensity. When the signal beam is in phase with the reference beam the intensity on the detector is given by:49
The light in the sample arm is modulated as it propagates through the sample, and the reference arm is modulated with two acousto-optic frequency shifters (Bragg cells) or acousto-optic modulators. The speckle pattern from the sample interferes with the reference beam on the CCD camera, which is placed under a small tilt angle . The reference light with the modulated part of the sample generates a static speckle pattern, which can be detected by the CCD camera. Interference of the nonmodulated sample light with the reference light leads to a dynamic speckle pattern, which in the CCD image will be smeared out. The spatial frequency domain (-space) of the recorded interference pattern reveals the level of the shot noise, the speckle decorrelation noise, and the tagged photons in different regions, due to spatial filtering associated with tilt angle . The sensitivity of this technique is limited by the shot noise of the reference light46,47 and setup limitations such as dynamic range and sensitivity of the camera.From the four frames of the CCD array, representing the four phases of the reference beam, a four-phase complex signal can be calculated. This complex signal in -space shows interesting properties. The heterodyne signal forms a narrow band of spatial frequencies at a position determined by tilt angle (region A in Fig 5), the decorrelation noise is relatively slow, and its fringes consists of the lower spatial frequencies (region B and C). The shot noise is represented by region D.
For example, the laser light is modulated at 85 MHz up and 80 MHz down with the use of two acousto-optic modulators. For practical reasons, two modulators are used instead of one at 5 MHz; that is the difference of both used modulators. The laser light inside the sample is modulated in the focus of a US transducer at 5 MHz.
This is known as single-phase detection; four-phase detection can be achieved by tuning one of the AO modulators slightly higher, 85.0000075 instead of 85 MHz. This way each next frame of the CCD array receives a reference beam with a phase shift of compared with the current frame assuming a frame rate of 30 frames per second.
The advantage of using multiple-phase detection is the higher SNR because the DC component is canceled out. All frames of the one to four phases should be taken within the speckle decorrelation time, making one or two phase imaging the better choice where this decorrelation time is short. Atlan et al.48 give an equation to calculate the complex signal for an arbitrary number of phases.
Another technique uses the same heterodyne setup with the addition of a photo refractive crystal (PRC) in front of the detector. The principle is to store intensity and phase information of the speckle pattern inside a hologram and readout this hologram with a CCD. The advantage of this self-adapting wave-front holography is a larger etendue than most CCDsand possibly a larger SNR. Bloningen et al.50,51 investigated the expected AO signal and phase shifts of the photons for photo-refractive-based detection by utilizing Monte Carlo simulations. They show that the average photon phase shift contains a DC phase shift, which is dominant over the usually used AC phase shift. This offers a possibility to compliment the AC data.
The crystal consists of a photorefractive material and thus has a refractive index that depends on the light intensity. The two beams, signal and reference, illuminate the crystal. The two interfering beams construct a three-dimensional (3-D) intensity profile in the crystal thus writing, using the photo refractive principle, a 3-D grating or hologram in the crystal within the response time. Typical response times are 100 µs for a multiple quantum well crystal (850-nm excitation wavelength),520.3 to 100 ms for GaAs (1064 nm),47,53188.8.131.52.–58 less than 10 ms for SPS:Te (790 nm)59 and 150 ms for .51,60,61 Changes in the speckle pattern slower than this response time are not recorded, making the detection insensitive for slow speckle decorrelation. The photo-refractive effect selects only optical frequencies from the signal beam, which are close to the reference-beam frequency.
During readout of the hologram, the reference beam is diffracted by the holographic grating. This way the wavefront of the reference beam is shaped in the signal wavefront, but now with a larger intensity. This constructed speckle pattern can be detected by a detector such as a CCD. Chi et al.62 describes this two wave mixing process in more detail. Gross et al.53 give a detailed description of detecting tagged light with a PRC.
This self-adapting wavefront holography is a promising technique, but the typical response time is too slow for in vivo use where the speckle decorrelation times () are in general much shorter than the crystal response time ().57
Spectral hole burning (SHB) is another crystal-based technique developed by Li et al.57,63 where the PRC is replaced by a SHB crystal. This technique has also a large etendue and is insensitive for speckle decorrelation. The crystal acts as a narrow optic bandpass filter where most optical frequencies are absorbed and only one specific frequency band is transmitted. The reference beam burns a spectral hole in the crystal at the same frequency as the tagged photons. The crystal is an inhomogeneously broadened optical absorber that can be modeled as a two-level system. The spectral line width is typically sub-MHz when cryogenically cooled and depends on the intensity of the reference beam. During the burning process only the crystal ions inside the crystal that are sensitive for this specific optical frequency are excited from the ground state. This makes the crystal transparent for this optical frequency, the nontagged photons are absorbed by the crystal.
After burning the spectral hole, the reference beam is switched off and the signal beam is switched on, the tagged photons can be detected by a photo detector such as a CCD. The cryogenic cooling () is the main disadvantage of this technique.
The lateral resolution of acousto-optic imaging is obtained by focusing the ultrasound, and its size is comparable with the width of the US focus when the weight of the unfocused part does not play a major role. The axial resolution is worse than the lateral resolution due to the elongated focus of the ultrasound transducer (Rayleigh length).
The axial resolution can be improved in several ways by changing frequency, amplitude, or phase in time in order to get a resolution equivalent to the pulse-echo in US tomography.
The first method that was implemented is the frequency sweep technique.64,65 In this method the ultrasound transducer generates a short chirped signal. The light along the acoustic axis is modulated with a different frequency for each depth. By using the Fourier transform on the detected signal, information in the axial direction is retrieved. Yao et al.65 combined this chirp technique with a parallel detection scheme. They also showed that this technique does not rely purely on ballistic light by using an obliquely incident laser beam. Wang et al.64 showed an axial resolution of approximately 1 mm;however, this resolution depends on the chirp rate. In chicken-breast tissue they achieved 4-mm resolution.65
Lesaffre et al.58 showed a different approach by making the phase of the US and the optical illumination random. Both use the same data set of typical 512 random phases with each of these phases set to either 0 or . This set of random phases has a typical period time 1 ms. The phase of the optical illumination is modulated with a delay relative to the modulation of the US. At a distance of times the speed of sound along the acoustic axis, both the optical illumination and US phase are correlated. Therefore only tagged photons from this position are detected using the heterodyne detection scheme with a photo refractive GaAs crystal. The resolution in the axial direction depends on the speed of sound and the time it takes to change phase, which in this case is () roughly 3 mm. Gross et al. give more theoretical background66 in regarding this technique.
A third method to increase the axial resolution uses short US pulses.44,48,51,67 These pulses consist of a few US cycles and thus modulate the light only locally in the axial direction. The US pulse reaches the point of interest after time , where is the distance from the transducer and the point of interest and the speed of sound (SOS). At this time () the sample is illuminated with pulse duration , and the modulated photons are recorded. The minimally resolvable size in the axial direction depends linearly on the SOS and pulse duration, assuming a short enough light pulse. Achieving high resolution thus requires short US pulses. The smaller the pulse duration, the bigger the spread of acoustic energy in the frequency domain. Detection schemes, which are only sensitive for a single frequency, require a longer US burst and thus a smaller US bandwidth to optimize the SNR.67 Atlan et al.48 used an US burst of 3 µs and the axial resolution is 5 mm.
Li et al.68 investigated the influence of amplitude modulated US on the speckle contrast in a parallel detection setup. In this case the US amplitude was modulated with a 250-Hz-to-16-KHz cosine envelope. They show that amplitude modulated US can outperform continuous wave US in signal strength. Weng69 investigated the relation between US pulse duration and amplitude on the AO signal and concluded that the imaging depth can be increased when reducing pulse duration and increasing the US peak amplitude.
Above methods all describe a way to improve the axial resolution by shaping the acoustic wave. The lateral resolution, however, is achieved by focusing the US transducer on the volume of interest. A different approach is suggested by Kuchment70 where an unfocused US transducer is used. An unfocused transducer has spherical wavefronts in the far field. When the US transducer emits a small US wave train, a small circular shaped region is tagged. By repeating these measurements for different positions of the US transducer, a synthetic focus is constructed on the positions where these circular regions intersect. Another example of what could be seen as synthetic focusing is used by Li et al.71 Sakadžić et al.44 used tomography to achieve very high resolutions in the order of 100 µm.
Let us underline that whatever is the technique that led to reduce the resolution it has to involve signal variation shorter than the decorrelation time. This is obviously realized for short pluses but must also be done for other approaches such as the chirpsor the random sequences.
Applications of Acousto-Optic Imaging
Measuring Optical Properties
While it is still unknown what the measured quantity is in AO, its signal probably holds information on the local fluence rate. Lev and Sfez12 showed the possibility of measuring the local fluence rate using AO between two optodes. This local fluence rate depends on position of illumination and detection and the optical properties distribution in the object such as the absorption coefficient. A map of the estimation of the optical absorption is obtained by scanning the US focus through the object and estimate the absorption coefficient for each position from the modulated fluence rate.42 The AO signal is affected not only by optical absorption and scattering in the US focus but also in the rest of the object. Therefore quantitative estimation of the absorption coefficient requires the number of tagged and untagged photons and an algorithm that solves the inversion problem as in DOT.30,33,54,7273.–74
Gunadi and Leung75 investigated the sensitivity of AO for the application of spectroscopy. From the spatial distribution of AO sensitivity they derived a penetration depth for AO of 14.8 mm for an intralipid solution with a µs’ of and absorption of µa of .
The absorption coefficient is in most cases dependent of the wavelength of the used light source, making functional imaging possible. For instance, chromophores such as hemoglobin and oxy-hemoglobin have distinct absorption spectra. This creates the possibility of determining the oxygenation levels of blood in biological tissues.33
Acousto-Optic Modulated Fluorescence
Fluorescence is often used to label particles and cells and enables the study of biological processes on the molecular and cellular level. However, light scattering makes it difficult to determine the origin of the fluorescent light in vivo. Efforts similar to DOT have to be made to perform optical tomography with fluorescent light.76 Kobayashi et al.,24 Hall et al.,77 and Yuan et al.76,7879.80.–81 investigated techniques to modulate and thereby localize this fluorescent light in turbid media.
The fluorescent light is noncoherent, and modulation has to come from modulation of optical properties of the sample. The US field causes variation in the density and, consequently, induces a gradient of the refractive index based upon the localized pressure variations in the medium. On these gradients the light is deflected, which leads to variations in photon density in the fluorophores (mechanism 6 in Fig. 2), therefore the intensity of the fluorescence is modulated. Also, the scattering coefficient is modulated (mechanism 5), which causes fluorescence intensity modulation by variation of the photon density distribution.24
Yuan et al.76,7879.80.–81 investigated the combination of fluorescence with AO and derived a mathematical model. They proposed two mechanisms that can explain the modulation of fluorescence. The first mechanism for low concentrations is similar to the explanation of Kobayashi et al.,24 the excitation light is modulated and thus also the fluorescent light. For high concentrations, they propose a modulation in quenching rate and thus in detected intensity. At high concentrations the variation in distance between fluorophores leads to variations in quenching, giving a modulation of fluorescence that is inverted compared to lower concentrations where quenching is absent. The signal strength found by Yuan80 shows a linear relation with the applied voltage on the ultrasound transducer, which is usually linear with the US pressure, and the signal was much weaker than found by Kobayashi et al.
Yuan et al. used a photo-multiplier tube to detect the modulation of the fluorescence and found that this signal depends quadratic on the acoustic pressure. Hall et al.77 showed a parallel detection scheme to detect the modulation in fluorescence by modulating the gain of the CCD with the US frequency. The phase difference between the US and CCD gain modulation was 0 and 180 deg.
Acousto-Optic Assisted Elastograpy
AO elastography is an imaging modality for quantifying mechanical properties. Kirkpatrick et al.82 stimulated tissue with a low frequency acoustic force (1 to 5 Hz), which induces a strain. The surface of the tissue, in their case the skin, will give a dynamic shift in the speckle pattern in the back-reflected laser light with this low frequency. The stiffer the material, the less shift in the speckle pattern is observed.
A high-intensity US burst is focused into the object on the region of interest. Inside the focus a shear wave is generated that radiates away from the focal point with a speed of typically 1 to , which depends on the viscoelastic properties of the medium. This low speed causes a speckle decorrelation on ms time scale in transverse direction. This technique has a lateral resolution in the order of mm and is not limited to use at the surface. Li et al.85 showed that by using the acoustic radiation force the resolution of a measurement can be increased by 40% to 110%. Sing et al.86 used two transducers, which caused a beat and shear wave frequency of 250 Hz. This low frequency shear wave was detectable with a speckle contrast measurement.
Acousto-Optic Assisted Light Focusing in Turbid Media
Xu et al.87 presented another detection technique called time-reversed ultrasonically encoded (TRUE) optical focusing based of a PRC used as a phase-conjugate mirror. The setup uses three arms: one signal arm where the object is located and two reference arms. The first reference arm is used for the creation of the hologram inside the PRC; the second illuminates the crystal from the opposite direction. By the hologram inside the crystal, the light is focused inside the object with the focal point inside the US focus region. This increases the local fluence at the US focus; in the future one can expect a significant improvement of the light level at the focus if the phase conjugate mirror exhibits a noticeable gain. This method makes AO the only technique that can create a guide star, which makes direct focusing of light possible in turbid media. Inside the US focus, this light is tagged and detected outside the medium. A second benefit from using this phase conjugate PRC setup is the improved spatial resolution with a factor of .
Ex Vivo and In Vivo Experiments
Various groups report measurements ex and in vivo with the use of AO.Kim et al.13 showed an image of chicken breast tissue with an embedded methylene-blue-dyed sentinel lymph node. Another ex vivo measurement on chicken breast is performed by Hisaka and Sasakura;88 they also used a reflection mode AO setup. Both experiments show the feasibility of a reflection mode AO setup on real tissue.
Kothapalli and Wang89 embedded mouse and rat blood vessels in tissue mimicking phantom material at a depth of 3 mm and imaged these vessels with an AO microscopy setup. They also tested the AO microscopy setup on a phantom.90 These experiments were performed ex vivo and show the possibility of measurements on tissue and have thus the correct scattering and absorption coefficients. Recently Lai et al.91 showed an ex vivo measurement on an HIFU induced lesion and compared the result with B-mode ultrasound. They concluded that AO sensing can follow the formation of the induced lesion in real time, and for noncavitating lesions AO gives a more robust signal compared with B-mode ultrasound. Murray et al.92 further investigated the changes in AO response of ex vivo chicken breast while elicited by a high-intensity ultrasound field. They concluded that with the use of AO it is possible to probe in real time the formation of a HIFU lesion.
In vivo experiments are performed by Lev et al.93,94 They demonstrate that ex vivo measurements suffer much less from speckle decorrelation than in vivo experiments. Further, they demonstrated the first AO tomography measurements on both mice and humans.
Discussion and Outlook
Tagged photons can be detected with several techniques. The photodetectors, which are fast enough to match the frequency of the US, have, in general, small detector surfaces. To avoid loss of SNR, averaging over multiple speckles must be avoided. These properties lead to a small etendue and thus a low count of tagged photons. CCD arrays can detect light from thousands of speckles in parallel and collect more light, resulting in an expected higher SNR. Unfortunately, these detectors are slow, therefore a lock-in detection is used, and only a few phases of the US are sampled.
Compared to untagged photons, tagged photons have a shifted optical frequency. Several spectrometer-based techniques are described in literature and have the advantage of greater stability. Small fluctuations in power, wavelength, and phase are canceled out assuming a long coherence length of the laser light. By using heterodyne detection it is possible to distinguish between speckle decorrelation by Brownian motion, shot-noise level, and the signal strength. Heterodyne detection also amplifies the signal optically to rise the signal above the level of the untagged photons.
PRC-based detection has a large etendue and the advantage of signal amplification. The disadvantage is the low response time of typical crystals, which makes it sensitive to speckle decorrelation. Spectral hole burning or confocal Fabry-Perot filters are insensitive to speckle decorrelation because they act as band pass filters and have the largest etendue.57 Crystal must be cryogenically cooled in order to achieve a narrow bandwidth in the MHz range, and Fabry-Perot need to be stabilized with a high precision feedback loop and isolated from mechanical vibrations.
To achieve high resolution in lateral direction, an ultrasound transducer with a small focal spot size is required. To achieve good axial resolution the frequency, phase and/or amplitude of the ultrasound can be made time dependent. This way only in a small part of the column, defined by the US beam, the light is tagged with these specific properties. The main advantage of giving a short US burst (amplitude) is the high acoustic peak pressure and therefore large modulation of light. The advantage of giving a chirped US pulse is the amount of information that is received in a single measurement; instead of measuring a single point, information can be extracted along the line in the US propagation direction.
The major challenges in in vivo AO applications are speckle decorrelation and low light levels. In general a detection method that detects as many tagged photons as possible is desired to obtain a good SNR. One strategy would be longer integration times; however, in vivo speckle decorrelation occurs on time scales of 1 ms. Most PRC have a response time orders of magnitude longer than 1 ms and are therefore not suitable for in vivo measurements. CCD-based techniques must acquire all light within the speckle decorrelation time, increasing shot noise and decreasing SNR. Spectral hole burning and Fabry-Perot interferometers are insensitive for speckle decorrelation and thus have none of these disadvantages; however, it is difficult to make sufficiently large etendue. Therefore these techniques are likely to be most suitable for in vivo imaging. Another more fundamental aspect in the introduction of acousto-Op. for in vivo application is the uncertainty regarding the exact information that is provided by AO. Although the spatial resolution of AO is potentially good, the property that it samples is related to the local fluence rate, which usually exhibits variations over a larger spatial scale. This aspect can be investigated using the new technologies that have been described in this review.
This work is performed under grant 09NIG01 of the Foundation for Fundamental Research on Matter (FOM), which is part of the Netherlands Organisation for Scientific Research (NWO).