
2
II. COARSE-GRAINED FREE ENERGY FOR
FERRONEMATICS
FNs are described in terms of two order parameters:
(i) the Q-tensor, which contains information about the
orientational order of the NLCs, and (ii) the magneti-
zation vector M, which gives the average orientation of
the magnetic moments of the MNPs. The Q-tensor is
symmetric and traceless, and is given by [19]:
Qij =Sninj+T kikj−(S+˜
T)δij
3.(1)
Here, the scalar order parameter Smeasures the uniaxial
degree of order about the leading eigenvector or the di-
rector n. Further, ˜
Tis the magnitude of the biaxial order
about the secondary director k. (A system with only uni-
axial order has ˜
T= 0. For such a system, the isotropic
phase corresponds to S= 0, and the nematic phase has
S 6= 0.) Taking into account the requirements of symme-
try and tracelessness, the Q-tensor can be expressed in
terms of five independent parameters as follows:
Q=
−q1+q2q3q4
q3−q1−q2q5
q4q52q1
.(2)
To obtain the nematic directors and S,˜
T, we choose
a frame of reference in which Qis diagonal. This pro-
vides us the three eigenvalues (λ3> λ2> λ1), and the
corresponding eigenvectors n,k,l. The largest eigen-
value λ3=S, and the corresponding eigenvector is the
primary direction of order n[19, 20]. We will use a
standard measure of biaxial order about the secondary
director k:T= (λ2−λ1)/λ3[7, 8, 20], which is pro-
portional to ˜
T. Naturally, λ1=λ2if the system is uni-
axial. The degree of biaxiality can also be defined as
B2={1−6Tr(Q3)2/[Tr(Q2)3]}[21, 22], where B2= 0
for the uniaxial state and B2= 1 for a state with maxi-
mum biaxiality. This definition of biaxiality also exploits
the difference between two eigenvalues to determine bi-
axial order, similar to T.
We use the Landau-de Gennes (LdG) approach to write
down the phenomenological free energy for this compos-
ite system. This is a functional of the order parameter
fields Q(r) and M(r) and has three contributions [23–28]:
G[Q,M] = ZdrA
2Tr(Q2) + C
3Tr(Q3) + B
4[Tr(Q2)]2
+L
2|∇Q|2+α
2|M|2+β
4|M|4+κ
2|∇M|2
−γµ0
2
3
X
i,j=1
Qij MiMj
.(3)
The first four terms in Eq. (3) represent the Ginzburg-
Landau (GL) free energy for the nematic component with
Landau coefficients A,B,C,Lhaving their usual mean-
ing. The next three terms correspond to the GL free
energy for the magnetic component. In the GL frame-
work, the gradient terms |∇Q|2and |∇M|2are essen-
tial to capture the effects of elastic interactions [29–33].
They penalise local variations in the order parameters
– this surface tension results in the motion of domain
boundaries in coarsening kinetics.
The magnitudes of the Landau coefficients determine
the scales of order parameter, length and time in the sys-
tem. For example, A=A0(T−TN) and α=α0(T−TM)
depend on the quench temperature Tand the critical
temperatures TN,TM. (Here, A0,α0are material-
dependent constants.) A direct estimate of the coeffi-
cients can be obtained from experimentally determined
quantities like the latent heat, order parameter magni-
tudes, susceptibilities, etc. [30, 33]. However, the cur-
rent experimental data on FNs is not adequate to pro-
vide accurate estimates of these coefficients. The utility
of the LdG framework lies primarily in predicting uni-
versal behaviors, e.g., power laws and their exponents,
scaling variables, etc.
The effect of dopant particles in LCs has been mod-
eled in several previous studies [34–37]. These models
describe the coupling of the dipole moment of ferroelec-
tric particles with the NLCs at a molecular level. The
induced field due to the impurity atoms acts like an
aligning field, and enhances orientational order in the
NLCs. On a similar footing, the last term in Eq. (3)
is the phenomenological magneto-nematic coupling de-
fined as a dyadic product of Qand Mand the parame-
ter γis the strength of the coupling. It is related to the
shape and size of the MNPs and their interaction with the
NLCs. This cubic magneto-nematic coupling term [24]
enforces the specific orientations of the magnetic and ne-
matic components essential for the emergence of biaxial
order in the system [15, 38]. A more accurate description
of the free energy can be obtained by incorporating dipo-
lar and quadrupolar interactions. This may be required
for studies of phase transitions and critical phenomena.
As discussed in Ref. [26], these terms may be ignored for
dilute ferronematic suspensions.
In their experiments, Liu et al. demonstrated that bi-
axial order emerges only when nand Mare tilted at
an angle. By manipulating the surface functionaliza-
tion, they could achieve a tilt angle up to 90◦. (Their
optical absorbance measurements to detect the biaxial
phase were carried out for a limited range from 10◦-65◦.)
Motivated by these experiments, we choose γ < 0 for
simplicity, which corresponds to a tilt angle of 90◦. In
principle, it is possible to modify the coupling term in
Eq. (3) such that nand Mare at an arbitrary angle,
but this makes the expression considerably more compli-
cated. The emergence of biaxiality (or the presence of
two distinguished directions) in the NLCs for non-zero
values of γ < 0 can be understood from the schematic
in Fig. 1: Choosing Malong the positive x-axis, the
LC molecules can align in two orthogonal directions, say
along the y-axis and z-axis.