Structure and Dynamics of Octamethylcyclotetrasiloxane Confined between Mica Surfaces

Correlation between ordering and shear thinning in confined OMCTS liquids

Despite the extensive research that has been conducted for decades on the behavior of confined liquids, detailed knowledge of this phenomenon, particularly in the mixed/boundary lubrication regime, remains limited. This can be attributed to several factors including the difficulty of direct experimental observations of the behavior of lubricant molecules under non-equilibrium conditions, the high computational cost of molecular simulations to reach steady state, and the low signal-to-noise ratio at extremely low shear rates corresponding to actual operating conditions. To this end, we studied the correlation between the structure formation and shear viscosity of octamethylcyclotetrasiloxane confined between two mica surfaces in a mixed/boundary lubrication regime. Three different surface separations corresponding to two-, three-, and five-layered structures were considered to analyze the effect of confinement. The orientational distributions with one specific peak for n=2 and two distributions, including a parallel orientation with the surface normal for n>2, were observed at rest. The confined liquids exhibited a distinct shear-thinning behavior independent of surface separations for a relatively low shear rate, \dot{\gamma}<10^{8} s^{-1}. However, the shear viscosities at \dot{\gamma}<10^{8} s^{-1} depended on the number of layered structures. Newtonian behavior was observed with a further increase in the shear rate. Furthermore, we found a strong correlation between the degree of molecular orientation and the shear viscosity of the confined liquids. The magnitude of the shear viscosity of the confined liquids can primarily be determined by the degree of molecular orientation, and shear-thinning originates from the vanishing of specific orientational distributions with increasing shear rate.

Structures of Nanoconfined Liquids Determined by Synchrotron X-ray Diffraction

We have successfully performed X-ray diffraction measurements of the liquids octamethylcyclotetrasiloxane (OMCTS, a quasi-spherical-shaped molecule) and n-hexadecane (a normal alkane) confined between mica surfaces at surface separation distances (D's) from 500 nm to the hard-wall thickness (1.9 nm for OMCTS and 1.0 nm for hexadecane). At all of the studied D's, we observed diffraction peaks corresponding to their mean intermolecular spacing at q = 8.6 nm^-1 (d = 0.73 nm) for OMCTS and q = 13.6 nm^-1 (d = 0.45 nm) for n-hexadecane. The peak intensity increased at D < ca. 50 nm for OMCTS even with the decreasing distance and exhibited a local maximum at D = 17-13 nm, indicating the sharp increase in the molecular order in this distance range. The peak intensities normalized by the D and I_normalized values of OMCTS and n-hexadecane were nearly constant at D's greater than 100 nm, though they appeared to increase slightly. The increase then became more significant with decreasing D below 100 nm, and finally the I_normalized values became 120 (for OMCTS) and 160 (for n-hexadecane) at the hard wall. These results clearly demonstrated the significant increase in the structural order of OMCTS and n-hexadecane under nanoconfinement, especially below 100 nm. The fwhm values of the peaks of OMCTS and n-hexadecane showed no significant change until small distances when the confinement effect was significant. These results indicated that the increase in the structural order should be mainly ascribed to the ordering of the molecules in the parallel plane in the enhanced layered structure formed under the confinement. The viscous parameters (b₂) of OMCTS and n-hexadecane obtained from the resonance shear measurement showed no increase at D's down to ca. 7 nm. This indicated that a certain ordering of the confined molecules was required for the observable increase in the viscosity.

Determination of the Free Energies of Mixing of Organic Solutions through a Combined Molecular Dynamics and Bayesian Statistics Approach

As new generations of thin-film semiconductors are moving toward solution-based processing, the development of printing formulations will require information pertaining to the free energies of mixing of complex mixtures. From the standpoint of in silico material design, this move necessitates the development of methods that can accurately and quickly evaluate these formulations in order to maximize processing speed and reproducibility. Here, we make use of molecular dynamics (MD) simulations, in combination with the two-phase thermodynamic (2PT) model, to explore the free energy of mixing surfaces for a series of halogenated solvents and high-boiling point solvent additives used in the development of thin-film organic semiconductors. Although the combined methods generally show good agreement with available experimental data, the computational cost to traverse the free-energy landscape is considerable. Hence, we demonstrate how a Bayesian optimization scheme, coupled with the MD and 2PT approaches, can drastically reduce the number of simulations required, in turn shrinking both the computational cost and time.