Abstract
Non-degenerate two-photon excitation (ND-TPE) has been explored in two-photon excitation microscopy. However, a systematic study of the efficiency of ND-TPE to guide the selection of fluorophore excitation wavelengths is missing. We measured the relative non-degenerate two-photon absorption cross-section (ND-TPACS) of several commonly used fluorophores (two fluorescent proteins and three small-molecule dyes) and generated 2-dimensional ND-TPACS spectra. We observed that the shape of a ND-TPACS spectrum follows that of the corresponding degenerate two-photon absorption cross-section (D-TPACS) spectrum, but is higher in magnitude. We found that the observed enhancements are higher than theoretical predictions.
© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
Two-photon microscopy [1] has had an enormous influence on biological studies providing a tool for high-resolution imaging hundreds of microns within live tissues. In neuroscience, two-photon imaging has been employed to trace neuronal connections [2], monitor biochemical signaling pathways [3], map protein distributions within individual cells and across cell populations [4], report gene expression [5,6], image oxygenation [7,8] and blood flow [9], trace neuronal circuits and image neuronal activity [10], and more. These studies relied on protein-based or synthetic fluorescent reporter molecules (hereafter referred to as fluorophores) suitable for two-photon excitation (TPE) [11].
In conventional two-photon microscopy, a fluorophore is excited by the simultaneous absorption of two photons of the same energy within the near-infrared (NIR) spectrum, derived from the same pulsed laser beam. This is called degenerate two-photon excitation, D-TPE. Alternatively, the same energy needed for the transition to the excited state can be delivered via absorption of two photons of different energy (i.e., different color) [12]. This non-degenerate two-photon excitation (ND-TPE) regime has been explored in two-photon excitation microscopy aiming to extend the excitation wavelength range [12–14], reduce out-of-focus fluorescence [15–18], increase spatial resolution [19,20], and increase penetration depth in scattering media [21,22].
The probability of two-photon excitation, quantified as the excitation cross-section, is an important parameter that, for a given fluorophore concentration and imaging conditions, translates into the laser power required for illumination [23]. However, a systematic study of non-degenerate two-photon absorption cross section (ND-TPACS) spectra of biologically relevant fluorophores, as a function of different combinations of excitation wavelengths, is missing. In addition, for deep imaging in biological tissue (e.g., live brain imaging), the choice of excitation wavelengths needs to be informed by considering the spectrum of water absorption and associated tissue heating.
Here, we performed ND-TPE spectroscopy of two fluorescent proteins (enhanced green fluorescent protein (EGFP) [24] and monomeric Kusabira orange 2 (mKO2) [25]) and three organic dyes (Fluorescein [26], Coumarin 343 [27] and Sulforhodamine 101 (SR101) [28]) using a sensitive fluorescence excitation technique described in our recent study [29]. Further, we explored water absorption and tissue scattering as additional considerations for the selection of the excitation wavelengths for efficient ND-TPE, and compared tissue heating under different excitation regimes using heat conduction simulations.
2. Results
We quantified the ratio of the non-degenerate two-photon absorption cross section (ND-TPACS) to the best achievable degenerate two-photon absorption cross section (D-TPACS) within the equivalent range of the total photon energy (referred to thereafter as peak D-TPAC). In this way, we mitigated the need to measure the collection efficiency of our system at the expense of quantifying the absolute value of the ND-TPACS (see Appendices A, B, and C for experimental setup, spectroscopy protocol, and sample preparation). Plotting ND-TPACS normalized by peak D-TPACS as a function of near-infrared (NIR) and infrared (IR) wavelengths, $\lambda _{\mathrm {NIR}}$ and $\lambda _{\mathrm {IR}}$ respectively, generated a 2-dimensional ND-TPACS spectrum for each examined fluorophore within the available tuning range of the system (Fig. 1(a)). The highest ND-TPACS was obtained when the sum of NIR and IR photon energies matched the ground-to-excited state transition energy. In Fig. 1(a), the excitation energy isocline, a line connecting points with equal total energy $\left(E_{\mathrm{NIR}}+E_{\mathrm{IR}}\right)$ on the 2-dimensional spectrum plot corresponding to the transition energy, is overlaid on the ND-TPACS spectrum for each fluorophore.
In Fig. 1(b), normalized D-TPACS of the five measured fluorophores along with ND-TPACS normalized to the peak D-TPACS are plotted as a function of the equivalent degenerate wavelength $\lambda _{\mathrm {D}}$, where $2hc/\lambda _{\mathrm {D}}=hc/\lambda _{\mathrm {NIR}}+hc/\lambda _{\mathrm {IR}}$ ($h$ is the Planck constant, $c$ is the speed of light in vacuum, and we have assumed equal refractive indices). In this way, we generated a linearized ND-TPACS spectrum that can readily be compared to the corresponding D-TPACS spectrum. Independently measured D-TPACS spectra were similar to the previously reported spectra of the fluorophores [23,30–32]. Several combinations of IR and NIR wavelengths correspond to the same equivalent degenerate wavelength (energy isoclines). Therefore, we report several ND-TPACS values (red dots) for each equivalent wavelength in Fig. 1(b). For each measured fluorophore, the shape of the linearized ND-TPACS spectrum was highly correlated with the D-TPACS spectrum for the longer wavelength region (Fig. 1(b)). However, we observed different amounts of enhancement, as indicated by different range of ND-TPACS values across different fluorophores. The measured fluorescent proteins have a larger spread in ND-TPACS values compared to the measured organic dyes. This might be because fluorescent proteins have a more complex molecular structure. Further studies are required to identify the mechanism behind this phenomenon.
In all examined cases, we observed an enhancement in ND-TPACS (red dots) compared to D-TPACS (black curves); a few red dots below the black curves fall well within our experimental measurement error ($15\%$ relative error, as shown in Appendix B). The observed enhancement in ND-TPACS, compared to D-TPACS, is partially due to the resonance enhancement effect [33–36]: when the energy difference between the high-energy incident photon ($\hbar \omega _{\mathrm {NIR}}$ in our experiment, where $\omega _{\mathrm {NIR}}=2\pi c/\lambda _{\mathrm {NIR}}$ is the frequency of the NIR beam) and that of the nearest one-photon allowed state is decreased, the probability of ND-TPE increases due to a longer lifetime of the reached virtual state as compared to the virtual state corresponding to D-TPE. In other words, the probability of two-photon absorption is proportional to the lifetime of the virtual intermediate excited state, which is inversely proportional to the energy difference between the closest one-photon allowed state and the virtual state [33–36]. This energy difference is called detuning energy $\Delta$, which is maximal in the case of D-TPACS, where the virtual state is in the middle between ground and first excited state. Therefore, for excitation to a specific excited state, D-TPA has a lower probability of occurrence compared to ND-TPA. The detuning energy, in case of excitation into the lowest excited state, is equal to the energy of the low-energy incident photon $\hbar \omega _{\mathrm {IR}}$, where $\omega _{\mathrm {IR}}$ is the frequency of the IR beam. Since all five measured fluorophores are non-centrosymmetric, two-photon absorption is allowed to the lowest excited state. Therefore, we are particularly interested in the case where one of the incident photons approaches the resonance with the lowest one-photon allowed state, known as intermediate state resonance enhancement (ISRE). Using the 2-level essential-state approximation of the sum-over-states (SOS) model, we find that the relation between $\Delta =\hbar \omega _{\mathrm {IR}}$ and ISRE$_2$ for non-centrosymmetric fluorophores is given by the following formula (see Appendix D):
The knowledge of ND-TPACS alone may be enough to guide the choice of wavelengths for imaging relatively thin and transparent samples. For deep imaging in live biological tissues, however, this choice is affected by other considerations, such as light scattering, water absorption, and the associated heating of the sample. In Fig. 2(a), we show the effects of water absorption, and tissue scattering for different wavelengths (adapted from [21] with permission). Longer wavelengths undergo lower scattering in tissue; however, as shown in Fig. 2(a), water absorption increases significantly at wavelengths near $1500$ nm and beyond $1800$ nm. Higher water absorption at longer wavelengths will result in greater tissue heating. In ND-TPE, the IR beam has longer wavelengths compared to D-TPE, and excessive tissue heating is a potential concern for ND-TPE, even when equal total excitation power is used. However, it has been shown that the increase in tissue temperature is much smaller than the increase of the water absorption coefficient [37] and therefore, the increase in tissue heating for longer wavelengths is not as drastic as expected. This is also evident in tissue heating simulations shown in Fig. 2(b), where we plot the maximum temperature change of the sample when illuminated by both NIR and IR beams, $\Delta T_{\mathrm {NIR+IR}}=\Delta T_{\mathrm {NIR}}+\Delta T_{\mathrm {IR}}$, relative to the maximum temperature change upon illumination with a single beam at the equivalent degenerate wavelength, $\Delta T_{\mathrm {D}}$, assuming equal total powers (see Appendix E for the description of the simulation). For example, the ratio of the water absorption coefficient at $1340$ nm to the absorption coefficient at $920$ nm is $\sim 30$ (red plot in Fig. 2(a) and Table 3), but the ratio of the temperature change between D-TPE with $\lambda _{\mathrm {D}}=920$ nm and ND-TPE with $\lambda _{\mathrm {NIR}}=700$ nm and $\lambda _{\mathrm {IR}}=1340$ nm and equal excitation power is $\sim 3$. This result agrees with the idea that most of the incident light is absorbed before escaping the thick tissue sample and therefore, the variation in the absorption coefficient of different wavelengths does not drastically affect the rate of tissue heating [37].
3. Conclusion
In conclusion, we showed that, within our measured excitation range, the ND-TPACS spectrum of all fluorophores tested follow their corresponding D-TPACS in shape and exhibit enhancement in magnitude. We showed that the amount of the observed enhancement is higher than the theoretical prediction given by the 2-level essential-state approximation of the SOS model. Different fluorophores exhibit different amounts of resonance enhancement and we observed up to $\sim 75\%$ enhancement in the ND-TPACS compared to the corresponding D-TPACS. We also showed that the increase in tissue heating, due to excitation with longer wavelengths, is not as large as water absorption increase for longer wavelengths. This increase in heating may be an acceptable trade-off for leveraging enhanced ND-TPACS and other advantages of ND-2PE [22] in deep tissue imaging applications.
Appendix A: Experimental setup
Detailed description of the spectroscopy setup is found in our recent work [29]. In brief, two temporally synchronized and spatially aligned pulsed femtosecond laser beams, one with near infrared (NIR) and another with infrared (IR) energies, were derived from a Ti:Sapphire laser and optical parametric oscillator (OPO), respectively. The beams, linearly polarized in the same plane, were combined on a dichroic mirror and focused by a microscope objective (Olympus, LMPLN-10X-IR, NA=$0.3$) into the sample solution (Fig. 3(a)). The NIR and IR wavelengths were scanned in $10$ nm increments in the range of $740-870$ nm and $1100-1400$ nm, respectively. We measured and recorded the fluorescence intensity and the following parameters for each laser beam: wavelength, spectral pulse bandwidth, average power, and repetition rate. The NIR optical path included a delay stage, so that the time-averaged intensity of the fluorescence could be measured as a function of the temporal delay between the NIR and IR pulses (Fig. 3(b)). All measurements were conducted at sufficiently low incident powers, ensuring that two-photon excitation was the only mechanism that generated the fluorescence signal and other non-linear phenomena such as stimulated emission and excited-state absorption are absent [29]. To show that the ND-TPE was the only fluorescence excitation mechanism, we plotted the fluorescence signal value at the peak of the temporal overlap (fluorescence generated by both ND-TPE and D-TPE) minus the fluorescence signal at no temporal overlap (fluorescence generated by D-TPE) as a function of the excitation power of each individual beam as shown in Fig. 4. The linear dependence of the fluorescence signal generated by ND-TPE on the excitation power of each beam is the characteristic signature of the ND-TPE [22,29].
Appendix B: Spectroscopy protocol
The fluorescence intensity measured as a function of the delay-line position (Fig. 3(b)), resulted in a cross-correlation curve that was generated for each combination of NIR and IR wavelengths (Fig. 3(b)). When NIR and IR pulses arrived at different times (Temporal offset $> 0$), the fluorescence signal was generated by degenerate two-photon excitation (D-TPE) due to absorption of two photons of equal energy coming from either the NIR or IR beam. An increase in the fluorescence signal was detected as laser pulses from the two beams overlapped in time indicating the additive effect of the non-degenerate two-photon excitation (ND-TPE) occurring due to simultaneous absorption of one photon from each NIR and IR beam. To isolate the ND-TPE signal, for each combination of NIR and IR wavelengths, we subtracted the D-TPE contribution from the fluorescence intensity at zero delay between NIR and IR pulses. Next, we compared the probability of two-photon absorption between D-TPE and ND-TPE. In D-TPE, the time-averaged fluorescence signal can be expressed as the following equation:
Measurement of the absolute ND-TPACS would require knowledge of the fluorophore quantum efficiency $\eta$ and the collection efficiency of the measurement system $\phi$. Our goal, however, was to compare ND-TPE to D-TPE and evaluate the possibility of enhancement of two-photon absorption upon ND-TPE. Therefore, for each fluorophore, we quantified the ratio of ND-TPACS to the peak (i.e., the best achievable) D-TPACS and not the absolute value of ND-TPACS. To that end, we normalized all ND-TPACS measurements by the peak D-TPACS, which occurred at the wavelength where the photon energy matched half of the total energy needed for the transition to the excited state. This is because under D-TPE, two photons are drawn from the same beam and have the same photon energy. D-TPACS was independently measured within the same range of the total energy of two photons by blocking the IR beam. This normalization mitigated the need to explicitly measure $\phi$ and $\eta$, because these terms disappear when dividing Eq. (3) by Eq. (2). For each fluorophore, the normalized ND-TPACS spectrum was computed by fitting cross-correlation curves (Fig. 3(b)) with Eq. (3) for each pair of NIR and IR wavelengths and normalizing the obtained valued by the peak D-TPACS that is a fluorophore-specific constant. For a more detailed description of the spectroscopy protocol please refer to our recent work on the spectroscopy technique [29]. We characterized the measurement error via repetitive measurement of the ND-TPACS along one isocline (Fig. 5). The main source of error was the uncertainty in the measurement of IR beam pulse width, fluorescence signal value, and laser powers.
Appendix C: Sample preparation
Fluorophores were dissolved in phosphate-buffered saline (PBS, pH $\sim 7.4$) or dimethyl sulfoxide (DMSO) to final concentrations listed in Table 2. The concentrations were measured using one-photon absorption spectroscopy (NanoDrop 2000 spectrometer, ThermoFisher) and published molar extinction coefficients. One photon-absorption spectra of inspected fluorophores are shown in Fig. 6.
Appendix D: Derivation of the intermediate state resonance enhancement (ISRE)
In this section we derive the formula for ISRE$_2$ (Eq. (1) in the main text). We start with the expression for TPA tensor $\mathbf {S}^{(f)}$ for the transition from ground state $|g\rangle$ to final state $|f\rangle$ which is given by:
Let’s assume, for simplicity that all transitions are polarized in the same direction and this direction is parallel to the two polarization vectors $\boldsymbol {e_1}$ and $\boldsymbol {e_2}$. This assumption will not affect the generality of the conclusions. Under this assumption, the $|\mathbf {S}^{(f)}|^2$ expression can be written in a simpler form:
Using the same procedure as for the 2-level system given above, here we derived the ISRE for a 3-level system with $|i\rangle$ as the additional intermediate state:
Appendix E: Heating comparison
To study potential heating under ND-TPE, we simulate tissue heating through scanned focused light using a MATLAB code, which was provided by K. Podgorski [37]. The code uses Monte Carlo simulations to simulate light spread in tissue and the finite differencing method for heat conduction, metabolic heating, and cooling by blood perfusion [37]. The code assumes a heterogeneous sample incorporating a $350$-$\mu$m glass window and immersion water above the sample [37]. The simulated focused beam is set to scan a $1$ mm$^2$ area at $250$ $\mu$m below the surface. We used previously published values for optical parameters of brain tissue [21,47] (Table S2). For simulation of heating under simultaneous excitation with NIR and IR beams, we used the code to simulate the heating under individual NIR and IR beams and then added the resulting tissue temperature changes to obtain the heating for the case that both beams are illuminating the sample at the same time. This is a viable approach since the superposition principle is valid for heat propagation which is described by a linear partial differential equation. Podgorski et al. showed that heating varies in different parts of the sample relative to the focal point and the temperature change is greatest deep below the brain surface and focal plane [37].
In Fig. 8(a) we show the spatial temperature profile of brain tissue in resting equilibrium and when illuminated using a NIR wavelength of $740$ nm, an IR wavelength of $1230$ nm and the equivalent degenerate wavelength of $920$ nm ($2/\lambda _{\mathrm {D}}=1/\lambda _{\mathrm {NIR}}+1/\lambda _{\mathrm {IR}}$). Light is focused and scanned $250$ $\mu$m below the cover glass over a $1$ mm$^2$ scan area. Power of the excitation beams were chosen such that $P_{\mathrm {NIR}}+P_{\mathrm {IR}}=P_{\mathrm {D}}=100$ mW and $P_{\mathrm {NIR}}=P_{\mathrm {IR}}$ at focal depth ($250$ $\mu$m below cover glass). For example, at $\lambda _{\mathrm {NIR}}=740$ nm we set $P_{\mathrm {NIR}}$ ($z=250$ $\mu$m) $=50$ mW at focal depth which means the input laser power at surface is $P_{\mathrm {NIR}}$ ($z=0$ $\mu$m) $=100\times (P_{\mathrm {NIR}}$ ($z=250$ $\mu$m)/$T(\lambda _{ \mathrm {NIR}}=740$ nm; $z=250$ $\mu$m)) $=250$ mW. Tissue exhibits the highest temperature change when illuminated with the IR beam (Fig. 8(b)) as expected given high absorption of water at this wavelength.
Funding
National Institutes of Health (MH111359, NS057198, S10RR029050, U01NS094232); National Science Foundation (CBET 1445158, ECCS 1405234, ECCS 1507146); Office of Naval Research (ECE4509); Cymer.
Acknowledgments
We thank Drs. Sergei A. Vinogradov and Luca Ravotto for their help in interpreting the experimental results and discussions of the theory of two-photon absorption in molecules. We also thank Dr. Kaspar Podgorski for kindly providing us with the tissue heating simulation code and helping us with the code implementation. We thank Drs. Andrew Dunn and David Miller for sharing numeric data from their publication.
References
1. W. Denk, J. H. Strickler, and W. W. Webb, “Two-photon laser scanning fluorescence microscopy,” Science 248(4951), 73–76 (1990). [CrossRef]
2. W.-C. A. Lee, J. L. Chen, H. Huang, J. H. Leslie, Y. Amitai, P. T. So, and E. Nedivi, “A dynamic zone defines interneuron remodeling in the adult neocortex,” Proc. Natl. Acad. Sci. 105(50), 19968–19973 (2008). [CrossRef]
3. Y. Okubo, S. Kakizawa, K. Hirose, and M. Iino, “Visualization of ip3 dynamics reveals a novel AMPA receptor-triggered ip3 production pathway mediated by voltage-dependent ca2+ influx in purkinje cells,” Neuron 32(1), 113–122 (2001). [CrossRef]
4. A. A. Oliva, M. Jiang, T. Lam, K. L. Smith, and J. W. Swann, “Novel hippocampal interneuronal subtypes identified using transgenic mice that express green fluorescent protein in GABAergic interneurons,” J. Neurosci. 20(9), 3354–3368 (2000). [CrossRef]
5. M. Mank, A. F. Santos, S. Direnberger, T. D. Mrsic-Flogel, S. B. Hofer, V. Stein, T. Hendel, D. F. Reiff, C. Levelt, and A. Borstet al., “A genetically encoded calcium indicator for chronic in vivo two-photon imaging,” Nat. Methods 5(9), 805–811 (2008). [CrossRef]
6. N. L. Rochefort and A. Konnerth, “Genetically encoded ca 2+ sensors come of age,” Nat. Methods 5(9), 761–762 (2008). [CrossRef]
7. S. Sakadžić, E. Roussakis, M. A. Yaseen, E. T. Mandeville, V. J. Srinivasan, K. Arai, S. Ruvinskaya, A. Devor, E. H. Lo, S. A. Vinogradov, and D. A. Boas, “Two-photon high-resolution measurement of partial pressure of oxygen in cerebral vasculature and tissue,” Nat. Methods 7(9), 755–759 (2010). [CrossRef]
8. A. Devor, S. Sakadžić, P. A. Saisan, M. A. Yaseen, E. Roussakis, V. J. Srinivasan, S. A. Vinogradov, B. R. Rosen, R. B. Buxton, A. M. Dale, and D. A. Boas, “Overshoot of o2 is required to maintain baseline tissue oxygenation at locations distal to blood vessels,” J. Neurosci. 31(38), 13676–13681 (2011). [CrossRef]
9. D. Kleinfeld, P. P. Mitra, F. Helmchen, and W. Denk, “Fluctuations and stimulus-induced changes in blood flow observed in individual capillaries in layers 2 through 4 of rat neocortex,” Proc. Natl. Acad. Sci. 95(26), 15741–15746 (1998). [CrossRef]
10. S. Nagayama, S. Zeng, W. Xiong, M. L. Fletcher, A. V. Masurkar, D. J. Davis, V. A. Pieribone, and W. R. Chen, “In vivo simultaneous tracing and ca2+ imaging of local neuronal circuits,” Neuron 53(6), 789–803 (2007). [CrossRef]
11. E. M. Callaway, “A molecular and genetic arsenal for systems neuroscience,” Trends Neurosci. 28(4), 196–201 (2005). [CrossRef]
12. P. Mahou, M. Zimmerley, K. Loulier, K. S. Matho, G. Labroille, X. Morin, W. Supatto, J. Livet, D. Débarre, and E. Beaurepaire, “Multicolor two-photon tissue imaging by wavelength mixing,” Nat. Methods 9(8), 815–818 (2012). [CrossRef]
13. S. Quentmeier, S. Denicke, and K.-H. Gericke, “Two-color two-photon fluorescence laser scanning microscopy,” J. Fluoresc. 19(6), 1037–1043 (2009). [CrossRef]
14. E. P. Perillo, J. W. Jarrett, Y.-L. Liu, A. Hassan, D. C. Fernée, J. R. Goldak, A. Bonteanu, D. J. Spence, H.-C. Yeh, and A. K. Dunn, “Two-color multiphoton in vivo imaging with a femtosecond diamond raman laser,” Light: Sci. Appl. 6(11), e17095 (2017). [CrossRef]
15. M. O. Cambaliza and C. Saloma, “Advantages of two-color excitation fluorescence microscopy with two confocal excitation beams,” Opt. Commun. 184(1-4), 25–35 (2000). [CrossRef]
16. C. M. Blanca and C. Saloma, “Two-color excitation fluorescence microscopy through highly scattering media,” Appl. Opt. 40(16), 2722–2729 (2001). [CrossRef]
17. C. Wang, L. Qiao, Z. Mao, Y. Cheng, and Z. Xu, “Reduced deep-tissue image degradation in three-dimensional multiphoton microscopy with concentric two-color two-photon fluorescence excitation,” J. Opt. Soc. Am. B 25(6), 976–982 (2008). [CrossRef]
18. D. Kobat, G. Zhu, and C. Xu, “Background reduction with two-color two-beam multiphoton excitation,” in Biomedical Optics, (Optical Society of America, 2008), p. BMF6.
19. S. Lindek and E. Stelzer, “Resolution improvement by nonconfocal theta microscopy,” Opt. Lett. 24(21), 1505–1507 (1999). [CrossRef]
20. C. Ibáïez-López, I. Escobar, G. Saavedra, and M. Martínez-Corral, “Optical-sectioning improvement in two-color excitation scanning microscopy,” Microsc. Res. Tech. 64(2), 96–102 (2004). [CrossRef]
21. D. R. Miller, J. W. Jarrett, A. M. Hassan, and A. K. Dunn, “Deep tissue imaging with multiphoton fluorescence microscopy,” Curr. Opin. Biomed. Eng. 4, 32–39 (2017). [CrossRef]
22. M.-H. Yang, M. Abashin, P. A. Saisan, P. Tian, C. G. Ferri, A. Devor, and Y. Fainman, “Non-degenerate 2-photon excitation in scattering medium for fluorescence microscopy,” Opt. Express 24(26), 30173–30187 (2016). [CrossRef]
23. C. Xu and W. W. Webb, “Measurement of two-photon excitation cross sections of molecular fluorophores with data from 690 to 1050 nm,” J. Opt. Soc. Am. B 13(3), 481–491 (1996). [CrossRef]
24. B. P. Cormack, R. H. Valdivia, and S. Falkow, “FACS-optimized mutants of the green fluorescent protein (gfp),” Gene 173(1), 33–38 (1996). [CrossRef]
25. A. Sakaue-Sawano, H. Kurokawa, T. Morimura, A. Hanyu, H. Hama, H. Osawa, S. Kashiwagi, K. Fukami, T. Miyata, H. Miyoshi, T. Imamura, M. Ogawa, H. Masai, and A. Miyawaki, “Visualizing spatiotemporal dynamics of multicellular cell-cycle progression,” Cell 132(3), 487–498 (2008). [CrossRef]
26. R. Sjöback, J. Nygren, and M. Kubista, “Absorption and fluorescence properties of fluorescein,” Spectrochim. Acta, Part A 51(6), L7–L21 (1995). [CrossRef]
27. G. Reynolds and K. H. Drexhage, “New coumarin dyes with rigidized structure for flashlamp-pumped dye lasers,” Opt. Commun. 13(3), 222–225 (1975). [CrossRef]
28. A. Nimmerjahn, F. Kirchhoff, J. N. Kerr, and F. Helmchen, “Sulforhodamine 101 as a specific marker of astroglia in the neocortex in vivo,” Nat. Methods 1(1), 31–37 (2004). [CrossRef]
29. S. Sadegh, M.-H. Yang, C. G. Ferri, M. Thunemann, P. A. Saisan, A. Devor, and Y. Fainman, “Measurement of the relative non-degenerate two-photon absorption cross-section for fluorescence microscopy,” Opt. Express 27(6), 8335–8347 (2019). [CrossRef]
30. M. Drobizhev, N. S. Makarov, S. E. Tillo, T. E. Hughes, and A. Rebane, “Two-photon absorption properties of fluorescent proteins,” Nat. Methods 8(5), 393–399 (2011). [CrossRef]
31. J. Mütze, V. Iyer, J. J. Macklin, J. Colonell, B. Karsh, Z. Petrášek, P. Schwille, L. L. Looger, L. D. Lavis, and T. D. Harris, “Excitation spectra and brightness optimization of two-photon excited probes,” Biophys. J. 102(4), 934–944 (2012). [CrossRef]
32. S. de Reguardati, J. Pahapill, A. Mikhailov, Y. Stepanenko, and A. Rebane, “High-accuracy reference standards for two-photon absorption in the 680–1050 nm wavelength range,” Opt. Express 24(8), 9053–9066 (2016). [CrossRef]
33. R. R. Birge and B. M. Pierce, “Semiclassical time-dependent theory of two-photon spectroscopy. the effect of dephasing in the virtual level on the two-photon excitation spectrum of isotachysterol,” Int. J. Quantum Chem. 29(4), 639–656 (1986). [CrossRef]
34. J. M. Hales, D. J. Hagan, E. W. Van Stryland, K. Schafer, A. Morales, K. D. Belfield, P. Pacher, O. Kwon, E. Zojer, and J.-L. Brédas, “Resonant enhancement of two-photon absorption in substituted fluorene molecules,” J. Chem. Phys. 121(7), 3152–3160 (2004). [CrossRef]
35. M. Drobizhev, N. S. Makarov, A. Rebane, G. de la Torre, and T. Torres, “Strong two-photon absorption in push- pull phthalocyanines: Role of resonance enhancement and permanent dipole moment change upon excitation,” J. Phys. Chem. C 112(3), 848–859 (2008). [CrossRef]
36. H. S. Pattanaik, M. Reichert, J. B. Khurgin, D. J. Hagan, and E. W. Van Stryland, “Enhancement of two-photon absorption in quantum wells for extremely nondegenerate photon pairs,” IEEE J. Quantum Electron. 52(3), 1–14 (2016). [CrossRef]
37. K. Podgorski and G. Ranganathan, “Brain heating induced by near-infrared lasers during multiphoton microscopy,” J. Neurophysiol. 116(3), 1012–1023 (2016). [CrossRef]
38. N. C. Shaner, P. A. Steinbach, and R. Y. Tsien, “A guide to choosing fluorescent proteins,” Nat. Methods 2(12), 905–909 (2005). [CrossRef]
39. K. M. Dean and A. E. Palmer, “Advances in fluorescence labeling strategies for dynamic cellular imaging,” Nat. Chem. Biol. 10(7), 512–523 (2014). [CrossRef]
40. G. McNamara, M. Difilippantonio, T. Ried, and F. R. Bieber, “Microscopy and image analysis,” Curr. Protoc. Hum. Genet. 46(1), 441–4434 (2005). [CrossRef]
41. W. L. Peticolas, “Multiphoton spectroscopy,” Annu. Rev. Phys. Chem. 18(1), 233–260 (1967). [CrossRef]
42. J. Bhawalkar, G. He, and P. Prasad, “Nonlinear multiphoton processes in organic and polymeric materials,” Rep. Prog. Phys. 59(9), 1041–1070 (1996). [CrossRef]
43. T. V. Esipova, H. J. Rivera-Jacquez, B. Weber, A. E. Masunov, and S. A. Vinogradov, “Two-photon absorbing phosphorescent metalloporphyrins: effects of $\pi$-extension and peripheral substitution,” J. Am. Chem. Soc. 138(48), 15648–15662 (2016). [CrossRef]
44. C. W. Dirk, L.-T. Cheng, and M. G. Kuzyk, “A simplified three-level model describing the molecular third-order nonlinear optical susceptibility,” Int. J. Quantum Chem. 43(1), 27–36 (1992). [CrossRef]
45. A. Rebane, M. Drobizhev, N. Makarov, E. Beuerman, S. Tillo, and T. Hughes, “New all-optical method for measuring molecular permanent dipole moment difference using two-photon absorption spectroscopy,” J. Lumin. 130(9), 1619–1623 (2010). [CrossRef]
46. A. Rebane, G. Wicks, M. Drobizhev, T. Cooper, A. Trummal, and M. Uudsemaa, “Two-photon voltmeter for measuring a molecular electric field,” Angew. Chem., Int. Ed. 54(26), 7582–7586 (2015). [CrossRef]
47. N. G. Horton, K. Wang, D. Kobat, C. G. Clark, F. W. Wise, C. B. Schaffer, and C. Xu, “In vivo three-photon microscopy of subcortical structures within an intact mouse brain,” Nat. Photonics 7(3), 205–209 (2013). [CrossRef]