‘Quasi-Tricritical’ and Glassy Dielectric Properties of a Nematic Liquid Crystalline Material

: Results of dielectric studies in the nematic and isotropic liquid phases of pentylcyanobiphenyl (5CB), a model rod-like liquid crystalline material, are shown. They are based on the discussion of the evolution of dielectric constant ( ε ), its changes under the strong electric ﬁeld (nonlinear dielectric e ﬀ ect, NDE), and ﬁnally, the primary relaxation time. It is shown that changes in ε ( T ) and NDE are entirely dominated by the impact of pretransitional ﬂuctuations (pre-nematic and pre-isotropic, respectively) which are associated with the weakly discontinuous character of the isotropic–nematic phase transition. This inﬂuence also extends for the low-frequency, ionic species dominated, region. Notable, that the derivative-based and distortions sensitive analysis revealed the tricritical nature of the I-N transition. Although the glass transition in 5CB occurs in the deeply supercooled state at T g ≈ − 68 ◦ C, the dynamics (changes of the primary relaxation time) follow a previtreous pattern both in the isotropic and in the nematic phase. Finally, the discussion of the ’ molecular ’ vs. ‘ quasi-critical ’ characterizations of the isotropic and nematic phases is presented. It shows the evident prevalence of the ‘quasi-critical-picture’, which o ﬀ ers the consistent temperature parameterization in the total tested temperature range.


Introduction
The name 'liquid crystals' may appear self-contradictory, but it describes the real nature of unique states locating between the isotropic liquid and the solid crystalline state. Liquid crystals (LCs) are built from molecules creating fluids in which elements of the long-range order of crystals with the mobility of liquids are combined. Liquid crystals are crucial for challenging chapters of physics, chemistry, and modern biology [1][2][3][4][5][6][7]. They are also omnipotent for a variety of applications, mainly displays used from smartphones to computers and TV sets [8,9]. As the most classical liquid crystalline material, one can indicate pentylcyanobiphenyl (5CB)-the rod-like compound with the permanent dipole moment (µ ≈ 5 D) approximately parallel to the long molecular axis [3]. It exhibits only the nematic mesophase, with the uniaxial orientational ordering while the translational disorder remains. The essential reference point for properties of 5CB are isotropic-nematic (I-N) and the nematic-solid (I-S) phase transitions. [3,4] Notable, that one of the most significant theoretical concepts of the Physics of Liquid Crystals and the Soft Matter Physics, Landau-de Gennes (LdG) model, was created to comment strong and surprisingly similar pretransitional changes of the Cotton-Mouton effect (CME, ∆n/B 2 ), Kerr effect (KE, ∆n/E 2 ) and the intensity of the scattered light (I L ) in the isotropic liquid phase [1,3,10]: where ∆n = n || − n ⊥ is the difference of refractive indexes for samples with molecules oriented in a parallel and perpendicular way in respect to the long molecular axis; E is the intensity of the electric field and B demotes the magnetic induction. T > T C = T * + ∆T * and T C = T I−N is for the I-N 'clearing' temperature, T * denotes the hypothetical, extrapolated, continuous phase transition and ∆T * is the measure of the discontinuity of the I-N transition. Equation (1) can be linked to the temperature evolution of compressibility described by the power/critical exponent γ = 1. For nematogenic liquid crystalline materials, particularly important for applications, such properties as dielectric constant (ε), refractive index (n) and electric conductivity (σ), primary relaxation time (τ) or viscosity (η) are particularly important [2][3][4]8]. There are five decades of studies indicating for these properties no pretransitional anomalies on approaching the I-N transition [2][3][4]8,[11][12][13]: n(T), ε(T) = A n,ε + B n,,ε T (2) η(T), σ(T), τ(T) = τ σ,η,τ 0 where Equation (3) is for the Arrhenius behavior and E σ,η,τ = const denotes the activation energy. Essentially important, regarding both fundamentals and applications, are changes of components of the refractive index (n ⊥ (T), n || (T)) and dielectric constant (ε ⊥ (T),ε || (T)) in the nematic phase [2][3][4]8]. Their evolutions are most often analyzed in terms of the Haller relation [3,4,[14][15][16][17][18]: and where the empirical exponent β = 0.1 − 0.22. It is expected that in the nematic phase below the clearing temperature T C = T I−N : ∆ε(T) = ε || (T) − ε ⊥ (T) = ∆ε = const and ∆n(T) = n || (T) − n ⊥ (T) = ∆n = const, i.e., constant values for the completely ordered sample are reached. They are assumed to reflect the molecular anisotropy of LC molecules [3]. Generally, no pretransitional behavior near the nematic-solid phase transition is expected. [1][2][3][4][5][6][7][8] Notwithstanding, there is a limited evidence of such behavior, in 2-cyano-4-pentylbiphenyl 4-(trans-4-pentylcyclohexyl) benzoate (the LC compound with the transverse location of the permanent dipole moment), [19] or 5CB and BaTiO 3 nanoparticles composites [20,21], for instance.
However, there are blots in the simple and broadly applied 'research landscape' sketched above [1][2][3][4]8]. First, both ∆n(T) and ∆ε(T) changes are considered as basic metrics of the order parameter in the Physics of Liquid Crystals and the Phase Transitions, including the LdG model [1][2][3]22]. The order parameter constitutes here the basic parameter, and its evolution has to be described by the universal, well-defined and constant value of the exponent β [22], in clear disagreement with the Haller relation [14]. Second, already in 1981 Bradshaw and Raynes [23] and in 1983 Thoen and Menu [24] presented the evidence for a notable pretransitional anomaly in the isotropic phase of n-cyanobiphenyls, including 5CB for dielectric constant. Starting from 1996, the author of this report tested the behavior of dielectric constant in the isotropic phase, first on approaching the I-N transition and subsequently also for the Isotropic-Smectic A (I-SmA), Isotropic-Smectic E (I-SmE) and the Isotropic-Chiral Nematic (I-N*) phase transitions. In each mentioned case, studies carried out even up to T C + 100 K validated the following pretransitional behavior [25][26][27]: where T > T C = T * + ∆T * , (ε * , T * ) are for the hypothetical, extrapolated, continuous phase transition; the power exponent α = 1/2 is linked to the specific heat pretransitional/critical anomaly c P (T). Following tests in the homogeneous phase of binary mixtures of limited miscibility, it was suggested Notable, that the pretransitional anomaly in the isotropic liquid phase described by Equation (6) is relatively strong and easy for detecting, but requiring high-resolution experiments which takes into account the general rule for 'critical phenomena tests': the temperature distance between subsequent tested temperatures should strongly decrease when approaching the singular temperature (T * in the given case).
In the last two decades, such in-deep studies taking into account Equation (6) have been carried out also by 'other' authors, extending the evidence confirming the appearance of the pretransitional effect in the isotropic liquid phase and the validity also for LC molecule beyond the basic rod-like pattern [28][29][30]. Notably, that the lack of ε(T) pretransitional 'anomaly' is also possible (see Equation (2)) but only for systems where the permanent dipole moment is perpendicular to the long molecular axis [19,31]. It is coupled to the exponent α = 0.
One should stress that direct experimental studies of the heat capacity/specific heat yielded a strong record for pretransitional anomalies associated with the exponent α ≈ 1/2 on both sides of the I-N transition, independently from the position or even absence of the permanent dipole moment [3,22]. In classical monographs regarding the Physics of Liquid Crystals the surrounding of the I-N transition in rod-like LC materials is indicated as the model case for which the basic mean-field theory (MFT) is valid [1][2][3][4]. From the experimental point of view power/critical exponents constitute the key check-point in testing theoretical predictions. For the basic MFT model [22]: α = 0, (T > T C = T * ) and α = 0 (T < T C = T * ), β = 1/2 and γ = 1, (T > T C = T * ), and γ = 1 (T < T C = T * ) for the compressibility related exponent.
Concluding, rod-like liquid crystalline materials with the nematic phase are considered as the experimental model-system for the fundamental insights and society-important applications [1][2][3][4][5][6][7][8][9][10][11]. Notwithstanding, some of their essential properties remain surprisingly puzzling, as indicated above. One can also add questions regarding the range and degree of the impact of phase transitions on basic features, or the still poorly tested relationship between pretransitional fluctuations and ionic dopants, which dominates dielectric properties in the low-frequency domain [11][12][13].
The target of this report is the experimental evidence and explanation of the cognitive gap indicated above. Results are based on broadband dielectric spectroscopy, and nonlinear dielectric effect (changes of dielectric constant under the strong electric field) studies carried out in isotropic and nematic phases of rod-like pentylcyanobiphenyl (5CB). Unique features, hidden so far, have been revealed due to the high-resolution measurements in the broad temperature range, as well as the innovative distortions-sensitive analytic methodology.

Materials and Methods
Experimental studies were conducted using the broadband dielectric spectrometer (Novocontrol) with the Quattro unit for temperature control. It enables complex impedance measurements up to f ∼ 3 GHz . The analysis explored the complex dielectric permittivity ε * = ε + iε . They were carried out in pentylcyanobiphenyl (5CB), with the mesomorphism [2,3]: Solid Crystal-(T m ≈ 280 K)-Nematic-(T C ≈ 308 K)-Isotropic liquid. 5CB molecules are terminated by the permanent dipole moment (µ ≈ 5 D), approximately parallel to the long molecular axis [3]. Samples of 5CB were deeply degassed immediately before measurements, basing on the sequence 'freeze and pump'. Such processing also significantly reduced the electric conductivity of tested samples. The analysis of experimental data explored innovative derivative-based tools, developed by the author [32][33][34]. Measurements of dielectric constant in the nematic phase of 5CB were carried out using the flat-parallel capacitors made from the gold-coated copper, with the gap d = 0.2 mm , and diameter 2r = 20 mm. The applied magnetic field B ≈ 1.2 T.
Dielectric studies (broadband dielectric spectroscopy BDS, nonlinear dielectric effect NDE) in the isotropic phase and for non-oriented samples in the nematic phase were carried out using the capacitor made from Invar. They also offer the insight into dynamics focused on the primary (structural, alpha) relaxation process. Details regarding this issue are given in the Appendix A.
Nonlinear dielectric effect (NDE) describes changes of dielectric constant due to the action of the strong electric field E: [35][36][37]. Such studies, associated with measurement of a tiny shift of the electric capacitance, ∆C(E)/ C ∼ 10 −6 , were performed using the dual-field resonant method, with the radio-frequency weak measuring field U peak−peak ≈ 1 V and the DC strong electric field U DC ≈ 300 V − 1000 V applied in the form of ca. ∆t ≈ 1 ms lasting pulses. The latter makes the impact of ionic-dopants, which can cause the parasitic heating effect, negligible. This report is based on the results of in-deep re-analyzed earlier findings of the author [21,33,38,39] and new measurements. The latter is related mainly to NDE studies and results presented in Figure 1 as well as in the Appendix A. Figure 1. The temperature evolution of dielectric constant (ε = ε ( f = 10 kHz)) in the isotropic liquid and nematic mesophases of pentylcyanobiphenyl (5CB). In the nematic phase, components related to the long (ε || ) and short-axis (ε ⊥ ) of the rod-like molecule are shown. The evolution of the mean dielectric permittivity in the nematic phase is calculated via the relation Pretransitional anomalies in the isotropic liquid and for the mean in the nematic region are described by Equations (6)- (8). The inset shows the 'focused' pretransitional effect in the isotropic liquid phase. Experimental data in blue are for the non-oriented sample: note supercooling well below the melting temperature T m . Figure 1 presents the behavior of dielectric constant in the isotropic and nematic phases of liquid crystalline 5CB. The inset in Figure 1 shows the clear manifestation of the pretransitional effect in the isotropic liquid, associated with the crossover: dε/ dT < 0 → dε/dT > 0 . The discontinuous I-N transition exhibits features of the weakly discontinuous phase transition, i.e., it is associated with the influence of the hidden continuous phase transition at temperature T * , located only ca. ∆T * ≈ 1.5 K below the clearing temperature T C = T I−N . Consequently, on cooling towards T C in the isotropic liquid phase, pretransitional fluctuations of the subsequent nematic phase appear. Their size and lifetime are described by Equation (1), and 5CB molecules within fluctuations are ordered in the way characteristic for the nematic phase. The latter means the equivalence of → n and − → n directors, i.e., vectors showing the preferred orientational arrangement of the long molecular axis of molecules [1][2][3]. This feature leads to the random orientation of permanent dipole moments within prenematic fluctuations, and then the cancellation of the related 'dipolar contribution' to dielectric constant characterizing fluctuations. Consequently, dielectric constant related to fluctuation is qualitatively smaller than for the 'chaotic, fludlike' surrounding: ε f luct. << ε surrounding . This leads to the appearance of the crossover temperature, dε(T = T cross. )/dT = 0 for T → T I−N , for which contributions to the total dielectric constant from prenematic fluctuations and the isotropic liquid surrounding are equal.

Results and Discussion
The behavior of dielectric constant in the isotropic phase is described via the relation [25][26][27]: where T > T I−N , T * = T I−N − ∆T * . Values of parameters obtained: ε * = 11.197, T * = 305.16 K, α = 0.51 ± 0.02, a = −0.00226 and A = 0.1236. This result is shown as the blue curve in Figure 1. Experimental data presented in Figure 1 make it possible to calculate the mean value of the dielectric constant in the nematic phase, due to the symmetry of 5CB molecule defined as It also exhibits a strong pretransitional anomaly, described by the relation [40,41]: Values of the fitted parameter ε * * = 11.878, T * * = 306.27 K, The experimental results presented in Figure 1 enables the analysis of the order parameter behavior in the nematic phase. According to the Physics of Critical Phenomena and Phase Transitions [22], temperature changes of the order parameter should follow the dependence [40,41]: This relation contains three adjustable parameters and has to be used for experimental data that terminate at T = T I−N , i.e., well above the singular temperature T * * . Thus, the domain of the most significant pretransitional changes of ∆ε(T) is excluded from the analysis, what has to increase the fitting uncertainty. However, one can consider the derivative of ∆ε(T) experimental data, reducing the number of adjustable parameters: The application of such analysis reduces the fitting solely to one parameter (T * * ). Results are shown in Figure 2.
It is particularly worth stressing that results presented in Figures 1 and 2 prove the dominated impact of the pretransitional behavior for the whole nematic and the isotropic liquid phases. This behavior is clearly linked to the impact of collective pretransitional fluctuations. Figure 1 also presents the temperature evolution of dielectric constant in the non-oriented supercooled nematic, for which supercooling well below T m was possible. In the given case, the nematic phase was reached on cooling from the isotropic liquid and the additional DC electric field (U=40 V) was applied only just at the phase transition temperature T I−N . This causes that in the nematic phase, changes of ε(T) follows the pattern close to ε || (T) behavior. Generally, it is expected that on cooling ε || (T) and ε ⊥ (T) approach horizontal asymptotes, which can estimate the 'molecular anisotropy' ∆ε = ε || − ε ⊥ = const for the perfectly ordered sample. For the non-oriented 5CB sample, the supercooling below T m was possible, and the mentioned tendency for ε || (T) qualitatively changed: the notable decrease of the value of ε N (T) is visible. In the opinion of the author, the link of this phenomenon to the growing impact of topological defects distorting the uniaxial arrangement can be indicated [42]. Such a possibility can show experimental studies and related modeling for 5CB + nanoparticles composites [20,21].
From the application point of view, the impact of ionic impurities on dielectric properties is meaningful, becoming significant for the low frequencies domain. The link between prenematic fluctuations and the impact of ionic contaminants is still poorly evidenced. [43][44][45] Results presented in Figures 1 and 2 show the dominated influence of pretransitional fluctuations on the dielectric constant in the broad surrounding of the I-N 'clearing' temperature. The question arises if such impact also extends into the low-frequency domain. Results presented in Figure 3 can respond to this question. In the plot as the reference, the value of dielectric permittivity in the clear static domain ε = ε ( f = 100 kHz), where the influence of the ionic contribution can be neglected, is assumed. It is visible that the presence of prenematic fluctuations leads to the linear behavior of the ionic-related contribution. This means that in the isotropic phase of 5CB the critical-like Equation (7) remains valid also for frequencies below the static domain, and the impact of ionic impurities appears only as an additional linear term. Such simple behavior can be linked to the impact of appearing prenematic fluctuations, what further suggest its reduction for very low frequencies. Figure 2. Results of the derivative-based analysis of the order parameter in the nematic phase of 5CB (see Equation (10). They strongly indicate the tricritical value of the order parameter exponent β ≈ 1/4, well beyond the mean-field value exponent β = 1/2 [3,22]. Figure 3. The experimental temperature dependence of the low-frequency, ionic-related, behavior of the real part of dielectric permittivity in the isotropic phase of 5CB. As the 'non-ionic' reference dielectric constant for the frequency f = 100 kHz was assumed [39].
Broadband dielectric spectroscopy studies also deliver essential messages regarding dynamic properties, amongst which the primary relaxation time plays the key role [46] (see the Appendix A). Figure 4 presents the re-analyzed experimental data, first introduced by the author in ref. [38], where the supercooled state of 5CB down to the glass temperature T g was reached. Such challenging target, in the non-oriented sample, required the special preparation of the sample, namely the deep cleaning leading to the very low electric conductivity (reached by the DC electric field action), deep degassing, and avoiding any 'dust contaminations" due to operation only is a special lab hood. Particularly noteworthy was the design of the measurement capacitor (made from invar and gold-coated), avoiding any 'edges'. To reduce the period of measurements in the non-equilibrium metastable state, tests were focused exclusively on the detection of the relaxation time. Experimental data are shown in the inset in Figure 4. The central part of the figure shows the proof for the preference of critical-like dynamics, however associated with the previtreous behavior and the glass transition.  (14)) focused on the validity of the critical-like behavior (Equation (13)), indicated by linear domains. The linear regression fits yield parameters for Equation (13) given in the plot. They are used for portraying via Equation (13) τ(T) experimental data presented in the inset. The magnitude H a (T) = H a (T)/R is for the normalized apparent activation enthalpy. MCT means 'mode-coupling theory' [47,48], and 'DSM' is for the 'Dynamical Scaling Model' [49]. These models are developed for glass-forming liquids in the high-and low-temperature dynamical domains. The glass temperature is estimated via the condition τ T g = 100 s [48].
The analysis presented in Figure 4 recalls patterns of dynamics developed for glass-forming systems. The preferred form of the description is revealed due to the application of the distortions-sensitive and derivative-based analysis via the transformation of experimental data introduced in ref. [32]: Primary relaxation time τ(T) →→-Activation Enthalpy H a (T)/R = d ln τ(T)/d(1/T). The latter dependence results directly from the general super-Arrhenius (SA) relation, with the temperature-dependent activation energy E a (T) [50][51][52]: The direct application of Equation (11) for describing experimental data is not possible because of the unknown general form of E a (T) evolution, and replacement relations have to be used. The most popular is the Vogel-Fulcher-Tammann (VFT) dependence [48]: where T > T g and T 0 < T g is the extrapolated VFT singular temperature. The amplitude D = D T T 0 = const. D T denotes the fragility strength coefficient. Comparing Equations (11) and (12) one obtains for the VFT approximation E a (T) = RD/(1 − T 0 /T). Notwithstanding, there is an extensive evidence for the limited fundamental validity of the VFT relation in supercooled liquids ( [51,52] and refs therein). For liquid crystalline compounds the critical like dependence seems to offer the most reliable description [35,38]: Comparing derivatives of the SA Equation (11) and the 'critical' Equation (13) one can show that in the domain where Equation (13) can be applied the following relation should obey for experimental data [32]: where

H a (T) = H a (T)/R = d ln τ(T)/d(1/T) and the magnitude H a (T) is for the (apparent) activation enthalpy.
The linear regression analysis for the linear domain found (see Figure 4) yields optimal values of parameters for Equation (13). Figure 4 shows that in 5CB for the ultraviscous domain, T g < T < T g + 80 K the critical like description is associated with the power exponent ϕ ≈ 9. Notable that such behavior was predicted by Colby [49] within the so-called dynamic scaling model (DSM) and suggested as the universal pattern for supercooled glass-forming low-molecular-weight liquids and polymers. Unfortunately, this model-indication was finally not accepted and rejected check-in/out studies by other researchers [53]. However, the DSM-like scaling pattern seems to describe perfectly the primary relaxation time changes in rod-like liquid crystalline glass-formers on approaching the glass temperature [33,34].
It is well established experimentally and theoretically, within the mode-coupling theory (MCT), that the critical-like Equation (13) can also be applied in the high temperature 'dynamical domain. Such MCT type behavior is associated with T C >> T g , τ(T C ) ∼ 0.1 µs , and the power exponent ϕ = 1.4 ÷ 3. Its ability to portray experimental data is expected for T > T C + 10 K [33,34,[50][51][52]. Exactly such MCT-type dynamics occur in the isotropic liquid phase of 5CB, as shown in Figure 4. Notably, the presence of the I-N transition manifests as a distortion visible from T I−N − 10 K in the nematic phase to T I−N + 5 K in the isotropic liquid phase. Figure 5 presents a strong electric field (E) related changes of dielectric constant-the nonlinear dielectric effect (NDE)-in the isotropic liquid phase. There are qualitatively different patterns for dielectric constant ε(T) (the inset in Figure 1) and NDE = (ε(E) − ε)/E 2 ( Figure 5). This is associated with the fact that dielectric constant detects the averaged arrangement of permanent dipole moments. The definition of NDE directly suggests that in a 'critical' system, it registers the contribution solely from pretransitional/precritical fluctuations. Generally, this contribution is described by the relation [36,37,54]: where the order parameter related susceptibility χ T = χ 0 (T − T C ) −γ and for the mean square of the order parameter fluctuations: ∆M 2 V ∝ (T − T C ) 2β ; T C denotes here the critical temperature which for the case of 5CB can be associate with T * .  (1) and (16)).
In the case of the isotropic phase of the nematic liquid crystal, the rod-like symmetry causes the increase of the number of neighbors what leads to the crossover to the mean-field type behavior for which generally the dimensionality d = 4 is required. For the tricritical case it reduces to d = 3 [22]. Hence, in the isotropic liquid phase: ∆M 2 V = ∆ε 0 ∆ε f and χ T = χ 0 /(T − T * ), where ∆ε 0 is for the anisotropy of dielectric permittivity for the DC electric field and ∆ε f at the applied to measure radio-frequency. In the case of the low-frequency (LF) NDE [54,55]: For the LF NDE the time-scale associated with the weak measuring field is qualitatively lower than the lifetime of pretransitional fluctuation, and then during one NDE measurement cycle the avere of few/several prenematic fluctuations is registered.
Results confirming the portrayal by the above relation is shown in Figure 5. Notable is the fair agreement with the temperature behavior for the Cotton-Mouton effect, Kerr effect, or the scattered light (Equation (1)) [1][2][3][4]. However, NDE is the only method for which the measurement frequency can be changed, which also leads to the possibility of the precise determination of ∆T * for phase transitions not only to nematic but also to different smectic or cholesteric mesophases [55].

Discussion
Pentylcyanobiphenyl (5CB) is one of the most classical thermotropic, rod-like, liquid crystalline materials commonly considered as the model LC system and also widely used in applications [3][4][5][6]. It exhibits the simplest liquid crystalline mesophase-nematic, for which solely the spontaneous orientational ordering takes place. The discussion presented in the given report is based on dielectric studies, the research method particularly well focused on detecting unique features of LC materials [1][2][3][4][5][6]. Notably, that the analysis of experimental data explored the innovative distortions-sensitive and derivative-based approach what enabled to reduce or avoid factors, and consequently misleading conclusions [32]. New experimental studies in the isotropic liquid and in the nematic phase yielded clear evidence for exponents α ≈ α ≈ 1/2 (specific heat-related values) and β ≈ 1/4 (order parameter related exponent). Such values indicate the tricritical nature [22] for the weakly discontinuous I-N transition in 5CB. One should stress that the tricritical-like description for dielectric constant obeys in the whole tested nematic phase, between T I−N and T m , as well as in the isotropic liquid phase.
It is often claimed that changes of dielectric tensor in the nematic phase are primarily the consequence of nematic ordering of molecules, and thus pretransitional fluctuations in the isotropic phase are also the consequence of that [2][3][4][5][6]12]. This 'molecular approach' means that dielectric properties are primarily described by the molecular analysis based on the Kirkwood-Froelich concepts, eventually extending by the Maier-Saupe (MS) model. This lead to dimers, trimers, . . . and to more complex as 'cybotactic groups' appearing in the LC material as the substitute of pretransitional (prenematic fluctuations) [3,4,[56][57][58][59][60]. There are also claims that the observed weakly first-order nature of the N-I transition must be traced at the level of molecular interactions because just before the phase transition to the nematic phase the relaxation rate (linked to the relaxation time) strongly decreases from several 10 kHz down to a few kHz. In fact, such behavior is reported as the phenomenon extending in the temperature range ∆T ≈ 1 K [57][58][59]. The 'molecular approach is also linked to the maximum of ε(T) in the isotropic liquid phase observed at ca. T I−N + 5 K as the particularly important reference point in explaining phenomena observer in the isotropic liquid phase [58,59].
Commenting on the discussion 'molecular' vs. 'quasi-critical' approach for describing dielectric properties of LC materials, the author would like to indicate the following issues: (i) Molecular' strategies are not able to portray temperature changes of dielectric constant, whereas the 'quasi-critical' approach describes experimental results from the Solid-Nematic melting temperature up to even T I−N + 100 K in the isotropic liquid phase [25][26][27][28][29][30][31][32][33]. (ii) The longest collective relaxations time associated with prenematic fluctuations in the isotropic liquid is equal τ f luct. (T I−N ) ≈ 1 µs and for the single-molecule related relaxation time primary relaxation time τ(T I−N ) < 50 ns. This is associated with frequencies ('relaxations rates')~1 MHz and~20 MHz, respectively [38,61]. Values of relaxation rates associated with frequencies <10 kHz for which 'special molecular processes' in the immediate vicinity of the clearing temperature was suggested in refs. [57,58] can be expected for relaxation times occurying only in the nematic phase, for τ || [62]. (iii) Notable is the lack of the additional relaxation process emerging in the immediate vicinity of the clearing temperature in the isotropic phase, heuristically indicated in refs. [57][58][59]. This issue is discussed extensively in the Appendix A. (iv) In refs. [57,58] the 'stretching' of the I-N transition in the domain ∆T ≈ 1 was evidenced, basic on the scan of dielectric constant.
The vast majority of reports show a sharp, discontinuous, I-N transition and no 'temperature stretching' for the I-N transition is reported [1][2][3][4]19,28,31,33,[38][39][40][41]. In fact, such 'stretching' means that the I-N ceases to be the discontinuous phase transition, against the recalled broad experimental evidence. Notwithstanding, such 'stretched behavior' can appear under the (very) strong electric field [3,63]-much higher than one used for NDE tests. It can also appear due to admixtures/impurities as the consequence of the biphasic region [22,64,65] between the isotropic liquid and the nematic phase, in ref. [57] the latter was declared to be absent. However, a definitive statement regarding the biphasic region is possible only after supplementary visual observations [64,65].
One can add here that the 'quasi-(tri)critical' approach does need poorly-defined 'cybotactic' groups [1][2][3]. It can explain and describe in a consistent way, temperature evolutions of dielectric constant, nonlinear dielectric effect, Kerr Effect, Cotton-Mouton effect, (see above). Finally, the 'quasi-critical' scenario explains the isomorphism of the pretransitional behavior observed in studies as the function of temperature and pressure [66][67][68][69][70]. Notable, that evolutions of ε(T) in the isotropic liquid and ε || (T), ε ⊥ (T), ε mean (T) in the nematic are completely described, be 'quasi-critical' relations. Their changes on approaching the I-N transition, or rather T * and T * * singular temperatures, can be well explained by prenematic (isotropic liquid phase) or pre-isotropic (nematic phase) fluctuations which correlation length (and the volume occupied by them in the sample) and the lifetime increase according to relations [22,41]: where T X = T * is for the isotropic liquid related extrapolated temperature of the hypothetical continuous phase transition (T > T I−N = T * + ∆T * ) and T X = T * * denotes the nematic phase-related extrapolated temperature of the hypothetical continuous phase transition (T < T I−N = T * * − ∆T * * ) Finally, worth addressing is the still popular [1][2][3][4] Haller relation [14]. It assumes the clearing temperature T C = T I−N as the reference Equations (4) and (5), instead of the temperature of the (hypothetical) continuous phase transition T * * . This is the reason of the 'strange and scattered' values of the 'exponent' β , which has to be influenced by the value of the discontinuity ∆T * * = T * * − T I−N .
Regarding dynamics, it follows patterns developed for supercooled and glass-forming liquids. Notably, that this also covers non-supercooled domains within the nematic phase and isotropic liquid. Pretransitional fluctuations distort the glassy dynamics patterns only in the surrounding of T I−N . Supercooling of 5CB down to the glass transition is a challenging experimental task: usually, it crystallizes at T m ≈ 15 • C [3,4]. However, also above T m the dynamics exhibit the previtreous character associated with the glass transition hidden in the supercooled state at T g ≈ −68 • C. These facts are worth stressing because, in basic monographs on the Physics of Liquid Crystals [1][2][3][4][5], the crucial role of T g and T I−N as origins of 'far impact singularities' governing key static and dynamic dielectric properties is omitted.

Conflicts of Interest:
The authors declare no conflict of interest. Figure A1 shows the evolution of the imaginary part of dielectric permittivity (ε ( f )) in the isotropic liquid phase of 5CB. The spectrum consists of the primary relaxation loss curve and for the low-frequency domain-the ionic DC conductivity part. The applied log-log scale in the inset shows characteristic linear domains. For the primary relaxation curve, they are associated with the distribution of relaxation times for the high-frequency ( f > f peak ) and the low-frequency ( f < f peak ) parts of the primary loss curve. Such behavior was first indicated in ref. [71]:

Appendix A
Exponents/coefficients m and n describe the distribution of relaxation times. The reference is related to the single relaxation time Debye process, for which m = n = 1. The DC electric conductivity is determined as σ = ωε ( f ). The central part of the plot is given in the semi-log scale, to show the temperature dependence of the loss curve peak height. It is described by the relation resembling Equation (7), for the pretransitional anomaly of dielectric constant, namely: Notable that in Figure A2, only a single relaxation process is visible, also near the clearing temperature. An additional (absent) relaxation process should manifest as the deformation of the primary loss curve. The lack of such extra-process confirms the distortions sensitive and derivative-based plot shown in Figure A2. Namely, recalling Equations (A1) and (A2) one obtains: The scenario based on Equations (A4) and (A5) offer a reliable way of determining the values of the distribution coefficients m and n. It also provides a precise estimation of the peak frequency, coupled to the primary relaxation time: d log 10 ε f peak d log 10 f = 0 (A6) Figure A2. The derivative-based and distortions-sensitive analysis of the primary relaxation loss curve (Equations (A4)-(A6)) in the isotropic phase of 5CB, the immediate vicinity of the clearing temperature. The plot explores experimental data from Figure A1. Notable is the clear manifestation of the Jonsher's scaling (Equations (A1)-(A4)), although in a limited range of frequencies. There are no hallmarks of an additional relaxation process.
The above discussion is focused on the primary relaxation process. The insight into the spectrum of the real (ε ( f )) and imaginary (ε ( f )) parts of dielectric permittivity, for the broad frequency range extending from 1 Hz to 1 GHz, is shown in Figure A3. The focus of these studies is searching for relaxation phenomena in the range 10 kHz-1 KHz suggesting in the immediate vicinity of the clearing temperature in refs. [57][58][59] Figure A3. Broadband dielectric spectra for the real and imaginary parts of the dielectric permittivity spectrum in the isotropic liquid phase of 5CB. The tested range of frequencies extends from f = 1 Hz (log 10 f = 0) to f = 1 GHz (log 10 f = 9). Results are for the immediate vicinity of the I-N clearing temperature (blue) and remote from the I-N discontinuous phase transition (red). For ε ( f ) data horizontal parts are for the frequency-independent dielectric constant ε ( f ) = ε. To reveal subtle features of spectra results of the derivative analysis of experimental data are also presented (solid curves, the right scale). Links between results/data and the left & right scales are indicated by arrows. Figure A3 are for T I−N + 0.1 K (blue) and T I−N + 100 K (red). The picture contains experimental data (left scale) and their derivatives (right scale). The latter can reveal even subtle feature of detection processes. Regarding the real part of dielectric permittivity, the horizontal section is for dielectric constant ε ( f ) = ε = const where the frequency shift does not change the value of ε ( f ). The distortions-sensitive derivative analysis of experimental data d log 10 ε ( f )/d log 10 f in the lower part of the plot (the right scale, solid curve) shows that in the immediate vicinity of the clearing temperature, the static domain extends from ca. from f ≈ 10 MHz to f ≈ 100 Hz. Only for lower frequencies, the gradual impact of ionic impurities becomes more significant.

Results presented in
The results presented in the upper part of Figure A3 shows that for lower frequencies ε ( f ) follows only the DC-conductivity related pattern, i.e., σ = ωε ( f ) = 2π f ε ( f ) = const. Such behavior is additionally confirmed by the derivative-based and distortions-sensitive analysis d log 10 ε ( f )/d log 10 f (solid curves, the right scale). In the immediate vicinity of T I−N it extends from ca. f ≈ 10 MHz down to even f ≈ 1 Hz.
Concluding, results presented above show no hallmarks of an additional low-frequency relaxation process in the range between 10 KHz and 1 kHz, indicated in the immediate vicinity of the I-N clearing temperature in ref. [57].