Computational chemistry calculations are broadly useful for guiding the atom-scale design of hard-soft material interfaces including how molecular interactions of single-component liquid crystals (LCs) at inorganic surfaces lead to preferred orientations of the LC far from the surface. The majority of LCs, however, are not single-component phases but comprise of mixtures, such as a mixture of mesogens, added to provide additional functions such as responsiveness to the presence of targeted organic compounds (for chemical sensing). In such LC mixtures, little is understood about the near-surface composition and organization of molecules and how that organization propagates into the far-field LC orientation. Here, we address this broad question by using a multiscale computational approach that combines density functional theory (DFT)-based calculations and classical molecular dynamics (MD) simulations to predict the interfacial composition and organization of a binary LC mixture of 4'-cyano-4-biphenylcarbolxylic acid (CBCA) and 4'-n-pentyl-4-biphenylcarbonitrile (5CB) supported on anatase (101) titania surfaces. DFT calculations determine the surface composition and atomic-scale organization of CBCA and 5CB at the titania surface, and classical MD simulations build upon the DFT description to describe the evolution of the near-surface order into the bulk LC. A surprising finding is that the 5CB and CBCA molecules adopt orthogonal orientations at the anatase surface and that, above a threshold concentration of CBCA, this mixture of orientations evolves away from the surface to define a uniform far-field homeotropic orientation. These results demonstrate that molecular-level knowledge achieved through a combination of computational techniques permits the design and understanding of functional LC mixtures at interfaces.
Keywords: anchoring energy; chemoresponsive systems; density functional theory; elastic energy; liquid crystals; molecular dynamics.