*in vivo*histological guidance during resection of brain tumors. The microscope utilizes a dual-axis confocal architecture to efficiently reject out-of-focus light for high-contrast optical sectioning. A biaxial microelectromechanical system (MEMS) scanning mirror is actuated at resonance along each axis to achieve a large field of view with low-voltage waveforms. The unstable Lissajous scan, which results from actuating the orthogonal axes of the MEMS mirror at highly disparate resonance frequencies, is optimized to fully sample 500×500 pixels at two frames per second. Optically sectioned fluorescence images of brain tissues are obtained in living mice to demonstrate the utility of this microscope for image-guided resections.

Surgical resection is a major component in the standard of care for the treatment of brain cancers, where clinical studies have demonstrated that improved outcomes correlate with more complete resections, especially in conjunction with adjuvant therapies such as radiochemotherapy.^{1, 2} Histological guidance provided by a miniaturized microscope could be used to improve the accuracy of resection at selected sites where tumor margins are difficult to identify, thereby improving survival and reducing neurological damage. We have been developing dual-axis confocal (DAC) microscopes for *in vivo* imaging with cellular resolution.^{3, 4} Unfortunately, the
$10\text{-}\mathrm{mm}$
-diam form factor of our current microscopes greatly exceeds their field of view (FOV) of less than
$0.5\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$
. The size of these instruments also surpasses that of conventional tools for tumor resection, which are typically
$1\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}2\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$
in diameter, rendering such a device inconvenient for guiding surgical procedures.

Here we demonstrate that the incorporation of a gradient-index (GRIN) lens at the tip of a microelectromechanical system (MEMS)-scanned DAC allows the image plane to be relayed forward by
$3\phantom{\rule{0.3em}{0ex}}\mathrm{cm}$
through a
$1.8\text{-}\mathrm{mm}$
-diam optic (see Fig. 1
). This form factor is designed to enable convenient microscopy in a clinical setting for optical biopsy and image-guided resection of questionable tumor-margin locations using fluorescent contrast agents. Previously, GRIN relay lenses have been utilized for conventional single-axis confocal^{5} and nonlinear microscopes.^{6} Here we show that two low-NA beams, with a crossing angle of
$\sim 45\phantom{\rule{0.3em}{0ex}}\mathrm{deg}$
, can also be relayed and scanned through such an optical element with minimal aberrations.

Studies have demonstrated the ability of a DAC to effectively reject out-of-focus scattered light for high-contrast optical sectioning in turbid media.^{7, 8, 9} The DAC utilizes low-power diode lasers along with simple low-NA optics to achieve deep tissue imaging. The long working distance afforded by the low-NA focused beams allows for postobjective beam scanning with a biaxial MEMS scanner.^{3, 4} This eliminates the need for compound lenses to mitigate the off-axis aberrations that result from scanning a beam through a focusing objective.

Our MEMS scanning device contains a total of four electrostatic actuators: a pair of actuators for each axis of the mirror to achieve tilt in opposing directions (pull-pull configuration).^{10} We have previously achieved raster-scanned imaging with these MEMS-scanned microscopes by actuating one axis of the mirror at its resonance frequency
$(\sim 1\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}3\phantom{\rule{0.3em}{0ex}}\mathrm{kHz})$
and actuating the orthogonal axis with a slow sawtooth waveform at a frame rate of
$4\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}30\phantom{\rule{0.3em}{0ex}}\mathrm{Hz}$
. The fast axis, because it is driven at resonance, can achieve large deflection angles (
$\pm 3\phantom{\rule{0.3em}{0ex}}\mathrm{deg}$
or more) with relatively low-voltage waveforms
$(\sim 20\phantom{\rule{0.3em}{0ex}}\mathrm{V})$
. Furthermore, only one actuator is necessary to drive an axis at its natural resonance frequency. On the other hand, relatively large voltages
$(\sim 100\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}120\phantom{\rule{0.3em}{0ex}}\mathrm{V})$
are applied to two actuators to achieve modest deflection angles
$(\pm 1\u20133\phantom{\rule{0.3em}{0ex}}\mathrm{deg})$
at off-resonance frequencies. Therefore, the ability to drive *both* axes of the mirror at resonance, in a Lissajous scanning pattern,^{11} enables large deflection angles to be achieved with fewer wired connections and lower voltages.

Previous implementations of Lissajous scanning have generally utilized stable scan patterns that self-repeat at a frame rate determined by the frequency ratio of the driving waveforms along each axis.^{12, 13} However, with arbitrary resonance frequencies, it is difficult to achieve a stable Lissajous scan pattern that repeats a rate of
$1\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}10\phantom{\rule{0.3em}{0ex}}\mathrm{Hz}$
and also provides coverage over a large FOV of
$500\times 500\phantom{\rule{0.3em}{0ex}}\text{pixels}$
. Here we demonstrate a general method to utilize an unstable (nonrepeating) “sliding” Lissajous pattern to fully sample
$500\times 500\phantom{\rule{0.3em}{0ex}}\text{pixels}$
at a frame rate of
$2\phantom{\rule{0.3em}{0ex}}\mathrm{Hz}$
. A similar strategy utilizing an unstable Lissajous scan of a MEMS mirror was demonstrated previously for two-photon microscopy.^{14} In that work, high frame rates of up to
$10\phantom{\rule{0.3em}{0ex}}\mathrm{Hz}$
were achieved with a FOV of
$256\times 256\phantom{\rule{0.3em}{0ex}}\text{pixels}$
. In our particular implementation, the scan frequencies are tuned to yield a Lissajous pattern that predictably slides at an optimized rate to ensure that each pixel is sampled. This eliminates the need for filtering or averaging to account for unsampled pixels.

Figure 1 illustrates the optical circuit of our microscope. This diagram includes a ray-trace simulation (Zemax, Bellevue, Washington) depicting the beam paths of the DAC as they are relayed through a nonmagnifying (1:1) GRIN triplet relay assembly consisting of a pair of ${\scriptstyle \frac{1}{4}}$ -pitch, $4.25\text{-}\mathrm{mm}$ -long, high-NA (0.51) GRIN lenses on either side of a ${\scriptstyle \frac{1}{2}}$ -pitch, $21.85\text{-}\mathrm{mm}$ -long, low-NA (0.19) GRIN lens. Figures 1 and 1 are design drawings of the basic components of the microscope. The parabolic reflector (AgilTek) insures that the two collimated beams intersect and focus to the same point, which defines the focal volume of the microscope. A pair of $0.1\text{-}\mathrm{deg}$ wedged alignment (Risley) prisms, shown in Fig. 1, are rotated to bring the collimated beams into parallel alignment with each other. Figure 1 is a photograph of a fully packaged microscope, including the piezoelectric actuator at the proximal end that translates the MEMS mirror in the axial direction for adjusting the focal depth.

The single-mode fiber-pigtailed collimators (GRINTECH Gmbh, Jena, Germany, focal
$\text{length}=4.5\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$
) used in our microscope yield a
$1\u2215{\mathrm{e}}^{2}$
beam diameter of
$0.9\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$
. As a result of the parabolic reflector having a
$6\text{-}\mathrm{mm}$
focal length and the collimators being spaced apart by
$4.7\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$
, the focusing half angle of each beam
$(1\u2215{\mathrm{e}}^{2})$
and the intersection half angle between the beams, respectively, are
$\alpha =0.075\phantom{\rule{0.3em}{0ex}}\mathrm{rad}$
and
$\theta =22.2\phantom{\rule{0.3em}{0ex}}\mathrm{deg}$
[see Fig. 1]. If the Gaussian illumination and collection beams are truncated (apodized) to allow 99% power transmission, diffraction theory predicts a full-width at half-maximum (FWHM) spatial resolution (in micrometers) of:^{7}

## Eq. 1

$$\Delta {x}_{d}=\frac{0.466\lambda}{n(\pi /2\cdot \alpha )\mathrm{cos}\phantom{\rule{0.2em}{0ex}}\theta}=2.1,\phantom{\rule{1em}{0ex}}\Delta {y}_{d}=\frac{0.466\lambda}{n(\pi /2\cdot \alpha )}=1.9,\phantom{\rule{1em}{0ex}}\Delta {z}_{d}=\frac{0.466\lambda}{n(\pi /2\cdot \alpha )\mathrm{sin}\phantom{\rule{0.2em}{0ex}}\theta}=5.1,$$As mentioned previously, an unstable Lissajous scan pattern is achieved by driving our MEMS mirror at resonance along both axes, enabling large-FOV imaging with relatively low-voltage waveforms. The resulting scan pattern is described by the following parametric equations as a function of time:

where ${\varphi}_{x}$ and ${\varphi}_{y}$ are phase shift constants (in radians), $a$ and $b$ are the slow and fast frequencies, respectively, in cycles/sec (Hz), and the FOV of the scan extends from 0 to $2A$ in the $x$ direction and 0 to $2B$ in the $y$ direction (in this case, $A=B=250$ ). The resonant frequencies of the two axes of our MEMS mirror are $985\pm 10\phantom{\rule{0.3em}{0ex}}\mathrm{Hz}$ and $2220\pm 10\phantom{\rule{0.3em}{0ex}}\mathrm{Hz}$ (FWHM bandwidth $\sim 20\phantom{\rule{0.3em}{0ex}}\mathrm{Hz}$ ).Figure 2 depicts the trajectory of a Lissajous scan pattern between $t=0$ and $4\phantom{\rule{0.3em}{0ex}}\mathrm{ms}$ (assuming ${\varphi}_{x}={\varphi}_{y}=0$ ). If the frequency ratio of the fast- and slow-axis waveforms is exactly $9\u22154$ , then a stable Lissajous pattern results, repeating once for every four cycles of the slow axis or every nine cycles of the fast axis ( $4\u2215a$ or $9\u2215b\phantom{\rule{0.3em}{0ex}}\mathrm{sec}$ ). If the ratio between the frequencies is slightly different than $9\u22154$ , an unstable Lissajous pattern is generated that appears to slide in time in the horizontal direction [Fig. 2]. The sliding speed, determined by the frequency ratio, is optimally chosen such that the Lissajous pattern shifts by a distance of less than $1\phantom{\rule{0.3em}{0ex}}\text{pixel}$ for every nearly repeated cycle ( $4\u2215a$ or $9\u2215b\phantom{\rule{0.3em}{0ex}}\mathrm{sec}$ ), ensuring that no pixels are missed over time. Figures 2 and 2 illustrate how an entire image is sampled over time by the sliding Lissajous pattern. Taking the derivative of Eq. 3 with respect to time reveals that the maximum velocity of the fast-axis scan is given by $2\pi \cdot b\cdot B\cong 3.5{\left(10\right)}^{6}\phantom{\rule{0.3em}{0ex}}\text{pixels}\u2215\mathrm{sec}$ . Thus, a data acquisition rate of $4\phantom{\rule{0.3em}{0ex}}\mathrm{MHz}$ is chosen to obviate dead (unsampled) pixels along the trajectory of the Lissajous scan.

Sinusoidal waveforms are generated with National Instruments (Austin, Texas) boards programmed in LabVIEW and sent through a voltage amplifier (Agiloptics, Albuquerque, New Mexico) to drive our biaxial MEMS mirror. A National Instruments digitizer board is used to record signals from a PMT detector at a rate of $4\phantom{\rule{0.3em}{0ex}}\mathrm{MHz}$ . All processes are synchronized for phase-locked operation, and the phase shifts ${\varphi}_{x}$ and ${\varphi}_{y}$ are determined via calibration. In our implementation, $a=986.0\phantom{\rule{0.3em}{0ex}}\mathrm{Hz}$ and $b=2219.0\phantom{\rule{0.3em}{0ex}}\mathrm{Hz}$ . For each detector value sampled at time $t$ , Eqs. 2, 3 are used to calculate the exact coordinates $x$ and $y$ of the Lissajous scan. These coordinates are rounded to the nearest integer to locate the most appropriate pixel for data entry within the $500\times 500\phantom{\rule{0.3em}{0ex}}\text{pixel}$ bitmap image (programmed in LabVIEW). Image degradation due to rounding errors is minimal due to the fact that we oversample the image. The spatial dimension of our target FOV is approximately $500\times 500\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$ . Therefore, with a pixel density of $500\times 500$ , and a lateral resolution of $\sim 4\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$ , our sampling $(1\phantom{\rule{0.3em}{0ex}}\text{pixel}\u2215\mu \mathrm{m})$ exceeds the Nyquist criterion by a factor of 2. Two million points are collected per image, resulting in a frame rate of $2\phantom{\rule{0.3em}{0ex}}\mathrm{Hz}$ . Redundantly sampled pixels contain an average of the sampled values at that pixel location.

Example images acquired with our GRIN-relayed DAC are shown in Figs.
3, 3, 3, 3, 3, 3
. Figure 3 is a reflectance image of a United States Air Force (USAF) bar target (group 7), in which the width of the smallest set of lines is
$2.2\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
per line. Due to the field curvature of the microscope, only the center of the FOV is in focus for this flat reflective target. In practice, field curvature is not a major concern, since our device is designed to image at least ten microns beneath the surface of a thick tissue, rather than imaging a flat surface. Figure 3 is a depth projection
$(\mathrm{FOV}=400\times 400\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m})$
of
$10\text{-}\mu \mathrm{m}$
-diam fluorescent beads embedded in Matrigel. The color map encodes depths from
$0\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}250\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
. Figure 3 is an *in vivo* image at a depth of
$50\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
$(\mathrm{FOV}=325\times 325\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m})$
of vasculature within the cerebrum of a mouse after a portion of the skull has been surgically removed (five-frame average). To highlight the vasculature in this image, fluorescein-dextran
$(\mathrm{MW}=\mathrm{70,000})$
was injected into the retro-orbital capillary plexus at a dose of
$\sim 1\phantom{\rule{0.3em}{0ex}}\mathrm{mg}$
(
$40\text{-}\mathrm{mg}\u2215\mathrm{kg}$
body weight) in a
$200\text{-}\mu \mathrm{L}$
volume. Figure 3 is of the same tissue as Fig. 3, except that it was acquired at a depth of
$100\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
(five-frame average). Figure 3 is an *ex vivo* image at a depth of
$25\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
$(\mathrm{FOV}=350\times 350\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m})$
of cerebellum from a transgenic mouse (
$10\phantom{\rule{0.3em}{0ex}}\text{days}$
of age) that expresses green fluorescent protein (GFP) under the Math1 promoter.^{15} After weaning, expression of this promoter is repressed in the cerebellum, except in medulloblastomas, a cerebellar brain tumor that spontaneously develops in this gene knock-out mouse (
$Pt{c}^{+\u2215-}p{53}^{-\u2215-}$
Math1-GFP). In Fig. 3, *in vivo* imaging of the tumor is performed after surgical removal of a portion of the skull. This image (
$\mathrm{FOV}=350\times 350\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
;
$\text{depth}=30\phantom{\rule{0.3em}{0ex}}\mu \mathrm{m}$
) reveals the GFP-expressing tumor in contrast to the dark normal brain tissue, and indicates that this instrument would have utility in detecting tumor margins when used in conjunction with appropriate exogenously applied fluorescent molecular probes. Figures 3 and 3 are single frames that have been filtered with a smoothing algorithm that replaces each pixel with an average of its adjacent pixels. Spatial resolution is minimally affected, since the image is oversampled.

We have presented the design and demonstration of a miniature DAC outfitted with a
$1.8\text{-}\mathrm{mm}$
-diam GRIN relay lens for convenient macroscopic positioning of the device at surgical sites. Our previous designs implemented raster-scanned MEMS mirrors, which were limited in their FOV because of the modest scan range attained by the nonresonant slow axis of the MEMS device. In addition, the raster-scanned MEMS mirrors required higher operating voltages that may not be preferred in clinical applications. Here we have described the use of resonant scanning along both axes of the MEMS mirror to achieve a large FOV with low actuation voltages and reduced numbers of actuators. We have developed the calibrated and synchronized instrumentation to decode the nonrepeating Lissajous pattern that results, and have demonstrated *in vivo* imaging of normal brain and tumor in living mice. Our goal is to use the GRIN-relayed DAC technology, along with exogenous fluorescent contrast agents, for spatially precise and molecularly specific resection of brain tumors

## Acknowledgments

This work was funded, in part, by grants from the National Institutes of Health, including U54 CA105296 and U54 CA136465 (NCI). These studies were also supported by the Stanford Center for Children’s Brain Tumors (CCBT) at the Lucile Packard Children’s Hospital. Liu was supported, in part, by a Canary Foundation/American Cancer Society postdoctoral fellowship for early cancer detection, and by an NIH Pathway to Independence Award, K99 EB008557 (NIBIB). We thank Ellis Garai for technical support, and Helen Rayburn and Matthew Scott for providing breeders for our medulloblastoma mouse colonies.

## References

*in vivo*imaging by a handheld dual-axes confocal microscope,” Opt. Express, 16 7224 (2008). https://doi.org/10.1364/OE.16.007224 1094-4087 Google Scholar

*In vivo*mammalian brain imaging using one-and two-photon fluorescence microendoscopy,” J. Neurophysiol., 92 3121 (2004). https://doi.org/10.1152/jn.00234.2004 0022-3077 Google Scholar

*in vivo*microscopy,” J. Microelectromech. Syst., 16 969 (2007). https://doi.org/10.1109/JMEMS.2007.892900 1057-7157 Google Scholar

*ARF*accelerates medulloblastoma in mice heterozygous for

*patched*,” Cancer Res., 61 513 –516 (2001). 0008-5472 Google Scholar