Nonlinear Diffuse fs-Pulse Reflectometry of Harmonic Upconversion Nanoparticles

Nonlinear diffuse femtosecond-pulse reflectometry is introduced as a powerful experimental tool for the unambiguous characterization of polar and non-polar point symmetry groups of harmonic upconversion nanoparticles. Using intense ultrashort 40 femtosecond laser pulses and an appropriate figure of merit (FOM), second and third harmonic emission serve for the structural characterization of polar Yb-doped lithium niobate and non-polar titanium dioxide nanoparticles. The tool is capable of differentiating these two samples by FOM values that differ by up to 13 orders of magnitude. The general applicability to harmonic upconversion nanoparticles over a broad range of intensities and wavelength spectrum, is discussed.


Introduction
The class of upconversion nanoparticles (UCNP) comprises nanoscopic crystals that are capable of converting two or more photons into a single one with higher photon energy, thus enabling light emission with large anti-Stokes shift with respect to the pump beam frequency.The most prominent examples are upconversion nanocrystals [1], where light emission is achieved via sequential absorption to intermediate states [2] featuring sharp emission spectra and long-lived excited states.These exceptional properties have already led to a diversity of applications and a strong focus in the realm of bioapplications such as imaging, sensing, and cancer therapy [3][4][5][6][7][8].
A new type of UCNPs-so-called harmonic upconversion nanoparticles-increasingly attract attention in recent years in material sciences in the framework of the search for UCNPs with novel properties and the exploitation of these properties for emerging applications.In harmonic UCNPs, the underlying mechanism for light emission is based on a (quasi-)instantaneous emission by means of the nonlinear optical principle of frequency conversion [9].Thereby, emission throughout the visible and near-infrared spectral range can be realized using the same nanoparticle by simply tuning the photon energy of the pump only.The possibility of even generating deep-ultraviolet light may potentially revolutionize bioapplications, particularly in the field of optogenetics [10][11][12].In order to clearly mark the difference to fluorescent mechanisms, the term "harmonic nanoparticles" is also frequently used [13,14].However, since all of the aforementioned optical processes yield an increase of emitted frequencies compared to those of incident waves, we will henceforth categorize them as "harmonic UCNP".
Accompanying this rapid development in material sciences, there is a need for an accurate characterization of harmonic UCNPs by means of appropriate spectroscopic tools capable of detecting respective nonlinear harmonic fingerprints.The first and most important step in material development is the verification of the nonlinear optical conversion mechanism, which is directly dependent on the structural feature of the as-synthesized nanoparticle sample.In particular, the mechanism of harmonic generation is highly sensitive to the presence of a polar or non-polar structural phase, so harmonic light emission can serve as a sensitive tool for structural characterization.For bulk crystals, appropriate tools were developed in the 1960s, and are successfully and widely applied in the framework of crystal growth with micro-and nanocrystalline materials.As a matter of fact, second harmonic generation (SHG)-where two pump photons are upconverted to a single photon-has been developed as an important complementary tool to XRD, especially if the structural analysis by the latter yields ambiguous results [15,16].In the small particle limit, however, present techniques may encounter difficulties due to diminishing nonlinear harmonic efficiencies [17], especially for higher harmonics beyond the second order, thereby enforcing intensities in the PW/m 2 -region.Combined with pulse durations of Q-switched lasers on the order of (1-10) ns, exposure leads to optothermal breakdown exceeding the damage threshold by several orders of magnitude [18], thus causing particles to melt, conglomerate, or become ablated.
A possible solution is to use pulse durations well below the nano-and even picosecond region [18]-in particular, femtosecond (fs) pulses, which are the current state-of-the-art .Nevertheless, the interaction of light pulses with nanocrystalline materials in a liquid solution as proposed for ns-setups [19] differs fundamentally for ultrashort lasers due to the particular properties of fs-pulse-matter interaction.For instance, pulse broadening via group velocity dispersion (GVD) may reduce the peak intensity and therefore the harmonic signal's intensity as well.Similar detriments stem from white light generation, harmonic generation, and multiphoton absorption from surrounding liquids, cuvettes, or cover slips.Still, these drawbacks may be circumvented in such transmittive setups [19,20] through various procedures; for example, GVD may be reduced through chirp compensation.In order to increase the overall nonlinear signal, though, large quantities ( 1 mg) of material are required, thereby favoring investigation on dry powders.
In the present paper, we propose a novel tool for the analysis of harmonic UCNP-powders that is able to overcome these disadvantages and focuses on the (quasi-)instantaneous upconversion through harmonic generation by means of nonlinear diffuse fs-pulse reflectometry.The tool is based on the concept of diffuse reflectometry, particularly with respect to the sample preparation.Using a tunable 40 femtosecond laser light source, it is shown that second and third harmonic signals can be obtained at high signal-to-noise ratios.The respective signals in turn are used to define a figure of merit that allows for the unambiguous assignment of the polar or nonpolar phase, respectively.Further discussion is devoted to the error susceptibility of such a figure of merit.The applicability and reliability of the introduced tool is demonstrated with a comparison of polar Yb-doped lithium niobate and non-polar titanium dioxide as an example.The results are discussed with respect to the impact of intensity and wavelength range, thus highlighting the general significance of harmonic upconversion nanoparticles.

Nonlinear Optics
Upconversion in nonlinear optical materials is commonly described through the coupled-wave theory, considering photon energy and momentum conservation respectively, with the summation carried out over the photon index i.Within the dipole approximation, expansion of the electric polarization P in a Taylor-series with respect to the electric field E yields with the susceptibilities χ (n) .The solution of the nonlinear coupled-wave equations aims at the description of the conversion efficiency η, which is the ratio of generated intensities to incoming intensities.According to Equation (1), it is mainly dependent on the susceptibility tensors χ (n) , whose even-orders-depending on the internal point symmetry group (PSG)-may vanish in the dipole approximation.A direct consequence is that centrosymmetric media do not exhibit significant even-order nonlinear effects (i.e., SHG), which are only weakly contributed by the symmetry-break at the crystal's surface [21][22][23].This signal break-down is the starting point of the unambiguous structural characterization as implemented in our introduced tool.
In case of bulk SHG, the nonlinear coupled-wave equations can be solved analytically under the non-depleting pump and slowly varying envelope approximations, yielding the well-known SHG-intensity I 2ω as a function of the fundamental intensity I ω [24]: Here, l is the interaction-length through the medium, n ω,2ω are the refractive indexes for fundamental (ω) and SHG (2ω) wave, and λ ω is the wavelength of the fundamental beam.The phase-mismatch is taken into account through ∆k.
An extension of the coupled-wave equations gives rise to third harmonic generation (THG), whose intensity follows the relation [25]: It should be noted that the THG phase-mismatch ∆k = k 3ω − 3k ω has already been neglected as will be motivated by the end of this section.An alternative to the direct third harmonic generation can be achieved in polar media through cascaded second-order processes where SHG is sum-frequency mixed with the remainder of the pump wave [25][26][27]: In contrast to even-order susceptibilities, the χ (3) -tensor will always contain nonzero elements, regardless of its PSG.This way, THG should be attainable from any material, yielding a third harmonic intensity I 3ω = I direct 3ω + I casc 3ω , whereas the detection of the opposite to an appreciable degree-that is, SHG in inversion symmetric materials-is strongly suppressed.
Introducing nanopowders in lieu of macroscopic bulk crystals allows significant simplificiations to Equations (2)-(4).Due to the statistical orientation of individual nanocrystals, angular relations as dictated through the nonlinear susceptibility-tensors will be averaged out, resulting in angle-independent scalar values for χ (n) .On the other hand, the phase-mismatch ∆k may be neglected due to the particles' dimensions with the interaction length l → 0, thus rendering sinc-terms approximately unity.

Introducing a Figure of Merit
Based on the strong relation between second harmonic signal and PSG, we introduce a figure of merit (FOM) composed of the measured signals S nω of SHG and THG in order to evaluate the crystal's polarity.This value is defined through the integrated spectral intensity distribution I nω (E ph ) recorded at the harmonic peak (n = 2 for SHG, n = 3 for THG) after baseline-correction.In other words, I nω (E ph ) is the emission spectrum of the n-th harmonic generation recorded by a spectrometer, whereas S nω is its integral.
Furthermore, it proves beneficial to introduce a figure of merit as the ratio between the SHG signal S 2ω and its THG counterpart S 3ω , in the following, termed "harmonic ratio": which yields a helpful indicator whether a material is polar or not.For nonpolar media, THG emission generated from the volume can be expected to exceed the surface-produced SHG by several orders of magnitude as a result of a higher number of dipoles involved in the former compared to the latter ( f R 1).In the opposite case where χ (2) does not diminish, SHG emission is predicted to be larger than or at least on the order of THG intensity, and therefore f R 1.
It is necessary to highlight the action of f R in more detail.In particular, it must be noted that this figure of merit inherently favors SHG over THG by a small margin, which can be typically neglected.This is the result of different peak widths of the harmonics and their respective exponents.More precisely, THG is broader than SHG on an energy scale.At ideal conditions where both peaks are perfectly Gaussian, the ratio between the peak widths of SHG and THG is √ 2 : √ 3.In case both peaks are equally tall with a peak spectral intensity A, the harmonic ratio would be f R ≈ 1.0036sA, with s being an arbitrary value describing the peak width.Regardless of the exact values for A and s, such a scenario already hints towards a polar response of the medium, but because A is a value that is typically specified in counts of the detection unit (thus well above unity), f R 1 as well.In terms of detection units, absolute values of f R are largely dependent on the evaluation of the measured signal.As shown in the previous paragraph, a sample may be falsely classified if the peak-counts are poorly chosen.Taking Gaussian spectral distributions into account with peak widths FWHM 2ω = √ 2s and FWHM 3ω = √ 3s (FWHM: full width at half maximum), Equation ( 5) becomes Clearly, even when the linear ratio between the peak intensities A 2ω /A 3ω remains constant (i.e., (N A 2ω )/(N A 3ω ) = const.),f R will not, and instead scale with N, thereby being dependent on the absolute number of counts or an equivalent set of data.For clarity, this issue is displayed on the left of Figure 1.Likewise, variation of the fundamental pulse's spectral width s caused by nonlinear processes such as self-or cross-phase modulation, four-wave mixing, Raman shifting, etc., affects f R in a similar fashion, as shown on the right side of Figure 1.
Due to the different slopes between equal linear ratios and f R , one might assume arising ambiguities.Mathematically, this possibility might be valid; however, from a physical point of view, this is rather unlikely.Typical spectrometers often possess a dynamic range of 2 16 values, as opposed to the approximately 2 67 values displayed in Figure 1 (left), thereby greatly reducing the likelihood of an intersection.Secondly, polar and nonpolar samples exhibit vastly different spectral characteristics with similarly large A 2ω /A 3ω -ratios, thereby reducing the chance of a false classification even further.6) with s = 1 (left) and as a function of the spectral peak width and the linear ratio A 2ω /A 3ω with max(A 2ω , A 3ω ) = 1 (right).The white dashed line marks value-pairs where f R = 1, whereas black lines illustrate equal ratios between second harmonic generation (SHG)-and third harmonic generation (THG)-peaks.Values scale from blue over red to yellow.
The third parameter that can vary in Equation ( 6) is the peak width s.Clearly, with rising peak widths, f R also increases, which makes exact comparisons with two different excitation sources difficult.Under normal circumstances, though, the experimental setup is unlikely to notably change, especially with regard to the pump laser.For a single system, s is therefore not expected to fluctuate significantly, thereby greatly reducing the parameter space depicted in Figure 1 (right).The cusp at (1, 1) results from the different exponents of the two peak values A 2ω and A 3ω , where the larger of the two is set to unity.This is why scaling a spectrum according to the larger peak of the nonlinear signal returns adverse results, especially in a fundamental intensity-dependent scan.Again, at ordinary experimental conditions, one peak typically outweighs the other throughout a large regime, except in certain specific scenarios [28] where the low-intensity-limit should be preferred.

Sample Preparation, Linear Optical Analysis, and Nonlinear Setup
Sample preparation follows standard work flows known from linear optical diffuse reflectance spectroscopy; i.e., tight compression of the powder [29].The nanopowder is placed in a steel sample plate (see Figure 2 (left)) and pressed with a planar tappet that is being struck with a hammer.Afterwards, additional material is added until a solid pellet can be obtained.In order to minimize specular reflexes, it is recommended to place a dried sheet of paper between tappet and powder during compression [30].For characterization, two crystallographically different nanopowders have been chosen; namely, LiNbO 3 :Yb and TiO 2 (anatase phase).The former belongs to the PSG 3 m, and is strongly polar [31][32][33], whereas the latter belongs to 4/mm [33,34] and is therefore inversion symmetric.The main part of both samples' size distributions is located at d = (300 ± 100) nm; this way, size-related effects may be neglected.
In order to define the spectral validity for the proposed method, a prior determination of the sample's band gap in a steady-state optical remission setup (FLS980, Edinburgh Instruments, Livingston, UK) seems expedient.As a white standard, the aluminium oxide ZWSK F2000 (Imerys Fused Minerals, Villach, Austria) has been chosen, so that the diffuse reflectance R ∞ is given by the ratio of the sample's against the white standard's remission.The absorption spectrum is subsequently calculated from the Kubelka-Munk relation F(R ∞ ) [35,36].Here, we assume an indirect allowed electronic transition, thereby defining the absorption (in arbitrary units) as α = (F(R ∞ )hν) 1/2 with respect to the incident photon energy hν [37].Lastly, the band gap is extracted from the corresponding Tauc plot [38,39].Figure 3   The naturally white appearance of the samples (e.g., LiNbO 3 :Yb is shown in Figure 2 (left)) is underlined by the measurements and offers a similar-if not better-reflectance to the white standard throughout the visible spectral range.However, a harsh drop-off of R ∞ towards the UV-region can be observed that is attributed to the absorption of radiation.In this case, the band gap amounts to E g = (3.85± 0.01) eV or λ g = (322 ± 1) nm.However, a safe estimate for an upper bound would be λ g ≈ 345 nm for the LiNbO 3 :Yb-nanopowder so that absorption becomes irrelevant.Similar measurements with the TiO 2 -nanopowder as shown in Figure 3 (right) lead to an energy gap of E g = (3.19± 0.01) eV or λ g = (388 ± 1) nm, and consequently to an upper bound of λ g ≈ 415 nm.For THG-measurements, these values should be used as a lower wavelength-limit where f R can still be reliably determined.
A specific experimental setup for nonlinear optical characterization of the sample is presented in Figure 2 (right).A regeneratively amplified and spectrally tunable femtosecond laser system (Astrella HE+ and TOPAS Prime, Coherent) is used to achieve high peak intensities necessary for harmonic generation.The excitation wavelength can be tuned all over the regime of λ ω = (800-2000) nm with a pulse length of τ 40 fs.The fundamental beam is loosely focused onto the sample at a low angle, with the focus point being behind the sample in order to suppress third harmonic generation from air [40].Intensity modulation can be achieved by placing an additional polarizing filter into the beam.Diffusely backscattered harmonic radiation is then collected perpendicular to the surface while avoiding the gathering of any specularly reflected light.The collected signal is spectrally filtered and detected-in our case, by a combination of a spectrograph and a back-illuminated air-cooled CCD-camera (IsoPlane SCT320 and PIXIS 2k, Princeton Instruments, Trenton, NJ, USA).

Harmonic Results
Initial proofs of principle are carried out by successive irradiation of the nanopowders with varying fundamental wavelengths for visual inspection.Results are presented in Figure 4, with the excitation wavelength shown at the bottom of each photograph.Evidently, the two samples show similar spectral characteristics at a first glance with the TiO 2 -nanopowder sample exhibiting a color shift of emitted radiation in comparison to LiNbO 3 :Yb-nanopowder at equal excitation wavelengths.In fact, coming from shorter fundamental wavelengths, TiO 2 -nanopowder shows no visible radiation to the naked eye up to λ ω ≈ 1200 nm when blue light gradually emerges and with further increase of λ ω traverses the entire visible spectrum from blue to red.On the other hand, LiNbO 3 :Yb-nanopowder displays such characteristics twice; one in the same wavelength region, and another at significantly lower wavelengths (as low as λ ω ≈ 800 nm).The transition-zone between these two is located at λ ω ≈ 1400 nm, where the red emission gently gives way to blue light.For further examination, emission spectra have been collected for both samples and are shown in Figure 5 alongside theoretical considerations in terms of nonlinear harmonic emission characteristics.
For analysis, each recorded spectrum in Figure 5 has been baseline-corrected: spectra have been shifted on the intensity-scale such that the intensity is zero at wavelengths corresponding to no harmonic emission.For low fundamental wavelengths λ ω < 900 nm, the influence of the backscattered fundamental radiation can be seen.The additional lobe in the UV region excited at fundamental wavelengths around 900 nm stems from internal reflections inside the spectrograph, and does not correspond to an actual signal.A comparison with the simulated spectra shows that light emission in the visible range from the investigated nanopowders stems from nonlinear harmonic generation that-depending on the fundamental wavelength-can be attributed to either second or third harmonic generation.
At an excitation wavelength of λ ω = 1400 nm, second and third harmonic radiation are inside the band gap of both materials so that subsequent absorption can be neglected.Undoubtedly, LiNbO 3 :Yb-nanopowder shows a significant amount of THG-intensity, so further investigation has been carried out at this wavelength.In Figure 6, harmonic emission peaks of both samples at this wavelength are shown in more detail.With suitable magnification, the weaker harmonic emission can be detected from both samples, with TiO 2 offering SHG-intensity several orders of magnitude lower than THG while being located at the flank of the THG emission peak.For calculation of the integrated spectral intensity, this behavior has been accounted for by firstly subtracting a linear term from the SHG-peak modeling the THG-flank, thus rendering the influence of the latter to zero.With the more intense peak of each spectrum scaled to unity (as shown in Figure 6), data can be extracted for calculation of the harmonic ratio f R which are presented in Table 1.
Table 1.Relevant values for the determination of the harmonic ratio f R extracted from Figure 6.FWHM is the full width at half maximum of each peak, the signal S is given by the area under the peak and the indexes 2ω and 3ω denote the second and third harmonic, respectively.In case of LiNbO 3 :Yb-nanopowder, the harmonic ratio f R is well above unity, thus supporting the assumption of a polar medium.On the other hand, the TiO 2 -nanopowder exhibits only minor SHG emission, and THG dominates the harmonic spectrum by several orders of magnitude, thus resulting in f R 1.This leads to the conclusion that measurements confirm the powder's inversion symmetry.

Intensity (In-)Dependence of the Harmonic Ratio
Due to the choice of exponents in Equation ( 5), the harmonic ratio is expected to be unaffected by intensity changes of the fundamental wave I ω .This functional (in-)dependence of f R on I ω is depicted in Figure 7 for both LiNbO 3 :Yb-and TiO 2 -nanopowder at a fundamental wavelength of λ ω = 1400 nm and λ ω = 1460 nm, respectively.The latter has been chosen for TiO 2 as a result of higher harmonic emission intensity, which is especially necessary for detection of the SHG-signal at lower fundamental intensities.In contrast to results presented in Figure 6, data has not been scaled to unity, thereby altering the magnitude of f R .However, the overall result remains largely the same: f R is well above unity for the polar medium and several orders of magnitude below unity for the nonpolar one.Furthermore, the value remains constant over almost an order of magnitude of the fundamental intensity.Deviations from the uniformity can be observed in the low-intensity-regime, where the signal to noise-ratio drastically diminishes and thereby lets the weaker peak fade into a quasi-constant noise-background.This behavior is indicated in Figure 7 by the dashed guidelines.It should be noted that although the absolute values of f R have increased compared to the ones from Table 1, they are still separated from each other by over 13 decades; in fact, a separation factor on the order of 10 15 is achieved.

Discussion
The results give strong evidence that our tool is appropriate to unambiguously distinguish between polar and nonpolar harmonic upconversion nanoparticles as exemplarily demonstrated for LiNbO 3 :Yb and TiO 2 .We note, however, that the absolute values of the harmonic ratio f R largely depend on the harmonic signals' strengths, and thus on the specific experimental setup-in particular on the detection unit.Therefore, either spectroscopic calibration of the detection unit or the use of the same nonlinear optical setup throughout a systematic study as shown in this article is recommended.We further note the importance of the specific units involved in the calculation of f R , especially in the context of comparability with other setups and materials.For this purpose, f R is to be determined from the spectral intensity detected by the photometer (in counts) as a function of the photon energy (in eV).
Additionally, with regard to f R , preliminary considerations in Section 2.2 have outlined the limits of its applicability.In fact, actual experimental conditions narrow the regime of physically realistic values even further, as shown in Figure 8.At a maximum peak height of 50,000 counts, a typical value for a detection system scaling signals in 2 16 steps, and with a linear ratio given by the scaling factors of Figure 6, LiNbO 3 :Yb-and TiO 2 -nanopowders can be easily distinguished.Furthermore, they are both well-displaced in relation to the f R = 1-line by approximately six to seven orders of magnitude.Taking the different scaling behaviors of the second and third harmonic intensities into account, weaker signals will run parallel to the f R = 1-line as intensity-independent constants, and thus-ideally-never cross it.Of course, with diminishing signals, noise will eventually overwhelm the signal and may cause inconstant harmonic ratios that might approach unity as visible for both samples in Figure 7.With the goal of measuring high-quality spectra, the detection system should be readjusted at such points.
The spectral peak width-dependency as illustrated in Figure 8 (right) once more clarifies the clear distinction of f R between the two material classes.The linear ratios are taken from Figure 6, and the larger peak is scaled to 50,000 counts.In this case, the harmonic ratios differ by roughly thirteen decades, and each remains either above or below unity with no crossings in the depicted regime.When the FWHM is measured on the energy scale, typical values with our setup will be on the order of (10-100) meV, which corresponds to the shaded area.In fact, a Fourier-limited 40 fs-pulse possesses a FWHM of 46 meV, which is marked by the vertical black dashed line [41].In order to attain a harmonic ratio of unity for LiNbO 3 :Yb-nanopowder with the same peak height yields a minimum pulse length τ min ≥ 450 ns.Therefore, with femtosecond-pulse setups, mix-ups with regard to the polarity-even with no reference standards for a given experimental arrangement-are highly unlikely.
Lastly, mixtures of particles with differing PSGs may be subject to investigation.In this case, even small fractions of polar samples will strengthen the otherwise minuscule SHG-peak, and thereby quickly push f R towards and above unity.In this case, the proposed method and evaluation can be used to trace fractions of polar media in an ensemble of nonpolar particles.For absolute determination of fractions, preliminary calibration with reference standards is required-ideally with pure samples whose mixtures are to be analyzed.
In fact, through calibration against a known material, this setup may also be capable of obtaining absolute values of nonlinear coefficients.The standard would be a nanocrystalline material of a well-defined size whose nonlinear coefficients-ideally both second-and third-order as a function of the wavelength-are known.Similar to linear optical characterization, measurements would then investigate upon the backscattered harmonic spectra, and thereby be able to examine material, size, and wavelength dependencies successively with a growing pool of data.
For LiNbO 3 :Yb-nanopowder, additional wavelength-dependent characterization has been conducted in the range λ ω = (1140-1580) nm.Within this range, slight fluctuations of the harmonic emission itself have been observed that could stem from the underlying absorption profile and especially the band gap itself [6].On the other hand, the harmonic ratio f R remains well above 10 5 at fundamental intensities I ω = (1-9) PW/m 2 .Thus, the applicability of arbitrary fundamental wavelengths can be demonstrated, which is ultimately only limited by the optical absorption characteristics of the material under investigation.Detrimental factors can affect either THG (typically at the band edge), SHG (color centers and other defects), or both.In such case, an appropriate wavelength-triplet between fundamental, second harmonic, and third harmonic wave should be determined prior to nonlinear investigation; i.e., by diffuse reflectance spectroscopy, outlined in Section 3.
In conclusion, the presented investigations have revealed that nonlinear diffuse fs-pulse reflectometry serves as a tool for quickly assessing (non-)linear optical properties of nanomaterials.Together with the introduction of the harmonic ratio f R as a figure of merit, structural analysis has been conducted on a polar and a nonpolar substance with clear results.By using femtosecond pulses, the risk of thermal degradation can be vastly reduced.A positive side-effect is the boost of nonlinear harmonic generation by the inherently large peak intensities of such laser sources.Yet, the most striking advantage in comparison to other approaches is that this method is specifically tailored to the analysis of powders, thereby making use of the sheer number of nanocrystals to augment upconverted signals to easily detectable levels.Since the particles will be isotropically aligned in such powder plaques, a nonlinear signal will be picked up due to orientational averaging as soon as a single tensor-element is nonzero.
To summarize the workflow for structural clarification through investigation of the polarity: (1) synthesis of the material to be analyzed; (2) compressing the powder into a pellet; (3) obtaining a remission-spectrum to find an appropriate wavelength triplet of fundamental, SHG, and THG with negligible absorption; (4) recording of harmonic spectra; (5) determination of harmonic signals S nω by integration of the respective baseline-corrected peak on an energy scale in eV; (6) calculation of the harmonic ratio f R according to Equation ( 5); (7) if f R 1, the material is likely polar, if f R 1, the powder exhibits nonpolar characteristics.
Furthermore, this technique is not limited to purely structural analysis.Additional cogitable subjects include higher harmonic generation beyond frequency-tripling, analysis of mixtures of two or more substances, or specific tailoring of illuminants for industrial purposes.It may further serve as an important additional tool for the determination of nonlinear optical coefficients when combined with polarization-dependent microscopy [42].Extension to a vast array of nanoscaled materials, each with their very own properties, can thus be easily realized and will therefore likely be subject to future investigations.

Figure 1 .
Figure 1.Logarithmized harmonic ratio log 10 ( f R ) as a function of spectral peak intensities according to Equation (6) with s = 1 (left) and as a function of the spectral peak width and the linear ratio A 2ω /A 3ω with max(A 2ω , A 3ω ) = 1 (right).The white dashed line marks value-pairs where f R = 1, whereas black lines illustrate equal ratios between second harmonic generation (SHG)-and third harmonic generation (THG)-peaks.Values scale from blue over red to yellow.

Figure 4 .
Figure 4. Harmonic emission of LiNbO 3 :Yb-(top) and TiO 2 -nanopowder (bottom) at different fundamental wavelengths in the range λ ω = (900-1900) nm photographed with a digital camera.A shortpass filter was placed in front of the camera lens to suppress residue infrared pump radiation.

2 Figure 5 .Figure 6 .
Figure 5. (Left): baseline-corrected emission spectra of LiNbO 3 :Yb-(top) and TiO 2 -nanopowder (bottom) as a function of the fundamental wavelength.Dashed lines indicate SHG-and THG-peak positions.Each spectrum's most significant peak has been scaled to unity; (Right): simulated spectra with their respective emission characteristics based on the spectral properties of the fundamental beam.All spectra are plotted on a logarithmic intensity scale with values exceeding unity clipped to one.Values scale from blue over red to yellow.

Figure 8 .
Figure 8. Left: logarithmized harmonic ratio log 10 ( f R ) as a function of spectral peak intensities according to Equation (6), with s = 46 meV; the colors scale from −11 (blue) to 13 (yellow).Harmonic ratios of LiNbO 3 :Yb (dotted) and TiO 2 (dash-dotted) are simulated for max(A 2ω , A 3ω ) = 50,000, their linear ratio given by the scaling factors of Figure 6, and extrapolated to low intensities.Right: f R of LiNbO 3 :Yb and TiO 2 as a function of the fundamental peak width s measured on an energy-scale with the linear ratio taken from Figure 6 and max(A 2ω , A 3ω ) = 50,000.The horizontal dashed black line marks f R = 1, whereas the vertical one marks s = 46 meV.Shaded in light gray are typically accessible experimental peak widths.Both ratios are separated by a factor of almost 10 13 .