Open access peer-reviewed chapter - ONLINE FIRST

Self-Assembly of Nanoparticles Decorated by Liquid Crystalline Groups: Computer Simulations

By Jaroslav Ilnytskyi

Submitted: August 13th 2019Reviewed: September 11th 2019Published: May 22nd 2020

DOI: 10.5772/intechopen.89682

Downloaded: 10


We present the results of the computer simulations for the self-assembly of decorated nanoparticles. The models are rather generic and comprise a central core and a shell of ligands containing terminal liquid crystalline group, including the case of the azobenzene chromophores. The simulations are performed using the coarse-grained molecular dynamics with the effective soft-core interparticle interaction potentials obtained from the atomistic simulations. The discussion is centred around the set of the self-assembled morphologies in a melt of 100–200 of such decorated nanoparticles obtained upon the change of the temperature, surface density of ligands, the type of the terminal group attachment, as well as the prediction of the possibility of photo-assisted self-assembly of the nanoparticles decorated by the azobenzene chromophores.


  • self-assembly
  • nanoparticles
  • molecular dynamics
  • azobenzene

1. Introduction

Nanoparticles gained extended research and technological focus due to their unique optical, electronic, magnetic and chemical properties [1]. Applications include: medical diagnostics, drug delivery, cancer therapy, nanoelectronics and information storage, sensors, (photo)catalysis, surface coatings [2]. Self-assembly of nanoparticles are in a core of many advanced materials developments coining the term ‘nanoarchitectonics’ [3]. Self-assembly process becomes more controllable by decoration of nanoparticles with functional ligands. Good example is decorated nanoparticles (DNP) coated by the liquid crystalline (LC) ligands [4, 5, 6, 7, 8, 9] that are considered in this chapter.

The structure of the self-assembled morphology depends on: (i) the details of molecular architecture, and (ii) external conditions. The group (i) includes the type of the core nanoparticle (metal/non-metal, magnetic/non-magnetic, etc.); grafting density, length, flexibility and chemical details of ligands; number and type of functionalisation groups, etc. The group (ii) includes the temperature, pressure/density, the presence of surfaces and external fields. Prediction of symmetry/structure/properties of the self-assembled morphology, especially by means of theoretical methods, is quite problematic.

Computer simulations are of great aid here, as these incorporate the relevant details of molecular architecture and tackle statistical behaviour of many-particle systems under specified external conditions. This approach faces a difficulty in terms of a wide span of the time- and length-scales typical for the supramolecular self-assembly. In particular, if more chemistry-specific details are taken into account—then relatively small system sizes and short simulation times can be covered. A good compromise can be achieved by using elements of a multiscale approach [10, 11], which builds a coarse-grained model based on the simulation data of a more chemically-detailed model. A coarse-grained model is of rather generic type capturing essential physical details of the atomistic system and allowing to reach required time- and length-scales of a self-assembly [12, 13]. Due to inevitable loss of specific chemical details, the comparison with particular chemical realisations is performed on a high level only—via the structure of observed morphologies and via the temperature- or density-driven phase transitions between them.

We cover the details of a coarse-grained modelling and self-assembly of DNPs containing LC groups (including the case of azobenzene chromophores). Computer simulations are performed via the molecular dynamics simulation combined with stochastic photoisomerisation events (in the case of azobenzene chromophores). As the result, this type of modelling enables to consider the role of grafting density and type of LC group attachment, as well as the role of the temperature and external fields (including illumination) on the process of self-assembly. Section 2 contains modelling and simulation details, in Section 3 we consider temperature related effects of the self-assembly, Section 4 covers the role of the details of molecular architecture, in Section 5 we cover the photo-aided self-assembly of DNPs containing azobenzene chromophores.

2. Modelling and simulation details

To study self-assembly of DNPs we use coarse-grained modelling, where relevant groups of atoms are replaced by single beads that interact via soft-core potentials. The model DNP is represented schematically in left frame of Figure 1. It is built from a central core and Nch ligands that are free to surf on its surface. Each ligand is terminated by a LC group. The model equally may represent a generation three carbosilane dendrimer [14, 15].

Figure 1.

Model DNP (left frame). Definition of sizes and distances for pairs of interacting beads (right frame).

The model contains spherical beads (designated thereafter as ‘sp’) and the spherocylinder (‘sc’) ones that mimic LC groups. The sphere-sphere interaction is soft repulsive of quadratic form


where rij=rij/σijis scaled distance between centres of ith and jth sphere and mixing rules σij=σi+σj/2are implied (see right frame of Figure 1). Umaxspspis the same for all interacting spheres.

Spherocylinders are of breadth Dand of elongation L/Dand interact via the soft anisotropic potential of Lintuvuori and Wilson [16]:


where dij=dij/Dis the dimensionless nearest distance between the cores of spherocylinders (see right frame of Figure 1), dcis the orientation dependent cutoff. Attractive interaction has the form:


and depend on orientations êiand êjof the long axes of spherocylinders and the unit vector r̂ijthat connects their centres [16]. P2x=1/23x21is the second Legendre polynomial, energy parameters Uattr, ϵ1and ϵ2are given below. The inclusion of the attractive contribution in Eq. (2) shifts the region for the LC stability towards smaller elongations L/D3compared to the case of purely repulsive interactions, where, typically, L/D68[14].

Mixed nonbonded interactions are evaluated in a similar way to Eq. (1)


where dij=dij/σijis a dimensionless distance between the centre of the ith sphere and the core of the jth spherocylinder (see right frame of Figure 1), with the scaling factor σij=σi+D/2.

Intramolecular interactions for the model include bond and angle interactions terms


where liis the instantaneous bond length, l0kis the set of effective bond lengths. Ligands stiffness is adjustable via the magnitude of ka, and θ0is set equal to π.

The force-field parameters are based on the coarse-graining of the LC dendrimer first performed and described in detail in Ref. [14]. Bulk behaviour of such macromolecules are very similar to that of the DNPs [17]. The soft-core diameters of spherical beads are: σ=21.37, 6.23 and 4.59 Å for the core, first ligand bead and the following ligans beads, respectively. Their respective masses: 62.441025, 2.201025and 0.701025kg. LC beads dimension are D=3.74Å, L/D=3, their mass is 3.941025kgand the moment of inertia: 6.001024kg. The energy parameters Umaxspsc, Umaxspscand Umaxspscare all equal to 701020J. LC-LC specific energy parameters are Uattr=15001020, ϵ1=1201020and ϵ2=1201020J. The bond lengths are: 14.9, 3.60 and 3.62 Å between the core-first spacer bead, first-second spacer bead and between following spacer beads, respectively. The spherocylinder is attached to the centre of its nearest spherical cap with the bond length of 2.98 Å. The bond interaction spring constant is 501020Jfor all bonds. The pseudo-valent angle spring constant is 201020J.

The simulations are carried out with the GBMOLDD program extended to the case of the NPxPyPzTensemble [18]. We use a single Nóse-Hoover thermostat in most cases for both translational and rotational degrees of freedom. In quenching or rapid heating runs the velocity rescaling was used instead. The timestep Δt20fswas used for velocity rescaling runs, whereas smaller timesteps Δt1015fswas required for the runs with the thermostat. To control pressure three barostats are used [18].

3. Temperature driven morphology changes

In this section we consider temperature driven transitions between ordered and disordered morphologies observed in the DNPs melt. The experimental evidence for such transitions are found in Refs. [4, 5, 6, 7, 8, 9, 17] and indicate a close relation between the average DNPs shape and the type of the ordered morphology. Therefore, we attempt to steer the self-assembly towards particular morphology by influencing the DNPs shape. To do so we impose the orientation field of a given symmetry, which acts on the LC beads. It is introduced via the following energy term


where eiis the orientation of ith LC bead, iis the direction of the field, whereas field strength F provides either uniaxial (F>0) or planar (F<0) preferred orientation of the LC beads. In turn, in these cases we expect, respectively, rod-like or disc-like conformations of DNPs.

The case of the uniaxial field is considered first. Applied to the isotropic morphology of 100 DNPs in bulk using the NPxPyPzTensemble at P=50atmand T=520K, the field of the strength F=21020Jand with icollinear to the Z-axis, it induces formation of a monodomain smectic A (SmA) morphology. It was studied then with the field switched off and the melt equilibrated for 2040nsat selected temperatures within a range of T350K650K. Its appearance is shown in Figure 2 indicating a lamellar structure with alternating layers of cores, ligand chains and LC beads (shown in left frame), in agreement with the experimental studies [5, 17]. The two-dimensional arrangement within layers is shown in the right frame.

Figure 2.

Snapshots of the SmA phase at T=470K including typical rod-like molecular conformation.

We examined spatial distribution of the DNPs cores next. The form of their radial distribution function gris similar within a range of T=350470Kindicating two maxima: one at r27Åfor the short-range order of cores within each layer, and another at r65Å, related to the interlayer distance (see, Figure 3, left frame). These characteristic distances are model dependent and can be examined by evaluting pair distribution functions gzrand gxyr. The former yields an interlayer distance at 59Å(Figure 3, right frame). The latter indicates some degree of local (but not long-ranged) positional order within the layers (Figure 3, middle frame), characteristic of a 2D liquid. Therefore the phase is identified as the smectic A (SmA) [5, 17].

Figure 3.

Radial distribution functions gr (left frame), gxyr (middle frame) and gzr (right frame) for the DNPs cores in the SmA phase at T=350K (solid lines) and at T=470K (dashed lines).

Application of a field with a planar symmetry (F<0) along the Zaxis, leads to the formation of a columnar (Col) morphology. The defect-free monodomain sample is achieved by using the field of a moderate magnitude, F=21020Japplied to the isotropic melt at P=50atmand T520K. It was subsequently studied in a series of runs at the temperatures in the range T350K650Kand is shown at T=490Kin Figure 4. It displays disc-like DNPs stacked into columns, whereas the columns theirselves are arranged hexagonally, in accordance with the experimental observations [4].

Figure 4.

Snapshots of the Col phase at T=490K including typical disc-like molecular conformation.

The structure of Col phase is analysed via the pair distribution functions, shown in Figure 5. The hexagonal arrangement of columns is clearly indicated in the form of gxyr, while the peaks in gzrallow the distance between dendritic cores within a column to be estimated as 18Åfor the Col phase at T=350Kand ∼20 Å at T=490K.

Figure 5.

Radial distribution functions gr (left frame), gxyr (middle frame) and gzr (right frame) for the DNPs cores in the Col phase at T=350K (solid lines) and at T=470K (dashed lines).

Upon heating up, both ordered morphologies, SmAand Col, undergo transitions to the spatially disordered phase as monitored via the order parameter SNof mesogens with respect to the axis i(nematic director for the SmAand the columns vector for the Col symmetry)


where P2xis the second Legendre polynomial and averaging is performed on all mesogens in a melt. As far as in this section the field is always collinear to the Z-axis, the notation Szis also used. Shape asymmetry of DNPs, a, is obtained from the components of the gyration tensor Gαβ


where α,βstay for the Cartesian axes, ri,αand Rαare the coordinates of the individual beads and of the DNP’s centre of mass, respectively, Rg2=ααGααis the squared radius of gyration. Each spherocylinder bead is replaced by a line of four centres. The shape anisotropy a is positive for the prolate shape and negative for the oblate one.

Sz, aand the system density ρare all shown as the functions of the temperature in Figure 6 for both SmA(left frame) and Col (right frame) morphologies. At lower temperatures both Szand a are non-zero, indicating the SmAand Col phases and quantifying the amount of orientation order and DNPs shape asymmetry at each temperature. Transition to the disordered isotropic morphology occurs at about T=490500K, where both Szand a simultaneously turn to zero.

Figure 6.

Evolution of the density, shape asymmetry and the order parameter at the SmA-I and Col-I phase transitions.

The simulations reproduce a strong link between the shape of a DNPs and the type of bulk morphology, as previously observed experimentally [4, 17]. In particular, rod shapes are compatible with the SmAsymmetry (Figure 2), discotic shapes are compatible with Col symmetry (Figure 4) and an spherical molecular shape is found in an isotropic state. This confirms the validity of the present model, which catches main features of the DNPs bulk assembly. The prolate-to-isotropic and oblate-to-isotropic shape transition occur simultaneously with vanishing the respective order parameter, as shown in Figure 6. The DNPs demonstrate shape bistability leading to the possibility to observe different symmetries at the same thermodynamics condition. As will be shown in the following section, this feature is dependent of the density of ligands defined via Nch.

4. The role of the details of the molecular architecture

4.1 Variation of a grafting density

In Section 3 the number of ligands in DNPs was kept fixed at Nch=32. The decoration density, however, can change the shape of the supermolecule leading to different types of self-organised structure [4, 6, 17]. On the first sight, the most favourable conformation at any Nchcould be estimated from space-filling considerations and then the type of self-assembled morphology could be predicted. This was attempted in [19] but showed that the temperature effects also plays a crucial role and the symmetry of the self-assembled morphology at each conditions is the result of a delicate balance between the enthalpic and entropic contributions to the free energy.

Based on experimental findings [5, 8, 9, 17, 20], we expect to observe the sequence of SmA, Col and cubic morphologies upon the increase of Nch. For each Nch, we performed aided self-assembly runs of duration 20nsat T=520Kwith the timestep of 20fsin the NPxPyPzTensemble. Both uniaxial and planar fields were used. Then, the field was switched off and the system was equilibrated at T=450K, about 50 K below the transition. Besides these, the spontaneous self-assembly runs were performed, too. In the latter, the temperature of the melt was reduced linearly from T=500Kdown to 450 K during first 20ns(cooling rate is 2.5 K/ns), followed by another run for 20nsat fixed T=450K. As the result, relatively defect-free smectic layers are found for the cases of Nch=12and Nch=20only, whereas at Nch24we obtained polydomain layered structures with either one type of domains (of smectic or discotic type) or a mixture of both. In the field-aided self-assembly the external field (6) was employed, similarly as described in Section 3. Both runs with uniaxial (F>0) and planar (F<0) fields were undertaken for each case of the number of ligands Nchbeing considered.

Due to a relatively small system size, firm identification of smectic and discotic domains turned to be quite problematic, the visual inspection indicates domains of just a few DNPs. Therefore, we opted to analyse the distribution of the DNPs shape asymmetries instead, assuming that each rod-shaped DNP is a part of a smectic domain and a disc-shaped one—of the discotic one, as demonstrated earlier in Section 1. To distinguish between two shapes, we introduce molecular ‘roddicity’ ar(always positive) and molecular ‘discoticity’ ad(always negative):


for each DNP. Here λmax, λmedand λminare the maximum, medium and minimum eigenvalues of the gyration tensor (8), respectively. If, for given DNP, ar>ad, then it is classified as a rod with the shape anisotropy a=ar, otherwise—as a disc with a=ad. At each time instance, the DNP melt splits into rods and discs subsystems, with their fractions termed as frand fd, respectively.

The histograms for the distribution paare built over all the DNPs in the system averaged over time trajectory and are shown in Figure 7 for selected values of Nch. Both the cases of unaided (left frame) and field-aided self-assembly (middle and right frames) are displayed. In the case of spontaneous self-assembly, rods and discs always coexist and the distributions of their shape asymmetry are relatively broad. With the increase of Nch, two maxima gradually merge into a spherulitic shape from both sides of a=0(at about Nch=64and higher). The histograms at the field-aided self-assembly are much narrower. In the case of uniaxial field, the discotic conformations are completely eliminated (except the case of Nch=48where smectic phase is not observed any more), as these are incompatible with the 1D symmetry of the aiding field. In the case of planar filed, which has a 2D symmetry, the rod-like conformations are not eliminated and do appear within XYplane, and are, in fact, the dominant ones at smaller values of Nch. With the increase of Nchabove 24the disc-like conformations are dominating.

Figure 7.

Histograms for the distributions of shape anisotropy pa′ at spontaneous self-assembly (left image), uniaxial field aided self-assembly (middle image) and planar field aided self-assembly (right image).

Resulting fractions of rods and discs, frand fd, are shown in Figure 8 as functions of Nchfor both types of runs. Left frame (unaided self-assembly case) indicates the broad region of the rod-disc coexistence at intermediate values of Nch. At Nch=64the system approaches symmetric case with both conformations transforming into a spherulitic shape. The right frame contains data for frfor uniaxial field aided self-assembly and data for fdfor planar field aided self-assembly, therefore, fr+fd1as both are obtained for different conditions. The shapes of both curves are much steeper as compared to the left frame plot indicating the possibility to control the molecular conformation by means of initial field of appropriate symmetry. Therefore, there is some interval of grafting density, where DNPs exhibit rod-disc shape bistability and, hence, the SmA-Col morphology bistability is observed in the melt. Spontaneous self-assembly yields the polydomain structure with both SmAand Col fragments. SmAor Col morphology can be made dominant or, at least, enhanced by an external field of appropriate symmetry [19]. One of such cases is discussed in detail in Section 5. These results reproduce the main trends of the DNPs self-assembly as previously seen experimentally [5, 8, 9, 17, 20].

Figure 8.

Fraction of rods and discs for unaided self-assembly (left frame, fr+fd=1). The same properties are shown on the right but fraction of rods is shown for uniaxial field aided runs and fraction of discs—for planar field aided runs.

4.2 Variation in liquid crystalline groups attachment

The model for DNP considered so far, is characterised by a longitudinal attachment of the terminal LC beads (Figure 1) and is found to self-assemble into the SmAand Col morphologies discussed above. With respect to the applications, the anisotropy in material mechanical properties depends on the spatial arrangement of the DNPs cores, whereas its optical properties—on the orientations of the LC beads. The SmAmorphology is lamellar with the layers normal a, and also optically uniaxial, characterised by the nematic director ncollinear to a(see Figure 2). The Col morphology is a supramolecular assembly of columns, where each column is a stack of disc-shape DNPs itself. The arrangement of columns is uniaxial and is described via vector a. However, there is no global nematic order of LC beads in this case, as their orientations are distributed radially in a plane perpendicular to a, see Figure 4. Therefore, optical response of the Col morphology is essentially different from that of the SmAone, bearing some analogy in the difference in opto-mechanical applications of the main- and side-chain LC architectures, see e.g. [21].

It is evident from Figure 4, that flat radial orientations of the LC beads in Col phase is the result of the energy penalty associated with bending of its host ligands. Lateral attachment of the LC bead, see, Figure 9, left frame, is also possible from the view of chemical synthesis [17], and may open up a possibility of greater orientation freedom of the LC beads. In particular, one may expect the optically uniaxial columnar morphology shown in Figure 9, right frame, or, some other morphologies as well. With this aim we modified the DNP model accordingly. The orthogonality of LC groups to the spacer-LC bond is maintained by employing two potentials: harmonic bond between the last monomer of the spacer and the centre of the spherocylinder (the bond length is 0.3 nm) and harmonic angular potential between the latter bond and the long axis of the LC bead (equilibrium angle is θ0=π/2).

Figure 9.

Model DNP with lateral attachment of LC beads (left frame) and uniaxial Col morphology (right frame).

Only the case of Nch=32ligands is considered and the self-assembly runs are all field-assisted. The field is directed along the Zaxis and its magnitude is set at F=241020J. The simulations are performed for the range of pressures P=50200atm, the duration of each run is 106molecular dynamics steps with a time step of Δt=20fs. For all values of the pressure within this interval the system indeed assemble into the uniaxial hexagonal columnar morphology depicted schematically in Figure 9 (right frame) and referred thereafter as uColh. We use this morphology as an initial state and perform a series of subsequent runs at selected temperatures with the field switched off are aimed on examining the temperature stability of this phase.

For the analysis of the structural changes we consider a number of order parameters. Besides the orientation order parameter SN(7) we also introduce the hexagonal order parameter SHwithin the XY plane and columnar order parameter SCdefined as follows


Here φkis the polar angle of the bond between kth and ith DNPs, where the summation is done over all kth DNPs that belong to a first coordination sphere of ith DNP (see Figure 10(a)), j=1, Nc,iis the number of DNPs such as their centres are found inside this cylinder of radius Rc=1nmdrawn around the core of the ith DNP along the Zaxis (see Figure 10(b)), and Nmaxis the normalisation factor introduced for the sake of convenience.

Figure 10.

(a) Flat hexagonal cluster of DNPs and the definition of bond angles φk. (b) A column of stacked DNPs. Order parameters SN, SH and SC upon heating the uColh.

Evolutions of all order parameters, SN, SHand SCupon the increase of the temperature are shown in Figure 10. At low temperatures, T<350K, the values of all three order parameters are essentially non-zero indicating the uniaxial hexagonal columnar phase uColhcharacterised by uniaxial nematic order, high columnarity and high hexagonal order of the columns. With the increase of the temperature all order parameters gradually decay, but while SNdecays very fast (almost linearly, in contrast with Figure 6), the essential delay is observed in the decrease of SHand SC. As the result, there is a temperature range around T=450K, where the value of SNdropped to about 0.1 (typical for the isotropic phase), but SHand SCare still almost the same as in the uColhmorphology at T=300400K. This is a columnar morphology with as a weak hexagonal order, wColhbut still with considerably ordered LC beads. The transition from uColhto wColhis gradual and the boundary between both shown in Figure 10 is for rather illustrative purpose. With further heating of the system, all order parameters SHand SCdrop to their minima at approximately T480Kindicating disordered morphology [22]. Therefore, the simulations predict two novel discotic morphologies for the DNPs with lateral attachment of LC groups: one characterised by an uniaxial, another—by random orientations of the LC groups.

5. Photo-assisted self-assembly for the azobenzene-decorated nanoparticles

Results presented in Sections 3 and 4 indicate that efficient self-assembly of model DNPs into monodomain morphologies faces certain difficulties. These are not the artefacts of the model or the simulation approach, but reflect physical properties of the LC polymers, which are relatively viscous and characterised by slow relaxation and strong tendencies for metastability attributed to the presence of a transient network [23]. One of the ways to stimulate/control the self-assembly of DNPs is to use the light [24, 25, 26, 27, 28]. To utilise this approach, the chromophores (e.g. azobenzene, cinnamoyl, diarylethene dithiophenols, etc.) are incorporated into the DNPs ligands. These effects share similarities with photoinduced structural changes in azobenzene-containing side-chain polymers.

Azobenzene chromophore is one of the most widely used, it exists in the trans- and cis-forms. The trans-isomer is prolate and has the LC properties, the cis-one is a non-LC. The photoisomerisation occurs between these forms when illuminated by a light with suitable wavelength. The trans-cis photoisomerization is angle-selective, the reverse one is not. If the absorption bands for both transformations overlap, the continuous trans-cis-trans isomerization cycles take place resulting in the reorientation of the trans-isomers predominantly perpendicularly to the light polarisation axis (‘orientation hole-burning’ or Weigert effect). Both polarised and unpolarized light can be used to photo-align the trans-isomers. In the latter case, as pointed out by Ikeda [29], ‘only the propagation direction is, in principle, perpendicular to the electric vector of the light. Thus, when unpolarized light is employed, it is expected that the azobenzene moieties become aligned only in the propagation direction of the actinic light’. This situation opens up a possibility for uniaxial alignment of azobenzenes in the direction of the light propagation and is exploited for our model DNPs melt.

We follow Ref. [31] in complementing the deterministic part of the simulations, provided by molecular dynamics, with stochastic part, which describes photoisomerisation of chromophores on a coarse-grained level. We represent trans-isomers via coarse-grained t-beads and their cis-counterparts—by c-beads. Both are considered of the same spherocylinder shape, but t-t interaction is of the LC type (described by a potential (2)), whereas the t-c and c-c interactions are non-LC soft repulsive (described by a potential of the similar form as (1) and (4)). The quantum mechanical nature of photoisomerization [32, 33] is accounted for implicitly, by applying the kinetic equations of a general form for the probabilities of the transitions between the t- and c-state of each ith bead:


where îis the unit vector collinear with the light polarisation axis, and ptand pcare the respective transition rates that depend on the chemical details of the chromophore group and the intensity and the wavelength of the illumination. The switch of the state is attempted for each chromophore at each MD step. Selective absorption of the light by the azobenzene chromophores is reflected in the angular dependence of the transition probability pitc. Photoisomerisation also involves the random change of the chromophore orientation [31, 32]. The choice being made for pt=0.001and pt/pc=0.5is justified in detail in [30] and leads to the concentration of c-beads of 5–10% in a photostationary state.

As was discussed in Section 3, the model DNP is found to exhibit rod-disc shape bistability at a wide interval of Nch=2456(see, Figure 8) resulting in a polydomain mixture of both SmAand Col domains upon cooling it down from isotropic state. To avoid this uncertainty we use the DNPs with Nch=12, in which case only SmAmorphology is found.

To monitor the level of SmAordering in a system, we consider a set of relevant characteristics. Each DNP is considered in terms of its equivalent ellipsoid provided by the gyration tensor components (8). The asphericity, A, is defined as an average A=3/2λmax+λmed+λmin/Rg21/2. Orientation order of DNPs is defined as S2=P2ÊN̂, where Êis the orientation of DNP long axis and the nematic director N̂is evaluated in a way usual for LC. Global smectic order in a system is linked to the level of ‘lamellarity’ in the arrangement of DNPs cores and is quantified via the amplitude of the density wave along the layer normal. It is evaluated by finding the maximum of the expression


as a function of p (here i=1). The maximum position pprovides the pitch of the SmA phase, whereas the smectic order parameter is Ss=Sp. To examine polydomain structure in a system we split it into separate SmA clusters (if more than one exist) and calculate their number per molecule, Nc, and the reduced maximum cluster size Mc, which is the number of DNPs in a largest cluster divided by the total number of DNPs. For the monodomain morphology one has Nc0and Mc1, whereas for a highly polydomain phase: Nc1and Mc0. The magnitude of all characteristics, A, S2, Ss, Ncand Mc, are, therefore, restrained to the interval between 0 and 1.

We start from the heating runs. The initial monodomain SmAmorphology was prepared with the aid of the orienting field (6). The changes in its properties are monitored then in a temperature interval from 400 to 550Kwith the field switched off. These changes are shown in Figure 11. As follows from there, the system undergoes sharp changes at T510K: both order parameters S2and Sssharply drop to zero indicating the presence of the order-disorder transition. The synchronicity in the evolution of A, S2and Ssindicates the absence of the purely nematic phase (S2>0and Ss=0) and a strong relation between the molecular shape and the symmetry of the ordered morphology, similarly to the case discussed in Section 3. Therefore, the transition that occurs at Tis the SmA-I transition. The behaviour of Ncand Mcindicates the monodomain morphology at T=400500Kwhich is transformed into a highly polydomain one in a course of a transition.

Figure 11.

A, S2 and Ss (left frame) and Nc and Mc (right frame) vs. T. Initial state: monodomain SmA, heated up to T indicated in each plot. Reprinted with permission from [30]. Copyright (2016) American Chemical Society.

We attempt next the reverse transition: a spontaneous self-assembly of the SmAmorphology out of the isotropic obtained upon cooling down the initial isotropic state. Various cooling rates ranging from 0.37to 4.5K/nsare used (note, that because of coarse-graining, the time scale is essentially contracted comparing with the real systems). The results are shown in Figure 12. Their behaviour demonstrate that at low enough cooling rate of 0.37K/nsA, Ncand Mcclosely follow their respective curves shown in Figure 11. The order parameter Ssdo not quite reach its respective value found in Figure 11 but is, nevertheless, close. At higher cooling rate of 1.12K/nsthe system is trapped in a polydomain state with Nc0.2(about 40 individual clusters). The maximum cluster size is still quite large, Mc0.6(about 120 DNPs). This indicates the coexistence of one large and about 40 smaller clusters. At higher cooling rate the system is split into a larger number of smaller clusters, as indicated by Ncand Mc.

Figure 12.

Temperature dependence of Ss, Nc and Mc. Initial state: isotropic system, cooled down from 490 K to a given T. Cooling rates: Black squares: 0.37K/ns, blue circles: 0.75K/ns, red triangles: 1.12K/ns, and orange diamonds: 4.5K/ns. No illumination is applied. Adapted with permission from [30]. Copyright (2016) American Chemical Society.

To examine the origin of this metastable state, we performed analysis of the translational and rotational mobility of DNPs in a temperature interval of 350500K. The system was quenched at each temperature for 30 ns and the initial 18nsof each the run were analysed by splitting it into time blocks at time instances tkof equal duration δt=tktk1=0.2ns. We define instantaneous translational and rotational diffusivities


at each tk, being averaged over initial 18ns; they provide short-time translational D=dtkand rotational R=rtkdiffusion coefficients of the DNPs at the early stage of a self-assembly. The temperature dependence of both is shown in Figure 13 as black legends. Both characteristics are high at 500 K and decay sharply as the temperature decreases, indicating a slowdown of translational and rotational mobility of DNPs. Hence, at lower temperatures one observes physically cross-linked domains, similarly to the case of side-chain LC polymers [23]. To quantify orientational arrest of chromophores at lower temperatures we estimated their rotation relaxation time as the function of the temperature from the time autocorrelation function êtê0for their orientations assuming its exponential decay

Figure 13.

Chromophores rotation relaxation time trot and short-time translational D and rotational R diffusion coefficients. Illuminated melt (red legends) vs. non-illuminated melt (black legends). Reprinted with permission from [30]. Copyright (2016) American Chemical Society.


The evaluation time interval is restricted again to the initial 18nsof each run. The dependence of troton temperature shows an essential increase of 11.5order of magnitude upon lowering the temperature from 500 down to 350K, see Figure 13.

Let us project now the effects of the illumination on the self-assembly of the DNPs. There are two known types of photomodulation in azobenzene-containing systems: (A) weakening (or elimination) of the LC order, and/or (B) order-order transition [29]. The effect (A) should weaken the interdomain links and increasing translation and rotation freedom of the DNPs. This is found in the performed simulations, as is seen from the Figure 13 (red legends). While the differences between illuminated and no illuminated cases are negligible around 500 K, trotunder illumination decreases by a factor of 2.5 at 350 K. Higher mobility of DNPs under illumination is also detected at higher valued of Dand R, especially at 460500K.

These effects of illumination manifest theirselves on the self-assembly in a course of cooling runs. The results for Ss, Ncand Mcunder illumination are shown in Figure 14. Similar behaviour is found for most characteristics as compared to the no illumination case, Figure 12, but at different respective cooling rates. The illumination reduces the self-assembly time-scale on average by a factor of 3–4 as compared to the no illumination case. Let us note again that the model time-scale is contracted comparing with the real units due to coarse-graining, therefore we emphasise on the relative speed-up of a self-assembly but not on the absolute cooling rates.

Figure 14.

Temperature dependence of Ss, Nc and Mc. Initial state: isotropic system, cooled down from 490 K to a given T indicated in the plot. Cooling rates: Blue circles: 0.75 K/ns, orange diamonds: 3.75 K/ns, pink triangles: 18 K/ns, green discs: 25 K/ns. The case with illumination. Adapted with permission from [30]. Copyright (2016) American Chemical Society.

To have additional proof on the ability of the illumination to aid self-assembly of DNPs, we performed quenching runs. With no illumination applied, such runs end up in a glassy-like metastable state with no evidence of the orientational or positional order. The situation is markedly different under illumination, as indicated in Figure 15 for a set of properties A, S2, Ss, Nc and Mc. Comparing with Figure 11 for heating runs, one finds the curves of similar respective shapes albeit all shifted to the lower temperatures. Therefore, computer simulation studies indicate the possibility for the SmA phase self-assemble by quenching the system in a broad temperature interval under illumination, otherwise impossible for the non-illuminated system. Although, these simulation findings cannot be compared straightaway with the particular experiments, they are very much in-tune with the general applications of the azobenzene chromophores to control LC alignment, aggregation and self-assembly in the systems of DNPs [24, 25, 26, 27, 28].

Figure 15.

A, S2 and Ss (left frame) and Nc and Mc (right frame) vs. T. Initial state: isotropic system quenched to T under illumination. Reprinted with permission from [30]. Copyright (2016) American Chemical Society.

6. Conclusions

Coarse-grained modelling technique maps chemical details of nanoparticles onto physical set of structure elements. In particular, the dimensions and the type of the core (metal or organic) are examined, the number and the properties of ligands (e.g. flexible or rigid, aliphatic or aromatic, etc.) are analysed, as well as the nature and properties of terminal functional groups (if any). These structure elements are replaced then by classical mechanics objects of appropriate shape that are connected via Hookean springs and interact via van der Waals forces. The parameters of the effective interaction potentials are found from related atomistic simulations using various techniques (e.g. force matching). While lacking chemical specificity of particular type of decorated nanoparticles, such models focused instead on important generic features of a whole class of underlying objects. More importantly, such simplification enables to reach required time- and length-scales of the self-assembly.

In this study we showed that such type of computer simulations is able to reproduce a wide spectra of experimentally observed features related to the self-assembly of decorated nanoparticles. Namely, we demonstrate self-assembly of model nanoparticles into lamellar smectic and columnar morphologies, reproduce and studied the shape-morphology relation, examine preference of ordered phase depending on the decoration density, study temperature driven order-disorder transitions. New optically uniaxial hexagonal phases are predicted for the case of laterally functionalised ligands.

Being viscous fluids, the melts of decorated nanoparticles are prone to slow relaxation and associated with this metastability. This restricts their application where fast regrowing/restructuring of the nanomaterial is needed in response on external stimulus. We demonstrate that the illumination of a suitable choice can be used for the speed-up of the self-assembly of the chromophore-containing nanoparticles.

The work can be extended in a number of ways. First, the use of the more accurately parameterised interaction potentials: if not for predicting the exact phase boundaries for new morphologies, then at least to show a right direction of where to search for them. Second, to consider specific (magnetic or non-metal) nanoparticles. Finally, the use of larger system sizes is always preferable. In conclusion: despite the increase in computing power, the coarse-grained picture always provides a valuable insight on the problem by focusing on its most relevant features.


The author thanks M.R. Wilson, J. Lintuvuori, M. Saphiannikova and A. Slyusarchuk for fruitful collaboration on the topic, as well as S. Sokołowski, S. Santer and E. Gorécka for stimulating discussions.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Jaroslav Ilnytskyi (May 22nd 2020). Self-Assembly of Nanoparticles Decorated by Liquid Crystalline Groups: Computer Simulations [Online First], IntechOpen, DOI: 10.5772/intechopen.89682. Available from:

chapter statistics

10total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us