Organic Crystals for THz Photonics

: Organic crystals with second-order optical nonlinearity feature very high and ultra-fast optical nonlinearities and are therefore attractive for various photonics applications. During the last decade, they have been found particularly attractive for terahertz (THz) photonics. This is mainly due to the very intense and ultra-broadband THz-wave generation possible with these crystals. We review recent progress and challenges in the development of organic crystalline materials for THz-wave generation and detection applications. We discuss their structure, intrinsic properties, and advantages compared to inorganic alternatives. The characteristic properties of the most widely employed organic crystals at present, such as DAST, DSTMS, OH1, HMQ-TMS, and BNA are analyzed and compared. We summarize the most important principles for THz-wave generation and detection, as well as organic THz-system conﬁgurations based on either difference-frequency generation or optical rectiﬁcation. In addition, we give state-of-the-art examples of very intense and ultra-broadband THz systems that rely on organic crystals. Finally, we present some recent breakthrough demonstrations in nonlinear THz photonics enabled by very intense organic crystalline THz sources, as well as examples of THz spectroscopy and THz imaging using organic crystals as THz sources for various scientiﬁc and technological applications.


Introduction
Terahertz (THz) waves have recently gained a considerable interest due to a wide variety of interesting applications ranging from fundamental scientific investigations to various industrial implementations. For example, THz waves are used to study fundamental carrier dynamics phenomena and phonon/vibrational modes of various materials systems [1] and even of single molecules [2]. On the other hand, there is a wide interest in THz imaging and spectroscopy for industrial materials testing, security, medical applications, food, construction, and the pharmaceutical industry [1]. This development has been made possible by the wide variety of THz sources that have become available during the last few decades, such as quantum cascade lasers, Gunn lasers, semiconductor antennas, as well as sources that are based on nonlinear optical (NLO) effects. While many of these sources cover only a limited THz range or can emit at particular THz frequencies only, the NLO-based sources are of particular interest because they can cover the complete THz frequency band (frequencies from 0.1 THz to beyond 10 THz). The NLO-based sources are either based on the NLO process of difference frequency generation or optical rectification. The NLO materials' requirements for both of these types are similar; they need to have a high second-order NLO susceptibility and should possibly satisfy the phase-matching condition between the fundamental optical pump waves and the generated THz waves. For this purpose, mainly three different material systems are considered, semiconducting NLO crystals such as zinc telluride (ZnTe) or gallium phosphide (GaP), the conventional inorganic NLO crystals such as lithium niobate (LiNbO 3 ) and organic NLO materials.
The purpose of this review is to discuss the progress in the development of THz sources and applications based on organic NLO crystals. Section 2 presents the origins of NLO response in organic materials, the intrinsic differences between organic and inorganic NLO materials, and the relations between the relevant microscopic and macroscopic properties of organic NLO crystals. We present and compare the linear and the nonlinear optical properties of some of the presently most widely used organic crystals for THz photonics, including DAST, DSTMS, OH1, HMQ-TMS, and BNA. Section 3 presents theoretical fundamentals of THz-wave generation via difference-frequency generation and via optical rectification. We summarize THz incoherent and coherent detection techniques that are most widely employed in combination with organic THz sources. The important phase-matching requirements for efficient THz-wave generation and detection based on NLO principles are also discussed. We compare THz refractive indices and THz absorption properties of most widely used organic crystals in relation to their phase-matching characteristics. Section 4 presents some of the most successful examples of organic-based THz sources in terms of power or bandwidth, both the narrowband tunable organic THz sources based on difference-frequency generation and broadband THz sources based on optical rectification. Section 5 summarizes various photonics application demonstrations based on organic crystalline THz sources. These are for example important fundamental demonstrations in the exciting new field of nonlinear THz photonics, as well as various non-destructive THz spectroscopy and THz imaging approaches, interesting for a variety of future scientific and industrial applications.

Second-Order Organic Nonlinear Optical (NLO) Crystals
The intrinsic NLO response of an organic material depends on the core molecular structure and microscopic nonlinearity of the constituting NLO molecules, their orientational distribution in a particular macroscopic system and the interactions between the molecules. For the particular NLO system configuration such as e.g., for THz-wave generation, the orientation of the molecules relative to the interacting electromagnetic waves and the macroscopic dimensions are essential. In this section, we discuss the involved physical processes and the state-of-the-art material parameters achievable in NLO organic crystals.

NLO Effects in Organic Materials
NLO organic materials are based on the nonlinearity of the constituting molecules. An electric field E induces a charge redistribution in an organic molecule. This induces a dipole moment p (in addition to the permanent dipole moment µ µ µ of a molecule), which is most commonly described by the power series as a function of the electric field E p i = 0 α ij E j + 0 β ijk E j E k + 0 γ ijkl E j E k E l + . . . , (1) where the usual Einstein summation convention holds and where 0 is the electric constant, α ij the linear polarizability, β ijk the first hyperpolarizability and γ ijkl the second hyperpolarizability of the molecule. Note that, in this review, we use, unless explicitly stated, SI units and the above defined power series expansion, which is the same as in Ref. [3]. Note that several different formulations are used for the above equation in both SI and cgs unit systems, so care must be taken when comparing the results from different sources [3][4][5]. While all molecules exhibit α ij and γ ijkl , the first hyperpolarizability β ijk is responsible for the asymmetric response to an external field and is only possible in asymmetric molecules. For molecules with non-zero β ijk , the second-order NLO effects such as second-harmonic generation, electro-optic effect, optical rectification, and THz-wave generation are possible and are of interest in this review. The conventional design of organic molecules to achieve a high first hyperpolarizability β ijk is based on a π-conjugated bridge with strong donor and acceptor substituents at the end of the molecule to promote the asymmetric response and has been extensively studied [6][7][8][9][10][11][12][13]. Figure 1a illustrates the simple physical mechanism responsible for the NLO response of an asymmetric molecule (β ijk = 0) to an external electric field: the electronic cloud favors the acceptor A over the donor D. As a consequence, the amplitude of the induced dipole moment depends on the polarity of the applied electric field. Note that the NLO response is usually significant only for very strong electric fields E. π-conjugated nonlinear optical (NLO) molecule (chromophore) upon application of an external sinusoidal electric field E(t); (b) NLO molecules in the simple oriented-gas model. The molecular-ordering angle θ kz is the angle between the charge-transfer axis z of the molecule and the polar axis k of the crystal; (c) an example of a centrosymmetric crystal structure with θ kz = 90 • , resulting in zero macroscopic second-order NLO susceptibility χ (2) . Examples of non-centrosymmetric structures with (d) perfectly parallel chromophores with θ kz = 0 • , maximizing the diagonal component of the χ (2) tensor and with (e) zig-zag-chromophore orientation with 0 < θ kz < 90 • with both diagonal and non-diagonal components of χ (2) ; (f) the acentric order parameter for the diagonal (cos 3 θ kz ) and non-diagonal (cos θ kz sin 2 θ kz ) components of χ (2) as a function of θ kz for simplest 1D-chromophores with one dominant first hyperpolarizability component β zzz .
general we need to consider the full tensorial nature of the first hyperpolarizability β ijk . The first hyperpolarizability of NLO chromophores can be evaluated theoretically using quantum chemical calculations [6,15] or measured experimentally, mostly using the Electric Field Induced Second Harmonic Generation (EFISH) [16] or the hyper-Rayleigh scattering (HRS) technique [17].
The macroscopic NLO material response is in general (also for inorganic materials) described using an analogous power series for the induced material polarization P as a function of the applied electric field E P i = 0 χ (1) ij E j + 0 χ (2) ijk E j E k + 0 χ (3) ijkl E j E k E l + . . . , (2) where χ (1) is the linear susceptibility, while χ (2) and χ (3) are the nonlinear susceptibilities of the second and the third order, respectively. The susceptibility tensor coefficients are all dependent on the frequency f (or the angular frequency ω = 2π f ) of the involved electric fields E.
As seen in Equation (2), the NLO effects of the second-order, which are of interest in this review, allow interaction of two electric fields E(ω 1 ) and E(ω 2 ) to induce a macroscopic polarization at a combination of their frequencies ω 3 , which is the source of the generated electric field E(ω 3 ). Figure 2 shows schematic illustrations of two basic second-order NLO processes: sum-frequency generation and difference frequency generation. In the sum-frequency generation process, the generated frequency is the sum of the two input frequencies ω 3 = ω 1 + ω 2 . From the quantum optics point of view, an NLO material allows interactions between photons, so that the total photon energy is conservedhω 3 =hω 1 +hω 2 and the total photon momentum is conservedhk 3 =hk 1 +hk 2 (h is the reduced Planck constant and k i the wave vector of the corresponding optical wave). The momentum conservation is in the field of NLO commonly known as the phase-matching condition. In Figure 2c, a general non-collinear phase-matching geometry is presented, while practically often a collinear geometry is preferred, in which all wave vectors are parallel. Any phase mismatch decreases the efficiency of the NLO conversion process. Phase matching requirements relevant for THz-wave generation are discussed in Section 3.4, while the reader is referred to general books on NLO for phase-matching approaches for other processes (see e.g., Refs. [18][19][20]). The schematic energy-level diagram in Figure 2b illustrates the photon energy conservation with the electronic ground state indicated by the solid line and the virtual states indicated by the dashed lines. The corresponding χ (2) NLO susceptibility depends on the frequencies of the involved fields, depending on how close these virtual levels are to the real electronic levels. In case one of the virtual levels is close to a real level, the NLO susceptibility is said to be resonantly enhanced. The frequency dependence of the NLO susceptibility for sum-frequency generation is here denoted as χ (2) (−ω 3 ; ω 1 , ω 2 ), where ω 3 is the frequency of the generated electric field. This description holds for the so-called parametric NLO processes, for which the photon energy is conserved and there is no electronic transition within the material involved. For non-parametric processes such as, e.g., the two photon-absorption, the imaginary part of the NLO susceptibility should also be considered [19]. A special case of sum-frequency generation is the well-known second-harmonic generation process with ω 1 = ω 2 = ω, where the generated frequency is twice the pump optical frequency ω 3 = 2ω and only a single-frequency pump optical field is required.

Sum-frequency generation
The difference-frequency generation process, also called an optical parametric amplification process, in which generated frequency ω 3 presents the difference of the pump frequencies ω 1 − ω 2 , can be explained by the absorption of the first photon at ω 1 to a virtual state and then a two-photon emission process, stimulated by the presence of the field at ω 2 (see Figure 2b). The field at ω 2 can be also spontaneously generated; in this case, the process is referred to as the optical parametric generation. An appropriate optical resonator can be added to an optical parametric generator to enhance the emission; this case is referred to as the optical parametric oscillator.
A special case of difference-frequency generation is the optical rectification process, for which ω 1 ≈ ω 2 ≈ ω and the generated field is a (quasi-)static field at zero (or close to zero) frequency ω 3 ≈ 0. THz-wave generation by difference-frequency generation or optical rectification results in a field at the THz frequency ω 3 ω 1 , ω 2 as discussed in more detail in Section 3. For THz-wave generation and electro-optic detection, the relevant processes involve static or quasi-static fields (including THz fields), i.e., oscillating at frequencies much lower than the optical frequencies. The electro-optic effect is most often described as a field-induced change of the refractive index n at ω as where the usual summation convention holds and where r ijk is the electro-optic tensor and can be related to the second-order susceptibility as where n i (ω) are the frequency-dependent refractive indices along the corresponding dielectric axis i. The linear susceptibility tensor χ (1) and the related relative permittivity tensor = χ (1) + 1 describe the common linear polarization response, which depends on the frequency of the applied field E(ω). Figure 3 schematically shows the frequency dependence of the relative permittivity (ω). The relative permittivity at low frequencies T (zero-stress relative permittivity) contains different contributions from the material polarization response: from acoustic phonon vibrations a , from optical phonon vibrations o , which also includes molecular vibrations in organic materials, and from electronic vibrations e T = a + o + e = a + S , where S = o + e is the clamped (zero-strain) relative permittivity and e determines the refractive index in the optical range as e = n 2 . Analogously, the nonlinear susceptibilities χ (2) and χ (3) also contain contributions from different possible material polarization responses [6,21,22]. The electrons can follow the fastest electric-field oscillations, therefore the electronic polarizability contributes to the linear and the nonlinear susceptibility up until the band-to-band transitions (electronic excitations). The optical phonon vibrations and the molecular vibrations contribute to the polarization response below about the THz/far-infrared range. The acoustic lattice vibrations can only follow relatively slow electric field oscillations below about the MHz range. Due to the different origin of the polarizability in organic and inorganic materials, the frequency dispersion of χ (i) (i = 1, 2 or 3) is essentially different, which is important for applications. For example, for the most well-known organic crystal DAST (4-N, N-dimethylamino-4 -N -methyl-stilbazolium tosylate), the optical refractive index along the polar axis a is about 2.1 and in the THz range it is about 2.3. The square root

√
T is also about 2.3 [22,23], which points to the very low dispersion from the low to the optical frequencies. For the most well-known inorganic NLO crystal LiNbO 3 , the optical refractive index along the polar c axis is about 2.2 and in the THz range it is about 5, which points to the much larger contribution of lattice vibrations to the polarizability of this crystal. This is also observed in the dispersion of the NLO susceptibility χ (2) : while the electro-optic coefficient (involving quasi-static fields) is of the same order of magnitude for DAST (r 111 = 48 pm/V at 1535 nm) and LiNbO 3 (r 333 = 30 pm/V at 1064 nm), the difference is much larger for the NLO coefficient for second-harmonic generation, involving only the optical fields (d 111 = 290 pm/V at 1535 nm for DAST and d 333 = 34 pm/V at 1064 nm for LiNbO 3 ) [21,22,24]. This shows that the origin of the linear and the nonlinear optical susceptibility in organic materials is essentially electronic and therefore extremely fast (femtoseconds). This is an important advantage for high-speed electro-optic modulation [6,25,26] and also for THz-wave generation and detection: due to the approximate matching of the THz and the optical refractive index, it is much easier to achieve the phase matching needed for highly efficient response, as discussed in Section 3.4.   . Schematics of the material linear polarization response as a function of frequency of the external electric field. For organic materials, the contributions of o and a to the static relative permittivity T are much less important compared to inorganic materials, which is an important advantage for high-speed electro-optic modulation and THz-photonics applications (see text). The scale below shows the THz frequency range in frequency f , wavelength λ, and spectroscopic wavenumber 1/λ.

Organic Single Crystalline Materials
To achieve a macroscopic second-order NLO response in organic materials with non-zero χ (2) , the constituting asymmetric NLO chromophores have to be distributed non-centrosymmetrically. This is by far not trivial because these molecules usually possess a high permanent dipole moment and they tend to orient antiparallel, as illustrated in Figure 1c, which cancels the macroscopic χ (2) response. There have been many possibilities proposed to achieve a macroscopic non-centrosymmetric order of NLO molecules. Most commonly, a non-centrosymmetric order is achieved by incorporating NLO chromophores in a polymer matrix and then pole the material by applying a very strong static electric field (just below the electric breakdown) above the glass transition temperature to align these molecules, based on their relatively high permanent dipole moment µ. The advantage of this approach is that it is in principle applicable for all chromophores with a permanent dipole moment, including those with the highest first hyperpolarizabilities β [6]. On the other hand, it is extremely challenging to achieve a desired long-term orientational stability of the non-centrosymmetric order at room temperature [6]. For THz photonics applications, the achievable material thickness of polymer-based materials is also very limited [27]. For highly-efficient THz-wave generation, the thickness of the THz generation material (organic or inorganic) in the range of 0.1-1 mm is required [28]. Another possibility to obtain an efficient macroscopic second-order active NLO organic material is by crystallization, offering the highest possible chromophore density, the best possible long-term orientational stability and sufficient sample thickness (∼ 0.1-1 mm range), which is the approach discussed in this review. Additional to the much superior temporal stability of crystals compared to poled polymers, their photochemical stability may be better by several orders of magnitude [29], which is particularly important for THz-wave applications, where high optical power densities are used and therefore a high damage threshold is required. On the other hand, only a very small ratio of known chromophores can crystallize in a non-centrosymmetric crystalline space group and there seems to be a limit in the conjugation length (i.e., the first hyperpolarizability strength β) that can still lead to a non-centrosymmetric alignment [15,30]. The chromophores, besides having a large β, need to be usually specially designed to promote a non-centrosymmetric alignment [15,30]. In addition, even for the chromophores that naturally self assemble into a non-centrosymmetric structure, the growth of high optical quality crystals large enough and with desired morphology may be too challenging to allow for practical applications. For all these reasons, although there are endless possibilities for molecular engineering to obtain organic NLO crystals, there have been only a few identified that lead to crystals that can be potentially used for practical NLO applications and some of these are presently commercially available. The commercial availability of these crystals has helped to considerably expand the research possibilities particularly for THz-wave generation with organic NLO crystals during the last couple of years, which has additionally motivated the continued research and development of novel and optimized organic NLO crystals for THz photonics applications.
The relation between the macroscopic and the molecular NLO coefficients is complex because the interactions between the neighboring molecules influence the first hyperpolarizability β in the crystalline form [11,31,32]. However, often the macroscopic second-order nonlinearities of organic materials can be reasonably well explained by the nonlinearities of the constituent molecules by the so-called oriented-gas model [14]. This relation is in particular simple for the chromophores in the 1D approximation, for which only one dominant tensor element β = β zzz can be considered, which is the case of most state-of-the-art organic NLO crystals. For example, the maximal diagonal element for the NLO susceptibility along the polar axis k of the crystal can be expressed as χ (2) kkk = N f local β eff kkk = N f local cos 3 θ kz β zzz .
The above relation reveals the most important contributions to the NLO effect in organic materials: the chromophore number density N, the local-field factor f local , the molecular nonlinearity β zzz along the main charge-transfer axis of the chromophore, and the orientational factor, the so-called acentric order parameter cos 3 θ kz , which is averaged over all molecules in a macroscopic system. For crystalline materials, the angle θ kz is the angle between the main charge transfer axis z of the chromophores and the polar axis k of the crystal, as illustrated in Figure 1b. Note that, for simplicity and practical relevance, we define here the polar axis k along the effective charge-transfer direction of the crystal, which may be, in low symmetry crystals (monoclinic and triclinic) slightly off (for highly nonlinear crystals, not more than a few degrees), the real polar axis defined by the effective permanent dipole moment of the crystal.
The effective first hyperpolarizability coefficients β eff ijk in the crystal can be introduced to consider the first hyperpolarizability tensor components β xyz of all single chromophores in the crystalline unit cell projected to a suitable Cartesian coordinate system of the crystal based on the oriented-gas model as [14,15] where n denotes the number of chromophores in the unit cell, s the site in the unit cell and θ s im are the angles between the dielectric axes i (i = 1, 2, 3) of the crystal and the molecular axis m (m = x, y, z) of the chromophore. The macroscopic second-order susceptibility tensor components can be then evaluated as Note that both the effective first hyperpolarizability components β eff ijk as well as the local-field factors f local strongly depend on the frequencies of the interacting electric fields.
The effective first hyperpolarizability components vanish in case of a centrosymmetric crystalline arrangement (β eff ijk = 0) like the alignment illustrated in Figure 1c. For perfectly parallel chromophore alignment along the k-axis as illustrated in Figure 1d and 1D chromophore approximation, the diagonal component along the polar axis is maximal β eff kkk = β zzz and all other tensor components of β eff ijk are zero. For an orientation with 0 < θ kz < 90 • and in 1D chromophore orientation as illustrated in Figure 1e, the maximal diagonal effective first hyperpolarizability component is β eff kkk = β zzz cos 3 θ kz and the maximal off-diagonal effective first hyperpolarizability component is β eff iik = β zzz cos θ kz sin 2 θ kz . The dependence of maximal diagonal and maximal off-diagonal components of β eff ijk are illustrated in Figure 1f. Note that in case the chromophores are additionally distributed in the azimuthal angle around the polar axis, another projection needs to be considered for different off-diagonal effective first hyperpolarizability components [14]. As seen in Figure 1f, to maximize the diagonal NLO susceptibility χ (2) kkk along the polar axis k, the projection angles θ kz should be close to zero, i.e., the charge transfer axes of the molecules should be close to parallel. To maximize the off-diagonal susceptibility elements χ (2) ikk , the molecules should make an angle of about 54.7 • with respect to the polar axis [14]. For electro-optic and THz photonics applications, for the majority of cases, the large diagonal NLO coefficient is desired and therefore the molecular orientation with small or zero θ kz is optimal.
The dispersion of the macroscopic NLO susceptibilities is usually modeled using Equation (8), by assuming a simple two-level model [33] for the dispersion of first hyperpolarizabilities β zzz and local field factors f local [8,9]. The dispersion depends on the frequencies of all involved interacting electric fields, therefore the NLO susceptibilities are different for different NLO processes. For the linear electro-optic effect, the dispersion of the dominant first hyperpolarizability component is in the two-level model given by and for the second-harmonic generation by where ω eg denotes the resonant angular frequency of the electronic transition from the ground (g) to the excited (e) state and a dispersion-free first hyperpolarizability β 0 presents the first hyperpolarizability extrapolated to infinite wavelengths (ω → 0), away from electronic resonances. The simple two-level model describes well the dispersion away from electronic resonances. More sophisticated models are available for modeling the dispersion closer or at electronic resonances such as the homogeneously damped two-level model [34,35], the homogeneously broadened vibronic model [34,36] and the inhomogeneously broadened vibronic model [36]. The dispersion-free first hyperpolarizability β 0 is a very useful parameter when comparing the nonlinearity of molecules with different resonant frequencies ω eg . In the conventional two-level model [33], the dispersion-free first hyperpolarizability β 0 depends on ω eg , the difference ∆µ = µ e − µ g between the excited µ e and the ground-state µ g dipole moments and the transition dipole moment µ eg as The inverse-square dependence of β 0 on ω eg is often referred to as the nonlinearity-transparency trade-off: the chromophores with a red-shifted absorption maximum (lower transmission in the visible wavelength range, lower ω eg ) typically present a higher β 0 . This is strongly related to the best phase-matching range for THz-wave generation with state-of-the-art organic materials, as discussed later on.
The macroscopic second-order susceptibility coefficients χ (2) ijk can be evaluated by the conventional Maker-fringe technique and electro-optic measurements [6,9,37,38]. For these measurements, high optical quality plane-parallel samples with polished surfaces and possibly with different orientations are required and the refractive index dispersion should be also measured with a relatively high precision (see an example in Ref. [39]), which is very time-consuming. For newly developed organic NLO crystals, simpler screening methods are therefore employed to evaluate their potential for NLO applications.
A relatively simple Kurtz and Perry second-harmonic generation powder test can be performed to screen the new materials [40]. The test can be done qualitatively only to determine if the material is non-centrosymmetric and estimate the order of magnitude for the χ (2) -effects, but also more qualitatively including the measurements of the second-harmonic generation efficiency as a function of powder size and as a function of pump wavelength, which may indicate a possible phase-matching enhancement [40]. For highly nonlinear organic crystals that typically present a red-shifted absorption edge in the visible (above ca. 600 nm, see the following section), the coherence length for second-harmonic generation is relatively short (in the order of 1 µm), much smaller than the typical powder size, which gives then a relatively good estimate for the macroscopic susceptibility based on the powder test (see, e.g., Figure 15 in Ref. [15]).
Another possibility to estimate the χ (2) -potential of novel materials is to evaluate first the molecular nonlinearities β xyz by quantum chemical calculations, determine the crystal structure by X-ray diffraction and then use the above-mentioned oriented-gas model according to Equations (7) and (8). Because of the additional influence of intermolecular interactions and local field corrections, the theoretical approximation with the effective first hyperpolarizability tensor β eff ijk only roughly explains the actually measured NLO and electro-optic coefficients but is still well suited to compare different materials and their potential for NLO applications [15].
More recently, several groups have used newly developed organic NLO crystals directly for THz-wave generation experiments and compared the generated THz fields with some standard material such as ZnTe; some examples are mentioned in Section 4.2. This method is also relatively convenient in case an appropriate THz system is available because the basic requirements for the sample quality for THz-wave generation are relatively low compared to the requirements for the optical characterization. This is because the (surface) quality requirements scale with the wavelength, which is about two orders of magnitude larger for THz compared to the optical waves. Figure 4 shows the crystal structures and the dominant sample morphology of some of the most well investigated organic NLO crystals for THz-wave generation during the last decade: the stilbazolium salts DAST and DSTMS (4-N,N-dimethylamino-4 -N -methyl-stilbazolium 2,4,6-trimethylbenzenesulfonate), the non-ionic configurationally locked polyene crystal OH1 (2-(3-(4-hydroxystyryl)-5,5-dimethylcyclohex-2-enylidene)malononitrile) and the quinolinium salt HMQ-TMS (2-(4-hydroxy-3-methoxystyryl)-1-methylquinolinium 2,4,6-trimethylbenzenesulfonate). The relevant crystallographic properties and the microscopic/macroscopic NLO properties of these crystals are listed in Table 1, which includes another widely investigated material BNA (N-Benzyl-2-methyl-4-nitroaniline). Examples of refractive index dispersion and electro-optic coefficient dispersion in the optical range are shown in Figures 5 and 6. In the following subsections, we discuss these properties in more detail for the selected ionic and non-ionic organic NLO crystals.  [32,44], the non-ionic configurationally locked polyene crystal OH1 (c) [45,46] and the quinolinium salt HMQ-TMS (d) [47]. The structures are available at the Cambridge Crystallographic Centre (CCDC) for the corresponding publications. The red arrows indicate the direction of the main charge-transfer axis of the NLO chromophores. The typical plate-like sample morphologies used in THz applications (as-grown for DAST, DSTMS, OH1 and cleaved for HMQ-TMS) have been drawn based on the information given in the corresponding references.

Ionic Organic Crystal DAST
DAST is the most well known and widely investigated organic electro-optic crystal. DAST was first reported already in 1989 by [41,42] and is still being recognized as one of the best state-of-the-art organic NLO crystals. High optical quality and large size DAST crystals can be grown from methanol solution by the slow cooling method [43,48]. DAST crystals have high second-order NLO and electro-optic coefficients, being respectively ten times and twice as large as those of the inorganic standard LiNbO 3 . Due to the relatively low zero-stress relative permittivity compared to inorganic alternatives, as discussed in the previous section, DAST is highly attractive for high-speed electro-optic applications and broadband THz-wave generation [22]. DAST is an organic salt that consists of a positively charged stilbazolium cation and a negatively charged tosylate anion. The stilbazolium cation is one of the most efficient NLO active chromophores that can crystallize in a non-centrosymmetric structure [15,30]. The counter ion tosylate is used to promote the non-centrosymmetric crystallization, but its contribution to the macroscopic nonlinearity is negligible compared to the cation [41,49]. Table 1. Crystallographic symmetry, wavelength of absorption maximum λ max in methanol solution, main first hyperpolarizability component β zzz (non-resonant) calculated by quantum-chemical calculations, the angle θ kz between the main charge-transfer axis z of the chromophores and the polar axis k of the crystal, the non-zero elements of the effective first hyperpolarizability tensor β eff ijk (non-resonant; i.e., Kleinman symmetry valid), and the measured electro-optic coefficients r ijk and nonlinear optical (NLO) coefficients for second-harmonic generation d ijk (ω) = 1 2 χ (2) (−2ω; ω, ω) at different wavelengths of some of the most well investigated organic NLO crystals for THz-wave generation. The results for first hyperpolarizabilities are given in cgs units because these are most often reported; 10 −30 esu corresponds to 4.18 · 10 −40 m 4 /V according to the definition of the first hyperpolarizability in Equation (1). [32,44,50,52] (m) [39,45,46] (mm2) [47,53,54] (m)

Crystal Name Space Group
The structure of DAST is shown in Figure 4a. The chromophores are packed with their main charge transfer axis oriented at about θ kz = 20 • with respect to the polar axis k = a, resulting in a high diagonal acentric order parameter of cos 3 θ kz = 0.83 (see also Figure 1f), which is close to optimal for the electro-optic and THz photonics applications. Due to the high diagonal acentric order parameter, DAST crystals are strongly anisotropic with a refractive index difference n 1 − n 2 > 0.5 in the visible and infrared wavelength range [22,23,48]. The dispersion of the refractive index is essential for phase-matched THz-wave generation and scales with the absorption band (electronic resonances) of materials according to the general Kramers-Kronig relation. DAST crystals are well suited in particular for applications at around telecommunications wavelengths, with a material absorption that is smaller than 0.1 mm −1 at 1300 nm and 1550 nm (see Figure 5b). As discussed later, this range is also best suited as pump wavelength range for THz-wave generation and detection with DAST.
DAST has very high NLO coefficients and electro-optic coefficients, as listed in Table 1. The highest low-frequency electro-optic tensor components r ijk measured in DAST are shown in Figure 6a. The measured dispersion has been modeled by the two-level model according to Equation (9) and is represented by solid or dashed curves in Figure 6a [23]; the deviation at shorter wavelengths is due to the limits of the simple two-level model close to the absorption edge as discussed below Equation (10). The particular orientation of the NLO stilbazolium chromophores shown in Figure 4a defines the main diagonal and non-diagonal components of the macroscopic NLO susceptibility.  DAST [22,23], DSTMS [52], OH1 [39], HMQ-TMS [53] and BNA [56].

Alternative Ionic Organic NLO Crystals
Theoretical evaluations additionally show that the upper limits of NLO susceptibilities in organic crystals have not been reached yet [11,31]. This is one of the reasons that the research and development of novel NLO organic crystals is continuing. Another reason for the reinforced development during the last couple of years is a requirement for different optical and THz properties that would allow applications in different wavelength regimes, as discussed in details in Section 3.4. Several salt-type organic crystals, i.e., with NLO cation and a counter-anion that promotes non-centrosymmetric packing, have been identified with similar or even superior NLO properties than DAST [32,49,50,[60][61][62][63][64][65]. For photonics applications, not only high nonlinear and appropriate linear optical properties are required. Many of the novel crystals may be very challenging to grow into suitable samples for applications or even simple characterization. Therefore, it is also important to develop organic NLO crystals with easier crystal growth procedures [44,66] or better crystal processing possibilities. For example, a possibility for cleaving as-grown crystals to produce samples with a desired thickness and optical-quality surfaces is very attractive because the expensive and time-consuming cutting and polishing may not be necessary anymore [47,67]. In addition, various crystal engineering approaches are being developed that can lead to crystals with a desired final morphology [68][69][70][71]. Two of the most successful organic NLO salt crystals developed during the last decade, considering their properties and growth/processing possibilities as well as reports on applications, are DSTMS and HMQ-TMS, which are discussed in the following subsections.
DSTMS [44] is structurally very similar to DAST: it has the same NLO cation chromophore as DAST, but the counter-anion is slightly modified (from 4-methylbenzenesulfonate to 2,4,6-trimethylbenzenesulfonate). The additional methyl groups on the anion of DSTMS do not considerably affect both the crystal structure and the nonlinear optical properties (see Figure 4 and Table 1). However, DSTMS has an important advantage for the crystal growth because the solubility in methanol is twice as large as in DAST [44]. In addition, the strongest THz phonon absorption peak at about 1 THz is in DSTMS only about half as that of DAST [72]. DSTMS has been employed by several research groups for THz photonics applications [72][73][74][75][76][77][78][79][80][81].
Another ionic crystal HMQ-TMS has a quinolinium-based ionic-cation core structure. The quinolinium derivatives have been recently identified as very attractive for NLO applications, resulting in non-centrosymmetric crystals with a relatively high probability [15]. HMQ-TMS has been first reported in 2013 [47] and has been already employed in different THz photonics schemes [53,54,71,[82][83][84]. The quinolinium chromophores in HMQ-TMS are oriented almost perfectly parallel (see Figure 4d) and therefore the diagonal NLO response is maximized. Although the NLO properties of HMQ-TMS may be comparable or slightly lower compared to DAST (see Table 1 [15,54]), a big advantage of this crystal is a possibility to obtain high optical-quality samples by a simple cleaving method [47]. HMQ-TMS has a blue-shifted absorption edge (see Figure 5) and different optical and THz properties compared to DAST and DSTMS, so it may be preferred for specific applications.

Non-Ionic Organic NLO Crystals
Besides organic salt crystals, the non-ionic organic NLO crystals have also gained considerable attention during the last decade for integrated electro-optic applications and THz photonics. The configurationally locked polyene crystal OH1 [45] has become very attractive due to the electro-optic coefficients as high as for ionic DAST, combined with a blue-shifted absorption edge and different THz properties [39,85]. OH1 exhibits a relatively high crystal symmetry (orthorhombic, see Figure 4c) compared to the majority of other highly nonlinear organic crystals (monoclinic or even triclinic), which simplifies the handling and characterization. OH1 has very good crystal processing possibilities from solution, including high optical quality bulk crystals [46,68,71] and large-area single crystalline thin films on substrates [86]. High-quality OH1 thin films can be also grown by the physical vapor deposition [87]. OH1 has been found particularly effective for THz-wave applications [85,[88][89][90][91][92][93][94][95][96][97][98].
Another non-ionic NLO organic crystal used for THz-wave generation by several groups is BNA (N-Benzyl-2-methyl-4-nitroaniline) [55]. BNA is based on the nitroaniline derivatives developed in the early days of organic NLO research [99]. Because of the yellow color of these crystals, due to the considerably blue-shifted absorption maximum (see Table 1 and Figure 5a), considerably lower NLO coefficients compared to the red-colored materials such as DAST are expected, based on the common nonlinearity-transparency trade-off (see Equation (11)). However, the measured NLO coefficients of BNA are surprisingly high, see Table 1, with the measured diagonal NLO coefficient of d 333 = 234 pm/V at the fundamental wavelength of 1064 nm [56], which is however resonantly enhanced. From the measured dispersion of the refractive indices and using a two-level dispersion model according to Equation (10), one can estimate d 333 ≈ 92 pm/V and d 333 ≈ 72 pm/V for BNA at 1500-nm and 1900-nm fundamental wavelength, respectively. Considering the two-level dispersion model for both the NLO and the electro-optic coefficient, we can estimate the highest electro-optic coefficient of BNA is r 333 = 24 pm/V at 1538 nm [100]. Although some experimental results suggest different values of the corresponding NLO and electro-optic coefficients [58,59,100], BNA crystal has been very successfully used for THz-wave generation applications, in particular in the pump wavelength regime at around 800 nm that is less optimal for high-nonlinearity red materials [101][102][103][104][105][106][107][108].
Several research groups worldwide identified many other interesting groups of ionic and non-ionic crystals recently, but have not yet been extensively used for THz applications. Some more crystal types are mentioned in Section 4, where we discuss different organic THz sources. For recent reviews on the development of novel organic materials see [15,30]. On the other hand, also some older organic NLO crystals with a lower nonlinearity may be interesting for THz photonics, which was not considered at the time of their development more than two decades ago. For a comprehensive review of these crystals, see, e.g., Ref. [11]. Only a few of these crystals, e.g., COANP (2-cyclooctylamino-5-nitropyridine) [109], DANPH (4-dimethylaminobenzaldehyde-4-nitrophenylhydrazone), MNBA (4 -nitrobenzylidene-3acetamino-4-methoxyaniline), MTTNPH (5-(methylthio)-thiophenecarboxaldehyde-4-nitrophenyl hydrazone), PNP (2-(N-prolinol)-5-nitropyridine) have been investigated for their THz properties [110].
The natural morphology of the crystals is also very important for applications because the NLO response critically depends on the orientation of the interacting fields with respect to the crystal axes and therefore an appropriate sample orientation is essential. For crystals with a natural morphology that does not allow a direct use for applications, a very cumbersome and low-yield cutting procedure is most often needed. A big part of the success encountered by DAST, DSTMS and OH1 crystals is due to the very favorable natural morphology [43,44,46] that allows for direct use in THz-applications, because of the plate-like morphology with the largest plane containing the polar axis, as illustrated in Figure 4. Another interesting case is the one of the HMQ-TMS crystal illustrated in Figure 4d, where a natural cleavage plane allows for relatively simple cleaving of samples with a desired orientation and thickness, which additionally results in high optical-quality parallel surfaces of the samples, without the need for additional polishing [47].

THz-Wave Generation and Detection via NLO Processes
In the early years of organic NLO crystals development, the most promising application considered was second-harmonic generation from the infrared radiation emitted e.g., from widely available infrared diode lasers to the visible radiation, motivated in particular by the very large NLO coefficients of organic materials, one order of magnitude larger compared to inorganic alternatives. Due to the huge progress in visible lasers development and the problems with photostability for second-harmonic generation, since the generated second-harmonic wavelengths are typically close to the absorption edge of these crystals, the interest has been moved to high-speed electro-optic applications [6,25,26] and in particular to THz photonics applications discussed in this review, for which photostability is usually high enough even for commercial applications. In this section, we summarize the theoretical basis of THz-wave generation relevant for organic crystals, the difference-frequency generation and the optical rectification and then present some examples of organic THz sources in Section 4.

Difference-Frequency Generation
To induce a nonlinear polarization oscillating at a THz (angular) frequency ω THz based on difference-frequency generation, the interacting electric fields E(ω 1 ) and E(ω 2 ) should oscillate at closely spaced frequencies ω 1 and ω 2 , so that their difference ω 3 = ω THz = ω 1 − ω 2 is in the THz frequency range. In order to efficiently generate ω THz , the phase matching condition ∆k ≈ 0 should be possibly well satisfied where k THz is the wave vector of the THz-wave with the magnitude of |k THz | = 2π n THz /λ THz and k i (i = 1, 2) are the wave vectors of the interacting optical fields with the magnitude of |k i | = 2π n i /λ i . n THz is the refractive index of the NLO material at the wavelength λ THz of the generated THz field, while n i (i = 1, 2) are the refractive indices of the NLO material at the wavelengths λ i (i = 1, 2) of the two interacting fields E(ω 1 ) and E(ω 2 ), respectively. The phase matching condition is often quantified by the so-called coherence length l c , which defines the maximum length of the NLO material, over which the generated field increases with the thickness In the simplest case of the collinear difference-frequency generation, the three interacting waves propagate in the same direction and the corresponding wave vectors are all parallel. The resulting coherence length is then given by This expression can be additionally simplified for relatively small differences in frequencies of the two pump fields (up to several THz) and in case the dispersion of the refractive index is relatively small, e.g., like in the infrared range for the organic crystals shown in Figure 5a. In this case, the dispersion in the optical range can be considered as n 2 = n 1 + (∂n/∂λ)(λ 2 − λ 1 ). This leads to the following coherence length for THz-wave generation where n g = n 1 − (∂n/∂λ)λ 1 is the group index of the optical wave. For efficient generation of THz waves by difference-frequency generation, besides a high χ (2) material susceptibility, the most important parameter is the low refractive index mismatch ∆n = n THz − n g between the generated THz and the pump optical waves. Organic materials have an important advantage compared to standard inorganic NLO materials such as LiNbO 3 because of the relatively weak dispersion of the refractive index between the optical and the THz frequency range, as discussed in Section 2.1. The phase-matching condition for THz-wave generation is in organics almost naturally satisfied, while for inorganic NLO crystals such as LiNbO 3 , more complex interacting-beam orientations are needed, e.g., tilted-pulse-front configuration [111]. The particular coherence lengths l c for organic NLO crystals are discussed in Section 3.4. For difference-frequency generation in case of perfect phase-matching and neglecting the pump-light absorption, the visible-to-THz conversion efficiency is given by [20] where is the NLO coefficient for THz-wave generation that scales with the electro-optic coefficient r, L the length of the NLO material, I 0 the pump intensity, α THz the absorption coefficient at the THz frequency, n o and n THz the refractive indices at the pump optical and the generated THz frequencies, respectively. Equation (16) shows that, besides phase-matching and minimal THz absorption α THz , the main material figure of merit for THz-wave generation FM THz can be defined as [6] FM THz = d 2 THz n 2 o n THz = n 6 o r 2 16n THz .
(18) Table 2 shows these parameters for several inorganic and organic NLO crystals at a wavelength that is as close as possible to phase matching. Conventional inorganic semiconducting crystals used for THz-wave generation and detection such as ZnTe and GaP have much lower electro-optic coefficients compared to best inorganic electro-optic crystals such as LiNbO 3 , but are anyway widely used because of the possibility for phase matching that is not directly possible with LiNbO 3 and other inorganic crystals with a high electro-optic coefficient. The organic crystals DAST, DSTMS, and OH1 show the largest figures of merit and they can be phase matched using pump lasers at telecommunication wavelengths in the range of 1300-1560 nm. The particular phase matching properties at different pump wavelengths and the generated THz frequencies are discussed for these crystals in details in Section 3.4 and examples of THz sources based on difference-frequency generation in organic crystals are given in Section 4.1. Table 2. Selection of organic and inorganic NLO materials that have been investigated for THz-wave generation and their most relevant parameters. Where possible, the parameters are given close to the optimal phase matching range: refractive index n o and the group index n g at the pump optical wavelength λ; refractive index n THz in the THz range; the electro-optic coefficient r; the NLO coefficient d THz for THz-wave generation according to Equation (17); the figure of merit FM THz for THz-wave generation according to Equation (18).

Optical Rectification
While for difference-frequency generation two pump-optical beams at closely spaced narrow single frequencies are used, THz-wave generation by optical rectification employs a single pump pulse of a very short pulse length. An ultrashort laser pulse in the femtosecond range induces a quasi-static polarization in χ (2) -NLO materials through optical rectification, which acts as a source of a THz field. A short laser pulse has an intrinsically broad bandwidth, i.e., a transform-limited laser beam with a pulse length of 10 fs or 200 fs has a bandwidth of roughly 40 THz or 2 THz, respectively. Different frequency components in such a pulse can mix with each other in an NLO material, producing a broadband THz pulse. In first approximation, considering non-depleted pump, negligible optical and THz absorption, plane-wave approximation and at perfect phase matching, the intensity of the generated THz wave is in frequency domain given by [6,112] I(ω THz , L) = 1 2 0 cn THz |E(ω THz , L)| 2 (19) ≈ d 2 THz I 2 0 ω 2 THz L 2 2 0 c 3 n 2 o n THz (20) and the conversion efficiency by where the same figure of merit FM THz as the one defined for difference-frequency generation in Equation (18) and listed in Table 2 is relevant. An analytical solution has been also derived for a more general case, considering both optical and THz-wave absorption and non-perfect phase matching, but still considering a non-depleted pump approximation [112,113]. In the plane-wave approximation and neglecting surface contributions, the generated electric field E THz (ω, λ, L) at the angular THz frequency ω THz = ω pumped by an optical pulse with the central wavelength λ is given by [112,113] where I(ω) is the Fourier transform of the pump optical pulse I(t) and L the material thickness.
Besides the NLO coefficient for THz-wave generation d THz , the following NLO material parameters are important: the refractive index n THz and the absorption coefficient α THz at the THz angular frequency ω, as well as the refractive index n, the optical group index n g and the absorption coefficient α at the optical wavelength λ. The second term in Equation (22) has a unit of length and its maximum is reached at exact phase matching (also called velocity matching) n THz = n g in case of negligible absorption. Its absolute value is referred to as the effective generation length L gen (ω, λ, L) and is a useful parameter to analyze the thickness dependence of the generated THz spectrum, which strongly depends on the material dispersion [112]. A simpler figure of merit resulting from the expression (22), which however neglects the optical and THz absorption effects, is again the coherence length l c , which is given by the same Equation (15)  The most common THz-wave generation scheme based on optical rectification employs a short fs pump pulse to produce a broadband THz pulse. Often, however, a single frequency THz source is desired, such as produced by difference-frequency generation for example. A tunable narrowband THz source can be also obtained by using a beating signal of two linearly chirped fs pulses with a variable delay as pump on a THz generator [114]. By controlling the length of the chirped pulse and the delay between the beating pulses, one can tune the frequency and the bandwidth of the generated narrowband THz pulses [114]. This method has been recently also demonstrated using organic THz sources [78,82,115].

THz-Wave Detection
The detection of the generated THz waves is often as challenging as THz generation. There are several possibilities used for THz-wave detection, which also depend on the type of the THz source. We mention here the approaches that are most often used for organic crystalline THz sources. For difference-frequency generation systems that produce narrowband THz pulses, most often incoherent detectors such as Golay cells, bolometers and pyroelectric detectors are used. These detectors are commercially available and they can be directly used to determine the intensity of the generated THz radiation. For systems based on optical rectification, one can additionally use coherent detection methods, which are able to detect the amplitude and the phase of the oscillating THz field, i.e., the exact time dependence of the generated electric field E THz (t). The coherent detection methods can give a lot of additional information and are therefore very attractive for applications. We briefly describe the most common methods that are used for coherent detection of THz fields produced by optical rectification in organic NLO crystals. Some of the often employed coherent detection methods for the optical rectification systems are illustrated in Figure 7.  Photoconductive antennas measure the transient photoconductivity induced by the probe optical field. When the generated THz electric field overlaps with the probe optical beam on a photoconductive antenna, a current is induced across the electrodes of the antenna and its magnitude and sign are proportional to the magnitude and the sign of the THz field. By delaying the probe pulse with respect to the THz pulse, a time-domain trace of the THz electric field E THz (t) is obtained. Most commonly, the response of the photoconductive antennas is limited to a few THz only, but, using special antenna designs, broadband detection has also been demonstrated [116].

Electro-Optic Sampling
When fs optical pulses with relatively high pulse energies (in the order of 1 mJ) are available, the detection is even possible in air by the so-called air-breakdown-coherent-detection (ABCD) [117,118].
Most commonly used technique for detecting the THz waves based on NLO crystals is the electro-optic sampling. This method employs electro-optic crystals with zero birefringence such as ZnTe or GaP [119]. When the THz field illuminates the detection electro-optic crystal, it induces a birefringence through the conventional electro-optic effect according to Equation (3). The induced birefringence changes the polarization state of the probe beam, which is usually detected using a combination of a quarter-wave plate, Wollaston prism, and two balanced photodiodes [119] (see Figure 7).
For highly birefringent materials such as organic NLO crystals, conventional electro-optic sampling detection cannot be used. Several alternative techniques have been proposed for birefringent materials [80,[120][121][122][123]. A relatively simple approach is the so-called THz-induced lensing technique, which also relies on the electro-optic effect [121]. In this method, the refractive index changes induced by a focused THz electric field through the electro-optic effect are exploited. The spatially varying (focused) THz field induces a spatially varying refractive index change ∆n(x, y, t) = − 1 2 n 3 rE THz (x, y, t).
This time-varying refractive index profile acts as a time-varying optical lens and leads to focusing/defocusing of the probe beam, which can be detected by various common optical detection techniques [121,124].
For the above mentioned electro-optic-based detection methods, the electro-optic sampling and the THz-induced lensing, the same phase-matching considerations are valid as for THz-wave generation with the corresponding material, considering the THz frequency limitations and the probe beam wavelength; i.e., for efficient detection, the THz refractive index n THz and the optical group index n g for the detection pulse should match. In the same way as for detecting with the photoconductive antennas, by delaying the probe pulse with respect to the THz pulse, a time-domain trace of the THz electric field E THz (t) is obtained by electro-optic sampling or THz-induced lensing (see Figure 7). Examples of the measured THz signals using different detection principles in systems based on optical rectification in organic crystals are presented in Section 4.2.

Phase-Matching Characteristics for THz-Wave Generation in Organic NLO Crystals
For efficient THz-wave generation by NLO processes (such as difference-frequency generation and optical rectification), it is essential to achieve phase matching (velocity matching) n THz n g between the fundamental optical and the generated THz waves, as discussed in Section 3. Because of the dispersion of the refractive indices in both ranges, the phase matching condition is strongly dependent on the wavelength λ of the pump optical beam(s) and the frequencies ω of the generated THz beams. The dispersion of the optical group indices n g = n − (∂n/∂λ)λ is shown in Figure 8a next to the dispersion of the refractive indices n(ω) in the THz range (see Figure 8b) for selected organic NLO crystals.
For different materials, there are different ranges where the velocity matching condition n THz n g is reasonably well satisfied and is roughly marked by the shaded thick lines in Figure 8. DAST for example has a stronger phonon resonance near 1.1 THz resulting in roughly two different regions: good velocity matching can be realized either for pump wavelengths between ca. 1000 nm and 1200 nm, yielding THz waves between 0.2 THz and 1 THz, or for pump waves with wavelengths between 1200 nm and 1700 nm for the generation of THz waves above 1.5 THz [112,125]. For crystals with a blue shifted visible absorption edge such as HMQ-TMS, the best velocity matching ranges are also blue shifted, which is very important when selecting an organic THz generator that should be optimized for a certain pump laser source. The optical phonon resonances and molecular vibrations that occur in the THz frequency range (see Figure 3) induce some absorption of the generated THz radiation in organic NLO crystals (see the graphs for α THz ( f ) in Figure 8b), which is also very important when optimizing these crystals for different desired frequency ranges. Note also that the properties in Figure 8 are only given along the polar axis of the selected NLO crystals. Several crystals also have notable off-diagonal χ (2) components (see Table 1). Although these coefficients may be lower compared to the diagonal susceptibility coefficients (employed when both optical and THz beams are polarized along the polar axis), they may allow other velocity matching ranges. For DAST, for example, using the pump beams polarized along the crystalline b axis (see Figure 4) generates THz radiation along the polar a axis based on the off-diagonal χ (2) 122 susceptibility, which is velocity matched for the pump wavelengths close to 800 nm [75,112], an important wavelength of available high-power Ti:Sapphire fs laser sources.
Best conditions for velocity matching can be straightforwardly evaluated if the relevant optical parameters given in Figure 8 are known. The simplest figure-of-merit for the velocity matching is the coherence length l c (λ, f ) given by Equation (15) and is shown in Figure 9 for selected organic materials obtained using the parameters in Figure 8 as a function of the pump optical wavelength λ and the generated THz frequency f = c/λ THz .  The corresponding properties are given for selected organic crystals along their polar axis, which is most often employed in THz applications. The n g (λ) are obtained from the dispersions n(λ) shown in Figure 5a. The THz properties n THz ( f ) and α THz ( f ) are according to the following references: DAST [126][127][128], DSTMS [129], OH1 [91], HMQ-TMS [53,83].
For highly nonlinear organic NLO crystals DAST and DSTMS presented in Figure 9, the best pump wavelength range is at around telecommunication wavelengths, which is very interesting given the availability of very compact and relatively low-cost fs fiber lasers [125]. The generation of THz waves with frequencies at around 1.1 THz is in DAST limited due to a strong optical phonon vibration [130]. In DSTMS (DAST derivative with an isomorphic structure), the optical phonon absorption near 1 THz is suppressed (see Figure 8b), most probably due to a heavier counter anion, and therefore this material is superior to DAST within this frequency range [72]. In a broader THz frequency range up to about 10 THz (for which the THz refractive index data is available), DSTMS shows similar THz velocity-matching properties compared to DAST.
OH1 crystal is a non-ionic crystal and has even higher figure of merit for THz-wave generation (see Table 2) compared to the stilbazolium salts DAST and DSTMS. OH1 has optimum velocity-matching between 1200 nm and 1460 nm for 0.3-2.5 THz, and has no absorption at about 1 THz but at about 3 THz [85]. Above 3 THz, the optimal phase matching conditions for OH1 are shifted to longer wavelengths, which is less optimal considering the available pump laser sources. Laser sources at around 1550 nm are not best suitable for OH1 due to a smaller absorption band in the optical range at this wavelength (see Figure 5b). In any case, OH1 has the highest efficiency among presently known organic crystals, at least in the lower THz frequency range below 3 THz [85]. For the quinolinium salt HMQ-TMS, see Figure 9d, the best phase-matching range is shifted towards shorter wavelengths compared to DAST and DSTMS, which is an advantage for using certain widely pump-wavelength sources at around 1000 nm [53,83,84].
Note that the measurements of material parameters required for evaluating l c may be very time consuming and often also have a limited precision, therefore these figures only give a first approximation for the coherence length and are only available for a few organic NLO crystals. Some other organic NLO crystals, for which phase matching has been studied in more detail, are OHQ-T [131], BNA [103] and COANP [109]. BNA and COANP are, for example, interesting, although having somewhat lower nonlinearity, because their absorption is considerably blue shifted (see Figure 5a for BNA) and therefore they allow for velocity matching at around the Ti:Sapphire laser wavelength of 800 nm [103,107,109].
The velocity-matching evaluation based on the coherence length l c as presented in Figure 9 does not take into account the influence of the optical and the THz-wave absorption. Based on the analytical solution for the THz electric field given by Equation (22), an equivalent parameter l opt has been defined, the so-called optimum crystal length l opt [112], which considers both optical/THz absorption and optical/THz group/refractive indices. The corresponding plots of l opt as a function of f and λ for DAST, DSTMS, OH1 and HMQ-TMS can be found in Ref. [15]. The main conclusions about the best velocity-matching ranges remain, however, very similar to those resulting from Figure 9.
It is usually assumed that the NLO coefficient for THz-wave generation d THz , which is important for the THz-wave generation efficiency beside phase matching (see Equations (16) and (22)), is approximately constant in the THz frequency range. This may not be the case due to the phonon and vibrational modes in this range and the NLO susceptibility in this range can be also considerably frequency dependent similarly to the linear optical properties shown in Figure 8b. A recent experiment shows that d THz may be even strongly enhanced at certain THz frequencies [132], which then additionally affects the final features in the generated THz spectrum in addition to phase matching.
There is a strong need to increase the variety of different NLO crystals that would allow for phase matching for various combinations of available pump optical and the desired generated THz frequency ranges. This is one of the reasons that the development of novel organic NLO crystals is being continued [15,30].

Organic Crystalline THz Sources
Over the last two decades, the interest in organic THz sources has been steadily increasing due to their high THz-wave generation efficiency and the ability to achieve ultra-high bandwidths as well as high THz electric fields. Zhang et al. [133] were the first to demonstrate the generation of THz-waves by optical rectification in organic DAST crystal. His group obtained results of one and two orders of magnitude larger than that from GaAs and LiTaO 3 inorganic crystals, respectively. This work triggered other researchers around the globe and in the following years, there were many reports of new achievements with organic THz sources. In this section, we briefly describe some of the most relevant research results, mostly from the past decade in the area of ultra-high bandwidth and high electric field THz-wave generation with organic crystals. Some of them involve new recently developed materials, while others employ well known organic crystals used in new studies.

THZ Sources Based on Difference-Frequency Generation
By using the process of difference-frequency generation, narrowband and in most cases frequency-tunable THz wave generation in various NLO crystals has been demonstrated. While the accessible THz wavelength range is intrinsically limited by the absorption of the NLO material employed, it is difficult to generate waves above 5 THz efficiently in inorganic crystals due to strong optical phonon resonances. The THz sources based on difference-frequency generation in organic NLO crystals are characterized by a wide tunability range, often beyond 20 THz (see Figure 10), a narrow single frequency output with linewidths that may be below 5 GHz [134] and high THz pulse energies beyond 500 nJ [97]. The typical dynamic range of such systems is about 30 dB (see Figure 10), limited by incoherent THz detection techniques employed (see Section 3.3).
A slightly non-collinear difference-frequency generation geometry may be used to achieve exact phase-matching conditions and improve the efficiency in a desired THz range [144]. Specially cut crystals may be used to optimize the phase matching at other pump wavelengths [147,148]. The possibility of quasi-phase matching has also been demonstrated in case the available dual-wavelength source is out of the best phase-matching conditions of a particular crystal [91]. Intrinsic gaps in the achievable spectrum are observed due to optical phonon and molecular vibrations (see Figure 8b), as shown in Figure 10a,b [106]. An interesting possibility to reduce these gaps is to combine the output generated by difference-frequency generation in different organic NLO crystals to obtain a more homogeneous THz spectrum, as demonstrated for DAST and BNA (see Figure 10c [106]).
In general, it is relatively straightforward to produce narrowband THz waves by differencefrequency generation in organic NLO crystals, which may be most often phase matched without the need for a special crystal cut at certain pump-wavelength ranges [135][136][137][138][139][140][141][142][143][144][145][146]. It is more challenging to develop an appropriate pump-beam pair with a desired separation in the THz frequency range. These systems are presently mostly based on high-power optical parametric oscillators and therefore relatively bulky and sensitive. More compact generation systems can be achieved by using optical rectification discussed in the following section.

THz Sources Based on Optical Rectification
The THz sources based on optical rectification in organic crystals are presently more widely employed compared to the sources based on difference-frequency generation. The THz electric field strength and its spectrum achieved by these sources depend on the characteristic properties of a particular THz generator, its nonlinearity and phase matching (see Section 3.4), as well as on the characteristics of the particular pump laser source through the source term I(ω) in Equation (22). The shorter the pump pulses, the broader the THz spectrum that can be achieved by optical rectification. For example, using 17-fs pump pulses, it was possible to generate a spectrum extending from 0.1 to beyond 25 THz in DAST crystal by optical rectification [153], and using 5-fs short pulses, the spectrum extended even up to 200 THz [154].
The final system specifications for applications such as e.g., needed for THz spectroscopy also strongly depend on the employed detection technique. A simplified example of a possible experimental set-up based on optical rectification is illustrated in Figure 7. The incoming pump laser pulse is split into a pump beam and a probe beam. The pump beam is used to generate a THz pulse in an organic NLO crystal and is then blocked by a filter transparent to THz radiation. For the optics to guide THz beams, parabolic or ellipsoidal mirrors are mostly used and different mirror combinations can be used to optimize the system for different THz photonics applications. The THz beam can be measured directly using incoherent THz detectors or coherently by overlapping THz pulses and probe optical pulses in a detector material, for which several different approaches are commonly used (see Section 3.3). The coherent detection, i.e., the measurement of the THz electric field as a function of time E THz (t), is most often achieved by varying the relative path length of pump with respect to probe pulses with a computer-controlled optical delay line. In some cases, e.g., when different materials are used for generation and detection with different optimal operation wavelength, probe pulses with a different wavelength can be used, which can be achieved e.g., by using a second-harmonic of the pump beam or a direct output of a Ti:Sapphire laser in case an additional optical parametric amplifier is used for the pump beam.
For table-top THz sources, particularly for those aiming at high-power THz pulses that can be used for nonlinear THz photonics (see Section 5.1), conventional femtosecond Ti:Sapphire lasers, combined with high-power optical parametric amplifier systems, are most often employed. Conventional Ti:Sapphire lasers emit at the wavelength of around 800 nm and they can be directly used only for organic crystals that allow for phase matching at this wavelength, such as BNA [107,108], COANP [109] and some other recently developed organic crystals, e.g., crystals based on the 2-(4-hydroxystyryl)-1-methylquinolinium (OHQ) chromophore [168] or the 4-(4-hydroxystyryl)-1methyl-pyridinium (OHP) chromophore [167]. For highly NLO crystals such as DAST, DSTMS (see Figure 9), wavelengths of above 1200 nm are optimal, which is why an additional high-power optical parametric amplifier seeded by a Ti:Sapphire laser is used for these sources to convert the output to longer wavelengths. With such a Ti:Sapphire-based pump system, THz electric fields exceeding 150 MV/m [93] or even beyond 8 GV/m using special focusing optimization techniques have been demonstrated [94]. Figure 11a shows time domain traces achieved with such a system for DSTMS and OH1, where the pump wavelength has been chosen as the optimum for each of the crystals, 1500 nm and 1350 nm, respectively. The spectral shape (see Figure 11b) such as the gaps in the spectrum for OH1, are due to the phonon and vibrational modes in this crystal, which corresponds well to the measured THz properties shown in Figure 8b and the theoretical calculations of phonon modes in OH1 [169]. Figure 11. Intense THz-wave generation using a high-power Ti:Sapphire/optical parametric amplifier system (65 fs pulse duration, 100 Hz repetition rate) in DSTMS (0.44-mm thick, 1500 nm pump wavelength and 3.8 mJ pump energy, blue curves) and in OH1 (0.48-mm thick, 1350 nm pump wavelength and 3.5 mJ pump energy, red curves without and black dotted curves with a 3-THz low-pass filter): (a) THz time-domain signals detected by air-breakdown-coherent-detection and (b) the corresponding amplitude spectra. The green curves in (a) represent the estimated generated electric field E THz (t) considering the mirror losses, with a maximum of 8.3 GV/m for DSTMS and 6.2 GV/m for OH1 [94]. Reprinted from Ref. [94] with permission from Springer Nature.
One interesting alternative to Ti-Sapphire-based pump systems is a high-energy Cr:Forsterite laser emitting at 1250 nm, which is still in the optimal pump range for the state-of-the-art organic NLO crystals (see Figure 9). With Cr:Forsterite femtosecond pump, THz electric fields higher than 1 GV/m in organic crystals DAST, DSTMS and OH1 have been generated [74,170,171].
With the variety of organic NLO crystals that are being developed, the possibility for phase matching should be also soon available for many other available high-power pump fs laser systems. For example, an ytterbium laser operating at 1030 nm [84] has been successfully used in combination with HMQ-TMS crystal, for which the best pump wavelength range is blue-shifted compared to DAST, DSTMS and OH1 (see Figure 9), also yielding high THz electric fields with a peak of about 20 MV/m [84].
In addition to high-power systems, relatively compact and ultra-broadband THz systems based on optical rectification in organic NLO crystals can be realized based on compact and relatively low-cost femtosecond laser sources operating at around 1560 nm, which is about the best phase-matching range for the state-of-the-art organic crystals such as DAST and DSTMS. Since the available pump pulse energies of these fs lasers are relatively small (in the nJ-range, which is six orders of magnitude smaller in comparison with the values up to several mJ for the high-power systems, see, e.g., the values in the caption of Figure 11), these systems do not provide an intense THz power, but they are still very interesting in combination with a sensitive detection technique, providing a high dynamic range needed e.g., for THz time-domain spectroscopy applications. Two of such examples based on a DSTMS organic crystal as THz generator are presented in Figures 12 and 13. In the first example, a highly sensitive photoconductive antenna optimized at 1560 nm wavelength has been used and in the second example, another DSTMS crystal has been used for the THz-induced lensing detection. Using organic crystals as both THz generation and detection materials enables achieving an ultra-broad THz spectral bandwidth beyond 15 THz, profiting from the phase-matching on both generation and detection part.  [79]. The red lines represent signals measured in air with a humidity of 44% and black lines are those obtained in dry-air environment. Reprinted from Ref. [79]. Licensed by Creative Commons CC BY 4.0.
One drawback of broadband systems based on optical rectification in organic NLO crystals is the modulation of the resulting THz spectra (see, e.g., the spectrum in Figure 11b for OH1 and in Figure 13b for DSTMS). This is due to the intrinsic phonon and vibrational modes of the generating materials, resulting in THz absorption peaks as shown in Figure 8b, leading to the absorption of the generated THz waves and varying phase-matching figures. To avoid gaps in the spectrum, a solution based on exchanging between different crystals at different THz frequencies as proposed for the systems based on difference-frequency generation (see Figure 10c) is here not applicable because the whole spectrum is generated and detected simultaneously through the time-domain traces E THz (t).
Therefore, one important challenge for future development of novel organic NLO crystals for THz photonics is to reduce the strength of these vibrations. This can be done through molecular engineering, but this is unfortunately still extremely difficult to control because even a slight chemical change of the constituting molecules can considerably change the crystalline packing and the crystal growth ability [172]. Another recently proposed approach combines two complementary organic NLO crystals in a combined structure as a tandem THz generator [54]. Using a tandem consisting of HMQ-TMS and OH1, it was shown that a certain compensation of phonon-absorption gaps in the generated THz spectra is possible for optimized crystal combinations [54].

Applications of Organic THz Sources
During the last decade, THz systems based on organic single crystalline THz sources reached many areas of applications from fundamental science to real-world applications. The research and practical applications such as THz spectroscopy and THz imaging in laboratories worldwide performed by organic NLO crystals are still more exceptional in comparison with competitive THz systems based on photoconductive antennas or inorganic NLO crystals that are more common, also because photoconductive antennas and inorganic NLO crystals are much more widely available. Nevertheless, the use of organic crystalline THz sources in fundamental science and development is moving to other interdisciplinary fields [174][175][176][177][178][179], which indicates a great potential for the future transfer of organic-crystal based THz spectroscopic and imaging systems into industrial applications. According to the THz science and technology roadmap prepared by Dhillon et al. [1], one can expect that organic THz sources will enter other fields of applications similarly as THz systems based on photoconductive antennas or inorganic NLO crystals. Furthermore, the last experiments showed that organic NLO crystals in a combination with photoconductive antennas provide good performance suitable for various applications [79]. Due to the possibility to reach very high THz electric fields with organic NLO crystals, several demonstrations of fundamental THz nonlinear photonics have been possible (see Section 5.1) using these sources, especially interesting due to the possibility to control matter with THz beams. The following sections summarize the past and current directions of use as well as provide an opinion on the future challenges and opportunities of organic crystalline THz sources.

THz Nonlinear Photonics
Recent progress in generation of intense THz pulses opened a range of new possibilities to study nonlinear physical phenomena. Interaction between mater and electric field produces various linear and nonlinear effects that can modify material properties. When a high external field is applied, the induced polarization inside the medium is a sum of the linear and nonlinear components (similarly as in the optical range described by Equation (2)), giving rise to many nonlinear effects [180,181], such as coherent THz manipulation of quantum states, high-order harmonic generation, nonlinear optical processes, and nonlinear transport phenomena in solids [182]. Various THz generation techniques are employed to achieve intense THz fields needed for nonlinear effects and most of them are based on optical rectification where either LiNbO 3 [182][183][184][185][186] or organic NLO crystals are used. Here we summarize some of the reported examples of nonlinear THz research based on organic crystalline THz sources.
DSTMS THz emitter crystal was used for generating THz pulses with a high peak electric field of about 300 kV/cm for a system enabling the THz Kerr effect spectroscopy and two-dimensional terahertz-terahertz-Raman spectroscopy [77,187,188], demonstrating a 2D experimental technique for measuring the nonlinear response of liquids in the THz field. The terahertz-terahertz-Raman pulse sequences that are sensitive to intermolecular modes (e.g., molecular orientational alignment) were employed to initiate intramolecular vibrational modes of liquids (see Figure 14). Vibrational coherence arising from intramolecular modes were observed in liquids such as diiodomethane at 3.66 THz and in simple halogenated liquid carbon tetrachloride at 6.50 THz [187]. These experiments open the door to THz spectroscopy that can study the dynamics of condensed-phase molecular systems such as the ability of liquids to solvate different compounds, the dynamics of biological macromolecules such as polypeptides, protein folding and DNA internal conversion as well as for amorphous solids and hydrogen-bonded liquids.
Time (ps) THz Kerr Effect Signal (a.u.) Figure 14. Measured THz Kerr-effect signals, generated by a terahertz-terahertz-Raman pulse sequence, of various liquids (black curves) [187]. Vibrational coherences are visible as oscillations in the measured signal. The bottom red trace is the square of the applied THz electric field and the blue trace below the tetrahydrofuran data shows the empty cuvette background response. Reprinted from Ref. [187] with permission from AIP Publishing.
Instantaneous Kerr nonlinearity and the retarded alignment of air molecules CO 2 , N 2 , and O 2 when triggered by intense THz pulses have been demonstrated by Shalaby et al [189], for which intense single-cycle THz pulses generated in a DSTMS crystal have been employed [94]. The same group demonstrated the visualization of broadband THz radiation in the 5-13 THz frequency range, based on nonlinear THz absorption mechanism in the standard silicon charge-coupled (CCD) sensor [190]. Intense THz pulses have been shown to induce heat-induced demagnetization in the ferromagnetic nickel [191] and to change dielectric properties and simultaneous electronic and magnetic response in a cobalt thin film [192]. Chefonov et al. investigated n-doped silicon with high-field THz transmission experiment where they observed huge transmittance enhancement of approximately 90 times under high THz field strengths generated by a DSTMS crystal [161]. The observable high increase in the transmission of n-doped silicon in the THz range by almost two orders can lead to the development of novel THz absorbers and THz field controlled electronic devices [161]. An ultrafast THz-induced optical modulator for pulse switching in a thin diamond window has also been demonstrated [193].
Vicario et al. reported on extreme laser spectral broadening and a change in optical transmission in GaP induced by a strong THz field [194]. DSTMS crystals were used in the experiment with thin aluminum and nickel films, where damage produced by a single and multiple THz pulses was investigated [171,195]. The observed damage resulted in ablation and through holes as well as delaminations and cracking near the damage threshold. The real-time room-temperature THz imaging developed by Fan et al. [196] was demonstrated on an aluminum foil by nonlinear frequency up-conversion in a DAST crystal. The results showed higher sensitivity in comparison to a commercial THz camera. Furthermore, analytical and numerical models [197] showed that the rapid heating of aluminum film with thicknesses of several tens of nanometers by the high-energy THz pulse causes a significant rise of effective electron collision frequency. Consequently, the increase in nonlinear transmission of the thin film as well as noticeable changes in reflectivity can also be observed. All of these recent examples demonstrate a very high potential of intense organic THz sources for nonlinear THz photonics.

General THz Spectroscopy and Imaging
Most often employed THz systems for THz spectroscopy and imaging applications are the so-called THz-TDS (THz time-domain spectroscopy) systems, which are based on optical rectification and can simultaneously measure the amplitude and the phase of the THz electric field and therefore, for example, the absorption coefficient and the refractive index in the THz range. THz-TDS systems can be also used for non-destructive sample-thickness or coating-thickness measurements [198]. One of the additional features of THz-TDS systems is the possibility of terahertz pulsed imaging. Terahertz pulsed imaging can be performed in transmission or in reflection system geometry by raster-scanning of the sample at the focal plane of a THz-TDS system. In the time-domain, terahertz pulsed imaging can provide 2D amplitude images based on the maximum amplitude or the peak-to-peak amplitude of the transmitted (or reflected) THz electric field as well as 2D phase images by considering the arrival time of the waveform. THz spectroscopic images are obtained by the Fourier transform of time-domain signals at each sample position. Spectroscopic THz imaging and spatial distribution mapping provide visualization and localization of individual chemical compounds within heterogeneous samples, e.g., pharmaceutical pellets. Thus, the complete spectroscopic information about each pixel is obtained, which allows besides detection also identification as well as classification of materials. The chemical analysis is possible at frequencies higher than 500 GHz [199], where THz spectral responses express significant characteristics that can be assigned to inter and intramolecular vibrations of individual chemical substances. Additionally, data of spectroscopic THz imaging can be further enriched by the use of different chemometrics methods. Among them, the most used is principal component analysis, where a series of spectra are simultaneously compared through measuring the covariance [200]. The method is also known as chemical mapping and is often applied in the pharmaceutical industry [201,202].
There have been already many examples of broadband THz spectroscopy measured with THz system based on organic NLO crystals. For example, optical properties in the THz range of organic NLO crystals themselves as those presented in Figure 8b and many other crystals have been measured by organic-based THz-TDS systems [53,83,91,109,126,129]. The refractive index and absorption coefficients have been also measured using THz systems based on difference-frequency generation in organic crystals [127,169].
A system based on organic NLO crystal DSTMS as a THz emitter and a photoconductive antenna as a THz receiver operating in the 0.2-8 THz range [79] was used for measuring a lactose film and two metamaterials, U-shape split-ring resonators and rectangle hole arrays. The obtained results coincide with previous results measured by Fourier-transform infrared spectroscopy or inorganic-based THz systems.
The refractive indices and the absorption of single crystals of aluminum nitride were measured in the 1-8 THz by THz-TDS based on DSTMS crystals [203]. The selected semiconductor has promising optical and optoelectronic properties that can be exploited for optical applications from the far infrared to deep ultraviolet spectral region.
Kawase et al. [204] performed several exploratory experiments with high-resolution time-of-flight THz tomography using DAST crystal for THz waves generation. They measured the thickness of Teflon films as a function of actual thickness. THz tomography images of multilayered structures of semiconductor device sample were recorded where GaAs layer was approximately 2 µm thick. This method should be valuable for car paint inspection. By using a metallic mesh, they tried to achieve imaging with high sensitivity in the THz region. They were able to detect 1 µm thickness differences in ultrathin PET film which is about 300 times thinner than the wavelength in the THz frequency range [204].
In the following subsections, we present further examples of THz spectroscopy and imaging using organic NLO crystals as THz sources in specific industrial sectors.

Pharmaceutical and Chemical Industry
The potential of various THz systems in the pharmaceutical industry was explored by several studies. For pharmaceutical use, mainly molecular structure, polymorphic forms [205], cocrystals [202], hydrates [206], dehydration processes [207], tablet coating thickness [208], concentration of active pharmaceutical ingredient (API) [178], spatial density distribution [209,210] and impurities [211] were analyzed. THz spectroscopic analysis based on organic NLO crystals was performed to distinguish between the absorption spectra of polymorphs I, II, and V of piroxicam, particularly in the regions below approximately 2 THz and above approximately 3.8 THz [177]. The organic compounds like piroxicam can come in different polymorphic forms that typically differ in their physicochemical properties, which also importantly affect their pharmaceutical performance. Since THz waves are sensitive to both intermolecular and intramolecular vibrations, THz-TDS can be used to investigate the crystalline state in different chemical species e.g., polymorphism. The same system based on DSTMS crystals was used to analyze hormone-based drug melatonin and its pharmaceutical product Circadin in the THz frequency range between 1.5 THz and 4.5 THz [178]. In the selected spectral range, melatonin showed characteristic spectral features at 3.21 THz and a weaker one at 4.20 THz allowing a quantitative analysis within the final product. Furthermore, spectroscopic THz imaging of different concentrations of Circadin and melatonin as an API within pellets was also performed, in order to verify the spatial recognition of different substances [178]. Another study with a DSTMS-based THz-TDS system was performed on paracetamol, where commercial pellets from two manufacturers were compared [212]. In both samples, the same characteristic spectral peaks were observed at 2.1 THz, 3.0 THz, 3.6 THz and 4.5 THz. The studies indicated that THz spectroscopy and imaging using organic THz sources could be a complementary tool to Raman and Fourier transform infrared spectroscopies, especially in applications of identification and quality control in the pharmaceutical industry.
In different chemical industries ranging from construction materials, paints, glasses, ceramics, adhesives, composites, paper products, polymers and plastics to animal feeds, food ingredients, pharmaceuticals, and cosmetics, powders are main process entering material used as fillers, pigments or additives [213]. Its usefulness and effectiveness in manufacturing processes are defined by several physicochemical properties such as granulation, particle size distribution, hardness, morphology, specific surface area, polymorphism and chemical purity. One of the most abundant raw material calcium carbonate (CaCO 3 ) was characterized by THz spectroscopic and imaging system operating with organic NLO crystals (see Figure 15) [174]. Several parameters such as CaCO 3 concentration in a polymer matrix, its granulation, and chemical treatment were analyzed based on the analysis of peak absorbance amplitude at around 3.2 THz, the appearance of additional spectral features such as shoulder and small peaks on the first derivative of absorbance spectra, and the frequency position and shift. In the last experiment, the calcite samples were treated by various chemicals. In that case, the frequency position of the characteristic spectral peak at 3.2 THz was shifted to higher frequencies although the powder granulation was the same. The less treated sample, i.e., natural calcite, had the most similar peak position as the pure calcite sample, whereas the samples with additives that modified the particle surface have the highest peak shift in comparison to the pure chemical. In addition, THz-TDS and amplitude THz pulsed imaging were also applied to distinguish CaCO 3 from other two calcium compounds, calcium oxide and calcium hydroxide, which appear during the carbonation/decarbonation production process. With spectroscopic THz imaging, a spatial distribution map was obtained, indicating regions with a various concentration of CaCO 3 in a polymer matrix. The results of this study demonstrated that the THz system based on DSTMS crystals has great potential as a rapid analytical tool for qualitative and quantitative analysis of powdered materials. The different size of particles is presented by circles of varying diameter and the error in the peak position is presented by a value of standard deviation (SD) next to each circle [174]. Reprinted from Ref. [174] with permission from Elsevier.

Defense and Security
Terahertz waves can, at least partially and at certain THz frequencies, penetrate various non-metallic materials such as fabrics, wood, paper, cardboard and plastics. Therefore, they can be used in defense and security screening to uncover the concealed illicit drugs or explosives within packages. Since many dangerous materials have unique spectral fingerprints in the THz range, it is possible to combine both, detection by THz imaging and subsequent THz spectral identification with THz-TDS. Organic crystalline THz sources can be used in the same way as inorganic sources for close-by-inspection of envelopes or suspicious powders hidden under clothes. THz time-domain spectroscopy and imaging system based on DSTMS was used to study the effects of various hiding techniques with different paper and textile barriers based on analyzing the spectral features for drugs and explosive simulants [214]. Besides deposited and impregnated powder samples, mixtures with polyethylene in a form of a pellet were also analyzed hidden behind or within textile or paper barriers simulant. The study confirmed that rapid detection and identification of concealed simulants is possible in the frequency range from 1.5 THz to 4.0 THz by using the spectral peak analysis method and a priori known characteristic spectral features from the database [214]. The DSTMS-based THz system was also used in a comprehensive study where different THz technologies were compared for detection and identification of 18 explosives in a pellet and powder forms [215]. The selected DSTMS crystals allowed access to higher frequencies compared to the inorganic-based systems (see Figure 16). In some cases, the results showed that minor absorption peaks, bands and shoulders that are visible in the respective frequency region in the Fourier-transform infrared spectroscopy data can appear more prominent in THz-TDS measurements [215].

Food and Agriculture
THz waves offer a non-ionizing non-destructive testing method to detect foreign bodies within food products with low water content such as chocolate or powdered eatables, e.g., sugar or flour. In addition, the chemical analysis of food samples such as edible oils is also possible with THz spectroscopy. THz interaction with food samples provides intramolecular information linked with the existence of certain characteristic bonds indicating the presence of certain molecules or groups of molecules as well as intermolecular information such as hydrogen bonding, alignment of molecules, amorphous or crystalline structure [216]. Because THz waves are very sensitive to water molecules, THz-TDS is an excellent tool to determine the water content within the food samples such as butter, vegetable oils or flour. A DSTMS-based THz spectroscopic system was used to obtain the absorption coefficient of five different edible oils from different plants in the frequency range from 1.5 THz to 3.5 THz [179]. By using additional chemometrics, also distinguishing between different edible oils was possible. The use of chemicals such as additives and hormones in food production is becoming a critical food safety issue not only for regulatory agencies but also for consumers. Some studies on hormones have already been performed using THz spectroscopy. One such example performed by a DSTMS-based THz system investigated the melatonin [178], which can be used as a dietary supplement. The experiment confirmed that the melatonin has a unique spectral feature at around 3.2 THz, which can be used to detect and quantify it within the dietary samples in order to control its authentication. Similarly, other chemicals such as pesticides and antibiotics can also be analyzed. Since pesticides have characteristic spectral response in the THz frequency range, the THz-TDS method can be used to detect pesticides residue in crop samples [216]. Furthermore, qualitative and quantitative analysis of food additives or antibiotics due to veterinary drugs can be performed by THz-TDS.

Building and Construction Sector
Although the reports on possible applications of THz technology in different sectors constantly increase, the construction and building industry lags behind them. However, several studies demonstrated that THz spectroscopy and imaging have promising potential and provide many opportunities for applications in construction and building materials characterization [217]. For instance, structural analyses of various thermal building insulation materials such as polymer foams, an aerogel composite, and a microencapsulated phase-change material by using a DSTMS-based THz-TDS and THz pulsed imaging were performed [175,176,218]. The insulating materials ensure low thermal conductivity and thermal transmittance to buildings. These characteristics are highly dependent on the thickness and structural properties of materials (e.g., homogeneity, porosity, humidity, density, inclusions, defects, etc.). Studies revealed that the absorption coefficient of polymer foams in the THz range is inversely related to the thermal conductivity of these foams [175]. The obtained THz amplitude images (see Figure 17) clearly visualized different manufacturing imperfections and internal structures such as beads structure, solvent distribution, air gaps and fibre direction in layers of aerogel blanket [175,176]. Based on the observation of additional spectral features during the solidification and liquefying of a phase-change material, THz spectroscopy was demonstrated to be a suitable method to characterize specific properties related to the freeze-thaw process [175]. This is especially important when phase-change materials are incorporated within a polymer matrix in order to increase the energy efficiency of insulation materials. In addition, a periodical fibrous structure of an aerogel composite was also observable in THz spectra, where several spectral peaks were present separated by a distance corresponding to the diameter size of fibers [218]. These examples show that THz-TDS and THz pulsed imaging based on organic THz sources have a potential to become a suitable non-destructive testing method when evaluating the existing building materials as well as when developing novel high-performance fibrous, cellular or granular materials. THz-based methods can provide important construction-material information, especially when combined with other non-destructive testing techniques using ultrasound, microwaves, infrared waves or X-rays. Figure 17. (a) signal processing for particle analysis of THz image determining the voids within a polymer foam [176]. Detection of imperfections (b) and internal structures (c) by THz imaging [175]. Reprinted from Ref. [176] and Ref. [175] with permission from Elsevier.

Future Challenges and Opportunities for Practical Applications
In the last decade, the rapid development of novel THz sources and detectors has opened additional possibilities of creating innovative imaging and sensing systems, allowing in situ measurements and real-time analysis of different materials. Among these, THz systems based on organic NLO crystals are very interesting due to their high THz-generation efficiencies, the possibility for high THz electric fields and a uniquely ultra-broadband nature. We can expect that the recent advances in the field will encourage new research in various disciplines with these systems. Taking into account also other studies performed by inorganic-based THz systems, we can soon expect some important studies in the field of biological and medical applications [219] because healthy and diseased tissue can be distinguished by THz waves due to the changes in water content, bacteria presence, structural and molecular changes. Therefore, THz spectroscopic and imaging analysis of small and macro biomolecules as well as cells and tissues is definitely one of the upcoming research opportunities for organic-crystals based systems. In addition, environmental applications are also an opportunity for laboratories operating with broadband THz systems such as based on organic NLO crystals. THz waves can penetrate even through rough environmental conditions like aerosol and heavily sooted flames [220]. Since THz radiation has the ability to gain vibrational information due to inter-and intra-molecular modes and rotational modes, the obtained information can be valuable for gaseous molecules quantification and identification especially in cases of air pollutants detection and atmospheric monitoring. One should also not forget other fields such as art conservation and archeology [221], where THz technology can play a key role in diagnostics before restoration of historical paintings and archaeological objects, and wood products industry [222]. Besides solids, which are already well explored in the THz region, more attention could also be given to liquids and aqueous solutions as well as gas spectroscopy.

Conclusions and Outlook
Organic NLO crystals have attracted a considerable interest within the fast-growing field of THz photonics. In this paper, we reviewed their unique properties that allow on the one hand for achieving state-of-the-art highly-intense THz fields and on the other hand can cover an ultrabroad THz frequency range of beyond 20 THz. These achievements give organic NLO sources an important advantage for the exciting fundamental and applied studies of THz light-matter interactions. We presented the basic material properties and requirements for optimized organic NLO crystals for THz photonics. These materials are composed of highly asymmetric electron-rich molecules with an ultra-fast response to the applied field, packed in a non-centrosymmetric crystalline lattice. Organic NLO crystals feature high long-term temporal as well as the required photochemical stability for high-optical-power applications. Few organic NLO crystals with record-high optical nonlinearities can be successfully grown and can be prepared with sufficient size, orientation, and optical quality for THz photonics applications. The development of novel materials is however being continued and is highly motivated by the wish of high THz power over a broad THz range without any gaps, which still present a challenge due to intrinsic phonon-mode resonances of NLO generation materials themselves. There is also a wish to increase the variety of organic NLO crystals optimized for certain available pump laser sources or compact high-power laser sources being currently developed, considering the optimal pump wavelength range and the desired THz frequency range. In parallel to the development of novel organic NLO crystals, there has been a significant progress in record high-power and ultra-broadband THz sources reported during the last few years. Organic crystalline THz sources became important in the development of the recently emerged field of nonlinear THz photonics, studying new fundamental aspects of matter enabled by intense THz waves. We also presented several examples of conventional THz spectroscopy and THz imaging based on organic crystalline THz sources, which show great potential for various scientific and industrial applications. We can therefore expect exciting further progress in research and applications of organic NLO crystals for THz photonics in the near future.