2.1. Helix Pitch and Cone Angle Measurement
The phase behavior of binary mixtures of CB7CB and 5CB was characterized as a function of weight percent 5CB,
x, in the range 12.5 <
x < 95, generating the phase diagram of
Figure 2a (obtained using differential scanning calorimetry, shown in
Supplementary Figure S1). The N–TB phase transition was probed using polarized transmission optical microscopy (PTOM) to observe the various mixtures in untreated and planar cells for
x up to 62.5 (
Supplementary Figure S2). In the mixtures with
x ≤ 25, the TB phase exhibits the typical stripe texture on cooling, which was shown to be due to a spontaneous undulation with displacement along
z of the planes of constant azimuthal orientation
φ [
42] caused by dilative stress on the helical structure due to the shrinking of
pH on cooling [
30]. For
x ≥ 37.5, we found that the stripes form briefly on cooling but that the texture then relaxes into a state of uniform birefringence, indicating that the addition of 5CB fluidizes the phase so that it can anneal into an undulation-free director field configuration [
43] (see
Supplementary Figure S3). This demonstrates the utility of making mixtures of 5CB with CB7CB to obtain a uniform and undulation-free, well-aligned TB phase. The isotropic (I) to nematic (N) transition temperature of the mixtures,
TIN, decreases monotonically as
x increases, while the N–TB phase transition temperature,
TNTB, decreases nearly linearly with the addition of 5CB until the transition becomes undetectable for
x > 62.5 (
Supplementary Figures S1 and S4). These behaviors are similar to those seen in mixtures of CB9CB/5CB [
19] and CB7CB/7OCB [
44].
Resonant soft X-ray scattering (RSoXS) at the carbon K-edge, used to carry out direct, in situ measurement of the bulk heliconical nematic structure in the mixtures [
30], exhibits diffraction arcs in the TB phase corresponding to periodicities in the 8–11 nm range. Extensive freeze-fracture studies (see
Supplementary Materials) also show the helix pitch
pH, the distance along
z for a 2π circuit of the helix, to be in the 8–11 nm range [
29]. These observations, discussed further below, lead us to assign the RSoXS diffraction to be from the fundamental periodicity of the helix, i.e., with
qH = 2π/
pH. In the mixtures, the diffraction arcs are smooth and are relatively narrow in wavevector
q at low temperatures, as in neat CB7CB [
30], indicating that the pitch is homogenous throughout the sample. As the temperature is raised, the scattering tends to broaden in
q into a distribution of individual arcs of differing
qH, including some with linewidths comparable to 0.0002 Å
−1, the wavevector resolution of the diffractometer (
Supplementary Figures S5 and S6, and [
30]). This behavior is taken as evidence for the development of domains in the sample with a distribution of values of average heliconical pitch. Detailed studies of neat CB7CB [
30] showed that the lower limit of the pitch in the measured distributions was a repeatable function of
T but that the upper limit varied erratically from scan to scan in
T, indicating that the TB helix pitch is much softer in response to stretching than to compression and that the stretching is due to non-uniform stress distributions that develop in the macroscopic textures of the helix axis [
30]. For this reason,
pH is taken to be the length corresponding to the lowest-
q half-height value in each of the measured distributions, as shown in
Supplementary Figure S6. The resulting
pH data are plotted vs. reduced temperature
TNTB −
T in
Figure 2b and vs.
T in
Supplementary Figure S7. As the TB phase of the
x = 0 and 12.5 mixtures is heated toward
TNTB, a coexistence range of
T is entered in which some of the TB domains melt. This causes the
pH distribution to narrow at the highest
T’s in the TB phase, with the upper limits of the pitch approaching the lower. However, the
x = 25 and 37.5 mixtures do not exhibit this behavior: the pitch range remains broad near the transition, indicating that the coexistence range is narrower at higher 5CB concentrations, consistent with what we observe using PTOM.
The optical cone angle of the TB helix,
θH, was determined for the mixtures from measurements of the birefringence, Δ
n, of the N and TB phases, as detailed in the
Supplementary Information (Supplementary Figure S8) [
45]. Published values of
θH from birefringence [
45] and NMR experiments [
46] are also available for neat CB7CB. The
θH data are plotted vs.
TNTB −
T in
Figure 2c and in
Supplementary Figure S9. The data show that Δ
n increases continuously on cooling through the N phase, and then near the N–TB transition abruptly begins to decrease. We take this change to indicate the onset of the collective heliconical ordering in the TB phase. The N phase birefringence is somewhat smaller than that of pure 5CB, presumably a consequence of the cyanobiphenyl groups in CB7CB being substantially tilted (by ~30°,
Supplementary Table S1) from the average dimer orientation in the N phase, which is substantial. The development of heliconical ordering then further reduces Δ
n in the TB phase.
2.2. Bend Deformation and the Geometry of the Twist-Bend Helix
If the heliconical ground state axis of the TB phase is taken to be along
z, then
n(
r) may be written as
where
φ(
z) is the azimuthal angle given by
as sketched in
Figure 1c–f. The heliconical structure can be represented as a rotation on a cone as in
Figure 1c or by the green director contour lines in
Figure 1e,f, representing the path (contour line) along which the incremental displacement is always along
n(
z). The local nematic order tensor is biaxial, with principal axes given by the director (
n), polarization (
p), and auxiliary (
m) unit vectors. We begin by considering the bend deformation of
n, given generally by the director rotation vector
where
s is the displacement along the contour. Since by definition of the pitch, we have
and, from the geometry of the helix in
Figure 1f,
we find
Figure 1f shows the geometry of reorientation on the director contour line of
Bn(
r), the bend rotation of
n(
r) about auxiliary vector
m(
r), and
Tp(
r), the twist rotation of the biaxial vector
p(
r) about
n(
r), wherein
where
BH is the magnitude of the bend
Bn of
n(
z) in the helix,
BH =
Bn(
θH).
The magnitudes of these deformations, BH and TH, are uniform in space and, under the constraint noted above, satisfy the condition BH2 + TH2 = 1/Rmol2, describing the exchange of director bend for biaxial twist as θH is decreased, accounting in a fundamental way for the possibility of having a short pitch in the limit of small θH.
BH values calculated using Equation (8) are plotted in
Figure 3a vs. sin
θH using the
qH and
θH data of
Figure 2 and vs. temperature in
Figure S10. The
BH values fall quite closely onto a straight line passing through the (
BH, sin
θH) origin, indicating that changing
T or the 5CB concentration
x just moves points along the line, a quite surprising result that enables an immediate prediction: if
BH(sin
θH) is indeed linear in
sinθH, then we must have
qH(
θH)cos
θH =
S, the (constant) slope of “the line”. At high temperatures in the TB phase where
θH is approaching zero, we will have cos
θH ≈ 1 and therefore
qH =
S. Thus, for data on the line, the slope of the line should give the limiting helix pitch
pHlim near the N–TB transition as
2π/S. Fitting the
BH data in
Figure 3 to a line through the origin yields
S = 0.64 nm
−1 and, therefore,
pHlim = 9.8 nm. This value is plotted as the yellow dot in
Figure 3a and is indeed close to the measured pitches at high
T in the TB phase, deviating from the maximum pitches of the different mixtures by less than 10%, characteristic of the deviation of the
BH(
θH) data from the fitted line, which has
pH increasing weakly with increasing 5CB concentration.
The observation that
qHcos
θH is nearly constant leads immediately to the question of how to interpret this fitted value of
S. That
BH =
Ssin
θH means that
S =
B(90°), the maximum achievable value of bend of
n(
z), obtained when
θ is extrapolated to 90°. However, in the helix, we also have a director twist of magnitude
which, if
qHcos
θH =
S, is given by
meaning that the twist elastic energy density
UT =
KTTH2/2 grows strongly with increasing
θH, effectively setting an upper limit of
θHmax ≲ 35° on the physically achievable range of
θH and making the limit
θH → 90° non-physical for the TB helix.
2.3. Pure Bend Regime (Hypothetical)
We propose, alternatively, that the extrapolation of
B(
θ) to
θ = 90° represents a completely different physical situation, the one exhibiting the maximum preferred bend of
n(
z). This must require a geometry: (i) in which there is only director bend (pure bend (PB)), (ii) in which this bend has its preferred value everywhere (constant magnitude of director bend), and (iii) in which
θ = 90°, that is
n ⊥
z. These conditions are uniquely realized in the geometry of
Figure 3c, in the system of CB7CB molecules in which their atoms are attracted to a cylindrical surface of variable radius, packed, and equilibrated. At low temperatures, this condition maximizes the number density, a condition explored by packing DFT-based (DFT/B3LYP/6-31G**) space-filling models of rigid all-trans CB7CB. Maximum density in a pure bend geometry is achieved when the molecules are arranged with
B(
r) and
n(
r) parallel to the
x-y plane (
θ = 90°) and on the cylinder of preferred radius 1/
Bmol. Since the shape of an extended CB7CB molecule matches a circle reasonably well, the preferred PB radius can be estimated from the construction shown in
Figure 3b, which minimizes the mean square atomic distance from a circle by varying the circle radius
Rmol. That is to say, for CB7CB, we take the preferred bend
Bmol. in
Figure 3c to be the inverse of the molecular radius of curvature
Rmol, indicated in
Figure 3b. The resulting effective molecular radius of curvature of CB7CB is found to be
Rmol = 1.58 nm. This corresponds to a molecular bend of
B = 1/
Rmol = 0.63 nm
−1, which is remarkably close to the slope
S = 0.64 nm
−1, independently derived in
Figure 3a from the
θH,
pH data. These data are plotted in
Figure 3a, with
S as the red half dot and
Bmol as the yellow half dot at
θ = 90°, where
B(
θ) extrapolates to
S.
Remarkably, we find that the
BH vs. sin
θH trajectory of the data in
Figure 3 apparently could have been predicted from the molecular quantity,
Rmol, even at small
θH where there is little director bend left in the structure, and, furthermore, that the limiting pitch at the N–TB transition can be given in terms of the molecular quantities as
pHlim = 2
π/S = 2
π/Bmol = 2
πRmol = 9.8 nm. That is, if the TB system
BH vs. sin
θH trajectory has slope
Bmol, meaning that
then the helix pitch at
θH = 0° is just the circumference of the circle in
Figure 3b,c, describing the molecular radius of bend curvature,
Rmol. This circumference accommodates about four CB7CB molecules of all-trans length (~2.8 nm) with a slight overlap of the CN groups (
Figure 3c).
Thus, if the rationale for the twist-bend phenomenon is based on the effects of molecular bend, then this observation suggests that the
qH(
θH) data respond to changes of temperature and concentration by moving on the trajectory
BH ≈
Ssin
θH, which is, in turn, being controlled through
S by molecular bend by way of
Rmol, even at small
θH. This result is also surprising because, for TB phases and especially near the N–TB transition where
θH tends to be the smallest, there is little director bend left in the structure, the bend magnitude
BH(
θH) being quite small, as shown in
Figure 3. Nevertheless, if the
BH vs. sin
θH data are on the line, the bend is still proportional to
Bmol. Thus, the data of
Figure 3 indicates that the structural preference for the TB ordering to give a well-defined
qH is not a preference for constant bend. Several recent studies analyzing the elasticity of the TB helix pitch have found the director curvature bend energy to be orders of magnitude too weak to account for the observed TB pitch compressional elastic constant,
C, measured as
θH becomes small, and were led to propose local lamellar smectic positional correlations as an alternate source of rigidity [
21,
47].
Figure 3 shows, however, that since the data are nearly on the line, the helix pitch appears to be controlled by molecular bend, even in the absence of director bend at small
θH.
2.4. Polygon Chain Model
These considerations lead next to the question of the geometrical meaning of “the line” and, in particular, the relationship
Thus, on the line, for different values of
θH,
qH is such that
s(Δ
φ = 2π), the net distance traveled along the director contour path for a 2π increase in
φ (one complete turn of the helix, shown as the black dashed line in
Figure 1e), is independent of
θH, and furthermore given by
for all
θH, including the PB regime
θ = 90° (the path through Δ
φ = 2π around the circle of circumference 2π
Rmol). This invariance does not appear to be a symmetry of the system because the TB ground state at small
θH is entirely different from that of the PB at
θ = 90°. However, this is clearly the condition that connects the two regimes.
We can shed light on this condition by developing a geometrical model, the rectangle-triangle (RT) polygon chain, sketched in
Figure 4, that, by design, exhibits
s(Δ
φ = 2π) = 2π
Rmol over the entire range of
θH. This model is an assembly of rectangular and triangular plates connected into a periodic chain where the lines representing shared edges of rectangles and triangles are bendable hinges. The rectangles are attached to, and constrained to be locally parallel to, a flexible rod in the form of a helical spring representing the contour line of the director (green lines,
Figure 1e,f and
Figure 4b) on a cylinder of tunable radius
R =
Rmolsin
θH. The corresponding tunable pitch,
pH = 2π
Rmol/cos
θH, guarantees that such a chain of length
s(Δ
φ = 2π) = 2π
Rmol (in the example shown in
Figure 4b, this length corresponds to eight rectangles, 8
L) always makes exactly one turn of the helix. The helical rod is inextensible, enforcing by construction the condition that
s(Δ
φ = 2π) = 2π
Rmol. In the basic structural unit of the chain, consisting of two rectangles and a 45º isosceles triangle, a 45º bend in the director is enforced by the triangular hinge when the triangle and rectangles are all in the same plane, the condition at
θ = 90° where the whole construction lies parallel to the
x-y plane (
Figure 4b). This directly models the molecular organization of the PB limit in the cylindrical shell packing of
Figure 3c.
The hinge angle
βo = 45° = 360°/8 was chosen because, as discussed below, a diffuse feature in the non-resonant X-ray scattering in the TB phase indicates that there are ~8 molecular half-lengths in the TB pitch at small
θH. Since the
BH values fall on “the line” in
Figure 3a, there must correspondingly be in the PB regime ~8 segments around the 2π
Rmol circumference. Indeed, as shown in
Figure 5a, the PB regime is well modeled by the arrangement of four 45° bent rod molecules. We propose that each ring in this structure is stabilized by neighboring rings in an arrangement where adjacent rings have a difference in azimuthal orientation of 45°, such that the flexible molecular centers in one ring are over the regions of fluctuating end-to-end molecular contacts in the neighboring ring, an entropically favored association. This makes a construction like a cylindrical brickwork chimney, as discussed in the next section. With this choice, eight rectangle, long edges must make a complete turn, so the rectangle length,
L, is chosen such that 8
L ~ 2π
Rmol. The corresponding magnitude of the director bend is then
The helix can be tuned by pulling the ends of the rod so that they become separated along a line parallel to
z; the separation being the pitch,
pH, as indicated by the black arrow in
Figure 1e, which decreases
θH and makes the rod less bent everywhere along its length. The bend angle,
β, of the local elements, thus decreases from the maximum of
βo = 45° causing them to buckle, the triangle swinging out of the plane of the rectangles to cause less bend, and in the process, inducing a local relative twist
τ of the rectangle planes, which are free to rotate about the rod axis, as sketched in
Figure 4a. With this geometry, if the tilt of
n(
z) relative to
z is
θH, then the angle between the triangle and rectangle planes will be –-
θH, the condition that keeps the triangle planes always parallel to the
x-y plane. If the separation of the rod ends is increased and
pH approaches the rod length 2π
Rmol, then
θH → 0,
BH → 0 as
BH = sin
θH/Rmol, and the rod becomes nearly straight, with the local geometry changing as shown in
Figure 4a. The triangle plane is eventually oriented normal to the rod, and its initial induced bend in the rod of Δ
φ = 45° is now completely converted to an induced local relative twist about the director of the rectangle plane normals through Δ
φ = 45° at each hinge, as illustrated in
Figure 4b. The bend angle,
β, twist angle,
τ, and
βmol are geometrically related, as indicated in
Figure 4a.
The rectangles also represent the principal axes of the local biaxial nematic ordering tensor of the director field (director
n, flexoelectric polarization direction
p, auxiliary unit vector
m), as in
Figure 4a. Thus, as
θH increases, the overall structure of a single pitch is converted from the
θH ~ 90° state: a series of eight steps of 45º rotation of director bend and of local biaxiality about its
m axis on the circumference of the circle of radius
Rmol; to the
θH ~ 0° state: a series of eight steps of 45° twist rotation of the local biaxiality about its
n axis, on a path along
z of length 2π
Rmol. This scenario precisely maintains
qHcos
θH =
BH throughout the range of
qH, i.e., puts
BH vs. sin
θH on “the line” (
Figure 3a). We denote these ranges of large and small
θH, respectively, as the pure bend (PB) regime (
θH ~ 90°) and the twisted biaxiality (TBX) regime. The RT model directly shows that the structural stability of the local elements through the transition from pure bend to twisted biaxiality is what is required to maintain the compressional elasticity of the pitch under the condition that
BH → 0 and director curvature elasticity drops out. Actual twist-bend phases typically have
θH ≲ 30°, so they are much closer to the twisted biaxiality limit than the pure bend. Thus, in “TB”, the twist should be taken to mean the twist of biaxiality. The
θH ~ 0° regime represents the state of the helix dominated by twisted biaxiality but having no macroscopic optical tilt. Such a state is achievable, as shown below. In the CB7CB mixtures, the TB phase appears to come in with a small but finite
θH, consistent with the optical, X-ray, and DSC evidence for a first-order N–TB transition.
The RT model can be made for any angle
βmol. If
βmol is small, then
with
β and
τ becoming the orthogonal projections of a vector of magnitude βmol, constrained to move on a circle (
Figure 4a). In the limit that βmol → 0 with L/βmol constant and assuming that the hinge bends remain highly flexible, the RT chain becomes like a sheet of paper bent into the accordion fold of a fan, with high bending rigidity in the radial direction, and low bending rigidity in the circumferential direction. Upon pulling the bend out, such a sheet will exhibit little elastic resistance against conversion from continuous bend to continuous twist. In the continuum limit
a result that is also derived in
Figure 1f from the projective geometry of the helix. The balance of bend and twist is controlled by
Bmol even in the limit of zero bend (
TBX ≅
Bmol −
BH2/2
Bmol).
2.5. Steric Oligomerization of Bent Molecular Dimers
In the RT model, the constraint that qH cosθH = Bmol, independent of θH, is built into the model by the fixed length of its chain of polygons, a condition that would seem most applicable to a system of locally bent flexible oligomer or polymer chains. In the dimer TB phases considered here, there are no chemical links between molecules, so it is necessary to understand, in the context of independent bent molecules, how such a similar, polymer-like condition could come about in both the PB and TBX regimes, how the PB and TBX regimes are linked, and, therefore, how biaxial twist in the absence of bend comes to be controlled by Bmol.
We propose that molecular bend and steric packing constraints of the condensed TB phase combine to stabilize oligomeric chains of molecules and that the brickwork packing motif, introduced in
Figure 3c and detailed in
Figure 5a, is the common structural feature that stabilizes the chains and connects the PB and TBX regimes. The brickwork packing of a pair of adjacent chains can be visualized as a string of segments, each containing a pair of oppositely directed molecular halves linked by interfaces, each containing the center of a molecule in one chain and the tails of two in the other. This motif has also been found in other molecular dimer liquid crystal structures [
48]. This assembly is stabilized by the well-known tendency for rigid and flexible molecular subgroups to nano-phase segregate [
49], with the flexible molecular centers most readily accommodating the fluctuations in a relative position of neighboring molecular ends or tails. We refer to a double helix chain formed in this way as a duplex helical tiled chain (DHT chain, DHTC). The intra-duplex tiled linking is responsible for the apparent fixed contour length,
s(Δ
φ = 2π) = 2π/
Bmol along
n(
r), and manifested in the construction of
Figure 4.
This proposal is supported by the observation of a diffuse, non-resonant X-ray scattering feature in the N and TB phases of pure CB7CB, having a peak on the
qz-axis at
qm ≈ 5.05 nm
−1. A similar peak is found in the N and/or TB phases of a variety of other bent molecular dimers, with
qm in the range 4 nm
−1 <
qm < 5 nm
−1 [
14,
15,
19,
22,
23,
31,
50,
51], as discussed in
Supplementary Figures S11–S13. The typical appearance of this TB phase feature is shown in
Figure S11, which plots the non-resonant X-ray structure factor of the TB phase in CB7CB, calculated from the molecular dynamics simulation in the TBX regime reported previously [
29]. The white ellipses indicate the on-axis peaks, which can also be seen in
Supplementary Figure S13, which plots
z-axis intensity scans
I(
qz) of CB7CB [
15] and of a DTC5C7/DTSe mixture [
31]. These scans indicate a periodic electron density modulation and, therefore, molecular positional ordering along the helix with a fundamental periodicity of
dm ≅ 1.25 nm, consistent with the presence of short-ranged periodic positional correlation of similarly structured molecular segments along
z. This finding supports the brickwork association proposal since this value of
dm is close to half of the molecular length
M = 1.4 nm of extended CB7CB, which is what required for the segment length in a brickwork tiling. In fact, a comparison of
dm with the extended molecular length for the bent molecular dimer systems for which
dm data is available noted earlier in this paragraph, shows that the condition
dm/
M ~ 1/2 appears to be a general trend, as illustrated in
Supplementary Figure S12. In CB7CB, since
dm is comparable to
pH/8 = 1.22 nm, it is close to the brickwork segment length in the PB regime on the
Rmol circle (
Figure 5a), an observation that can be taken as evidence for there being similar segments at small
θH.
We made an initial evaluation of whether the peak at
qm (the chain segment scattering) could be understood on the basis of a model in which the self-assembly of a pair of molecular chains is described as a periodic chain of half molecule-long segments, each connected to an adjacent segment by a nearest-neighbor harmonic spring. The structure factor of this standard model for 1D positional ordering exhibits only short-range order at finite temperature [
52,
53], as described by the monotonic increase in the mean square of fluctuations in the dynamic separation of pairs of elements of the chain with increasing mean separation:
where
n counts the segments along the chain. The 1D structure factor
I1D(
qz) fits the data for CB7CB [
15] and the DTC5C7/DTSe mixtures [
31] quite well (
Supplementary Figure S13), giving, in both cases, a distance along the chain of ~15 segments for translational order to be lost, i.e., where √<
δun2> becomes equal to the segment interval
dm.
Another common feature of the chain segment scattering,
I(
qz,
q⊥), in the materials listed in
Supplementary Figure S12 is that the width of the diffuse peak in the direction normal to the helix axis
z,
δq⊥, is significantly larger than
δqz, its extent along the helix, as seen in
Supplementary Figure S11. In some cases, this appears to be mosaically broadening due to alignment defects, but in the TB phase in the DTC5C7/DTSe mixtures [
31], for example, the narrow angular width in
q of the resonant
qH peaks shows that the sample is well aligned, and therefore that the broadening of
I(
qz,
q⊥) in the
q⊥ direction is intrinsic. The corresponding correlation lengths,
ξ⊥ and
ξz, have the inverse relationship, implying that the correlation volumes giving the diffuse non-resonant scattering are extended along
z, i.e., a chain-like periodicity along
z rather than layer-like correlations.
2.6. Duplex Helical Tiled Chain (DHTC) Structure of the Twist-Bend Phase
The challenge then is to develop a model of the TB helix with small
θH in which it is made up of at least pairs of molecular chains in a brickwork tiling with subsections along
z of pairs of antiparallel half-molecules, in which, for CB7CB, the structural twist between segments is 45°. To this end, we considered the organization of single-stick and space-filling molecular models consistent with the above requirements. The PB regime is readily modeled by the packing of all-trans space-filling models of CB7CB and, as in
Figure 5, by either two- or three-segment bent stick models having 45° or 30° bend, respectively. In the PB limit, brickwork tiling of either stick model gives a bend of 45° per segment and four molecules per ring (
Figure 5a) so that the change in azimuthal orientation
φ is 45° per segment.
The required structures are shown in various representations in
Figure 5b,c and
Figure 6, and
Supplementary Figure S14. The basic structural associations are of three molecules like that of the green, cyan, and yellow groups inside the black elliptical rings in
Figure 5a,b, wherein terminal groups of the cyan and yellow tuck into the volume of hard-to-fill space vacated by the bending of the green and can associate with the flexible central aliphatic linkers. This scenario is repeated for the next segment along
z, among a group rotated through 45° relative to the initial one and having the cyan molecule in the center, and so on for all
z. The stick models in
Figure 5b,c and
Figure 6 show that this structure is double-helixed, made up of two identical right-handed helical chains of molecules, each transforming into the other by a translation of a single segment length followed by a 45° rotation in azimuthal angle. The paired assembly of two chains is stabilized in both the PB and TBX cases by a combination of a constraint, in the former to be on the cylinder or in the latter to be in a tube created by neighboring chains and by the pressure exerted by the neighboring molecules. In the pairing of the single-strand chains, the overlaps stabilize the structure and the interlocking bends promote the filling of space. In the bent stick representation of
Figure 5c, the half-molecular rods can be taken to represent the principal axes of the halves of the molecule. Taking the half-molecule polarizability to be uniaxial, the effective optical anisotropy of a segment of the double helix can be obtained using the construction in
Figure 5c. Here, the white square at each level is marked with a black dot that marks the midpoint between the intersections of the two chains with the square. The dark green line connects the midpoints from square to square. Thus, in a given segment, the green line construction will give the orientation of the local principal axis of the average dielectric tensor with the largest refractive index, which we take to be the local director. Thus, the green line represents the trajectory of the optical
n(
r), which is also a right-handed helix. This construction shows that, in a given segment of the DHT chain, the tilts of the half-molecular optic axes away from
z, in this case by ~22°, tend to cancel one another, leaving a much smaller effective optical heliconical cone angle, in this case
θH ~ 11°. The magenta labels in
Figure 5c indicate the handedness of the various helices, with the single chains and the director helix being right-handed (RH). Interestingly, the double helix is left-handed (LH).
For clarity, the molecules in
Figure 5 are positioned with more symmetry than they will actually have in the typical case. Generally, the planes of the bent molecules in the helices of
Figure 5b,c will be tilted away from
z through an angle,
ψ, as shown in
Figure 6 and in
Supplementary Figure S14 for both signs of tilt from
z. The untilted case could occur at some particular temperature, like the unwinding of the helix in a chiral nematic at a particular temperature.
Figure 6 presents the fully formed DHTC structure in the pure TBX limit for which the optical director tilt
θH = 0°. In this structure, the projection of the halves of a given tilted molecule onto the
x-
y plane has an opening angle between them of 45° (
Figure 6c), the same as the rotation Δ
φ = 45° per segment. This condition requires a tilt of the molecular plane from
z of 9.9°. In this case, the two molecular halves in a given segment have parallel projections onto the
x-
y plane (
Figure 6a,c). Since they also have equal and opposite tilts, there must be a principal axis along
z of their average biaxial contribution to the dielectric tensor (
Figure 6d). Starting from this structure, a heliconical director field of finite
θH can be generated by changing the molecular tilt (
Supplementary Figure S14) or by helical deformation of the DHT chain (
Supplementary Figure S14). Introducing a director bend into the DHT chain reduces biaxial twist, following the geometric projection scenario of
Figure 3f and of the RT model in
Figure 4. This comes about as illustrated in
Figure 6d,e, showing that, on the boundary between the two duplex chain segments containing the halves of the red molecule (denoted by a black circle), the projections of the halves of its cyan and yellow molecular neighbors make a 45° angle to one another. As indicated in
Figure 6d, this corresponds to twist
τ = 45° for
β = 0 at a yellow triangle in
Figure 4a. In the presence of director (heavy green line) bend,
β, the rotation of these neighbors relative to the red molecule is of opposite sign (+
β/2, −
β/2) and applied on the projections, as on the edges of the yellow triangle, causing the black disc plane, with application of bend, to rotate about
p, remaining, as in
Figure 6b, coplanar with the yellow triangle as it reorients (
Figure 4a). Elastic deformation of the DHT chain then satisfies the RT model constraints, which put
BH(
θH) on “the line” in
Figure 3a. With Equation (10), we have
and
For
BH small, then, the reduction in biaxial twist is controlled by
Bmol, with
In the limit of large bend, twist is eliminated, and the structure evolves toward the PB limit (
Figure 6e).
It is useful here to recall Ref. [
16], in which we reported the TB phase in the mesogen UD68, a rigid-core molecule in the form of a ~120° bent rod with unequal-length arms. This phase appears to be in the TBX regime, in which case the structure of
Figure 6 may be applicable, with one small change, to describe this TB phase, namely that the spacing between black discs would alternate from one pair to the next, between the short-arm length and the long-arm length.
2.7. Three-Dimensional Heliconical State
The bulk TB phase is a 3D space-filling packing of DHTCs. The overall orientational ordering with uniaxial positive birefringence means that the DHT chains are generally running parallel to one another, making the TB a hierarchical nematic self-assembly of anisotropic, self-assembled oligomeric chains. In the packing of cylindrical objects that are helically modulated, the helical contours on adjacent facing cylinders cross each other (like the stripes on a pair of parallel barber poles of the same handedness if put into contact). This geometry tends to suppress melding of the chains and to maintain the cylinders as distinct entities in the packing. Each DHT chain is then effectively confined to an on-average cylindrical hole in the fluid by its neighboring chains, which exert an effective pressure like that coming from osmotic pressure in a depletion interaction. This picture is supported by the experimental finding that the
BH(
θH) data of all the mixtures lie on the same line in
Figure 3a, indicating that they behave as if they all have the same
Bmol (at
x = 37.5 we might have expected a significant dilution effect leading to a smaller
Bmol). The constancy of
Bmol suggests that in the structures determining the pitch, the DHT chains, in the case of the TBX, are comprised dominantly of the bent dimers and that the 5CB is a filler in between. The 5CB dilution lowers the phase stability and reduces
θH, but this all occurs with
qHcos
θH =
Bmol, implying chains under the same constraint: d
φ/d
s =
Bmol.
Next, we consider the steric packing of the DHT chains that make up the bulk phase. The fact that the resonant X-ray scattering from the bulk TB exhibits diffraction spots from oriented domains that are 3D smectic-like, that is, having resolution-limited width in
δqz, indicates that the long-range ordered pseudo-layer scattering objects are arrays of lamellar sheets extended in the in-plane direction [
30,
31]. This means that in the bulk TBX packing, the phase
φ of the twist in a DHT chain must become coherent with that of its neighbors, a condition that has been observed in nematic phases made by packings of chiral particles internally structured as a steric repulsive helical object, realized, for example, in suspensions of helical flagella [
40] and in the extensive simulations of steric helices of Kolli et al. [
41]. The existence of the DHT chains opens up the possibility of a number of other arrangements of these objects such that, for instance, it is feasible to observe a transverse twisting of these DTH chains akin to a TGB* phase in some range of parameters. Another example relevant to the TB phase is the helical nanofilament phase found in neat bent-core systems [
9,
54] in which chiral filamentous bundles of a few smectic layers achieve macroscopic phase coherence of their twist solely by interacting through their periodic biaxiality.
The Kolli simulations appear to be particularly applicable to describe the interaction of, and the potential of, long-range phase ordering for, the DHT chains for finite
θH in the TBX regime.
Supplementary Figure S15 shows an example of the systems of interacting particles employed by Kolli et al. [
41], composed of rigid helical chains of contour length,
L, made of truncated hard spheres of diameter,
D. Comparison of the Kolli particles with the steric shape of the DHT chains of CB7CB, made in
Supplementary Figure S12 for the
θH = 10° case in
Figure 5b, shows particles with helical radius
r/D ≈ 0.2 and pitch,
p ≈
L ≈ 10
D match the CB7CB DHT chains quite well. The Kolli et al. phase diagram for
r/D = 0.2, also reproduced in
Supplementary Figure S15, shows that particles having
r/D = 0.2 systematically give I, N, TB, and smectic phases, with the TB range decreasing as the pitch becomes comparable to and longer than the particle contour length. Thus, the single-pitch duplex CB7CB chains should be able to order into a 3D TB phase if sufficiently long and rigid. The TB range in
Supplementary Figure S15 is limited with increasing volume fraction by the appearance of smectic ordering, corresponding to the positional ordering of the particles into smectic layers of thickness comparable to their length. In the case of living polymer chains, like what we propose for CB7CB, the effective particles will be transient and polydisperse, the latter condition well known to strongly suppress smectic ordering [
55], an effect which may expand the TB range.
For
p ≈
L and
r/D = 0.2, the Kolli et al. helical particles behave as if they are smooth, like those of Barry et al. [
40], which have helical glide symmetry. In these cases, if the steric helical interaction is reduced, for example, by reducing
r/D or making
pitch/L large, the system will revert to a simple nematic or smectic phase. However, the DHTCs are not smooth but are periodically structured with a local biaxial shape, so the role of variations of the steric shape along the DHT chain must also be considered.
Figure 4b,
Figure 5c,
Figure 6,
Figures S11, S15 and S14 all exhibit aspects of the biaxiality of the DHT chains.
Figure 4b for the
θH = 0° case and
Figure 5c show that the projection of the segments onto the drawing plane varies in effective shape along the chain, with a period equal to half that of the helix. This variation is also evident in the projection of the steric shape of a duplex chain in
Figure S15. Generally, each segment is biaxial, with a steric cross-section in the
x-
y plane that has the symmetry of an ellipse. This elliptical shape rotates in azimuthal angle φ along the chain (biaxial twist), as in a twisted ribbon of zero net local curvature, with a period of
pH/2, equal to four segment lengths, as is clear in
Figure 4b. In the TBX limit, this is the only periodicity of the DHTC. In a dense packing of the DHTCs, steric variations in shape, especially periodic ones, will lead to the development of correlations between chain positions along
z. This will be an especially strong effect if the oligomerization has substantially reduced the translational entropy for displacement along
z. In the helical nanofilament case [
10,
54], where the structural periodicity is the helix half-pitch, and the filament steric profile is almost circular, such that the neighboring filaments only weakly sense each other’s grooves, the filaments have a strong tendency to order with their biaxial twist in phase. In the DHTC case, sufficiently large ellipticity and packing density will lead to a 3D structure in which adjacent duplexes will align out of phase to facilitate packing. Twisted ribbons, for example, pack best when shifted by a quarter of their pitch.
2.8. Resonant Soft X-ray Scattering as a Probe of the Duplex Helical Tiled Chain Model
Given that we now have a fairly detailed structural model, we reconsider RSoXS as a probe of the heliconical structure of the TB phase. Interestingly, the first application of resonant scattering to LCs was to probe the heliconical molecular orientational ordering in chiral tilted smectic phases, in which the molecules are confined to layers, sorting out layer-by-layer sequences of azimuthal orientations of tilted molecules [
56,
57]. In this context, the general theory of resonant scattering was applied to the smectic case [
58]. This formalism has recently been applied in a comprehensive analysis of RSoXS scattering from the TB phase by Salamończyk et al. [
59,
60], which supports the duplex helical model.
RSoXS at the carbon K-edge (incident wavelength,
λ = 4.4 nm) gives a range of scattering vectors
q < 2π/2.2 nm, probing length scales through the nanometer range down to ~2 nm: molecular, but not atomic, size. In this
q range, molecular subcomponents such as the biphenyls in CB7CB act nearly as composite entities in the scattering process, being describable by second-rank molecular polarizability tensor scattering cross-sections, as in deGennes’ formulation of light scattering by fluctuations in director orientation [
36]. In analogy with visible light microscopy, RSoXS could even be used to visualize patterns of birefringence of LC phases and textures with X-ray resolution using depolarized transmission.
In probing the DHTC model of the TB phase, we first consider the RSoXS from individual filaments, illustrated in
Supplementary Figure S16, focusing on the essential qualitative features of the scattering in the simplest geometry. This figure shows about 1.5 pitches of the
ψ = 0 DHT chain in
Figure 5b, represented by space-filling models of CB7CB. In a typical experiment, the sample cell with the LC between silicon nitride windows is oriented parallel to the Figure plane, the TB helix axes of the LC are aligned parallel to the windows, and we consider illuminating a domain with the local helix axis vertical, as shown in the figure. Incident X-rays pass through the image plane and are forward-scattered onto a 2D detector behind. The incident and scattered directions can be chosen so that the scattering vector
q is parallel to the DHTC
z-axis. In this example, we take the incident X-ray polarization,
i, to be horizontal. The helical winding of the filament is apparent in the figure, with the director giving the orientation of a principal axis of the polarizability tensor following a helical trajectory, as shown in
Figure 1 and
Figure 4. According to deGennes, the depolarized scattering field amplitude probing director orientation is approximately
where
f is the outgoing polarization, nearly parallel to
z, and
δn(
z) is the angular deviation from
z. The key feature of this relation is that
E is linear in
δn(
z), so that the sinusoidal projection of the helix structure onto the
i-
z plane gives a sinusoidally varying scattering amplitude as
δn(
z) = sin(
qHz), which, in turn, produces scattering at
q =
qH. This is the basis of the claim that the depolarized scattering peak determines the helix pitch
pH = 2π/
qH, which clearly should be applicable to scattering from a single DHTC. In addition to the helical undulation, the DHTC exhibits smaller-scale roughness, a result of the precessing biaxiality discussed in the previous section. Inspection of the DHCT shows that, like in
Figure 4b, there are, in general, four distinct projections of the biaxial order on any vertical plane. Scattering from these variations has amplitude
where
p(
z) is the biaxial orientation vector. Here, the scattered amplitude is independent of the sign of
p(
z) so that the projection of
p(
z) onto a vertical plane
δn(
z) = sin(
qHz) generates a polarized scattering amplitude
Ep(
z) ∝ cos(2
qHz), the second harmonic of the scattering from the helix, explicitly showing that the periodicity of the biaxiality is a half-pitch: flipping the rectangles around Δ
φ = π does not change their biaxial polarizability. All the DHT chains presented here share this property. Salamończyk et al. have pointed out that in the scattering from columns of helical precessing tilted rods, averaging together pairs of columns shifted relatively in phase by a half-period of their biaxial polarizability renders the net polarizability the same in every quarter period and the second harmonic disappears. In the few experiments where the second harmonic might have been seen, it has not been observed [
30,
31,
59], indicating that such averaging may be taking place in the TB phase. In the case of the DHTCs, shifting a pair of chains by two segments and averaging will eliminate the second harmonic. However, achieving this in arrays of DHTCs may be problematic since frustration effects come into play on the closest-packed 2D hexagonal lattices.
2.10. Model Systems of Bent, Rigid Molecules
We sought to explore the role of molecular bend in other TB systems. The only others of which we are aware and for which data sets of
pH vs.
θH are available are the mean-field theoretical model for bent rods of Greco et al. [
61] and the Monte Carlo simulation of Greco et al. of hard spheres assembled to make steric circular arcs [
62]. These models are of particular interest because they treat collective TB behavior for bent objects that are rigid.
Figure 7a and b show plots of
BH vs. sin
θH calculated from
pH(
θH) for the arcs and bent rods, respectively. The black line in each plot is drawn through the origin and
BH(
θH) for the smallest
θH. The general behavior of
BH(sin
θH) is similar to that of the CB7CB mixtures in
Figure 3a, but with a tendency to increase relative to the black line with increasing sin
θH, which is also seen weakly in neat CB7CB (
Figure 3a).
In the case of the hard, circular arcs, we carried out the
Rmol construction of
Figure 3b, with the result shown in
Figure 7a, finding
Rmol = 12.6 in units of the sphere diameter, σ. The corresponding
Bmol = 1/
Rmol = 0.08/σ is comparable to the
B = 0.1/σ extrapolation of the black line, indicating a relation between the PB and TBX limits similar to that in the CB7CB mixtures. This makes the hard arcs a very interesting system for exploration of the DHTC structure.
Turning to the bent rod case,
Figure 7b shows that
BH(
θH) obeys
BH =
Ssin
qH rather well for
θH < 15°, with a slope
S = 0.56/
L, where
L is the length of one of the arms in the bent rods. In
Figure 7b, we have used the small-angle value of
S to extrapolate to
θH = 90° in order to determine the radius of the cylinder
R = 1/
S in the PB regime. The resulting construction using the shape of the simulated bent rods shows a quite reasonable PB limit. This is an exciting result because this model approaches understanding the TBX regime from a mean-field statistical mechanical approach that is entirely different from the geometrical model building that we have employed here. That it captures the essence of the geometry of “the line” offers an opportunity to understand in detail the evolution of the local geometry to keep the system on “the line” in the absence of changes in molecular conformation.
Both models consider rigid molecules and yet seem to exhibit the same essential geometrical behavior as the CB7CB system, which was rationalized based on nanophase segregation of flexible (the central alkyl linker and the tail ends) and rigid molecular subcomponents (cores): the molecular ends find entropic freedom by associating with the flexible cores. We propose that in systems of rigid, bent, hard particles, the analogous association is between the particle ends and the free volume available in the pockets of difficult-to-fill space created by molecular bend.