Scholarly article on topic 'Spontaneous liquid crystal and ferromagnetic ordering of colloidal magnetic nanoplates'

Spontaneous liquid crystal and ferromagnetic ordering of colloidal magnetic nanoplates Academic research paper on "Nano-technology"

Share paper
Academic journal
Nature Communications
OECD Field of science

Academic research paper on topic "Spontaneous liquid crystal and ferromagnetic ordering of colloidal magnetic nanoplates"



Received 23 Oct 2015 | Accepted 8 Dec 2015 | Published 28 Jan 2016^^HDOI:iai038/ncomms1039J OPEN

Spontaneous liquid crystal and ferromagnetic ordering of colloidal magnetic nanoplates

M. Shuai1, A. Klittnick1, Y. Shen1, G.P. Smith1, M.R. Tuchband1, C. Zhu2, R.G. Petschek3, A. Mertelj4, D. Lisjak4, M. Copic4'5, J.E. Maclennan1, M.A. Glaser1 & N.A. Clark1

Ferrofluids are familiar as colloidal suspensions of ferromagnetic nanoparticles in aqueous or organic solvents. The dispersed particles are randomly oriented but their moments become aligned if a magnetic field is applied, producing a variety of exotic and useful magnetomechanical effects. A longstanding interest and challenge has been to make such suspensions macroscopically ferromagnetic, that is having uniform magnetic alignment in the absence of a field. Here we report a fluid suspension of magnetic nanoplates that spontaneously aligns into an equilibrium nematic liquid crystal phase that is also macroscopically ferromagnetic. Its zero-field magnetization produces distinctive magnetic self-interaction effects, including liquid crystal textures of fluid block domains arranged in closed flux loops, and makes this phase highly sensitive, with it dramatically changing shape even in the Earth's magnetic field.

1 Department of Physics and Soft Materials Research Center, University of Colorado, Boulder, Colorado 80309, USA. 2 Advanced Light Source, Lawrence

Berkeley National Laboratory, Berkeley, California 94720, USA. 3 Department of Physics, Case Western Reserve University, Cleveland, Ohio 44106, USA.

4Jozef Stefan Institute, SI-1000 Ljubljana, Slovenia. 5Faculty of Mathematics and Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia.

Correspondence and requests for materials should be addressed to N.A.C. (email:

Colloidal suspensions of nanoparticles offer exciting and ever-expanding opportunities for creation and application of novel hybrid particle/fluid phenomena. A classic example is the paramagnetic ferrofluid, a colloidal suspension of ferromagnetic nanoparticles dispersed in isotropic solvent1'2. The ferrofluid particles disperse with random orientations but readily align and move in magnetic fields, resulting in dazzling hybrid behaviours, including radical changes in shape3'4 and magnetically controlled transport, which enable a variety of technical and biomedical applications5'6. A longstanding challenge has been to make fluids that are ferromagnetic, such that they exhibit spontaneous, equilibrium magnetic ordering in the absence of an applied field7-11. Such phases are conceivable because in such hybrid suspensions the dipole-dipole interaction energy between the supermolecular magnetic moments of neighbouring particles can be many times the thermal energy at room temperature10,11.

One approach to achieve ferromagnetism in colloidal fluids is to create stable, fluid suspensions of well-dispersed magnetic nanoparticles in isotropic solvent, and to design their mutual interactions to produce equilibrium, zero-field magnetization. Since ferromagnetic ordering implies orientational ordering of the particles, such fluid suspensions are also liquid crystals (LCs). Suspensions of magnetic nanoplates with moments normal to their planes are of particular interest to this study. One reason is that anisotropic steric interactions can lead to spontaneous LC ordering12, which have been demonstrated in various colloidal suspensions of nanoplates and nanosheets13,14. Meanwhile, the plate-like (as opposed to spherical or rod-like) shape promotes parallel rather than antiparallel orientation of the moments of neighbouring particles, making the interplate magnetic dipole interactions sufficient to stabilize ferromagnetic ordering, as recently demonstrated in dilute suspensions of nanoplates orientationally ordered by a thermotropic nematic LC host


In our experiments, barium hexaferrite (BF) nanoplates are suspended in isotropic n-butanol and surfactant-stabilized to produce a system of functionalized nanoplates with weak electrostatic repulsion, strong and anisotropic steric repulsion and magnetic interaction. The electrostatic repulsion enables stable suspensions at essentially any concentration, including at high volume fractions where spontaneous LC ordering may occur. We observe nematic ordering of the nanosheets at volume fractions in the range where the Onsager model and subsequent theory and simulations predict an isotropic (Iso)-nematic transition for polydisperse, plate-like colloids119. This nematic is distinctly ferromagnetic, with a magnetization density two orders of magnitude larger than that reported in the previously found thermotropic nematic LC host (4-cyano-40-pentylbiphenyl, 5CB) fluid system15.


Isotropic to ferromagnetic nematic phase transition.

Suspensions of BF nanoplates (overall thickness t = 7 nm, diameter D = 48 ± 21 nm, Supplementary Fig. 1) in n-butanol (BF/ BuOH) are loaded in capillaries and observed at various magnifications using both unpolarized and depolarized transmitted light microscopy, which enabled visualization of their basic phase behaviour as a function of the volume fraction of the functionalized nanoplates, f = p<D2t>v/4, where v is the number density (Fig. 1a). For f <0.28, the suspensions are optically isotropic, giving excellent extinction for any sample orientation between crossed polarizer and analyser, and uniform transmission with only the polarizer, indicating that the particles are well dispersed and that the suspensions are optically

Figure 1 | Liquid crystal ordering of barium ferrite nanoplatelets suspended in n-butanol. (a) Barium hexaferrite nanoplatelet suspensions viewed in transmitted light with optical polarization conditions indicated (Polarizer: magenta, P; analyser: cyan, A). Low-volume fraction suspensions (f <0.25) are isotropic (Iso), appearing dark between crossed polarizers. The orange/red colour is due to optical absorption by the nanoplatelets. At higher concentrations (f >0.28), a birefringent ferromagnetic nematic (NF) phase appears in the lower part of the cell. (b) An applied in-plane magnetic field Bext induces birefringence in the isotropic phase, with the principal axes of the optical dielectric tensor along and normal to Bext and the induced macroscopic magnetization density M parallel to Bext. (c) The Nf phase is separated gravitationally from the isotropic region by a sharp, horizontal interface. Equilibrium Iso and NF structures deduced from birefringence and dichroism measurements are illustrated, where mo is the nanoplate magnetic moment and n the director, indicating the local, mean nanoplate normal. (d) The Iso phase is magnetized and the Iso-NF interface becomes continuous under 0.5-mT applied magnetic field. Samples are sealed in rectangular glass capillaries of thickness L = 50 mm and width W = 1 mm. The boundaries of the cells are indicated by the solid thin white lines. The air-liquid interfaces are indicated by the dashed yellow lines.

homogeneous and isotropic on micron and larger scales. The characteristic orange-red colour is due to nanoplate absorption (Supplementary Fig. 2). The isotropic suspensions are stable, with only weak gravitational concentration gradients developing over periods of hundreds of days (Supplementary Fig. 3), consistent with an estimated Perrin length /P = kBT/ApLS(K<D2t>/4)g~ 1 cm (DpLS = mass density difference between nanoplates and solvent, g = 9.8 m s " 2).

With the application of an external magnetic field Bext, the Iso suspensions exhibit induced birefringence, with significant optical transmission between crossed polarizers when B ext b 10 mT is

applied in the plane of the cell at an angle of 45° to the polarizer, P (Fig. 1b). This induced birefringence is spatially homogeneous, further evidence that the nanoplate dispersion is uniform on the optical scale. The direction with lowest refractive index is found by using a Berek optical compensator to be parallel to Bext. The Iso suspensions also exhibit magnetically induced optical dichroism, with the absorbance largest for optical polarization parallel to the planes of the nanoplates (Supplementary Fig. 2). Thus, the magnetized isotropic phase has field-induced, uniaxial nematic ordering of the plates, with the nematic director n(r), the unit vector along the local, the mean orientation direction of the nanoplate normals, aligned parallel to Bext, along the induced magnetization density M. The induced nematic has negative birefringence An(Bext) = иц — n> <0 and a visible light-dichroic ratio Z(Bext) = OD||/OD><1 (where || (>) indicates orientation parallel (normal) to n).

At higher concentrations (f >0.28), fresh suspensions begin to sediment higher-density domains at the bottom of the capillary that are birefringent even in the absence of applied magnetic field. Over periods of hours to days, a well-defined, horizontal, planar interface between the birefringent region and the upper isotropic phase develops. The birefringent region, if left unperturbed, gradually develops LC-like textures with uniform domains that appear to grow in size and become more ordered with the passage of time (Fig. 1a,c, below the interface, and Supplementary Fig. 4). We identify this behaviour as a two-phase coexistence, with the denser phase sedimented from the less-dense Iso phase by gravity.

Further insight into the ordered phase can be obtained by observing the Iso-NF interface with an external field applied, where we see that an initially sharp interface formed in the absence of field (Fig. 1c) begins to deform and move into the Iso region as Bext is increased. For sufficiently large field (Bext~0.5mT), the ordered phase reorients to have the director along Bext, nematic order induced in the Iso phase becomes comparable to that of the ordered phase and the interface becomes continuous (Fig. 1d), qualitatively as described by the mean-field model of a field-induced isotropic/nematic transition of Fan and Stephen20. Induction of nematic order requires larger Bext as the interface is pushed further up in the cell because the number density of nanoplates in the Iso region decreases with increasing height (Supplementary Fig. 3). If the applied field is removed, the interface reforms and returns to the state of Fig. 1c. These responses to changing field take several minutes to hours to complete because the distribution in the number density of plates must also change to re-establish a uniform chemical potential of the plates in the colloidal system. The field-induced continuity provides evidence that the birefringent phase below the interface is identical to that induced in the Iso phase, a uniaxial, magnetically polar nematic, having the first- and second-rank nanoplate orientational order parameters Q1 and Q2 simultaneously non-zero (ferromagnetic nematic, NF). If Bext is instead reduced to a small but non-zero value, the Iso-NF interface was also found, unexpectedly, to exhibit a field-induced spiking instability, discussed in detail below.

To assess the nature of the magnetic ordering in the vicinity of the Iso-NF transition, we measured and analysed An induced by an external magnetic field Bext, in effect probing the field-induced, uniaxial orientational order in the Iso and NF phases (Supplementary Fig. 5). Analysis including the effects of polydispersity of the response of extremely low-concentration suspensions (f = 0.005) gives an estimate of the mean nanoplate magnetization of mo = 2 x 10 —18 A m2 (Supplementary Note 1), which is very close to the value obtained from aligned and dried BF nanoplates with a mean diameter of 70 nm (ref. 15). As f is increased, the initial (low-field) magnetic susceptibility, w(f), of the Iso phase, derived from the measurements of Q2 = An/Ansat,

where Ansat is the saturated birefringence, shows a substantial growth. Comparison shows that the standard Langevin-Weiss mean-field model21 significantly overestimates the growth of ) and predicts the Iso-NF transition near f = 0.09, a much lower value than observed. However, treating the orientationally diffusing nanoplates of the Iso phase effectively as a system of dipolar hard spheres and applying models that include interparticle correlations11 gives a reasonable qualitative description of the Iso phase susceptibility (Supplementary Note 2).

The LC order of the NF phase was characterized by synchrotron X-ray diffraction. A high-concentration sample (f = 0.28) was kept vertical for several days until it developed a clear Iso-NF interface. As shown in Fig. 2a, with the beam intersecting the upper part of the Iso part of the sample (Location Iso), the diffraction pattern is circularly symmetric, demonstrating the isotropic nature of the sample. When the beam is incident on the two-phase interface (Location Iso-NF), a pair of diffuse arcs appears, with a noticeable isotropic background. The higher-density phase region (Location NF) gives multiple diffraction arcs. If the diffracted signal is integrated along concentric circles around the beam, we obtain the dependence of scattering intensity versus wavevector q shown in Fig. 2b. The diffraction arcs are broad, indicating short-ranged lamellar correlations characteristic of a nematic phase22, with the peaks at almost the same locations in both regions of the sample. Fitting the structure factor S(q) = q2I(q) with the function S(q) = sh(qa)/[ch(qa) — cos(qd)] (a is the average of the magnitude of the orientation difference between neighbouring plates), describing lamellar correlations in an array of plates, gives peaks in the pair correlation function at d = 20.0, 20.9 and 22.1 nm for the NF, Iso-NF and Iso regions, respectively (Fig. 2c), showing that the difference in the number density between the two phases is small. The azimuthal distributions of the first-order diffraction intensities at Location Iso-NF and Location NF are plotted in Fig. 2d. Analysis of the azimuthal distributions of the diffraction intensity gives a second-rank order parameter Q2 = 0.8 (Supplementary Note 3).

Equilibrium magnetic domain structure in a thin capillary. The

response of the higher-density phase to an applied magnetic field was studied in some detail. This phase is a low-viscosity fluid that flows readily if the capillary is tilted and anneals if left undisturbed into domains that are several millimetres across and span the 50-p.m thickness of the capillary (Fig. 3 and Supplementary Fig. 4). Cells were equilibrated over periods of up to several days with external static fields Bequil. If Bequil = 0, obtained by locally cancelling the Earth's magnetic field and other stray fields, a characteristic texture of large, distinctly dichroic and birefringent domains, which divide the cell into almost uniform regions of different orientations separated by sharp domain boundaries as shown in Fig. 3, consistently develops. Measurement of the dichroism and birefringence of these domains and applying the relationship between n, An and Z discussed above in connection with the Iso phase show that the nanoplate director n(r) is oriented as indicated in Fig. 3a. If Bequil is applied parallel to n in the left- and right-hand domains in Fig. 3a, then one of these domains (in this case the left domain) shrinks and the other expands (Fig. 3b). If the sign of Bequil is reversed, then the domain growth/shrinkage is also reversed. The field-induced displacement of the inversion walls is plotted in Supplementary Fig. 6. These walls behave as a system of lines with inherent tension that adopt a geometry minimizing their net length while mediating the changes in director orientation. The field-induced wall displacements are a consequence of

-S ¿S

105 104 103 102 101

0.1 0.3 0.5 0.7 0.9 1.1 1.3 1.5 q (nm-1)

(x105) * 8


Iso-NF • Iso

120 180 270 330 360 в (degrees)

0.4 0.8 q (nm-1)

Figure 2 | Synchrotron X-ray diffraction of BF/BuOH magnetic nanoplatelet suspensions. A suspension with f = 0.28 was filled into a flat capillary of thickness L = 50 mm and width W = 1 mm and sedimented for 1 day, after which it developed an Iso-NF interface. (a) 2D diffraction patterns measured at the three heights in the sample indicated at left by the white boxes representing the beam size and location. Location Iso: immediately below the air-liquid interface; Location Iso-NF: immediately above the Iso-NF interface; Location NF: in the NF phase. (b) Diffraction intensity versus q obtained by circular integration of the 2D images in a. The integrated diffraction from a f = 0.12 sample is shown for comparison. This suspension is uniformly isotropic and the scattering is featureless. (c) Structure factor S(q) = q2I(q) at the three locations probed in the f = 0.28 cell. The position of the first diffuse peak moves from qmax = 0.30 nm "1 in the Iso-NF region to 0.25 nm "1 in the Iso region, corresponding to an ~15% density difference between the two locations, consistent with Supplementary Fig. 3. The solid curves are a fit to the equation shown, modelling the layer correlations in a suspension of thin plates, where a is the average of the magnitude of the orientation difference between neighbouring nanoplates. (d) Diffraction intensity I(d) along the circle passing through the first diffuse peak for the Iso-NF and NF locations, respectively. This angular distribution is fitted by the function I(0) = £n p2ncos2n0 (solid lines), where d is the azimuthal angle as illustrated in the inset.

stresses that tend to change the domain areas according to the sign of the magnetic energy density UM = — M • Bequil. The linear response of the displacement at low fields is a definitive evidence for a magnetization, M, that is finite and independent of applied field strength around Bequil = 0. The field-induced motion of the domains walls is maximal for Bequil in the n direction, showing that M is along n and uniquely determining its sign, as indicated in Fig. 3a. Equilibration in a finite field, such as the Earth's magnetic field, BequilB0.05mT, produces arrays of similarly structured domains that are more complex (Supplementary Fig. 4).

The domain magnetization is also evident from the internal response of the domains to applied fields Btrans changing over periods of several seconds. The big domains in the cell are observed to have a static, grainy texture of subdomains в 10 mm in size in which the in-plane orientation varies around the mean by a few degrees. Applying Btrans parallel to M within a domain produces a continuous enhancement of the extinction of that domain when viewed between crossed polarizers as the magnitude of the applied field is increased (Fig. 3 c). The field applies torques to the nanoplates that tend to reduce the orientational fluctuations of the director field n(r), orienting the plates more parallel on average to Btrans and thereby increasing Q2. Applying Btrans antiparallel to M, on the other hand, generates torques that tend to enhance the orientational disorder of the nanoplates, producing, above a certain threshold BtransB0.005 mT, an instability marked by the appearance of a speckle-like pattern of birefringence generated by an array of subdomains of M(r) that have rotated away from the field direction in random directions and through a distribution of angles. Since these

subdomains are substantially smaller than the cell thickness, they must be a response to static fluctuations and thus to quenched disorder in the local M,n fields. Strong magnetic self-interactions lead to the local expulsion of magnetic charge and the creation of subdomains responding as independent units to Bext. Since these changes take place at values of Bext much smaller than the internal field (BM~40mT, see Supplementary Note 4), the response appears to be a soft deformation mode in the frozen transverse fluctuations of M. This behaviour is quite different from the Freedericksz transitions found in the BF/5CB system23.

This dependence of the domain response on the sign of Btrans, combined with imaging of the dichroism and birefringence-based optical textures, enables detailed mapping of M(r) and n(r) in any domain pattern. M(r) is found to be always parallel to n(r), which is consistent with the nanoplate magnetic moments being along their normals, and with the observed domain wall movements in applied fields. Importantly, the domain magnetization is non-zero in the absence of any applied magnetic field, as an equilibrium property of the phase.

The textural features and response to an applied field of the BF/BuOH ferromagnetic nematic phase may be understood as phenomena resulting from the combined effects of magnetostatic and nematic elastic energies. The interplay of Frank elasticity and magnetostatic self-interactions is expressed through the characteristic length £M = \JK/(moM2) <1 mm (K is the LC Frank elastic constant in the one-constant approximation, see Supplementary Fig. 7 and Supplementary Note 4), which for the BF/BuOH NF phase is small compared with any cell dimension. Generally, the imperative produced by this condition is to reduce magnetostatic energy by reducing the excess internal

JZ? + ® - - > —1. - »

Figure 3 | Equilibrium ferromagnetic nematic flux loop texture in a thin capillary. (a) Orientational domain textures in the NF phase in a cell equilibrated in the absence of applied magnetic field viewed in transmission with optical polarization conditions [P, A] indicated. Uniform domains, with magnetization M(r) orient in a discontinuous loop to minimize magnetic charge at the domain boundaries. The domains extinguish wherever the director field n(r) is parallel or normal to crossed polarizers. Rotating the polarizers confirms that the domains are birefringent. Optical dichroism is revealed with P only, the nanoplate absorbance being larger for P parallel to the nanoplate planes (normal to n). (b) Applying BequN (yellow arrows; from left: BequN = 0.03 mT, 0.04 mT, 0.08 mT and 0.10 mT) induces domain wall displacement, allowing the domain with M parallel to BequN to grow in area. (c) A small transient field (white arrows; Btrans~0.02mT) applied parallel to M makes the orientation within that domain more uniform (darker). Btrans applied antiparallel to M produces a random texture of brighter, reorienting subdomains. (a-c) Scale bar, 200 mm. The cells are 1 mm in width, with the boundaries indicated by thin white lines. (d) An illustration of the twist wall passing along the cell centreline from one cell surface to the other. (e) Structure of the splay/bend/splay domain boundaries interior to the cell, along the line indicated by the olive arrow in a, stabilized by a balance of LC bend elasticity and magnetic charge attraction. (f) Magnified image of the +1 escaped twist wall in n(r) and M(r) imposed by the loop structure of M(r). Scale bar, 50 mm. L = 50 mm, W = 1mm.

and external magnetic fields generated. A distortion such as splay of M(r) that produces bulk magnetic charge of a particular sign [rm(r) = — V • M(r)] must necessarily also produce the same quantity of opposite charge elsewhere in the bulk or at the interfaces with non-magnetized material [ps(r) = M(r) • s(r)]. Therefore, the drive for energy reduction leads either to non-splayed regions of M(r) tangent to boundaries, which produce minimal magnetic charge, or to orientational variation of M(r) such that the positive and negative magnetic charges generated are located as close as possible to each other, a separation that is at least £m because of the LC elasticity (Supplementary Fig. 7). The resulting structure of uniform domains separated by distinct walls is a general feature of the BF/BuOH textures reported here. However, in systems where XM is much larger than the cell dimensions, such as in the case of BF/5CB, the magnetic self-interactions affect the texture only weakly (Supplementary Note 5).

Our observations indicate that the NF region in the flat capillary cells is magnetized with M preferentially parallel to the largest cell face (Fig. 3d). This is expected, considering that, for a slab of infinite area, the magnetic energy of a uniform M field-free to orient in any direction is lowest when M is in the plane of the slab (Supplementary Note 4). The response of the capillary cells to an applied field confirms this confinement effect (Supplementary Fig. 8), showing that, while fields of less than 0.2 mT were needed to induce a complete in-plane alignment of the entire sample, in the case of perpendicular field there is little observable response until Bext> reaches ~10mT, and

extinguishing regions first appear at ~20mT, indicating a domain-mediated transition orienting n normal to the cell face.

The equilibrium in-plane texture in the absence of external applied fields can be understood by approximating the NF region as a two-dimensional (2D) slab the thickness of the cell, in which n(r) is a 2D field, uniform along the cell normal, c (Fig. 3d), and assuming that boundary conditions force n(r) to be parallel to the edges of the cell, which eliminates magnetic charge on the surfaces. M(r), thus constrained to be parallel to the a,b plane in Fig. 3d, adopts an equilibrium in-plane texture that is determined only by the in-plane sample shape and the magnetoelastic behaviour. In the case of rectangular cells with planar boundaries, the equilibrium configuration is a loop of uniform domains of M separated by sharp splay/bend/splay walls, such as that shown in Fig. 3e and Supplementary Fig. 7, forced into the texture of Fig. 3 a by the bends of M(r) at the corners of the rectangle. These walls are observed to have an optical-resolution (micron scale or smaller) linear core and a diffuse width comparable to the cell thickness. Figure 3e shows that the deformation of M(r) through a splay/bend/splay wall generates lines of opposite magnetic charge, the attraction of which favours narrow walls. The Frank elastic energy associated with bend of the director n(r), on the other hand, favours broad walls, to reduce the curvature of the director reorientation at the wall centre24. Using KS b 6kBT/D = 5 x 10 " 13 N from a numerical simulation25, this balance yields an equilibrium wall core structure of width w~£M<1 mm, which is smaller than the optical resolution. This scenario is an example of 'orientational fracture', similar to that

found in — 1 topological defects in the director field structure and textures of thermotropic ferroelectric smectic C LCs with large permanent polarization26, as detailed in Supplementary Note 6.

The loop geometry leads to an overall topological disclination line defect of strength + 1 in the M,n fields of the sample. In the centre of the cell we observe a linear boundary between domains of antiparallel M (Fig. 3a-c), which higher magnification imaging shows to be a + 1 disclination line singularity in the M,n fields (Fig. 3f) that starts on one cell surface at a point where three domains meet (blue square) and ends on the other cell surface at the other such point, as shown in Fig. 3d. This defect can be constructed by starting with a + 1 bend disclination line (blue box) running along c from one cell face to the other (required by the boundary conditions on the cell edges) and then sliding the intersection points on the cell surfaces to the locations where the three domains meet (blue squares). This defect will 'escape' into the third dimension, forming a continuous twist wall in the cell centre (orange square) confined by the mutual attraction of the opposite-signed magnetic space charge.

The equilibrium, field-free domain pattern of M(r) in Fig. 3 resembles that of the closed flux loops found in solid

ferromagnetic single crystals27, known to be an organization of their magnetization field M(r) that reduces magnetic charge rm(r) = — V • M(r) and thereby the magnetic field and energy external to the sample. In crystals, such patterns are achieved with the aid of the crystal anisotropy, which pins M(r) to certain 'easy' symmetry directions. In ferromagnetic LCs, however, there is no underlying lattice to stabilize discontinuous reorientation of M(r). Equilibrium domain patterns are established by the sample shape and by magnetoelastic interactions described in Supplementary Note 4. An interesting special case where the magnetization can vary continuously with no need for domain walls is the ring cell—for example, one shaped like a flat washer or a racetrack. The equilibrium M(r) here is everywhere parallel to the cell boundary, a geometry with continuous bend depending only on radius.

Spikes at the Iso-NF interface. A dramatic manifestation of the fluid ferromagnetic ordering in the BF/BuOH system is the interfacial spiking behaviour shown in Fig. 4 and Supplementary Fig. 4. In fresh cells, the Iso-NF interface is initially flattened by

0 0.01 0.02 0.03 0.04 0.05 0.06

6equil (mT)

e„,. = 0.05 mT

6..., = 0.5 mT

37 min

1 min 2 min

= 0.05 mT

10 min

Figure 4 | Ferromagnetic nematic spikes at the BF/BuOH isotropic/nematic interface. (a) Spikes obtained by equilibrating a f = 0.28 suspension in a rectangular capillary in an external field (yellow arrow, Bequil = 0.05 mT) for 2 days. The optical polarization conditions [P, A] are indicated. The reaction of the spikes to transient fields (white arrows, Btrans~0.02 mT) indicates that n and M in the spikes are in the same direction as Bequil. (b) Spike height versus Bequii. The spike growth is limited by their penetration into the lower-susceptibility Iso phase in the top of the cell. (c) Sketch of ferromagnetic spike at the Iso-NF phase boundary showing magnetization and magnetic charges. (d) Spike dynamics. Application of a large external field (Bext = 0.5 mT) induces saturated M, except near the Iso-air interface. On reducing Bext to 0.05 mT, the NF-Iso interface reforms initially with an undulation instability, which then evolves into spikes. After application of an external field Bext = 0.05 mT to an initially smooth interface formed at Bext = 0, on interfaces instability appears and grows into spikes. The cells are 1 mm in width with boundaries indicated by the thin white lines. Scale bar, 200 mm.

gravitational sedimentation of domains of NF (Fig. 1c). However, applying an external magnetic field of Bequil >0.005 mT normal to this interface causes a quasi-periodic array of birefringent spikes to grow, which equilibrate to a field-dependent height above the interface over the course of a day (Fig. 4a,b). Exposure of the spikes to transient fields of few-second duration shows that they are nematic and ferromagnetic, with magnetization parallel to the spike axis along Bequil, with a field response very similar to that of the bulk NF domain.

The spike geometry appears to be a soliton-like response of the interface that enables the magnetization of a ferromagnetic fluid to align with Bext. The extended shape of the spikes, with M(r) oriented along the spike axis, minimizes the generation of magnetic charge, with the remnant opposite magnetic charges on opposite ends of the spikes, indicated in Fig. 4c, decreasing as the spikes become longer and narrower. Each spike is effectively a small bar magnet, making the spikes mutually repulsive, as is the case for similarly aligned bar magnets placed next to one another. The spikes are roughly conical in shape, becoming narrower in both the a and c directions with increasing distance from the interface, the latter case evident from the birefringent interference extinctions at places where the phase shift is an odd multiple of p.

Spikes tend to grow in an external field if the magnetic energy gain is larger than the cost of increasing the interface area, for example, for a cylindrical spike of radius RB25 mm when Bext>2aIN/RM~0.02mT. As discussed in connection with the Iso phase, application of a sufficiently large Bext normal to the interface of an equilibrated sample induces nematic order everywhere, expanding the strong ordering of the bulk nematic well into the Iso region for Bext>0.5mT (Fig. 1d). On the removal of Bext from such a sample previously equilibrated with spikes, both the interface and the spikes evolve quickly over a period of a few minutes, as shown in Fig. 4d, to an instability of interface position. This experiment shows that the spike length at a given field is limited by the decrease in the susceptibility to magnetic ordering with distance from the Iso-NF interface.

We propose that spike formation can be understood qualitatively as a novel type of nonlinear capillary instability of the sort produced in fluid interfaces when they become electrically or magnetically charged by the application of an external field3,5,28. As discussed above, in a ferromagnetic fluid, the permanent magnetization M is parallel to the interface when Bext = 0, but, as M rotates in response to an applied field normal to the interface, it develops a normal component M> = Bext/mo. This, in turn, produces magnetic charge Bext/mo, and an effective magnetization energy density for driving the instability. Near the threshold, the instabilities are periodic with wavelength 1G = 2p • lG, where lG = \Ja/(Apg) is the gravitational capillary length (a = surface tension, Ap = mass density difference between Iso and NF phases). The associated increase in capillary energy density3 is UG = \J a/(Apg). Taking sIN = 0.1kBT/D2~10"7Jm"2 as an estimate of the Iso-NF interfacial tension29,30, we find 1G ~ 200 mm, comparable to the observed spike spacing, and UG~0.002Jm_ 3, energy that is provided by the applied field.


A variety of observations discussed above indicate that the Iso and Nf phases are equilibrium dispersions of individual nanoplates. First, the magnitude of the magnetization of the orientable unit producing the magnetic birefringence in the Iso phase is near that of a single plate, which must therefore act as thermalized individuals. Second, the Iso-NF phase coexistence appears at f~0.28, equivalent to a scaled density of v<D3> = 3.3, a volume fraction value where nematic ordering

due to the steric and charge interactions of single nanoplates is expected18,19,31. Third, X-ray diffraction from oriented monodomains of the NF phase reveals several diffuse quasi-Bragg reflections from short-ranged lamellar ordering with a layer spacing of d = 20 nm (Fig. 2). Assuming that the nanoplates form nearly continuous sheets, that is, the volume fraction f = 0.28 is determined by the spacing between the neighbouring plates, it gives a plate thickness t = fd = 5.6 nm. This is comparable to the 7-nm thickness of a single nanoplate, providing further evidence that the BF/BuOH system is a suspension of single plates dispersed in equilibrium and distinguishing the ferromagnetism observed here from that of the non-ergodic aggregation of magnetic particles8,9. Finally, observations above suggest that magnetic dipolar and steric quadrupolar interactions, tending to stabilize Q1 and Q2, respectively, contribute to the ordering transition.

Several theoretical treatments of ordering transitions in systems with both dipolar and quadrupolar interactions have explored the interplay of the ordering of Q1 and Q2, in particular the Sivardiere and Blume (SB) Ising model in the molecular field approximation32, and the Baus and Colot (BC) density functional theory of a system of dipolar hard ellipsoids33. Despite the very different nature of these two systems, the SB and BC phase diagrams, shown in Supplementary Fig. 9, exhibit some key qualitative common features relevant to the BF/BuOH system. Both the SB and BC show that the two types of ordering are surprisingly weakly coupled, with the isotropic/nematic quadrupolar transition independent of the dipole interaction until the latter becomes strong enough for dipole ordering to occur at a smaller density than the quadrupolar ordering. On the other hand, the dependence of the (isotropic or nematic)/NF transition on the dipole strength is only weakly affected by whether the quadrupolar ordering has taken place or not. Our independent assessment above of the proximity of the two transitions at fB0.28 leads us to conclude that on these phase diagrams the BF/BuOH system is in the vicinity of the tricritical point where, with an increasing dipole strength, the nonferromagnetic nematic phase has just disappeared. This is a place where the remnant Iso-NF transition will be first order, and where there will be significant short-ranged quadrupolar ordering (nonpolar alignment of the nanoplates) and thus an enhanced susceptibility for inducing Q2.

Those familiar with smectic LCs may recognize the similarity of the phase diagram in Supplementary Fig. 9 to McMillan mean-field isotropic/nematic/smectic A model, with ferromagnetic ordering replacing the appearance of smectic A layering34. McMillan's distinct mean-field approach is thus yet another route to the generic dipolar/quadrupolar phase diagram. The Landau model of Palffy-Muhoray, Lee and Petschek also produces a similar Iso, nematic, NF phase diagram35.


Preparation of the BF suspensions. BF nanoplatelets were suspended in n-butanol to form a stable colloidal suspension following a previously established method36,37. Briefly, BF nanoplatelets with the nominal composition BaFe115Sc0.5O19 were synthesized hydrothermally at 240 °C, with dodecylbenzenesulphonic acid added as a surfactant. The nanoplatelets were subsequently washed with water and nitric acid. To increase the adsorption of the dodecylbenzenesulphonic acid on the nanoplatelets, they were suspended in an HNO3 solution (pH = 1.5) and heated at 100 °C for 2.5 h. After this, the nanoplatelets were washed again with water and acetone, dried in air at 60 °C and dispersed in n-butanol using probe sonication. The final volume fraction of platelets f = 0.003 was obtained by concentrating the initial suspension (f = 0.001) in a rotary evaporator. The functionalized nanoplatelets have an average magnetic core thickness of 5 nm, with 1 nm of surfactant on each side, and a mean diameter of D = 48 nm with a s.d. of 21 nm, which gives an aspect ratio of B7. Suspensions prepared at the volume fraction obtained in the synthesis of the platelets, f = 0.003, formed an isotropic liquid. Higher concentrations were prepared by compression, using ultracentrifugation at 160,000g for 40 min. The

clear supernatant was then removed and the remaining concentrated suspension, at f b0.1, homogenized. Repeated addition of low-concentration suspension followed by recompression eventually yields suspensions with desired higher concentrations. The suspensions, loaded into the 50 mm thick x 1 mm wide x 25 mm long sample cavity of rectangular capillaries, flame-sealed and held vertically, are stable for months under ambient conditions, showing no sedimentation or aggregates visible in the optical microscope.

Synchrotron X-ray experiment. The orientational order of the nanoplates in a phase-separated suspension was studied using synchrotron X-ray diffraction performed on Beamline 7.3.3 of the Advanced Light Source, Berkeley, CA, USA. The suspensions were contained in rectangular (0.05 x 1.0 mm2) borosilicate glass capillaries. The sample was probed at multiple locations in the capillary with a 10-keV synchrotron X-ray beam with a Mo/B4C double-multilayer mono-chromator and a beam at the sample ~700 mm wide by 300 mm high.


1. Papell, S. S. Low viscosity magnetic fluid obtained by the colloidal suspension of magnetic particles. US Patent 3215572 (1965).

2. Rosensweig, R. E. Ferrohydrodynamics (Cambridge University Press, 1985). 3

Cowley, M. D. & Rosensweig, R. E. The interfacial stability of a ferromagnetic fluid. J. Fluid Mech. 30, 671-688 (1967).

Bacri, J. C., Cebers, A. & Perzynski, R. Behavior of a magnetic fluid microdrop in a rotating magnetic-field. Phys. Rev. Lett. 72, 2705-2708 (1994). Odenbach, S. Colloidal Magnetic Fluids: Basics, Development, and Application of Ferrofluids, Lecture Notes in Physics (Springer, 2009). Berkovski, B. & Bashtovoy, V. Magnetic Fluids and Applications Handbook (Begell House Inc., 1996).

Brochard, F. & de Gennes, P. G. Theory of magnetic suspensions in liquid crystals. J. Phys. 31, 691-708 (1970).

Rault, J., Cladis, P. E. & Burger, J. P. Ferronematics. Phys. Lett. A 32, 199-200 (1970).

Ivanov, A. O. & Kantorovich, S. S. Chain aggregate structure and magnetic birefringence in polydisperse ferrofluids. Phys. Rev. E 70, 021401 (2004).

10. Wang, Z. W., Holm, C. & Muller, H. W. Molecular dynamics study on the equilibrium magnetization properties and structure of ferrofluids. Phys. Rev. E 66, 021405 (2002).

11. Huke, B. & Lucke, M. Magnetic properties of colloidal suspensions of interacting magnetic particles. Rep. Prog. Phys. 67, 1731-1768 (2004).

12. Onsager, L. The effects of shape on the interaction of colloidal particles. Ann. N Y Acad. Sci. 51, 627-659 (1949).

13. van der Kooij, F. M., Kassapidou, K. & Kekkerkerker, H. N. W. Liquid crystal phase transitions in suspensions of polydisperse plate-like particles. Nature 406, 868-871 (2000).

14. Gabriel, J. C. P. & Davidson, P. New trends in colloidal liquid crystals based on mineral moieties. Adv. Mater. 12, 9-20 (2000).

15. Mertelj, A., Lisjak, D., Drofenik, M. & Copic, M. Ferromagnetism in suspensions of magnetic platelets in liquid crystal. Nature 504, 237-241 (2013).

16. Clark, N. A. Soft-matter physics: ferromagnetic ferrofluids. Nature 504, 229-230 (2013).

17. Frenkel, D. & Eppenga, R. Monte Carlo study of the isotropic-nematic transition in a fluid of thin hard disks. Phys. Rev. Lett. 49, 1089-1092 (1982).

18. Bates, M. A. & Frenkel, D. Nematic-isotropic transition in polydisperse systems of infinitely thin hard platelets. J. Chem. Phys. 110, 6553-6559 (1999).

19. Bates, M. A. Influence of particle shape on the nematic-isotropic transition of colloidal platelet systems. J. Chem. Phys. 111, 1732-1736 (1999).

20. Fan, C. P. & Stephen, M. J. Isotropic-nematic phase transition in liquid crystals. Phys. Rev. Lett. 25, 500-503 (1970).

21. Reitz, J. R., Milford, F. J. & Christy, R. W. Foundations of Electromagnetic Theory (Addison-Wesley, 1979).

22. Michot, L. J. et al. Liquid-crystalline aqueous clay suspensions. Proc. Natl Acad. Sci. USA 103, 16101-16104 (2006).

23. Mertelj, A., Osterman, N., Lisjak, D. & Copic, M. Magneto-optic and converse magnetoelectric effects in a ferromagnetic liquid crystal. Soft Matter 10, 9065-9072 (2014).

24. Zhuang, Z., Maclennan, J. E. & Clark, N. A. Device applications of ferroelectric liquid-crystals - importance of polarization charge interactions. Proc. SPIE 1080, 110-114 (1989).

25. O'Brien, P. A., Allen, M. P., Cheung, D. L., Dennison, M. & Masters, A. Elastic constants of hard thick platelets by Monte Carlo simulation and virial expansion. Soft Matter 7, 153-162 (2011).

26. Link, D. R., Chattham, N., Maclennan, J. E. & Clark, N. A. Effect of high spontaneous polarization on defect structures and orientational dynamics of tilted chiral smectic freely suspended films. Phys. Rev. E 71, 021704 (2005).

27. Kittell, C. Physical theory of ferromagnetic domains. Rev. Mod. Phys. 21, 541-583 (1949).

28. Zelazo, R. E. & Melcher, J. R. Dynamics and stability of ferrofluids - surface interactions. J. Fluid Mech. 39, 1-24 (1969).

29. Reich, H., Dijkstra, M., van Roij, R. & Schmidt, M. Entropic wetting and the free isotropic-nematic interface of hard colloidal platelets. J. Phys. Chem. B 111, 7825-7835 (2007).

30. van der Beek, D. et al. Isotropic-nematic interface and wetting in suspensions of colloidal platelets. Phys. Rev. Lett. 97, 087801 (2006).

31. Rowan, D. G. & Hansen, J. P. Salt-induced ordering in lamellar colloids. Langmuir 18, 2063-2068 (2002).

32. Sivardie, J. & Blume, M. Dipolar and quadrupolar ordering in S = 3/2 Ising systems. Phys. Rev. B 5, 1126-1134 (1972).

33. Baus, M. & Colot, J. L. Ferroelectric nematic liquid-crystal phases of dipolar hard ellipsoids. Phys. Rev. A 40, 5444-5446 (1972).

34. McMillan, W. L. Simple molecular model for the smectic A phase of liquid crystals. Phys. Rev. A 4, 1238-1246 (1971).

35. Palffy-Muhoray, P., Lee, M. A. & Petschek, R. G. Ferroelectric nematic liquid crystals: realizability and molecular constraints. Phys. Rev. Lett. 60, 2303-2306 (1988).

36. Lisjak, D. & Ovtar, S. Directed assembly of BaFe12O19 particles

and the formation of magnetically oriented films. Langmuir 27, 14014-14024 (2011).

37. Lisjak, D. & Drofenik, M. Chemical substitution—an alternative strategy for controlling the particle size of Barium Ferrite. Cryst. Growth Des. 12, 5174-5179 (2012).


This work was supported by the Soft Materials Research Center under NSF MRSEC Grants DMR-0820579 and DMR-1420736, by the Institute for Complex Adaptive Matter Postdoctoral Fellowship Award OCG5711B and by Slovenian Research Agency Grants P1-0192 and P2-0089-4. The synchrotron X-ray experiments were supported by the Director of the Office of Science and Office of Basic Energy Sciences of the US Department of Energy under Contract No. DE-AC02-05CH11231.

Author contributions

M.S. and N.A.C. designed the research. M.S., A.K., Y.S., G.P.S., M.R.T., C.Z. and N.A.C. performed the experiments. A.M., D.L. and M.C. supplied the materials. M.S., R.G.P., M.A.G., J.E.M. and N.A.C. analysed and interpreted the data. M.S., J.E.M. and N.A.C. wrote the manuscript. All the authors read and commented on the manuscript.

Additional information

Supplementary Information accompanies this paper at naturecommunications

Competing financial interests: The authors declare no competing financial interests.

Reprints and permission information is available online at reprintsandpermissions/

How to cite this article: Shuai, M. et al. Spontaneous liquid crystal and ferromagnetic ordering of colloidal magnetic nanoplates. Nat. Commun. 7:10394 doi: 10.1038/ncomms10394 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit