SYSTEMS AND METHODS FOR STIMULATED BRILLOUIN MICROSCOPY
20220326080 · 2022-10-13
Inventors
Cpc classification
G01N21/636
PHYSICS
H01S3/108
ELECTRICITY
International classification
Abstract
Systems and methods are provided for analyzing a biomechanical property of a medium using stimulated Brillouin scattering microscopy. The method can include a first step of applying a probe beam and pulsed pump beam to a target section of the medium, wherein the pump beam interacts with the probe beam to generate at least one acoustic wave in the medium and at least one Brillouin signal is produced as a result of the generated acoustic wave. The method can also include a second step of receiving the produced Brillouin signal and a third step of determining, using a processor and the
Claims
1. A method of analyzing at least one mechanical property of a medium using Brillouin microscopy, the method comprising: applying a probe beam and pulsed pump beam to a target section of the medium, wherein the pump beam interacts with the probe beam to generate at least one acoustic wave in the medium and at least one Brillouin signal is produced as a result of the generated acoustic wave; receiving the produced Brillouin signal; and determining, using a processor and the Brillouin signal, information associated with at least one mechanical property of the target section.
2. The method of claim 1, wherein the acoustic wave is amplified based on the interaction between the pump beam and the probe beam.
3. The method of claim 1, wherein the probe beam is pulsed.
4. The method of claim 3, wherein the probe beam pulses are applied at twice the rate that the pump beam pulses are applied to the target section.
5. The method of claim 1, wherein the total time averaged power of both the pulsed probe beam and the pump beam is less than 10 mW.
6. The method of claim 1, wherein the ratio of the time averaged power of the pulsed probe beam to the time averaged power of the pump beam is between 1:3 and 1:5.
7. The method of claim 1, wherein the pulse width of the pump beam is between 2 and 5 nanoseconds.
8. (canceled)
9. (canceled)
10. The method of claim 1, wherein determining the at least one biomechanical property of the medium comprises compensating for the relative intensity noise of the probe beam.
11. The method of claim 1, wherein the at least one biomechanical property includes the longitudinal modulus the target section.
12. The method of claim 1, wherein the integration time of the Brillouin signal is less than 100 μs.
13. (canceled)
14. A Brillouin microscopy system for analyzing at least one mechanical property of a medium, the system comprising: a probe optical source configured to apply a probe beam to a target section of the medium; a pump optical source configured to apply a pulsed pump beam to the target section of the medium, wherein the pulsed pump beam interacts with the probe beam to generate at least one acoustic wave in the medium and at least one Brillouin signal is produced as a result of the generated acoustic wave; a sensor configured to receive the produced Brillouin signal; and a processor configured to determine information associated with at least one mechanical property of the target section using the Brillouin signal.
15. The system of claim 14, wherein the at least one mechanical property includes the longitudinal modulus of the target section.
16. The system of claim 1, wherein sensor is configured to measure a transmitted form of the probe beam after it has been transmitted through the target section, and wherein the Brillouin signal is contained within the transmitted probe beam.
17. The system of claim 1, further comprising: at least one quarter wave plate configured to alter the polarization of the probe beam; and a polarizing beam splitter configured to redirect the transmitted probe beam to the sensor.
18. The system of claim 14 further comprising: a filter configured to prevent back-reflected pump beam light from contacting the sensor.
19. (canceled)
20. The system of claim 14, wherein the probe optical source comprises a continuous wave laser and an electro-optic modulator.
21. (canceled)
22. (canceled)
23. The system of claim 14, wherein the probe optical source is configured to tune the frequency of the probe beam using a feedback signal provided from a second photo-diode sensor.
24. (canceled)
25. The method of claim 14, wherein the sensor integrates the Brillouin signal and wherein the integration of the electrical signal from photo-diode sensor is performed using a boxcar integrator.
26. The system of claim 14, wherein the pump optical source and the probe optical source are configured to apply a pump beam and a probe beam that are approximately counter-propagating with one another.
27. (canceled)
28. (canceled)
29. (canceled)
30. The system of claim 1, further comprising: a spatial pinhole to select Brillouin signal from a confocal point in the target section.
31. A method of measuring shear modulus information of a medium using Brillouin microscopy, the method comprising: applying a probe beam and pump beam to a target section of the medium, wherein the pump beam interacts with the probe beam to generate at least one acoustic wave in the medium and at least one Brillouin signal is produced as a result of the generated acoustic wave; receiving the produced Brillouin signal; and determining, using a processor and the Brillouin signal, information associated with the shear modulus of the target section.
32. The method of claim 31, wherein a single laser is used as an initial light source for both the probe beam and pump beam.
33. (canceled)
34. The method of claim 31, wherein the pump beam is formed of light modulated with a square waveform.
35. The method of claim 31, wherein the polarization state of the probe beam is orthogonal to the polarization state of the pump beam when applied to the target section.
36-42. (canceled)
43. A Brillouin microscopy system for measuring shear modulus information of a medium, the system comprising: a probe optical source configured to apply a probe beam to a target section of the medium; a pump optical source configured to apply a pump beam to a target section of the medium, wherein the pump beam interacts with the probe beam to generate at least one acoustic wave in the medium and at least one Brillouin signal is produced as a result of the generated acoustic wave; a sensor configured to receive the produced Brillouin signal; and a processor configured to determine information associated with the shear modulus of the target section using the Brillouin signal.
44-46. (canceled)
47. The system of claim 43, wherein the probe optical source comprises at least one acousto-optic frequency shifter configured to modify the frequency of the probe beam light.
48. (canceled)
49. The system of claim 43, wherein sensor is configured to measure a transmitted form of the probe beam after it has been transmitted through the target section, wherein the Brillouin signal is contained within the transmitted probe beam.
50-53. (canceled)
54. A method of measuring longitudinal and shear moduli information of a medium using Brillouin microscopy, the method comprising: applying a first probe beam and a pulsed first pump beam to a target section of the medium, wherein the first pump beam interacts with the first probe beam to generate at least one first acoustic wave in the medium and at least one first Brillouin signal is produced as a result of the generated acoustic wave; applying a second probe beam and a second pump beam to the target section of the medium, wherein the second pump beam interacts with the second probe beam to generate at least one second acoustic wave in the medium and at least one second Brillouin signal is produced as a result of the generated acoustic wave; receiving the produced first and second Brillouin signals; and determining, using a processor and the first and second Brillouin signals, information associated with both the longitudinal modulus and shear modulus of the target section.
55. The method of claim 54, wherein the first probe beam and first pump beam are not applied to a target section of the medium at the same time that the second probe beam and second pump are applied to the target section of the medium.
56. The method of claim 54 further comprising: switching between applying the first probe and pump beams and applying the second probe and pump beams to the target section.
57. The method of claim 54, wherein the first Brillouin signal is used when determining the information associated with the longitudinal modulus of the target section.
58. The method of claim 54, wherein the second Brillouin signal is used when determining the information associated with the shear modulus of the target section.
59. A Brillouin microscopy system for measuring the longitudinal and shear moduli information of a medium, the system comprising: a longitudinal configuration comprising: a first probe optical source configured to apply a first probe beam to a target section of the medium; a first pump optical source configured to apply a pulsed first pump beam to the target section of the medium, a shear-wave configuration comprising: a second probe optical source configured to apply a second probe beam to the target section of the medium; a second pump optical source configured to apply a second pump beam to the target section of the medium, wherein the first pump beam interacts with the first probe beam or the second pump beam interacts with the second probe beam to generate at least one acoustic wave in the medium and at least one Brillouin signal is produced as a result of the generated acoustic wave; a sensor configured to receive the produced Brillouin signal; and a processor configured to determine information associated with at least one of a longitudinal modulus or a shear modulus of the target section.
60. The system of claim 59, wherein the first probe beam and first pump beam are not applied to a target section of the medium at the same time that the second probe beam and second pump beam are applied to the target section of the medium.
61. The system of any one of the preceding claims, further comprising: an adjustment system configured to switch between applying the first probe and pump beams and applying the second probe and pump beams to the target section.
62-63. (canceled)
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0023] The present disclosure will hereafter be described with reference to the accompanying drawings, wherein like reference numerals denote like elements.
[0024]
[0025]
[0026]
[0027]
[0028]
[0029]
[0030]
[0031]
[0032]
[0033]
[0034]
[0035]
DETAILED DESCRIPTION
[0036] The systems and methods provided herein are not limited to the particular embodiments described. The terminology used herein is selected for the purpose of describing particular embodiments only, and is not intended to be limiting. The scope of the present invention will be limited only by the claims. As used herein, the singular forms “a”, “an”, and “the” include plural embodiments unless the context clearly dictates otherwise.
[0037] It should be apparent to those skilled in the art that many additional modifications beside those already described are possible without departing from the inventive concepts. In interpreting this disclosure, all terms should be interpreted in the broadest possible manner consistent with the context. Variations of the term “comprising” should be interpreted as referring to elements, components, or steps in a non-exclusive manner, so the referenced elements, components, or steps may be combined with other elements, components, or steps that are not expressly referenced. Embodiments referenced as “comprising” certain elements are also contemplated as “consisting essentially of” and “consisting of” those elements.
[0038] The present disclosure overcomes the aformentioned drawbacks of prior systems by using stimulated Brillouin scattering (SBS). SBS is a nonlinear optical phenomenon used in the physical sciences and waveguide devices. Thus far, SBS has been considered unsuitable for biomedical applications because it requires too high optical power (exceeding 270 mW) for biological samples.
[0039] Spontaneous Brillouin scattering is generated by acoustic waves that arise from thermal motion of molecules. In continuous medium, these vibrations can be viewed as collective vibrational modes at all frequencies propagating in all directions as acoustic waves, equivalent to acoustic phonons in solid state matter. The speed of acoustic propagation in soft tissues is about 2,000 to 2,500 m/s (V.sub.L) for longitudinal waves and 3 to 30 m/s for shear waves (v.sub.S). The acoustic waves form refractive index gratings via the photo-elastic effect. Among these thermodynamically generated acoustic modes, those with acoustic wavelengths approximately equal to half of the optical wavelength (λ.sub.o/n) in the medium (with a refractive index, n) give rise to phase-matched Bragg reflection of light. This phenomenon is spontaneous Brillouin scattering. The frequency of the scattered light is shifted by the Doppler effect. The Brillouin frequency shift is given by: f.sub.B=(2n/λ.sub.o)v.sub.L,S cos(θ/2) where θ is the scattering angle: 180° for backward scattering and 90° for the cross-axis geometry. The linewidth of the Brillouin scattering spectrum is the inverse of the lifetime of acoustic waves and therefore is related to the viscosity of the medium at the acoustic frequency. The scattering mechanism of SBS is identical except that the acoustic waves are generated by interference of a pump and a probe beam via electrostriction (i.e. electric-field dependent refractive index change).
[0040] The systems and methods described herein can use pump and probe light in order to produce stimulated Brillouin scattering. The techniques described herein can be integrated with standard confocal microscope platforms. The resulting systems can readily incorporate fluorescence contrast for labelling specific cells, subcellular compartments, and specific tissues in living organisms. Thus, the SBS microscopy techniques described herein offer a leap forward from traditional spontaneous Brillouin microscopy methods, and can be used as high-resolution, comprehensive tools for mechanical investigations of biomaterials, transparent ocular tissues, fresh tissue slices, organoids, cells, small animal models, and the like.
[0041] As further described, for longitudinal SBS, an intensity modulation scheme can be used to generate time-synchronized, nanosecond pulses for pump and probe light. Calculations show that this can generate >1,000 times higher signals than spontaneous Brillouin scattering at equivalent optical powers. The pulsed operation can represent ˜100 times higher efficiency than is possible using continuous-wave (CW) schemes.
[0042] Shear-wave SBS has not yet been developed. To date, only longitudinal SBS has been successfully demonstrated. In some aspects, the systems and methods of the present disclosure allow for both longitudinal and shear moduli to be quasi-simultaneously obtained from samples at microscopic resolution. Thus, the SBS microscopy techniques described herein can in principle determine the complete stiffness matrix: that is, both the Lame parameters, λ and μ, of isotropic materials; and all 4 coefficients for fibrillary structures with hexagonal symmetry.
[0043] The primary output of Brillouin spectroscopy is the Brillouin frequency (the other output is acoustic loss via the linewidth). SBS microscopy can map the Brillouin frequency by scanning the focal volume at which the pump and probe beams meet. As this technique essentially measures the acoustic speed, it is analogous to ultrasound shear wave imaging. One difference, however, is that SBS microscopy uses focused light to generate acoustic waves whereas the ultrasound technique excites ultrasound using a transducer that requires physical contact. Another major difference is that the phase-matching acoustic frequency involved in SBS is several GHz for longitudinal waves and MHz for shear waves. In both cases, the phase-matching acoustic wavelength can be equal to a half of the optical wavelength, which is, for example, 280 nm for λ.sub.o=780 nm and n=1.4. This submicron resolution contrasts with longer ultrasound wavelengths in the range of ˜200 μm to 1 mm. The acoustic frequency range used in Lamb-wave OCT elastography is 0.3-1 kHz, and the corresponding acoustic wavelength is rather long, 2-5 mm. The shear modulus of isolated single cells has been measured recently using vibrational pipettes, but this technique is applicable only to isolated cells or small transparent objects with low modulus (<<100 Pa). By contrast, the SBS microscopy techniques described herein may be uniquely suited for non-contact, microscopic mapping in hydrogels, 2D and 3D cell cultures, tissue slices, and organoids.
[0044]
[0045] In the method 100, the at least one biomechanical property may include the longitudinal modulus of the target section. The probe beam may also be pulsed. If the probe beam is pulsed, the probe beam pulses may be applied at twice the rate that the pump beam pulses are applied to the target section. The total time averaged power of the pulsed probe beam and the pulsed pump beam may be less than 10 mW. The pulsed probe beam and the pump beam may both have a total time averaged power or energy density below the maximum permissible exposure limits for a human tissue. The pulsed probe beam and the pump beam may both have a total time averaged power or energy density below the maximum permissible exposure limits for human ocular tissue. The maximum permissible exposure limit for the pulsed probe beam and the pump beam may be calculated using the American National Standards Institute (ANSI) Z136 guidelines. The ratio of the time averaged power of the pulsed probe beam to the time averaged power of the pulsed pump beam may be between 1:2 and 1:6, between 1:3 and 1:5, or about 1:4. The pulse width of the pump beam may be between 2 and 5 nanoseconds.
[0046] Although not depicted, the method 100 may comprise a fourth step of displaying the at least one biomechanical property of the target section of the medium. The medium may be tissue, or more specifically, mammalian tissue or mammalian cells. The method step 106 may comprise compensating for the relative intensity noise of the probe beam. The acoustic wave may be amplified based on the interaction between the pump beam and the probe beam. The pump beam may be in counter-propagation with respect to the probe beam. The Brillouin signal may have a signal-to-noise ratio of at least 50 dB, at least 55 dB, at least 60 dB, or at least 65 dB. The Brillouin signal may be contained within a transmitted form of the probe beam.
[0047]
[0048] The at least one biomechanical property analyzed by the system 200 may include the longitudinal modulus of the target section 209. As shown, the sensor 206 may be configured to measure a transmitted form of the probe beam 222 after it has been transmitted through the target section 209, wherein the Brillouin signal is contained within the transmitted probe beam 222. In one form, both the pump beam 214 and the probe beam 212 may be vertically polarized. In order to achieve this vertical polarization, system 200 further include at least one quarter wave plate configured to alter the polarization of the probe beam 212. For instance, a pair of quarter wave plates may be positioned in the optical paths of both the pump beam 214 and probe beam 212. Both the pump beam 214 and probe beam 212 may be circularly polarized at the medium 208. After the medium, the probe beam 212 may continue toward the shared path of pump beam 214 and convert into horizontally polarized light after passing through a quarter wave plate. The system may include at least one polarizing beam splitter configured to redirect the transmitted probe beam 222 to the sensor. For instance, the probe beam 212 may be redirected to the sensor 206 through a polarizing beam splitter between the path of pump beam 214 and probe beam 222. Although not depicted, the system 200 may further comprise a filter positioned between the medium 208 and the sensor 206 and configured to prevent back-reflected pump beam light from contacting the sensor 206.
[0049] The probe beam 212 may be pulsed. In one form, the probe optical source 202 may comprise a continuous wave laser and an electro-optic modulator. Similarly, the pump optical source 204 may comprise a continuous wave laser and an electro-optic modulator. The sensor 206 may be a photo-diode sensor. In one form, the electrical signal provided by the photo-diode sensor may be a pulse train, where the amplitude of the pulse represents the addition of the probe power and the Brillouin signal. The integration of the electrical signal from photo-diode sensor may be performed using a boxcar integrator. The boxcar integrator may enhance the signal to noise ratio of the transmitted prove beam 222. The boxcar integrator may be specifically triggered at the pulse repetition rate and its gating window may be configured to coincide with the pulse duration of the transmitted probe beam 222. The integrated signal may be subsequently provided to an analog to digital converter before reaching the processor 210. The probe optical source 202 may be configured to tune the frequency of the probe beam 212 using feedback signal provided from a second photo-diode sensor. The probe optical source 202 may be configured to provide the pulsed probe beam with a total time averaged power or energy density below the maximum permissible exposure limits for human tissue. The probe optical source 202 may be configured to provide the pulsed probe beam with a total time averaged power or energy density below the maximum permissible exposure limits for human ocular tissue. The maximum permissible exposure limit for the pulsed probe beam may be calculated using the American National Standards Institute (ANSI) Z136 guidelines.
[0050] In the depicted system 200, the pump optical source 204 and the probe optical source 202 are configured to apply a pump beam 214 and a probe beam 212 that are counter-propagating with one another. The pump optical source 204 may have a peak power of less than 1 Watt. The pump optical source 204 may be configured to provide the pulsed probe beam with a total time averaged power or energy density below the maximum permissible exposure limits for human tissue. The pump optical source 204 may be configured to provide the pulsed probe beam with a total time averaged power or energy density below the maximum permissible exposure limits for human ocular tissue. The maximum permissible exposure limit for the pump beam may be calculated using the American National Standards Institute (ANSI) Z136 guidelines.
[0051] Although not depicted, the system 200 may further comprise a display in electrical communication with the processor 210 and configured to display the information associated with the at least one biomechanical property of the target section 209. As shown, a positioning system 215 may be configured to move the target section 209 of the medium with respect to the probe optical source 202 and pump optical source 204.
[0052]
[0053] In the method 300, a single laser may be used as an initial light source for both the probe beam and pump beam. The pump beam may be formed of light modulated with a square waveform. The polarization state of the probe beam may be orthogonal to the polarization state of the pump beam when applied to the target section. The total power of the probe beam and the pump beam may be less than 10 mW. The probe beam and the pump beam may both have a total time averaged power or energy density below the maximum permissible exposure limits for a human tissue. The probe beam and the pump beam may both have a total time averaged power or energy density below the maximum permissible exposure limits for human ocular tissue. The maximum permissible exposure limit for the probe beam and the pump beam may be calculated using the American National Standards Institute (ANSI) Z136 guidelines. The ratio of the power of the probe beam to the power of the pump beam is between 1:2 and 1:6, between 1:3 and 1:5, or about 1:4.
[0054] Although not depicted, the method 300 may comprise a fourth step of displaying the information associated with the shear modulus of the target section. The medium may be tissue, or more specifically, mammalian tissue or mammalian cells. The method step 306 may comprise compensating for the relative intensity noise of the probe beam. The acoustic wave may be amplified based on the interaction between the pump beam and the probe beam. The Brillouin signal may be contained within a transmitted form of the probe beam.
[0055]
[0056] The probe optical source 402 and the pump optical source 404 may share a single laser as an initial light source for both the probe beam 412 and pump beam 414. If this is the case, the single laser may be a continuous-wave laser. The pump optical source 404 may comprise an electro-optic modulator configured to modulate the pump beam light to have a square waveform. The probe optical source 402 may comprise at least one acousto-optic frequency shifter configured to modify the frequency of the probe beam light. The sensor 406 may be a photo-diode sensor. The sensor may be configured to measure a transmitted form of the probe beam 422 after it has been transmitted through the target section 408, wherein the Brillouin signal is contained within the transmitted probe beam 422. Although not depicted, the system 400 may further comprise a display in electrical communication with the processor 410 and configured to display the information associated with the shear modulus of the target section. As shown, the system 400 may comprise a positioning system 415 configured to move the target section 409 of the medium with respect to the probe optical source 402 and pump optical source 404.
[0057]
[0058] In the method 500, the first probe beam and first pump beam may not be applied to a target section of the medium at the same time that the second probe beam and second pump are applied to the target section of the medium. Although not depicted, the method may further comprise an additional method step of switching between applying the first probe and pump beams and applying the second probe and pump beams to the target section. The first Brillouin signal may be used when determining the information associated with the longitudinal modulus of the target section. Likewise, the second Brillouin signal may be used when determining the information associated with the shear modulus of the target section. The multiple probe beams and pump beams of the method 500 may all have a total time averaged power or energy density below the maximum permissible exposure limits for a human tissue, or more specifically human ocular tissue. The maximum permissible exposure limit for the multiple probe beams and pump beams may be calculated using the American National Standards Institute (ANSI) Z136 guidelines.
[0059]
[0060] The first probe beam 611 and first pump beam 613 may not be applied to a target section 609 of the medium at the same time that the second probe beam 612 and second pump beam 614 are applied to the target section 609 of the medium. As shown, the system 600 may comprise an adjustment system 618 configured to switch between applying the first probe and pump beams (611, 613) and applying the second probe and pump beams (612, 614) to the target section 609. As shown, the adjustment system may include flip mirrors. Also as shown, the system 600 may further comprise a positioning system 615 configured to move the target section 609 of the medium with respect to the longitudinal and shear-wave configurations (620, 630).
EXAMPLES
[0061] The following examples are provided in order to demonstrate and further illustrate certain embodiments and aspects of the present disclosure and are not to be construed as limiting the scope of the disclosure.
Example 1
Longitudinal Stimulated Brillouin Scattering
[0062] The following experimental and prophetic remarks relate to the development and characterization of a system and method for longitudinal stimulated Brillouin scattering.
[0063] Without being bound by theory, longitudinal-wave Brillouin spectroscopy determines f.sub.L which is proportional to v.sub.L. The acoustic speed is related to the compliance coefficient along the acoustic propagation direction: if the direction is in the x-axis, c.sub.11=ρv.sub.L.sup.2 (if z-axis, c.sub.33=ρv.sub.L.sup.2). For isotropic materials, c.sub.11=λ+2μ (=c.sub.33), and c.sub.11 is equal to longitudinal modulus M=ρ.sub.0v.sub.L.sup.2. M is related to bulk modulus K and shear moduli G: M=K+4/3G. ρ.sub.0v.sub.S.sup.2=G for shear waves. In soft tissues, G K, so M≈K. K=1/β where β=−(1/V)(∂V/∂.sub.p) is volume compressibility, where V is volume and p hydrostatic pressure; therefore, f.sub.L≈2 (2n/λ.sub.o)√{square root over (1/(ρβ))}(Eq. 1). Water-containing soft materials, such as hydrogels and tissues, can be modeled as a two-part mixture of liquid water plus a “solid” part. The latter is mostly hydrated proteins including collagen fibers. If the contribution of solid-liquid interactions is neglected, the liquid and solid parts additively contribute to the physical properties of the tissue. Let V.sub.w denote the volume of water in tissue, and V.sub.s the volume of the other constituents. Since the total volume is the sum of the two parts, V=V.sub.w+V.sub.s, the volume fraction of water is expressed as f.sub.w=V.sub.w/(V.sub.w+V.sub.s), and the volume fraction of the solid part is 1−f.sub.w. In terms of the volume fraction of water, the refractive index, density, and compressibility can be described as: n=n.sub.wf.sub.w+n.sub.s(1−f.sub.w), ρ=ρ.sub.wf.sub.w+ρ.sub.s(1−f.sub.w); and β=β.sub.wf.sub.w+β.sub.s(1−f.sub.w), where n.sub.w, ρ.sub.w and β.sub.w are the refractive index, density, and compressibility of water, and n.sub.s,ρ.sub.s and β.sub.s are these properties of the solid components.
[0064] From Eq. (1), =Δf.sub.L/f.sub.L=Δn/n−0.50ρ/ρ−0.5Δβ/β. This can be expressed in terms of the change of hydration H, which is defined as H=(total tissue weight−dry weight)/dry weight, and the compressibility of the solid part, β.sub.s: Δf.sub.L/f.sub.L≈−0.08ΔH/H−0.1Δβ.sub.s/β, when f.sub.w≈0.75−0.8 and H≈3.2. The equation indicates that the longitudinal Brillouin frequency provides information about hydration and the compressibility β.sub.s of the protein constituents. For corneal stroma, collagen fibrillary structure is the dominant contributor to β.sub.s.
[0065] For cells, longitudinal Brillouin frequency is sensitive to changes in solid-liquid volume ratio and the compressibility of the cytoplasmic macromolecules and cytoskeleton. While cell volume can vary as much as 100%, the density variability of cells is small, <1%. A small difference in cell density can affect cell biology and is universally related to shear modulus. The interaction between the liquid and solid components within the cytoplasm regulates the rheological behavior of cells and can play a prominent role in how cells deform and move. In turn, the hydro-mechanical properties of cells can influence intracellular water transport, cytoskeletal network modulation and cell volume regulation.
[0066] The Brillouin scattering cross section in isotropic liquids is given by
where K is bulk modulus, k.sub.B the Boltzmann constant, T temperature, V volume of the interaction region, θ the scattering angle, and γ=ρ∂ϵ/∂ρ the electrostriction constant where ε is dielectric constant and ρ is density. For isotropic materials, γ is related to the photo-elastic coefficient p.sub.12 for longitudinal acoustic waves and p.sub.44 for shear acoustic waves. For longitudinal waves, γ.sup.2≈1. The theoretical maximum collection efficiency of Brillouin scattering in the on-axis confocal configuration is given by
In a lossless system, S.sub.c≈5.4×10.sup.−11. The theoretical shot-noise-limited SNR is: SNR (sponatneous, longitudinal)=0.5 ρα.sub.sS.sub.cP.sub.0≈2.7×10.sup.−11 ηα.sub.sP.sub.0 (Eq. 2) where η is the photon-to-electron conversion efficiency of a photodetector or CCD camera, α.sub.s is optical loss in the system, and P.sub.0 is the input optical power. Experimentally, we measured α.sub.s≈0.006, which includes a ˜10-dB optical loss in a spectrometer and 4 dB loss in fiber coupling. We found that current spontaneous Brillouin microscopy has reached the theoretical SNR limit within a factor of <10.
[0067] We have considered the SBS process between a pump and a probe beam. The SBS gain coefficient is g.sub.0=γ.sup.2ω.sup.2/(nv.sub.L,Sc.sup.3ρ.sub.0Γ.sub.B) (Eq. 3), where γ.sub.B is the Brillouin linewidth, typically Γ.sub.B/2π=˜400 MHz. We find g.sub.0=3.5×10.sup.−5 (0.78/λ.sub.0).sup.3 μm/W. This theoretical value agrees well with the experimentally measured gain value of water, g.sub.0=3.4×10.sup.−5 μm/W at λ.sub.0=780 nm. The SBS gain is given by g=g.sub.0 I.sub.1Δz, where I.sub.1 is the pump beam intensity and Δz is the acousto-optic interaction length. The probe beam intensity is described by: I.sub.2(output)=I.sub.2(input)*(1+g). For focused beams with a Gaussian profile, we can define an effective beam waist area, A=Δx.sup.2, where Δx=0.4λ.sub.oNA and NA is numerical aperture in a medium with refractive index n. The effective confocal interaction length is Δz≈0.8 z.sub.R, where z.sub.R=nλ.sub.o/(NA).sup.2. Using I.sub.1=P.sub.1/A we obtain g=3.1×10.sup.−4 P.sub.1 [W.sup.−1] at λ.sub.o=0.78 μm (Eq. 4). In a lossless system, the optical signal at a photodetector is then P.sub.2(out)=P.sub.2(in) {1+g cos(2πf.sub.1t)}, where f.sub.1 is the modulation frequency of the probe (˜1 MHz, see below). The shot-noise-limited SNR is SNR (longitudinal SBS)=0.5 g.sup.2ηP.sub.2 (Eq. 5).
[0068] With Eq. (2) and Eq. (5), we can compare the SNR of the spontaneous and stimulated processes, using P.sub.0=P.sub.1
+
P.sub.2
, where the bracket indicates time average. For a given total power, it can be shown that shot-noise-limited SNR in SBS is maximized when I.sub.1=2I.sub.2 and, thus, P.sub.0=3
P.sub.2
. The condition for SBS to have higher gain than spontaneous scattering is: g.sup.2P.sub.2>6×10.sup.−11 P.sub.0 and thusg>1.3×10.sup.−5, from Eq. (4), the pump power required to produce this much gain is P.sub.1>42 mW, and the total power,
P.sub.1
+
P.sub.2
, should be greater than 64 mW. P.sub.0 is usually limited to maximum 10 mW (at 780 nm) to avoid thermal damage to samples. Then, max
P.sub.1
=6.7 mW and max
P.sub.2
=3.3 mW. Therefore, when CW light is used for SBS, spontaneous Brillouin microscopy has higher SNR unless the total pump power exceeds the maximum permissible exposure (MPE) level of ˜32 mW at λ.sub.o=0.78 μm for skin and corneal tissues (and cells). Consequently, stimulated Brillouin scattering microscopy using continuous wave light is of limited value.
[0069] In our attempt to overcome the limitations of spontaneous Brillouin microscopy, we discovered that stimulated Brillouin scattering microscopy with a pulsed pump and probe light could yield improved results while still operating within the maximum permissible exposure levels of skin and corneal tissues. As further explained below, using optical pulses with a judicious choice of pulse widths and duty cycles, the optical powers P.sub.1 and P.sub.2 can each be two orders of magnitude greater than their time-average values P.sub.1
and
P.sub.2
; therefore, the SNR can be greatly enhanced compared to CW SBS.
[0070]
[0071] Laser frequency control. Accurate measurement of the frequency difference of the pump and probe is required for reliable, long-term experiments. Although direct locking of two CW laser sources is possible, we can initially lock both lasers to a .sup.85Rb vapor transition line using a lock option (e.g. DLS-Pro, Toptica), and while the pump laser remains locked, the probe laser can be tuned, using a PZT built in the laser source, away from the pump frequency (lowered by ˜5.7 GHz). The beat frequency between the pump and probe can be measured using a high-speed photodiode. Tuning of the probe laser can be performed in a feedback subsystem where the electrical signal from the high-speed photodiode is mixed with a tunable local oscillator signal. The frequency of the mixed signal can represent the error signal and can be compensated by a proportional-integral-derivative controller through PZT control to match the beat frequency between the pump and probe beams to the configured local oscillator frequency. Alternatively, the beat frequency between the pump and probe beams can be measured using a high-speed photodiode and a microwave counter and used to control the tuning of the probe laser.
[0072] Filtering out the pump. We can use a high-pressure Rb vapor filter (60° C.) in front of the final photo-detector to remove back-reflected pump light and transmit the SBS-modulated probe light. The frequency lock of the pump laser can be automatically aligned the laser frequency to the filter. The effectiveness of Rb cells in Brillouin spectroscopy has been previously demonstrated. Without being bound by theory, a high extinction over >100 dB can in principle be obtained.
[0073] SNR in pulsed operation. The signal current of the detector can be expressed as i.sub.s(t)
=eη
αP.sub.2
g cos(2πf.sub.1t), where α is the optical loss in the path from the sample to the detector. Typical silicon detectors have a responsivity of 0.55 A/W and η≈3.4×10.sup.18 at 780 nm and a thermal noise current of i.sub.th≈2.2 pW/✓Hz. The total noise power is a sum of thermal, shot, and intensity noises:
i.sub.n.sup.2(t)
=(i.sub.th.sup.2+2e.sup.2η
αP.sub.2
+e.sup.2RIN η.sup.2
αP.sub.2
.sup.2, where RIN is relative intensity noise (RIN) of the probe laser. For
αP.sub.2
>>1 nW, the thermal noise is negligible compared to the other two terms. Because the probe light is pulsed, the peak power can be greater than the average power by a ratio m P.sub.2/
P.sub.2
, and RIN was increased by a factor m compared to CW operation at the same average power. For an integration bandwidth of BW [Hz], SNR=0.5 g.sup.2η
αP.sub.2
/(2 BW+mRIN η
αP.sub.2
BW). Since SBS gain g is proportional to m, g.sup.2∝m.sup.2 and SNR ∝˜m.
[0074] Balanced detection. Although this analysis already shows the SNR benefit of pulsed operation, SNR could be increased further by reducing RIN. Unlike thermal and shot noise, which are of quantum origin, RIN noise is classical noise and can be measured by tapping a small portion of the probe pulses as depicted in
The experimental laser sources can have relatively low RIN below ˜145 dB/Hz at frequencies above 100 kHz. The condition that the shot noise is equal to the RIN noise is: ηαP.sub.2
=2/RIN. For α=0.25 (6 dB loss), we get
P.sub.2
=0.74 mW. It can be shown that the optimal power is
P.sub.2
/P.sub.0=(√{square root over (9+4a)}−3)/2a, where a RIN ηαP.sub.0. For α=0.25 and P.sub.0=5 mW, we get α=13.4 and
P.sub.2
/P.sub.0=0.18. For a in a wide range of 1 to 100, the optimum
P.sub.2
/P.sub.0 ranges between 0.1 and 0.3. So, as a general guideline, we choose
P.sub.2
=0.2P.sub.0, i.e.,
P.sub.1
:
P.sub.2
=4:1.
[0075] Expected SNR. As initial parameters, we can use P.sub.1
=4 mW,
P.sub.2
)=1 mW and repetition rates of 1 MHz for the pump and 2 MHz for the probe, as shown in
[0076] This case is summarized in Table 1. To determine the peak Brillouin frequency shift, the probe laser frequency can be tuned with a sinusoidal waveform over a range of about 5.6 to 6.2 GHz (or as necessary) relative to the pump laser frequency. The photodetector signal can be digitized at a sampling rate of 10 MS/s. A Brillouin spectrum can be constructed from the data collected as a function of the pump-probe frequency difference. Prior spontaneous Brillouin systems have been limited to operating at an integration time of about 0.2 s. At the same SNR, the experimental SBS system and method can be operated with a pixel integration rate of 10 kHz, roughly 2,000 times faster than spontaneous Brillouin microscopy. Multiple spectra could be averaged to improve the frequency sensitivity at the expense of the pixel rate. The Brillouin frequency measurement error can be <10 MHz.
TABLE-US-00001 TABLE 1 Experimental operation parameters and the expected SNR gain and data acquisition speed compared to a spontaneous Brillouin microscopy system. Pump Pump Pump Probe pulse Pump peak SNR Integration Mode power power width rate power (Hz.sup.−1) time Spontaneous 5 mW — CW — 5 mW 2,500 0.2 s (exp.) Stimulated 4 mW 1 mW 5 ns .sup.(1) 1 MHz 0.8 W 5 × 10.sup.6 100 μs
Example 2
Shear-wave Stimulated Brillouin Scattering
[0077] The following experimental and prophetic remarks relate to the development and characterization of a system and method for shear-wave stimulated Brillouin scattering.
[0078] Without being bound by theory, in soft-matter physics and engineering, as well as tissue and cell biomechanics, arguably the most important mechanical property is shear modulus G or μ. Most biomaterials are viscoelastic, so both the real (elastic) and complex (viscous) components, G′ and G″, are of interest. G″ is responsible for stress relaxation of ECM, which affects cell spreading and differentiation. Shear modulus is of purely entropic origin related to the configurational free energy of polymeric chains, whereas bulk modulus is due largely to the elastic energy stored in both intramolecular and intermolecular chemical bonds. So, shear and longitudinal moduli provide different, complementary information.
[0079] The shear Brillouin frequencies of most biological samples range between ˜1 to 100 MHz, which we estimate from published values of G measured mechanically and acoustically. Unfortunately, optical filters capable of reliably discriminating MHz frequency differences are not readily available. Spontaneous shear Brillouin measurement has been applied to solid samples with GHz-scale frequency shifts but has never been demonstrated for biological samples. SBS could overcome this difficulty because the Brillouin scattered light is phase-coherent with the probe and pump beams, especially when they come from a single-frequency laser. The shear SBS signal can be detected with heterodyne detection using a standard photodiode.
[0080] For shear SBS in hard solid materials, shear electrostriction γ.sup.2 is 100 times smaller, acoustic speed V.sub.S is 500 times smaller, and linewidth Γ.sub.B is 100 times smaller than for longitudinal SBS. Thus, from Eq. (3), we find that the shear-wave SBS gain coefficient would be ˜500 times greater than longitudinal-wave SBS. For soft materials, this estimation should be corrected for three factors. First, the gain is reduced by the presence of water because liquid water does not support shear waves and it produces zero SBS gain. So, the SNR gain should be multiplied by the solid volume fraction, which is 0.2-0.3. Second, the pump and probe beams must be crossed)(˜90°, so the interaction length becomes ˜2Δx, slightly reduced from Δz by a factor of NA/n. Lastly, because of small Γ.sub.B pulsed operation does not make sense and the optimal CW pump and probe powers are close to P.sub.1=4 mW and P.sub.2=1 mW. Taking all these factors into account, the expected shear-wave SBS gain is about the same as that of CW longitudinal-wave SBS. Considering the lower optical loss in the SBS setup compared to VIPA Brillouin spectroscopy, we estimated that the acquisition speed of shear SBS microscopy would be equal to or 10-times faster than spontaneous longitudinal systems.
[0081] Recently, we developed a spontaneous shear Brillouin microscope using an off-axis confocal beam geometry and a two-stage VIPA spectrometer. The polarization state of shear-induced Brillouin light is orthogonal to the input polarization state, unlike longitudinal Brillouin scattering where the original polarization state is maintained. Therefore, polarization-based filtering discriminates shear from longitudinal Brillouin signals. The Brillouin spectra of polycarbonate, fused quartz, and borosilicate glass were recorded with an acquisition time of 1-10 s at an input power of 10 mW with a spatial resolution of ˜1 μm. The spectrometer-based system is able to resolve is >1 GHz, so it cannot be applied to soft materials with MHz-range Brillouin frequencies.
[0082] Design of a shear-wave SBS microscope: For shear SBS, we have devised an alternative scheme based on interference detection. The schematic is shown in
[0083] If the diffraction efficiency variation over the frequency tuning range of ±40 MHz causes too large power variation, we would employ a double-pass, cat's eye configuration for each FS device to reduce the intensity variation and eliminate any beam diffraction-induced deflection of the probe light after the devices.
[0084] The measurable modulus range would be about 1 Pa (f.sub.S=80 kHz) to 1 MPa (f.sub.S=80 MHz), which covers most soft biomaterials. The lower limit comes from the RIN noise increasing sharply below 1 kHz. The upper limit is due to the finite tuning range of frequency shifters. For much harder materials, shear SBS can in principle be obtained with a setup modified from the longitudinal-wave SBS microscope such that the pump and probe beams are now cross-polarized and cross-propagating; the pulsed operation provides a SNR gain compared to CW shear SBS.
The system parameters and expected performance are summarized in Table 2.
TABLE-US-00002 Pump Probe Pump & SNR Pixel Frequency Frequency Mode power power probe (Hz.sup.−1) rate resolution range Shear 4 mW 1 mW CW 2.5 × 10.sup.4 ~10 Hz 100 kHz 0.1-80 MHz SBS
Example 3
Combined Longitudinal-Shear Stimulated Brillouin Microscopy
[0085] The following experimental and prophetic remarks relate to the development and characterization of a system and method for combined longitudinal-shear stimulated Brillouin microscopy.
[0086] The measurement of both longitudinal and shear Brillouin frequency shifts simultaneously will enable us to determine the compliance matrix of a microscopic volume within a sample. Following the inventions described above, we can combine the two systems into a common sample stage. This allows the acquisition of both longitudinal and shear moduli information nearly simultaneously. The longitudinal and shear moduli provide complementary information about the polymeric content, degree of crosslinking, and water-to-solid ratio of samples. Therefore, we expect that the combined system can be used to characterize the mechanical properties of hydrogels, cells, and tissues more comprehensively than individual instruments that measure either of the two. Furthermore, the combined system will allow us to investigate the phenomenological correlation between longitudinal and shear moduli, which has been controversial.
[0087] Without being bound by theory, Hooke's law in 3D is defined as σ=Cε where σ is the stress tensor, ε is the strain tensor, and C is the compliance matrix of the material. The stiffness matrix S is defined by ε=Sσ; therefore, S=C.sup.−1, and S and C are each a 6×6 symmetric matrix with 21 independent components. In most materials with certain symmetry, the number of independent components is considerably reduced. Materials with fibrillary structures—for example, corneal stroma with organized collagen fibrils—have transverse-isotropic mirror-symmetry. This type of material belongs to the space group D.sub.∞h and has an elastic tensor in hexagonal symmetry with four independent parameters, C.sub.11, C.sub.13, C.sub.44, and C.sub.66 (∵C.sub.12=C.sub.11C.sub.66). The longitudinal acoustic wave speed, v.sub.L, and shear wave speed, v.sub.S, are directly related to these parameters. The acoustic speeds propagating along the cylindrical axis are v.sub.L,3.sup.2=C.sub.33/ρ and v.sub.S,3.sup.2=C.sub.66/ρ, where subscripts L and S represent longitudinal and shear waves, respectively, and p is the density of the material. For waves propagating in the transverse plane, v.sub.L,1,2.sup.2=C.sub.11/ρ and v.sub.S,1,2.sup.2=C.sub.44/ρ. The simplest case is that of a material that is isotropic in all directions, such as polymers and hydrogels consisting of randomly-orientated polymeric chains. The elastic tensor of such an isotropic material is described by two Lamé parameters, λ and μ: C.sub.44=μ, C.sub.12=C.sub.13=λ, C.sub.11=C.sub.33=λ+μ, and thus
Recall that Brillouin frequencies are given by: f.sub.L=nv.sub.L/λ.sub.o and f.sub.S=nv.sub.S/λ.sub.0. For most materials, n and ρ are known (˜1.41 and ˜1,050 kg/m.sup.3 respectively). Therefore, by measuring f.sub.L and f.sub.S (in three directions), the full compliance matrix can in principle be determined. For isotropic materials, a single measurement along one axis is sufficient to determine all elastic moduli. For transverse-isotropic materials, C.sub.11, C.sub.13, C.sub.44, and C.sub.66 can be computed from the measurement of longitudinal and shear Brillouin frequencies at two different sample orientations.
[0088] To determine both longitudinal and shear moduli simultaneously, we modify the longitudinal SBS setup and combine it with the shear SBS microscope. The major modification required is to re-direct the pump and probe beams for longitudinal SBS so that they both enter the sample from the below, making a 90-deg angle to one another, as illustrated in
[0089] The longitudinal SBS pulse modulation scheme is identical to that described earlier. A motorized sample stage is used as before. This allows the probe and sample beams in both longitudinal and shear SBS measurements to be fixed, ensuring stable foci overlap. The pump and probe beams for longitudinal SBS are aligned to the same paths as the pump and probe beams for shear SBS, such that the system can switch between the two modalities simply by using flip mirrors. This combined system provides an unprecedented opportunity to non-invasively measure the compliance matrix at each location within a sample with microscopic resolution.
Example 4
[0090] In order to validate the effectiveness of the longitudinal stimulated Brillouin scattering systems and methods described herein, an exemplary experiment was conducted.
[0091] The goal of this project is to build a high speed subcellular stimulated Brillouin microscope for mapping the elastic properties of cells and tissues. Imaging the elasticity of cells is important for understanding how cells respond to the stiffness of the extracellular matrices. Stiffness has been found to affect wound healing and metastasis, further investigation on this effect requires high resolution imaging instrument that could capture the dynamics of elasticity in cells. The existing Brillouin microscope had been used to study elastic properties of cells. However, the spontaneous scattering process by which it operates has small response signal thus requires a long integration time. In this project, we overlap two laser beams at the sample to stimulate Brillouin scattering. The two lasers are configured to have frequency separation that matches the acoustic vibration of the sample, thus enhances the generation of phonons and allows faster imaging speed. Stimulated Brillouin scattering (SBS) requires high laser power for the enhanced signal, to maintain a safe laser exposure for cells, we use pulse scheme to reduce the average power-. Depending on the configuration of the lasers, either the longitudinal or shear acoustic waves can be stimulated. We began by building the Brillouin microscope system for the longitudinal wave, optimizing the detection signal level. Subsequently, we aim to alter the microscope system to include shear wave measurement, followed by testing the microscope system with various materials.
Experimental Results
[0092] The primary components of a stimulated Brillouin microscope system for longitudinal waves have been constructed. The system includes two lasers, configured to counter-propagate and overlap each other in the sample to stimulate Brillouin scatterings. The schematics of the microscope system are shown in
[0093] Each of the pump and probe laser beam passes through an electro-optic modulator (EOM), the laser beam experiences a π-phase shift through varying the Pockels cell's birefringence with high voltage fast switching. Then the laser beam passes through a polarizer, forming optical pulses that retain the peak power of the incoming continuous wave light. The EOMs are controlled by the respective EOM driver and an EOM controller that configures and synchronizes the pulse timing for both pump and probe lasers. In a preliminary test, the pulse periods were set to 2.5 μs and 2.0 μs for pump and probe lasers, respectively. For calibration purposes, the duty cycles for both laser pulses were set to 50% to enhance the signal level. Each laser beam passes through a quarter wave plate, focuses down to a spot of ≈1.0 μm diameter by a microscope objective and spatially overlaps with each other to stimulate Brillouin scattering in the sample. The probe laser carrying Brillouin signal experiences half wave shift when passing through the two quarter wave plates. Having half wave difference with the pump beam, the probe beam is directed toward the detection path by the polarizing beam splitter.
[0094] The detection path consists of a custom-made 150 mm Rubidium vapor cell filter heated up to 70° C. to suppress the pump reflection up to 50 dB. The probe laser with the Brillouin signal was detected by a photoreceiver (Newfocusmodel 2031). The detected signal was fed into a lock-in amplifier where the reference signal (pulses with repetition rate corresponds to the frequency difference between the two laser pulses) was generated by the EOM controller. The integration time of the lock-in amplifier was set to 10 ms, the maximum sensitivity level was set to 10 mV, the reference signal has a frequency of 100 kHz (The frequency difference between pump and probe pulses). The output signal of the lock-in amplifier was acquired by a digital acquisition board and processed by a computer. In the data acquisition process, the local oscillator frequency was stepping sequentially at 1 MHz interval over the frequency range of interest. The reading from lock-in amplifier was recorded for each frequency step, at 1 kHz rate.
[0095] The peak pump and probe power measured at the sample was 0.8 W and 32 mW, respectively. The microscope objective has NA=0.6, the calculated laser focusing spot size is 2w≈1.0 um and the Rayleigh length 2z≈1.4 μm. The effective beam area is A ≈0.5 μm.sup.2. Given pump power of 0.8 W, the pump intensity is Ip≈1500 GW/m.sup.2. Consider ethanol, which has Brillouin gain factor of g=0.12 (m/GW), Brillouin gain is calculated as G=g I.sub.p 2z≈2.5×10.sup.−4. With the probe power measured at the detection ˜0.5 mW, and the conversion gain of the photoreceiver given as 1.2×10.sup.3, the expected signal strength is ˜150 μV. In the case of water, the Brillouin gain factor is 0.048 (m/GW) and the expected signal strength is 60 μV. We tested the stimulated Brillouin microscope with water and ethanol contained inside a shallow well of 1 mm depth with a transparent glass base, the Brillouin gain spectra obtained experimentally are shown in
[0096] While the above detailed description has shown, described, and pointed out novel features as applied to various embodiments, it will be understood that various omissions, substitutions, and changes in the form and details of the devices or algorithms illustrated can be made without departing from the spirit of the disclosure. As will be recognized, certain embodiments of the disclosures described herein can be embodied within a form that does not provide all of the features and benefits set forth herein, as some features can be used or practiced separately from others. The scope of certain disclosures disclosed herein is indicated by the appended claims rather than by the foregoing description. All changes which come within the meaning and range of equivalency of the claims are to be embraced within their scope.