Jordan Austin-Frank Wilson,
Allen Sun
and
Zhiwei Li
*
Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, USA. E-mail: zlinano1@umd.edu
First published on 15th August 2025
Magnetic north and south poles attract, leading to the well-known end-on attachment between two magnets. In this work, we leverage distinct size-dependent dipole–dipole interactions between nanorod magnets-magnetized nanorods-to assemble hierarchical photonic stripes and provide an effective mechanism to prepare complex photonic materials with precisely tunable structural colours. Under a magnetic field, Fe3O4@SiO2 nanorods are magnetized with an induced parallel magnetic dipole, favour a size-dependent equilibrium bonding state, and then assemble along a critical angle into an offset packing instead of the end-on attachments. Our experimental results and computational model indicate that magnetic nanorods exhibit a distinct bonding mode compared to magnetic nanospheres, underpinning the important role of magnetic anisotropy in determining the nanoscale magnetic interactions, structural symmetry, and photonic properties. This scale-dependent interaction produces tetragonal photonic crystals, which further assemble into hierarchical photonic stripes driven by magnetic repulsion between crystals at a high concentration under the assistance of a template. This work provides an assembly approach to fabricating photonic materials and devices that have orders across different length scales based on unique interactions between nanostructured materials.
Among various physical interactions, capillary force is one of the most broadly used physical interactions to assemble nanoparticles into photonic crystals.15–17 During this physical process, the nanoparticle building blocks are first dispersed into a good solvent and then assembled on pre-designed substrates through solvent evaporation.18–20 The strong capillary force assembles the nanoparticles into close-packed photonic crystals.8,12,21 Because of the strong connection between the packed nanoparticles, it is practically difficult to dynamically tune the crystal periodicity and structural colours on the substrates. Therefore, exploring new driving forces that allow one to actively tune the connection between the building blocks is important to extend the applications of photonic crystals.13 Instead of the evaporation-assembly approach to solid-state materials, magnetic assembly was developed to assemble magnetic nanoparticles into ordered structures in a solution phase.22–27 The dipole–dipole coupling between neighbouring magnetic nanoparticles drives the assembly of two-dimensional or three-dimensional photonic crystals with vivid structural colours in a colloidal solution.22,23,28–32 Because the strength of the induced dipoles of individual nanoparticles and the dipolar coupling between them is dependent on the strength of the external magnetic field, it is possible to tune the particle interactions and therefore the structural periodicity by simply changing the strength or direction of the external fields.32–36 Therefore, magnetic assembly has been developed as a dynamic approach to responsive photonic crystals that show reversible colour switching in response to external stimuli.1,36,37 For example, different types of nanoparticles, including nanospheres, nanocubes, and nanorods, have been assembled into responsive photonic crystals in solutions, which exhibit precisely controlled structural colours in response to either magnetic field strength or direction changes.26,27,30,38,39 However, further processing these photonic materials into devices needs additional mechanisms to pattern the crystal materials into macroscale periodic structures that have order and complexity in different length scales. Here, we report direct magnetic assembly of Fe3O4@SiO2 nanorods into hierarchical photonic crystals with orders across nano to macroscale and develop a reliable approach to photonic patterns with dynamic colour switching in response to a magnetic field. We studied the shape anisotropy of magnetic nanoparticles by using nanospheres and nanorods as model systems and discussed how the nanorod shape leads to a size-dependent bonding angle between two interacting nanorods as well as the formation of tetragonal photonic crystals in a magnetic field. This work elucidates new insights into the interaction dynamics between nanoparticles and provides a reliable approach to complex and responsive photonic patterns that can be used in device fabrication, sensing, anticounterfeiting, and information encryption.
The structure of the tetragonal photonic crystals is shown in the TEM image in Fig. 2c, which demonstrates high-quality periodicity and well-defined positional and orientational order of nanorods inside the single crystal. This TEM image is a projection of the single tetragonal crystal along its 〈100〉 crystal direction, leading to a centred rectangular phase. One unit cell is shown in Fig. 2c to illustrate the crystal orientation in the TEM image. The position of the nanorods is highlighted by red dots in Fig. S1a to elucidate the symmetry and positional order of the tetragonal photonic crystals. To analyse the unit cell, two parameters are introduced in Fig. S1, including c for periodicity along vertical direction and a for periodicity along other two directions. Based on the TEM image in Fig. S1b, we measured the pixel intensity of this image along the a- and c-axis, which is highlighted by green and red dashed lines in Fig. S1b. As demonstrated in Fig. S1c and d, unit cell parameter c is 476 ± 14 nm while a is 255 ± 23 nm, demonstrating high positional order and high-quality periodicity of the tetragonal photonic crystals. Under an optical microscope, a high-concentration solution of the magnetic nanorods does not show obvious structural colours (Fig. S2a). Under an external magnetic field, however, the magnetic nanorods assemble into tetragonal microcrystals and exhibit bright structural colours. In a horizontal magnetic field (0°), the crystals show green colours. When the field increases from the horizontal to vertical direction, the structural colours change from green to red under an optical microscope (Fig. 2d–f). Similarly, the photonic crystals have brilliant structural colours in a capillary and show colour changes from blue to red by changing the field directions (Fig. S2b) or viewing angles (Fig. S2c). This redshift of structural colours is due to the periodicity increases as the orientation changes of the microcrystals. Due to magnetic shape anisotropy, the nanorods and the anisotropic photonic crystals prefer to parallelly align their long axes to the external magnetic field directions, making it possible to dynamically tune the crystal orientation in solution phases. As a result, under a vertical incident light and a horizontal magnetic field (Fig. S1d), the diffraction occurs along the short axes of the nanorods, and the small periodicity produces green structural colours. When the magnetic field changes from horizontal to vertical directions, the periodicity will increase driven by the preferred alignment of crystals to the external magnetic field, leading to the redshift of the structural colours (Fig. S1d). The viewing angle-dependent structural colours can also be explained by the periodicity increases in Fig. S2e. When the viewing angle is zero, the diffraction occurs along the short axes of the nanorods, leading to green colour. Increasing the viewing angle to 60° leads to an increase in periodicity as well as the redshift of the observed colours.
To elucidate the formation mechanism of the tetragonal photonic crystals, we developed a three-dimensional model and computed the magnetic interactions between two nanorods using finite element analysis. Consider the shape anisotropy of the nanorods and image the magnetic interactions between the induced magnetic dipoles in each nanorod. It is difficult to calculate their magnetic interactions using the super magnetic dipole approximation, in which one assumes a magnetic dipole in the whole volume and neglects the shape effect of the nanorod. We consider the interactions between nanospheres and between nanorods in two models to elucidate how the shapes determine the interaction directionality and magnitudes. To do so, we use the finite element method and divide the magnetic nanospheres and nanorods into many, tiny volumes, with each volume having an induced magnetic dipole. The magnetic interactions are first calculated between any volume pairs of the two nanoparticles, and the force vectors are generated by summing over the whole nanoparticle volume to compute the overall magnetic interactions between the two nanoparticles. In Fig. 3a, we start with a two-nanosphere system and compute the magnetic force by changing one particle position from 0° to 90° under the presence of a horizontal magnetic field. The dependence of the magnetic interactions on the angles and the nanosphere separation is illustrated in Fig. 3b. Because the magnetic interaction between two dipoles is inversely proportional to the fourth power of their separation,22,23 it is reasonable that the magnetic force magnitude decreases consistently as interparticle separation increases. The change in force direction is more important but difficult to analyse based on the force vector plots in Fig. 3b. We therefore decompose the magnetic forces into tangent and normal components in Fig. 3c and d, respectively. In the horizontal magnetic field, the tangent magnetic force first increases and then decreases as the angle changes from 0° to 90°. While there is a maximum magnitude, its direction is always pointing to the 90° position within the simulation domain, indicating that there is a consistent force to assemble the two nanospheres into linear chains along the 90° into an end-on configuration (Fig. 3c). We observed a simultaneous change in the normal component, whose interaction direction is highly dependent on the position of the nanospheres as defined by the angle. At 90°, there are strong attractions between the two nanospheres as the normal component is pointing to the primary nanosphere, leading to assembly of two nanospheres into the chains. At 0°, however, the normal component is pointing away, demonstrating a strong repulsion between the two nanospheres. This force analysis based on the finite element method predicts that the magnetic interaction between two identical nanospheres are directional, leading to repulsion and attraction domains depending on their separation, and that nanospheres assemble into linear chains along the direction of the external magnetic field. To experimentally verify this theoretical model, we used magnetic nanospheres as building blocks (Fig. S3) and assembled them by applying a magnetic field. As shown in the TEM image in Fig. 3e, we observed the formation of long chains made of the magnetic nanospheres. Therefore, this theoretical consideration is consistent with the long chains observed in TEM image as well as previous results on magnetic assembly of nanospheres into 1D chains.23,43
To study the magnetic interactions between magnetic nanorods, we used the same finite element analysis on two identical nanorods in a horizontal uniform magnetic field (Fig. 4a). The interparticle separation remains constant when the position of the secondary nanorod changes from 0° to 90°. The two magnetic nanorods are both horizontally oriented in the horizontal magnetic field due to the Zeeman coupling between the induced magnetic dipoles and the horizontal field.26,44
In an analogy to the case of the magnetic nanospheres, we divide the nanorod into tiny volumes and assume one induced magnetic dipole in each volume. Then, we compute the magnetic interactions between each volume pair in the two nanorods and sum over the volume pair interactions to calculate the overall magnetic forces. Fig. 4b shows the dependence of the magnetic interactions between the two nanorods on position and separation. We further plot the normal and tangent components of the magnetic interactions by mathematically decomposing the magnetic interactions. Based on the directionality of the normal force, we observed the presence of attraction and repulsion domains (Fig. 4c). Within the attraction domain, we noticed gradual decrease and then increase of the tangent component, demonstrating the presence of a critical point when the tangent component becomes zero (Fig. 4d). Further analysis on this critical point implies that there is no relative shift between the two nanorods and therefore defines a thermodynamic equilibrium state. At 90°, the tangent component also turns zero and the normal component demonstrates a strong attraction between the nanorods. However, a tiny thermal fluctuation would disassemble the nanorods from this position at 90° because, once the secondary nanorod is off this point driven by the fluctuation, the tangent force increases and is directing the nanorod away from the 90° position. One sequence of this event is the movement of the secondary nanorod to the position where the tangent force is zero within the attraction domain, which is driven by the normal force and leads to the equilibrium of the two nanorods at the critical angles. This theoretical analysis is important to identify the thermodynamics of the nanorod assembly, demonstrating a quasi-stable state at 90° and the thermodynamic equilibrium state at the critical angle. More importantly, this model predicts that the magnetic nanorods prefer to assemble along a size-dependent critical angle. This new understanding explains the formation of the tetragonal photonic crystals (Fig. 4e). Due to the shape anisotropy of the nanorods, the resulting equilibrium status at the critical angle is different from our understanding of the bulky magnet interactions, which favour end-on assembly due to the north–south pole attraction. The packing of magnetic nanorods along the critical angles leads to the formation of the tetragonal photonic crystals. The predicted equilibrium state of the nanorods can be clearly verified by the positional orders in the TEM (Fig. 4f) and SEM (Fig. 4g) images. Considering the unit cell of the body-centered tetragonal crystal, we can experimentally determine the critical angle based on the unit cell parameters a and c (Fig. S4a). In the {100} plane, we measured the a and c to be 255 ± 23 nm and 476 ± 14 nm, respectively (Fig. S1). On the other hand, the critical angle is the angle (θ) between the line (blue arrow in the bottom left in Fig. S4a) connecting the two nanorods and the a-axis (green arrow in the bottom right in Fig. S4a) in the {110} plane. Overlapping multiple {100} and {110} planes produces different configurations, which can be used to identify the crystal orientation in TEM images. As a result, overlapping {110} and {100} planes lead to the formation of rectangular (Fig. S4b) and body-centered rectangular (Fig. S4c) phases, respectively. Based on this prediction, we could conclude that the projection image in Fig. S4d represents the multilayers of {100} planes, which is consistent with our understanding of Fig. 2c. Based on the geometry relation between the rod packing in the {100} and {110} planes, we find that there exists a geometry relation and that θ is calculated to be 52.85° by using the unit cell parameters (a and c) measured from the TEM image in Fig. S4d and plotted in Fig. S1c and d. This experimental measurement can be verified by our theoretical model that is illustrated in Fig. 4. To do so, we considered different silica thickness from 5 nm to 50 nm and computed the tangent and normal magnetic forces between the magnetic nanorods. The results are plotted in Fig. S5a and b for the tangent and normal forces, respectively. Specifically, the zero points in the normal forces separate the repulsion domain from the attraction domain. When silica thickness increases from 5 nm to 50 nm, the angle where the normal magnetic force is zero decreases (Fig. S5b). By comparing the two plots, we can conclude that the critical angles are located within the attraction regime and that it is possible to predict the critical angles for nanorods of different silica thickness. In our experiments, the silica thickness was measured to be 50 nm. The critical angle when the tangent force is zero for 50-nm thick silica shells is shown to be 53.3° in Fig. S5a. This theoretically predicted angle by our analytical model is consistent with the experimentally calculated critical angle of 52.85° as derived from the TEM images and the periodicity plots in Fig. S1.
To observe the long-range order of the tetragonal photonic crystals, we could selectively etch away the Fe3O4 by HCl in the crystals, producing tetragonal crystals of the same structures made of hollow silica shells. A low-magnification TEM image of the SiO2 nanoshells shows the long-range order of the tetragonal photonic crystals (Fig. S6). This observation combined with above structural analysis demonstrates the long-range order and high-quality periodicity of the formed photonic crystals. To further illustrate the positional order of the nanorods in the tetragonal crystals, two TEM images in Fig. S7a and b represent the projection of the (110) and (100) planes, respectively. Their relative crystal orientation can be explained by the crystal unit cell in the insets, which are consistent with the analysis presented in Fig. S4. Overall, our theoretical and experimental studies on magnetic interactions and assembly of nanorods provide new understanding of magnetic interactions between nanorods and new insights into the anisotropic interaction dynamics between anisotropic nanostructures, which further inspires the development of emerging photonic and electronic materials.
The presence of the size-dependent critical angles leads to the formation of offset packing of the nanorods into tetragonal crystals. Such periodic superstructures explain the observation of structural colours in an optical microscope (Fig. 1d–f). When prepared in a high concentration with the assistance of a template (Fig. S8), we observed well-defined hierarchical photonic stripes under the optical microscope (Fig. 5a). The photonic strips were aligned parallel to the external magnetic field and exhibit different colours if the field direction was changed. Such a simultaneous red shift of the structural colours is consistent with the increases in periodicity during the field direction changes. We then measured the reflection spectra of the photonic stripes during this process, and the strong diffraction peaks occur in the visible range due to the constructive interference of the light with the tetragonal crystals (Fig. 5b). As the field angle increases, the tetragonal photonic crystals rotate according to the external field. Under normal light incidence, the periodicity will increase, leading to the redshift of the structural colours. To verify the hierarchical photonic stripe structures and formation, the assembled photonic crystals were fixed by silica coating under the presence of the external magnetic field. While the nanorods were assembled into the hierarchical photonic stripes, the ongoing SiO2 coating simultaneously occurred on the nanorod surface, which stabilized the hierarchical photonic stripes by forming a conformal coating. The SEM image in Fig. 5c indicates the strip-like morphologies of the assembled photonic crystals in the external magnetic field. Crystals have a very high aspect ratio in the magnetic field. Under high concentration and strong magnetic field, the nanorods first assembled in tetragonal photonic crystals with short fiber-like morphologies. Due to the parallel alignment of the fibers in the magnetic field, they experienced strong repulsion between their neighbouring fibers. One consequence of this event is the formation of longer photonic stripes through fusion of multiple fibers along the longitudinal direction, which is further regulated by the presence of the template. The strong repulsion along the transverse direction under the presence of the channel template produces the long-range periodic superstructures made of individual photonic stripe. The tetragonal packing of nanorods in the photonic stripes can be further verified in the TEM image in Fig. 5d, which confirms the formation of tetragonal photonic crystals and the long-range hierarchical photonic stripes in high concentration nanorod dispersion and a strong magnetic field.
Polydimethylsiloxane (PDMS) templates were used to assemble the photonic stripes under the presence of a magnetic field. The structure of the templates is shown in Fig. S8. The PDMS films were prepared by mixing silicone elastomer base and curing agent at room temperature with a ratio of 10:
1 (base to curing agent). The trapped bubbles during mixing were removed under vacuum for ∼1 hour. The prepolymer solution was then transferred to the surface of a commercial silicon wafer with stripe patterns in a Petri dish. Afterwards, the Petri dish was placed in an over at 80 °C for 2 hours to cure the PDMS template films. Through this soft lithography process, the PDMS templates can be prepared for assembly. To assemble the magnetic nanorods into photonic stripes, about 20 μL of the magnetic nanorods was added to the PDMS templates sitting on a glass subtract. A cover glass was placed on top of the liquid droplets. The sample was observed under an optical microscope while a magnetic field was applied to assemble the nanorods. The formed photonic stripes could be observed under the optical microscope.
This journal is © The Royal Society of Chemistry 2025 |