Hexagonal Bond-Orientational Order (HBOO)

• Hexagonal bond‐orientational order (HBOO) is defined by sixfold symmetry in particle bonding, separating orientational alignment from translational order.
• It is characterized using order parameters like Ψ6 and Q6 and measured via correlation functions that exhibit algebraic decay in the hexatic phase.
• HBOO informs studies in liquid crystals, colloids, nucleation, and supercooled liquids through both experimental techniques (e.g., X-ray Fourier analysis) and simulation.

Searching arXiv for recent and foundational papers on hexagonal bond-orientational order. Search query: "hexagonal bond orientational order hexatic smectic XCCA Steinhardt SymBOP" Hexagonal bond-orientational order (HBOO) denotes sixfold orientational symmetry in the arrangement of bonds linking a particle or molecule to its neighbors. In the literature on two-dimensional melting, liquid crystals, colloids, nanoparticle lattices, and supercooled liquids, HBOO is the orientational counterpart to positional order: a system may retain liquid-like or short-range translational correlations while sustaining short-range, quasi-long-range, or long-range sixfold bond correlations. In two dimensions this distinction is central to the hexatic phase, for which the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario predicts algebraically decaying orientational correlations and exponentially decaying positional correlations; in three-dimensional layered systems and other complex materials, HBOO is analyzed through local $\Psi_6$-type fields, spherical-harmonic invariants such as $Q_6$, X-ray Fourier harmonics $C_{6m}$, and symmetry-specific bond-order projections (Agosta et al., 2017, Tanaka, 2013).

1. Definitions and order parameters

The standard two-dimensional hexatic order parameter is the local complex field

$$\Psi_6 = \frac{1}{N_b} \sum_{j=1}^{N_b} \exp(i6\theta_{ij}),$$

where $\theta_{ij}$ is the angle of the bond from particle $i$ to neighbor $j$, and $N_b$ is the number of neighbors. This representation makes explicit that HBOO is a sixfold orientational quantity: its magnitude measures the local degree of hexagonal alignment, while its phase encodes bond orientation (Tanaka, 2013).

In three dimensions, bond-orientational order is commonly represented with the Steinhardt order parameters. For a particle $i$,

$$q_{\ell m}(i) = \frac{1}{N_b(i)} \sum_{j=1}^{N_b(i)} Y_{\ell m}(\hat{\mathbf{r}_{ij}}),$$

and the associated rotational invariant is

$Q_6$0

The $Q_6$1 channel, $Q_6$2, is the conventional scalar measure most closely associated with hexagonal or close-packed orientational environments in three-dimensional particulate systems (Russo et al., 2012). A closely related rotationally invariant form is

$Q_6$3

and its averaged Lechner-Dellago variant

$Q_6$4

These definitions are used to distinguish liquid-like, crystal-like, cubic, hexagonal, and hydrate-like local environments in simulations (Zerón et al., 2024).

For liquid-crystal hexatics, HBOO is also described by a coarse-grained order field

$Q_6$5

where $Q_6$6 is the local bond angle relative to a reference axis. This notation is especially common in descriptions of the hexatic-B phase and in analyses of scaling among higher sixfold harmonics (Zaluzhnyy et al., 2017).

2. Correlation functions, Fourier harmonics, and scaling laws

The defining observable for HBOO is the bond-orientational correlation function. In a particle-based two-dimensional description it is written as

$Q_6$7

or, in a layer-resolved form used for simulated smectics,

$Q_6$8

In the KTHNY framework, $Q_6$9 decays algebraically in the hexatic phase,

$C_{6m}$0

with $C_{6m}$1, whereas positional correlations decay exponentially (Agosta et al., 2017).

In X-ray studies of liquid-crystal hexatics, the azimuthal intensity on the diffraction ring is expanded as

$C_{6m}$2

and for sixfold symmetry only harmonics with $C_{6m}$3 contribute significantly. The normalized bond-orientational order parameters are then

$C_{6m}$4

with $C_{6m}$5 the peak wavevector (Zaluzhnyy et al., 2014).

Angular X-ray cross-correlation analysis (XCCA) extracts these quantities through the two-point angular correlation function

$C_{6m}$6

whose Fourier coefficients satisfy

$C_{6m}$7

This route avoids assuming a particular line shape for the azimuthal intensity and was used to determine successive sixfold harmonics directly in freely suspended hexatic films (Zaluzhnyy et al., 2017).

A central scaling result for three-dimensional hexatics is the multicritical scaling theory (MCST) relation

$C_{6m}$8

with

$C_{6m}$9

The parameter $$\Psi_6 = \frac{1}{N_b} \sum_{j=1}^{N_b} \exp(i6\theta_{ij}),$$0 is reported as $$\Psi_6 = \frac{1}{N_b} \sum_{j=1}^{N_b} \exp(i6\theta_{ij}),$$1 for 3D hexatics and $$\Psi_6 = \frac{1}{N_b} \sum_{j=1}^{N_b} \exp(i6\theta_{ij}),$$2 for 2D hexatics, while the next correction in 3D is $$\Psi_6 = \frac{1}{N_b} \sum_{j=1}^{N_b} \exp(i6\theta_{ij}),$$3 (Zaluzhnyy et al., 2017). This formalism organizes HBOO as a hierarchy of coupled sixfold harmonics rather than a single scalar.

3. Hexatic phases, quasi-long-range order, and layered systems

In two dimensions, the hexatic phase is characterized by quasi-long-range orientational order and short-range positional order. A notable three-dimensional realization of this two-dimensional phenomenology was reported in a molecular-dynamics simulation of a smectic phase formed by a one-component system of particles interacting via a spherically symmetric pair potential. In that system, the in-layer bond-orientational correlations decayed algebraically as

$$\Psi_6 = \frac{1}{N_b} \sum_{j=1}^{N_b} \exp(i6\theta_{ij}),$$4

in quantitative agreement with the KTHNY prediction for the hexatic-isotropic transition in two dimensions (Agosta et al., 2017).

The same study emphasized that previously observed three-dimensional hexatic smectics exhibited truly long-range in-layer bond-orientational order, rather than algebraic decay. The contrast had remained unexplained and had sometimes been heuristically attributed to interlayer interactions or molecular anisotropy. The simulated layered phase therefore provided the first observation of theoretically predicted two-dimensional hexatic order in a three-dimensional system, and it did so with isotropic interactions rather than particle anisotropy or explicitly directional forces (Agosta et al., 2017).

The simulated system contained 50,000 identical particles in a cubic box with periodic boundaries, at fixed density $$\Psi_6 = \frac{1}{N_b} \sum_{j=1}^{N_b} \exp(i6\theta_{ij}),$$5. Upon cooling from an isotropic liquid, a first-order transition to a hexatic smectic phase was observed, followed at lower temperature by a transition to a smectic B crystal. The hexatic smectic retained liquid-like in-layer diffusion and short-range positional order, as indicated by an exponentially decaying radial distribution function $$\Psi_6 = \frac{1}{N_b} \sum_{j=1}^{N_b} \exp(i6\theta_{ij}),$$6 (Agosta et al., 2017).

Experimental liquid-crystal hexatics display a related, but not identical, layered phenomenology. In freely suspended films of 75OBC and 3(10)OBC, the angular profile of a single diffraction peak is close to Gaussian near the hexatic-B–smectic-A transition and is better fitted by a Voigt function at lower temperatures in the hexatic-B phase. The effective Hamiltonian introduced by Aharony and Kardar predicts finite angular fluctuations in many coupled layers, and the corresponding theoretical estimates were reported to be in good agreement with the X-ray data (Zaluzhnyy et al., 2016). This suggests a distinction between stacked three-dimensional hexatics stabilized by interlayer coupling and quasi-two-dimensional hexatic order realized within layers.

4. Experimental quantification in liquid-crystal films

Spatially resolved X-ray studies of freely suspended hexatic films of 3(10)OBC revealed substantial spatial inhomogeneity of bond-orientational order near the hexatic-smectic transition and the formation of large-scale hexatic domains at lower temperatures. At the transition, the magnitude of the sixth-order Fourier component $$\Psi_6 = \frac{1}{N_b} \sum_{j=1}^{N_b} \exp(i6\theta_{ij}),$$7 varied strongly across the sample, indicating coexistence of regions with stronger and weaker bond-orientational order. Upon cooling to $$\Psi_6 = \frac{1}{N_b} \sum_{j=1}^{N_b} \exp(i6\theta_{ij}),$$8, the system formed extended single hexatic domains with uniform bond orientation over areas up to at least $$\Psi_6 = \frac{1}{N_b} \sum_{j=1}^{N_b} \exp(i6\theta_{ij}),$$9 (Zaluzhnyy et al., 2014).

A distinctive result of that XCCA study was the direct determination of up to 25 successive sixfold bond-orientational order parameters $\theta_{ij}$0 deep in the hexatic phase. This unusually large number of resolved harmonics allowed a stringent test of MCST across the full temperature range of hexatic stability. The analysis confirmed the scaling law $\theta_{ij}$1 and directly determined the first-order correction term, with $\theta_{ij}$2 and $\theta_{ij}$3, consistent with theory. At the highest orders, $\theta_{ij}$4, a small but systematic deviation required a second-order correction $\theta_{ij}$5 (Zaluzhnyy et al., 2014).

A related X-ray study on three liquid-crystal compounds—3(10)OBC, 75OBC, and PIRO6—reported that the MCST description remains valid over the whole temperature range of the hexatic phase for all three compounds. Near the smectic-hexatic transition, the temperature dependence of the first bond-orientational harmonic follows a conventional power law with an unusually small critical exponent,

$\theta_{ij}$6

Higher harmonics scale as powers of the first harmonic,

$\theta_{ij}$7

which was interpreted as evidence for nonlinear coupling among bond-orientational order parameters of different order (Zaluzhnyy et al., 2017).

The same study reported compound-dependent but closely clustered MCST parameters, including $\theta_{ij}$8 for 3(10)OBC, $\theta_{ij}$9 for 75OBC, and $i$0 for PIRO6. This supports the claim, made explicitly there, that compounds of various composition and with different low-temperature crystalline phases display the same thermodynamic behavior in the hexatic phase and in the vicinity of the smectic-hexatic transition (Zaluzhnyy et al., 2017).

5. Methodological refinements and symmetry-specific generalizations

A longstanding methodological issue in HBOO analysis is the dependence of conventional $i$1 parameters on the definition of the local neighborhood. For a particle $i$2,

$i$3

It was argued that the procedure used to select $i$4 can influence both the numerical values and the qualitative trend of $i$5 more strongly than changes in physical control parameters such as packing fraction, and that the discrete nature of $i$6 makes $i$7 discontinuous with respect to particle coordinates (Mickel et al., 2012).

To remedy these shortcomings, a morphometric or Minkowski version was introduced: $i$8 Here the contributions are weighted by Voronoi facet areas rather than counted uniformly. This construction is parameter-free, continuous under smooth particle displacements, and mathematically equivalent to Minkowski tensors in the cited formulation (Mickel et al., 2012).

A related refinement addresses robustness under noise. Using Voronoi-facet weighting with exponent $i$9,

$j$0

it was reported that $j$1 produces robust and continuous bond-orientational order parameters that distinguish noisy FCC and BCC environments much more reliably than threshold, Delaunay, or $j$2 schemes. Since $j$3 is the symmetry channel relevant to hexagonal arrangements, this refinement directly affects HBOO detection in noisy soft-matter data (Haeberle et al., 2019).

A more symmetry-selective development is the Symmetry-specific Bond Order Parameter framework. In two dimensions, $j$4 remains the reference object for sixfold order; in three dimensions, standard multipoles are projected onto symmetry-invariant subspaces: $j$5 These SymBOPs and Symmetrized Bond Order Parameters can be assigned locally to individual bonds or averaged globally, and they preserve orientational information that is absent from rotational invariants such as $j$6. They were shown to identify coherent crystalline domains with different orientations, detect dislocations and grain boundaries, and distinguish individual sublattices such as interpenetrating FCC lattices within a cubic diamond (Logan et al., 2023, Logan et al., 2021). A plausible implication is that symmetry-specific projection is especially valuable when HBOO must be separated from other orientational motifs in polycrystalline or defect-rich samples.

Rotationally invariant local bond order parameters have also been revisited for hydrate structure identification. In systems where clathrate hydrates coexist with liquid water, the combination $j$7 and $j$8 was reported to outperform the more conventional $j$9-$N_b$0 pairing for distinguishing hydrate-like from liquid-like molecules across carbon dioxide, methane, tetrahydrofuran, nitrogen, and hydrogen hydrates (Zerón et al., 2024). This is not a statement about HBOO alone, but it shows that the $N_b$1 channel is not universally optimal once local geometry departs from the close-packed environments for which $N_b$2 is most commonly used.

6. Roles of HBOO in nucleation, slow dynamics, and emergent ordered states

In hard-sphere crystallization, bond-orientational ordering was reported to precede positional ordering. Computer simulations of a metastable hard-sphere melt found no evidence for a two-step mechanism in which crystals form inside dense precursor regions. Instead, nucleation was driven by fluctuations of orientational order rather than density, and $N_b$3 fluctuations were spatially extended while density fluctuations remained local. The density at the nucleus center extrapolated to the fluid density as nucleus size went to zero, which was interpreted as evidence for the absence of dense precursors (Russo et al., 2012). This places HBOO at the center of a controversy over whether crystallization in simple fluids is density driven or symmetry driven.

An independent dynamical controversy concerns supercooled water. In simulations of ST2, TIP5P, TIP4P/2005, and mW water, the Steinhardt-Nelson-Ronchetti $N_b$4 order parameter was used to monitor bond-orientational relaxation near reported liquid-liquid-transition conditions. Across all models and state points considered, density relaxed more slowly than or at least as slowly as bond-orientational order, contradicting the artificial polyamorphism hypothesis, which requires $N_b$5 (Palmer et al., 2015). In this context, HBOO is not merely a structural classifier but a dynamical observable used to test whether apparent polyamorphism could instead be a crystallization artifact.

HBOO also appears in electronically ordered solids. In the magnetic kagome metal FeGe, a dimerization-driven two-dimensional hexagonal charge-diffuse precursor was reported above the multiple-$N_b$6 charge density wave transition. In the regime $N_b$7, the diffuse scattering around the $N_b$8 wavevector is anisotropic, indicating quasi-long-range bond-orientational order. The local order field was defined as

$N_b$9

with correlation function

$i$0

In the ordered phase $i$1 is constant, while in the intermediate regime it decays algebraically, paralleling the phenomenology of liquid-crystal hexatics and defect-mediated melting (Subires et al., 2024).

Finally, recent work on two-dimensional colloids linked HBOO to dynamical heterogeneity as reported by dumbbell probes. Across the liquid-hexatic transition, the host system’s rotationally resolved probe dynamics changed from Brownian in the isotropic liquid to non-Gaussian in the hexatic and solid phases. In mobile domains, the dumbbells undergo rotational jumps of $i$2, consistent with the sixfold symmetry of HBOO; in immobile domains, they librate within cages formed by surrounding discs. The disappearance of this non-Gaussianity upon reentrant melting driven by size polydispersity was presented as evidence for a close connection between HBOO and probe dynamics (Kim et al., 2 Oct 2025). This suggests that HBOO can be inferred not only from static bond geometry but also from symmetry-locked transport signatures.

HBOO therefore occupies several distinct but connected roles: it is the defining order of the hexatic phase, a measurable field in layered liquid crystals, a precursor variable in crystallization, a test variable in supercooled-liquid dynamics, a symmetry marker in nanoparticle and hydrate analysis, and, in frustrated electronic materials, an intermediate order that survives after translational coherence is lost. The common feature across these settings is sixfold bond alignment decoupled, to varying degrees, from full positional crystallinity (Agosta et al., 2017, Zaluzhnyy et al., 2014).