Abstract

The second harmonic generation microscope (SHGM) constructs images of intensity distributions of SH waves produced by the interaction of fundamental waves with a polar material. We have developed this nonlinear optical microscope in order to make possible nondestructive, three-dimensional (3D) observations of various kinds of inorganic and organic materials. The SHGM can disclose also inverted domain structures of antiparallel spontaneous polarizations using the interference with the reference SH waves. The observation principle and several applications to structural characterizations of LiNbO3 and LiTaO3 quasi-phase matching devices, domain structure analyses of a relaxor/ferroelectric solid solution Pb(Zn1/3Nb2/3)O3-9%PbTiO3 at the morphotropic phase boundary, development of order parameter in a quantum paraelectric relaxor Li-doped KTaO3, and antiphase polar domain structures of muscle fibers and myofibrils are surveyed by stressing the high effectiveness of the SHGM as a tool of material diagnosis.

1. Introduction

Optical second harmonic generation (SHG) is a phenomenon where coherent light with twice the fundamental frequency is generated through the nonlinear interaction between a laser and a material. Since SHG is sensitive to the break of spatial inversion (SI) symmetry or time inversion (TI) symmetry, it provides a powerful experimental technique for investigating ferroelectricity related to the SI symmetry break [13], magnetic ordering related to the TI symmetry break [4, 5], and multiferroics related to both phenomena [69].

Certain kinds of materials are classified as ferroics where spontaneous electric polarization (ferroelectrics), spontaneous magnetic order (ferromagnetics), or spontaneous strain (ferroelastics) can be reversed by their conjugate external fields , , or , respectively. Ferroics have a phase transition temperature . Above these properties disappear, while below , the spontaneous order parameters develop in macroscopic or mesoscopic regions having different orientations but same magnitude of spontaneous values. These regions are termed domains or variants. Spontaneous order parameter displays a hysteresis loop upon applying these conjugate external fields and these phenomena are the process of changing in domain states from multiple to single domain and vice versa. Recently much attention has been paid to cross-correlations of and in multiferroics and a lot of studies have been performed in this promising field.

In the case of ferroelectrics, the development of can be detected by SHG in a wide dynamic range to determine the critical index at the vicinity of the phase transition point [10, 11]. SHG has also been used to observe domain structures by exploiting the spatial distributions of the SHG tensor component [1215]. Based upon this advantage of SHG for material characterizations, we have developed the SHG microscope (SHGM) for material characterization, in particular, nondestructive and three-dimensional (3D) observations of domain structures of inorganic, organic, and biological materials. This paper reviews the principle of our SHGM and its application to observations of domain structures in certain kinds of ferroelectrics and biological tissues of muscles as well. Section 2 describes the principle of SHG (Section 2.1), the principle of SHG interference microscopy for observing antiparallel domain structure (Section 2.2), an intrinsic issue ( issue) for realizing 3D domain structures using the strongly focused scanning-type SHGM (Section 2.3), and the polar diagram representation of SHG images and related problems (Section 2.4). Section 3 explains the optical system of our microscope. Section 4 shows some examples of SHG images and their analyses picking up 4 cases: 3D observations of periodically poled domain (PPD) structures in quasi-phase matching (QPM) devices for wavelength converters of lasers (Section 4.1), complex domain structures of a relaxor/ferroelectric mixed solution Pb(Zn1/3Nb2/3)O3-9%PbTiO3 (PZN/9PT) at the morphotropic phase boundary (MPB) (Section 4.2), a quantum paraelectric relaxor TaO3 (KLT) (Section 4.3), and muscle fibers and myofibrils (Section 4.4). Section 4 sums up the results and discusses future trends of SHG microscopy.

2. Theoretical Treatment of SHG

2.1. Principle of SHG

As a result of the nonlinear coupling of fundamental wave fields and materials, the following nonlinear wave propagation equation (1) is derived from the Maxwell equation. SHG is a case of the second-order nonlinear optical process where is proportional to the square of the fundamental wave fields. As a result, coherent waves with the frequency twice, therefore, the wavelength half of the fundamental waves are generated: Here the right hand terms express the light source of harmonic waves: is the nonlinear polarization induced by the fundamental laser beam, is the nonlinear magnetization, is the electric quadrupole moment, and and are dielectric and magnetic susceptibilities in a sample, respectively. In general cases, the first term is dominant in the SHG, and we consider only this term in the following. The left hand terms express the propagation of waves generated by the right hand terms.

and the electric field of the fundamental wave are related through the SHG tensor component (-constant) as is expressed in the following equation: Here and mean the third-rank tensors and is invariant and is variant with respect to the time inversion operation. These -constants are related to and magnetic order parameter as The relation of these constants to the existence of spatially inversion and magnetic inversion symmetries is summarized in Table 1. Since SHG intensity is proportional to the square of -constant, the behavior of structural and magnetic order parameters and the correlation between them can be obtained from SHG measurements. Hereafter we describe phenomena related to in nonresonant wavelength region, that is, the case of noncentrosymmetric and nonmagnetically ordered transparent materials.

For getting the SH wave field from (1), we consider the collinear case where the fundamental wave with TEM00 Gaussian mode with a beam radius at the focal position and the Rayleigh length and SH waves propagate along the axis. The complex field amplitude of the SH wave generated in a sample is obtained as [16] Here is a measure of the wave number offset between the fundamental () and SH waves () and is expressed by is the sample thickness and is the distance of the focal position from the incident surface of the sample. is an absorption parameter expressed as , and and are the intensity absorption coefficients of the fundamental and SH waves, respectively. In a sample with multidomain structures, , , , and are dependent on the spatial coordinates ().

2.2. Principle of SHG Interference Microscopy for Observing Antiparallel Domain Structures

A simple and general method of observing domain structures by SHG is to exploit the anisotropy of  -constants. For example, 90° domain structures with directed along the , and axes in ferroelectric tetragonal BaTiO3 can be observed using the difference in magnitude of and as shown in Figure 1 [12]. However, 180° domain with antiparallel cannot be distinguished, because SH waves have same amplitude for both domains. Other optical method cannot also detect this domain structure. The nondestructive observation of 180° domain structure has been required in particular to characterize precisely QPM devices. This motivated us to develop the SHG interference microscope (SHGIM) that makes the observation of 180° domain structure possible [12, 13, 17, 18].

The SHG interference microscope (SHGIM) exploits the interference of SH waves generated in a sample and a reference plate to observe 180° domain structures in ferroelectrics. The principle is schematically illustrated in Figure 2, where a sample contains two domains with opposite polarities. In this case, the signs of -constants are opposite and the phase difference of SH waves generated in these two domains is . Therefore, with uniform reference waves, the interference produces an intensity contrast. This principle of SHGIM can be described analytically as below in a strongly focusing optical system, such as the present optical system.

Equation (4) is modified when up and down domains exist in a sample as The interference intensity of   and the reference SH wave are described by Here, is the phase difference between the fundamental and SH waves. Because is proportional to , the interference intensities, the third term in the right-hand side, of antiparallel domains are different, and the 180° domain structure can be observed as the intensity contrast. The phase is varied to obtain the maximum intensity contrast by a phase modulator.

2.3. Principle of Observations of Internal Domain Structures: Issue

Our SHGIM enables us to observe polar domain structures in cross-sections of a sample. For this 3D observation, we adopt a confocal optical system composed of a pair of objective lenses with high numerical aperture (NA) and a pinhole (see the optical system of SHGIM in Section 3). The optical diffraction limit at the focal position and the Rayleigh length are measures of the spatial resolutions along the lateral and axial directions, respectively, and both characteristic lengths are related to NA as

Here, is the wavelength of an incident wave, and is the refractive index of a sample at . Since and are proportional to (NA)−1 and (NA)−2, respectively, a higher NA yields higher spatial resolutions. To describe the degree of focus, a nondimensional parameter is introduced: It should be noted that the strongly focused optical system with high NA yields phenomena substantially different from a weakly focused one, which affects the visibility of domain structures on the inside of a sample [19].

In the weakly focused case , the plane wave approximation is valid and the SH intensity is expressed as Equation (10) indicates that does not depend on the sign of . However, when becomes larger, cannot be approximated by (10) and depends strongly on the sign of and the focus position as shown in Figure 3 [19]. Let us call this phenomenon “ issue,” stressing the important observation criterion of the internal domain structures using SHGM. This simulation result is calculated using (4) with similar magnitudes as in the present experiments. The SHG intensity decreases drastically inside the sample when   changes from negative to positive. On the other hand, when the focus position is located in the vicinity of end surfaces, remains large and decreases slowly with . This phenomenon induces the degradation of the visibility of inside structures in the case of positive [16] and the use of -constants with negative is strongly requested.

The physical origin of the issue is explained as follows. Generally, the phase of Gaussian wave changes from to upon propagation. This phase is known as the Gouy phase and is largely changed at the focus point: the more strongly focused is the wave, the larger is the change in the phase at the focus position. This phase change makes the effective phase   of SH waves toward the positive direction and makes away the phase-matching condition. On the other hand, negative of the sample pulls the effective phase toward which can satisfy the phase match condition to obtain large SHG intensities.

2.4. Polar Diagram Representation of SHG Images

For quantitative analyses of SHG images, 2D maps of polarization dependences of SHG intensity in the polar diagram (PolD) are measured for each area in mesh under the condition that the polarization directions of the incident and SH waves are kept parallel to each other [20]. The result is fitted with calculated values that are obtained from the domain configurations. The following equations show the case of PZN/9PT under the assumption of the monoclinic space group Pm. In this calculation, we adopt the laboratory coordinates () parallel to the cubic axes with the fundamental wave incident along the axis and the SH intensities produced in the plane. The SHG intensities are expressed as a function of the polarization angle :Type 1:Type 2:Type 3:

In the monoclinic system, the twenty-four spontaneous polarization states exist and the observed PolDs are classified into three groups, expressed by (11)–(13). Types 1–3 correspond to the classification of Table 2, and the form of PolD is determined by the direction of projected onto the plane. The results of PZN/9PT are discussed in Section 4.2.

3. Optical System of SHGIM

The optical system of the SHGIM is illustrated in Figure 4. The light source is an Nd : YVO4 laser with a wavelength of 1064 nm, a pulse width of 10 ns, a repetition rate of 40 kHz, and a maximum power of 6 W. The fundamental beam is expanded by a couple of lenses (2a and 2b) and divided into two beams by a separator . One is used for changing the phase of the fundamental wave by a phase modulator , and the other is for producing the reference SH waves by a QPM wave converter . Intensities of the fundamental and reference SH waves can be changed independently by attenuators composed of a half wave plate and a thin-film polarizer (4 and 4′). The polarization direction of the incident waves is changed by a half-wave plate . Fundamental and reference SH waves are focused to a sample by an objective (12a) with a numerical aperture (NA) of 0.7. Average fluence of the laser beam at the sample position is less than 30 J/cm2 so as to produce no optical damage in samples. Total SH waves generated from the sample and the reference plate pass through a rear objective (12b), and the polarization direction of the SH waves is selected by an analyzer . The light with wavelength 532 nm is selected by a spectrometer , and the intensity is detected point by point by a photomultiplier tube that is synchronized with incident laser pulses by a lock-in amplifier . The sample stage is driven by the combination of stepping motors and piezoactuators along the three mutually perpendicular directions. Sectional SHG images are displayed by a computer . With this microscope we can select the noninterference geometry simply by shutting off the second path to the reference plate.

We also constructed an SHGIM equipped with an image-intensified CCD. In this system, the 2D SHG image of a sample integrated along the incident direction is observed. This has not high spatial resolution such as the scanning microscope but has a merit of higher temporal resolution of several seconds. We applied this ICCD SHGIM for observing development of the order parameter of KLT (see Section 4.3).

4. SHG Images and Domain Structure Analyses

4.1. Periodically Poled Domain Structures in Quasi-Phase Matching Devices

As a first example of SHGM observations, we show here nondestructive and 3D characterizations of PPD structures in LiNbO3 and LiTaO3 QPM devices for wavelength converters of laser. The wavelength conversion technique is highly demanded in the optoelectronics for obtaining short-wavelength laser sources or optical parametric oscillators. To obtain high conversion efficiency, it requires the phase coincidence of optical waves generated in a nonlinear medium. This is termed phase matching techniques, among which the QPM is most promising because of flexibility of the incident wavelength and possible use of the maximum SHG tensor component. Since the efficiency depends strongly on the exact periodicity of the PPD structure in a designed volume area, nondestructive and 3D observations have been required in place of a chemical etching technique which damages devises and provides only domain structure at sample surfaces. However, only two methods are reported: one is the two-wave mixing technique based on the photorefractive effect [21] and another uses four SHG interference images with different phases [22]. Both methods have drawbacks: the former method depends on the specific optical property of a sample which is undesirable as a stable wavelength converter, and the latter needs complicate process of analyses.

We apply the strongly focused scanning SHGIM to the observation of cross-section images of the PPD structure in an MgO: stoichiometric LiTaO3 (MgSLT) QPM [17]. In this study, however, we observe clear images of the PPD structure with the period of 8 μm only near surfaces. This is because we use the component whose is positive. On the other hand, when we observed the PPD structure of a QPM devise fabricated into MgO-doped LiNbO3 using the component with negative , the PPD structure of a designed 19.5 μm period is clearly observed inside the sample. The result is shown in Figure 5.

4.2. Domain Structure Analyses of PZN/9PT

High piezoelectric constants and mechanical coupling coefficients are observed in solid solutions of relaxors and ferroelectrics at specific composition ratios. This compositional phase boundary is termed the MPB [23]. The symmetry of the MPB has been found to be lower than those of the two end members (trigonal in relaxor-rich side and tetragonal in ferroelectric rich side) [24, 25]. This symmetry lowering to the monoclinic phase makes the rotation of more facile, which induces greater piezoelectricity [26]. The symmetry of its MPB has been investigated using several kind of experimental methods without the definite conclusion. This is the motive of our study of domain structure analyses of PZN/9PT at the MPB using SHG microscopy. We observed the domains with the spatial resolution of the optical diffraction limit of about several hundred nm, and the symmetry was determined by analyzing the 2D polar diagram maps [20]. In this study, we cannot observe 180° domain structures because we do not use the SHG interference technique.

Figure 6(a) shows an SHG image of the PZN-9PT (010)c plate. Bright and dark parts represent the area generating strong and weak SH intensities, respectively. A photograph of the same area taken by the polarization microscope in the crossed Nicols is also shown in Figure 6(b) for comparison. Similar domain structures are observed in both figures. However, the SHG image provides more clearly the domain structure as bright and dark stripe patterns as shown in this figure, and furthermore, more quantitative information can be obtained from SHGM as described below.

Figure 7(a) shows SHG images of the region A taken with finer scanning step and Figure 7(b) shows the PolD maps corresponding to Figure 7(a). It is obtained with parallel P and A configuration. Two patterns exist indicating blue and red lines, and these two patterns almost exactly coincide with the SHG image of Figure 7(a). PolDs calculated with (11)–(13) assuming monoclinic    are shown in Figure 7(c). Experimental results are in good agreement with the theory. Furthermore, the direction of can be estimated by the polar diagram and are shown in Figure 7(c). It is also revealed that the direction of is not located in this plane and the arrows represent the projection of on the (100) plane.

The same procedure is applied to region B. In this case, only one pattern was observed in the polar diagram, and the form can be also explained with monoclinic symmetry. Full analyses of the PolD maps reveal the complicate domain structures as shown in Figure 8. Thus we conclude that the macroscopic symmetry of the MPB of PZN/9PT is monoclinic . We would stress that the mapping of the polar diagrams can provide useful and quantitative information on the complex domain configuration.

4.3. History Dependence and Time Evolution of the Order Parameter in KLT

KTaO3 is known to be quantum paraelectrics where a ferroelectric phase transition does not take place although a dielectric constant increases with decreasing temperature and attains high value at low temperature. With a partial substitution of Li ions at K sites, a dielectric peak appears at low temperature and the peak temperature shifts to higher temperature with increasing Li concentration. However, the polar state at low temperature had not been well recognized [27, 28]. Yong et al. carried out dielectric and neutron diffuse scattering experiments and claimed that KLT is a relaxor because of the dielectric dispersions and diffuse scatterings related to the existence of polar nanoregions (PNRs) [29]. However, more convincing experimental results are needed on this point. To clarify the relaxor nature, we performed precise SHGM observations of KLT [2, 3].

We carried out the experiments using ICCD SHGM which enables to observe the 2D polarization distributions in a sample from the SH wave intensity distribution. The measurements were performed in the following five processes; (i) zero field cooling (ZFC), (ii) zero field heating after ZFC (ZFH/ZFC), (iii) field heating after ZFC (FH/ZFC), (iv) field cooling (FC), and (v) FH after FC (FH/FC). The electric field of 100 V/mm was applied along the direction and obtained results are shown in Figure 9. In ZFC and ZFH/ZFC processes, no marked SH intensity is observed in the measurement temperature region. In FC, SH intensity starts to increase drastically at a certain temperature. The temperature at which SH intensity begins to increase shifts to higher temperature with increasing Li concentration. In FH/FC process, SH intensity generated from the specimen at low temperature decreases with raising temperature and vanishes around a temperature . In FH/ZFC, SH intensity dose not appear at low temperature even applying an electric field. With heating the sample, a weak SH intensity observed at low temperature starts to increase abruptly at certain temperature (42 K for KLT-2.6% and 83 K for KLT-6.8%), shows a peak and decreases, and finally vanishes at a temperature . coincides with where SH intensity disappears in FH/FC process. Based upon these results, we confirm the marked history dependence of the order parameter in KLT, which is a typical feature of relaxors.

In ferroelectrics, the domain switching time is known to be several  sec. However, in relaxor, a relaxation time to the equilibrium state under external stimuli is several ten minutes. To observe this slow kinetics of the order parameter in KLT, SHGM observations were carried out in the following thermal procedure [30]. First the specimen was cooled down to 24 K without an electric field and heated up without an electric field to a certain temperature. When a temperature becomes stable, the electric field was applied and the SH intensity was measured as a function of time. Figure 10 shows the time dependence of the SH intensity at several temperatures observed in KLT-2.8%. The present results indicate that the slow kinetics observed in SH intensity strongly depends on the temperature where an electric field is applied. The time evolution of SH intensity shows a drastic difference around 40 K. Below 40 K, SH intensity develops after a long incubation time. On the contrary, SH intensity shows a fast evolution followed by the slow process above 40 K. This phenomenon could be explained by introducing a long-range strain field proposed by Chandra, who explained a long domain switching time observed in BaTiO3 [31]. In the case of KLT, off-centered Li ion could produce the long-range strain field. Thus above-mentioned two experimental results concerning the behavior of the order parameter, that is, the history dependence and the time evolution with a long relaxation time, confirm that KLT is relaxor.

4.4. Muscle Fibers

Muscles are functional body tissues and have structural hierarchy as shown in Figure 11: a muscle consists of bundles of muscle fibers of several ten μm diameter and a muscle fiber is an assemblage of myofibrils of several μm diameter. A myofibril is composed of a one-dimensional sequence of sarcomere with the period of about 2 μm along the fiber direction. A sarcomere is a unit structure of muscle and contains thick myosin filaments, actin filaments, and titins that connect myosin filaments with sarcomere wall (-line). Contraction and elongation of muscles are caused by relative sliding of actin and myosin fibers. X-ray and electron microscope revealed that a myosin filament is formed by coiled coil dimmers characterized by two myosin heads. The myosin chain is a coiled portion of two alpha helices, and SH waves are generated by the asymmetric charge distribution in hydrogen bonds in alpha helices. It is interesting to note that this spiral structure is opposite with respect to M-band which separates the sarcomere into two parts. Thus it is expected to have antiparallel polarization structure, or we can say, anti-ferroelectric structure in mesoscopic scale. However, no observation has been performed on this possible polar structure as yet. This is our observation motivation of myofibril using our SHG interference microscope.

For SHGM observation, we used white skeletal muscles, or psoas muscle of rabbit. Figure 12(a) shows an SHG image of muscle fibers. It is shown that muscle fibers are formed by regular assembly of the bright part. Figure 12(b) shows an SHG image of a myofibril extracted from muscle fibers by centrifuge. The periodicity of the bright parts is exactly same as that of muscle fiber, and finer scanning reveals 2 μm periodicity of the bright parts. This fact indicates that the bright part corresponds to the A-band in sarcore. It is also interesting that muscle fibers shown in Figure 12(a) are a coherent bundle of myofibrils.

To reveal the polarity of sarcomere, we exploit the SH interference. The result is shown in Figure 13. Here the relative phase of SH waves from the sample and the reference plate is changed by a phase modulater. The most left image shows SHG image without interference, the next one from the left is an image with phase difference 0, then phase shift  , and  . The position of the brightest part periodically shifts depending on the phase with periodicity of  . This means that the sarcomere consists of two parts with opposite polarities, and as a total, the myofibril has a periodically inverted polarization structure. It should be pointed out that the dimension of inverted part is about 1 μm, which is almost the limit of spatial resolution. Therefore the two parts with opposite polarization cannot be observed separately.

It is interesting that nature produces PPD structure in muscle fibers. If we consider this fact from an engineering perspective, it is really the QPM device, when the laser beam is introduced along the fiber axis. The periodicity of 2 μm could provide the QPM device at the violet wavelength region. Actually the diameter of muscle fiber of 100 μm is almost similar to an optical fiber and that of myofibril of several μm to a core of an optical fiber.

5. Summary: Future Trend of the SHG Microscopy

Due to the high sensitivity to the breaks of spatial inversion and time reversal symmetries, it is more and more recognized that the SHGM becomes a powerful method to investigate the nature of ferroic materials, that is, the behavior of order parameters, domain structures, symmetry changes, and so forth. It should be also pointed out that SHGM has been adopted not only for the material science but also for the biological and medical field [32, 33].

One of the bottle necks of SHGM is its spatial resolution. Using the conventional optical geometry, the spatial resolution is limited by the diffraction limit that is about the wavelength (several hundred nm) of the incident waves. One way to overcome this limitation is the use of a near-field microscope (SNOM). The SNOM has a long research history and the basic idea was already proposed by Synge in 1928 [34]. Ash and Nicholls demonstrated experimentally the super-resolution at microwave frequencies [35]. Since then, the SNOM has been extended to the visible wavelength region and the resolution reaches 25 nm [36].

The techniques of SNOM can be divided into two categories dependent on a probe. One is aperture mode whose probe is used either to deliver the incident wave or collect the scattered light. The spatial resolution of this system is determined by the size of the aperture. Another way is called as an apertureless mode. This technique utilizes the near-field waves scattered by a probe and the scattered wave is detected by a photomultiplier.

In recent years, SNOM has been also applied to the SHG. The first experimental result was reported by Smolyaninov et al. concerning the domain structure on Ni single crystal [37]. They used the aperture near-field scanning technique with an uncoated tapered fiber tip as a probe and achieved the resolution around 150 nm. The apertureless SNOM SHG system was also constructed based on the atomic force microscope and realized fast-time resolution measurement [38]. The measurement of local variations of -constants components is also reported [39]. The SNOM-SHGM is promising technique with high spatial resolution to investigate nanosize domain structures in thin films and the structural and physical properties of domain boundaries. However, extremely low intensity of the signal causes the low signal to noise ratio. Further improvements are required to establish the SNOM-SHGM as a universal probe microscope.

Acknowledgments

The authors are grateful to Professor S. Ishiwata, Waseda University for preparing living samples of muscle fibers and myofibrils and a lot of valuable comments on the SHGM observations of them. This work is partly supported by the Grants-in-aid for Scientific Research (C) of MEXT, Japan and the special grants for research (A) and (B) in Waseda University. J. Kaneshiro is grateful for the support of Grants-in-Aid for JSPS Fellows. H. Yokota is thankful to the support of Career-Support Program for Woman Scientists at Chiba University.