This Letter investigates the defect structures that arise between two colloidal spheres immersed in a nematic liquid crystal. Molecular simulations and a dynamic field theory are employed to arrive at molecular-level and mesoscopic descriptions of the systems of interest. At large separations, each sphere is surrounded by a Saturn ring defect. However, at short separations both theory and simulation predict that a third disclination ring appears in between the spheres, in a plane normal to the Saturn rings. This feature gives rise to an effective binding of the particles. The structures predicted by field theory and molecular simulations are consistent with each other.

We report a methodology that permits quantitation of the azimuthal anchoring energy of the nematic liquid crystal 4-cyano-4'-pentyl-biphenyl on surfaces patterned with oligopeptides. The oligopeptide (IYGEFKKKC), an optimized substrate for the Src protein kinase, was covalently immobilized via the terminal cysteine to monolayers of amine-terminated tetra(ethylene glycol) formed on gold films. The measurements of anchoring energies, which were based on a torque-balance method, revealed a systematic decrease in anchoring energy from 3.7 +/- 0.6 microJ/m2 with increasing surface density of oligopeptide. We calculate that a mass density of oligopeptide of less than 1 ng/cm2 can lead to a measurable change in the anchoring energy of the nematic liquid crystal. These results suggest that measurements of anchoring energies of liquid crystals on surfaces may offer the basis of quantitative and label-free methods for detecting biomolecules on surfaces.

We report measurements of the orientations and azimuthal anchoring energies of the nematic liquid crystal 4-cyano-4'-pentylbiphenyl (5CB) on polycrystalline gold films that are deposited from a vapor at an oblique angle of incidence and subsequently decorated with organized monolayers of oligomers of ethylene glycol. Whereas the gold films covered with monolayers presenting tetra(ethylene glycol) (EG4) lead to orientations of 5CB that are perpendicular to the plane of incidence of the gold, monolayers presenting tri(ethylene glycol) (EG3) direct 5CB to orient parallel to the plane of incidence of the gold during deposition of the gold film. We also measure the azimuthal anchoring energy of the 5CB to be smaller on the surfaces presenting EG3 (3.2 +/- 0.8 microJ/m2) as compared to EG4 (5.5 +/- 0.9 microJ/m2). These measurements, when combined with other results presented in this paper, are consistent with a physical model in which the orientation and anchoring energies of LCs on these surfaces are influenced by both (i) short-range interactions of 5CB with organized oligomers of ethylene glycol at these surfaces and (ii) long-range interactions of 5CB with the nanometer-scale topography of the obliquely deposited films. For surfaces presenting EG3, these short- and long-range interactions oppose each other, leading to small net values of anchoring energies that we predict are dependent on the level of order in the EG3 SAM. These measurements provide insights into the balance of interactions that control the orientational response of LCs to biological species (proteins, viruses, cells) on these surfaces.

We report results for the potential of mean force (PMF) and the defect structures that arise when spherocylindrical nanoparticles are immersed in a nematic liquid crystal. Using a dynamic field theory for the tensor order parameter Q of the liquid crystal, we analyzed configurations, including one, two, and three elongated particles, with strong homeotropic anchoring at their surfaces. For systems with one nanoparticle, the most stable configuration is achieved when the spherocylinder is placed with its long axis perpendicular to the far-field director, for which the defect structure consists of an elongated Saturn ring. For systems with two or three nanoparticles with their long axes placed perpendicular to the far-field director, at small separations the defect structures consist of incomplete Saturn rings fused with new disclination rings orthogonal to the original ones, in analogy to results previously observed for spherical nanoparticles. The shape of these orthogonal rings depends on the nanoparticles' configuration, i.e., triangular, linear, or parallel with respect to their long axis. A comparison of the PMFs indicates that the latter configuration is the most stable. The stability of the different arrays depends on whether orthogonal disclination rings form or not, their size, and the curvature effects in the interparticle regions. Our results suggest that the one-elastic-constant approximation is valid for the considered systems; similar results were obtained when a three-constant expression is used to represent the elastic free energy. The attractive interactions between the elongated particles were compared to those observed for spheres of similar diameters. Similar interparticle energies were observed for linear arrays; in contrast, parallel and triangular arrays of spherocylinders yielded interactions that were up to 3.4 times stronger than those observed for spherical particles.

The structure of polymers in the vicinity of spherical colloids is investigated by Monte Carlo simulations and integral equation theory. Polymers are represented by a simple bead-spring model; only repulsive Lennard-Jones interactions are taken into account. Using advanced trial moves that alter chain connectivity, depletion and packing effects are analyzed as a function of chain length and density, both at the bond and the chain level. Chain ends segregate to the colloidal surface and polymer bonds orient parallel to it. In the dilute regime, the polymer chain length governs the range of depletion and has a negligible influence on monomer packing in dense polymer melts. Polymers adopt an ellipsoidal shape, with the larger axis parallel to the surface of the particle, as they approach larger colloids. The dimensions are perturbed within the range of the depletion layer.

Monte Carlo simulations and dynamic field theory ͑DyFT͒ are used to study the interactions between dilute spherical particles, dispersed in nematic and isotropic phases of a liquid crystal. A recently developed simulation method ͑expanded ensemble density of states͒ was used to determine the potential of mean force ͑PMF͒ between the two spheres as a function of their separation and size. The PMF was also calculated by a dynamic field theory that describes the evolution of the local tensor order parameter. Both methods reveal an overall attraction between the colloids in the nematic phase; in the isotropic phase, the overall attraction between the colloids is much weaker, whereas the repulsion at short range is stronger. In addition, both methods predict a new topology of the disclination lines, which arises when the particles approach each other. The theory is found to describe the results of simulations remarkably well, down to length scales comparable to the size of the molecules. At separations corresponding to the width of individual molecular layers on the particles' surface, the two methods yield different defect structures. We attribute this difference to the neglect of density inhomogeneities in the DyFT. We also investigate the effects of the size of spherical colloids on their interactions.

We analyze the response of a nematic liquid-crystal film, confined between parallel walls, to the presence of nanoscopic particles adsorbed at the walls. This is done for a variety of patterns of adsorption ͑random and periodic͒ and operational conditions of the system that can be controlled in experimental liquid-crystal-based devices. We compute simulated optical textures and the total optical output of the sensor between crossed polars, as well as the correlation function for the liquid-crystal tensor order parameter; we use these observables to discuss the gradual destruction of the original uniform orientation. For large concentrations of particles adsorbed in random patterns, the liquid crystal at the center of the sensor adopts a multidomain state, characterized by a small correlation length of the tensor order parameter, and also by a loss of optical anisotropy under observation through crossed polars. In contrast, for particles adsorbed in periodic patterns, the nematic at the center of the cell can remain in a monodomain orientation state, provided the patterns in opposite walls are synchronized.

Monte Carlo simulations and dynamic field theory are used to study spherical particles suspended in a nematic liquid crystal. Within these two approaches, we investigate the binding of the defects to the particles, the adsorption of a particle at a solid surface, and two particles interacting with each other. Quantitative comparisons indicate good agreement between the two approaches. A Monte Carlo method based on the combination of canonical expanded ensemble simulations with a density-of-state formalism is used to determine the potential of mean force between one particle and a hard wall. On the other hand, the potential of mean force is evaluated using a dynamic field theory, where the time-dependent evolution of the second rank tensor includes two major aspects of liquid crystalline materials, namely the excluded volume and the long-range order elasticity. The results indicate an effective repulsive force that acts between the particle and the wall. Layer formation at the surface of the hard wall gives rise to local minima in the potential of mean force. The director profile for a particle at contact with a solid surface is characterized by a disclination line distorted and attracted towards the wall. The structure of the nematic for two particles at short distances is also investigated. Our results indicate a structure where the two particles are separated by a circular disclination line. The potential of mean force associated with this configuration indicates an effective attractive interaction between the two particles.