Landau–de Gennes Model for the Isotropic Phase of Nematogens: The Experimental Evidence Challenge

TL;DR

This paper challenges the Landau–de Gennes model by analyzing experimental evidence of pretransitional effects in liquid crystals, revealing limitations of the model.

Contribution

The study provides a new model-based explanation for pretransitional effects in nematogens that the traditional Landau–de Gennes model cannot explain.

Findings

Pretransitional effects vary depending on molecular structural features like dipole moment loci and steric hindrance.

The interplay between observation and pretransition fluctuations time scales is essential for explaining experimental results.

A new coherent explanation is proposed for phenomena unaccounted for by the canonical Landau–de Gennes model.

Abstract

The Landau–de Gennes model is one of the most significant fundamental frameworks in The Physics of Liquid Crystals and Soft Matter Physics. It is validated by the universal parameterisation of the Cotton–Mouton effect, the Kerr effect, and light scattering in the isotropic phase of nematogens. However, as early as 1974, de Gennes identified the first two puzzling problems of this model. Over the following decades, this list has expanded. This report presents the first comprehensive analysis of these issues, with the explicit experimental reference. It focuses on the hardly coherently discussed pretransitional changes in the dielectric constant and the extension in a strong electric field, specifically the nonlinear dielectric effect (NDE). Notably, there are uniquely different pretransitional forms of pretransitional effects, depending on molecular structural features such as permanent dipole moment loci or a steric hindrance. It is tested for 5CB, 5*CB, and MBBA: nematogenic liquid crystalline materials that differ in the above features. The obtained specific pretransitional effects and the evidence for the essential importance of the interplay between observation and pretransition fluctuations time scales led to a new coherent, model-based explanation of all the discussed problems, which cannot be explained within the canonical Landau–de Gennes model.

Figures

Full text

1. Introduction

The isotropic phase of rod-like nematic liquid crystals (LCs) is a gateway to liquid crystal (LC) mesophases whose extraordinary properties inspired the research validating the Landau–de Gennes model, one of the crucial concepts in The Physics of Liquid Crystals [1,2,3,4,5,6,7]. Its significance includes the Soft Matter category [8,9,10,11,12], comprising liquid crystals, near-critical systems, polymers, colloids, micellar systems, supercooling glass formers, …. and model foods or biosystems [9,10,11,12,13,14,15]. It demonstrates the interdisciplinary nature of the field, which brings together physicists, chemists, and biologists [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15]. In recent years, soft matter models have also entered the fields of socio-economics [16] and topology-governed systems [17]. Despite qualitative differences, these systems share common scaling patterns. The isotropic phase of nematogens is also important in Critical Phenomena and Phase Transitions Physics due to a melting/freezing unusual discontinuous phase transition related to only one element of symmetry [18,19,20,21]. The isotropic phase of rod-like nematogenic LC has also proven to be significant for studying previtreous properties in glass-forming systems, i.e., for the Glass Transition Physics [22,23,24,25,26].

A unique significance of the isotropic phase of the nematogens also stems from the fact that it constitutes an experimental model system, facilitating the following cognitive feedback process: $[eqn]$ .

The Landau–de Gennes (LdG) model is a key reference in all of the aforementioned fields. Pretransitional effects studies in the isotropic liquid phase of rod-like nematic LC constitute the canonical reference for LdG model validation, most often recalled in monographs [2,3,4,5,6,7,8,9,10,12,19,20]. However, when compared with experimental results, several puzzling issues emerge. This report presents this peculiar experimental evidence and proposes coherent model explanations of emerging challenges.

Five decades ago, Pierre Gilles de Gennes published seminal papers discussing the anomalous increase in the Cotton–Mouton effect (CME) in the isotropic liquid phase of rod-like liquid crystalline (LC) compounds, as they approach the isotropic–nematic (I-N) transition [27,28]. To explain the surprisingly simple and common pattern of changes, de Gennes considered the expansion of the Landau model free energy ( $[eqn]$ ) using the uniaxial quadrupolar order parameter [27,28]:

$$[eqn]$$

Substituting the tensorial uniaxial order parameter $[eqn]$ , where $[eqn]$ refers to the component of the director vector with ‘nematic’ $[eqn]$ and $[eqn]$ and $[eqn]$ is the scalar metric of ordering, de Gennes obtained the expansion with the scalar metric [1,2,27,28]:

$$[eqn]$$

De Gennes introduced the extrapolated temperature of a hypothetical continuous phase transition in order to take into account the long-range pre-transitional effect associated with discontinuous phase transitions: $[eqn]$ , where $[eqn]$ is the I-N ‘clearing’ temperature and $[eqn]$ is the I-N discontinuity metric. Consequently, the amplitude of the second-order term in the above relation is given by $[eqn]$ . The last term in the above equations reflects the impact of the external field. Refs. [27,28] recalled Cotton–Mouton effect (CME: refractive index birefringence induced by a strong magnetic field).

In subsequent years, apart from the CME [27,28,29], such modelling was employed for the Optical Kerr effect (OKE) [30,31,32], Electro-Optic Kerr effect (EKE) [32,33,34], Rayleigh light scattering (LS, $[eqn]$ ) [35], turbidity ( $[eqn]$ ) [36], compressibility ( $[eqn]$ ) [37], and nonlinear dielectric effect (NDE, $[eqn]$ ) [38]. For each of these methods, the same temperature pattern for the strong and long-range pretransitional effect was reported [1,2,3,4,19,20,29,30,31,32,33,34,35,36,37,38]:

$$[eqn]$$

where $[eqn]$ is the amplitude related to the given method; parameters $[eqn]$ .

In the basic monographs for Liquid Crystals Physics, the linear changes in CME, KE, or LS reciprocals are presented as the crucial validation of the Landau–de Gennes (LdG) model, as well as a more general mean-field description validation in the isotropic liquid phase [1,2,3,4,5,6,7].

The grand success of the LdG model was a significant contribution to honouring Pierre G. de Gennes’ Nobel Prize in 1991 [39].

However, already in the first edition of the grand monograph The Physics of Liquid Crystals (1974, [1]), de Gennes pointed out some puzzling problems of the model:

• (Q1) ‘For the majority of second-order phase transitions the Landau model approximation does not occur (magnetic systems, superfluid Helium, …) showing for response functions the behaviour $[eqn]$ , with the exponent $[eqn]$ , but not $[eqn]$ ‘ [1].

It becomes even more puzzling when including later results for Electro-Optic Kerr effect (EKE) and NDE upon approaching the critical consolute point, with anomalies described by the exponents $[eqn]$ and $[eqn]$ , respectively [40,41].

• (Q2) ‘Why is the temperature singular $[eqn]$ is so close to $[eqn]$ , i.e.,  $[eqn]$ ?’ [1].

It means that the experimentally detected discontinuity via properties encountered in Equation (3) is equal to $[eqn]$ [1,4,19,21,29,30,31,32,33,34,35,36,37,38]. Notably, the mean-field analysis via basic Landau–de Gennes or Maier–Saupe models gives $[eqn]$ [1,4,42,43,44].

The results obtained in subsequent decades showed that this list can be significantly longer, namely, the following:

• (Q3) For $[eqn]$ , $[eqn]$ , and $[eqn]$ , pretransitional changes and characteristic deviations from the LdG Equation (3) appear in the immediate vicinity $[eqn]$ . They are commented only heuristically via the impact of not-defined ‘cybotactic groups’ [3,4,31,32,45,46,47,48,49,50,51,52].
• (Q4) Heat capacity/specific heat shows the pretransitional anomaly $[eqn]$ , with the exponent α = 1/2, in the isotropic liquid and in the nematic phases [4,19,47,48,49,50,51]. The basic LdG model yields in the isotropic liquid phase α = 0 [4,5,6,53,54,55,56]
• (Q5) The puzzle regarding the explanation of the pretransitional anomaly of the dielectric constant $[eqn]$ existed [57,58,59,60,61,62,63,64]. It includes the explicit dependence on the molecular structure [63,64].
• (Q6) Nonlinear dielectric effect (NDE) is the strong electric field-related extension of dielectric constant [65,66,67]:

$$[eqn]$$

Generally, the infinite rise in NDE upon approaching a critical point is expected: $[eqn]$ . However, for some exceptional LC materials, the ‘negative’ pretransitional anomaly was found: $[eqn]$ [68].

The unusual, universalistic pretransition changes described by Equation (3) served as the inspiration and origin for the Landau–de Gennes model. Hence, significant experimental problems may raise both concerns and a call for clarification.

Notably, the Landau–de Gennes model has remained the subject of intensive theoretical research for decades, aimed at testing its peculiarities and the limits of its applicability [1,2,3,4,5,6,7,8,9,15,69,70,71,72,73]. However, this report refers to the puzzling experimental evidence (Q1–Q6) that has not been explained so far, and is essential for the ultimate validation of the model.

To address this challenge, focused studies of the dielectric constant and NDE are essential. These results are presented below and are followed by model considerations in the Section 3. They comment on and explain the puzzling experimental features recalled in questions Q1–Q6.

2. Results

The results presented below are for three nematogenic rod-like liquid crystalline compounds, with the following mesomorphism:

• (i)4-cyano-4-n-pentylbiphenyl (pentylcyanobiphenyl, 5CB), with the following mesomorphism [4,52]:

$$[eqn]$$

• (ii)4-(2-methylbutyl)-4′-cyanobiphenyl (isopentylcyanobiphenyl, 5*CB [69,74,75]:

$$[eqn]$$

• (iii)p-methoxy-benzylidene p-n-butyl-anilline, MBBA [4,21,52]:

$$[eqn]$$

5CB and 5CB share the same chemical formula C_18_H_19_N, but they differ slightly in molecular structure, as depicted in Figure 1 and Figure 2 below. 5CB essentially crystallizes on cooling. Supercooling is only possible when fast cooling is matched with a specially prepared sample [74]. 5CB can be supercooled at an arbitrary cooling rate down to the glass temperature that occurs in the chiral nematic phase at $[eqn]$ [68,76,77]. Both 5CB and 5*CB are associated with the permanent dipole moment $[eqn]$ [4,15], approximately parallel to the long molecular (cyanobiphenyl) part of molecules. It leads to the positive anisotropy of dielectric permittivity: $[eqn]$ for 5CB [4].

For MBBA, the permanent dipole moment is almost perpendicular to the long molecular axis, $[eqn]$ , and it is located at the angle $[eqn]$ with respect to the long molecular axis. It leads to the negative anisotropy of dielectric permittivity, $[eqn]$ [75]. The tested compounds were purchased from Fluka and carefully degassed before the experiments.

In conclusion, the tested compounds are characterized by qualitatively different positions of the permanent dipole moment (5CB and MBBA) and the presence of structural steric hindrance (5CB). These factors have a significant impact on the form of pretransitional effects in the isotropic liquid phase, as noted in the literature and as presented in the results below. This is a unique situation, as the standard properties studied in the isotropic phase, such as compressibility (usually tested via KE or IR measurements) and heat capacity (specific heat), consistently exhibit the same pretransitional patterns. This includes 5CB, 5CB, and MBBA, tested in the given report.

Figure 1 and Figure 2 present examples of registered spectra for the real part of dielectric permittivity in 5CB and 5*CB, focusing on features relevant to the physical properties discussed in the given report. High-resolution spectra for MBBA have only recently been reported in ref. [64] Both figures show the same pattern, with a long horizontal domain where the frequency shift does not change the value of $[eqn]$ . This is the static domain, defining the dielectric constant, i.e., $[eqn]$ [65].

When increasing the frequency above the static domain $[eqn]$ drops, which is the hallmark that the long molecular axis gradually ceases to follow changes in the electric field. This is the onset of the dynamic domain used for testing relaxation processes, which are not the focus of the given report. On lowering the frequency below the static domain, values of $[eqn]$ ‘dramatically’ rises. It is the low-frequency region, and the rise is generally linked to ‘ionic impurities’, whose impact can be detected for low enough excitation frequencies [60,61,65]. For the authors, an alternative explanation via translational shifts in the basic molecules can yield a similar registered output [64].

Generally, experimental studies of dielectric constant temperature evolution are related to the real part of dielectric permittivity scans at a selected, constant frequency, most often in the range from 10 kHz to 100 kHz, and are carried out [4,21,56,57,58,59,60,61,62,63,64]. Such a scenario can be successful for 5CB and also for MBBA, with a similar $[eqn]$ spectrum. However, in the isotropic phase of 5CB, the static domain strongly shifts on cooling (see Figure 2), which is related to the enormous temperature range possible for testing and approaching the glass temperature. It causes the above constant frequency scenario to fail. For the 5CB dielectric constant studies, the focus should be on the middle of the static domain, shifted towards lower frequencies upon cooling. As visible in Figure 2, for 5*CB, the static domain extends from ca ~3 kHz to well below 1 Hz just above the clearing temperature. For 5CB and also for MBBA, it is almost always related to the frequency range $[eqn]$ .

2.1. Pretransitional Behaviour of Dielectric Constant

Figure 3 shows dielectric constant changes in the isotropic liquid phase of MBBA, up to $[eqn]$ . There is an explicit linear behaviour, with no hallmarks of a pretransitional effect for $[eqn]$ , i.e., upon approaching the I-N transition. Nevertheless, the derivative analysis of high-resolution data in ref. [64] revealed weak hallmarks of the pretransitional behaviour for $[eqn]$ tests.

Figure 4 presents changes in the dielectric constant in the isotropic liquid phase of 5CB and 5*CB, using the normalized scales. The plot shows the explicit manifestations of pretransitional effects, namely

$$[eqn]$$

where the clearing temperature $[eqn]$ is related to I-N transition for 5CB and I-N* for 5*CB.

Such pretransitional behaviour is particularly evident in the inset in Figure 5, for the focused behaviour in the close vicinity of $[eqn]$ .

The curves portraying pretransitional changes are related to the following relation [56,57,58,59,60]:

$$[eqn]$$

where $[eqn]$ is for parameters describing the hypothetical continuous phase transition, at the tested temperatures $[eqn]$ ; the exponent is related to the one describing the mean-field related heat capacity anomaly, namely φ = 1 − α = 1/2 [57,58].

For 5CB and 5CB, dielectric constant tests enable determining the values of the discontinuity metric: for both compounds within the limit of the experimental error, $[eqn]$ . Worth stressing is the range of the extreme range of pretransitional effects and significantly lesser ‘strength’ of the pretransitional effect for 5CB in comparison with 5CB. It is reflected by the relation of amplitudes $[eqn]$ , as results from the data given in Table 1 and visible in Figure 4.

There are essential problems with the direct derivation of the relation for $[eqn]$ directly from the Landau–de Gennes model, and includes the molecular structure impact.

To portray $[eqn]$ anomaly in the isotropic phase, Drozd-Rzoska [58,63] recalled Mistura’s [78] model reasoning for precritical changes in the dielectric constant in the homogeneous phase above the critical consolute point (CCP) and gas–liquid–critical point (GLCP). Mistura indicated the link to the critical behaviour of the specific heat/heat capacity $[eqn]$ [78]:

$$[eqn]$$

where $[eqn]$ denotes the critical temperature.

For systems with CCP and GLCP, the critical exponent α ≈ 0.115. It is the value for the $[eqn]$ (space dimensionality) and $[eqn]$ (order parameter dimensionality) universality class [79]. The Mistura model concept [78] served as the essential reference for Sengers et al.’s [80] model analysis, leading to an equation that portrays $[eqn]$ changes upon approaching CCP and GLCP, solving theoretically the challenge that has existed since 1934 [81]. The functional form of the model relation derived by Sengers et al. [80] agrees with Equation (6), although it additionally includes critical corrections-to-scaling terms, important as we move away from the CCP, starting only from $[eqn]$ above $[eqn]$ [82,83].

Drozd-Rzoska [58,63] noted that Mistura’s model [78] analysis does not explicitly refer to systems associated with the critical universality class $[eqn]$ , but addresses generally a supercritical liquid (above the critical point) with the dominant impact of pretransitional, critical fluctuations. For the isotropic phase of nematogens, prenematic fluctuations are observed, and the critical temperature can be linked to the extrapolated singular temperature $[eqn]$ , for instance, of a pseudospinodal type [84]. Drozd-Rzoska further recalled the mean-field character of the isotropic phase, leading to the critical exponent α = 1/2 for the specific heat anomaly. Subsequently, the integration of Equation (7) can directly yield Equation (6) when additionally taking into account the ‘background’ term reflecting linear temperature changes in the nondipolar polarization component [58,63]. Finally, for the mean-field description, correction-to-scaling terms, appearing in Sengers et al.’s [80] model, are inherently absent.

In 1983, Thoen and Menu [57] carried out high-resolution dielectric constant measurements in the isotropic liquid phase on n-octyl-cyanobiphenyl, a rod-like LC compound with $[eqn]$ mesomorphism. They demonstrated a fair portrayal using Equation (6), with the heuristic explanation indicating some informal similarity to the behaviour in the homogeneous phase as it approaches the critical consolute point. The Drozd-Rzoska model [58,63] provided a formal justification for this preliminary result.

2.2. Pretransitional Behaviour of Nonlinear Dielectric Effect

This section presents the essential evidence for the dielectric constant counterpart of a strong electric field—the nonlinear dielectric effect (NDE) [65,66,67], defined by Equation (4). The pattern of NDE temperature changes in the isotropic phase of nematogenic LC compounds can be directly derived from the Landau–de Gennes model, similarly to other properties encountered in Equation (3). Recalling the analysis from refs. [85,86], one obtains the following:

$$[eqn]$$

Consequently, the NDE reciprocal should follow a linear pattern, as indicated in Equation (3). It essentially simplifies the analysis of experimental results.

In the above relation, $[eqn]$ . The coefficient ‘a’ is related to the second-order amplitude in Equation (2) and defines the amplitude of compressibility, or, alternatively, orders parameter-related susceptibility $[eqn]$ with $[eqn]$ ; $[eqn]$ is the local field-related factor reflecting the fact that the electric field acting on a molecule within a dielectric material differs from the externally applied one; $[eqn]$ and $[eqn]$ are metrics of molecular anisotropies of ‘dielectric constant’ for the hypothetical zero-frequency limit and the applied measurement frequency. This is related to the action of the DC strong electric field and the weak scanning field, respectively.

Recalling experimental conditions, values of molecular anisotropies in Equation (8) should be associated with the real part of dielectric anisotropy determined for frequencies in the static domain for LC samples in the nematic phase perfectly oriented by a strong magnetic field in directions perpendicular and parallel to the scanning electric field. Hence, $[eqn]$ in the static domain and $[eqn]$ , i.e., it can also cover frequencies beyond the static domain.

Figure 5 shows the NDE pretransitional rise in the isotropic phase of 5CB and MBBA for two scanning frequencies of the weak electric field. For MBBA, the pretransitional effect is ca. 50× weaker than for 5CB: it means it is of the order $[eqn]$ for MBBA and $[eqn]$ for 5CB.

This difference does not influence the resolution of registered NDE values. For the given apparatus concept, one can detect even the lowest possible contribution to NDE from statistical polarization fluctuations in the nondipolar liquid, of the order $[eqn]$ [65,66]. The inset in Figure 5 is for the reciprocal of NDE changes from the main part of the plot, focused on the visual validation of the LdG model-related temperature pattern (Equations (3) and (8)). Notable is the range of such descriptions extending even above $[eqn]$ .

For MBBA, the shift in the scanning frequency does not influence the phenomenon. However, such an impact significantly manifests in 5CB for the scanning frequency $[eqn]$ . It should be stressed that both frequencies are related to the static domain of the dielectric constant, as shown in Figure 1 for 5CB and in ref. [64] for MBBA.

Figure 6 presents the pretransitional behaviour in the isotropic liquid phase of glass-forming 5CB. The striking feature is the anomalous pretransitional effect, which tends toward negative values. The impact of the scanning frequency is also visible. The lower of the tested frequencies corresponds to the static domain of the dielectric constant, as shown in Figure 2. Generally, for obtaining the ‘negative’ pretransitional anomaly for NDE within the Landau–de Gennes model, opposite signs of anisotropies $[eqn]$ and $[eqn]$ are required in the context of Equation (8). However, it is not possible to obtain these single-molecule-related properties for 5CB. Nevertheless, the question regarding the parameterization arises.

The standard way of testing NDE, KE, $[eqn]$ ….pretransitional patterns is the ‘reciprocal plot’, as shown in the inset in Figure 5. However, for 5*CB, the crossover $[eqn]$ excludes such analysis.

To overcome this problem, the following transformation of experimental data has been carried out:

$$[eqn]$$

where constant parameters are related to the following linear behaviour: $[eqn]$ and $[eqn]$ .

The inset in Figure 6 shows the results of such data transformation, based on Equation (9). The fair linear behaviour validates the mean-field type evolution $[eqn]$ , with the negative amplitude $[eqn]$ .

As for the behaviour well remote from the clearing temperature, one may expect that the influence of pretransitional fluctuation diminishes. The form of the temperature evolution in this domain, shown by the brown line, allows us to recall the classic Piekara model for NDE changes in molecular liquids with the tendency to pair dipole–dipole coupling [65,87]:

$$[eqn]$$

where $[eqn]$ is for the number of permanent dipole moments in a unit volume, $[eqn]$ denotes the permanent dipole moment metric, $[eqn]$ is the Boltzmann factor, and $[eqn]$ denotes Piekara’s factor describing the possible dipole–dipole coupling and for the antiparallel arrangement $[eqn]$ ; $[eqn]$ is the local electric field factor.

Originally, Equation (8) was introduced to explain the positive NDE value in nitrobenzene, with the explicit antiparallel dipole–dipole local arrangements [65,87].

3. Discussion

This section addresses the experimental challenges presented in questions Q1–Q6, with emphasis on the experimental results discussed above. Section 3.1 provides some cognitive background comments for the Physics of Critical Phenomena, which are significant for the subsequent model discussion.

3.1. Critical Fluctuations and Mean-Field Behaviour

The isotropic–nematic transition in rod-like molecular systems exhibits unique features of an experimental model system, i.e., a system that can be relatively simply mirrored by a conceptual/theoretical model, offering exceptional $[eqn]$ cognitive feedback analysis [4,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,42,43,44,45,46,47,48,49,50,51].

For this path, of particular importance is the molecular field interaction model for uniaxial rods considered by Onsager in 1946 [88], later developed by Flory [89] and followers [90,91,92,93,94]. These theoretical studies, later supplemented by simulations, showed that a specific I-N solidification process is associated only with uniaxial orientation, i.e., a single element of symmetry, freezing. The discontinuous character, with step changes in relevant physical properties, can be considered a generic feature of such systems, as indicated already by Onsager [88].

It is essentially a comment on the discontinuous character of the I-N transition in rod-like LC systems. However, a question arises as to why this discontinuity is so small, as pointed out by de Gennes in his Q2 question.

It seems that a simple heuristic answer might be the ‘freezing’ limited to only one element of symmetry. Two further essential experimental facts should also be considered: (i) the extremely long-range pretransitional effects, which are generally absent for standard discontinuous melting/freezing transitions, and (ii) the mean-field nature of these pretransitional effects.

Regarding the first issue, for weakly discontinuous phase transitions, the nature of the phenomenon often causes it to be associated with a critical-like anomaly ‘hidden’ just below the transition. An example is pseudospinodal behaviour in the close surroundings of a critical point [84]. Such a ‘pseudospinodal’ picture can resemble the behaviour observed for the I-N phase transition [84,95,96].

Generally, pretransitional effects are caused by the appearance of multimolecular pretransitional fluctuations, which exhibit locally features of the next, approaching phase. For the given case, they are prenematic fluctuations. The size (correlation length, $[eqn]$ ) and lifetime ( $[eqn]$ ) of pretransitional fluctuations critically increase upon approaching the singular temperature, critical $[eqn]$ , or the pseudospinodal $[eqn]$ temperatures [79,84]:

$$[eqn]$$

where for the cases of GLC and CCP, the exponent $[eqn]$ , and the dynamic exponent $[eqn]$ for the conserved and $[eqn]$ for the non-conserved order parameter.

For the I-N transition, one can assume the following: $[eqn]$ : $[eqn]$ and $[eqn]$ [81].

Notably, values of critical exponents are described by space $[eqn]$ and order parameter $[eqn]$ dimensionalities. Hence, all near-critical systems can be arranged in ( $[eqn]$ universality classes [79]. Remote from the critical/singular temperature, the correlation length of critical fluctuations reduces to as few as two, or even fewer, molecules. It is matched with their extremely short lifetime. Consequently, the impact of standard intermolecular interactions may become dominant again. This crossover is described by the Ginzburg criterion [97] in the cases of GLC and CCP. It is also the boundary between the non-classical description and the ‘classic’ mean-field type behaviour. However, pretransitional fluctuations remain relevant because the system’s properties are still described in relation to the critical/singular point. It is often emphasized that in the classic domain, the long-range intermolecular interactions, in comparison to critical fluctuations, cause the influence of the latter to be averaged out [85,86]. However, there are systems where a mean-field description occurs in the presence of large multi-element pretransitional fluctuations. Ferroelectrics in the paraelectric phase can be a canonical example. They are associated with long-range Coulomb interactions, as well as a dominance of uniaxiality, which appears to be essentially important for the appearance of the mean-field behaviour [98].

In practice, the ‘infinite’ range of interactions means the possibility of interacting with a larger number of molecules (or species) than is possible via a direct, ‘geometric’ environment of the same type of molecules. This can be achieved by a sufficiently high effective dimension of space. For a single, isolated critical point, this is related to $[eqn]$ , with the inherent mean-field description. The related critical type is universal for any system, regardless of its microscopic features. It is expressed by critical exponents: for compressibility γ = 1, correlation length $[eqn]$ , and order parameter $[eqn]$ . For the specific heat (heat capacity), α = 1/2 in the low-temperature phase (below $[eqn]$ ) and α = 0 (no anomaly) in the high-temperature phase [19,98].

Another type of mean-field behaviour and pretransitional effects appears in the vicinity of the tricritical point (TCP) where three lines of critical points intersect. In some systems, a simplified symmetric TCP phenomenon can appear. It manifests via a smooth crossover from the line of discontinuous to the line of continuous phase transitions [99,100]. It was noted for some magnetic systems or phase transitions associated with the superfluid behaviour [99,100]. For TCP, the border dimensionality reduces to $[eqn]$ , and the pretransitional effects are coupled with the following values of exponents: γ = 1, $[eqn]$ , $[eqn]$ , and α = 1/2, both above and below the TCP singular temperature [100,101,102,103,104,105].

3.2. I-N Transition Mean-Field Nature

The evidence for universalistic long-range anomalous changes in the Cotton–Mouton effect, Kerr effect, and Rayleigh light scattering, obtained at the turn of the 1960s and 1970s (see Equation (3)) [27,28,29,30,31,32,33,34,35,36,37], has become the crucial evidence for validating the mean-field nature of the I-N transition and the isotropic phase. The hallmark feature is the exponent γ = 1 for describing critical changes in compressibility, linked to properties recalled in Equation (3). It should be noted that at that time, only Rayleigh light scattering could be explicitly linked theoretically to compressibility. The model relationship between the Kerr effect, NDE, and a simple extension to CMEs was demonstrated only in 1993 [40].

The behaviour of these properties can be phenomenologically described by the basic LdG model formulation, as defined in Equations (1) and (2). However, it also yields the exponent α = 0 in the high-temperature phase, which means there is no anomaly for the specific heat and related properties in the isotropic phase, and is in explicit disagreement with the experimental evidence for the specific heat [4,19,53,54,55,56] or dielectric constant (see the Section 2).

The problem was solved by implementing the more general Landau & Ginzburg model expansion, which takes into account the impact of supplementary factors, particularly reflecting weak local density fluctuations [4,19,42,43,44]. Nevertheless, it introduces additional parameters to relations describing pretransitional changes, making the fitting of experimental data tricky. To reflect the collective behaviour expressed by Equation (11), the gradient term for local changes has been included. This is the extended Landau–Ginzburg–de Gennes (L-G-dG) model [4,19].

The agreement with the specific heat can also be reached by extending the Landau–de Gennes expansion up to the term proportional to $[eqn]$ in Equation (2), which implies the near-tricritical nature of the I-N transition [4,19,101,102,103,104,105] rather than TCP pseudospinodal-type behaviour. It naturally yields the required form of the pretransitional anomaly for the specific heat/heat capacity.

The discussion regarding the basic mean-field vs. TCP-type origin of the I-N transition has a long history [4,19,50,101,102,103,104,105], and only unambiguous experimental results can be decisive. Ultimate estimations of the order parameter exponent value in the nematic phase can be essential since the value $[eqn]$ for TCP and $[eqn]$ for the basic mean field [19,79]. However, matching experimental data, most often based on refractive index or dielectric constant measurements, appears to be a challenging task due to the discontinuity of the I-N phase transition, i.e., experimental data are available only well beyond the singular temperature. These biassing factors can be minimized by the linearized, distortion-sensitive analysis, which, for nematogenic 8OCB and 5CB, yielded α ≈ 0.25, in agreement with the TCP hypothesis [56,63].

Another ultimate test of the ‘TCP hypothesis’ can be the experimental demonstration of the above-mentioned crossover from a discontinuous phase transition line to a continuous phase transition line. Several experimental studies have suggested that this behaviour can be achieved for the I-N transition under a sufficiently strong external electric field [101,102,103,104,105]. The authors of this report have extensive experience in NDE and Electro-Optic Kerr effect studies, which are inherently associated with strong electric fields. In our opinion, the results reported in refs. [101,102,103,104,105] and related cannot be considered as a definitive validation. First, they involve electric fields up to $[eqn]$ [101,102,103,104,105]. Such extreme electric fields can introduce numerous parasitic impacts, for instance, related to electrical conductivity and local heating. Second, doubt is related to the visible ‘stretching’ of the I-N transition in the range of even 1–2 K, rather than the strong reduction in the discontinuity metric.

For the authors, another factor can be decisive for the smooth crossover from the line of discontinuous to continuous phase transitions: the length of the rod-like molecule, i.e., the primary structural factor. This may be indicated by the variability of the discontinuity value $[eqn]$ in the homologous series of n-cyanobiphenyls to which 5CB belongs. They range from ~10 K for 14CB to ~0.7 K for 4CB, as shown in the Appendix A. See also refs. [106,107] for similar evidence in other LC homologous series. The rising length of rod-like molecules, matched with the rising ability to undergo a uniaxial ordering, can be a decisive parameter defining the reduction in discontinuity, and a continuous transition can appear below some minimal molecular length. This interpretation supports the recently obtained continuous phase transition, within the experimental error limit, for approaching I-N transition in 5CB on compressing, with a frustrated impact of a small addition of paraelectric BaTiO_3_ nanoparticles [108].

In standard temperature tests for a given LC compound, only one point from such near-TCP line is empirically available. For 5CB, it is associated with a $[eqn]$ discontinuity.

3.3. Discontinuity of the I-N Transition

De Gennes’ question (Q2) indicates the unique ‘weakness’ of the I-N transition, which reflects small values of $[eqn]$ . It can be precisely determined using properties encountered in Equation (3). This issue was partially discussed at the end of the last subsection. A notable commentary on de Gennes’s concern may be the qualitatively large model estimates of the discontinuity metric. For the basic LdG model, it is $[eqn]$ , and for the alternative Maier–Saupe model, it is even $[eqn]$ [1,4,42,43,44]. Approximate agreement with experiment was obtained by Mukherjee in 1998 [44], assuming the TCP character of the I-N transition.

The discontinuity metric $[eqn]$ value is considered a particularly significant experimental checkpoint for theoretical models [1,2,3,4,5,6,19,42,43,44]. It is also a recognized ‘material’ characteristic of LC materials [4]. However, it is rarely noted that this meaning of $[eqn]$ has some limitations. In nematogenic LC materials, the discontinuity value strongly increases with compression; for example, in 5CB, it can reach ~7 K at a relatively moderate pressure of $[eqn]$ [63]. As indicated in refs. [63,106], at least in the homologous series nCB, pressure changes lead to the appearance of the isotropic–nematic–smectic A triple point associated with the only $[eqn]$ value, which is specific, i.e., only for a given LC material.

3.4. Dielectric Constant: Model Explanation

In the Section 2, Drozd-Rzoska’s model analysis [58,63] was presented to describe the pretransitional behaviour of the dielectric constant in the isotropic phase of nematogens. It employed Mistura’s (1974) [78] and Sengers et al.’s (1980) [80] concepts, developed to explain the critical anomaly of the dielectric constant as one approaches the critical consolute point (CCP) in mixtures of limited miscibility.

However, one can also consider the application of the droplet model introduced by Oxtoby for the approximate description of critical effects in liquids (1977) [109], and next implemented to describe anomalies in dielectric properties, namely the dielectric constant, Kerr effect, and NDE in critical binary mixtures (1979) [110]. It led to the following relation of dielectric constant excess change on approaching the critical consolute point in binary mixtures [80,83]:

$$[eqn]$$

The Oxtoby model [109,110] considers exclusively the dielectric constant critical anomaly in a system, which arises due to the dielectric constant’s local excess associated with precritical fluctuations, approximated as droplets with a radius described as the correlation length (Equation (11)). Such a picture can also be directly implemented in the isotropic–nematic phase, if system-specific features are taken into account. It includes elements of the inherent uniaxial anisotropy specific to the next nematic phase, as well as a semi-continuous and mean-field type nature of the transition. Following this reasoning, Equation (12) directly transforms into the isotropic phase-related Equation (6), with $[eqn]$ .

The characteristic feature of continuous/critical phase transition is the appearance of pretransitional fluctuations with local characteristics recalling the next approaching phase. For the isotropic phase of nematogens, it is related to prenematic fluctuations. They have to show the basic feature of the nematic phase, i.e., the uniaxial ordering, whose dominant direction is indicated by the director. It is significant that directors $[eqn]$ and $[eqn]$ , indicating the dominant uniaxial ordering, are equivalent [4]. The latter feature directly leads to the approximate ‘cancellation of contribution from permanent dipole moments, if they are parallel to the long molecular axis. Notably, the cancellation does not require intermolecular dipole–dipole coupling and can be purely ‘statistical’. Consequently, the dielectric constant within fluctuations is essentially less than that of the ‘chaotic’ surrounding. A ‘contrast factor’ between fluctuations and their liquid-like surrounding appears. It yields explicit conditions for the above-considered modelling, which recalls the droplet model background.

The same framework excludes the appearance of the dielectric constant pretransitional effect in the isotropic phase of MBBA, specifically with respect to the transverse permanent dipole moment. In the given case, the ‘cancellation’ of the permanent dipole moment has to be absent. There is no ‘contrast factor’, and therefore no pretransitional anomaly.

Regarding the significant difference in the manifestation of pretransitional effects between the isotropic phases of 5CB and 5*CB, it is worth noting that they are isomeric compounds, sharing the same permanent dipole moment with respect to the long molecular axis, but showing a slight structural difference, as illustrated in Figure 1 and Figure 2.

For 5CB, a steric hindrance associated with the molecular structure distorts the prenematic arrangement within fluctuations, as illustrated in Figure 7. For 5CB, the contribution of permanent dipole moments can be almost totally cancelled. For 5CB, one can expect a decoupling: $[eqn]$ . Only the parallel component can be subject to prenematic cancellation, and then a non-vanishing transverse component of the permanent dipole moment ( $[eqn]$ ) for the prenematic fluctuations remains. The limited cancellation of the permanent dipole moment leads to a lesser contrast factor in 5CB than in 5CB. It further leads to a lesser manifestation of dielectric constant pretransitional anomaly in 5CB than in 5C, as visible in Figure 4 and Table 1.

Dielectric constant scans register the effective changes in this quantity in the total volume of the sample, which in a given case resembles a specific ‘critical colloid’ consisting of the isotropic volume part with chaotically arranged molecules and a part occupied by prenematic (locally ordered) fluctuations: $[eqn]$ .

The volume occupied by the fluctuations increases on cooling $[eqn]$ , which is related to the simultaneous decrease in the $[eqn]$ contribution. Both contributions to the registered dielectric constant equilibrate at some temperature for which $[eqn]$ . Below this crossover temperature, in the immediate vicinity of the I-N clearing temperature, the impact of fluctuations dominates.

This is a case of the indirect detection of the impact of critical fluctuations.

3.5. Nonlinear Dielectric Effect: Timescale Meaning and the Model Explanation

The Kerr effect (KE), the Cotton–Mouton effect (CME), light scattering (LS), and the nonlinear dielectric effect (NDE) are methods that directly register the response from multimolecular fluctuations. NDE introduces the issue of the interplay between the timescale associated with the scanning electric field employed by the method and the relevant system timescale. The latter is associated with the prenematic fluctuations lifetime, which critically rises on cooling, as shown in Equation (11). For 5CB, it terminates at $[eqn]$ , which is ca. 3 decades less than the molecular primary relaxation time $[eqn]$ [63,74].

KE, CME, and LS are methods employing light for scanning properties. For the ‘red’ light limit, it is related to frequency $[eqn]$ , and then, the scanning timescale $[eqn]$ . Hence, for the mentioned light-related methods, $[eqn]$ at any temperature in the isotropic liquid phase. NDE detects the system’s properties using a radio-frequency weak electric field, in the $[eqn]$ range [66]. It yields the scanning timescale $[eqn]$ from ca. $[eqn]$ to 0.1 μs. The NDE scanning timescale can coincide with the system timescale ( $[eqn]$ ) in the isotropic phase.

Experimental evidence for OKE, EKE, CME, and LS shows the (apparent) agreement with the LdG model (Equation (3)), regarding both the temperature behaviour pattern and the method-related amplitude $[eqn]$ . The latter includes significant molecular properties, particularly anisotropies of the refractive index $[eqn]$ or dielectric constant $[eqn]$ . Following the discussion in the given report, the detection was carried out by OKE, EKE, CME, and LS methods, which were used to ‘observe’ fluctuations as an average from individual molecules constituting the average fluctuation.

The scanning timescale of NDE causes the entire pretransitional (prenematic) fluctuation to be ‘observed’/scanned. In the particular case of low-frequency (LF) NDE, it is possible that $[eqn]$ , and the effective average response from an assembly of fluctuations was tested.

In 1993, Rzoska [40] proposed a model for explaining the mystery of NDE and EKE pretransitional effects in the homogeneous phase of critical binary mixtures. It led to the output equation [40]:

$$[eqn]$$

where $[eqn]$ is the mean of the local ordered parameter fluctuations square and $[eqn]$ is order parameter related susceptibility, i.e., the compressibility; $[eqn]$ is the physical magnitude related model constant.

In ref. [41], it was extended for the supercritical domain above the gas–liquid critical point.

Significant are the unique conditions for a ‘standard’ supercritical liquid under a strong electric field [40]:

• Initially, spherical/isotropic critical fluctuations become elongated and semi-uniaxial under the strong electric field, so the correlation length becomes uniaxial $[eqn]$ , and additionally shows a ‘mixed criticality: $[eqn]$ . That is, it follows the standard for critical mixtures non-classical pattern, and $[eqn]$ , i.e., it follows the mean-field pattern.
• The different experimental definitions of NDE are significant: $[eqn]$ , and EKE, $[eqn]$ .

The above model can be extended for the isotropic phase above the I-N transition [40,41]. First, one should consider the substitution $[eqn]$ in Equation (13). Next, one should take into account that for deriving Equation (3), only the term $[eqn]$ is taken into account. It means that in the LdG series, given by Equations (1) and (2), other terms can be negligible for the methods encountered in Equation (3). It can be explained by the phase transition discontinuity. The authors of this report indicate that a strong external electric (or magnetic) field can change the tested system, creating a ‘uniaxial’ semi-critical colloid. Consequently, under such conditions, the system may return to its basic mean-field state.

Following Equation (13), it leads to the following dependence in the isotropic liquid phase of nematic LC:

$$[eqn]$$

The mean-field character causes the decoupling of the order parameter fluctuations and the temperature behaviour governed by the susceptibility: $[eqn]$ .

For NDE, local order parameter fluctuations are associated with local dielectric constant fluctuations in the static domain in a hypothetical zero-frequency limit. It is linked to the DC strong electric field, inducing the anisotropy ( $[eqn]$ ), i.e., zero-frequency extrapolation of the static domain. It reduces the detection of a fluctuation $[eqn]$ to a single-molecule-related anisotropy of the dielectric constant. For $[eqn]$ , it is not possible, and the average single fluctuation features are detected.

Equation (14) appears to resemble the LdG Equation (8). However, the parameters in Equation (14) are related to the ‘whole’ fluctuations, detected via the relevant timescale to the scanning frequency. For LdG Equation (8), there are explicit molecular properties.

Following the information above, one can consider the nature of pretransitional NDE anomalies presented in Figure 5 and Figure 6:

For MBBA, $[eqn]$ , and one can consider the uniaxial ordering of rod-like molecules with the transverse permanent dipole moment, which allows for a free rotation and then the common orientation under a strong electric field, yielding for the fluctuations $[eqn]$ and $[eqn]$ , and then $[eqn]$ , in agreement with Figure 5.

• For 5*CB, $[eqn]$ ; however, for fluctuations, only the component of the permanent dipole moment parallel to the director can be cancelled due to the prenematic arrangement, and a significant transverse component of the permanent dipole remains (See Figure 7). It can be freely oriented under the strong leading electric field, leading to $[eqn]$ and then $[eqn]$ , i.e., the negative pretransitional anomaly, in agreement with the results presented in Figure 6.
• For 5CB both $[eqn]$ and $[eqn]$ , leading $[eqn]$ , as evidenced in Figure 5.

For a sufficiently low scanning frequency, the averaged response from a few fluctuations is detected, thus eliminating the detection of disturbing impacts, such as pre-Smectic arrangements within fluctuations. It can explain the validity of Equation (13) for the I-SmA or I-SmE transition, for homologous series of n-isothiocyanatobiphenyls and n-thiocyanatobihenyls [106,107]. New experimental insights can facilitate the development of NDE devices with adjustable scanning frequencies.

3.6. Distortions from Landau–de Gennes Model Pattern Close to I-N Transition

A significant and still-inconclusive discussion concerns small deviations from the LdG temperature pattern (Equation (3)) that occur in the immediate vicinity of the I-N transition for KE, CME, or LS. They are particularly evident when testing reciprocals of the above properties, via a ‘bending down’ from the expected linear behaviour. It is heuristically explained as the influence of undefined ‘cybotactic groups’ [4,19,30,31,32,33,34,35,36,45,46,47,48,110,111].

For the authors, a possible explanation of this ‘inherent’ distortion is provided when analyzing the phenomenon via the model Equations (13) and (14). The extreme detection frequency for light–related methods, in comparison with the lifetime of fluctuations, causes the detection of fluctuations to be reduced to an apparent average of a single-molecule-related property. The crucial meaning is related to the molecular anisotropy of the refractive index. However, this is an average value seen over the entire uniaxial order within a prenematic fluctuation. The image of the prenematic fluctuation for 5CB in Figure 7 is an idealization. In reality, the molecules are slightly disoriented, and the director indicates only the dominant averaged prenematic direction. Uniaxial ordering improves when approaching the I-N transition, which is also related to the large increase in the number of molecules in the fluctuation. The mentioned methods (KE, CME, and LS), in fact, also detect (anisotropic) fluctuations. Therefore, it can be assumed that the detection using these methods is essentially related to the nondipolar component of dielectric polarization component in the direction indicated by the director $[eqn]$ . As $[eqn]$ approaches, this component increases on $[eqn]$ due to the improved ordering. The CME, OKE, and EKE amplitudes in Equation (3) also increase in comparison to the behaviour remote from the I-N transition. Consequently, the mentioned ‘bending down’ occurs.

The indirect evidence for this picture could be explained by the anomalous changes in the higher-frequency NDE pretransitional effect shown in the inset in Figure 5.

For NDE, it can be explained by the joint impacts of the measurement frequency and, first and foremost, the imperfect uniaxial ordering, changing the amplitude $[eqn]$ . On cooling towards $[eqn]$ , a registerable effective transverse component for large-enough fluctuations appears, which decreases $[eqn]$ , and starts the route towards the negative NDE (as for 5*CB). It is the ‘bending up’ that is visible in the inset in Figure 5. However, on approaching the singular temperature, the compressibility critically rises, which, under the strong electric field, improves uniaxiality. The initial trend diminishes, and the ‘bending down’ in the inset in Figure 5 appears. $[eqn]$ rises towards the value related to the perfect uniaxial arrangements.

4. Materials and Methods

Liquid crystalline compounds were purchased from Fluka (Seelze, Germany), and a few [freezing, vapour removal, melting and heating up to $[eqn]$ , and then freezing] were carried out. It degassed samples and removed parasitic contaminations. The mesomorphism of the tested compound is given in Section 2.

Broadband dielectric spectroscopy studies were carried out using a Novocontrol impedance analyzer, for frequencies from $[eqn]$ to $[eqn]$ , with the voltage of the measurement (scanning) electric field $[eqn]$ , which yields the optimal resolution of measurements. Samples were placed in a flat-parallel capacitor with the gap $[eqn]$ , which yields $[eqn]$ for the intensity of the scanning electric field. The capacitor was placed in the Quattro Novocontrol unit for temperature control, enabling stabilization of $[eqn]$ . Temperature was monitored using a Pt100 thermistor (UnipressEquipment, Warszawa, Poland), placed within one of the capacitor’s plates.

For the nonlinear dielectric effect (NDE), the single-resonant circuit concept was employed, where a strong electric field pulse is applied to the tested sample, causing a frequency shift. The latter was monitored via a modulation domain analyzer, allowing for frequency versus time scans. The design of the apparatus is given in refs. [66].

The strong electric field was applied in the form of ‘rectangular (DC)’ pulses lasting $[eqn]$ , with the voltage changing from $[eqn]$ to $[eqn]$ . For the applied capacitor, it led to electric field intensities ( $[eqn]$ ) from $[eqn]$ to $[eqn]$ . During the experiment, the condition related to Equation (4) $[eqn]$ was permanently tested, and this condition was tested during experiments.

During the experiment, the condition related to Equation (4) $[eqn]$ was permanently tested, and this condition was tested during experiments.

For the weak, radio-frequency-related scanning electric field, the intensity $[eqn]$ was applied. It was associated with the voltage $[eqn]$ applied to the measurement capacitor, which enables an optimal resolution. Note that NDE measurements are associated with extremely small dielectric constant changes $[eqn]$ , which show the technological challenges of this method [64,65,112,113].

The macroscopic gap of the measurement capacitor is notable, as it avoids surface-related effects that can occur when micrometric gaps are used. The latter are often employed in BDS studies. A notable feature of the NDE technique used in this study is the application of a ‘rectangular’ pulse of the DC electric field, which yields a comparable response from the activated sample. This enables any parasitic impacts that could distort the measurements to be identified during the experiment. These impacts can result from gas bubbles if the tested sample is not satisfactorily degassed or from solid contamination, such as dust particles. These factors all contribute to making the applied NDE measurement concept a reliable experimental solution. Another advantage of the NDE measurement concept is its sensitivity, enabling research in the range of $[eqn]$ to $[eqn]$ and the detection of statistical fluctuations in the polarizability of non-dipole liquids [66]. The only limitation is the lack of commercially available solutions despite the existence of proven advanced solutions, such as those in the laboratory of the authors of the given report [67].

5. Conclusions

The isotropic phase of rod-like nematic compounds is a unique experimental model system that provides a significant reference point for theoretical and simulation analyses in The Physics of Liquid Crystals and the broader category of Soft Matter systems. It is also important for Critical Phenomena Physics, due to the emergence of strong, long-range critical-like pretransitional effects despite the discontinuous nature of the I–N transition.

The Landau–de Gennes model (1969/1970) [1,27,28] plays a key role here, as it was inspired by pretransition effects exhibiting the parameterization shown in Equation (3). Numerous monographs and review reports have validated the Landau–de Gennes model by recalling experimental results related to this equation (Equation (3)). However, already in 1974, the model’s creator, Pierre-Gilles de Gennes, in his famous monograph, pointed out troubling problems with the LdG model for the I-N transition, which are recalled as questions Q1 and Q2 [1] in the Introduction. Over time, more such issues have emerged and are noted in questions Q3–Q6 in the Introduction. This report presents the first comprehensive response to all these perplexing problems.

Particular attention is paid to evidence associated with the ‘linear’ (dielectric constant) and ‘nonlinear’ (NDE) dielectric studies, where the manifestation of unusual and diverse features of pretransition effects in the isotropic phase is pervasive and crucial for explaining Q1–Q6 quests.

This work demonstrates that a fundamental understanding of various aspects of pretransitional effects, which manifest as apparent deviations from the mean-field description, is possible in the context of Critical Phenomena Physics [18,19,79]. However, it is essential to acknowledge factors such as the role of exceptional prenematic and pretransitional fluctuations, the impact of molecular structure on these fluctuations, and the fundamental importance of specific measurement methods, emphasizing the significance of the scanning timescale. This is where the model introduced in refs. [58,63] becomes particularly important, as demonstrated by the discussion of Equations (13) and (14). It enables a consistent description of pretransitional effects, including KE, LS, and NDE. For the latter, it includes the unusual negative NDE shown in Figure 6. It also provides a common interpretation of the pretransitional effects in the isotropic phase of nematogens and in the supercritical domains near the critical consolute and gas–liquid critical points. Finally, it answers question Q1, which particularly troubled de Gennes.

The report shows that it is possible to provide a coherent explanation for questions Q1–Q6, indicating puzzling problems related to the isotropic phase of nematogens.

A significant part of the report constitutes the results of focused dielectric constant and nonlinear dielectric effect studies in the isotropic phase of selected LC compounds. Notably, a similar, critical-like pattern has been evidenced for nanocolloidal systems composed of nCB liquid crystals and nanoparticles, such as C_60_ fullerenes and BaTiO_3_. These reports also cover LC mesophases and the solid phase, and in each case, the dominant impact of pretransitional fluctuations was demonstrated. These results are particularly important when considering the emerging applications of such materials in innovative renewable energy devices. The report does not address issues related to the dynamics of nematogens in the isotropic phase, which remain a significant challenge for model analysis within Landau–de Gennes theory. However, it provides an extremely important reference for glass transition physics. In this context, the universal critical-type characteristics of changes in the primary relaxation time are worth mentioning, as these describe both dynamics in the isotropic phase as it approaches the nematic phase and the glass temperature in supercooled liquids, as demonstrated by Drozd-Rzoska in ref. [26]:

$$[eqn]$$

where $[eqn]$ is the relative distance for the extrapolated singular temperature.

Numerous reports have discussed the implementation of the Landau–de Gennes model for liquid crystalline systems, particularly in the context of the isotropic liquid phase. The latter represents the initial consideration of the system in the history of Model I, which remains its hallmark reference. However, none of these reports have yet addressed a fundamental issue from the perspective of the scientific method [114]: coherently addressing the numerous problems that arise when comparing fused experimental evidence with the basic LdG model. This report is the first of its kind, demonstrating that many experimental phenomena that are surprising in the context of the LdG model can be coherently explained using the model. However, further challenges arise, e.g., related to the issue of dynamic phenomena in the isotropic phase, as indicated above.

The question arises as to whether the model explanations proposed and discussed in the given report fall within or beyond the framework of the Landau–de Gennes model. In our opinion, they are a crucial supplementation to the canonical Landau–de Gennes model and should be considered in any of its future developments, particularly in light of the clear experimental evidence.