Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# A case of antiferrochirality in a liquid crystal phase of counter-rotating staircases

## Abstract

Helical structures continue to inspire, prompted by examples such as DNA double-helix and alpha-helix in proteins. Most synthetic polymers also crystallize as helices, which relieves steric clashes by twisting, while keeping the molecules straight for their ordered packing. In columnar liquid crystals, which often display useful optoelectronic properties, overall helical chirality can be induced by inclusion of chiral chemical groups or dopants; these bias molecular twist to either left or right, analogous to a magnetic field aligning the spins in a paramagnet. In this work, however, we show that liquid-crystalline columns with long-range helical order can form by spontaneous self-assembly of straight- or bent-rod molecules without inclusion of any chiral moiety. A complex lattice with Fddd symmetry and 8 columns per unit cell (4 right-, 4 left-handed) characterizes this “antiferrochiral” structure. In selected compounds it allows close packing of their fluorescent groups reducing their bandgap and giving them promising light-emitting properties.

## Introduction

A large part of condensed matter physics has been devoted to the study of the spontaneous long-range orientational ordering of electronic spins, leading to ferromagnetism and antiferromagnetism. Another intriguing self-ordering phenomenon is spontaneous synchronization of chirality, i.e., breaking the mirror symmetry and favoring either right- or left-handedness. A much-studied problem in biology is understanding how it happened that most amino acids in nature are strictly left-handed, causing the alpha-helix in the proteins they make to be right-handed. It is not accidental that such questions have occupied the minds of biologists and solid-state physicists alike1.

In liquid crystals (LC), long-range order can generally be achieved more easily and quickly than in crystals, making them interesting for functional materials where order is required for their function. Such long-range order is normally in molecular orientation and/or position of phase-separated domains (e.g., aromatic and aliphatic). However, in exceptional cases order of the third type, chirality, is also self-propagating in LCs, without the molecules themselves being chiral. This can happen e.g., in layered phases of bent (banana-shaped) molecules. While a stack of upright bananas on a shelf has mirror planes, a stack of tilted bananas has none, hence is chiral, and the domino tilting effect can ensure long-range chiral order2.

It has also been discovered recently that two common types of bicontinuous 3D network LC phases are always chiral3. These are the triple-network cubic4,5 and the tetragonal “Smectic-Q” phase6. They form in “polycatenar” rod-like molecules bearing more than one flexible end chain. Their phases consist of column-like segments with the molecules arranged in rafts of two or three lying normal to the column axis. Such segments are joined at three- or four-way junctions. To avoid clashing of the end chains, the successive dumbbell-shaped rafts twist resulting in the segments being helical. The question is why the same twist sense is maintained throughout the entire infinite network. It was proposed3,5 that this long-range chiral order is enforced at junctions, where all three or four merging segments must be homochiral, i.e., have the same twist sense, to achieve close packing.

While the bicontinuous network phases are relatively rare7, columnar LC phases, with infinite parallel columns and no junctions, are common. Most often they have 2D-hexagonal symmetry (Colhex) and form in compounds with disc- and wedge-shaped molecules8,9,10,11,12,13,14, in honeycomb-forming rod-like amphiphiles with side-chains15,16,17, and in the polycatenars18,19,20. Columnar phases of many compounds are 1D semiconductors and have desirable electro-optic properties when aligned. They can be used e.g., in light harvesting and light emission, in sensors, ionic conductors, etc.21. Combining these features with chirality could enhance their versatility as functional materials, e.g., as emitters of circularly polarized light or as membrane materials for separation of left- and right-handed isomers (enantiomers). As columnar mesogens normally have multiple pendant chains, it is likely that the adjacent molecular rafts, or strata, are mutually rotated. However, it is an open question whether uniform twist sense can propagate over long range in junction-free columnar phases as it does in bicontinuous networks. What would eliminate helix reversals and ensure that the same twist sense is maintained?

There have been a number of claims of helical columnar LCs of achiral compounds22, some of them ingeniously designed with high barrier to conformational flip (e.g., “propeller-blade” bipyridines23. The claims are often based on observed lines perpendicular to the fiber in fiber X-ray patterns (layer lines), but their long-range isochirality has not been proven convincingly. Layer lines, although a feature of diffraction on helices, are just as likely in fiber patterns of completely randomly twisted chains or columns24. In the extreme, placing a brick rotated on top of another creates a portion of a helix, but no long-range chirality can be achieved as a random reversal of rotation direction cannot be prevented. Even when the molecules have a chiral group or a chiral dopant is added, while there is amplified preponderance of one helical hand, the length of homochiral sequences is limited and dopant-dependent. This short-range order is characteristic of a paramagnet, with the dopant acting as the external field. However, the question is whether the twisted sense can propagate to true long range on its own in columns of true LCs, including those of achiral compounds. Can column helicity be of ferromagnetic type? Can we have examples of ferrochirality25 or even antiferrochirality?

Helix reversals cannot be prevented in an isolated 1D column. In contrast, in column-containing crystals helical sense can propagate as reversals are inhibited by interactions with “correctly” positioned molecules on surrounding lattice; prime examples are various soft-crystal forms of dendronized perylene bisimides (PBIs) bearing wedge-shaped end groups26, different “ordered” discotic phases27 or in numerous achiral crystalline synthetic polymers28. Many helical soft crystals transform to Colhex at higher temperatures, e.g., PBIs and some polymers such as Teflon. Teflon’s chains lose their regular helicity on the transition from crystal to columnar LC at 30 °C29,30, probably with soliton-like helix reversals31.

Here we report a 3D LC phase where we show that homochirality of helical columns can indeed propagate to long-range even without molecular chirality or network junctions. The phase consists of counter-rotating twisted columns with elliptical or star-like cross-sections. There is long-range order between helical columns but not between molecules, hence it is a liquid crystal and not a crystal. Two of the compounds involved contain a chromophore; the helical configuration enables their efficient π–π stacking, making them interesting for optoelectronic applications. A simple theory based on interacting quadrupoles is developed confirming that, compared to the alternatives, the observed structure is energetically the optimal packing of helices of linear polycatenar dumbbells.

## Results and discussion

### The materials

The Fddd phase is found in two types of polycatenars with three chains at each end, one with a straight (FCN16 and FO16) and the other with a bent-core (IC3/n—see Fig. 1). The high-T trigonal columnar phase of IC3/n compounds was described in ref. 32. For synthesis, see “Methods” and Supplementary Methods 1.

All compounds form a uniaxial columnar phase at higher temperatures, IC3/n with trigonal (Coltri)32, and FCN16 and FO16 with hexagonal symmetry (Colhex), see Table 1, Fig. 1, Figures further below and Supplementary Figs. 1, 7, 8, 10, 11, and 16. Both phases are also referred to collectively as Col. Below Col all except IC3/10 also form a 3D LC phase, which is the focus of this study. Its structure is determined by a detailed X-ray study and by atomic force microscopy (AFM), see the next section. Fan (“spherulitic”) textures, typical of columnar developable domains33, appear in polarized optical micrographs (POM, Fig. 2 and Supplementary Figs. 24). Comparing insets in Fig. 2a, b, we see that the slow axis, hence the aromatic cores, is roughly perpendicular to the column axis. The increase in birefringence following the phase transition is most pronounced in FCN16 (Fig. 2c, d) where, from Michel-Levy color chart, optical retardation increases by 1/3 (from ~150 to ~200 nm), indicating molecular alignment closer to column normal in the low-T phase.

### Structure of the mesophase

The absence of second harmonic (SHG) in the high-T phase of FCN16/FO16 means that the hexagonal Laue symmetry of the small-angle X-ray (SAXS) pattern (Fig. 3a and Supplementary Figs. 8 and 11) comes from a Colhex phase, unlike in IC3/n where the presence of SHG indicates trigonal symmetry32. As in the Col phase, only a diffuse wide-angle X-ray scattering (WAXS) maximum at 4.5 Å (q = 1.40 Å−1) is observed for IC3/n in the low-T phase (Supplementary Fig. 5b). This confirms the LC nature of the phase. An additional narrower but still diffuse peak appears at 3.4 Å (1.85 Å−1) in the low-T phase of FCN16 and FO16 (Fig. 3d and Supplementary Figs. 7 and 10). Its near-meridional position as a vertical streak close to the horizon in the grazing incidence (GIWAXS) pattern of a sheared film (Fig. 4b–e) identifies it as arising from good π–π stacking along z (column) axis of parallel aromatic planes.

The transmission powder SAXS patterns of low-T phase of the compounds are rather complex (Fig. 3a and Supplementary Figs. 6a and 12a). The observed diffraction peaks can be indexed to an orthorhombic lattice (a~180 Å, b~110 Å and c~40 Å, Table 1) with the help of grazing incidence SAXS (GISAXS) on surface-oriented thin films (IC3/12, Fig. 3b; IC3/14, Supplementary Fig. 6b; FCN16, Fig. 3c, Supplementary Fig. 9; FO16, Supplementary Fig. 12b). The GISAXS patterns confirm the reflection conditions h + k, h + l, k + l even, and h + l, h + k equal to 4n for h0l and hk0. This narrows the choice of spacegroup to F2dd (No. 43) or Fddd (No. 70, see Supplementary Table 2). As no 0kl reflection was observed, the deciding condition k + l = 4n could not be tested. However, since Fddd has an inversion centre and F2dd has not, we performed a second harmonic generation (SHG) test. As can be seen in Supplementary Fig. 15, while the noncentrosymmetric Coltri phase of IC3/n generates the second harmonic, isotropic liquid and the orthorhombic mesophase do not. This restricts the choice of spacegroup to Fddd alone.

Electron density (ED) maps were reconstructed from SAXS intensities (Supplementary Section 5 and Supplementary Tables 1 and 5). According to the map of the low-T phase of IC3/12 (Fig. 5a, c and Supplementary Movie 1) there are eight columns in the unit cell. As in the Col phase32, the calculated number, μ, of molecules in each column stratum of thickness 4.6 ± 0.1 Å is 3 (Table 1, Supplementary Table 8, and Supplementary Discussion 2). The cross-section of an averaged column is rounded triangular, similar to that in the Col phase, three bent-core molecules back-to-back forming a three-arm star32. However, in contrast to the high-T phase, the orientation of the stars changes monotonically with increasing z-elevation, resulting in four left- and four right-handed helical “staircases”. Assuming uniform twist, each consecutive three-arm stratum, or step, is rotated by 13.2° around the column axis—see the model in Figs. 5b, d and 6a. A molecular model of a layer with periodic boundary conditions, using the experimental lattice constants and subjected to 30 cycles of molecular dynamics (MD) annealing is shown in Fig. 6c; note efficient space-filling.

The ED map of the Fddd phase of FCN16 (Fig. 5e, g and Supplementary Movie 4) shows even more clearly the twist of the columns, here appearing like twisted ribbons. The column section is now oval, with close to μ = 2 parallel linear molecules per column stratum. This is referred to below as “dimer” (Table 1 and Supplementary Table 8); see also Fig. 5f and the MD snapshot in Fig. 6d. Even though the strata are thinner than in IC3/n (3.4–3.5 Å), the twist angle between strata in the linear compounds is larger (18–19°) than that in IC3/n (13°) (see Supplementary Table 10), making the pitch (=2c) considerably shorter—compare models in Fig. 5d, h.

The discontinuous drop in μ in a thermal transition, as observed here from 3 to 2 in FO16 and FCN16, is rare but not unique in columnar phases. A discontinuous drop from μ = 4 to 3.5 has been observed in pizza-slice-shaped molecules of trialkoxybenzoate salts (“dendron ejection transition”)34, a process that can even result in unusual columnar superlattices35.

All columns in (100) rows along b axis co-rotate, either clockwise or anticlockwise, but along {110} diagonals they counter-rotate, alternating between left and right helicity (Fig. 5a–h). xy sections of the model at different z-elevations (Fig. 6a, b and Supplementary Figs. 17 and 18) show that the linear dimers of FO16 and FCN16 succeed in evading clashes at every level (see Supplementary Movie 6), while the co-rotating stars of IC3/n unavoidably clash along b at z = 4.5 Å and z = 13.5 Å (see Supplementary Movie 3). Interestingly, while the a/b ratio for FCN16, FO16, and IC3/14 is exactly √3 (186.9/107.9, 186.0/107.4, 186.8/107.8), meaning that the column axes sit precisely on a hexagonal lattice, for IC3/12 the ratio is 1.635, i.e., the unit cell is stretched along the clash direction b (Fig. 6a). The deviation from hexagonal grid also causes separation of (400) and (220) reflections (Fig. 3a). Probably the “padding” by the C12-chains is insufficient to fully mitigate the clashes. Significantly, IC3/10 does not form Fddd at any temperature (Table 1).

The alternating helices in Fddd phase are also observed directly by AFM on FCN16 (Fig. 5j and Supplementary Figs. 20 and 21). Columns without helicity are observed in the Col phase in Fig. 5k, matching the ED map and model in Fig. 5l, m. The ~53 Å distance between helical columns in Fig. 5j and the ~35 Å helical pitch match almost exactly the ED map and the proposed structural model of Fddd. Edge dislocations, stacking faults of mismatched chiral columns, and helix reversal defects are also seen in Fig. 5j.

An interesting consequence of the different column cross-sections in IC3/n and FCN16/FO16 appears to be a higher bending modulus of the star-like IC3/n columns compared with that of the bendable ribbons of FCN16/FO16. The difference is evident in the type of disclinations in POM images in Fig. 2; only FCN16/FO16 columns can tolerate the high curvature at the center of a s = +1 disclination (Fig. 2a, inset 1). In contrast, the stiffer IC3/n columns prefer ends to sharp bends. Thus the central column circle splits into two semicircles with a larger radius (two s = +1/2 disclinations), with a bundle of straight but finite length columns filling the central part (Fig. 2a, inset 2).

### 3D order of helices, not molecules

The high surface alignment in GISAXS/GIWAXS experiments allows some further structure refinement and clarification of the true nature of the mesophase. Significantly, the 3.4 Å peak in powder WAXS (Fig. 3d) corresponds to the streak at qz = 2π/3.5 Å in GIWAXS of columns aligned parallel to substrate (Fig. 4b, Supplementary Fig. 13b). The 0.1 Å difference comes from the 8° tilt of the aromatic planes from the xy plane, causing a displacement of the intensity maximum away from the z-axis (horizon); see I(qxy) intensity profile, Fig. 4c. The tilt is taken into account in the model in Fig. 5i, which is more realistic than the simplified version in 5 h.

Notably, while there are many small-angle Bragg reflections with resolution-limited width, i.e., with correlation lengths of at least 1800 Å (domain size >5000 Å), there are no Bragg reflections at wide angles. This means that Fddd is a true LC and not a crystal. A more precise distinction is provided by the 3.5 Å streak. From the qz-profile the correlation length along the columns is $${{\xi }}_{z}=2/\triangle {q}_{z}$$ = 67 Å (Fig. 4d). However, from the qxy-profile along the streak, we get $${{\xi }}_{{xy}}=2/\triangle {q}_{{xy}}$$=5.6 Å, or “crystal size” of πξxy = 18 Å, which is smaller even than the width of one single column (>50 Å). These measurements confirm that the Fddd is a true LC and not a crystal, soft or otherwise, as in many reported helical column-like soft structures. Thus, while there is true 3D long-range order in packing of helical columns, there is no intercolumnar correlation between positions of individual molecules. One may think of the Fddd phase as consisting of orderly interlocked helical tubes filled with mobile molecules with no preferred position and no lattice. A somewhat similar situation has been observed previously in a chiral helicene, but with substantially lower order in “tube” packing36.

According to our model, there is a fundamental difference between helical columns assembled from rod-like molecules on the one hand and those from disc- and plate-like molecules on the other. For rod-like molecules such as those studied here, which include both straight and bent rods, the twist angle ϕ between successive rafts, or strata, in a column is a balance between the steric repulsion R(ϕ) between the alkyl brushes at rod ends (red dashed curve in Fig. 7a) and the attraction A(ϕ) between the π-conjugated rods (blue dashed curve). R(ϕ) has a maximum when the rods are parallel (maximum clash, ϕ = 0) and a minimum when the rods are perpendicular (ϕ = 90°). A(ϕ) has a sharp minimum when the rods are parallel (ϕ = 0°) as there is maximum π-π overlap along the entire length of the rods. The resulting sum energy function S(ϕ) (black solid line) therefore has two minima at ±ϕ1 where ϕ1 is a small angle. In the Fddd phase the interaction with neighboring columns will break the degeneracy, making the two minima unequal, thus favoring one twist sense over the other—see the separate section on column packing below. Due to the molecular motion in the LC, the averaged columns are smooth twisted ribbons.

In contrast, the situation is different where discs or plate-like mesogens two or more benzene rings wide are involved. Here A(ϕ) has minima additional to that at ϕ = 0, because by rotating the successive platelets around the column axis further aromatic rings will come in contact producing local minima (Fig. 7b). The twist angle ϕ2 will then be fixed by the global minimum in the wavy sum function S(ϕ) = R(ϕ) + A(ϕ). ϕ2 will typically be considerably larger than ϕ1 in narrow rod molecules, as seen in disc-like37 or elongated plate-like mesogens such as PBI26,38, or in hat-like molecules like cyclotriveratrylene39. Instead of being a continuous helix, such knobbly columns are discontinuous helices that tend to interlock with neighboring columns resulting in crystals with 3D long-range positional order of molecules.

### Antiferrochirality

Spontaneous long-range alignment modes of electric dipoles in SmC* LCs (layers of tilted chiral molecules) and in some bent-rod LCs, are well known, respectively, as ferroelectric or antiferroelectric40. In the same vein, it is useful to regard the present Fddd phase as “antiferrochiral”. This helps distinguishing it from the reported cases of helical columns, where the range of homochirality is short and dependent on the possible presence of chiral substituents or dopants, which act as external field acts on a paramagnet. The well-known frustration of antiferromagnets on a triangular lattice is also present here, causing symmetry-breaking of hexagonal to orthorhombic.

### Comparison with cubic LCs and Fddd phase in block copolymers

A phase with Fddd symmetry, albeit with a distinctly different structure, has been reported previously in block copolymers (BCPs)41,42. The Fddd lattice parameters (a is defined as the largest lattice parameter as used in this paper, and in the following parameters from literature have been changed accordingly for easier comparison.) in poly(styrene-b-isoprene) or poly(isoprene-b-styrene-b-ethyleneoxide) are naturally an order of magnitude larger, a varying from 750 to 2125 Å43. However, the a:b ratio is similar to that in our LCs, about $$\sqrt{3}\!:$$1 (Supplementary Table 11). Unlike our LC phase which consists of columns, the Fddd in BCPs is a single network bicontinuous phase with 3-way junctions, as shown in orange, scaled, in Fig. 6e, f and Supplementary Fig. 19, together with our IC3/12 model. In Fddd phase of PS-b-PI, PS forms the network embedded in PI34. In a crude analogy, the aromatic cores in our model can be said to occupy PI-rich areas in the BCP phase. Evidently, the IC3/n helical columns fill the space between network segments, with the twist sense of the columns matching that of the surrounding network. The Fddd phase in BCPs is found between lamellar and Cubbi phases ($${Ia}\bar{3}d$$ or I4132)42, while our Fddd appears below the columnar.

The chiral columns in Fddd phase remind us of helical segments of networks forming the bicontinuous phases mentioned in the introduction. These are the double-gyroid cubic $${Ia}\bar{3}d$$ phase7, the triple-network cubic I234 or the tetragonal “Smectic-Q” (I4122)6. While the latter two are chiral even in achiral compounds, the two networks in the $${Ia}\bar{3}d$$ have opposite hands, resulting in no overall chirality. The gyroid can thus be considered antiferrochiral and the triple-network and the SmQ ferrochiral. In all three phases the twist between successive molecular strata is 8–10°. As in the Fddd, this twist balances the attraction between the cores A(ϕ) and the repulsion R(ϕ) of the end chains. Since the compounds forming the above bicontinuous phases contain only 3–4 end chains, their smaller twist than the ~18° in the six-chain FCN16/FO16 is understandable.

In bicontinuous phases, the homochirality of each network is enforced at junctions, where all three or four merging helical segments must be homochiral to minimize steric clash. However in the Fddd there are no junctions and the columns are parallel. What propagates their homochirality is the 3D order of interlocking helices which ensures that helix reversals, such as that in Fig. 5j2, are eventually corrected by the surrounding lattice. In columnar phases of covalent or self-assembled discs columns have cylindrical symmetry, hence such 3D interlock is unlikely. The present results thus suggest that a noncircular column cross-section, such as oval or star-like, is beneficial in facilitating 3D interlock that could sustain homochirality.

### A simple theory of packing of helices: Fddd and alternative models

The interaction between two FCN16/FO16 dumbbell dimer rafts from neighboring columns at the same elevation can be described quantitatively as that between two linear quadrupoles oriented at β1 and β2 (Fig. 6g, h). The repulsion of equal charges at the ends of the rod-like quadrupole represents the steric repulsion of the pendant flexible chains. Considering each helix as a series of such twisting quadrupoles, the interaction energy E between two neighboring columns is the sum of interactions between their corresponding quadrupoles at all elevations. Our calculations (Supplementary Discussion 1) show that E is determined by constant α, defined as $$\alpha ={\beta }_{2}-{\beta }_{1}$$ for co-rotating and $$\alpha ={\beta }_{1}+{\beta }_{2}$$ for counter-rotating columns (Fig. 6i, j).

For a row of columns, the minimum E, Emin, is found for counter-rotating columns with $$\alpha ={\beta }_{1}+{\beta }_{2}=\pi /2$$ (Fig. 6j). If all columns must co-rotate, a higher Emin has $$\alpha ={\beta }_{2}-{\beta }_{1}=\pi /2$$ (Fig. 6i). In $$\frac{3{\phi }^{2}}{32{r}_{12}^{5}}$$ units, where ϕ is the quadrupole moment and r12 the distance between quadrupoles, Emin per dimer/quadrupole is −29 in a counter-rotating, 3 in a co-rotating, and 6 in a random row. However, the condition to keep all neighboring columns counter-rotating cannot be satisfied on a 2D-hexagonal lattice. The obvious choice then is to keep as many counter-rotating neighbors as possible, as shown in Fig. 6k1. A 2 × 2 supercell is assumed, as both counter- and co-rotating minimum-energy rows have a two-column repeat (Fig. 6i, j). We have two left- and two right-handed helices in the 2 × 2 supercell, and each has four counter-rotating and two co-rotating neighbors. Energy minimization shows that it is impossible to have $$\alpha =\pi /2$$ for all rows simultaneously; instead the best solution has $$\alpha =\pi /2$$ for co-rotating and $$\alpha =\frac{5\pi }{12}$$ for counter-rotating columns (Fig. 6k1), with interaction energy per dimer about −45.6 (Supplementary Table 13). Co-rotating columns having $$\alpha =\pi /2$$ means that neighboring helices are vertically shifted by quarter-pitch (Fig. 6i). This results in a 3D orthorhombic cell with Fddd symmetry, fitting almost exactly our observed dimer orientations at different elevations (cf. Fig. 6k1 and b/3.5 Å).

Other possible arrangements of helical columns have also been explored (Fig. 6k2–k4). The energies for models 1–4 were −45.6, −36.75, 15, and 15, respectively. The energy of Model 1 with two LH and two RH rotations, our proposed Fddd structure, was clearly the lowest. In comparison, the calculated energy for a completely uncorrelated Col phase is 18. For more details, see Supplementary Discussion 1. It is therefore concluded that the antiferrochiral Fddd structure is a result of packing optimization of helical columns with a noncircular cross-section.

It is also interesting to note that the homochiral Model 3, with four ribbons in a unit cell, all twisted in the same direction and which represents a ferrochiral structure, has high energy and is therefore unlikely to form spontaneously in non-chiral compounds. At the same time, however, the superlattice Model 2 with 3 LH and 1 RH column (or vice versa), having an enantiomeric column ratio 3:1, and thus being chiral, has an energy not too much higher than the Fddd. There is thus a possibility that it might be observed in some systems in which other structural factors would override the energy disadvantage calculated by our simple model. Being equivalent to a ferrimagnet, its structure could be regarded as ferrichiral. In fact, it mirrors precisely the prototype ferrimagnet Fe3O4, which has three Fe moments parallel and one antiparallel.

### Spectra and π–π stacking

With conjugated fluorophore cores, FCN16 and FO16 are potential electroptic materials44,45. Their UV–vis and fluorescence spectra are shown in Fig. 8 together with HOMO and LUMO orbitals calculated by density functional theory (DFT). The bandgaps measured from FO16 and FCN16 spectra in Col phase are 2.15 and 1.52 eV, smaller than calculated (Fig. 8 and Supplementary Table 12) due, at least partly, to the neglect of intermolecular conjugation in the calculation. The bandgaps narrowed by a further 0.07 and 0.17 eV on Col-Fddd transition, causing red shifts of 20 and 50–60 nm in FO16 and FCN16. These results indicate an increase in intermolecular conjugation upon the transition. In the Col phase 4.5 Å thick three-molecule strata orient randomly in the xy plane (Figs. 4a and 5l, m and Supplementary Fig. 13a). On Col-Fddd transition the number of molecules drops to about two and the inter-strata spacing to 3.5 Å.

As mentioned in Supplementary Discussion 2, the difference between the 4.5 Å and 3.5 Å spacing in columnar phases is usually associated with the change from the T-type stacking mode (aromatic planes on successive molecules close to perpendicular) to the parallel mode (aromatic ring planes nearly parallel). In FO16 and FCN16, this change upon the Col-Fddd transition is also consistent with the change to a more strongly negative birefringence (rings closer to perpendicular to the column; Fig. 2c, d) and a narrowing of the bandgap. These changes are absent in IC3/n compounds which maintain the T-type stacking in both phases.

The particularly large bandgap narrowing in FCN16 is consistent with the potent malononitrile acceptor (Hammett constant for CN σp = 0.66 (http://www.wiredchemist.com/data/hammett-sigma-constants) interacting strongly with the electron-rich adjacent aromatic cores46 and stabilizing the π-stack47. FCN16 stands out also by its large transition enthalpy (17J/g) and increase in birefringence (Fig. 2c, d). The near-IR fluorescence, expected in FCN16, could be useful in achieving greater light penetration through opaque tissue48.

In summary, a complex 3D liquid-crystal phase is discovered in bent- and straight-core compounds, having orthorhombic Fddd symmetry and consisting of counter-rotating helical columns. Unlike previous reports of helical columnar LCs, in spite of being non-crystalline, helical periodicity is long-range and spontaneous. As antichiral near-neighbor interactions are favored, the structure is equivalent to an antiferromagnet, hence termed antiferrochiral. Its discovery confirms the universality of the principle of twisting polycatenar columns as building blocks for complex 3D self-assembly, previously established in bicontinuous network LCs, albeit here without network junctions. The long-range twist sense is shown to propagate through mere steric inter-helical interaction. Comparison of different helical packing models based on linear quadrupoles confirms that in achiral compounds the antiferrochiral Fddd is the lowest energy solution, while the ferrochiral alternative is energetically the least favored.

Other helical column phases could be designed on similar principles. As the chiral Model 2 structure, with a 3:1 ratio of column enantiomers, has an energy not too high above that of the Fddd, by modifying the molecule it may be possible to obtain it by spontaneous assembly. With suitable fluorophores this could lead to new organic circularly polarized LEDs, enantioselective filters or sol-gel templates for chiral porous ceramics.

## Methods

### Synthesis

IC3/n compounds were synthesized by a Cu(I)-catalyzed click reaction32. Synthesis and chemical characterization of the straight-core compounds FO16 and FCN16 is described in Supplementary Methods 1.

Reactions requiring an inert gas atmosphere were conducted under nitrogen and the glassware was oven-dried (105 °C). Tetrahydrofuran (THF) was distilled from sodium prior to use. Commercially available chemicals were used as received. 1H-NMR and 13C-NMR spectra were recorded on a Bruker-DRX-300 spectrometer and a Bruker-DRX-400 spectrometer. The mass of the compounds was characterized by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) on a Bruker autofleX max instrument. Elemental analysis was performed using an Elementar VARIO EL elemental analyzer. Thin-layer chromatography was performed on aluminum plates precoated with 5735 silica gel 60 PF254 (Merck). Column chromatography was carried out on Merck silica gel 60 (230–400 mesh).

Transmission powder SAXS/WAXS experiments were carried out at station I22 of Diamond Light Source, U.K. Powder samples in 1-mm glass capillaries were held in a modified Linkam hot stage. Pilatus 2 M detector (Dectris) was used and the X-ray energy was 12.4 keV. GISAXS/GIWAXS experiments were done at BM28 of European Synchrotron Radiation Facility, France, and I16 of Diamond Light Source. 2D diffraction patterns were collected using a MAR165 CCD camera at BM28 and Pilatus 2 M at I16. Thin-film samples were prepared from melt on the silicon substrate. n-tetracontane was used to calibrate the sample to detector distance.

Electron density maps were calculated by inverse Fourier transformation using the standard procedure as described in International Tables for Crystallography. Integral intensities of all peaks were measured using Gaussian peak fitting. GISAXS intensities were used to help resolve overlapping peaks. More details are in Supplementary Information.

Molecular models were built using Materials Studio (Accelrys). Geometry optimization and molecular dynamic annealing were performed using Forcite Plus module with Universal Force Field. NVT annealing dynamics was performed through 30 cycles between 300 and 600 K, with a total annealing time of 30 ps. Frontier orbital distributions of fluorene derivatives were calculated using the DFT at level B3LYP, 6–31 G(d) basis set.

AFM imaging was done in tapping mode on a Bruker Multimode 8 instrument with Nanoscope V controller. The sample was dissolved in toluene and spin-coated on a highly ordered pyrolytic graphite substrate.

Samples of ca. 50-μm thickness for frequency-doubled light emission, or second harmonic generation (SHG) measurements were prepared on the rough side of silicon wafer and examined using a Zeiss LSM 510 Meta upright laser-scanning confocal microscope with a ×40/0.75NA objective. The temperature was controlled by a Linkam hot stage. A Chameleon Ti:Sapphire femtosecond pulsed laser (Coherent, California), tuned to 800 nm, was attached to the microscope and focused onto the sample resulting in a SHG signal detectable at 400 nm. In order to exclude possible surface effects, first the beam was focused at the bottom of the LC film in contact with the Si substrate, then at 20 μm height. The SHG ratio was calculated by dividing the intensities at 400 nm by the averaged background around the peak (see also ref. 32).

For UV–vis and fluorescence emission spectroscopy, the samples were melted and spread evenly between quartz plates (1 mm thickness). UV–vis spectra were recorded using Lambda 900 (Perkin Elmer). Fluorescence spectra of FO16 were excited at 420 nm and recorded by Fluoromax4 (Horiba).

DSC thermograms were recorded on a DSC 200 F3 Maia calorimeter (NETZSCH) with heating/cooling rates as specified.

Polarized optical micrographs were recorded using an Olympus BX-50 equipped with a Mettler HS82 hot stage.

## Code availability

All derived mathematical equations are given in Supplementary Information and no special code was required in their use.

## References

1. Frank, F. C. On spontaneous asymmetric synthesis. Biochim. Biophys. Acta 11, 459–463 (1953).

2. Link, D. R. et al. Spontaneous formation of macroscopic chiral domains in a fluid smectic phase of achiral molecules. Science 278, 1924–1927 (1997).

3. Dressel, C. et al. Dynamic mirror-symmetry breaking in bicontinuous cubic phases. Angew. Chem. Int. Ed. 53, 13115–13120 (2014).

4. Zeng, X. B., Ungar, G. & Imperor-Clerc, M. A triple-network tricontinuous cubic liquid crystal. Nat. Mater. 4, 562–567 (2005).

5. Zeng, X. B. & Ungar, G. Spontaneously chiral cubic liquid crystal: three interpenetrating networks with a twist. J. Mater. Chem. C 8, 5389–5398 (2020).

6. Lu, H. J. et al. The solution of the puzzle of Smectic-Q: the phase structure and the origin of spontaneous chirality. Angew. Chem. Int. Ed. 57, 2835–2840 (2018).

7. Kutsumizu, S. Recent progress in the synthesis and structural clarification of thermotropic cubic phases. Isr. J. Chem. 52, 844–853 (2012).

8. Oswald, P. & Pieranski, P. Smectic and Columnar Liquid Crystals: Concepts and Physical Properties Illustrated by Experiments (CRC Press, 2005).

9. Sun, H. J., Zhang, S. D. & Percec, V. From structure to function via complex supramolecular dendrimer systems. Chem. Soc. Rev. 44, 3900–3923 (2015).

10. Vera, F., Serrano, J. L. & Sierra, T. Twists in mesomorphic columnar supramolecular assemblies. Chem. Soc. Rev. 38, 781–796 (2009).

11. Wohrle, T. et al. Discotic liquid crystals. Chem. Rev. 116, 1139–1241 (2016).

12. Kato, T., Yasuda, T., Kamikawa, Y. & Yoshio, M. Self-assembly of functional columnar liquid crystals. Chem. Commun. 7, 729–739 (2009).

13. Fall, W. S. et al. Molecular ejection transition in liquid crystal columns self-assembled from wedge-shaped minidendrons. Soft Matter 15, 22–29 (2019).

14. Li, Y. X. et al. New type of columnar liquid crystal superlattice in double-taper ionic minidendrons. Chem. Eur. J. 25, 13739–13747 (2019).

15. Cheng, X. et al. Calamitic bolaamphiphiles with (semi) perfluorinated lateral chains: polyphilic block molecules with new liquid crystalline phase structures. J. Am. Chem. Soc. 125, 10977–10996 (2003).

16. Chen, B. et al. Liquid crystalline networks composed of pentagonal, square, and triangular cylinders. Science 307, 96–99 (2005).

17. Ungar, G. et al. Self-assembly at different length scales: polyphilic star-branched liquid crystals and miktoarm star copolymers. Adv. Funct. Mater. 21, 1296–1323 (2011).

18. Destrade, C. et al. Disc-like mesogen polymorphism. Mol. Cryst. Liq. Cryst. 106, 121–146 (1984).

19. Gharbi, M., Gharbi, A., Nguyen, H. T. & Malthete, J. Polycatenar liquid crystals with long rigid aromatic cores: a review of recent works. Curr. Opin. Colloid Interface Sci. 7, 312–325 (2002).

20. Gainar, A., Tzeng, M. C., Heinrich, B., Donnio, B. & Bruce, D. W. Incompatibility-driven self-organization in polycatenar liquid crystals bearing both hydrocarbon and fluorocarbon chains. J. Phys. Chem. B. 121, 8817–8882 (2017).

21. Wöhrle, T. et al. Discotic liquid crystals. Chem. Rev. 116, 1139–1241 (2016).

22. Fontes, E., Heiney, P. A. & de Jeu, W. H. Liquid-crystalline and helical order in a discotic mesophase. Phys. Rev. Lett. 61, 1202–1205 (1988).

23. Metzroth, T. et al. Unravelling the fine structure of stacked bipyridine diamine-derived C3-discotics as determined by X-ray diffraction, quantum-chemical calculations, Fast-MAS NMR and CD spectroscopy. Chem. Sci. 2, 69–76 (2011).

24. Zagrovic, B. Helical signature motif in the fibre diffraction patterns of random-walk chains. Mol. Phys. 105, 1299–1306 (2007).

25. Baumgarten, J. L. Ferrochirality: a simple theoretical model of interacting, dynamically invertible, helical polymers, 2. molecular field approach: supports and the details. Macromol. Theory Simul. 4, 1–43 (1995).

26. Wang, L. et al. Extraordinary acceleration of cogwheel helical self-organization of dendronized perylene bisimides by the dendron sequence encoding their tertiary structure. J. Am. Chem. Soc. 142, 9525–9536 (2020).

27. Heiney, P. A. et al. Frustration and helicity in the ordered phases of a discotic compound. J. Physique 50, 461–483 (1989).

28. Corradini, P., Auriemma, F. & De Rosa, C. Crystals and crystallinity in polymeric materials. Acc. Chem. Res. 39, 314–323 (2006).

29. Clark, E. S. & Muus, L. T. Partial disordering and crystal transitions in polytetrafluoroethylene. Z. Kristallogr. 117, 119–127 (1962).

30. Ungar, G. Thermotropic hexagonal phases in polymers: common features and classification. Polymer 34, 2050–2059 (1993).

31. Savin, A. V. & Manevitch, L. I. Topological solitons in spiral polymeric macromolecules: a chain surrounded by immovable neighbors. Phys. Rev. B 63, 224303 (2001).

32. Cheng, H. F. et al. Trigonal columnar self-assembly of bent phasmid mesogens. Chem. Commun. 54, 156–159 (2018).

33. Kleman, M. Developable domains in hexagonal liquid-crystals. J. Physique 41, 737–745 (1980).

34. Fall, W. S. et al. Quantised molecular ejection transition in liquid crystal columns self-assembled from wedge-shaped minidendrons. Soft Matter 15, 22–29 (2019).

35. Yen, M. H. et al. Added alkane allows thermal thinning of supramolecular columns by forming superlattice—an X-ray and Neutron Study. J. Am. Chem. Soc. 138, 5757–5760 (2016).

36. Shcherbina, M. A. et al. Hollow six-stranded helical columns of a helicene. Angew. Chem. Int. Ed. 48, 7837–7840 (2009).

37. Feng, X. et al. Towards high charge-carrier mobilities by rational design of the shape and periphery of discotics. Nat. Mater. 8, 421–426 (2009).

38. Percec, V. et al. Self-assembly of dendronized perylene bisimides into complex helical columns forming thermodynamically controlled 2D and kinetically controlled 3D periodic arrays. J. Am. Chem. Soc. 133, 12197–12219 (2011).

39. Roche, C. et al. Homochiral columns by chiral self-sorting during supramolecular helical organisation of hat-shaped molecules”. J. Am. Chem. Soc. 136, 7169–7185 (2014).

40. Takezoe, H. Polar liquid crystals—ferro, antiferro, banana, and columnar. Mol. Cryst. Liq. Cryst. 646, 46–65 (2017).

41. Meuler, A. J., Hillmyer, M. A. & Bates, F. S. Ordered network mesostructures in block polymer materials. Macromolecules 42, 7221–7250 (2009).

42. Tyler, C. A. & Morse, D. C. Orthorhombic Fddd network in triblock and diblock copolymer melts. Phys. Rev. Lett. 94, 208–302 (2005).

43. Wang, X.-B., Lo, T.-Y., Hsueh, H.-Y. & Ho, R.-M. Double and single network phases in polystyrene-block-poly(l-lactide) Diblock copolymers. Macromolecules 46, 2997–3004 (2013).

44. Chakraborty, C. et al. Vice versa donor acceptor fluorene–ferrocene alternate copolymer: a twisted ribbon for electrical switching. Chem. Commun. 51, 13123–13126 (2015).

45. Kim, D. et al. Solution-processable fluorene derivative for organic thin-film transistors. Org. Electron. 76, 105464 (2020).

46. Hansch, C., Leo, A. & Taft, R. W. A survey of Hammett substituent constants and resonance and field parameters. Chem. Rev. 91, 165–195 (1991).

47. Nuckolls, C. et al. Synthesis and aggregation of a conjugated helical molecule. J. Am. Chem. Soc. 121, 79–88 (1999).

48. Ravetz, B. D. et al. Photoredox catalysis using infrared light via triplet fusion upconversion. Nature 565, 343–346 (2019).

49. Tauc, J., Grigorovici, R. & Vancu, A. Optical properties and electronic structure of amorphous germanium. Phys. Stat. Solidi (b) 15, 627–637 (1966).

## Acknowledgements

For support with experiments at Diamond Light Source, we thank Dr. Olga Shebanova and Prof. Nick Terill (I22) and Dr. Gareth Nisbet and Prof. Steve Collins (I16), as well as Drs. Oier Bikondoa and Paul Thompson at XMaS beamline (BM28) at ESRF. We thank Prof. Julia Weinstein for help with fluorescence spectroscopy. For performing mass spectrometry, we also thank the application scientist Qinghe Wang from Bruker Daltonics China. Financial support is acknowledged from the 111 Project 2.0 of China (BP2018008), from EPSRC (EP-P002250, EP-T003294), from NSFC (21865038), and Yunnan Joint Fund Key Project (2018FY001(−012)). Y.X.L. thanks CSC for the stipend and UoS for a fee waiver.

## Author information

Authors

### Contributions

X.H.C., X.B.Z., and G.U. supervised the study. H.F.G. and Q.C. synthesized and characterized the compounds chemically, H.F.G. performed DFT calculations. Y.X.L. carried out X-ray, UV–vis, fluorescence experiments, reconstructed ED maps, and performed MD simulation. R.B.Z. and Y.X.L. performed AFM experiments. K.G., G.A.G., and X.B.Z. developed the packing theory for helices. G.U., X.B.Z., and Y.X.L. wrote the paper with contributions from all authors. All authors discussed the results and commented on the manuscript.

### Corresponding authors

Correspondence to Xiao-hong Cheng, Xiang-bing Zeng or Goran Ungar.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Peer review

### Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Li, Yx., Gao, Hf., Zhang, Rb. et al. A case of antiferrochirality in a liquid crystal phase of counter-rotating staircases. Nat Commun 13, 384 (2022). https://doi.org/10.1038/s41467-022-28024-1

• Accepted:

• Published:

• DOI: https://doi.org/10.1038/s41467-022-28024-1