Simulation of the gyroid phase in off-lattice models of pure diblock copolymer melts

Emergence of layered nanoscale mesh networks through intrinsic molecular confinement self-assembly

Block copolymer self-assembly is a powerful tool for two-dimensional nanofabrication; however, the extension of this self-assembly concept to complex three-dimensional network structures is limited. Here we report a simple method to experimentally generate three-dimensional layered mesh morphologies through intrinsic molecular confinement self-assembly. We designed triblock bottlebrush polymers with two Janus domains: one perpendicular and one parallel to the polymer backbone. The former enforces a lamellar superstructure that intrinsically confines the intralayer self-assembly of the latter, giving rise to a mesh-like monoclinic (54°) M¹⁵ network substructure with excellent long-range order, as well as a tetragonal (90°) T¹³¹ mesh. Numerical simulations show that the spatial constraints exerted on the polymer backbone drive the assembly of M¹⁵ and yield T¹³¹ in the strong segregation regime. This work demonstrates that intrinsic molecular confinement is a viable path to bottom-up assembly of new geometrical phases of soft matter, extending the capabilities of block copolymer nanofabrication.

Effect of the counterion size on microphase separation in charged-neutral diblock copolymers

In this work, the question of the influence of the counterion size on the self-assembly in melts of diblock copolymers with one charged block was studied using coarse-grained molecular dynamics simulations. It was assumed that the blocks were fully compatible, i.e., the Flory-Huggins parameter χ between them was equal to 0. Due to the presence of correlation attraction (electrostatic cohesion) between the charged species, the systems with all types of counterions underwent transitions to ordered states, forming various morphologies, including lamellae, perforated lamellae, and hexagonally packed cylinders. Phase diagrams were constructed by varying the chain composition f_c and locating the order-disorder transition positions in terms of the electrostatic strength parameter λ (dimensionless Bjerrum length). Despite having a rather large ion size mismatch, the systems with smaller counterions demonstrated an even better tendency to form microphase separated states than the systems with larger ones. It was found that the differences between the phase diagrams of the systems with different counterions can be roughly rationalized by using coordinates (volume fraction of the charged block φ_c-modified interaction parameter λ*). The latter parameter assumes that the electrostatic energy is simply inversely proportional to the characteristic distance between the ions of different signs. Such an approach appeared to be rather effective and allowed the diagrams obtained for different counterion sizes to almost coincide. The results of this work suggest that the counterion size can be used as a tool to control the system morphology as well as the effective incompatibility between the blocks.

Experimental and Computational Evaluation of Self-Assembled Morphologies in Diblock Janus Bottlebrush Copolymers

Diblock Janus-type "A-branch-B" bottlebrush copolymers (di-JBBCPs) consist of a backbone with alternating A and B side chains, in contrast to the side chain arrangement of conventional bottlebrush copolymers. As a result, A and B blocks of di-JBBCPs can microphase-separate perpendicular to the backbone, which is located at the interface between the two blocks. A reparametrized dissipative particle dynamics (DPD) model is used to theoretically investigate the self-assembly of di-JBBCPs and to compare with the experimental results of a range of polystyrene-branch-polydimethylsiloxane di-JBBCPs. The experimentally formed cylinder, gyroid, and lamellar morphologies showed good correspondence with the model phase diagram, and the effect of changing volume fraction and backbone length is revealed. The DPD model predicts a bulk-stable perforated lamella morphology together with two unconventional spherical phases, the Frank-Kasper A15 spheres and the hexagonally close-packed spheres, indicating the diversity of morphologies available from complex BCP molecular architectures.

Obtaining and Characterizing Stable Bicontinuous Cubic Morphologies and Their Nanochannels in Lyotropic Liquid Crystal Membranes

Amphiphilic monomers in polar solvents can self-assemble into lyotropic liquid crystal (LLC) bicontinuous cubic structures under the right composition and temperature conditions. After cross-linking, the resulting polymer membranes with three-dimensional (3D) continuous uniform channels are excellent candidates for filtration applications. Designing such membranes with the desired physical and chemical properties requires molecular-level understanding of the structure, which can be obtained through molecular modeling. However, building molecular models of bicontinuous cubic structures is challenging due to their narrow regime of stability and the difficulty of self-assembly of large unit cells in molecular simulations. We developed a protocol for building stable bicontinuous cubic unit cells involving both parameterization and assembly of the components. We validate the theoretical structure against experimental results for one such LLC monomer and provide insight into the structure missing in experimental data, as well as demonstrate the qualitative nature of water and solute transport through these membranes.

Finding the bulk periodicity of lamellar and cylindrical structures using the pressure tensor

It has been known for nearly thirty years that in all molecular simulations of periodic ordered morphologies (such as those formed by block copolymers), when the periodic boundary conditions of the simulation box do not match the bulk periodicity L₀ of the morphology, they change the structure and even the stability of the morphologies obtained in the simulations. Few studies, however, have focused on finding L₀ in simulations. Taking dissipative particle dynamics simulations of asymmetric diblock copolymer melts as an example, we found a simple way of using the pressure tensor, which can be readily calculated in molecular simulations in the continuum, to find L₀ of lamellae and (regular-hexagonally packed) cylinders regardless of their orientation in (cuboid) simulation boxes. Variation of the diagonal and off-diagonal elements of the pressure tensor with the orientation of lamellae and cylinders in the box is explained by coordinate system rotation and confirmed by our simulation results. We also showed that the pressure tensor cannot be used to find L₀ of 3D periodic ordered structures with a cubic symmetry such as the double gyroid.

Stabilizing a Double Gyroid Network Phase with 2 nm Feature Size by Blending of Lamellar and Cylindrical Forming Block Oligomers

Molecular dynamics simulations are used to study binary blends of an AB-type diblock and an AB₂-type miktoarm triblock amphiphiles (also known as high-χ block oligomers) consisting of sugar-based (A) and hydrocarbon (B) blocks. In their pure form, the AB diblock and AB₂ triblock amphiphiles self-assemble into ordered lamellar (LAM) and cylindrical (CYL) structures, respectively. At intermediate compositions, however, the AB₂-rich blend (0.2 ≤ x _AB ≤ 0.4) forms a double gyroid (DG) network, whereas perforated lamellae (PL) are observed in the AB-rich blend (0.5 ≤ x _AB ≤ 0.8). All of the ordered mesophases present domain pitches under 3 nm, with 1 nm feature sizes for the polar domains. Structural analyses reveal that the nonuniform interfacial curvatures of DG and PL structures are supported by local composition variations of the LAM- and CYL-forming amphiphiles. Self-consistent mean field theory calculations for blends of related AB and AB₂ block polymers also show the DG network at intermediate compositions, when A is the minority block, but PL is not stable. This work provides molecular-level insights into how blending of shape-filling molecular architectures enables network phase formation with extremely small feature sizes over a wide composition range.

Periodicity and global order parameter of hexagonally packed cylinders in a periodic box

In all molecular simulations of periodic ordered morphologies (such as those formed by block copolymers), the periodic boundary conditions (PBCs) of the simulation box usually do not match the bulk periodicity L₀ of the morphology, thus changing the structure and even the stability of the morphologies obtained in the simulations. To address this problem for hexagonally packed cylinders, we first proposed a general method of calculating the periodicity of such cylinders in a cuboid simulation box with the PBCs applied in all directions, which further allows one to enumerate all possible orientations and periodicities of such cylinders within an estimated range that can fit into a cuboid box of given lengths. We then showed how to choose the lengths of a cuboid box such that regular-hexagonally packed (RHP) cylinders with given intercylinder distance and orientation can fit into the box. Next, taking as an example the dissipative particle dynamics (DPD) simulations of a cylinder-forming diblock copolymer melt, we showed that L₀ of RHP cylinders oriented along the body diagonal of a cubic box is found when all the off-diagonal elements of the pressure tensor vanish. Finally, based on our general method of calculating the periodicity of hexagonally packed cylinders, we designed a global order parameter for such cylinders, which takes into account their ordering only for the orientations that can fit into the simulation box. Using again the DPD simulations, we showed that our global order parameter can be used to quantify the formation of hexagonally packed cylinders in each collected configuration and to monitor their orientation (thus periodicity) during the simulation run.

Microphase separation in helix-coil block copolymer melts: computer simulation

By means of molecular dynamics simulation, the process of the microphase separation in the melts of diblock helix-coil copolymers comprising a flexible and a helical block was studied. The resulting microstructures were examined, and the spatial distribution of the blocks and molecular packing were investigated. The phase diagram was built in terms of the fraction of the helical block and the incompatibility parameter of the blocks. The comparison of the diagrams for helix-coil and the classic coil-coil copolymer blends was carried out. It was shown that the total region where the ordering into distinctive microstructures takes place is similar for both diagrams. But for the helix-coil copolymers the area of the cylinders splits into the region of those with circular and elliptical cross-sections; the bicontinuous phase area is much wider; in the lamellar phases, the helical blocks were oriented precisely perpendicular to the lamellar interface, forming a cohesive interlocked structure of densely packed helices.

Phase Behavior of Gradient Copolymer Melts with Different Gradient Strengths Revealed by Mesoscale Simulations

Design and preparation of functional nanomaterials with specific properties requires precise control over their microscopic structure. A prototypical example is the self-assembly of diblock copolymers, which generate highly ordered structures controlled by three parameters: the chemical incompatibility between blocks, block size ratio and chain length. Recent advances in polymer synthesis have allowed for the preparation of gradient copolymers with controlled sequence chemistry, thus providing additional parameters to tailor their assembly. These are polydisperse monomer sequence, block size distribution and gradient strength. Here, we employ dissipative particle dynamics to describe the self-assembly of gradient copolymer melts with strong, intermediate, and weak gradient strength and compare their phase behavior to that of corresponding diblock copolymers. Gradient melts behave similarly when copolymers with a strong gradient are considered. Decreasing the gradient strength leads to the widening of the gyroid phase window, at the expense of cylindrical domains, and a remarkable extension of the lamellar phase. Finally, we show that weak gradient strength enhances chain packing in gyroid structures much more than in lamellar and cylindrical morphologies. Importantly, this work also provides a link between gradient copolymers morphology and parameters such as chemical incompatibility, chain length and monomer sequence as support for the rational design of these nanomaterials.

Correlation between morphology and anisotropic transport properties of diblock copolymers melts

Molecular simulations of coarse-grained diblock copolymers (DBP) were conducted to study the effect of segregation strength and morphology on transport properties. It was found that in the strong segregation limit (i.e., high χN, where χ is the Flory-Huggins parameter and N is the degree of polymerization), the presence of the DBP interfaces imposes topological constraints similar to those of entanglements as manifested in the rheological signature of the polymer (i.e., a plateau modulus). Furthermore, compared to the behavior of isotropic melts, the crossover from Rouse to reptation scaling of the self-diffusion coefficient (D) parallel to the DBP interface takes place at a smaller N, an effect that depends on temperature and is more pronounced in the Lamellae morphology than in the hexagonal cylinder morphology. Additionally, it is shown that for an entangled melt (i.e., N ≫ Ne where Ne is the entanglement length) block retraction is instrumental for chains to diffuse parallel to the interface of lamellar layers. Lastly, it is found that the anisotropic viscosity of different morphologies is mostly affected by the orientation of the chains relative to the shear flow direction, exhibiting reduced values when chains align in the neutral or flow directions.

Diffusion in Lamellae, Cylinders, and Double Gyroid Block Copolymer Nanostructures

We study transport of penetrants through nanoscale morphologies motivated by common block copolymer morphologies, using confined random walks and coarse-grained simulations. Diffusion through randomly oriented grains is 1/3 for cylinder and 2/3 for lamellar morphologies versus an unconstrained (homopolymer) system, as previously understood. Diffusion in the double gyroid structure depends on the volume fraction and is 0.47-0.55 through the minority phase at 30-50 vol % and 0.73-0.80 through the majority at 50-70 vol %. Thus, among randomly oriented standard minority phase structures with no grain boundary effects, lamellae is preferable for transport.

Phase behavior and orientational ordering in block copolymers doped with anisotropic nanoparticles

A molecular field theory and coarse-grained computer simulations with dissipative particle dynamics have been used to study the spontaneous orientational ordering of anisotropic nanoparticles in the lamellar and hexagonal phases of diblock copolymers and the effect of nanoparticles on the phase behavior of these systems. Both the molecular theory and computer simulations indicate that strongly anisotropic nanoparticles are ordered orientationally mainly in the boundary region between the domains and the nematic order parameter possesses opposite signs in adjacent domains. The orientational order is induced by the boundary and by the interaction between nanoparticles and the monomer units in different domains. In simulations, sufficiently long and strongly selective nanoparticles are ordered also inside the domains. The nematic order parameter and local concentration profiles of nanoparticles have been calculated numerically using the model of a nanoparticle with two interaction centers and also determined using the results of computer simulations. A number of phase diagrams have been obtained which illustrate the effect of nanoparticle selectivity and molar fraction of the stability ranges of various phases. Different morphologies have been identified by analyzing the static structure factor and a phase diagram has been constructed in coordinates' nanoparticle concentration-copolymer composition. Orientational ordering of even a small fraction of nanoparticles may result in a significant increase of the dielectric anisotropy of a polymer nanocomposite, which is important for various applications.

Why Is Gyroid More Difficult to Nucleate from Disordered Liquids than Lamellar and Hexagonal Mesophases?

Block copolymers, surfactants, and biomolecules form lamellar, hexagonal, and gyroid mesophases. Across these systems, the nucleation of lamellar from the disordered liquid is the easiest and the nucleation of gyroid the most challenging. This poses the question of what are the factors that determine the rates of nucleation of the mesophases and whether they are controlled by the complexity of the structures or the thermodynamics of nucleation. Here, we use molecular simulations to investigate the nucleation and thermodynamics of lamellar, hexagonal, and gyroid in a binary mixture of particles that produces the same mesophases as those of surfactants and block copolymers. We demonstrate that a combination of averaged bond-order parameters q̅₂ and q̅₈ identifies and distinguishes the three mesophases. We use these parameters to track the microscopic process of nucleation of each mesophase and investigate the existence of heterogeneous nucleation (cross-nucleation) between mesophases. We estimate the surface tensions of the liquid/mesophase interfaces from nucleation rates using classical nucleation theory and find that they are comparable for the three mesophases with values that are about a third of those expected for liquid-crystal interfaces. The driving forces for nucleation, on the other hand, are quite different and increase in the order gyroid < hexagonal < lamellar at any temperature. We find that the nucleation rates of the mesophases follow the order of their driving forces. We conclude that the difficulty to nucleate the gyroid originates in its lower temperature of melting and extremely low entropy of melting compared to those of the hexagonal and lamellar mesophases.

Molecular dynamics simulation of thermotropic bolaamphiphiles with a swallow-tail lateral chain: formation of cubic network phases

T-shaped bolaamphiphiles (TBA) with a swallow-tail lateral chain have been found to provide a fertile platform to produce complex liquid crystalline phases that are accessible through changes of temperature and lateral chain length and design. In this work, we use molecular simulations of a simple coarse-grained model to map out the phase behavior of this type of molecules. This model is based on the premise that the crucial details of the fluid structure stem from close range repulsions and the strong directional forces typical of hydrogen bonds. Our simulations confirm that TBAs exhibit a rich phase behavior upon increasing the length of their lateral chain. The simulations detect a double gyroid phase and an axial-bundle columnar phase which bear some structural resemblance to those found in the experiment. In addition, simulations predict two cocontinuous phases with 3D-periodicity: the "single" diamond and the "single" plumber's nightmare phase. Our analysis of energetic and entropic contributions to the free energy of phases formed by TBA with either swallow-tail or linear side-chains suggest that the 3D-periodic network phases formed by the former are stabilized by the large conformation entropy of the side-chains.

Universality between Experiment and Simulation of a Diblock Copolymer Melt

The equivalent behavior among analogous block copolymer systems involving chemically distinct molecules or mathematically different models has long hinted at an underlying universality, but only recently has it been rigorously demonstrated by matching results from different simulations. The profound implication of universality is that simple coarse-grained models can be calibrated so as to provide quantitatively accurate predictions to experiment. Here, we provide the first compelling demonstration of this by simulating a polyisoprene-polylactide diblock copolymer melt using a previously calibrated lattice model. The simulation successfully predicts the peak in the disordered-state structure function, the position of the order-disorder transition, and the latent heat of the transition in excellent quantitative agreement with experiment. This could mark a new era of precision in the field of block copolymer research.

Phase diagrams of block copolymer melts by dissipative particle dynamics simulations

Phase diagrams for monodisperse and polydisperse diblock copolymer melts and a random multiblock copolymer melt are constructed using dissipative particle dynamics simulations. A thorough visual analysis and calculation of the static structure factor in several hundreds of points at each of the diagrams prove the ability of mesoscopic molecular dynamics to predict the phase behavior of polymer systems as effectively as the self-consistent field-theory and Monte Carlo simulations do. It is demonstrated that the order-disorder transition (ODT) curve for monodisperse diblocks can be precisely located by a spike in the dependence of the mean square pressure fluctuation on χN, where χ is the Flory-Huggins parameter and N is the chain length. For two other copolymer types, the continuous ODTs are observed. Large polydispersity of both blocks obeying the Flory distribution in length does not shift the ODT curve but considerably narrows the domains of the cylindrical and lamellar phases partially replacing them with the wormlike micelle and perforated lamellar phases, respectively. Instead of the pure 3d-bicontinuous phase in monodisperse diblocks, which could be identified as the gyroid, a coexistence of the 3d phase and cylindrical micelles is detected in polydisperse diblocks. The lamellar domain spacing D in monodisperse diblocks follows the strong-segregation theory prediction, D∕N(1∕2) ~ (χN)(1∕6), whereas in polydisperse diblocks it is almost independent of χN at χN < 100. Completely random multiblock copolymers cannot form ordered microstructures other than lamellas at any composition.

Lattice model of oligonucleotide hybridization in solution. I. Model and thermodynamics

A coarse-grained lattice model of DNA oligonucleotides is proposed to investigate the general mechanisms by which single-stranded oligonucleotides hybridize to their complementary strands in solution. The model, based on a high-coordination cubic lattice, is simple enough to allow the direct simulation of DNA solutions, yet capturing how the fundamental thermodynamic processes are microscopically encoded in the nucleobase sequences. Physically relevant interactions are considered explicitly, such as interchain excluded volume, anisotropic base-pairing and base-stacking, and single-stranded bending rigidity. The model is studied in detail by a specially adapted Monte Carlo simulation method, based on parallel tempering and biased trials, which is designed to overcome the entropic and enthalpic barriers associated with the sampling of hybridization events of multiple single-stranded chains in solution. This methodology addresses both the configurational complexity of bringing together two complementary strands in a favorable orientation (entropic barrier) and the energetic penalty of breaking apart multiple associated bases in a double-stranded state (enthalpic barrier). For strands with sequences restricted to nonstaggering association and homogeneous pairing and stacking energies, base-pairing is found to dominate the hybridization over the translational and conformational entropy. For strands with sequence-dependent pairing corresponding to that of DNA, the complex dependence of the model's thermal stability on concentration, sequence, and degree of complementarity is shown to be qualitatively and quantitatively consistent both with experiment and with the predictions of statistical mechanical models.

Design rule for colloidal crystals of DNA-functionalized particles

We report a Monte Carlo simulation study of the phase behavior of colloids coated with long, flexible DNA chains. We find that an important change occurs in the phase diagram when the number of DNAs per colloid is decreased below a critical value. In this case, the triple point disappears and the condensed phase that coexists with the vapor is always liquid. Our simulations thus explain why, in the dilute solutions typically used in experiments, colloids coated with a small number of DNA strands cannot crystallize. We understand this behavior in terms of the discrete nature of DNA binding.

A systematic coarse-graining strategy for semi-dilute copolymer solutions: from monomers to micelles

A systematic coarse-graining procedure is proposed for the description and simulation of AB diblock copolymers in selective solvents. Each block is represented by a small number, n(A) or n(B), of effective segments or blobs, containing a large number of microscopic monomers. n(A) and n(B) are unequivocally determined by imposing that blobs do not, on average, overlap, even if complete copolymer coils interpenetrate (semi-dilute regime). Ultra-soft effective interactions between blobs are determined by a rigorous inversion procedure in the low concentration limit. The methodology is applied to an athermal copolymer model where A blocks are ideal (theta solvent), B blocks self-avoiding (good solvent), while A and B blocks are mutually avoiding. The model leads to aggregation into polydisperse spherical micelles beyond a critical micellar concentration determined by Monte Carlo simulations for several size ratios f of the two blocks. The simulations also provide accurate estimates of the osmotic pressure and of the free energy of the copolymer solutions over a wide range of concentrations. The mean micellar aggregation numbers are found to be significantly lower than those predicted by an earlier, minimal two-blob representation (Capone et al 2009 J. Phys. Chem. B 113 3629).

Simulations of the gyroid phase in diblock copolymers with the Gaussian disphere model

Pure melts of asymmetric diblock copolymers are studied by means of the off-lattice Gaussian disphere model with Monte-Carlo kinetics. In this model, a diblock copolymer chain is mapped onto two soft repulsive spheres with fluctuating radii of gyration and distance between centers of mass of the spheres. Microscopic input quantities of the model such as the combined probability distribution for the radii of gyration and the distance between the spheres as well as conditional monomer number densities assigned to each block were derived in the previous work of F. Eurich and P. Maass [J. Chem. Phys. 114, 7655 (2001)] within an underlying Gaussian chain model. The polymerization degree of the whole chain as well as those of the individual blocks are freely tunable parameters thus enabling a precise determination of the regions of stability of various phases. The model neglects entanglement effects which are irrelevant for the formation of ordered structures in diblock copolymers and which would otherwise unnecessarily increase the equilibration time of the system. The gyroid phase was reproduced in between the cylindrical and lamellar phases in systems with box sizes being commensurate with the size of the unit cell of the gyroid morphology. The region of stability of the gyroid phase was studied in detail and found to be consistent with the prediction of the mean-field theory. Packing frustration was observed in the form of increased radii of gyration of both blocks of the chains located close to the gyroid nodes.