Document Detail

Dispersion-cancelled biological imaging with quantum-inspired interferometry.
Jump to Full Text
MedLine Citation:
PMID:  23545597     Owner:  NLM     Status:  In-Data-Review    
Quantum information science promises transformative impact over a range of key technologies in computing, communication, and sensing. A prominent example uses entangled photons to overcome the resolution-degrading effects of dispersion in the medical-imaging technology, optical coherence tomography. The quantum solution introduces new challenges: inherently low signal and artifacts, additional unwanted signal features. It has recently been shown that entanglement is not a requirement for automatic dispersion cancellation. Such classical techniques could solve the low-signal problem, however they all still suffer from artifacts. Here, we introduce a method of chirped-pulse interferometry based on shaped laser pulses, and use it to produce artifact-free, high-resolution, dispersion-cancelled images of the internal structure of a biological sample. Our work fulfills one of the promises of quantum technologies: automatic-dispersion-cancellation interferometry in biomedical imaging. It also shows how subtle differences between a quantum technique and its classical analogue may have unforeseen, yet beneficial, consequences.
M D Mazurek; K M Schreiter; R Prevedel; R Kaltenbaek; K J Resch
Related Documents :
24901977 - A new dusts sensor for cultural heritage applications based on image processing.
23390117 - Mid-hindbrain malformations due to drugs taken during pregnancy.
23202167 - A coded aperture compressive imaging array and its visual detection and tracking algori...
21276917 - In-vivo imaging of femoral artery nitinol stents for deformation analysis.
22132507 - Visual chimaeras obtained with the riesz transform.
16018247 - Controlling the occlusion of single restorations with cad/cam technology.
Publication Detail:
Type:  Journal Article    
Journal Detail:
Title:  Scientific reports     Volume:  3     ISSN:  2045-2322     ISO Abbreviation:  Sci Rep     Publication Date:  2013 Apr 
Date Detail:
Created Date:  2013-04-02     Completed Date:  -     Revised Date:  -    
Medline Journal Info:
Nlm Unique ID:  101563288     Medline TA:  Sci Rep     Country:  England    
Other Details:
Languages:  eng     Pagination:  1582     Citation Subset:  IM    
1] Institute for Quantum Computing and Department of Physics & Astronomy, University of Waterloo, Waterloo, ON N2L 3G1 CANADA [2].
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Sci Rep
Journal ID (iso-abbrev): Sci Rep
ISSN: 2045-2322
Publisher: Nature Publishing Group
Article Information
Download PDF
Copyright © 2013, Macmillan Publishers Limited. All rights reserved
Received Day: 04 Month: 02 Year: 2013
Accepted Day: 18 Month: 03 Year: 2013
Electronic publication date: Day: 02 Month: 04 Year: 2013
collection publication date: Year: 2013
Volume: 3E-location ID: 1582
PubMed Id: 23545597
ID: 3613801
Publisher Item Identifier: srep01582
DOI: 10.1038/srep01582

Dispersion-cancelled biological imaging with quantum-inspired interferometry
M. D. Mazurek14
K. M. Schreiter14
R. Prevedel12
R. Kaltenbaek13
K. J. Rescha1
1Institute for Quantum Computing and Department of Physics & Astronomy, University of Waterloo, Waterloo, ON N2L 3G1 CANADA
2Research Institute of Molecular Pathology (IMP) & Max F. Perutz Laboratories, University of Vienna, Dr. Bohr Gasse 7-9, 1030 Vienna, Austria
3Vienna Center for Quantum Science and Technology, Faculty of Physics, University of Vienna, Vienna, Austria
4These authors contributed equally to this work.

Quantum information science promises powerful and unconventional capabilities across a broad range of technologies. An important example relates to the imaging technology, optical coherence tomography (OCT). OCT can noninvasively reconstruct the 3-dimensional structure of tissue with micron resolution1; it is emerging as an important clinical tool with diverse medical applications. OCT can diagnose retinal diseases such as glaucoma, analyze artherosclerotic tissues within arteries, and detect early-stage cancerous lesions in breast tissue2, 3. In addition, OCT has found application in precision laser machining4. Since OCT relies on low-coherence interferometry, its axial resolution is limited by the coherence length of the light, inversely proportional to the bandwidth. The coherence length determines ultimate resolution, however, material dispersion can limit the practical one.

Fundamental studies in quantum optics showed that interference with energy-time entangled photon pairs5 exhibits inherent robustness against unbalanced dispersion6, 7. Even-order effects of dispersion, including the dominant group-velocity dispersion, are automatically cancelled, effectively solving the dispersion problem. This dispersion cancellation is automatic since one does not need to precisely measure and compensate the dispersion. When energy-time entanglement is strong, but not perfect, the effect is more accurately described as automatic even-order dispersion reduction since the dispersion is dramatically reduced, not cancelled8. OCT based on entangled-photon interferometry was proposed to harness automatic dispersion cancellation9. This quantum-optical coherence tomography (QOCT) has two significant barriers to practical implementation. Firstly, the reliance on producing and resolving individual photon pairs places stringent limits on the QOCT signal using state-of-the-art systems10, 11. Secondly, in samples with multiple interfaces, QOCT produces a signal for each interface and an additional, artifact feature for each unique pair of interfaces9, 12, 13; the number of artifacts grows quadratically with the number of interfaces, cluttering the image of complex samples. QOCT was first used to perform an axial scan of a coverslip13 and later applied to measure the surface topography of a gold-coated onion sample14. The gold coating was essential to increase reflectivity but rendered the technique impractical for in vivo applications; furthermore, it prevented imaging the sample's internal structure, which is one of the main benefits of OCT and necessary for most medical applications.

Recently, several different approaches have shown that dispersion cancellation does not require entanglement, but can also be observed in classical systems15, 16, 17, 18, 19. While all of these methods could, in principle, solve the low-signal problem of QOCT, each suffers from unwanted artifacts. Here, we focus on one of these techniques, chirped-pulse interferometry (CPI). A method for identifying artifacts in CPI has been demonstrated, but it requires multiple axial scans of a sample and is thus inherently slow25.

In the present work, we describe and demonstrate a new method for CPI using a single beam of shaped laser pulses. Our method produces background-free, dispersion-cancelled signals, completely free of artifacts without the need for multiple scans. We apply this technique to image a biological sample, demonstrating dispersion cancellation and observing the sample's internal structure. CPI overcomes both limitations of QOCT while retaining its advantages, demonstrating its potential for future practical application.

Theoretical description

CPI uses classical light with strong frequency anti-correlations. To create these anti-correlations, we constructed a 4-F pulse shaper with a spatial light modulator (SLM)20, 21. In CPI, this method has distinct advantages over pulse-stretching techniques with bulk optics22, 23, including: straightforward optimization of the chirp parameter8, better stability and efficiency, and more complex pulse shapes. We apply a frequency-dependent phase shift to the laser pulses, ϕ(ω) = −A(ωω0)|ωω0|, where A is a positive constant. The absolute value distinguishes this from the quadratic phase leading to linear chirp: ϕ(ω) applies a linear chirp to red-shifted frequencies (ω < ω0) and an equal, opposite chirp (antichirp) to blue-shifted frequencies (ω > ω0). The resulting pulse has frequency ω0 at its lagging edge, and instantaneous frequencies in the preceding part of the pulse obey the function . These are the frequency anti-correlations needed for CPI.We refer to this as a Blue-Antichirped-Red-Chirped (BARC) pulse.

CPI can be understood by considering the schematic in Figure 1a. Light in the upper arm of the interferometer travels through a dispersive material of length L as well as a distance L1L through free space; light in the lower arm travels a distance L2 through free space. The dispersive material in the upper arm has a wavevector of light that can be expanded about frequency ω0 as k(ω) = k(ω0) + α(ωω0) + β(ωω0)2 + …, where α and β describe the group delay and group velocity dispersion, respectively. At any time, two laser beams with frequencies ω0 + Δ and ω0 − Δ enter the interferometer with corresponding amplitudes E(Δ) and E(−Δ). After travelling through the interferometer they overlap at the nonlinear crystal for sum-frequency generation (SFG). Two paths produce SFG with frequency 2ω0; either blue-shifted light travels the upper arm and red-shifted light travels the lower arm, or vice versa. The amplitudes for these two paths interfere to give a signal S(Δ, τ) = |E(Δ)E(−Δ)|2 (1 + cos [ϕ+(Δ, τ) − ϕ(Δ, τ)]), where τ = (L2L1 + L)/c is the time delay between the two paths18. To second order in the wave-vector and ignoring a global phase, the respective phases of the two paths are, . The final signal results from integrating Δ over the pulse bandwidth, forming a peak at τ = αL26. Since unbalanced dispersion contributes the same phase βΔ2 to each path, the effect cancels out of the final signal; this is automatic dispersion-cancellation. With imperfect anti-correlations, dispersion cancellation persists if the unbalanced dispersion is much less than the chirp parameter, A (see ref. 8).

Ideally, each interference signal would correspond to an interface in the sample. However, in both QOCT13 and CPI25 an additional artifact signal appears halfway between the signals arising from each pair of interfaces. In complex samples, artifacts can outnumber the features from real interfaces, seriously impeding reconstruction and interpretation. The origin of the artifacts in these two techniques are subtly, but importantly different (see Supplementary Fig. S1 online). Artifacts in CPI are a result of interference in the SFG light at frequencies blue-shifted or red-shifted from the operating frequency, 2ω0, by an amount Δω = Δτ/(4A), where Δτ is the time-delay difference between the two sample interfaces. Spectral filtering of the SFG can remove all artifacts arising from pairs of interfaces separated by more than some minimum delay. The analogous method of artifact removal in QOCT requires coincidence detection with tens of femtoseconds time-resolution, which is extremely difficult in practice.

Experimental setup and characterization

Our experimental setup is shown in Figure 1b. We used this system to image a stack of two microscope coverglass slips (Figure 1c) and an onion (Figure 1d). First, we focus on the coverglass sample to benchmark our system performance. We measure the SFG power as a function of time delay, τ, with two types of pulses, a transform-limited (TL) pulse and the BARC pulse. Using the TL pulse is equivalent to OCT with background-free autocorrelation24, which does not suffer from artifacts but is not dispersion cancelling. In order to remove artifacts from signals acquired with BARC pulses, the SFG signal is sent through a 0.35-nm-bandwidth filter. See Methods for more details.

The data measured using TL and BARC pulses are shown in Figures 2a and b, respectively. Black (red) lines show data without (with) unbalanced dispersion from 3-mm BK7 in the sample arm, where the light passes twice through the glass. The left-most four peaks correspond to the front and back surfaces of the first and second coverglass pieces. The average delay between the first (second) two peaks is 248 ± 4 μm (236 ± 4 μm). Dividing by the group index ng = 1.517 of the coverglass gives thicknesses of 163 μm and 156 μm for the two slides, in good agreement with 164 ± 3 μm and 157 ± 3 μm measured with a caliper. We subtracted a constant 1.6 mm from the delay-arm motor position in the unbalanced-dispersion data to compensate for the group delay from the additional glass so corresponding peaks could be overlayed for comparison.

Figures 2a and b show that no additional artifact features arise going from a TL to a BARC pulse. In the online Supplementary Information, we show the maximum layer separation giving rise to artifacts is 3 μm, narrower than the peak widths, 4.2 μm. As expected, any resolvable artifacts are filtered out of our signal.

The average signal peak width for the TL pulses is 3.6 ± 0.2 μm (FWHM) which broadened by 106% to 7.4 ± 0.6 μm by the dispersion. Uncertainties represent the standard deviations of the five peak widths. In contrast, the BARC-pulse produced average signal widths of 4.2 ± 0.1 μm that broadened by only 7% to 4.5 ± 0.2 μm, demonstrating dispersion cancellation. The measured peak widths for the TL and BARC pulses are in good agreement with the theoretical calculations of 3.4 μm and 4.1 μm, respectively, using the 60 nm acceptance bandwidth of our system. The peak width of the BARC-pulse signal is broadened compared to the TL pulses as a result of the narrow filtering of the SFG; if a broader bandwidth was measured instead, the widths would become equal, but artifacts can reappear. The dispersion cancellation observed cannot be explained by the slightly broader BARC-pulse signal in the balanced dispersion case. If one used a TL pulse yielding a 4.2 μm peak width in the absence of dispersion, equal to the width to our BARC-pulse signal, we calculate that two passes through 3-mm BK7 glass would broaden the signal to 5.8 μm, a 38% increase. Our BARC-pulse signal is broadened by 7%, thus the dispersion cancellation is significant even when compared with this conservative benchmark.

Dispersion-cancelled biological imaging

We prepared a sample of onion as depicted in Figure 1d and took a set of axial (z) scans, moving the sample in the y-direction between scans. The data are displayed in Figure 3, where the four panels a, b, d, and e show cases with and without dispersion (3-mm BK7) for both TL and BARC pulses. The vertical axes are the delay-arm motor positions and the horizontal axes show the transverse y-positions. The images show the cellular structure of the onion deep into the sample and are artifact-free. Beside each image, we show a single axial scan taken at the y-position marked by the red line in each plot. The TL-pulse signal peaks are dramatically broadened when unbalanced dispersion is added, but the BARC-pulse peak widths are unchanged, directly demonstrating automatic dispersion cancellation in our image. In order to compare the effect of dispersion over the entire images produced by the TL pulses and BARC pulses, we extracted the widths of signal peaks throughout the images and display corresponding histograms in Figures 3c and f. The TL-pulse peaks broadened by 61% from an average of 3.6 ± 0.5 μm to 5.7 ± 1.2 μm, where uncertainties are the standard deviations of each distribution. The BARC-pulse peak widths increased by only 4% from 4.2 ± 0.8 μm to 4.3 ± 0.6 μm, less than the standard deviation. Hence, dispersion is cancelled througout the BARC-pulse images.

To further show the capability of our method, we took a set of axial (z) scans over a grid of x- and y-positions of a different onion sample prepared in the same manner as before, using the BARC pulses with no unbalanced dispersion. From this data, we extracted 2D cross-sections of the sample shown in Figures 4a–c and the 3D cell wall structure of the top layer of cells is shown in Figure 4d. Thus CPI, with its inherent dispersion cancellation, is suitable for practical 3D imaging of biological samples.


We have demonstrated dispersion-cancelled, artifact-free, optical-coherence-tomography imaging of a biological sample. For future work, incorporating a nonlinear material with larger nonlinearity and acceptance bandwidth24, 27 will improve the system, increasing acquisition rates and image resolution. Dispersion-cancelled OCT with chirped-pulse interferometry draws upon insights from quantum information science. Exploiting subtle differences in the analogous roles played by different physical parameters between the techniques allowed problems inherent to the quantum scheme to be solved in the classical technology; extraction of artifacts is very hard, or even technologically impossible, in the quantum device, yet becomes straightforward in CPI. Our results remove the technological barriers to dispersion-cancelled biological imaging and underscore the importance of understanding classical analogues to quantum mechanical effects.


Pulses from a titanium:sapphire laser (808 nm, 90 nm FWHM) pass through a 4-F pulse shaper incorporating a spatial light modulator (CRi SLM-640-D-VN)20, 21, also see ref. 28 for details. The SLM served two purposes, compressing the pulses by compensating for balanced dispersion in the setup, and applying the BARC-phase ϕ(ω) = −A(ωω0)|ωω0|, where A = 2500 fs2 and λ0 = 2πc/ω0 = 809.60 nm. The shaped pulses were split on a beamsplitter with 16 mW sent to the sample and 24 mW into a variable-delay line (a retroreflector on a motorized stage). A 3-mm-thick BK7 glass window could be inserted into the sample arm to introduce unbalanced dispersion of β = 132 fs2. A small onion sample was placed inside a 1″ lens tube, submerged in water, covered with a 1-mm-thick BK7 window (Figure 1d), and mounted on a motorized x-y stage. A washer separated the onion and BK7 window by 2.7 mm. A 19-mm achromatic lens on a motorized stage focused the beam inside the sample to a spot size of approximately 6 μm. The delay-arm and sample-arm beams were focused using a 75-mm achromatic lens onto a 0.5-mm BBO crystal, cut for type-I SFG. The beam separation at the lens was 14 mm. The SFG signal was collimated and sent through a monochromator (Princeton Instruments Acton Advanced SP2750A) and detected with a single-photon-counting photomultiplier (Hamamatsu H10682-210). The monochromator was centred at 404.80 nm when the TL-pulse was used, and 404.64 nm when the BARC-pulse was used, to account for a small shift in the signal frequency induced by the added dispersion (see Supplementary Fig. S2 online). The acceptance bandwidth of the monochromator was 0.35 nm in both cases. This filtering of the SFG light was performed in order to remove artifacts from signals acquired with BARC pulses.

Axial depth scans were taken by moving the retro-reflector and recording the photomultiplier signal every 0.4 μm. For the coverglass samples, the delay-stage speed was 0.5 mm/s and the delay stage scanned a range of 700 μm. For the 2D onion data, axial scans were taken over an 800 μm range in the y-direction with one scan every 4 μm. The delay-stage speed was 0.1 mm/s and it was scanned over a total range of 600 μm. The acquisition time per image was 1 hour.

For the histograms in Figure 3, all peaks from each axial scan were fitted provided their amplitudes were between 250 and 6000 counts per delay-stage step, so that their widths were not obscured by noise or detector saturation. Each histogram was fit with a Gaussian peak to estimate the mean and variation of the peak widths in each image.

The 3D onion data was taken over a range of 300 μm, 500 μm, and 350 μm in the x, y, and z directions, respectively. One axial scan was taken every 10 μm in the x and y directions. Every five data points in the z direction were binned to provide a point every 2 μm. Smoothing and threshold algorithms were applied to the raw data to create the 2D images. The 3D structure was visualized with the Imaris (Bitplane, Inc.) software after an FFT bandpass filter was applied. The delay-stage speed was 0.3 mm/s and data acquisition took 2.5 hours.

Author Contributions

K.J.R. and R.K. conceived of the experiment. K.M.S., R.K. and R.P. set up the pulse shaper and cross correlator and performed preliminary experiments. K.M.S. wrote the control software and calibrated the SLM. M.D.M. modified the setup and performed the experiments and simulations presented here, and along with R.P. analyzed the data. M.D.M. and K.J.R. wrote the first draft of the manuscript. All authors contributed to the final version.

Supplementary Material Supplementary Information

Dispersion-cancelled biological imaging with quantum-inspired interferometry-Supplementary Information

We thank Jonathan Lavoie, Kostadinka Bizheva, Donna Strickland, and Rolf Horn for helpful discussions. Sellmeier coefficients for the coverglass slides were provided by Corning Inc. We are grateful for financial support from NSERC, OCE, CFI, QuantumWorks, Industry Canada, and MRI ERA. R.P. acknowledges support from the FWF (J2960-N20), MRI, the VIPS Program of the Austrian Federal Ministry of Science and Research and the City of Vienna as well as the European Commission (Marie Curie, FP7-PEOPLE-2011-IIF). R.K. acknowledges support from the Austrian Academy of Sciences (APART) and support from the European Commission (Marie Curie, FP7-PEOPLE-2010-RG).

Povazay B.,et al.Submicrometer axial resolution optical coherence tomography. Opt. Lett.27, 1800–1802 (Year: 2002).18033368
Fercher A. F.,, Drexler W.,, Hitzenberger C. K., & Lasser T.,Optical coherence tomography — principles and applications. Rep. Prog. Phys.66, 239–303 (Year: 2003).
Zysk A. M.,, Nguyen F. T.,, Oldenburg A. L.,, Marks D. L., & Boppart S. A.,Optical coherence tomography: a review of clinical development from bench to bedside. J. Biomed. Opt.12, 051403 (Year: 2007).17994864
Webster J. L. P.,, Muller M. S., & Fraser J. M.,High speed in situ depth profiling of ultrafast micromachining. Opt. Express15, 14967–14972 (Year: 2007).19550777
Hong C. K.,, Ou Z. Y., & Mandel L.,Measurement of Subpicosecond Time Intervals between Two Photons by Interference. Phys. Rev. Lett.59, 2044–2046 (Year: 1987).10035403
Steinberg A. M.,, Kwiat P. G., & Chiao R. Y.,Dispersion cancellation and high-resolution time measurements in a fourth-order optical interferometer. Phys. Rev. A45, 6659–6665 (Year: 1992).9907789
Steinberg A. M.,, Kwiat P. G., & Chiao R. Y.,Dispersion Cancellation in a Measurement of the Single-Photon Propagation Velocity in Glass. Phys. Rev. Lett.68, 2421–2424 (Year: 1992).10045393
Resch K. J.,, Kaltenbaek R.,, Lavoie J., & Biggerstaff D. N.,Chirped-pulse interferometry with finite frequency correlations. Proc. of SPIE7465, 74650N (Year: 2009).
Abouraddy A. F.,, Nasr M. B.,, Saleh B. E. A.,, Sergienko A. V., &Teich M. C.,Quantum-optical coherence tomography with dispersion cancellation. Phys. Rev. A65, 053817 (Year: 2002).
Altepeter J.,, Jeffrey E., & Kwiat P.,Phase-compensated ultra-bright source of entangled photons. Opt. Express13, 8951–8959 (Year: 2005).19498929
Fedrizzi A.,et al.High-fidelity transmission of entanglement over a high-loss free-space channel. Nature Phys.5, 389–392 (Year: 2009).
Strekalov D. V.,, Pittman T. B., & Shih Y. H.,What we can learn about single photons in a two-photon interference experiment. Phys. Rev. A57, 567–570 (Year: 1998).
Nasr M. B.,, Saleh B. E. A.,, Sergienko A. V., & Teich M. C.,Demonstration of Dispersion-Canceled Quantum-Optical Coherence Tomography. Phys. Rev. Lett.91, 083601 (Year: 2003)14525237
Nasr M. B.,et al.Quantum optical coherence tomography of a biological sample. Opt. Commun.282, 1154–1159 (Year: 2009).
Erkmen B. I., & Shapiro J. H.,Phase-conjugate optical coherence tomography. Phys. Rev. A74, 041601(R) (Year: 2006).
Banaszek K.,, Radunsky A. S., & Walmsley I. A.,Blind dispersion compensation for optical coherence tomography. Opt. Commun.269 152–155 (Year: 2007).
Resch K. J.,, Puvanathasan P.,, Lundeen J. S.,, Mitchell M. W., & Bizheva K.,Classical dispersion-cancellation interferometry. Opt. Express15, 8797–8804 (Year: 2007).19547215
Kaltenbaek R.,, Lavoie J.,, Biggerstaff D. N., & Resch K. J.,Quantum-inspired interferometry with chirped laser pulses. Nature Phys.4 864–868 (Year: 2008).
Le Gouët J.,, Venkatraman D.,, Wong F. N. C., & Shapiro J. H.,Experimental realization of phase-conjugate optical coherence tomography. Opt. Lett.35, 1001–1003 (Year: 2010).20364197
Weiner A. M.,Femtosecond pulse shaping using spatial light modulators. Rev. Sci. Inst.71, 1929–1960 (Year: 2000).
Légaré F. L.,, Fraser J. M.,, Villeneuve D. M., & Corkum P. B.,Adaptive compression of intense 250-nm-bandwidth laser pulses. Appl. Phys. B74(Suppl.), S279–S282 (Year: 2002).
Martínez O. E.,Matrix Formalism for Pulse Compressors. IEEE J. Quantum Electron.234, 2530–2536 (Year: 1988).
Pessot M.,, Maine P., & Mourou G.,1000 times expansion/compression of optical pulses for chirped pulse amplification. Opt. Commun.62, 419–421 (Year: 1987).
Pe'er A.,, Bromberg Y.,, Dayan B.,, Silberberg Y., & Friesem A.,Broadband sum-frequency generation as an efficient two-photon detector for optical tomography. Opt. Express15, 8760–8769 (Year: 2007).19547211
Lavoie J.,, Kaltenbaek R., & Resch K. J.,Quantum-optical coherence tomography with classical light. Opt. Express17, 3818–3825 (Year: 2009).19259223
Kaltenbaek R.,, Lavoie J., & Resch K. J.,Classical Analogues of Two-Photon Quantum Interference. Phys. Rev. Lett.102, 243601 (Year: 2009).19659004
Nasr M. B.,et al.Ultrabroadband Biphotons Generated via Chirped Quasi-Phase-Matched Optical Parametric Down-Conversion. Phys. Rev. Lett.100, 183601 (Year: 2008).18518370
Schreiter K.,Optical pulse shaping for chirped pulse interferometry and bio-imaging. M. Sc. Thesis. University of Waterloo, Canada. (Year: 2011).

Article Categories:
  • Article

Previous Document:  Clinical outcomes and considerations of the lumbar interbody fusion technique for lumbar disk diseas...
Next Document:  The Mental Activity and eXercise (MAX) Trial: A Randomized Controlled Trial to Enhance Cognitive Fun...