Influence of birefringence on polarization resolved nonlinear microscopy and collagen SHG structural imaging

We analyze the influence of the anisotropy of molecular and biological samples on polarization resolved nonlinear microscopy imaging. We show in particular the detrimental influence of birefringence on Second Harmonic Generation (SHG) and Two-Photon Excited Fluorescence (TPEF) polarization resolved microscopy imaging, which, if not accounted for, can lead to an erroneous determination of the sample properties and thus to a misinterpretation of the read-out information. We propose a method to measure this birefringence and account for this effect in nonlinear polarization resolved experiments. © 2010 Optical Society of America OCIS codes: (180.4315) Nonlinear microscopy; (190.1900) Diagnostic applications of nonlinear optics. References and links 1. J. Gannaway and C. J. R. Sheppard, “Second harmonic imaging in the scanning optical microscope,” Opt. Quantum Electron. 10, 435–439 (1978). 2. I. Freund, M. Deutsch, and A. Sprecher, “Connective tissue polarity. Optical second-harmonic microscopy, crossed-beam summation, and small-angle scattering in rat-tail tendon,” Biophys. J. 50, 693–712 (1986). 3. M. Flörsheimer, C. Radüge, H. Salmen, M. Bösch, R. Terbrack, and H. Fuchs, “In-situ imaging of Langmuir monolayers by second-harmonic microscopy,” Thin Solid Films 284, 659–662 (1996). 4. L. Moreaux, O. Sandrea, and J. Mertz, “Membrane Imaging by Second Harmonic Generation Microscopy,” J. Opt. Soc. Am. B 17, 1685–1689 (2000). 5. P. J. Campagnola, A. C. Millard, M. Terasaki, P. E. Hoppe, C. J. Malone, and W. A. Mohler, “3-Dimesional HighResolution Second Harmonic Generation Imaging of Endogenous Structural Proteins in Biological Tissues,” Biophys J. 81, 493–508 (2002). 6. D. A. Dombeck, K. A. Kasischke, H. D. Vishwasrao, M. Ingelsson, B. T. Hyman, and W. W. Webb, “Uniform polarity microtubule assemblies imaged in native brain tissue by second-harmonic generation microscopy,” Proc. Natl. Acad. Sci. U.S.A. 100, 7081 (2003). 7. M. Strupler, M. Hernest, C. Fligny, J. L. Martin, P.-L. Tharaux, and M. C. Schanne-Klein, “Second Harmonic Microscopy to Quantify Renal Interstitial Fibrosis and Arterial Remodeling,” J. Biomed. Opt. 13, 054041 (2008). 8. O. Nadiarnykh, R. LaComb, M. Brewer, and P. J. Campagnola, “Second Harmonic Generation Imaging Microscopy of Ovarian Cancer,” Biophys. J. 94, 4504–4514 (2008). 9. S. W. Teng, H. Y. Tan, J. L. Peng, H. H. Lin, K. H. Kim, W. Lo,, Y. Sun, W. C. Lin, S. J. Lin, S. H. Jee, P. T. C. So, and C.Y. Dong, “Multiphoton autofluorescence and second-harmonic generation (SHG) imaging of ex-vivo porcine eye,” Invest. Ophthalmol. Vis. Sci. 47, 1216–1224 (2006). #127860 $15.00 USD Received 3 May 2010; accepted 26 May 2010; published 28 Jun 2010 (C) 2010 OSA 5 July 2010 / Vol. 18, No. 14 / OPTICS EXPRESS 14859 10. E. Ralston, B. Swaim, M. Czapiga, W. L. Hwu, Y. H. Chien, M. G. Pittis, B. Bembi, O. Schwartz, P. Plotz, and N. Raben,“Detection and imaging of non-contractile inclusions and sarcomeric anomalies in skeletal muscle by second harmonic generation combined with two-photon excited fluorescence,” J. Struct. Biol. 162, 500–508 (2008). 11. D. Axelrod, “Carbocyanine dye orientation in red cell membrane studied by microscopic fluorescence polarization,” Biophys. J. 26, 557–573 (1979). 12. J. Borejdo and S. Burlacu, “Measuring Orientation of Actin Filaments within a Cell: orientation of Actin in Intestinal Microvilli,” Biophys. J. 65, 300–309 (1993). 13. M. Flörsheimer, M. Bösch, C. Brillert, M. Wierschem, and H. Fuchs, “Second-harmonic microscopy a quantitative probe for molecular surface order,” Adv. Mater. 9, 1061-1065 (1997). 14. C. Anceau, S. Brasselet, and J. Zyss, ”Local orientational distribution of molecular monolayers probed by nonlinear microscopy,” Chem. Phys. Lett. 411, 98–102 (2005). 15. P. Stoller, B. M. Kim, A. M. Rubenchik, K. M. Reiser, and L. B. Da Silva, ”Polarization-dependent optical second harmonic imaging of a rat-tail tendon,” J. Biomed. Opt. 7, 205–214 (2000). 16. P. Stoller, K. M. Reiser, P. M. Celliers, and A. M. Rubenchik, “Polarization-Modulated Second Harmonic Generation in Collagen,” Biophys. J. 82, 3330–3342 (2002). 17. T. Yasui, K. Sasaki, Y. Tohno, and T. Araki, “Tomographic imaging of collagen fiber orientation in human tissue using depth-resolved polarimetry of second-harmonic-generation,” Opt. Quantum Electron. 37, 1397–1408 (2005). 18. R. M. Williams, W. R. Zipfel, and W. W. Webb, ”Interpreting Second-Harmonic Generation Images of Collagen I Fibrils,” Biophys. J. 88, 1377-1386 (2005). 19. S. Yen, W. L. Chen, H. Y. Wu, F. C. Li, W. Lo, S. J. Lin, S. H. Jee, Y. F. Chen, P. T. C. So, and C. Y. Dong, “Collagen Thermal denaturation study with high resolution second harmonic generation imaging and polarization optical microscopy,” Biophys. J. 91, 2620–2625 (2006). 20. F. Tiaho, G. Recher, and D. Rouède, “Estimation of helical angles of myosin and collagen by second harmonic generation imaging microscopy,” Opt. Express 15, 4054–4065 (2007). 21. C. Odin, Y. Le Grand, A. Renault, L. Gailhouste, and G. Baffet, “Orientation fields of nonlinear biological fibrils by second harmonic generation microscopy,” J. Microsc. 229, 32–38 (2008). 22. S. Psilodimitrakopoulos, S. I. C. O. Santos, I. Amat-Roldan, A. K. N. Thayil, D. Artigas, and P. Loza-Alvarez, “In vivo, pixel-resolution mapping of thick filaments’ orientation in nonfibrilar muscle using polarization-sensitive second harmonic generation microscopy,” J. Biomed. Opt. 14, 014001 (2009). 23. A. Gasecka, T.-J. Han, C. Favard, B. R. Cho, and S. Brasselet, “Quantitative imaging of molecular order in lipid membranes using two-photon fuorescence polarimetry,” Biophys. J. 97, 2854–2862 (2009). 24. A. Gasecka, L.-Q. Dieu, D. Bruehviler, and S. Brasselet, “Probing molecular order in zeolite L inclusion compounds using two-photon fluorescence polarimetric microscoy,” J. Phys. Chem. B 114, 4192–4198 (2010). 25. V. Le Floc’h, S. Brasselet, J.-F. Roch, and J. Zyss, “Monitoring of orientation in molecular ensembles by polarization sensitive nonlinear microscopy,” J. Phys. Chem. B 107, 12403–12410 (2003). 26. S. Brasselet, V. Le Floc’h, F. Treussart, J.-F. Roch, J. Zyss, E. Botzung-Appert, and A. Ibanez, “In situ diagnostics of the crystalline nature of single organic nanocrystals by nonlinear microscopy,” Phys. Rev. Lett. 92, 207401– 204405 (2004). 27. K. Komorowska, S. Brasselet, J. Zyss, L. Pourlsen, M. Jazdzyk, H. J. Egelhaaf, J. Gierschner, and M. Hanack, “Nanometric scale investigation of the nonlinear efficiency of perhydrotriphynylene inclusion compounds,” Chem. Phys. 318, 12–20 (2005). 28. S. Brasselet, and J. Zyss, “Nonlinear polarimetry of molecular crystals down to the nanoscale,” C. R. Physique 8, 165–179 (2007). 29. P. Schön, F. Munhoz, A. Gasecka, S. Brustlein, and S. Brasselet, “Polarization distortion effects in polarimetric two-photon microscopy,” Opt. Express 16, 20891–20901 (2008). 30. B. Richards, and E. Wolf, “Electromagnetic Diffraction in Optical Systems. II. Structure of the Image Field in an Aplanatic System ,” Proc. R. Soc. London Ser. A. 253, 358–379 (1959). 31. E. Y. S. Yew, and C. R. J. Sheppard, “Effects of axial field components on second harmonic generation microscopy ,” Opt. Express 14, 1167–1174 (2006). 32. N. Sandeau, L. L. Xuan, D. Chauvat, C. Zhou, J. F. Roch, and S. Brasselet, “ Defocused imaging of second harmonic generation from a single nanocrystal,” Opt. Express 15, 16051–16060 (2007). 33. P. Schön, M. Behrndt, D. Ait-Belkacem, H. Rigneault, and S. Brasselet, “Polarization and Phase Pulse Shaping applied to Structural Contrast in Nonlinear Microscopy Imaging,” Phys. Rev. A 81, 013809 (2010). 34. D. Lara and C. Dainty,“Axially resolved complete Mueller matrix confocal microscopy,” Appl. Opt. 45, 1917– 1930 (2006). 35. C. E. Bigelow, and T. H. Foster, “Confocal fluorescence polarization microscopy in turbid media: Effects of scattering-induced depolarisation,” J. Opt. Soc. Am. A 23, 2932 (2006). 36. I. Ledoux, C. Lepers, A. Prigaud, J. Badan, and J. Zyss, “Linear and nonlinear optical properties of N-4nitrophenyl L-prolinol single crystals,” Opt. Commun. 80, 149–154 (1990). 37. F. P. Bolin, L. E. Preuss, R. C. Taylor, and R. J. Ference, “Refractive index of some mammalian tissues using a #127860 $15.00 USD Received 3 May 2010; accepted 26 May 2010; published 28 Jun 2010 (C) 2010 OSA 5 July 2010 / Vol. 18, No. 14 / OPTICS EXPRESS 14860 fiber optic cladding method,” Appl. Opt. 28, 2297–2303 (1989). 38. D. T. Poh, “Examination of refractive index of human epidermis in-vitro and in-vivo.,” Proc. Inter. Conf. Lasers. 96, 118–125 (1996).


Introduction
Since its first developments [1][2][3] and its introduction in bio-imaging [4][5][6], second harmonic generation (SHG) microscopy is now widely used to image ordered biomolecular assemblies in complex samples at depths reaching a few hundreds of micrometers.Coherent SHG occurring naturally in non-centrosymmetric structures such as collagen [2], skeletal muscles [5] and microtubules [6], is today exploited as a functional contrast [7,8], possibly in conjunction with Two-Photon Excited Fluorescence (TPEF) [9,10], with the ultimate goal of developing diagnostics of pathological effects related to tissues and cell architecture.
In addition to their unique imaging capabilities, these contrasts are dependent on the incident light polarization state, providing an interesting way to probe molecular orientation and disorder.Structural information, from membrane architecture and proteins aggregates to biopolymers and tissue assemblies, is a key parameter in a large variety of biological phenomena.Polarization probe imaging has been exploited a long time ago in fluorescence [11,12] and SHG [13,14] in ordered molecular samples.Its extension to SHG in tissues leads to a large amount of current studies [15][16][17][18][19][20][21][22].
Recent works in polarization resolved nonlinear microscopy (or nonlinear polarimetry) have pointed out the need to rely on a variable state of incident polarization to retrieve relevant information on the sample structure.In fluorescence, measuring the signal from a continuous variation of the incident linear polarization state circumvents the limitations of the traditional fluorescence anisotropy imaging, which explores a limited number of polarization states and therefore requires a priori knowledge information on the studied system [11,12].TPEF polarimetry has shown in particular its ability to reveal non-ambiguous information on molecular angular distributions in cell membranes and crystals [23,24], with a typical 300nm lateral resolution.
Although polarimetry microscopy imaging is a powerful technique, relevant information can be retrieved only if particular care is brought on the possible polarization distortions occurring in the instrument.The degree of ellipticity appearing from reflections on mirrors and dichroic mirrors can lead to a deformation of the retrieved information, in particularly when using the epi-geometry [29].In addition, the excitation and detection polarization are scrambled due to the use of high numerical aperture objectives [11].Finally an extra coupling along the propagation direction appears to be non-negligible when using high aperture focusing [30].These different polarization distortion factors have been thoroughly addressed and can be readily accounted for in microscopy polarization analysis [25,[31][32][33].
Apart from these instrumental effects, other polarization distortions can originate from the sample itself, especially when propagating through micrometric scale depths.Scattering involves more specifically depolarization, which can be identified and eventually corrected for by measuring depolarization properties using for instance the Mueller matrix formalism [34,35].In this work we focus on birefringence, which can be non-negligible and even sometimes significant in crystalline [36] and biological samples, in particular from fibrillar structures such as collagen even at a few micrometer depths [37,38].Birefringence, which introduces a phase shift in the input field polarization state between optical axes of a priori unknown orientations, can lead to erroneous deductions on the measured properties such as SHG nonlinear tensorial components.We show that this effect in collagen, which has been so far evoked in only a few works in nonlinear polarized microscopy [16], and ignored in most of the previous polarization resolved studies [17,[20][21][22], can be detrimental to the measured signals.In this work, we analyze the effect of birefringence in polarimetric TPEF and SHG microscopy applied to anisotropic samples widely studied in molecular and biological sciences.Although these contrasts are chosen for the possibility to combine them on the same microscope, the present study can be easily adapted to one-photon fluorescence or other nonlinear processes.

Effect of birefringence on polarimetric microscopy imaging
The nonlinear polarimetric microscope used in this work has been described in previous works [23,26].The excitation light source is a tunable Ti:Sapphire laser (Chameleon, Coherent) that delivers 150 fs pulses at a repetition rate 80 MHz.The incident wavelength is set at 800nm (for SHG) and 840nm (for TPEF) with a typical averaged power of a few mWs.The laser beam is reflected by a dichroic mirror and focused on the sample by a relatively low Numerical Aperture (NA) objective (×40, 0.6 NA), in order to discriminate high aperture focusing effects from the present birefringence effect study.The backward emitted signal is collected by the same objective and directed to a polarization beamsplitter that separates the beam towards two avalanche photodiodes.Images are performed by scanning the sample on a piezoelectric stage, which allows precise location of polarimetric measurement points.The linear polarization of the incident laser beam is continuously rotated in the sample plane by an achromatic half waveplate mounted on a step rotation motor at the entrance of the microscope.For each value of the polarization angle α (relative to X), the emitted signal is recorded on the two perpendicular directions X (I X ) and Y (I Y ) defining the sample plane [Fig.1(a), 1b)].In order to measure locally the birefringence properties of the studied sample, we implemented on the same microscope a forward polarized detection identical to the backward set-up, the detected signal being the incident fundamental beam at the input ω frequency.Upon rotation of the incident linear polarization, the ellipticity occurring from the sample birefringence is expected to strongly modify the polarimetric dependence of the measured signal through the whole sample thickness L [Fig.1(c)].At last, polarization distortions dichroism and ellipticity occurring from the mirrors reflection are preliminarily characterized following a calibration procedure using a fluorescent solution as described in [29].
For both the TPEF and SHG contrasts, varying the incident polarization angle α allows retrieving information on the unknown molecular angular distribution denoted f (θ , φ , ψ), with (θ , φ , ψ) the Euler set of angles defining the molecular frame orientation in the sample.In what follows, we discard the ψ angle since cone-shape distributions are considered where only the (θ , φ ) angle components are relevant by cylindrical symmetry, as often applied to biological media such as in cell membranes [11] and fibrillar structures [2].Nevertheless this model can be extended to more general distribution functions.
The TPEF polarimetric signal measured from a molecular ensemble within an orientational distribution f (θ , φ ) (with (θ , φ ) the molecular dipole orientation angle), analyzed along a given polarization state I = (X,Y ), can be expressed as the incoherent sum over individual fluorescent intensities [23,25]: where μ abs is the molecular absorption dipole, and μ em is its emission dipole, both fixed relative to the molecular microscopic frame.E(α) the incident linear polarization state vector oriented with an angle α relative to X.This expression, written in the planar wave approximation, assumes a small collection angle of the fluorescence radiation.In the case of a high aperture collection angle, the | μ em (θ , φ ) • I| 2 emission contribution in Eq. ( 1) has to be replaced by the J I (θ , φ ) corrected function, previously introduced to account for the broad range of collection angles within a high numerical aperture (NA) objective [11].In the present work, the use of a low-NA objective (NA=0.6)allows considering the planar wave approximation to be relevant [29,33].
In the case of SHG, the signal results from the coherent radiation of nonlinear molecular induced dipoles, induced in the focal spot of the objective.This leads to the following SHG intensity measured along the I = (X,Y ) polarization analysis direction: with χ (2) IJK the macroscopic nonlinear SHG tensorial coefficients with (I, J, K) = (X,Y, Z).This expression is written in the planar wave approximation but can nevertheless be completed for high NA [31,33].The χ IJK coefficients are related to the angular molecular distribution in the focal spot of the microscope f (θ , φ ) following expressions developed previously [14,[20][21][22]25].
Note that due to the small size of the focal volume relative to the nonlinear coherent lengths, the phase matching relations occurring usually at a large propagation distance are not accounted for here.
In both expressions [Eqs.( 1) and ( 2)], the optical field E(α) is expressed at the dipoles positions (either fluorescent or nonlinear induced) in the focal spot of the microscope objective, at a given depth d from the sample surface [Fig.1(c)].In anisotropic samples where this incident field undergoes a birefringence retardation, its expression has to be re-written in order to account for the consequent polarization distortions.In the present work, we consider that the object projection in the sample plane is uni-axial, which is consistent with a cylindricalsymmetry distribution lying in the (X,Y ) plane as mentioned above.This is also relevant in most of the systems imaged in nonlinear microscopy.In the other cases, only a picture of the sample projection is possible since the optical coupling is limited to the X and Y polarization directions.Denoting Θ b the angle between the X macroscopic axis and the fast optical axis X b of this object [Fig.1(b)], and Φ b the phase shift between its fast and slow optical axes, then the new expression of the optical field polarization state at the focal depth distance d is, in the planar wave approximation: with E 0 X,Y the optical field polarization components in the macroscopic (X,Y ) frame at the sample surface (Z = 0) [Fig.1(c)], and Φ λ .Δn.d, with λ the incident wavelength and Δn the refractive index difference between the fast and slow axis of the object in the sample plane.
In addition to its effect on the incident excitation field, the birefringence also affects the detected signal, which propagates back in the sample in the case of a backward detection.A similar approach as above can be implemented to account for this effect, assuming that the same Δn value applies to both incident and emitted wavelengths (this is the case if no resonance is involved).The relation between the macroscopic emission dipole components (either occurring from a TPEF process or an induced nonlinear SHG process) at the focal depth d and at the exit of the sample (Z = 0) follows the same equation as in Eq. ( 3), introducing the detection wavelength in the expression of Φ Glass substrate E (Z=d)

In situ characterization of the sample local birefringence
Equation (3) shows that for a varying angle α of the incident polarization E(α), the sample birefringence parameters can be deduced from a polarimetric measurement performed in the forward direction for the ω frequency of the excitation field [Fig.1(a)].In such geometry, the field propagates through the whole thickness of the sample L and therefore Φ b (L) is measured.Theoretical polarimetric responses of the incident intensity I ω X (α) are depicted in Fig. 2(a  with a π periodicity, the periodicity being consistent with the fact that this technique does not discriminate the fast or slow axes of the system (this indetermination does not affect the results for the detected signals, as seen below).Note that accounting for the ellipticity and dichroism introduced by the dichroic mirror as mentioned in [29] leads to a slight modification of the fit solutions [Fig.2(d)] which has to be ultimately introduced in the models for TPEF or SHG contrasts as detailed below.Figure 2(d) shows in particular that introducing the instrumental distortions in the data analysis is crucial before any sample analysis.The birefringence parameters Φ b (L) and Θ b determined from the fit of [Fig.2(b)] shows that (i) Θ b lies roughly along the observed collagen fiber orientation (here 45 • ) and (ii) the measured birefringence phase shift Φ b (L) = 97 • can be used to deduce the approximate sample thickness L at the measured sample location.Assuming Φ b ≈1.35 • /μm (Δn≈0.003[37,38]) between the long axis of collagen fibers and their orthogonal direction, this measurement leads to L≈72μm, which is reasonable in the studied sample.
In the following sections, both TPEF and SHG contrasts are studied with known instrumental factors, Θ b and Φ b being introduced as fitting parameters in anisotropic molecular assemblies.Such systems can be viewed as an ensemble of excitation/emission dipoles lying within a cylindrical symmetry angular distribution f (θ , φ ).This distribution is most generally modeled by a cone shape of given aperture and mean direction angles.This model has been applied to a large number of situations in polarization resolved fluorescence such as in labeled lipid membranes [11], labeled Actin fibers [12] or disordered molecular crystals [28,33] as well as for SHG in collagen where the traditionally used C 6 crystal point group model leads to a nonlinear tensor equivalent to that of a cone [2,[20][21][22].In such uni-axial samples the molecular order imposes an anisotropy along the main axis that can lead to non-negligible birefringence.

TPEF polarization resolved microscopy in crystalline samples
In order to visualize the effect of birefringence on TPEF, a model is applied to a cone lying in the (X,Y ) plane, which main orientation ϕ 0 relative to X also corresponds to the optical axis of the uni-axial object.Figure 3 shows the effect of birefringence on the I X and I Y TPEF polarimetric response in a system of large cone aperture.An increasing birefringence leads to a deformation of the polarimetric polar plots, in particular when the optical axis is away from the experimental projection axes X or Y .This would clearly lead to a misinterpretation of the polarimetric data, even for slight birefringence phase shifts.In particular for large values of Φ b , the polarimetric data resemble those of a much smaller cone aperture than what is is actually.Experimental TPEF polarimetric data are shown in Fig. 4 on a 1D macroscopic uni-axial molecular co-crystal of Perhydrotriphenylene(PHTP)-4-Dimethylamino-40nitrostilbene(DANS) described in a previous work [27].In this crystal the birefringence is nonnegligible, which is particularly visible when the crystal is tilted [Fig.4(c)].The birefringence phase shift deduced from the fit of I ω Y is only indicative since a large thickness (about 500 μm) is crossed in the forward direction, and therefore many π phase shift periods are expected.To fit the TPEF polarimetric data, the molecular 1D distribution is modeled by a cone of infinitely small aperture angle (dirac function) oriented along the main axis of the crystal by an angle ϕ 0 relative to X.All angles ϕ 0 , Θ b and Φ b , are used as fitting parameters.The fit, which is clearly improved by accounting for birefringence [Fig.4(e), 4(f)], leads to values consistent with both the tilt angle of the crystal (ϕ 0 ≈135 • ) and a large birefringence.At a penetration depth of 5 μm, the fitted birefringence phase shift leads to Δn ≈ 0.04, which is relevant in such system where Δn values ranging between 0.001 and 0.8 have been reported in molecular crystals [36].Note that the measured birefringence value is averaged over the excitation volume of the objective, which introduces an error margin on the penetration depth value and therefore on the Δn determination.In addition, the difference between ϕ 0 and Θ b found when fitting I ω and I T PEF might be due to a slight heterogeneity of the crystal main axis orientation through its large thickness.Nevertheless, if the birefringence effect is not accounted for in the TPEF polarimetric fit, the quality of the fit is considerably reduced [Fig.4(f)].
Fit of the TPEF polarimetric response with no birefringence included.All the fits include the dichroism and ellipticity parameters of the dichroic mirror (dichroism γ = 0.009, ellipticity factor δ = 13 • with the mirror used for this experiment [29]).

SHG polarization resolved microscopy in collagen
In this part the effect of birefringence is explored on the SHG polarimetric data on collagen type I fibers from rat-tail tendons, attached at both ends on a coverslip and immersed in 0.15M NaCl.To prepare the fibers, Adult Sprague Dawley rats were euthanased for purposes unconnected with the present research.Tails were removed and immediately snap frozen in liquid nitrogen cooled isopentane.At the time of use the tissue was thawed and the tendon exposed.Individual fibers were teased out by microdissection and either examined immediately or stored frozen until required.Control Raman spectra were identical in either case, contained none of the peaks characteristic of proteoglycans, and were indistinguishable from those obtained from fibers purified by enzymatic extraction.
SHG imaging has been widely studied on similar systems.Analysis of SHG polarization resolved experiments is traditionally done assuming the collagen fibers to be of C 6 symmetry, equivalent to considering a collection of nonlinear individual dipoles within a cone distribution of typical half apertures angles ranging from 50 • to 65 • [15,16,18,21,22].Calculated SHG polarimetric responses are depicted in Fig. 5(a) for different tilt angle ϕ 0 of theoretical collagen fibers.The effect of a theoretical birefringence is seen to be detrimental to the polarimetric responses, introducing extra lobes in the polar plot shape comparing to the absence of birefringence.Again in this case, analyzing SHG polarimetry not accounting for birefringence leads to erroneous tensorial analysis of the SHG response.
Experimental analyzes performed in collagen type I fibers of thickness ranging from 70μm to 100μm confirm the strong deformation induced on the polarimetric responses, even for a few μm penetration depths in the sample [Fig.6 the X direction, two lobes appear clearly in the I SHG X polarimetric response, that cannot be explained in the absence of birefringence.The fit of the I ω X response [Fig.6(d)] confirms that the birefringence originates from a uni-axial system oriented close to the X axis (Θ b = 18 • ), with a phase shift Φ b = 119 • .This birefringence is consistent with Δn ≈ 0.003 [37,38] in a ≈ 88μm thick fiber.The fitting of the SHG polarimetric responses at different penetration depths are performed using Θ b and Φ b as fitting parameters for both I X and I Y .Excellent fitting quality is obtained for a cone distribution of 52 • half aperture angle, oriented at ϕ 0 = 10 • , which corresponds to the fiber direction observed in the SHG image (Fig. 6).The SHG fitting results show that the orientation of the optical axis Θ b = 2 − 9 • , measured locally, lies close to ϕ 0 .The difference with the averaged Θ b = 18 • measured over the whole sample thickness is most probably due to the sample heterogeneity in the Z propagation direction.At last, the measured birefringence phase shifts increase when further penetrating in the sample, as expected.Supposing a birefringence Δn ≈ 0.003, the deduced theoretical penetration depth (see values in brackets given in the Fig. 6 legend) are consistent with the expected depths of focus d, including error margins due to the focal volume averaging over several micrometers in the Z direction.
In the data shown in Fig. 6, the direction of the optical axis Θ b , measured locally, is still relatively close to its averaged value over the whole sample thickness.This is not always the case as shown in another area of the same sample (Fig. 7).In this situation the SHG polarimetric fit gives Θ b = (−23) • − (−12) • in the explored depths of focus [with consistent birefringence phase shifts, see Fig. 7(a)-7(c)], whereas Θ b = 36 • ≡(−54) • on average.This is a clear signature of a progressive rotation of the collagen fiber / optical axes orientation when penetrating in the sample depth.Note that still accounting for the variations of ϕ 0 seen on the SHG images does not notably change these results.
These results show finally, that particular care has to be brought on the local properties of the sample, which can be very distinct from averaged birefringence since the optical axis of the system can typically rotate progressively through the whole sample thickness.At last, a complete fitting procedure on both the cone shape and birefringence parameters at each measurement point is therefore necessary to lead to reliable data on the nonlinear tensor of the structure.These experimental SHG polarimetric data, performed at distances still relatively close to the sample surface, indeed clearly evidence the presence of birefringence and demonstrate the im- portance of accounting for such effect to retrieve sample structural information.

Conclusion
We have demonstrated that birefringence from molecular and biological samples is a detrimental factor in polarization resolved information in microscopy imaging applied to anisotropic samples.We have shown that this birefringence can be however fully characterized and accounted for.In a context where SHG microscopy constantly evolves towards the exploration of complex tissues where birefringence is present at increasing depths, this analysis is more and more crucial to distinguish the local nature (symmetry, orientational disorder) of molecular and biomolecular assemblies.This study can be extended to other contrasts such as Third Harmonic Generation (THG) or Coherent Anti-Stokes Raman Scattering (CARS) for which polarization-resolved microscopy has also been implemented.

Fig. 1 .
Fig. 1.(a) Experimental set-up.D: dichroic mirror, APD: avalanche photodiodes, P: polarizer, M: mirror, S : sample, FLP : long pass filter, FBP : band pass filter.(b) Definition of the input polarization angle α and the birefringence slow axis direction θ b relative to X. (c) Experimental configuration: an incident linear polarized field at the ω frequency is focused at a distance d from the sample surface in the optical axis Z direction.The transmitted fields at ω at both the distance d and the sample output distance L are elliptical if the sample is birefringent.The emitted field (fluorescence or 2ω for SHG) is measured in the epi direction.
) in a polar representation, for a sample tilt angle of the optical axis of Θ b = 30 • , and for different birefringence phase shifts.If no birefringence occurs in the sample, the polarization response is a cos 2 α two-lobes pattern with a maximum intensity along X (α = 0 • ) and a vanishing intensity along Y (α = 90 • ).In the presence of birefringence, this pattern tends to open (no extinction occurs anymore) and rotates progressively when Φ b increases.Similar responses are expected for I ω Y , rotated by 90 • with respect to I ω X .As seen in Fig.2(b), experimental I ω X,Y polarimetric measurements performed at a given position on a collagen type I fiber of roughly 70-100μm thickness exhibit such behavior (the sample preparation details can be found in the SHG section below).

Figure 2 (Fig. 2 .
Fig. 2. (a) Theoretical linear polarimetric responses I ω X through a birefringent sample at different values of Φ b , for an optical axis tilt of Θ b = 30 • .(b) Experimental linear polarimetric responses I ω X,Y (dots) through a collagen type I fiber from a rat-tail tendon, oriented at 45 • relative to X.The fit (black lines) uses the model detailed in the text.(c,d) Maps of the mean square error between fitted and experimental intensities : (c) with no dichroic mirror polarization distortion, (d) with polarization distortions (dichroism factor γ = 0.043, ellipticity factor δ = 55 • , with the notations of [29]).Fit solutions : Θ b = 44 • and Φ b = 97 • .

Fig. 5 .
Fig. 5. Theoretical SHG polarimetric responses (I SHG X (red) and I SHG Y (green)) for a collagen fiber which nonlinear tensor corresponds to a cone half angle aperture distribution of the nonlinear molecular dipoles of 50 • .The orientation of the fiber ϕ 0 also corresponds to the optical axis orientation Θ b . (a) ϕ 0 = Θ b = 0 • and (b) ϕ 0 = Θ b = 45 • , with different values of the birefringence phase shift Φ b .

Fig. 6 .
Fig. 6.Experimental SHG polarimetric responses for collagen type I fibers oriented close to the macroscopic X axis.The SHG images (image scale is in counts/s) represent a 30μm×30μm area of the sample scanned at different depth d: (a) 2μm, (b) 20μm, (c) 60μm.Both SHG experimental intensities I X (red dots) and I Y (green dots) are depicted (after normalization for more visibility), together with the corresponding fit according to the model detailed in the text.The parameters obtained from the fits are 52 • and ϕ 0 = 10 • for respectively the cone distribution orientation and aperture, and : (a) Θ b= 7 • , Φ b = 19 • [14μm];(b) Θ b = 9 • , Φ b = 29 • [22μm];(c) Θ b = 2 • , Φ b = 88 • [65μm].(d)Birefringence measurement and fit on the I ω X component of the incident field, leading to birefringence parameters of Θ b = 18 • , Φ b = 119 • [88μm], when the sample is propagated through.The found birefringence axes are represented in the inset scheme.