The present invention relates to an optical device that requires optical resolution, and in particular, to a CARS microscope for focusing an optical beam and acquiring a response signal by changing the relative positions of the optical beam and an observation object irradiated with the optical beam.
Raman spectroscopic microscopes are quite effective to observe biologically-related samples. In a Raman spectroscopic microscope, an observation target is irradiated with a focused laser beam to detect Raman scattered light generated from the observation target. Raman scattered light has a shifted frequency from the frequency of the excitation light, and a Raman spectrum is measured with a spectrometer or the like. Scanning an observation target with an irradiation beam while changing their relative positions can obtain an optical spectrum at each position, and an image can be formed on the basis of such spectrum. A Raman spectrum at each observation position reflects the vibrational excited state of a molecule at the position, and thus is characteristic of the molecule. Using the characteristics of such a spectrum can know, if living cells are observed, a distribution of biomolecules in the cell tissue.
Measurement of the aforementioned Raman scattering takes a long time as the intensity of the obtained scattered light is weak. As a method that can obtain intense scattered light, there is known spectroscopy that uses nonlinear Raman scattering called CARS (Coherent Anti-Stokes Raman Scattering). Using such a method can also obtain a Raman spectrum and know the molecular vibrational state. To generate CARS, pulsed laser with high peak power is used. CARS is generated from such a pulsed laser beam due to the nonlinear effect, and the intensity of the CARS becomes orders of magnitude higher than that of Raman scattering as the peak power is higher. Accordingly, it is possible to obtain a signal with a high signal-noise ratio and significantly reduce the measurement time.
CARS is based on the third-order polarization. In order to generate CARS, a pump beam, a Stokes beam, and a probe beam are required. Typically, the pump beam is substituted for the probe beam in order to reduce the number of light sources. In that case, the induced third-order polarization is represented as follows.
P
AS
(3)(ωAS)=|χr(3)(ωAS)+χnr(3)|EP2(ωP)E*S(ωS) [Formula 1]
Herein, χr(3)(ωAS) is a resonant term of a vibration of a molecule with the third-order electric susceptibility, and χnr(3), which has no frequency dependence, is a nonresonant term. In addition, the electric fields of the pump beam and the probe beam are represented by EP, and the electric field of the Stokes beam is represented by ES. In Formula (1), the asterisk that appears in ES represents the complex conjugate. The intensity of a CARS beam is represented as follows.
I
CARS(ωAS)∝|PAS(3)(ωAS)|2 [Formula 2]
A mechanism by which a CARS beam is generated will be described using a molecular energy-level diagram (
Since Raman scattering was first discovered in 1928, a spectrum of a variety of molecules has been researched, and data thereon has been accumulated. Thus, it is desirable to identify molecules with reference to such spectral data. A CARS beam is represented by Formulae (1) and (2), and Im[χr(3)(ωAS)] is a portion corresponding to the Raman scattering spectrum. This is the complex portion of the resonant term, and interferes with the nonresonant term χnr(3) as described above. Thus, the shape of the spectrum obtained from CARS differs from that of the Raman scattering spectrum Im[χr(3)(ωAS)]. Therefore, it would be difficult to directly analyze a CARS spectrum with reference to the Raman scattering spectrum.
Development of a method for extracting a Raman scattering spectrum from a CARS spectrum is an important challenge to be addressed, and a variety of methods has been developed (Non Patent Literature 1). For example, the maximum entropy method, which is a method for restoring a phase spectrum from an intensity spectrum, includes determining a complex portion of a resonant term through mathematical computation. Alternatively, a method that uses interference is also known (Non Patent Literature 2).
As a spectral region that is sensitive to the molecular structure, there is a Raman scattering spectral region (of from 1800 to 800 cm−1) called a fingerprint region. For detection of a CARS beam, a spectrum in a similar region is desirably obtained. In the method introduced in Non Patent Literature 1, the spectral bandwidth of a Stokes beam for excitation is about 140 cm−1, which cannot cover such region. Non Patent Literature 3 introduces a method that uses a photonic fiber for a light source to cover such deficiency. Specifically, the method includes irradiating a photonic fiber with ultrashort pulsed laser to generate a broadband beam called a supercontinuum beam, and using it as a Stokes beam.
A method for observing living cells using a CARS beam is advantageous in that it is “noninvasive.” In order to generate a CARS beam, a measurement target is simultaneously irradiated with a pump beam and a Stokes beam that are ultrashort pulsed beams. Typically, as the wavelengths of the two excitation beams used for CARS, wavelengths that are not absorbed by living cells are used. Thus, the method can be said to be “noninvasive” under low peak power conditions and thus will hardly damage cells. However, if the peak power is increased too much, a multi-photon process may occur, which in turn may influence the cells. Thus, even though the method is “noninvasive,” the peak power of beams that irradiate living cells is desirably low. An object of the present invention is to improve the signal-to-noise ratio of a weak CARS signal that is generated under suppressed peak power conditions of excitation beams.
A CARS microscope for solving the above problem includes a first laser beam with a frequency ωP; a normal dispersion nonlinear optical fiber excited by the first laser beam; a second laser beam with a frequency ωST in a supercontinuum beam (hereinafter referred to as a SC beam) generated from the nonlinear optical fiber; a third laser beam as a reference beam with a frequency ωAS=2·ωP−ωST in the SC beam generated from the nonlinear optical fiber; an optical unit configured to align the first beam and the second beam on the same axis; a mechanism for adjusting the phases of the first beam and the second beam; an objective lens configured to focus the first and second laser beams; a scanning mechanism for scanning the observation sample; an objective lens configured to detect a CARS beam generated from the observation sample; an interference optical unit configured to cause the CARS beam and the third beam to interfere with each other; spectrometers each configured to disperse the interference beam; photodetectors each configured to detect the dispersed beam; a computing unit configured to process signals from the photodetectors; and a display device configured to display an image on the basis of information of the computing unit.
As a method that can distinguish between a change in the number of molecules in a measurement region and a change in a spectrum that occurs due to a change in the molecular structure, a method that uses a broadband beam as a Stokes beam is adopted. Although the present invention uses a SC beam as a broadband beam, such SC beam has not been used as a light source for measurement of homodyne interference, which may be able to improve the signal-to-noise ratio of a CARS signal, so far. This is because, in order to obtain a broadband, coherent SC beam, which is required to generate a CARS beam and is required for homodyne measurement, an excitation light source with a pulse width of several femtoseconds would be needed, which is not very realistic. In the present invention, a highly coherent SC beam is used as a light source to improve the signal-to-noise ratio of a CARS signal.
According to the present invention, it is possible to integrate interference beams of a CARS beam generated for each laser pulse and a reference beam on each detector, and thus improve the signal-to-noise ratio. In order to maintain the coherence for each pulse of a SC beam, which is generated without using a normal dispersion nonlinear optical fiber, it would be necessary to perform excitation using an ultrashort pulse with a pulse width on the order of 10 femtoseconds. In order to generate such an ultrashort pulsed laser beam, an expensive solid-state laser would be required at present. Therefore, CARS microscopes become expensive and the spread of CARS microscopes for purposes other than research becomes difficult. In contrast, if a normal dispersion nonlinear optical fiber is used, it becomes possible to maintain the coherence for each pulse of a laser beam with a pulse width of greater than or equal to 100 femtoseconds. In such a case, although the peak power of the laser beam needs to be higher than when optical fibers of other type dispersion are used, this can be addressed by using a fiber laser. Thus, a less expensive optical device configuration can be provided.
Other problems, configurations, and advantageous effects will become apparent from the following description of embodiments.
Hereinafter, embodiments for carrying out the present invention will be described with reference to the drawings.
The frequency range in which a SC beam is generated desirably reaches a maximum of ±3000 cm−1 of the frequency top of the first beam. This is because if an organic substance in a living cell is to be observed, a stretching vibration of 2930 cm−1 of CH3 or a stretching vibration of 2850 cm−1 of CH2 may be a target to be observed, and thus that such bands should be included in the frequency range. However, as a spectral range that is effective in identifying molecules is in the range of about 800 to 1800 cm−1 (i.e., fingerprint region), a SC beam in a narrow frequency range may be sufficient in some cases. The photonic crystal fiber may be provided with a polarization-maintaining optical fiber. In that case, the polarization direction is stabilized, and thus, the intensity and the spectral shape of a SC beam are also stabilized.
The SC beam is split in two by a dichroic mirror 257, with the frequency ωP as the boundary. That is, a beam with a low frequency shown by a region 163 in
A laser beam in the high-frequency region reflected by the dichroic mirror 257 is also a SC beam, and is used as a local beam including the frequency (2·ωP−ωST), that is, a reference beam (i.e., third beam). The reference beam is collimated by lenses 258 and 206, and passes through a polarization beam splitter 216 and a Fresnel rhomb waveplate 217 having the effect of a λ/4 plate, and is then returned to the Fresnel rhomb waveplate 217 by a mirror 218. The mirror 218 is used to adjust the optical path length. A laser beam that has passed through the Fresnel rhomb waveplate 217 is a p-polarized beam, which is then reflected by the polarization beam splitter 216, and travels toward the half beam splitter 213.
It follows that beams polarized in different directions enter the half beam splitter 213 from two directions, and the beams are split in two directions, so that interference beams are emitted in two directions. A method called phase-diversity detection is used to detect |EAS(ω)|. A Fresnel rhomb waveplate 221 having the effect of a λ/2 plate whose optical axis is tilted by 22.5 degrees is disposed for interference beams that are emitted to the right of the half beam splitter 213 on the paper surface. Then, the interference beams are focused onto spectrometers disposed at the focus positions by a condensing lens 215. A polarization beam splitter 223 is disposed in the optical path before the spectrometers, so that the interference beams are decomposed into components in the s direction and the p direction, which are then detected by spectrometers 106 and 108, respectively. Herein, it is assumed that the observation object is a point object in the optical axis on the focal plane, and the complex amplitude of a CARS beam with a frequency ω from the observation object and the complex amplitude of the reference beam are represented by EAS(ω) and ELO(ω), respectively. Provided that a differential signal of the spectrometers 106 and 108 at the respective wavelengths is IC(ω), the differential signal IC(ω) is represented as follows.
I
C(ω)=α|EAS(ω)|·|ELO(ω)|cos Φ(ω). [Formula 3]
Symbol a represents a coefficient including signal amplification, the efficiency of the spectrometers, and the like, and symbol Φ(ω) represents the phase difference between the CARS beam from the observation object and the reference beam. A Fresnel rhomb waveplate 222 having the effect of a λ/4 plate whose optical axis is tilted by 45 degrees is inserted for interference beams that are emitted in the upward direction of the half beam splitter 213 on the paper surface. The interference beams focused by a condensing lens 214 are detected by spectrometers 105 and 107. Specifically, the interference beams are separated into s-polarized beams and p-polarized beams by a polarization beam splitter 224 disposed in the optical path, which are then detected by the respective spectrometers. Herein, provided that a differential signal of the spectrometers 105 and 107 at the respective wavelengths is IS(ω), the differential signal IS(ω) is represented as follows.
I
S(ω)=α|EAS(ω)|·|ELO(ω)|sin Φ(ω). [Formula 4]
Only interference components are detected in IC(ω) and IS(ω).
A computing unit 109 performs computation represented as follows.
I(ω)=√(IC2(ω)+IS2(ω))=α|EAS(ω)|·|ELO(ω)| [Formula 5]
I(ω) is proportional to the amplitude of the CARS beam from the observation object and the amplitude of the reference beam. Thus, if the wavelength dependence of |ELO(ω)| is small, increasing |ELO(ω)| can obtain I(ω) with amplified |EAS(ω)|. Typically, the spectrum of a SC beam is not flat. Thus, in order to obtain a more accurate spectrum I(ω), it is necessary to perform correction using the amplitude spectrum of the SC beam.
Next, the complex components of the resonant term of the CARS beam are extracted to extract the Raman scattering spectrum. Φ(ω) that is the phase difference between the local beam and the CARS signal beam is represented by Φ(ω)=ωτ+θS(ω)+θinst(ω). Symbol ωτ represents the optical path difference between the two beams, θS(w) represents the phase difference due to a resonant beam, and θinst(ω) represents the phase difference derived from the device. Herein, it is assumed that the local beam has no frequency dependence. tan Φ(ω) is determined from Formulae (3) and (4), and Φ(ω) can also be determined. First, as an observation sample, a sample that generates only a nonresonant CARS beam is measured to determine ωτ+θinst(ω) Next, an observation sample that generates resonant CARS is measured. Accordingly, θS(ω) can be determined. Thus, the complex number portion of the resonant components can be determined as I(ω)sin θS(ω). Accordingly, a portion corresponding to the Raman scattering spectrum can be obtained. Detectors such as CCDs may be used for detection of beams with the spectrometers. The display device 110 displays the scanned position and the display position of the observation object 202 in association with each other. Displaying the complex components of a resonant beam at a frequency position that is characteristic of a molecular vibration can know a distribution of molecules.
In this embodiment, a number of pulses are integrated. If the coherence between pulses is lost, the phase of Φ(ω) in IC(ω)=α|EAS(ω)|·|ELO(ω)|cos Φ(ω) or IS(ω)=α|EAS(ω)∥ELO(ω)|sin Φ(ω) becomes random, so that the value of IC(ω) or IS(ω) obtained by integrating a number of pulses becomes zero. However, as a SC beam generated from a normal dispersion nonlinear optical fiber is used in this embodiment, the coherence between pulses is maintained. Thus, there is no possibility that IC(ω) or IS(ω) may become zero, which would otherwise occur if there is no coherence between pulses.
Although the embodiment shown in
In the aforementioned embodiment, a method of causing linearly polarized beams, which are orthogonal to each other, to interfere with each other is adopted to perform phase diversity detection. As an alternative method, it is also possible to use a method of converting a CARS beam from an observation object and a reference beam into a right-handed circularly polarized beam and a left-handed circularly polarized beam that are orthogonal to each other, and causing the beams to interfere with each other. For example, a CARS beam is converted into a right-handed circularly polarized beam, and a reference beam is converted into a left-handed circularly polarized beam. When the optical axes of analyzers for detecting the interference beams are set to 0, 45, 90, and 135 degrees, respectively, it is possible to obtain beams with phase differences of 0, 90, 180, and 270 degrees. Combining such signals can obtain a signal represented by Formula (5) and obtain effects that are similar to those in the aforementioned embodiment.
When the optical shutters 219 and 209 are open, a signal represented by the following formula is output from the spectrometer 105,
S
C(ω)=|ELO|2+|EAS(ω)|2+2|ELOEAS(ω)|cos Φ(ω), [Formula 6]
and
a signal represented by the following formula is output from the spectrometer 107.
S
S(ω)=|ELO|2+|EAS(ω)|2+2|ELOEAS(ω)|sin Φ(ω). [Formula 7]
|ELO|2 and |EAS(ω)|2 in Formula (6) and Formula (7) can be determined by shutting off one of them. When the optical shutter 219 is closed and the optical shutter 209 is opened, |ELO(ω)|2 is output to the spectrometers 105 and 107. To the contrary, when the optical shutter 219 is opened and the optical shutter 209 is closed, |EAS(ω)|2 is output to the spectrometers. The computing unit 109 computes |ELO(ω)EAS(ω)| from the outputs.
Further spectral correction can be performed by measuring the Stokes beam (i.e., second beam). Though not shown, a plane mirror is inserted immediately before the objective lens 201, and an optical spectrum is measured by the spectrometers 105 and 107 with the optical shutter 219 open and the optical shutter 209 closed. The optical spectrum includes the optical spectrum of the Stokes beam ES(ω). Thus, the influence of the spectral distribution of the Stokes beam can be corrected by taking into Formula (1) into consideration.
Using the results, a phase difference generated by the resonant beam is computed with the aforementioned method that uses interference, so that [|ELO(ω)EAS(ω)|sin θS(ω)] that is the complex component of the resonant term is extracted to obtain a result equivalent to that of Raman spectroscopy.
In this embodiment, a bandpass filter 220 is inserted in the optical path of a pump beam that is the first beam reflected by the beam splitter 251. As a laser beam emitted from the pulsed light source 141, a laser beam with a narrow pulse width is used to maintain the coherence between pulses. Therefore, the spectral bandwidth of the first beam is wide. However, if the first beam is used as it is to generate a CARS beam, a desired spectral resolution may not be obtained in some cases. To address this, a bandpass filter for narrowing the spectral bandwidth was inserted.
The present invention is not limited to the aforementioned embodiments, and includes a variety of variations. For example, although the aforementioned embodiments have been described in detail to clearly illustrate the present invention, the present invention need not include all of the structures described in the embodiments. It is possible to replace a part of a structure of an embodiment with a structure of another embodiment. In addition, it is also possible to add, to a structure of an embodiment, a structure of another embodiment. Further, it is also possible to, for a part of a structure of each embodiment, add, remove, or substitute a structure of another embodiment.
According to the present invention, it is possible to acquire a high-resolution image using a CARS beam, and provide a noninvasive optical device for measuring a distribution of biomolecules or a change in the distribution.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2013/072349 | 8/22/2013 | WO | 00 |