Next Article in Journal
Total Synthesis of Nicrophorusamide A and Structural Disproof of the Proposed Noursamycin A
Next Article in Special Issue
Electrically Induced Structural Transformations of a Chiral Nematic under Tangential-Conical Boundary Conditions
Previous Article in Journal
Natural Products–Pyrazine Hybrids: A Review of Developments in Medicinal Chemistry
Previous Article in Special Issue
Multilayer, Broadband Infrared Reflectors Based on the Photoinduced Preparation of Cholesteric Liquid Crystal Polymers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Impact of Intermolecular Interactions on the Glassy Phase Formation of Twist-Bend Liquid Crystal Dimers: Insights from Dielectric Studies

1
Institute of Materials Engineering, Faculty of Science and Technology, University of Silesia, 75 Pułku Piechoty 1a, 41-500 Chorzów, Poland
2
Faculty of Electrical Enginesering, Automatics, Computer Science and Biomedical Engineering, AGH University of Science and Technology, 30-059 Krakow, Poland
3
Department of Applied Chemistry and Life Science, Graduate School of Engineering, Toyohashi University of Technology, Toyohashi 441-8580, Japan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2023, 28(21), 7441; https://doi.org/10.3390/molecules28217441
Submission received: 5 October 2023 / Revised: 30 October 2023 / Accepted: 31 October 2023 / Published: 6 November 2023
(This article belongs to the Special Issue Liquid Crystals II)

Abstract

:
The formation of the nematic to twist-bend nematic (NTB) phase has emerged as a fascinating phenomenon in the field of supramolecular chemistry, based on complex intermolecular interactions. Through a careful analysis of molecular structures and dynamics, we elucidate how these intermolecular interactions drive the complex twist-bend modulation observed in the NTB. The study employs broadband dielectric spectroscopy spanning frequencies from 10 to 2 × 109 Hz to investigate the molecular orientational dynamics within the glass-forming thioether-linked cyanobiphenyl liquid crystal dimers, namely, CBSC7SCB and CBSC7OCB. The experimental findings align with theoretical expectations, revealing the presence of two distinct relaxation processes contributing to the dielectric permittivity of these dimers. The low-frequency relaxation mode is attributed to an “end-over-end rotation” of the dipolar groups parallel to the director. The high-frequency relaxation mode is associated with precessional motions of the dipolar groups about the director. Various models are employed to describe the temperature-dependent behavior of the relaxation times for both modes. Particularly, the critical-like description via the dynamic scaling model seems to give not only quite good numerical fittings, but also provides a consistent physical picture of the orientational dynamics in accordance with findings from infrared (IR) spectroscopy. Here, as the longitudinal correlations of dipoles intensify, the m1 mode experiences a sudden upsurge in enthalpy, while the m2 mode undergoes continuous changes, displaying critical mode coupling behavior. Interestingly, both types of molecular motion exhibit a strong cooperative interplay within the lower temperature range of the NTB phase, evolving in tandem as the material’s temperature approaches the glass transition point. Consequently, both molecular motions converge to determine the glassy dynamics, characterized by a shared glass transition temperature, Tg.

Graphical Abstract

1. Introduction

The simplest liquid crystalline mesophase, the uniaxial nematic (N) phase, is characterized by the orientation of the molecules in a preferred direction, but with low positional order. The introduction of chirality, either through a chiral center in the mesogen or by adding a chiral dopant to the N phase of the host, results in the formation of a chiral nematic (N*) phase that has a helical structure. Independently, Meyer [1] and Dozov [2] predicted that liquid crystal phases with local chirality could be formed with bent achiral mesogens. The helical phase could occur spontaneously as a result of the simultaneous bending and twisting of a local director in an array of non-chiral molecules. The resulting twist-bend nematic (NTB) phase is present in an equal number of degenerate domains with opposite twist directions and the director is tilted relative to the helical axis. This new nematic-type phase is distinguished by unusual periodic patterns observed with polarizing microscopy and very fast electro-optic switching in the microsecond regime [3,4]. However, an unambiguous determination of the phase structure appears to be very difficult, since it exhibits no electron density modulation. Extensive studies [5,6,7] suggested that this structure corresponds to a phase with a spontaneous conical twist-bend director distortion, which was theoretically predicted to be driven by an elastic instability with the sign inversion of the bend elastic constant K33 for bent-shaped molecules [2].The NTB phase was first discovered for 1,7-bis-4-(4′-cyanobiphenyl)heptane, CB7CB, which consists of two mesogenic units connected with a flexible spacer [5,6,7,8]. Other dimeric mesogen [9,10,11,12,13,14,15,16,17] bent-core species [18,19,20] have been reported to manifest the NTB phase. Despite the significant number of materials investigated, a general structure–property relationship for the NTB phase remains elusive; however, it is found that a spatially uniform curvature of the molecule is necessary to form the NTB phase regardless of the underlying chemical groups [21,22].
The primary motivation for this work is the analysis of the orientational dynamics of molecules in the nematic phases: N and NTB formed with the polar symmetric dimers CBSC7SCB and the similar asymmetric dimers CBSC7OCB. Dielectric relaxation measurements were carried out for both samples in the entire temperature range up to the glass transition. The observed relaxation processes, as in earlier dielectric relaxation studies [5,20,23,24], were interpreted using different molecular motions defined with theoretical models of dielectric relaxation [25,26,27]. In both nematic phases, the rotation of dipolar groups associated with the terminal cyanobiphenyl groups lead to two relaxation modes related to the rotational dynamics of the molecules. The low-frequency mode, denoted m1, is caused by “end-to-end” motion of dipolar groups, which are parallel to the mesogen axis. The high-frequency mode, denoted m2, is the result of the precessional rotation of the dipolar groups around the director. These modes contribute to the complex dielectric permittivity differently depending on the orientation of the dielectric relative to the measured electric field.
The important part of the analysis concerns the study of the glass transition identified in both longer dimers clearly below the N-NTB transition. The glass transition can be characterized by the temperature Tg and the rate at which various properties of the material change with temperature as the liquid phase approaches the glass transition. For systems exhibiting a certain number of structural relaxations, it is possible to identify dielectric glass transitions corresponding to each of the two types of molecular motion independently. It is interesting whether changes in molecular interactions observed in previous spectroscopic studies [28,29] can be related to the observed glass transitions. The results for the CBSC7SCB and CBSC7OCB dimers presented here are discussed for information on which molecular movements are frozen during glass formation. The results for the symmetric and asymmetric liquid crystal dimers presented here are discussed to provide information on how the changes in molecular interactions observed near the N-NTB transition and below affect the aggregation process and how they are frozen during glass formation.

2. Results

2.1. Dielectric Relaxation

The simplest dielectric absorption model applicable to most materials is the Debye model, which describes the single dipole relaxation process. The temperature dependence of the relaxation time can be represented by the phenomenological Arrhenius equation, which introduces the activation energy Ea for the reorientation of a molecular dipole in a dielectric environment.
τ = τ 0   exp [   E a / k B T ] ,
Although the Arrhenius equation describes the dielectric relaxation of many simple fluids well, there are several materials for which the equation fails and it is necessary to modify the equations in describing experimental results. One of the phenomenological equations most commonly used to describe the temperature dependence of relaxation time data (τ) is the Vogel–Fulcher–Tammann (VFT) equation [30]:
τ = τ 0   exp [   B / T T 0 ] ,
Despite some concerns about the validity of the VFT equation, as it assumes dynamic divergence of the relaxation time at some finite temperature T0, some theories are based on the VFT equation such as the Adam–Gibbs entropy model (AG) [31] and some more recent theoretical models [32,33,34]. The basic idea behind these theories is that glass formation is associated with highly cooperative movements in a structurally fluctuating sample, with cooperativity growing as Tg is approached. Combining the above ideas of glass formation with a mean-field description of the virtual phase transition, the dynamic scaling (DS) model [35,36] has been proposed such that
τ = τ 0   [ T T c / T c ] Φ ,
where the pre-factor τ0 is defined as the relaxation time at 2Tc; the temperature Tc is the temperature of the virtual phase transition (the critical temperature), usually located slightly below Tg. For the high-temperature dynamic domain, well above the glass transition where the coupling mechanisms can be disregarded, Mode Coupling Theory (MCT) provides a power law function similar to the above eq. but with fitting parameters of the different physical meanings to those reported by the DS model. In that case, Tc accounts for the crossover temperature from the ergodic to the non-ergodic domain. The critical temperature Tc seems to correlate with the caging temperature TA [37].
An alternative non-divergence description that has recently been considered is the Waterton–Mauro (WM) parameterization, derived empirically by Waterton [38] in the 1930s and recently derived theoretically by Mauro and colleagues [39]; it is now considered a promising function for representing relaxation times, defined as follows:
τ = τ 0   exp [ K / T   exp C / T ] ,
where τ0 has the same meaning as in the VFT equation; K and C parameters are related with effective activation barriers, being both defined as thermal activation fitting parameters.
The problem is more complex in dielectric studies of liquid crystals, because of their macroscopic anisotropy and the presence of a nematic potential. The solution of the rotational diffusion equation for a rigid dipolar molecule in the presence of a nematic potential predicts that the frequency dependence of each component of the electric permittivity can be characterized by two exponential decays [25,26]. Rotation of the molecule about a short molecular axis in the presence of the nematic potential gives rise to the low-frequency relaxation mode m1 detected in the parallel component of the permittivity, which is classified as τ00 in terms of the spherical harmonics. Relaxation of the longitudinal component of the molecular dipole with precession about the director axis contributes to a high-frequency mode, m2, detected in the perpendicular component of the permittivity and is classified as τ10 in terms of the spherical harmonics. This is for a nematic potential assumed to be of the form
U / k B T = σ cos 2 θ ,
where σ is proportional to the uniaxial order parameter S (for S < 0.6), while approximate expressions for the relaxation times can be derived [26] as
σ = 3 2 S 5 π S / 1 S 2 ,
The nematic potential barrier parameter, σ, is defined as σ = q/RT where q is the height of the barrier separating two minima along the n direction.
τ 00 τ D = e σ 1 σ 2 σ δ / π 1 + σ + 2 σ 1 ,
τ 10 τ D = 1 S 1 + S / 2 ,
The relaxation time τD is that for rotational diffusion in the isotropic phase. Longitudinal, μl, and transverse, μt, components of the molecular dipole moment, μ, contribute to the dielectric permittivity differently and they relax at different frequencies of the probe field. In solving the rotational diffusion equation, a uniaxial nematic potential has been assumed. An inconsistency arises here, since molecules that have a dipole inclined to the long molecular axis are intrinsically biaxial.

2.2. Molecular Modes in the Nematic Phase

The dynamic dielectric responses of CBSCnSCB and CBSCnOCB (n = 5, 7) were measured in metal cells with a 50 μm spacer, in the frequency range of 5 Hz to 1 GHz [40]. The dielectric response of CBSCnSCB and CBSCnOCB in the nematic phase shows two relaxation processes, the contribution of which to the dielectric spectrum depends on the orientation of the director in the measurement cell. These dielectric relaxations can be related to the rotational diffusion of dimer molecules: the reorientation of end-to-end dipolar groups parallel to the director at low frequencies (m1), according to the theoretical model of dielectric relaxation in nematic dimers [25,26], and the precessional movement of dipolar groups around the director for the high-frequency branch of the spectrum (m2).
As already stated, the ε″ spectra at higher frequencies are dominated with two maxima. These correspond to the two molecular relaxation modes m1 and m2 of the symmetric CBSCnSCB and asymmetric CBSCnOCB dimers, as also made for CBCnCB by Cestari et al. [5], Merkel et al. [24], and López et al. [41]. We analyzed the two observed processes, in the N phase, and we were enabled to describe them unequivocally using an orientational order, S, and so-called Debye relaxation time, which corresponds to the relative relaxation in the absence of orientational order. This can be simply interpreted as an extension of the isotropic relaxation time into the temperature range of the nematic phase. Following Equation (4), the contribution to the perpendicular component of permittivity originates from the precession rotation of the transverse component of the dipole moment, i.e., rotation of the bow axis of the bent-core conformation contributes significantly to the perpendicular component.
This is the assignment of mode m2. The higher frequency relaxation process, m2, can arise from the rotations of a segment of the molecule, i.e., the internal rotation of each monomer with the spacer anchored involves the fluctuations of the cyanobiphenyl dipolar moment. As the length of the spacer in between the two mesogenes in the dimer is large enough, such an independent internal rotation of each monomer of the dimer is highly feasible. The temperature dependencies of δε and of the relaxation time τ 10 are suggestive of the precessional rotation of the longitudinal component of each tio-cyanobiphenyl (SCB) dipole moment around the director in a planar-aligned cell.
Relaxation times of the modes, however, seem to be well defined solely with the local orientation order in the nematic phase. All dimer dynamics behavior in the nematic phase is quite well reproduced with the molecular dynamic model (Figure 1a,b). The experimental relaxation times, τ00 and τ10, for mode 1 and 2, respectively, were well fitted using Equations (6)–(8), assuming S and τD as unknown values. Results are shown in Figure 1a,b. Temperature dependences of the orientational parameters correspond quite well to the results obtained with infrared spectroscopy [29], but their values are slightly smaller than referenced because dielectric ones are related to the S-CB or (O-CB) dipole moments but IR ones correspond to the molecular axis.
In the NTB phase, the temperature dependences of the τ10 of the m2 mode keep their trends without any step at the transition temperature, irrespective of if the director becomes tilted with respect to the symmetry axis. It seems this process is exclusively determined with the local orientational order (relative to the local director), which is continuously growing on decreasing temperature. This is not the case for end-over-end relaxation time, τ00, of the mode m1. First of all, relaxation time dependence clearly shows a kink at the transition temperature, which indicates an increased relaxation rate at the transition temperature. This is more likely due to critical fluctuations of the director in the vicinity of transition [42]. Then, in the NTB phase, the temperature trend of the τ00 relaxation time suddenly increases, with respect to the trend in the nematic phase. We can simply interpret this fact as an apparent growth of the potential barrier, q, separating two minima along the n direction, upon entering the NTB phase [26,38]. This finding corresponds well with an orientational correlation effect of the longitudinal dipoles ( g | | 1 ), which was revealed with IR spectroscopy [28,29].
Models for the interpretation of the low-frequency relaxation in liquid crystals are often based on a single particle relaxation process, but spectroscopic probes of molecular motion such as magnetic resonance, neutron scattering, and time-resolved fluorescence depolarization suggest that reorientation times for mesogens are of the order of 10−9 s to 10−10 s in isotropic, nematic, and disordered smectic phases. Thus, interpretation of dielectric relaxation processes at MHz or even kHz frequencies in terms of single molecule rotation is not likely to be correct. The low-frequency relaxations observed in liquid crystals are the result of collective molecular motion, although the models outlined above are useful in analyzing results and comparing materials.

3. Discussion

Dynamic Characterization on Glass Forming

We will now focus on the temperature range below the transition N-NTB. Only two longer dimers, with seven carbons in the link (i.e., CBSC7SCB and CBSC7OCB), show typical glass-forming behavior on cooling [40]. We used phenomenological equations and the Vogel–Fulcher–Tammann formula, Equation (2), to describe the temperature dependence of relaxation time data (τ). The results of fitting are shown in Figure 1a,b as dashed lines. VFT parameters obtained for Equation (2) are listed in Table 1. It should be stressed, as can be observed in Table 1, that two different dielectric glass transition temperatures are obtained, though rather close to one another, within about 1 K. It is also interesting to note that the pre-factor τ0 for the m2 mode is of the order of 10−11 s.
A more detailed method to study the dynamics of a glass-forming behavior is to analyze the derivatives of the relaxation time with respect to the inverted temperature [43,44]. Application of this procedure by Rzoska and Drozd-Rzoska leads to the relationship
d   l n τ d 1 / T   1 / 2 = H A T R   1 / 2 = B 1 / 2 1 T 0 / T ,
where HA is denoted as the apparent enthalpy of activation. In terms of the applicability of the VFT equation,   H A T / R   1 / 2 is predicted to have a linear dependence on inverse temperature. The results of applying this equation to both modes are shown in Figure 2a for CBSC7SCB and Figure 2b for CBSC7OCB.
In the case of CBSC7SCB, and mode m2, the activation enthalpy (almost constant in N phase, HA = 42 kJ/mol) begins to increase in the NTB phase showing a region of linearity (of H A T / R   1 / 2 ) down to 325 K (1000/T ≅ 3.08 K−1). On the contrary, the activation enthalpy of the mode m1 jumps at transition N-NTB (from 54 kJ/mol to 113 kJ/mol) then increases, but less rapidly, to join m2 below 325 K, and from this temperature, both modes have the same enthalpy of activation, probably due to cooperative behavior. Similarly, mode m2 for CBSC7OCB gradually increases activation enthalpy in the NTB phase (from 53 kJ/mol at transition temperature up to 116 kJ/mol, additionally showing linear region ( H A T / R   1 / 2 ) down to 330 K (1000/T ≅ 3.03 K−1)). In this case, mode m2 shows two linear domains: first from the transition temperature down to 330 K and the second below it. This coincides with the temperature when lateral interaction begins to be important [29].
The activation enthalpy, of the mode m1, again jumps at transition N-NTB (from 53 kJ/mol to 101 kJ/mol) then grows but less rapidly to join the m2 mode at about 320 K and since then, both mode have the same activation enthalpy (~116 kJ/mol) due to the cooperative behavior. Continuous lines in both figures are linear fits according to Equation (9).
We now consider the temperature-derivative procedure to be applied to the critical-like description through Equation (3). A linear dependence with temperature is obtained in the region of validity of the critical-like descriptions according to the relationship
T 2 d   l n τ d 1 / T   1 = H A T T 2 R   1 = Φ 1 T T c ,
The TC in the above MCT critical-like equation accounts for the crossover temperature from the ergodic to the non-ergodic domain. The critical temperature Tc seems to correlate with the caging temperature TA. Experimental data for both modes are shown in Figure 3 as a plot of [T2R/HA(T)] vs. temperature. It should be stressed that for both dimers, the data of the m1 mode exhibit a linear behavior (i.e., follow the DS model) over the entire temperature range of the NTB mesophase. The temperature TC is the temperature of the virtual phase transition (the critical temperature), usually located slightly below Tg. The exponent Φ ≈ 6–15 was suggested as universal for glass-forming polymers. In contrast, m2 data in the NTB mesophase clearly show two linear domains, one at low temperatures (DS model) from Tg up to about 320–324 K (denoted as TA in the figures); for both modes, they seem to be indistinguishable from behavior of the m1 mode and another domain at high temperatures (MCT-like description) from TA up to the vicinity of the NTB-N phase transition. The crossover temperature coincides well with IR spectroscopy observation for the temperature when lateral interaction increases their strength [29]. Below that temperature, it seems, both modes become coupled and can be described in the same way.
The subsequent linear fittings according to Equation (10) yield the values of the parameters for each mode and dynamic domain, namely, TC and the exponent Φ. The final fitting of the relaxation data according to Equation (3) is focused on the pre-factor τ0. All these derived parameters are listed in Table 3 and results of the fittings are drawn in Figure 3.
One of the most noticeable results in Table 3 is related to both the TC and TA temperatures corresponding to the m1 and m2 modes. It seems that a common critical temperature (≈273 K) for the virtual phase transition and also the same glass transition temperature of about 276 K are obtained for both modes. As for the exponent, in the DS domain, Φ is 7.0 (m1 mode) and 7.4 (m2 mode), typical values usually reported for spin-glass-like systems. In the MCT domain for the m2 mode, Φ is within the usual range of values.
The application of the temperature-derivative procedure to Equation (4) does not allow for a similar simple linearization methodology as in the previous eq. Instead, an enthalpy function is described with two uncorrelated variables (K and C) in the form
    d   l n τ d 1 / T   = H A T R   = K 1 + C / T e C / T
Now, we can obtain K and C parameters and, consequently, the glass transition temperatures for both relaxation modes. The fitting procedure is unable to fit to only one dynamic domain for each mode as observed in Figure 4, but it allows us to better distinguish the different dynamic domains. It should be stressed that according to Equation (10) and the results of the fits for the C and K parameters (see Table 4), ln[HA(T)/R] exhibits a nearly linear dependence with the inverse of temperature. As suggested recently [44], to proceed further with a linearized expression, Equation (11) can be written as
l n H A T R 1 + C / T   = l n K + C / T
Different crossover temperatures for each set of relaxation modes, one at about 300 K for m2 data and another at about 330 K for m1 data, can be observed. The former has not been appreciated for VFT or critical-like analyses of m2 data. However, the latter has already been observed in the VFT analysis of m1 data and apparently coincides with the estimation of the caging temperature, TA, shown with the m2 data. The final fitting of the relaxation time data according to Equation (4) leads to the pre-factor τ0. All the (WM)-fitting parameters, C, K, and τ0, for each dynamic domain and relaxation mode are listed in Table 4. The final fits are also drawn in Figure 4. If we focus on the results for τ0 and Tg (see Table 4), a strong parallelism with the VFT results can be observed. Again, τ0 for the m2 mode is of the order of 10−10 s (high-temperature dynamic domain), far from 10−14 s. As for the glass transition temperature, two different values are obtained, but less than ca. 3K apart.

4. Materials and Methods

In this study, we investigated both symmetrical and asymmetrical liquid crystal dimers featuring cyanobiphenyl (CB) mesogenic groups. The symmetrical dimers, denoted as CBSCnSCB, incorporate a thioether bridge (S-Cn-S) with varying linker lengths, specifically, either five or seven methylene groups (n = 5, 7). On the other hand, the asymmetrical dimers, referred to as CBSCnOCB (n = 5, 7), connect the mesogens to an alkyl chain with either five or seven methylene groups on one side via a thioether bridge and on the other side via an ether bridge. The molecular structures and typical textures in the nematic and twist-bend nematic phases for investigated compounds are illustrated in Figure 5, and for comprehensive details regarding the synthesis of thioether/ether dimers, we refer to the cited references [45,46].
Dielectric spectroscopy investigations were conducted over a broad frequency range spanning from 10 Hz to 1 GHz. These measurements were performed using two impedance/network analyzers, specifically the HP4I92A and HP4I95A, and encompassed samples aligned in both planar and homeotropic orientations. Diverse cell configurations with gap sizes ranging from 1.6 to 12 μm were employed. Homeotropic cells were achieved by using a commercial solution of the AL 60,702 polymer (JSR Korea, Gongju-City, Republic of Korea) and planar cells equipped with SE-130 polymer alignment layers from Nissan Chemical Industries, Ltd. For high-frequency dielectric measurements, we utilized gold-plated cells. Quality of alignment was tested using the ITO cells of the same geometry and alignment agents. Geometry of planar and homeotropic sandwich cells and molecular frame of reference are shown in Figure 6.
For better deconvolution of the modes in the dielectric spectra, we used permittivity derivatives vs. ln(f) [20,24].
d ε d ( ln f ) = d ε d ( ln ω ) = j = 1 n Re δ ε j α ( i ω τ j ) α [ 1 + ( i ω τ j ) α ] 2 ,
Here, δεj, τj, αj are the fitting parameters of the equation of the derivative of ε′. Figure 7 shows the fitting examples of the ′/d(lnf) at two temperatures (370 K, 350 K) and how they reproduce the experimental curve. The errors are less than 5% for frequency and 7% for strength of the mode. The shape parameters were found to be quite small (α < 0.08).

5. Conclusions

The dielectric relaxation spectra for symmetric CBSC7SCB and asymmetric CBSC7OCB dimers are resolved in terms of two relaxation processes, m1 and m2, as predicted with the rotational diffusion model. The high-frequency relaxation process originates from the precession rotation of the longitudinal component of the monomer dipole moment of each thiocyanobiphenyl (SCB) dipole moment around the director. The relaxation rate is accelerated with respect to the isotropic relaxation because τ 10 is shorter as compared to τ D . The low-frequency relaxation corresponds from the “end-over-end” rotation of the monomer dipole moment. The relaxation rate of that mode is retarded with respect to the isotropic one. It is clear that dynamics of the two relaxation process modes m1 and m2 are quite different in the two temperature regions of the NTB mesophase. In the first region, which is extended almost 30 K from N-NTB transition, for mode m1 assigned to “end-over-end” rotation, activation enthalpy initially jumps almost twice on entering NTB then increases continuously over almost 30 K. This turns out to be well associated with the appearance of longitudinal dipole correlations discovered in IR studies. In the same region, the activation enthalpy of mode m2 (assigned to precession of dipoles) grows continuously but much faster than m1 and finally they meet at the same value of HA(T). This meeting temperature corresponds well to the appearance of transverse interactions as seen in IR measurements. The description of the critical behavior of dynamic domains related to cooperative molecular motions is quite sufficient to adequately describe the dynamics associated with the two main molecular motions in the NTB mesophase of the symmetric and asymmetric dimer. Both types of molecular motion appear to strongly cooperate at low temperatures, changing in a coordinated manner as the temperature of the material approaches the glass transition point. As expected, it turns out that both molecular motions determine glassy dynamics with the same glass transition temperature Tg.

Author Contributions

Conceptualization: A.K.; Synthesis: Y.A.; Methodology: A.K. and K.M.; Investigation: A.K., K.M. and M.C.; Formal Analysis: A.K. and K.M.; Writing—review and editing: A.K. and K.M.; Visualization: A.K. and K.M. The manuscript was written through contributions of all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Science Centre, Poland, for Grant Nos. 2018/31/B/ST3/03609 and 2020/39/O/ST5/03460.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

RepOD will be the main platform for data storage after the project will be finished. https://repod.icm.edu.pl/dataset.xhtml?persistentId=doi%3A10.18150%2FRXRLON (accessed on 1 January 2024).

Acknowledgments

Authors A.K. and M.C. thank the National Science Centre for funding through grant no. 2018/31/B/ST3/03609. K.M. thanks the National Science Centre for funding through Project No. 2020/39/O/ST5/03460. All DFT calculations were conducted with the Gaussian09 program using the PL-Grid Infrastructure on the ZEUS and Prometheus clusters in Academic Computer Center CYFRONET AGH (AGH University of Science and Technology) in Cracow, Poland.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Meyer, R.B. In Les Houches Summer School in Theoretical Physics; Balian, R.G., Weil, G., Eds.; Gordon and Breach: New York, NY, USA, 1976; pp. 273–373. [Google Scholar]
  2. Dozov, I. On the Spontaneous Symmetry Breaking in the Mesophases of Achiral Banana-Shaped Molecules. Europhys. Lett. 2001, 56, 247–253. [Google Scholar] [CrossRef]
  3. Panov, V.P.; Balachandran, R.; Nagaraj, M.; Vij, J.K.; Tamba, M.G.; Kohlmeier, A.; Mehl, G.H. Microsecond Linear Optical Response in the Unusual Nematic Phase of Achiral Bimesogens. Appl. Phys. Lett. 2011, 99, 261903-3. [Google Scholar] [CrossRef]
  4. Panov, V.; Balachandran, R.; Vij, J.K.; Tamba, M.G.; Kohlmeier, A.; Mehl, G.H. Field-Induced Periodic Chiral Pattern in the Nx Phase of Achiral Bimesogens. Appl. Phys. Lett. 2012, 101, 234106-4. [Google Scholar] [CrossRef]
  5. Cestari, M.; Diez-Berart, S.; Dunmur, D.A.; Ferrarini, A.; de la Fuente, M.R.; Jackson, D.J.B.; Lopez, D.O.; Luckhurst, G.R.; Perez-Jubindo, M.A.; Richardson, R.M.; et al. Phase Behavior and Properties of the Liquid Crystal Dimer 1″,7″-bis(4-cyanobiphenyl-4′-yl) heptane: A Twist-Bend Nematic Liquid Crystal. Phys. Rev. E 2011, 84, 031704-20. [Google Scholar] [CrossRef] [PubMed]
  6. Chen, D.; Porada, J.H.; Hooper, J.B.; Klittnick, A.; Shen, Y.; Tuchband, M.R.; Korblova, E.; Bedrov, D.; Walba, D.M.; Glaser, M.A.; et al. Chiral Heliconical Ground State of Nanoscale Pitch in a Nematic Liquid Crystal of Achiral Molecular Dimers. Proc. Natl. Acad. Sci. USA 2013, 110, 15931–15936. [Google Scholar] [CrossRef]
  7. Borshch, V.; Kim, Y.K.; Xiang, J.; Gao, M.; Jákli, A.; Panov, V.P.; Vij, J.K.; Imrie, C.T.; Tamba, M.G.; Mehl, G.H.; et al. Nematic Twist-Bend Phase with Nanoscale Modulation of Molecular Orientation. Nat. Commun. 2013, 4, 2635–2638. [Google Scholar] [CrossRef]
  8. Meyer, C.; Luckhurst, G.; Dozov, I. Flexoelectrically Driven Electroclinic Effect in the Twist-Bend Nematic Phase of Achiral Molecules with Bent Shapes. Phys. Rev. Lett. 2013, 111, 067801-5. [Google Scholar] [CrossRef]
  9. Panov, V.P.; Nagaraj, M.; Vij, J.K.; Panarin, Y.P.; Kohlmeier, A.; Tamba, M.G.; Lewis, R.A.; Mehl, G.H. Spontaneous Periodic Deformations in Nonchiral Planar-Aligned Bimesogens With a Nematic-Nematic Transition and a Negative Elastic Constant. Phys. Rev. Lett. 2010, 105, 167801-4. [Google Scholar] [CrossRef]
  10. Henderson, P.A.; Imrie, C.T. Methylene Linked Liquid Crystal Dimers and the Twist-Bend Nematic Phase. Liq. Cryst. 2011, 38, 1407–1414. [Google Scholar] [CrossRef]
  11. Adlem, K.; Čopič, M.; Luckhurst, G.R.; Mertelj, A.; Parri, O.; Richardson, R.M.; Snow, B.D.; Timimi, B.A.; Tuffin, R.P.; Wilkes, D. Chemically Induced Twist-Bend Nematic Liquid Crystals, Liquid Crystal Dimers, and Negative Elastic Constants. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2013, 88, 022503-8. [Google Scholar] [CrossRef]
  12. Mandle, R.J.; Davis, E.J.; Lobato, S.A.; Vol, C.-C.A.; Cowling, S.J.; Goodby, W.J. Synthesis and Characterisation of an Unsymmetrical, Ether-Linked, Fluorinated Bimesogen Exhibiting a New Polymorphism Containing the NTB or ‘twist-bend’ Phase. Phys. Chem. Chem. Phys. 2014, 16, 6907–6915. [Google Scholar] [CrossRef] [PubMed]
  13. Mandle, R.J.; Davis, E.J.; Archbold, C.T.; Cowling, S.J.; Goodby, W.J. Microscopy Studies of the Nematic NTB Phase of 1,11-di-(1″-cyanobiphenyl-4-yl)undecane. J. Mater. Chem. C 2014, 2, 556–566. [Google Scholar] [CrossRef]
  14. Mandle, R.J.; Davis, E.J.; Archbold, C.T.; Vol, C.-C.A.; Andrews, J.L.; Cowling, S.J.; Goodby, W.J. Apolar Bimesogens and the Incidence of the Twist–Bend Nematic Phase. Chem. Eur. J. 2015, 21, 8158–8167. [Google Scholar] [CrossRef] [PubMed]
  15. Gorecka, E.; Vaupotic, N.; Zep, A.; Pociecha, D.; Yoshioka, J.; Yamamoto, J.; Takezoe, H. A Twist-Bend Nematic (NTB) Phase of Chiral Materials. Angew. Chem. Int. Ed. 2015, 54, 10155–10159. [Google Scholar] [CrossRef] [PubMed]
  16. Paterson, D.A.; Gao, M.; Kim, Y.-K.; Jamali, A.; Finley, K.L.; Robles-Hernández, B.; Diez-Berart, S.; Salud, J.; de la Fuente, M.R.; Timimi, B.A.; et al. Understanding the Twist-Bend Nematic Phase: The Characterisation of 1-(4-cyanobiphenyl-4′-yloxy)-6-(4-cyanobiphenyl-4′-yl)hexane (CB6OCB) and Comparison with CB7CB. Soft Matter 2016, 12, 6827–6840. [Google Scholar] [CrossRef]
  17. Paterson, D.A.; Abberley, J.P.; Harrison, W.T.A.; Storey, M.D.J.; Imrie, C.T. Cyanobiphenyl-Based Liquid Crystal Dimers and the Twist-Bend Nematic Phase. Liq. Cryst. 2017, 44, 127–146. [Google Scholar] [CrossRef]
  18. Chen, D.; Nakata, M.; Shao, R.; Tuchband, M.R.; Shuai, M.; Baumeister, U.; Weissflog, W.; Walba, D.M.; Glaser, M.A.; Maclennan, J.E.; et al. Twist-Bend Heliconical Chiral Nematic Liquid Crystal Phase of an Achiral Rigid Bent-Core Mesogen. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2014, 89, 022506-5. [Google Scholar] [CrossRef]
  19. Sreenilayam, S.P.; Panov, V.P.; Vij, J.K.; Shanker, G. The NTB Phase in an Achiral Asymmetrical Bent-Core Liquid Crystal Terminated with Symmetric Alkyl Chains. Liq. Cryst. 2017, 44, 244–253. [Google Scholar] [CrossRef]
  20. Merkel, K.; Kocot, A.; Vij, J.K.; Shanker, G. Distortions in Structures of the Twist Bend Nematic Phase of a Bent-Core Liquid Crystal by the Electric Field. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2018, 98, 022704-8. [Google Scholar] [CrossRef]
  21. Walker, R.; Pociecha, D.; Salamończyk, M.; Storey, J.M.D.; Gorecka, E.; Imrie, C.T. Supramolecular Liquid Crystals Exhibiting a Chiral Twist-Bend Nematic Phase. Mater. Adv. 2020, 1, 1622–1630. [Google Scholar] [CrossRef]
  22. Mandle, R.J. The Dependency of Twist-Bend Nematic Liquid Crystals on Molecular Structure: A Progression from Dimers to Trimers, Oligomers and Polymers. Soft Matter 2016, 12, 7883–7901. [Google Scholar] [CrossRef] [PubMed]
  23. Merkel, K.; Welch, C.; Ahmed, Z.; Piecek, W.; Mehl, G.H. Dielectric Response of Electric-Field Distortions of the Twist-Bend Nematic Phase for LC Dimers. J. Chem. Phys. 2019, 151, 114908–114920. [Google Scholar] [CrossRef] [PubMed]
  24. Merkel, K.; Kocot, A.; Welch, C.; Mehl, G.H. Soft Modes of the Dielectric Response in the Twist–Bend Nematic Phase and Identification of the Transition to a Nematic Splay Bend Phase in the CBC7CB Dimer. Phys. Chem. Chem. Phys. 2019, 21, 22839–22848. [Google Scholar] [CrossRef] [PubMed]
  25. Nordio, P.L.; Rigatti, G.; Segre, U. Dielectric Relaxation Theory in Nematic Liquids. Mol. Phys. 1973, 25, 129–136. [Google Scholar] [CrossRef]
  26. Coffey, W.T.; Kalmykov, Y.P. Rotational Diffusion and Dielectric Relaxation in Nematic Liquid Crystals. Adv. Chem. Phys. 2000, 113, 487–551. [Google Scholar]
  27. Stocchero, M.; Ferrarini, A.; Moro, G.J.; Dunmur, D.A.; Luckhurst, G.R. Molecular Theory of Dielectric Relaxation in Nematic Dimers. J. Chem. Phys. 2004, 121, 8079–8097. [Google Scholar] [CrossRef]
  28. Kocot, A.; Loska, B.; Arakawa, Y.; Merkel, K. Structure of the Twist-Bend Nematic Phase with Respect to the Orientational Molecular Order of the Thioether-Linked Dimers. Phys. Rev. E 2022, 105, 044701-10. [Google Scholar] [CrossRef]
  29. Merkel, K.; Loska, B.; Arakawa, Y.; Mehl, G.H.; Karcz, J.; Kocot, A. How Do Intermolecular Interactions Evolve at the Nematic to Twist–Bent Phase Transition? Int. J. Mol. Sci. 2022, 23, 11018. [Google Scholar] [CrossRef]
  30. García Colin, L.S.; del Castillo, L.F.; Goldstein, P. Theoretical Basis for the Vogel-Fulcher-Tammann Equation. Phys. Rev. B 1989, 40, 7040–7044, Erratum in Phys. Rev. B 1990, 41, 4785. [Google Scholar] [CrossRef]
  31. Adams, G.; Gibbs, J.H. On the Temperature Dependence of Cooperative Relaxation Properties in GlassForming Liquids. J. Chem. Phys. 1965, 43, 139–146. [Google Scholar] [CrossRef]
  32. Cohen, M.H.; Grest, G.S. Liquid-glass transition, a free-volume approach. Phys. Rev. B 1979, 20, 1077–1098. [Google Scholar] [CrossRef]
  33. Bouchaud, G.P.; Biroli, G. On the Adam-Gibbs-Kirkpatrick-Thirumalai-Wolynes Scenario for the Viscosity Increase in Glasses. J. Chem. Phys. 2004, 121, 7347–7354. [Google Scholar] [CrossRef] [PubMed]
  34. Lubchenko, V.; Wolynes, P.G. Theory of Structural Glasses and Supercooled Liquids. Ann. Rev. Phys. Chem. 2007, 58, 235–266. [Google Scholar] [CrossRef] [PubMed]
  35. Colby, R.H. Dynamic Scaling Approach to Glass Formation. Phys. Rev. E 2000, 61, 1783–1792. [Google Scholar] [CrossRef] [PubMed]
  36. Erwin, B.M.; Colby, R.H. Temperature Dependence of Relaxation Times and the Length Scale of Cooperative Motion for Glass-Forming Liquids. J. Non-Cryst. Solids 2002, 307–310, 225–231. [Google Scholar] [CrossRef]
  37. Drozd-Rzoska, A.; Rzoska, S.J.; Pawlus, S.; Martínez-García, J.C.; Tamarit, J.-L. Evidence for Critical-Like Behavior in Ultraslowing Glass Forming Systems. Phys. Rev. E 2010, 82, 031501-8. [Google Scholar] [CrossRef]
  38. Waterton, S.C. The Viscosity-Temperature Relationship and Some Inferences on the Nature of Molten and of Plastic. J. Soc. Glass. Technol. 1932, 16, 244–249. [Google Scholar]
  39. Mauro, J.C.; Yue, Y.; Ellison, A.J.; Gupta, P.K.; Allan, D.C. Viscosity of Glass-Forming Liquids. Proc. Natl. Acad. Sci. USA 2009, 106, 19780–19784. [Google Scholar] [CrossRef]
  40. Kocot, A.; Czarnecka, M.; Arakawa, Y.; Merkel, K. Dielectric Study of Liquid Crystal Dimers: Probing the Orientational Order and Molecular Interactions in Nematic and Twist-Bend Nematic Phases. Phys. Chem. B 2023, 127, 7082–7090. [Google Scholar] [CrossRef]
  41. López, D.O.; Sebastian, N.; de la Fuente, M.R.; Martínez-García, J.C.; Salud, J.; Pérez-Jubindo, M.A.; Diez-Berart, S.; Dunmur, D.A.; Luckhurst, G.R. Disentangling Molecular Motions Involved in the Glass Transition of a Twist-Bend Nematic Liquid Crystal Through Dielectric Studies. J. Chem. Phys. 2012, 137, 034502-10. [Google Scholar] [CrossRef]
  42. Pociecha, D.; Crawford, C.A.; Paterson, D.A.; Storey, J.M.D.; Imrie, C.T.; Vaupotič, N.; Gorecka, E. Critical Behavior of the Optical Birefringence at the Nematic to Twist-Bend Nematic Phase Transition. Phys. Rev. E 2018, 98, 052706-5. [Google Scholar] [CrossRef]
  43. Stickel, F.; Fischer, E.W.; Richert, R. Dynamics of Glass-Forming Liquids. I. Temperature-Derivative Analysis of Dielectric Relaxation Data. J. Chem. Phys. 1995, 102, 6251–6257. [Google Scholar] [CrossRef]
  44. Martínez-García, J.C.; Tamarit, J.L.; Rzoska, S.J. Enthalpy Space Analysis of the Evolution of the Primary Relaxation Time in Ultraslowing Systems. J. Chem. Phys. 2011, 134, 024512-7. [Google Scholar] [CrossRef]
  45. Arakawa, Y.; Komatsu, K.; Tsuji, H. Twist-Bend Nematic Liquid Crystals Based on Thioether Linkage. New J. Chem. 2019, 43, 6786–6793. [Google Scholar] [CrossRef]
  46. Arakawa, Y.; Ishida, Y.; Tsuji, H. Ether- and Thioether-Linked Naphthalene-Based Liquid-Crystal Dimers: Influence of Chalcogen Linkage and Mesogenic-Arm Symmetry on the Incidence and Stability of the Twist–Bend Nematic Phase. Chem. Eur. J. 2020, 26, 3767–3775. [Google Scholar] [CrossRef]
Figure 1. Plots of relaxation times for modes m1 and m2 (5 μm planar cell): (a) symmetric dimer (CBSC7SCB), (b) asymmetric dimer (CBSC7OCB). Symbols: red circles ()—m1 mode, dark blue squares ()—m2 mode, rotational diffusion model fitting: black solid line—m1 mode, dash blue line—τD relaxation time, solid blue line—VFT model [30], solid red line—fitting of experimental data by Equation (2), purple/burgundy triangles (, )—order parameter, S (as an insert).
Figure 1. Plots of relaxation times for modes m1 and m2 (5 μm planar cell): (a) symmetric dimer (CBSC7SCB), (b) asymmetric dimer (CBSC7OCB). Symbols: red circles ()—m1 mode, dark blue squares ()—m2 mode, rotational diffusion model fitting: black solid line—m1 mode, dash blue line—τD relaxation time, solid blue line—VFT model [30], solid red line—fitting of experimental data by Equation (2), purple/burgundy triangles (, )—order parameter, S (as an insert).
Molecules 28 07441 g001
Figure 2. Plots show a linear dependence H A T / R   1 / 2 on the inverse temperature in NTB phase. (a)—CBSC7SCB dimer, (b)—CBSC7OCB. The red arrows show the jump of the HA that omits the critical fluctuation at the N-NTB transition, where HA is denoted as the apparent enthalpy of the activation. Symbols: red circles ()—m1 mode, dark blue squares ()—m2 mode, solid line: model fitting according to Equation (9). The fitting parameters are listed in Table 2.
Figure 2. Plots show a linear dependence H A T / R   1 / 2 on the inverse temperature in NTB phase. (a)—CBSC7SCB dimer, (b)—CBSC7OCB. The red arrows show the jump of the HA that omits the critical fluctuation at the N-NTB transition, where HA is denoted as the apparent enthalpy of the activation. Symbols: red circles ()—m1 mode, dark blue squares ()—m2 mode, solid line: model fitting according to Equation (9). The fitting parameters are listed in Table 2.
Molecules 28 07441 g002
Figure 3. Results of the temperature-derivative analysis (Equation (10)) in NTB phase applied to both m1 and m2 modes in which linear dependences indicate domains of validity of the critical-like description; (a) symmetric dimer (CBSC7SCB), (b) asymmetric dimer (CBSC7OCB). Symbols: red circles ()—m1 mode, dark blue squares ()—m2 mode, solid lines: model fitting according to Equation (10). The fitting parameters are listed in Table 3.
Figure 3. Results of the temperature-derivative analysis (Equation (10)) in NTB phase applied to both m1 and m2 modes in which linear dependences indicate domains of validity of the critical-like description; (a) symmetric dimer (CBSC7SCB), (b) asymmetric dimer (CBSC7OCB). Symbols: red circles ()—m1 mode, dark blue squares ()—m2 mode, solid lines: model fitting according to Equation (10). The fitting parameters are listed in Table 3.
Molecules 28 07441 g003
Figure 4. Results of the temperature-derivative analysis (Equation (11)) in NTB phase applied to both m1 and m2 modes of (a) CBSC7SCB and (b) CBSC7OCB. Dashed lines (red and black) indicate dynamic domains of validity of the critical-like description. Symbols: red circles ()—m1 mode, black squares ()—m2 mode.
Figure 4. Results of the temperature-derivative analysis (Equation (11)) in NTB phase applied to both m1 and m2 modes of (a) CBSC7SCB and (b) CBSC7OCB. Dashed lines (red and black) indicate dynamic domains of validity of the critical-like description. Symbols: red circles ()—m1 mode, black squares ()—m2 mode.
Molecules 28 07441 g004
Figure 5. (a) Molecular structure of symmetric (CBSCnSCB) and asymmetric (CBSCnOCB) dimers. (b) POM textures of the 5 μm planar cell for the CBSC7OCB sample in the N phase (375 K) and the NTB phase (355 K) (cooling rate of 5 °C/min). The scale bar equals 50 μm.
Figure 5. (a) Molecular structure of symmetric (CBSCnSCB) and asymmetric (CBSCnOCB) dimers. (b) POM textures of the 5 μm planar cell for the CBSC7OCB sample in the N phase (375 K) and the NTB phase (355 K) (cooling rate of 5 °C/min). The scale bar equals 50 μm.
Molecules 28 07441 g005
Figure 6. Laboratory and molecular frame of reference. Planar (top) and homeotropic cell (bottom) (a) in the N phase, (b) in the NTB phase, where k is the helical axis, (c) molecular frame of reference.
Figure 6. Laboratory and molecular frame of reference. Planar (top) and homeotropic cell (bottom) (a) in the N phase, (b) in the NTB phase, where k is the helical axis, (c) molecular frame of reference.
Molecules 28 07441 g006
Figure 7. Fitting of the derivative of permittivity vs. ln(f) showing the deconvolution of the molecular modes: m1 and m2.
Figure 7. Fitting of the derivative of permittivity vs. ln(f) showing the deconvolution of the molecular modes: m1 and m2.
Molecules 28 07441 g007
Table 1. Fitting parameters according to Equation (2) for the different dimers and the calculated glass transition temperature for the m1 and m2 modes.
Table 1. Fitting parameters according to Equation (2) for the different dimers and the calculated glass transition temperature for the m1 and m2 modes.
SampleModelog10 [τ(s)]T0 (K)B (K)Tgl (K)
CBSC7SCBm1−10.1266.6682.6289.8
m2−11.9265.6631.9286.3
CBSC7OCBm1−9.2261.3596.0284.3
m2−11.5260.5621.7280.5
Table 2. Fitting parameters according to Equations (2) and (9) for the different dimers.
Table 2. Fitting parameters according to Equations (2) and (9) for the different dimers.
Samplelog10 [τ(s)]Equation (2)Equation (9)Range
[1000/T(K−1)]
CBSC7SCB−10.32692892.96–3.3
−11.62642893.15–3.3
CBSC7OCB−8.92652663.08–3.3
−11.52612703.08–3.3
Table 3. Fitting parameters according to Equation (10) for the different dimers.
Table 3. Fitting parameters according to Equation (10) for the different dimers.
SampleTAΦ (T < TA)TcΦRange (K)Description
CBSC7SCB324 K287.4 7.32937.3305–347DS
3211.3330–351MCT
CBSC7OCB320 K282.5 7.0287.47.0295–350DS
303.51.3320–365MCT
Table 4. Fitting parameters according to Equation (11) for the different dimers.
Table 4. Fitting parameters according to Equation (11) for the different dimers.
Samplelog10 [τ(s)]K (K)C (K)Tgl (K)Range
(1000/T (K−1))
CBSC7SCB−7.780.07893016268.22.97–3.28
−9.080.00823653267.93.02–3.22
CBSC7OCB−7.150.0583028263.93.03–3.30
−9.050.0313198260.62.98–3.35
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kocot, A.; Czarnecka, M.; Arakawa, Y.; Merkel, K. Exploring the Impact of Intermolecular Interactions on the Glassy Phase Formation of Twist-Bend Liquid Crystal Dimers: Insights from Dielectric Studies. Molecules 2023, 28, 7441. https://doi.org/10.3390/molecules28217441

AMA Style

Kocot A, Czarnecka M, Arakawa Y, Merkel K. Exploring the Impact of Intermolecular Interactions on the Glassy Phase Formation of Twist-Bend Liquid Crystal Dimers: Insights from Dielectric Studies. Molecules. 2023; 28(21):7441. https://doi.org/10.3390/molecules28217441

Chicago/Turabian Style

Kocot, Antoni, Małgorzata Czarnecka, Yuki Arakawa, and Katarzyna Merkel. 2023. "Exploring the Impact of Intermolecular Interactions on the Glassy Phase Formation of Twist-Bend Liquid Crystal Dimers: Insights from Dielectric Studies" Molecules 28, no. 21: 7441. https://doi.org/10.3390/molecules28217441

APA Style

Kocot, A., Czarnecka, M., Arakawa, Y., & Merkel, K. (2023). Exploring the Impact of Intermolecular Interactions on the Glassy Phase Formation of Twist-Bend Liquid Crystal Dimers: Insights from Dielectric Studies. Molecules, 28(21), 7441. https://doi.org/10.3390/molecules28217441

Article Metrics

Back to TopTop