Liquid-Crystal Ordering and Microphase Separation in the Lamellar Phase of Rod-Coil-Rod Triblock Copolymers. Molecular Theory and Computer Simulations.

A molecular model of the orientationally ordered lamellar phase exhibited by asymmetric rod-coil-rod triblock copolymers has been developed using the density-functional approach and generalizing the molecular-statistical theory of rod-coil diblock copolymers. An approximate expression for the free energy of the lamellar phase has been obtained in terms of the direct correlation functions of the system, the Flory-Huggins parameter and the Maier-Saupe orientational interaction potential between rods. A detailed derivation of several rod-rod and rod-coil density-density correlation functions required to evaluate the free energy is presented. The orientational and translational order parameters of rod and coil segments depending on the temperature and triblock asymmetry have been calculated numerically by direct minimization of the free energy. Different structure and ordering of the lamellar phase at high and low values of the triblock asymmetry is revealed and analyzed in detail. Asymmetric rod-coil-rod triblock copolymers have been simulated using the method of dissipative particle dynamics in the broad range of the Flory-Huggins parameter and for several values of the triblock asymmetry. It has been found that the lamellar phase appears to be the most stable one at strong segregation. The density distribution of the coil segments and the segments of the two different rods have been determined for different values of the segregation strength. The simulations confirm the existence of a weakly ordered lamellar phase predicted by the density-functional theory, in which the short rods separate from the long ones and are characterized by weak positional ordering.

Solvent Vapor Annealing of a Diblock Copolymer Thin Film with a Nonselective and a Selective Solvent: Importance of Pathway for the Morphological Changes.

Diblock copolymer thin films of polystyrene-block-poly(dimethyl siloxane) (PS-b-PDMS) featuring PDMS cylinders in a PS matrix are investigated during solvent vapor annealing with mixtures of n-heptane (which is strongly selective for PDMS) and toluene (which is close to nonselective for both blocks). Swelling in the vapor of one of the pure solvents and exchanging it stepwise by the vapor of the other solvent is compared to swelling in a given binary solvent vapor mixture for a prolonged time. The resulting structural changes, such as ordering of the cylinders on a hexagonal lattice and their transition into lamellae, are followed using in situ, real-time grazing-incidence small-angle X-ray scattering (GISAXS). In three runs, the sequence of solvent vapor swelling and vapor exchange is varied. Compiling the resulting morphologies in a diagram of states in dependence on the solvent content in the film and in the minority nanodomains allows insight into the role of the glass transition of the PS matrix for the ordering processes and their time scales. Based on these findings, a protocol is suggested to efficiently obtain an order-to-order transition from the cylindrical morphology with random domain orientation to the oriented lamellar state.

Ordering of anisotropic nanoparticles in diblock copolymer lamellae: Simulations with dissipative particle dynamics and a molecular theory.

Local distribution and orientation of anisotropic nanoparticles in microphase-separated symmetric diblock copolymers has been simulated using dissipative particle dynamics and analyzed with a molecular theory. It has been demonstrated that nanoparticles are characterized by a non-trivial orientational ordering in the lamellar phase due to their anisotropic interactions with isotropic monomer units. In the simulations, the maximum concentration and degree of ordering are attained for non-selective nanorods near the domain boundary. In this case, the nanorods have a certain tendency to align parallel to the interface in the boundary region and perpendicular to it inside the domains. Similar orientation ordering of nanoparticles located at the lamellar interface is predicted by the molecular theory which takes into account that the nanoparticles interact with monomer units via both isotropic and anisotropic potentials. Computer simulations enable one to study the effects of the nanorod concentration, length, stiffness, and selectivity of their interactions with the copolymer components on the phase stability and orientational order of nanoparticles. If the volume fraction of the nanorods is lower than 0.1, they have no effect on the copolymer transition from the disordered state into a lamellar microstructure. Increasing nanorod concentration or nanorod length results in clustering of the nanorods and eventually leads to a macrophase separation, whereas the copolymer preserves its lamellar morphology. Segregated nanorods of length close to the width of the diblock copolymer domains are stacked side by side into smectic layers that fill the domain space. Thus, spontaneous organization and orientation of nanorods leads to a spatial modulation of anisotropic composite properties which may be important for various applications.

Vertical vs Lateral Macrophase Separation in Thin Films of Block Copolymer Mixtures: Computer Simulations and GISAXS Experiments.

Mixtures of two diblock copolymers of very different lengths may feature both macro- and microphase separation; however, not much is known about the mechanisms of separation in diblock copolymer thin films. In the present work, we study thin films of mixtures of two compositionally symmetric block copolymers, both in the one-phase and in the two-phase state, combining coarse-grained molecular simulations (dissipative particle dynamics, DPD) with scattering experiments (grazing-incidence small-angle X-ray scattering, GISAXS). We reveal that the film thickness and selective adsorption of different blocks to the substrate control the distribution of macrophases within the film as well as the orientation of the lamellae therein. In thick films, the mixtures separate in the vertical direction into three layers: Two layers being rich in short copolymers are formed near the film interfaces, whereas a layer being rich in long copolymers is located in the film core. The lamellar orientation in the layers rich in short copolymers is dictated by the surface selectivity, and this orientation only weakly affects the vertical orientation of lamellae in the film core. This provides the opportunity to control the domain orientation in the copolymer films by mixing block copolymers with low-molecular additives instead of relying on a more complicated chemical modification of the substrate. In thinner films, a lateral phase separation appears.

Effect of Cross-Linking on the Structure and Growth of Polymer Films Prepared by Interfacial Polymerization.

Interfacial polymerization of tri- and bifunctional monomers (A3B2 polymerization) is investigated by dissipative particle dynamics to reveal an effect of cross-linking on the reaction kinetics and structure of the growing polymer film. Regardless of the comonomer reactivity and miscibility, the kinetics in an initially bilayer melt passes from the reaction to diffusion control. Within the crossover period, branched macromolecules undergo gelation, which drastically changes the scenario of the polymerization process. Comparison with the previously studied linear interfacial polymerization (Berezkin, A. V.; Kudryavtsev, Y. V. Linear Interfacial Polymerization: Theory and Simulations with Dissipative Particle Dynamics J. Chem. Phys. 2014, 141, 194906) shows similar conversion rates but very different product characteristics. Cross-linked polymer films are markedly heterogeneous in density, their average polymerization degree grows with the comonomer miscibility, and end groups are mostly trapped deeply in the film core. Products of linear interfacial polymerization demonstrate opposite trends as they are spontaneously homogenized by a convective flow of macromolecules expelled from the reactive zone to the film periphery, which we call the reactive extrusion effect and which is hampered in branched polymerization. Influence of the comonomer architecture on the polymer film characteristics could be used in various practical applications of interfacial polymerization, such as fabrication of membranes, micro- and nanocapsules and 3D printing.

Linear interfacial polymerization: theory and simulations with dissipative particle dynamics.

Step-growth alternating interfacial polymerization between two miscible or immiscible monomer melts is investigated theoretically and by dissipative particle dynamics simulations. In both cases the kinetics for an initially bilayer system passes from the reaction to diffusion control. The polymer composed of immiscible monomers precipitates at the interface forming a film of nearly uniform density. It is demonstrated that the reaction proceeds in a narrow zone, which expands much slower than the whole film, so that newly formed polymer is extruded from the reaction zone. This concept of "reactive extrusion" is used to analytically predict the degree of polymerization and distribution of all components (monomers, polymer, and end groups) within the film in close agreement with the simulations. Increasing the comonomer incompatibility leads to thinner and more uniform films with the higher average degree of polymerization. The final product is considerably more polydisperse than expected for the homogeneous step-growth polymerization. The results extend the previous theoretical reports on interfacial polymerization and provide new insights into the internal film structure and polymer characteristics, which are important for membrane preparation, microencapsulation, and 3D printing technologies. A systematic way of mapping the simulation data onto laboratory scales is discussed.

Hybrid approach combining dissipative particle dynamics and finite-difference diffusion model: simulation of reactive polymer coupling and interfacial polymerization.

A novel hybrid approach combining dissipative particle dynamics (DPD) and finite difference (FD) solution of partial differential equations is proposed to simulate complex reaction-diffusion phenomena in heterogeneous systems. DPD is used for the detailed molecular modeling of mass transfer, chemical reactions, and phase separation near the liquid∕liquid interface, while FD approach is applied to describe the large-scale diffusion of reactants outside the reaction zone. A smooth, self-consistent procedure of matching the solute concentration is performed in the buffer region between the DPD and FD domains. The new model is tested on a simple model system admitting an analytical solution for the diffusion controlled regime and then applied to simulate practically important heterogeneous processes of (i) reactive coupling between immiscible end-functionalized polymers and (ii) interfacial polymerization of two monomers dissolved in immiscible solvents. The results obtained due to extending the space and time scales accessible to modeling provide new insights into the kinetics and mechanism of those processes and demonstrate high robustness and accuracy of the novel technique.

Simulation of heterogeneous end-coupling reactions in polydisperse polymer blends.

The influence of polydispersity on the interfacial kinetics of end-coupling and microstructure formation in the melt of immiscible polymers was studied using dissipative particle dynamics simulations. The irreversible reaction started at a flat interface between two layers, each of which contained polymer chains of two different lengths with functionalized or unreactive end groups. As in the case of fully functionalized monodisperse reactants [A. V. Berezkin and Y. V. Kudryavtsev, Macromolecules 44, 112 (2011)], four kinetic regimes were observed: linear (mean field coupling at the initial interface), saturation (decreasing the reaction rate due to the copolymer brush formation or reactant depletion near the interface), autocatalytic (loss of the initial interface stability and formation of a lamellar microstructure), and terminal (microstructure ripening under diffusion control). The interfacial instability is caused by overcrowding the interface with the reaction product, and it can be kinetically suppressed by increasing chain length of the reactants. Main effects of polydispersity are as follows: (i) the overall end-coupling rate is dominated by the shortest reactive chains; (ii) the copolymer concentration at the interface causing its instability can be not the same as in the lamellas formed afterwards; (iii) mean length of the copolymer product considerably changes with conversion passing through a minimum when a microstructure is just formed.

Template copolymerization near a patterned surface: computer simulation.

We perform a Monte Carlo simulation of irreversible template copolymerization near a chemically heterogeneous surface with a regular distribution of discrete adsorption sites that selectively adsorb from solution one of the two polymerizing monomers and the corresponding chain segments. In the polymerization model, the chain propagation process is simulated by adding individual monomers to the end of growing macroradical. We focus in this paper on the influence of polymerization rate, adsorption energy, and the distance between adsorption sites on the chain conformation and the primary sequence of the resulting two-letter (AB) copolymers and, specifically, on the coupling between polymerization and adsorption. The conditions for the realization of conformation-dependent copolymerization are formulated. For this regime, we observe the formation of a quasiregular copolymer with two types of alternating sections. One of them contains randomly distributed A and B segments. The second one consists mainly of strongly adsorbed A segments. It is found that the average length of the random sections is proportional to the distance between the nearest neighbor adsorption sites. The average length of the A-rich sections is determined by the "adsorption capacity" of adsorption site. By varying the strength of the effective monomer-substrate interaction and the distribution of adsorption sites on the substrate, the copolymers with different surface-induced primary sequences can be designed and synthesized in a controlled fashion.