The Characterization and Application of Two Liquid Crystal Mixtures in the Low THz Region

: In the previous work, two new nematic liquid crystal (NLC) mixtures, E7-2 and S200-2, were produced by adding eight LC monomers to two commercial LCs S200 and E7, respectively. At λ = 589 nm, the birefringence ( ∆ n ) characteristics of the two LC nematic mixtures E7-2 ( ∆ n = 0.260) and S200-2 ( ∆ n = 0.298) are greater than those of the commercial LC E7 ( ∆ n = 0.224) and S200 ( ∆ n = 0.266). The properties ( T N-I , ε // , ∆ (cid:15) , K 11 , and K 33 ) of these four NLCs were measured. A double-layer metal loop arrays modulation structure based on metamaterial (MM) metal–dielectric–metal (MDM) was designed and fabricated for use in the THz frequency range. The results show that the LC mixtures E7-2 and S200-2 have greater modulation depth ( MD ) and less modulation insertion loss ( IL ) than E7 and S200 at THz frequencies. The results show that LC mixtures have signiﬁcant potential for designing active tunable LC-based devices in the THz and visible light range. and E7-2, S200-2 LC mixtures. The experimental results show that the metrics ∆ T , MD , and IL of these mixture-based LCs can be signiﬁcantly higher than those of commercial LCs in their working frequency band. This suggests that two high birefringence NLCs can also be used to fabricate highly e ﬃ cient and functional THz devices, such as phase shifters and switches. the dielectric anisotropy of mixture LC E7-2 is higher than that of commercial LC E7. The results indicate that mixture LC E7-2 can signiﬁcantly increase the maximum transmission spectra of double-layer metal loop arrays’ modulation. The relative phases of E7 and E7-2 have a small change ( ∆ ϕ < 90 ◦ ) at di ﬀ erent voltages.


Introduction
In recent years, terahertz (THz) technologies have been used for THz wireless communications [1], THz imaging [2], biomedical research [3], and environmental monitoring [4]. It is critical to control THz waves, with low-loss and high modulation, for applications using THz systems. Therefore, tunable THz functional devices have gradually become a major research topic [5]. There are many kinds of modulation mechanisms in the relevant materials that have been reported, such as graphene [6,7], VO 2 [8], electron gas [9], and liquid crystals (LCs) [10]. In 2014, the nematic liquid crystal (NLC)-based metamaterial was designed, fabricated, and used in THz tunable devices [11]. The transmittance tunability of the system is approximately 19% by changing the magnitude of the AC bias voltage from 0 to 300 V. In 2015, Kowerdziej et al. proposed a tunable LC THz MM device. As the applied voltage increased to 200 V, a maximum tuning value of 22% occurred at 0.82 THz. However, the modulation depth of the LC-based metamaterial is still very narrow and the bias voltage is too high.
Liquid crystals, in which the anisotropy can be changed by electric or magnetic means, have been successful at controlling and propagating electromagnetic waves. Wavelengths have included millimeter waves, visible light, and THz waves [12][13][14]. High birefringence values are especially important for devices working in the millimeter or THz ranges of electromagnetic waves. In recent years, 2 of 9 there has been increasing interest in characterizing LC materials as a tunable materials unit for THz applications, such as modulators [15], filters [16], phase shifters [17,18], switches [19], and absorbers [20]. High birefringence nematic liquid crystals have been broadly researched to satisfy the needs of the extensive variety of optical infrared range photonics devices [21]. Compared to NLCs for use in the visible region, these NLCs show weak birefringence in THz regions. Some papers demonstrated that the birefringence exhibited by LCs is relatively low in the THz band [22,23]. The relatively low birefringence can lead to performance limitations for THz electronically tunable devices. An effective way to solve these problems is to synthesize high birefringence LC compound materials. Recently, Vieweg et al. reported the birefringence of the LC mixture BL037, from ∆n = 0.16 at 0.5 THz to ∆n = 0.2 at 2.5 THz [24]. Wang et al. proposed a new LC mixture, called NJU-LDn-4, with a large birefringence of 0.306 from 0.4 to 1.6 THz [25]. Reuter et al. proposed the high birefringence LC mixtures 1852 (∆n = 0.33 at 0.5 THz) and 1825 (∆n = 0.38 at 0.5 THz) [26].
In the previous work, two new NLCs mixtures (E7-2 and S200-2) were proposed [27], which can be used to construct tunable THz devices due to tunable dielectric properties. In this paper, the more physical properties of these mixture LCs were measured and analyzed. Currently, to develop effective terahertz modulators and switches, electromagnetic metamaterials, an artificial composite material, are essential components. They consist of two-or three-dimensional arrays of tailored metallic unit cells [28]. In order to directly manipulate THz waves, we demonstrated THz LC modulators based on electromagnetic metamaterials and E7-2, S200-2 LC mixtures. The experimental results show that the metrics ∆T, MD, and IL of these mixture-based LCs can be significantly higher than those of commercial LCs in their working frequency band. This suggests that two high birefringence NLCs can also be used to fabricate highly efficient and functional THz devices, such as phase shifters and switches.

Materials and Methods
In order to obtain NLCs with high birefringence, eight LC monomers were designed and synthesized based on multi-ring rigid cores and polar terminal groups in Figure 1 [29]. We produced two nematic liquid crystals mixtures, E7-2 and S200-2, by adding eight LC monomers to two commercial NLCs: E7 (Jiangsu Hecheng Technology Co., Ltd. Nanjing, Jiangsu province, China.), and SLC103014-200 (S200) (Shijiangzhuang Chengzhi Yonghua Display Material Co., Ltd. Shijiangzhuang, Hebei province, China.) [27]. The NLC mixture of E7-2 and S200-2 were composed of 84 wt % of two commercial LC materials (E7 and S200), and 2 wt % of each new liquid crystal monomer. Dielectric properties of these two LC mixtures were measured using a frequency-selective surface in the sub-THz band. The dielectric anisotropy of the nematic E7-2 and S200-2 LC mixtures was 16.3% and 16.4% higher than that of the commercial E7 and S200 LCs, respectively [27].
Crystals 2020, 10, x FOR PEER REVIEW  2 of 10 THz applications, such as modulators [15], filters [16], phase shifters [17,18], switches [19], and absorbers [20]. High birefringence nematic liquid crystals have been broadly researched to satisfy the needs of the extensive variety of optical infrared range photonics devices [21]. Compared to NLCs for use in the visible region, these NLCs show weak birefringence in THz regions. Some papers demonstrated that the birefringence exhibited by LCs is relatively low in the THz band [22,23] [26].
In the previous work, two new NLCs mixtures (E7-2 and S200-2) were proposed [27], which can be used to construct tunable THz devices due to tunable dielectric properties. In this paper, the more physical properties of these mixture LCs were measured and analyzed. Currently, to develop effective terahertz modulators and switches, electromagnetic metamaterials, an artificial composite material, are essential components. They consist of two-or three-dimensional arrays of tailored metallic unit cells [28]. In order to directly manipulate THz waves, we demonstrated THz LC modulators based on electromagnetic metamaterials and E7-2, S200-2 LC mixtures. The experimental results show that the metrics ΔT, MD, and IL of these mixture-based LCs can be significantly higher than those of commercial LCs in their working frequency band. This suggests that two high birefringence NLCs can also be used to fabricate highly efficient and functional THz devices, such as phase shifters and switches.

Materials and Methods
In order to obtain NLCs with high birefringence, eight LC monomers were designed and synthesized based on multi-ring rigid cores and polar terminal groups in Figure 1 [29]. We produced two nematic liquid crystals mixtures, E7-2 and S200-2, by adding eight LC monomers to two commercial NLCs: E7 (Jiangsu Hecheng Technology Co., Ltd. Nanjing, Jiangsu province, China.), and SLC103014-200 (S200) (Shijiangzhuang Chengzhi Yonghua Display Material Co., Ltd. Shijiangzhuang, Hebei province, China.) [27]. The NLC mixture of E7-2 and S200-2 were composed of 84 wt % of two commercial LC materials (E7 and S200), and 2 wt % of each new liquid crystal monomer. Dielectric properties of these two LC mixtures were measured using a frequency-selective surface in the sub-THz band. The dielectric anisotropy of the nematic E7-2 and S200-2 LC mixtures was 16.3% and 16.4% higher than that of the commercial E7 and S200 LCs, respectively [27]. The physical properties (no, Δn, TN-I, ε//, Δɛ, K11, and K33) of the eutectic mixtures (E7-2 and S200-2) were measured using two commercial LC materials (E7 and S200) as benchmarks for comparison. The measurement results are shown in Table 1. The birefringence of NLC mixtures have been measured using the abbe refractometer (NAR-4T) at λ = 589 nm. The measured birefringence of E7 is ∆n = 0.224 and E7-2 is ∆n = 0.260. The birefringence measured of S200 is ∆n = 0.224. The S200-2 is ∆n The physical properties (n o , ∆n, T N-I , ε // , ∆ , K 11 , and K 33 ) of the eutectic mixtures (E7-2 and S200-2) were measured using two commercial LC materials (E7 and S200) as benchmarks for comparison. The measurement results are shown in Table 1. The birefringence of NLC mixtures have been measured using the abbe refractometer (NAR-4T) at λ = 589 nm. The measured birefringence of E7 is ∆n = 0.224 and E7-2 is ∆n = 0.260. The birefringence measured of S200 is ∆n = 0.224. The S200-2 is ∆n = 0.260. The birefringence values of the LC mixtures E7-2 and S200-2 are 16.1% and 12.0% higher than the commercial LC materials E7 and S200, respectively. It means that eutectic mixtures with more benzene rings, NCS terminal, and triple bridge bond exhibit larger birefringence values when compared to commercial LCs. The clearing point (T N-I ) was measured using differential scanning calorimetry. The values of T N-I for E7 and E7-2 are 58.2 • C and 72.2 • C, respectively, and the values of T N-I for S200 and S200-2 are 101.1 • C and 114.9 • C, respectively. The T N-I of the eutectic mixtures is much larger than it is for commercial NLCs. This means eutectic mixtures have wide nematic phase temperature ranges. The dielectric properties (ε // and ∆ ), and the splay (K 11 ) and bend (K 33 ) elastic constants of these NLC mixtures, were measured using Instec's automatic liquid crystal test equipment (ITO cell thickness: 11 µm). Here, ε // is the maximum dielectric constant of the LC cell and ∆ (∆ = ε // − ε ⊥ ) is the LC's dielectric anisotropy at 1 kHz. The values of ε // , ∆ , K 11 and K 33 in eutectic mixtures (E7-2 and S200-2) are smaller than that in commercial LCs E7 and S200 at 1 kHz. Table 1. Physical properties of nematic liquid crystals (NLCs) E7, E7-2, and S200 and S200-2.

Applications of Devices for Terahertz Continuous Wave Modulation
The application frequency range of interest for a LC for a tunable metamaterial-LC-metamaterial (MLM) structure is about from 290 GHz to 320 GHz. The unit cell structure of the LC-based tunable THz MLM is described in Figure 2a. The double-layered metal MM structure is a periodic surface with identical two-dimensional arrays of sub-wavelength circular air loops, using the ultraviolet lithography process shown in Figure 2b. This structure was embedded in the LC layer contained within an MLM structure, as is shown in Figure 2. It is composed of two parallel fused silica substrates that are separated by a 45 µm spacer. A layer of arrayed metal metamaterial unit cells was deposited on the inner surface of the quartz substrate. Polyimide (Pi) is an organic insulating material with a low dielectric constant. It has been widely used to align LC molecules. Both electrodes are spin-coated with sulfonic Pi as the alignment layers and are photo-aligned. A cell is infiltrated with the NLCs. The optical axis of the LC is initially parallel to the connecting wires. These LC materials are composed of anisotropic molecules whose direction may change due to an external electrical field between the copper patterned layer and the metal ground.
The cell is driven by applying a 1 kHz square wave current signal through copper tapes connecting the electrodes to the electric power source. As the applied voltage increased, the optical axes of the LC molecules gradually change along the electric field vector as shown in Figure 2c. Without an electric field, the long axes of the LC molecules lie parallel to the surface of each quartz plate. The dielectric constant of the LC layer is ε ⊥ . When a saturating bias voltage is applied to the LC layer, the long axes of the LC molecules will be parallel to the electric field. The dielectric constant of the LC layer is ε // . The structure demonstrates a pronounced radiant mode that can be used to develop THz emitters. The cell is driven by applying a 1 kHz square wave current signal through copper tapes connecting the electrodes to the electric power source. As the applied voltage increased, the optical axes of the LC molecules gradually change along the electric field vector as shown in Figure 2c. Without an electric field, the long axes of the LC molecules lie parallel to the surface of each quartz plate. The dielectric constant of the LC layer is ε⊥. When a saturating bias voltage is applied to the LC layer, the long axes of the LC molecules will be parallel to the electric field. The dielectric constant of the LC layer is ε//. The structure demonstrates a pronounced radiant mode that can be used to develop THz emitters.
The transmittance spectrum can be changed to different NLC dielectric constants. The transmission characteristics and the LC parameters are calculated and analyzed for THz electromagnetic waves in the frequency range from 220 to 330 GHz. The EM performance of the LCbased MLM structure was analyzed and optimized using the finite difference time domain (FDTD) method. In the simulation, assuming that the unit cell is located in a periodic environment, and illuminated by a linear polarization and normal incidence plane wave. Period boundary conditions are used in the x-axis and y-axis directions, and the boundary condition in the z-axis direction is set to open. So, the EM properties of the THz MDM structure were obtained by FDTD. The simulated transmittance spectra of the modulator for different dielectric constants is shown in Figure 3. The transmittance spectrum can be changed to different NLC dielectric constants. The transmission characteristics and the LC parameters are calculated and analyzed for THz electromagnetic waves in the frequency range from 220 to 330 GHz. The EM performance of the LC-based MLM structure was analyzed and optimized using the finite difference time domain (FDTD) method. In the simulation, assuming that the unit cell is located in a periodic environment, and illuminated by a linear polarization and normal incidence plane wave. Period boundary conditions are used in the x-axis and y-axis directions, and the boundary condition in the z-axis direction is set to open. So, the EM properties of the THz MDM structure were obtained by FDTD. The simulated transmittance spectra of the modulator for different dielectric constants is shown in Figure 3. The cell is driven by applying a 1 kHz square wave current signal through copper tapes connecting the electrodes to the electric power source. As the applied voltage increased, the optical axes of the LC molecules gradually change along the electric field vector as shown in Figure 2c. Without an electric field, the long axes of the LC molecules lie parallel to the surface of each quartz plate. The dielectric constant of the LC layer is ε⊥. When a saturating bias voltage is applied to the LC layer, the long axes of the LC molecules will be parallel to the electric field. The dielectric constant of the LC layer is ε//. The structure demonstrates a pronounced radiant mode that can be used to develop THz emitters.
The transmittance spectrum can be changed to different NLC dielectric constants. The transmission characteristics and the LC parameters are calculated and analyzed for THz electromagnetic waves in the frequency range from 220 to 330 GHz. The EM performance of the LCbased MLM structure was analyzed and optimized using the finite difference time domain (FDTD) method. In the simulation, assuming that the unit cell is located in a periodic environment, and illuminated by a linear polarization and normal incidence plane wave. Period boundary conditions are used in the x-axis and y-axis directions, and the boundary condition in the z-axis direction is set to open. So, the EM properties of the THz MDM structure were obtained by FDTD. The simulated transmittance spectra of the modulator for different dielectric constants is shown in Figure 3. Numerical results for resonant frequencies were obtained from the numerical FDTD, and fitted to the measured results. The conductivity of the copper layer was set to 5.8 × 10 7 S/m. The dielectric constants are set from 2.54 to 3.2. The transmittance spectrum can be changed to different LC dielectric constants. As the dielectric constants increase, the lower transmittance spectrum will shift toward a Crystals 2020, 10, 99 5 of 9 low resonance frequency. This reveals that the LC-based MLM structure can change from its maximum to its minimum transmittance value.
The modulation depth (MD), insertion loss (IL), and ∆T are used to describe the performance of the LC-based modulator device. The MD, IL, and ∆T of the modulator structure at a frequency point can be calculated using: where the parameters T max and T min are maximum and minimum transmittances of the structure at a certain frequency point. The variables ϕ max and ϕ min are the maximum and minimum phases of the structure at a certain frequency point, respectively.

Experimental Results and Discussions
In order to experimentally verify the performance of the proposed device, the MLM structure was fabricated with a two-layer metal MM array using the ultraviolet lithography process. The manufactured device and the microscopic image of the unit cell are shown in Figures 3b and 4a. The EM characteristics of the MLM structure were measured at room temperature using a vector network analyzer system and the transmission spectra were obtained at different bias voltages as shown in Figure 4c. The effective dielectric constant of the LCs molecules depends on the bias voltage. The Agilent N5224A vector network analyzer and the VDI mm-wave module extender (from 220 GHz to 325 GHz) were used to measure the transmittance spectra of the LC samples. Two pairs of lenses were used to focus the THz wave. The maximum voltage (20 V) with a 1 kHz square wave was applied to the patterned copper layer and the copper ground plane in the experiment. Numerical results for resonant frequencies were obtained from the numerical FDTD, and fitted to the measured results. The conductivity of the copper layer was set to 5.8 × 10 7 S/m. The dielectric constants are set from 2.54 to 3.2. The transmittance spectrum can be changed to different LC dielectric constants. As the dielectric constants increase, the lower transmittance spectrum will shift toward a low resonance frequency. This reveals that the LC-based MLM structure can change from its maximum to its minimum transmittance value.
The modulation depth (MD), insertion loss (IL), and ΔT are used to describe the performance of the LC-based modulator device. The MD, IL, and ΔT of the modulator structure at a frequency point can be calculated using: where the parameters Tmax and Tmin are maximum and minimum transmittances of the structure at a certain frequency point. The variables φmax and φmin are the maximum and minimum phases of the structure at a certain frequency point, respectively.

Experimental Results and Discussions
In order to experimentally verify the performance of the proposed device, the MLM structure was fabricated with a two-layer metal MM array using the ultraviolet lithography process. The manufactured device and the microscopic image of the unit cell are shown in Figure 3b and Figure  4a. The EM characteristics of the MLM structure were measured at room temperature using a vector network analyzer system and the transmission spectra were obtained at different bias voltages as shown in Figure 4c. The effective dielectric constant of the LCs molecules depends on the bias voltage. The Agilent N5224A vector network analyzer and the VDI mm-wave module extender (from 220 GHz to 325 GHz) were used to measure the transmittance spectra of the LC samples. Two pairs of lenses were used to focus the THz wave. The maximum voltage (20 V) with a 1 kHz square wave was applied to the patterned copper layer and the copper ground plane in the experiment.  For applied voltages increasing from 0 V to 20 V, the transmission characteristics of the transmissive MLM structure are shown in Figures 4 and 5, respectively. When the bias voltage is varied from 0 to 4 V, samples have rapid spectral transmittance changes. The transmittance of the samples slowly rises and eventually stabilizes when the bias voltage is higher than 4 V. For a zero applied bias voltage, the black curve (0 V) represents the maximum transmission spectrum. When the applied voltage increased to its maximum, the pink curve (20 V) is the minimum transmission spectra. The ∆T and MD values increase as the external applied bias is increased at a frequency point.
transmissive MLM structure are shown in Figures 4 and 5, respectively. When the bias voltage is varied from 0 to 4 V, samples have rapid spectral transmittance changes. The transmittance of the samples slowly rises and eventually stabilizes when the bias voltage is higher than 4 V. For a zero applied bias voltage, the black curve (0 V) represents the maximum transmission spectrum. When the applied voltage increased to its maximum, the pink curve (20 V) is the minimum transmission spectra. The ΔT and MD values increase as the external applied bias is increased at a frequency point.
To further investigate the performance of the demonstrated sample at multiple frequency points, from Figure 5, we extracted the MD, IL, ΔT, and Δφ for specific frequency points. Table 2 lists the values for from frequencies ranging from 296 GHz to 312 GHz.   Figure 5a, the resonance frequency point of the E7-2 shift toward 304.8 GHz in Figure 5b, and the ΔT and MD values of the LC E7-2 are larger than that of the commercial LC E7 between 296 GHz and 305 GHz. The IL of the LC E7-2 sample is less than that of the commercial LC E7 from 296 GHz to 312 GHz. For example, the values of ΔT and MD of the LC E7-2 are 0.56 and 0.76 at 300 GHz, respectively, and ΔT and MD of the LC E7 are 0.47 and 0.69, respectively. The value of IL (1.51 dB) of the LC E7-2 sample is less than the IL (1.54 dB) of the commercial LC E7. These comparisons show that the dielectric anisotropy of mixture LC E7-2 is higher than that of commercial LC E7. The results indicate that mixture LC E7-2 can significantly To further investigate the performance of the demonstrated sample at multiple frequency points, from Figure 5, we extracted the MD, IL, ∆T, and ∆ϕ for specific frequency points. Table 2 lists the values for from frequencies ranging from 296 GHz to 312 GHz. Compared to the resonance frequency point (311.6 GHz) of the low transmittance spectrum of the E7 LC characteristic shown in Figure 5a, the resonance frequency point of the E7-2 shift toward 304.8 GHz in Figure 5b, and the ∆T and MD values of the LC E7-2 are larger than that of the commercial LC E7 between 296 GHz and 305 GHz. The IL of the LC E7-2 sample is less than that of the commercial LC E7 from 296 GHz to 312 GHz. For example, the values of ∆T and MD of the LC E7-2 are 0.56 and 0.76 at 300 GHz, respectively, and ∆T and MD of the LC E7 are 0.47 and 0.69, respectively. The value of IL (1.51 dB) of the LC E7-2 sample is less than the IL (1.54 dB) of the commercial LC E7. These comparisons show that the dielectric anisotropy of mixture LC E7-2 is higher than that of commercial LC E7.
The results indicate that mixture LC E7-2 can significantly increase the maximum transmission spectra of double-layer metal loop arrays' modulation. The relative phases of E7 and E7-2 have a small change (∆ϕ < 90 • ) at different voltages. Figure 6 shows the transmission spectra results of NLC S200 and S200-2. These transmission spectra curves were measured for different voltages. The values of MD, IL, ∆T, and ∆ϕ are listed in Table 3.
Crystals 2020, 10, x FOR PEER REVIEW 7 of 10 increase the maximum transmission spectra of double-layer metal loop arrays' modulation. The relative phases of E7 and E7-2 have a small change (Δφ < 90°) at different voltages. Figure 6 shows the transmission spectra results of NLC S200 and S200-2. These transmission spectra curves were measured for different voltages. The values of MD, IL, ΔT, and Δφ are listed in Table 3. Compared to the resonance frequency point (308.8 GHz) of the low transmittance spectrum (the pink curve) of the S200 in Figure 6a, the resonant frequency point of the S200-2 shifts toward 302.7 GHz in Figure 6b. The ΔT and MD of the LC S200-2 mixture are larger than those of the commercial LC S200 from 296 GHz to 305 GHz. The IL of the LC S200-2 sample is less than that of the commercial LC S200. The phase of the S200 and S200-2 have small changes (Δφ < 90°) at different voltages in the working frequency band.
For example, the ΔT (0.6) and MD (0.83) of S200-2 are larger than the corresponding ΔT (0.49) and MD (0.75) values of the commercial LC S200 at 300 GHz. The IL (1.36 dB) of the LC S200-2 sample is less than the IL (2.2 dB) of the commercial LC S200 at 300 GHz. According to numerical FDTD simulation data, the values show that the dielectric anisotropy of the LC S200-2 mixture is higher than that of commercial LC S200.   Compared to the resonance frequency point (308.8 GHz) of the low transmittance spectrum (the pink curve) of the S200 in Figure 6a, the resonant frequency point of the S200-2 shifts toward 302.7 GHz in Figure 6b. The ∆T and MD of the LC S200-2 mixture are larger than those of the commercial LC S200 from 296 GHz to 305 GHz. The IL of the LC S200-2 sample is less than that of the commercial LC S200. The phase of the S200 and S200-2 have small changes (∆ϕ < 90 • ) at different voltages in the working frequency band.
For example, the ∆T (0.6) and MD (0.83) of S200-2 are larger than the corresponding ∆T (0.49) and MD (0.75) values of the commercial LC S200 at 300 GHz. The IL (1.36 dB) of the LC S200-2 sample is less than the IL (2.2 dB) of the commercial LC S200 at 300 GHz. According to numerical FDTD simulation data, the values show that the dielectric anisotropy of the LC S200-2 mixture is higher than that of commercial LC S200.
These results indicate that the LC mixtures E7-2 and S200-2 can significantly increase the maximum transmission spectra of double-layer metal loop arrays' modulation capabilities. The dielectric anisotropy of LC mixtures can be significantly increased by adding eight LC monomers. The birefringence these NLC materials exhibit high into the visible range would present comparatively high dielectric anisotropy in THz range as well.

Conclusions
In this paper, two new NLC mixtures (the LC S200-2 and the LC E7-2) were proposed and the birefringence of the nematic LC E7-2 and S200-2 mixtures is higher by 16.1% and 12.0% than the birefringence of the commercial LC products E7 and S200 at visible light, respectively. The eutectic mixtures (E7-2 and S200-2) have wide nematic phase temperature ranges. The values of ε // , ∆ , K 11 , and K 33 in eutectic mixtures are smaller than those of the commercial LCs E7 and S200 at 1 kHz. Based on the proposed LC and metamaterials, a THz MLM structure was fabricated to control terahertz waves using the ultraviolet lithography process. The experimental results are shown that LC mixtures E7-2 and S200-2 have larger ∆T, MD, and less IL for modulation than those E7 and S200 at THz frequencies for normal incidence of THz waves. The relative phases of LC mixtures have small change (∆ϕ < 90 • ) at different voltages. It also means that the ∆ of the commercial liquid crystals (E7 and S200) is improved by adding eight LC monomers. Our work may provide an effective solution for the design of new high birefringence LC materials in the THz range.