Ivan Lyubimov

Tunable molecular orientation and elevated thermal stability of vapor-deposited organic semiconductors

Significance Glasses are solids that lack the regular order of crystals. Organic glasses, when produced by deposition from the vapor, can exhibit high levels of molecular orientation that improve performance of devices such as organic light-emitting diodes. We show here that molecular orientation in such glasses is primarily controlled by the substrate temperature during deposition, suggesting that the performance of almost any device based upon amorphous organic materials might be systematically optimized by this route. We explain molecular orientation in the glass in terms of the orientation present near the surface of the corresponding liquid. The highly oriented glasses formed here also exhibit high density and improved thermal stability. These features will likely further enhance the performance of organic electronics devices. Physical vapor deposition is commonly used to prepare organic glasses that serve as the active layers in light-emitting diodes, photovoltaics, and other devices. Recent work has shown that orienting the molecules in such organic semiconductors can significantly enhance device performance. We apply a high-throughput characterization scheme to investigate the effect of the substrate temperature (Tsubstrate) on glasses of three organic molecules used as semiconductors. The optical and material properties are evaluated with spectroscopic ellipsometry. We find that molecular orientation in these glasses is continuously tunable and controlled by Tsubstrate/Tg, where Tg is the glass transition temperature. All three molecules can produce highly anisotropic glasses; the dependence of molecular orientation upon substrate temperature is remarkably similar and nearly independent of molecular length. All three compounds form “stable glasses” with high density and thermal stability, and have properties similar to stable glasses prepared from model glass formers. Simulations reproduce the experimental trends and explain molecular orientation in the deposited glasses in terms of the surface properties of the equilibrium liquid. By showing that organic semiconductors form stable glasses, these results provide an avenue for systematic performance optimization of active layers in organic electronics.

Model vapor-deposited glasses: growth front and composition effects.

A growing body of experimental work indicates that physical vapor deposition provides an effective route for preparation of stable glasses, whose properties correspond in some cases to those expected for glasses that have been aged for thousands of years. In this work, model binary glasses are prepared in a process inspired by physical vapor deposition, in which particles are sequentially added to the free surface of a growing film in molecular dynamics simulations. The resulting glasses are shown to be more stable than those prepared by gradual cooling from the liquid phase. However, it is also shown that the composition of the resulting glass, which is difficult to control in physical vapor deposition simulations of thin films, plays a significant role on the physical characteristics of the material. That composition dependence leads to a re-evaluation of previous results from simulations of thinner films than those considered here, where the equivalent age of the corresponding glasses was overestimated. The simulations presented in this work, which correspond to films that are approximately 38 molecular diameters thick, also enable analysis of the devitrification mechanism by which vapor-deposited glasses transform into the supercooled liquid. Consistent with experiments, it is found that this mechanism consists of a mobility front that propagates from the free interface into the interior of the films. Eliminating surface mobility eliminates this route of transformation into the supercooled liquid.

Thermodynamic consistency in variable-level coarse graining of polymeric liquids.

Numerically optimized reduced descriptions of macromolecular liquids often present thermodynamic inconsistency with atomistic level descriptions even if the total correlation function, i.e. the structure, appears to be in agreement. An analytical expression for the effective potential between a pair of coarse-grained units is derived starting from the first-principles Ornstein-Zernike equation, for a polymer liquid where each chain is represented as a collection of interpenetrating blobs, with a variable number of blobs, n(b), of size N(b). The potential is characterized by a long tail, slowly decaying with characteristic scaling exponent of N(b)(1/4). This general result applies to any coarse-grained model of polymer melts with units larger than the persistence length, highlighting the importance of the long, repulsive, potential tail for the model to correctly predict both structural and thermodynamic properties of the macromolecular liquid.

First-principle approach to rescale the dynamics of simulated coarse-grained macromolecular liquids.

We present a detailed derivation and testing of our approach to rescale the dynamics of mesoscale simulations of coarse-grained polymer melts (I. Y. Lyubimov, J. McCarty, A. Clark, and M. G. Guenza, J. Chem. Phys. 132, 224903 (2010)). Starting from the first-principle Liouville equation and applying the Mori-Zwanzig projection operator technique, we derive the generalized Langevin equations (GLEs) for the coarse-grained representations of the liquid. The chosen slow variables in the projection operators define the length scale of coarse graining. Each polymer is represented at two levels of coarse graining: monomeric as a bead-and-spring model and molecular as a soft colloid. In the long-time regime where the center-of-mass follows Brownian motion and the internal dynamics is completely relaxed, the two descriptions must be equivalent. By enforcing this formal relation we derive from the GLEs the analytical rescaling factors to be applied to dynamical data in the coarse-grained representation to recover the monomeric description. Change in entropy and change in friction are the two corrections to be accounted for to compensate the effects of coarse graining on the polymer dynamics. The solution of the memory functions in the coarse-grained representations provides the dynamical rescaling of the friction coefficient. The calculation of the internal degrees of freedom provides the correction of the change in entropy due to coarse graining. The resulting rescaling formalism is a function of the coarse-grained model and thermodynamic parameters of the system simulated. The rescaled dynamics obtained from mesoscale simulations of polyethylene, represented as soft-colloidal particles, by applying our rescaling approach shows a good agreement with data of translational diffusion measured experimentally and from simulations. The proposed method is used to predict self-diffusion coefficients of new polyethylene samples.

Molecular modeling of vapor-deposited polymer glasses

We have investigated the properties of vapor-deposited glasses prepared from short polymer chains using molecular dynamics simulations. Vapor-deposited polymer glasses are found to have higher density and higher kinetic stability than ordinary glasses prepared by gradual cooling of the corresponding equilibrium liquid. In contrast to results for binary Lennard-Jones glasses, the deposition rate is found to play an important role in the stability of polymer vapor-deposited glasses. Glasses deposited at the slowest deposition rate and at the optimal substrate temperature are found to correspond to the ordinary glasses that one could hypothetically prepare by cooling the liquid at rates that are 4-5 orders of magnitude slower than those accessible in the current simulations. For intermediate-length polymer chains, the resulting vapor-deposited glasses are found to be highly anisotropic. For short chains, however, the glasses are isotropic, showing that structural anisotropy is not a necessary condition for formation of stable glasses by physical vapor deposition.

Effective Soft-Core Potentials and Mesoscopic Simulations of Binary Polymer Mixtures

Mesoscopic molecular dynamics simulations are used to determine the large-scale structure of several binary polymer mixtures of various chemical architecture, concentration, and thermodynamic conditions. By implementing an analytical formalism, which is based on the solution to the Ornstein−Zernike equation, each polymer chain is mapped onto the level of a single soft colloid. From the appropriate closure relation, the effective, soft-core potential between coarse-grained units is obtained and used as input to our mesoscale simulations. The potential derived in this manner is analytical and explicitly parameter dependent, making it general and transferable to numerous systems of interest. From computer simulations performed under various thermodynamic conditions the structure of the polymer mixture, through pair correlation functions, is determined over the entire miscible region of the phase diagram. In the athermal regime mesoscale simulations exhibit quantitative agreement with united atom simulations....

Age and structure of a model vapour-deposited glass

Glass films prepared by a process of physical vapour deposition have been shown to have thermodynamic and kinetic stability comparable to those of ordinary glasses aged for thousands of years. A central question in the study of vapour-deposited glasses, particularly in light of new knowledge regarding anisotropy in these materials, is whether the ultra-stable glassy films formed by vapour deposition are ever equivalent to those obtained by liquid cooling. Here we present a computational study of vapour deposition for a two-dimensional glass forming liquid using a methodology, which closely mimics experiment. We find that for the model considered here, structures that arise in vapour-deposited materials are statistically identical to those observed in ordinary glasses, provided the two are compared at the same inherent structure energy. We also find that newly deposited hot molecules produce cascades of hot particles that propagate far into the film, possibly influencing the relaxation of the material.

Theoretical reconstruction of realistic dynamics of highly coarse-grained cis-1,4-polybutadiene melts

The theory to reconstruct the atomistic-level chain diffusion from the accelerated dynamics that is measured in mesoscale simulations of the coarse-grained system, is applied here to the dynamics of cis-1,4-polybutadiene melts where each chain is described as a soft interacting colloidal particle. The rescaling formalism accounts for the corrections in the dynamics due to the change in entropy and the change in friction that are a consequence of the coarse-graining procedure. By including these two corrections the dynamics is rescaled to reproduce the realistic dynamics of the system described at the atomistic level. The rescaled diffusion coefficient obtained from mesoscale simulations of coarse-grained cis-1,4-polybutadiene melts shows good agreement with data from united atom simulations performed by Tsolou et al. [Macromolecules 38, 1478 (2005)]. The derived monomer friction coefficient is used as an input to the theory for cooperative dynamics that describes the internal dynamics of a polymer moving in a transient regions of slow cooperative motion in a liquid of macromolecules. Theoretically predicted time correlation functions show good agreement with simulations in the whole range of length and time scales in which data are available.

Influence of Vapor Deposition on Structural and Charge Transport Properties of Ethylbenzene Films

Organic glass films formed by physical vapor deposition exhibit enhanced stability relative to those formed by conventional liquid cooling and aging techniques. Recently, experimental and computational evidence has emerged indicating that the average molecular orientation can be tuned by controlling the substrate temperature at which these “stable glasses” are grown. In this work, we present a comprehensive all-atom simulation study of ethylbenzene, a canonical stable-glass former, using a computational film formation procedure that closely mimics the vapor deposition process. Atomistic studies of experimentally formed vapor-deposited glasses have not been performed before, and this study therefore begins by verifying that the model and method utilized here reproduces key structural features observed experimentally. Having established agreement between several simulated and experimental macroscopic observables, simulations are used to examine the substrate temperature dependence of molecular orientation. The results indicate that ethylbenzene glasses are anisotropic, depending upon substrate temperature, and that this dependence can be understood from the orientation present at the surface of the equilibrium liquid. By treating ethylbenzene as a simple model for molecular semiconducting materials, a quantum-chemical analysis is then used to show that the vapor-deposited glasses exhibit decreased energetic disorder and increased magnitude of the mean-squared transfer integral relative to isotropic, liquid-cooled films, an effect that is attributed to the anisotropic ordering of the molecular film. These results suggest a novel structure–function simulation strategy capable of tuning the electronic properties of organic semiconducting glasses prior to experimental deposition, which could have considerable potential for organic electronic materials design.