Computational Design of High-χ Block Oligomers for Accessing 1 nm

Apr 16, 2018 - Molecular dynamics simulations are used to design a series of high-χ block oligomers (HCBOs) that can self-assemble into a variety of ...
0 downloads 4 Views 12MB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

Computational Design of High-# Block Oligomers for Accessing 1-Nanometer Domains Qile P. Chen, Leonel Barreda, Luis E. Oquendo, Marc A. Hillmyer, Timothy P. Lodge, and J. Ilja Siepmann ACS Nano, Just Accepted Manuscript • Publication Date (Web): 16 Apr 2018 Downloaded from http://pubs.acs.org on April 16, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Computational Design of High-χ Block Oligomers for Accessing 1-Nanometer Domains Qile P. Chen,†,‡,¶ Leonel Barreda,‡ Luis E. Oquendo,‡ Marc A. Hillmyer,‡ Timothy P. Lodge,∗,†,‡ and J. Ilja Siepmann∗,†,‡,¶ Department of Chemical Engineering and Materials Science, University of Minnesota, 421 Washington Avenue SE, Minneapolis, Minnesota 55455-0132, United States, Department of Chemistry, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States, and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States E-mail: [email protected]; [email protected]

∗ To

whom correspondence should be addressed

† CEMS ‡ CHEM ¶ CTC

1

ACS Paragon Plus Environment

χMF (Tsim = 400 K or Texpt = 300 K)

Graphical TOC Entry mono-oligomers HO

y-1 OH

x-2

di-oligomers (x-3)/2 HO

y-1 (x-3)/2 OH

10 mono-oligomers di-oligomers Expt mean-field 1

0.1

0

2

5 10 15 20 Domain period d [nm]

Page 3 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Abstract Molecular dynamics simulations are used to design a series of high-χ block oligomers (HCBOs) that can self-assemble into a variety of mesophases with domain sizes as small as 1 nm. The exploration of these oligomers with various chain lengths, volume fractions, and chain architectures at multiple temperatures reveals the presence of ordered lamellae, perforated lamellae, and hexagonally-packed cylinders. The achieved periods are as small as 3.0 nm and 2.1 nm for lamellae and cylinders, respectively, which correspond to polar domains of approximately 1 nm. Interestingly, the phase behavior of these oligomers is distinct from that of either solvent-free surfactants or block polymers in detail. The simulations reveal that the behavior of these HCBOs is a product of interplay between both “surfactant factors” (head group interactions, chain flexibility, and interfacial curvature) and “block polymer factors” (χ, chain length N, and volume fraction f ). This insight promotes the understanding of molecular features pivotal for the mesophase formation at the sub-5 nm length scale, which facilitates the design of HCBOs tailored towards particular desired morphologies.

Keywords high-χ block polymer, solvent-free surfactant, block polymer phase behavior, molecular simulations, computational design

In the microelectronics industry, the growing need for microprocessors with faster speed and storage devices with larger storage densities demands the miniaturization of the associated features below the current practical limit of about 10 nanometers. 1 Further feature size decreases using traditional lithographic processes faces significant challenges going forward. 2 Nanolithography utilizing directed self-assembly of block polymers is considered an attractive and cost-effective approach to achieve the desired single digit nanometer length scales. 1,3 For block polymers, the morphology and the separation between self-assembled domains can be tailored by tuning the block polymer composition ( fi ), the overall degree of polymerization (N), and the interaction parameter 3

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

between the blocks (χ). Pioneering research has been devoted to exploring the materials design space and identifying, synthesizing, and implementing new high-χ block polymers (oligomers) that can be readily used for such applications. 4–8 In addition, a variety of processing techniques including thin-film coating, 9 defect removal, 10–12 and domain alignment 13–15 have been developed to make block polymers more practically viable. However, typical domain periods d (d is defined as the domain spacing for lamellae or the 100 interplanar distance for hexagonally packed cylinders) for self-assembled block polymers (d ≈ 15–100 nm) are still higher than the desired periods (d ≤ 5 nm) for nanolithography 1,3,16 and other applications. 17–20 Since d increases with N, while χN must be large enough to access the microphase separation region, the key to reducing d is to find “high χ–low N” block polymers (oligomers) which contain more incompatible blocks. 1 There are numerous examples in the literature aimed at practical approaches to such highχ–low N systems in attempts to drive down d. 21–44 For example, Kennemur and co-workers 35 have demonstrated that domain sizes of sub-5 nm can be achieved without the incorporation of silicon, fluorine, or metal species (i.e., “all-organic” block polymers), the presence of which is undesirable for convenience of post-processing and sustainability. 1 In addition, the use of highlyhydroxylated blocks has recently been proved a promising method. In particular, Otsuka et al. utilized maltoheptaose to fabricate block polymers with d as small as 7.7 nm. 39 Nowak et al. designed a sugar-olefin conjugate that achieved the nanostructure formation with d = 6.6 ± 0.7 nm. 42 Similarly, Jeong et al. utilized a solid-state transition to self-assemble block oligomers that contain multiple hydroxyl groups into 5.6-nm-spaced lamellae. 43 Further reduction in domain separation calls for predictive design of high-χ block oligomers (HCBOs). This, in turn, creates a demand for more comprehensive understanding of the phase behavior at this length scale, which lies at the boundary between self-assembled block polymers and solvent-free surfactants. 1 The microphase separation of block polymers is driven by the interplay between segment incompatibility and the enforced block connectivity. Thus, the associated complex phase behavior can be largely attributed to the chain flexibility that can induce preferred curvature at the interface, while the precise molecular geometry (e.g., head group area, tail group

4

ACS Paragon Plus Environment

Page 4 of 34

Page 5 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

volume) plays a less important role. 5 On the other hand, low molar mass surfactants typically adopt only layered smectic or nematic phases in their solvent-free liquid crystalline state and not other complex morphologies. 45,46 This liquid crystalline-type phase behavior stems largely from anisotropic interactions such as hydrogen bonding and π-π interactions, where the self-assembly behavior is heavily influenced by molecular geometry. 47 While the general principles for both block polymer and liquid crystalline surfactant self-assembly are well-established, it is still a matter of discussion as to whether these theories are applicable to sub-5 nm HCBOs. 5,30 Therefore, the complication resulting from this interplay motivates a more rigorous understanding of the phase behavior at the sub-5 nm length scale, and molecular simulation offers a promising approach to address this challenge. Earlier works in simulating block polymers (in fact, polymers in general) often utilize coarsegrained force fields for which the “resolution” of the model is relatively low (e.g., one monomer is modeled as a spherical interaction site and its chemical structure is obscured). These pioneering studies provide valuable insight into the phase behavior of block polymers, 48–54 but coarse-grained force fields are insufficient for the predictive design of new block oligomers, because the detailed molecular features play a more important role at shorter length scales. Therefore, we use atomisticlevel force fields with higher resolution to model HCBOs. This allows for more refined descriptions of molecular-level interactions, and higher transferability of the force field parameters compared to coarse-grained molecular models. As a consequence, a direct relationship between molecular structure and phase behavior can be established without any extraneous fitting parameters. This provides a convenient route to identify HCBOs that can self-assemble into structures with sub-5 nm domains, and to generate further insight into critical molecular features governing self-assembly at the block polymer/surfactant boundary. We present a molecular dynamics (MD) simulation aimed at designing HCBOs that can selfassemble into various complex ordered mesophases with domain periods ranging from 2 to 4 nm. These amphiphilic molecules consist of hydrophilic polyol head groups and alkanes as hydrophobic tail(s). The molecular structures of these oligomers are shown in Figure 1a. It is noteworthy

5

ACS Paragon Plus Environment

a)

mono-x-y

di-x-y

tri-x-y

(x-4)/3

(x-3)/2 HO

b)

y-1 OH

x-2

Narm

HO

y-1 (x-3)/2 OH

tri-13-4 DIS-MIC tri-16-4

DIS-MIC

di-13-4

CYL-2.1 di-21-4

DIS di-17-4

mono-12-4 L-3.0 mono-16-4

tri-22-4

PL-3.5 mono-20-4

(x-4)/3 HO

y-1 (x-4)/3 OH

DIS-MIC tri-28-4

DIS-MIC

CYL-2.2 di-29-4

CYL-2.6

PL-3.9 mono-28-4

fnon-polar Figure 1: (a) Molecular structures of the HCBOs investigated. x and y indicate the number of hydroxyl groups and the number of CHn groups in the alkyl tail(s), respectively. (b) Simulation snapshots of representative morphologies observed at Tsim = 493 K. The caption labels the morphology followed by the domain period d (in nm). “L” = lamellae, “PL” = perforated lamellae, “CYL” = hexagonally-packed cylinders, and “DIS(-MIC)” = disordered(-micellar) phase. Atoms in the polar block are rendered in red (including the carbon atoms) and atoms in the non-polar block are rendered in cyan. For the two PL snapshots, one lamellae layer is depicted to show the holes. For both tri-16-4 and tri-28-4 snapshots, the second images show different micelles with different colors.

6

DIS

Page 7 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

that simulated oligomers are a mixture of different stereoisomers since the tacticity is not preserved in the generation of initial molecular structures. The design space (chain length, volume fraction, chain architecture, and temperature) is explored in silico via molecular simulations. Morphologies of these oligomers are identified through both visualization and quantitative analysis. We also probe key factors that govern the self-assembly behavior such as chain orientation and hydrogen bonding from the analysis of simulation trajectories to understand the interplay between “surfactant factors” (e.g., head group interactions, chain flexibility, interfacial curvature) and “block polymer factors” (χ, N, and f ). We provide predictions for HCBOs that exhibit exciting features for nanotechnological applications and elucidate pivotal factors that contribute to the mesophase formation at the nanometer length scale.

Results and Discussion Molecular models and experimental validation. The oligomers studied in this work are represented by the Transferable Potentials for Phase Equilibria united-atom (TraPPE–UA) force field. 55–58 Details about the molecular models and simulation techniques are provided in the Simulation Methodology section, as well as in the Supporting Information. The key to accurately predicting block oligomer phase behavior is to quantitatively reproduce the χ parameter between the two blocks, i.e., the interactions between the polyol head group and the alkyl tail(s). Among the experimentally measurable quantities, the solubility of an alkane in a polyol is very sensitive to χ and, hence, it serves as the target property for the molecular model validation. In particular, the solubilities of propane in 1,2-ethanediol at various temperatures and pressures are calculated from Gibbs ensemble Monte Carlo simulations 59,60 and compared to the experimental data, 61 due to the unavailability of solubility data for longer alkane and polyols (likely because they are strongly immiscible). The simulation results shown in Figure S1 suggest that the solubilities are overestimated due to the lack of polarizability in the TraPPE–UA model, but the solubilities at p ≥ 10 MPa (where both phases are highly incompressible and “liquid-like”) can be reproduced within 10% er-

7

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ror if the simulation temperatures are scaled by a factor of 0.7–0.8 (see Supporting Information for complete simulation details and results). Therefore, relatively high temperatures are used in the subsequent simulations to account for the underestimation of the unlike interactions. Characterization of the complex mesophases observed. The morphologies of the solvent-free oligomers are examined via MD simulations. After an initial exploration of the design space, we focus on oligomers that are composed of four hydroxyl groups (y = 4) and up to three alkyl tails at Tsim = 493, 523, and 573 K. The molecular structures and representative simulation snapshots of the mesophase at Tsim = 493 K are shown in Figure 1. The morphologies and volume fractions of all the oligomers can be found in Table S5. A variety of mesophases are observed, including lamellae, perforated lamellae, hexagonally-packed cylinders, disordered micellar phases, and disordered phases. As illustrated in Figure 1b, one-arm (mono-) oligomers with a short hydrophobic block (x = 12) can phase-separate into lamellar domains. Similarly, mono-oligomers having a longer alkyl block (x = 14–26) can form perforated lamellae. When the length of the alkyl tail further increases to x = 28, the long-range ordering vanishes and the morphology transitions into the disordered phase. For two-arm (di-) oligomers, hexagonally-packed cylinders are observed in the simulations when x = 17–29. When the number of arms Narm further increases to three (tri-), disordered spherical micelles emerge from the simulations for all four chain lengths investigated (x = 13–28). The presence of these mesophases is confirmed by the peak patterns in the structure factors shown in Figures 2, S4–S10. In addition, the domain period d can be calculated from the peak positions in the structure factor through the relationship d = 2π/q∗ , where q∗ is the characteristic peak position determined from the peak pattern. For lamellae, d is equal to the center-to-center spacing of the lamellae, L (d = L = Lpolar + Lnonpolar , in which Lpolar and Lnonpolar are the widths of each lamella, respectively), and for hexagonally-packed cylinders, d is the 100 interplanar distance. √ The center-to-center spacing of cylinders L can be calculated via d = 3L/2. From the structure factor data, d for (perforated) lamellae structures at Tsim = 493 K varies from 3.0 to 3.9 nm when 12 ≤ x ≤ 20, and the d for hexagonally-packed cylinders ranges from 2.1 to 2.6 nm for 17 ≤

8

ACS Paragon Plus Environment

Page 8 of 34

Page 9 of 34

d [nm] 54 3

2

1 mono-12-4 LAM q/q* = 1, 2, 3, ...

mono-16-4 PL q/q* = 1, 2, 3, 4, ...

S(q) [relative unit]

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

mono-20-4 PL q/q* = 1, 2, 3, 4, ...

di-17-4 HEX-CYL 1/2 1/2 q/q* = 1, 3 , 4 , ...

tri-16-4 DIS-MIC

0.0

0.2

0.4

-1 0.6 q [Å ]

0.8

1.0

Figure 2: Structure factors for representative systems at Tsim = 493 K in logarithmic scale. The data are shifted vertically for clarity. The shadow indicates uncertainties of the data. The arrow symbols indicate the peak patterns corresponding to ordered structures. x ≤ 29. The domain periods d of 3.0 nm for lamellae and 2.1 nm for cylinders are the smallest separations reported in the experimental literature for “all-organic” block polymers/oligomers. 1,41 The corresponding domain sizes (defined as the lamellar width Lpolar for lamellae, and the cylinder diameter Dcyl for hexagonally-packed cylinders) estimated from d and volume fraction f (see Supporting Information for calculation details) are as small as 1 nm. The phase behavior of these HCBOs is unique as it stands in between that of solvent-free nonionic surfactants and block polymers. On the one hand, unlike diblock polymers, mono-oligomers can self-assemble into lamellae and perforated lamellae, but without the formation of cylindrical or spherical structures (but note that long-chain primary alcohols, such as 1-octanol, are known

9

ACS Paragon Plus Environment

to form cylindrical micelles 62 ). On the other hand, compared to commonly-used solvent-free nonionic surfactants, the transition temperatures from the liquid crystalline state to the isotropic liquid state of these oligomers are much higher, which implies that χ is much greater for these oligomers. 5,46,62 The unique aspects of block oligomer phase behavior suggests that molar mass and chain architecture both play a pivotal role; the remainder of this paper is focused on further understanding of their phase behavior. a)

b)

mono-16-4 (PL)

180 150 3.0

di-17-4 (CYL)

120

Probability density

2.5

1.5 1.0

60

r2

2.0

tri-16-4 (SPH)

φ [deg] 90

Φ

r1

30

0

r2

r1

mono-16-4 (PL) di-17-4 (CYL) tri-16-4 (SPH)

0.5 0.0 -1.0

-0.5

0.0 cos(φ)

0.5

1.0

Figure 3: (a) Simulation snapshots of selected alkyl tails (cyan) conformations at the polar-apolar interface (black dashed lines). (b) Angular distributions of the neighboring alkyl arm end-to-end vectors for three sets of oligomers at Tsim = 493 K. The schematic drawings depict the anti-parallel, perpendicular, and parallel orientation of the neighboring arms. The effects of chain architecture on interfacial curvature and hydrogen bonding. To elucidate why high-χ oligomers with multiple arms can self-assemble into cylindrical and spherical structures, we first focus on how Narm shapes the phase behavior through its influence on interfacial curvature. As illustrated in the simulation snapshots in Figure 3a, the polar-apolar interface is more curved as Narm increases, likely resulting from a bulkier alkyl tail at the interface. The shape of the interface can be quantitatively described by the angle between the end-to-end vectors of the neighboring alkyl arms (the alkyl arm is defined as one linear alkyl chain which stems from

10

Page 11 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

the linker carbon at the interface, e.g., each di-oligomer has two alkyl arms) as indicated in Figure 3b. For three sets of oligomers that have approximately the same alkyl chain length (x = 16 for mono- and tri-, and x = 17 for di-oligomers) but different values of Narm , the angular distributions are clearly distinct. For mono-oligomers, a maximum at φ = 180◦ suggests that the parallel orientation is dominant, which is in agreement with a flat polar-apolar interface. For di-oligomers, the peak height at φ = 180◦ decreases dramatically and the angle between two neighboring arms are more randomly distributed. This packing results in a more curved polar-apolar interface and thus the formation of cylinders, since the random orientation corresponds to neighboring arms around the cross section of the cylinder and the parallel orientation can likely be attributed to the arms along the longest axis of the cylinder. For tri-oligomers, the parallel orientation is further suppressed and the distribution becomes nearly uniform with a slightly higher bump at around φ = 100◦ . The flat orientational distribution corresponds to a more spherical geometry as there is no preferred orientation for chain packing. Therefore, the increase in Narm can effectively engender a greater curvature at the polar-apolar interface and thus a more curved morphology is favored when Narm > 1 (cylinders for Narm = 2 and spheres for Narm = 3). This result is consistent with the trend predicted from the molecular packing parameter, 63 as well as the mean-field prediction of miktoarm block polymer phase diagram. 39,64–67 In addition to the interfacial curvature engendered by the bulkiness of multiple arms, the second crucial factor that facilitates or frustrates the formation of ordered mesophases is the hydrogen bonding capability (equivalently, head group interactions in surfactant theory, or χ in block polymer theory). Figure 4a shows the number of hydrogen bonds per molecule Nhb at Tsim = 493 K. To ensure that the observed trends are not caused by the geometric hydrogen-bond criteria being too strict or too loose, we performed additional hydrogen-bonding analysis using a looser criterion and similar trends can be found in the data presented in Figure S16 and Table S6. Stronger hydrogen bonding capability implies weaker head-tail interactions relative to head-head interactions, and equivalently, a higher effective χ parameter between two blocks. However, it is noteworthy that the hydrogen bonding is the result of complexation events, and it is essential to use more

11

ACS Paragon Plus Environment

ACS Nano

a) 2.6

Nhb [per mlcl]

2.5

2.4

mono-oligomers di-oligomers tri-oligomers

2.3

2.2 10

15

20 x

b) [Nhb(mono-) - Nhb(tri-)] / Nhb(mono-) [%]

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 34

25

1

25

30

2 3 4

20

15

10

5 x = 16 x = 28 0 1

2 3 Hydroxyl group label

4

Figure 4: (a) The number of hydrogen bonds per molecule Nhb for three types of oligomers at Tsim = 493 K. The filled and open symbols correspond to oligomers in the ordered and disordered phases, respectively. Snapshots of perforated lamellae for mono-14-4, mono-20-4, and mono-26-4 are shown. (b) The relative decrease in Nhb from mono-oligomers to tri-oligomers as a function of hydroxyl group position for oligomers with two different alkyl tail lengths at Tsim = 493 K. complicated models to quantitatively capture the effect of hydrogen bonding on the effective χ parameter. 68 One can see from Figure 4 and Figure S16 that Nhb values for mono-oligomers using both of the hydrogen-bond criteria show a continuous decrease as a function of alkyl chain length x. In contrast, Nhb values for di-oligomers show no significant chain-length dependence. As for tri-oligomers, Nhb decreases slightly as x increases. This is likely due to the different morphologies of the mesophases present for mono-, di-, and tri-oligomers. For mono-oligomers, as x increases, the polar lamellae contain more holes as the volume fraction of the polar block decreases (see simulation snapshots of polar lamellae in Figure 4). The increase in interfacial area leads to a larger

12

ACS Paragon Plus Environment

Page 13 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

fraction of the hydroxyl groups being at the interface and, hence, results in a significant decrease in Nhb . In contrast, di-oligomers form cylindrical structures, where the change in interfacial area as a function of x is relatively small (the diameter of the cylindrical polar region changes from 1.0 nm for di-17-4 to 1.2 nm for di-29-4). This small change in diameter results in the insensitivity of Nhb to the chain length. As for tri-oligomers, they form the disordered micellar phase, which leads to a chain-length dependence in between mono- and di-oligomers. Increasing the tail length for the tri-oligomers leads to a decrease in the aggregation number and, hence, an increase in the interfacial area (see Figures S12 and S13). In addition, the fraction of free tri-oligomers that are not involved in hydrogen-bonded aggregates increases with increasing x (see Figure S12 and S13). Together, the two effects lead to a slight decrease in Nhb . However, the effect of chain architecture on Nhb becomes more pronounced when Narm further increases to three. Nhb drops by 12% for tri-16-4 oligomers compared to that of the mono-16-4 oligomers, more than a factor of two greater than the decrease from mono- to di-oligomers. The decrease is likely due to a combination of phase change and steric hindrance imposed by the three bulky arms that suppresses the presence of other hydroxyl groups near the interface. The steric hindrance works against the enthalpic driving force towards phase separation, and the resulting disordered phase further lowers Nhb . This hypothesis is supported by the most significant decrease in Nhb for the hydroxyl group closest to the alkyl tail(s) among all four hydroxyl groups, as shown in Figure 4b. For oligomers with two different alkyl tail lengths, the hydroxyl group closer to the alkyl tail(s) suffers more from the loss of hydrogen bonds. The first hydroxyl group in the trioligomers has more than 20% fewer hydrogen bonds compared to that in the mono-oligomers with the same chain length, whereas the decrease in Nhb for the last two hydroxyl groups stays below 6%. Therefore, the reduction in hydrogen bonding for tri-oligomers results in their inability to form any ordered mesophase. In conclusion, hydrogen bonding is one of the main driving forces towards phase separation, and tuning the chain architecture is one effective approach to control hydrogen bonding and thus engineer the mesophase formation. The effect of chain length on molecular conformation and packing. The key to interpreting

13

ACS Paragon Plus Environment

a)

c)

mono-12-4

0.61

20

d)

0.6

Ree [Å]

0.0

mono-16-4 0.79

493 K 573 K

ν ≈ 0.6

16

x

20

28

32

ν ≈ 0.8

Weight fraction

0.73

12

Lnonpolar

0.2

0.88

b)

Lpolar

0.4

0.61 0.69

15

10

0.8

0.65

0.77

1.0

0.8 0.6 0.4 0.2 0.0

mono-20-4

0.8 0.6 0.4

polar apolar

0.2 0.0 -20

-10 0 10 Distance to bilayer center [Å]

Figure 5: The effect of alkyl chain length on molecular conformation and packing in the (perforated) lamellae mesophase. (a) End-to-end chain distance of the alkyl arm (Ree ) as a function of x for mono-oligomers, on logarithmic scale. The slopes for the two or three neighboring data points, which correspond to the Flory exponent ν, are given. (b) Schematic drawing of how the chain conformation affects packing, where red and blue lines correspond to polar and apolar blocks, respectively, and dashed lines indicate the interface. (c) Simulation snapshots of representative (perforated) lamellae for three oligomers at Tsim = 493 K. Atoms in the polar block are rendered in red and magenta for two leaflets (chains that have their head groups in the same cluster form one leaflet), respectively, and apolar groups that stem from two leaflets are rendered in cyan and green, respectively. (d) Composition profile along the axis perpendicular to the lamellae plane corresponding to the three oligomers shown in (c). The composition profiles for apolar blocks are shown separated for each leaflet (dashed lines) and for both leaflets combined (solid lines). The cyan-colored region indicates the interdigitation region where apolar weight fractions of both leaflets are above 0.2.

14

20

Page 15 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

the lack of cylindrical and spherical mesophases for mono-oligomers is the understanding of how the alkyl chain length alters the molecular conformation and thus the stability of mesophases. Figure 5a illustrates the scaling behavior of the end-to-end distance of the alkyl arm as a function of its chain length Ree ∝ xν , where ν is the Flory exponent. Note that the end-to-end carbon distance of the polar block is always about 3.3 Å (a predominantly gauche conformation, see Table S7), and thus, it plays a less important role in affecting the bulk morphology. In the disordered state at Tsim = 573 K, ν is found to decrease as x increases. It ranges from 0.8 for shorter chains to 0.6 for longer ones, which is greater than the scaling exponent of 0.5 for a polymer in its pure melt. This indicates that the oligomers are more rigid compared to their polymer counterpart, which aligns well with results from our earlier studies on olefin oligomers. 69 At Tsim = 493 K, the Flory exponent decreases continuously as a function of x, from 0.9 to 0.6, which reveals that particularly the shorter chains become more rigid at higher temperatures. It is noteworthy that the Flory exponents of the short oligomers (x ≤ 20) are greater than those predicted by the mean-field theory for block polymers in the ordered state (ν = 0.5 from the weak segregation theory 70,71 and ν = 0.67 from the strong segregation theory 71 ). The greater persistence of the chains implies a transition from flexible coils to rigid rods as the alkyl tail chain length decreases. For flexible block polymers, perforated lamellae are considered thermodynamically unstable. 72,73 However, the presence of perforated lamellae has been reported by both theoretical 74 and experimental 75 studies for rod-coil block polymers, which is consistent with the mono-oligomer phase behavior. The enhanced stability of perforated lamellae is schematically explained in Figure 5b. When the chains are sufficiently long, the scaling behavior (ν ≈ 0.6) approaches that of the polymer melt (ν ≈ 0.5). In this case, the chains can be modeled as cones that can effectively pack into hexagonal structures. In contrast, when the oligomer chains are more rigid (ν ≈ 0.8), the rod-like chains lack the ability to form hexagonally-packed structures and instead prefers layered structures. 74 At this range of chain length, the oligomer phase behavior resembles that of a liquid crystal, where cylindrical and spherical structures can result in packing defects and thus are energetically not favored. 76

15

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 34

The effect of chain conformation on the layered structure packing is further confirmed by the simulation snapshots of the bilayer structure and its composition profile depicted in Figures 5c and 5d, respectively. From the simulation snapshots, alkyl chains tend to stretch out from the interface to the center of bilayer without much interdigitation between the chains around the center of bilayer. This observation is supported by the composition profile along the axis perpendicular to the lamellae plane. As x increases from 12 to 20, the thickness of the alkyl layer Lnonpolar grows from 2.1 to 3.0 nm (Table S8). However, the change in the thickness of the alkyl interdigitation region (shown in cyan in the figure) is less pronounced, from 0.35 nm for x = 12 to 0.45 nm for x = 20. This implies that the chains attempt to avoid interdigitation, which would potentially perturb their extended conformation and parallel alignment. This observation confirms the rod-like chain conformation and its effect on the phase behavior. It is intriguing whether the patterning of such small domain will result in larger line edge roughness (LER), or even defects in the lamellae plane due to the Plateau-Rayleigh instability. 77 LER (3σ ) is computed for the bulk self-assembled mono-12-4 oligomer at Tsim = 493 K and a value of 2.0 ± 0.3 Å is obtained (see Supporting Information for calculation details). The LER is around 20% of the domain width, which is sufficient for most nanolithography applications. 78–80 Figure 6a illustrates the phase map that includes all the oligomers investigated in this study. The mean-field χ parameter, χMF , is obtained from fitting our simulation results to the order-disorder phase boundary from self-consistent field theory 64,65,73,81,82 (see Supporting Information for more details), which results in the following:

χMF (Tsim ) =

2150 − 1.1 Tsim

(1)

χMF (Tsim ) =

3440 − 2.8 Tsim

(2)

for mono-oligomers, and

for di-oligomers. For both cases, a typical reference volume of 118 Å3 is used. There is insufficient data to fit χMF for tri-oligomers since all the tri-oligomers covered in this study are in the disordered

16

ACS Paragon Plus Environment

Page 17 of 34

a) 40 30

mono-DIS mono-LAM mono-PL mean-field boundary

χMFN

20 10 40 30

di-DIS di-CYL mean-field boundary

20 10 0.5

b)

χMF (Tsim = 400 K or Texpt = 300 K)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

0.6

0.7 0.8 fnon-polar

0.9

1.0

10 mono-oligomers di-oligomers Expt mean-field 1

0.1

0

5

10 15 Domain period d [nm]

20

Figure 6: (a) The phase map for the oligomers investigated shown in symbols. χMF , the meanfield χ, is fitted to the order-disorder phase boundary from the self-consistent mean-field theory 64,65,73,81,82 represented by the dashed lines. (b) The χMF as a function of domain period d for mean-field boundary (χMF > 13.7/(d/nm)2 , the red dashed line), available “all-organic” materials (black symbols), 21,29,32–38 and oligomers from this work (triangles). χMF values are from extrapolation to Texpt = 300 K or Tsim = 400 K, and the simulation domain period values are from the Tsim = 493 K data.

17

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

phase. The comparison of χ and the achieved domain period d for the mean-field boundary (assuming a typical statistical segment length of b = 0.7 nm for convenience of comparison), a few selected experimentally available block polymers, 21,29,32–38 and the oligomers from this work are shown in Figure 6b. χMF can be as large as 4.3 and 5.8 at Tsim = 400 K (Texpt ≈ 300 K) for mono-oligomers and di-oligomers, respectively. These oligomers have χMF values that are one to two orders of magnitude greater than those of other “all-organic” block oligomers. For example, the χMF of poly(styrene)-block-poly(methyl methacrylate), one of the most widely studied diblock polymer for self-assembly, is around 0.03 at T = 300 K, 83 while the largest χMF value for “all-organic” diblock polymers reported experimentally is less than 0.5 at T = 300 K. 32,35 The greater χ for these oligomers allows the decrease in N and thus enables the formation of 1-nm domains. One should note, however, that use of χMF is a convenient but flawed approach. Nevertheless, presenting some of the data and discussion in terms of the χ parameter allows for comparison of our results to those from the standard diblock copolymer model applicable to blocks consisting of at least a few Kuhn segments and not involving strong directional interactions. There are two major factors that are not included in the fitting of χMF . First, complexation interactions such as hydrogen bonding are not considered in the Flory–Huggins theory formulation. 68 Consequently, χMF is a function of y due to the cooperative nature of hydrogen bonding. For example, we performed simulations for mono-6-2 and mono-9-3 at Tsim = 493 K. Both are found to form homogeneous disordered phases, which leads to the conclusion that χMF (y = 4) > χMF (y = 3) (Table S5). Second, the effects of short chain length and semiflexibility are not taken into account in the fitting of mono- and di-oligomers data using the mean-field AB diblock copolymer model and the meanfield AB2 star polymer model, respectively. 84,85 A recent computational work by Mao et al. 85 has discussed the effects of finite chain length and semiflexibility on the order-disorder transition point χNODT for symmetric diblock polymers. According to their results, 85 χNODT differs by around 20% from the mean-field theory for symmetric diblock copolymers that have similar chain length and semiflexibility as HCBOs. Despite that the quantitative determination of χ is hampered by the

18

ACS Paragon Plus Environment

Page 18 of 34

Page 19 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

lack of χNODT ( f ) theory data for oligomeric semiflexible chains, these valuable data demonstrate that the χMF values from the semi-quantitative fitting in our manuscript are only for comparison purposes and not accurate reflections of only the interactions between the two blocks in standard block polymer models. Despite the loss of predictive power from the mean-field theory, the deviations from the Flory– Huggins theory framework can be more of a blessing than a curse. Potential applications of HCBOs with directional hydrogen-bonding interactions 86 and fairly rigid backbone can arise from their ability to form domain morphologies, such as anisotropic lamellae phases 74,87 and complex spherical packing phases, 88,89 that are difficult to access for flexible diblock polymers. Therefore, further simulation, experimental, and theoretical studies are needed to elucidate the phase behavior and additional applications of these HCBOs.

Conclusions In this study, we have designed a series of high-χ block oligomers with various chain lengths, volume fractions, and chain architectures. From MD simulations, these oligomers can self-assemble into a variety of mesophases including lamellae, perforated lamellae, and hexagonally-packed cylinders. The smallest domain periods in lamellae and cylinders are 3.0 nm and 2.1 nm, respectively, and the smallest polar domain sizes are around 1 nm. These are the smallest reported in the literature so far for “all-organic” block oligomers. At this length scale, the phase behavior stands in between that of block polymer and that of solvent-free surfactant, which arises from the interplay between both “block-polymer factors” (χ, N, and f ) and “surfactant factors” (head group interactions, chain flexibility, and interfacial curvature). The chain architecture can greatly affect the interfacial curvature and hydrogen bonding, which consequently facilitates or frustrates the phase separation and the resulting mesophases. The chain conformation, on the other hand, is more stretched than that of its polymer counterpart, which is believed to be the cause for the unusually high stability of perforated lamellae phase. These insights provide a route to engineer

19

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

similar HCBOs at the sub-5 nm length scale to the desired morphology and domain separation for applications such as nanolithography.

Simulation Methodology Molecular models. The TraPPE–UA force field is used to model these HCBOs. 55–58 The model represents the hydrogen and oxygen atoms of hydroxyl groups as individual interaction sites, and each CHn (n = 0–3) group as one pseudo-atom located on the position of the carbon atom. The nonbonded interactions among interaction sites include pairwise-additive Lennard-Jones (LJ) and Coulomb potentials, which are calculated for all intermolecular pairs and for intramolecular pairs that are separated by four or more bonds. The potential is truncated at a spherical cutoff of 14 Å 90 and is shifted to ensure the continuity of the potential at the cutoff radius. Lorentz–Berthelot combining rules are used to describe the interactions between unlike interaction sites. 91 As implemented in the TraPPE–UA force field, a short range repulsive term between the oxygen and hydrogen atom that are separated by four bonds is used to prevent the unphysical overlap of the two sites (since the proton is not protected by the LJ potential and the 1–4 oxygen-oxygen LJ potential is not included). In the standard TraPPE–UA model, the interaction sites are joined by fixed bond lengths, but bending and dihedral angles are allowed to vary. In this work, the bond length is still fixed for bonds that contain protons to allow the use of larger time step, while harmonic stretching potentials are used for all the other bonds. The equilibrium bond length is taken from the TraPPE–UA force field and the stretching force constant is adopted from the OPLS force field. 92 All the force field parameters are listed in Tables S1 and S2. Additional simulations are performed to validate the molecular model, and details are included in the Supporting Information. Simulation details. MD simulations in the N pT ensemble are performed via the GROMACS 5.0.0 software 93,94 with a time step of 2 fs. System sizes of 500–1000 molecules are used for the simulations. The temperature and pressure are controlled via the Nosé–Hoover thermostat 95,96 with a coupling time constant of 0.4 ps and the Parrinello–Rahman barostat 97 with a 2 ps time constant,

20

ACS Paragon Plus Environment

Page 20 of 34

Page 21 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

respectively. The Particle mesh Ewald method is used for the electrostatics calculation. 98 All bond lengths with hydrogen atoms are kept constant using the P-Lincs algorithm. 99 The initial structures are generated via the Monte Carlo (MC) simulations in the canonical ensemble (Tsim = 3000 K, ρ = 0.7 g/cm3 ) for 500 MC cycles (one MC cycle consisting of N randomly chosen MC steps, with N being the number of molecules in the system) using MCCCS–MN (Monte Carlo for complex chemical systems–Minnesota) software program. 100 The types of MC moves employed are similar to those used in our earlier work on oligomers, namely translation, rotation, and configurational moves. 69 In addition, the k-d tree data structure is used to accelerate the simulation. 101 After the disordered configuration is obtained, GROMACS is used to equilibrate the system at Tsim = 573 K and p = 1 bar. For subsequent lower temperature runs, the system is annealed from equilibrated higher temperature configuration with a rate of 2 K/ns. Each simulation is performed for at least 150 ns, during which each molecule moves at least three times its radius of gyration (Table S4). The first 50 ns of the trajectory is discarded and the remainder is used for the analysis. The uncertainties are calculated from four independent simulations with a block size of 20 ns for each run, and are reported as 95% confidence interval. VMD and Taychon are used to visualize the simulation trajectories and make the snapshot images. 102,103 Some analyses are enabled via the TRAVIS trajectory analyzer. 104 The details for the analysis can be found in the Supporting Information. To avoid the incommensurability effect, 105,106 all of the simulations are performed for each system using an anisotropic orthorhombic simulation box, where the three dimensions of the simulation box are allowed to fluctuate independently to accommodate mesophases that are anisotropic in at least two dimensions (e.g., lamellae, perforated lamellae, and cylinders). In addition, simulations initiated from ordered mesophases are performed to assess the inaccuracy in determining the phase transition due to hysteresis. It turns out that the hysteresis is less than 30 K for the largest system size used in this study (Figure S2). Furthermore, to avoid being trapped in a local minimum, we performed four independent simulations starting from different configurations, and all of them converge to the same morphology with the same domain spacing. Figures S19 and S20 show the individual structure factors obtained from each independent simulation for the mono-12-4

21

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

and di-17-4 systems. The excellent agreement between independent simulations indicates that the simulations reproducibly lead to the same morphology. The average lengths for the orthorhombic box obtained from these independent simulations are provided in Table S9. It should be noted that the relative box lengths are not the same for all simulations; that is, the ordered morphology can be achieved through different arrangements with respect to the periodic simulation cell. As a further test to assess whether the simulation results are free from major artifacts that may be caused by the relatively small system sizes accessible to simulations with atomistically detailed force fields, we performed additional simulations for the mono-12-4 oligomer using a larger system size (N = 1500, in addition to the N = 1000 previously reported). The computed structure factor is included in Figure S18, and the peak positions are a very close match for the two system sizes but, again, the relative box lengths are different.

22

ACS Paragon Plus Environment

Page 22 of 34

Page 23 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Acknowledgement This work was supported by the National Science Foundation through the University of Minnesota MRSEC under Award Number DMR-1420013. The authors acknowledge the Minnesota Supercomputing Institute (MSI) at the University of Minnesota for providing computational resources that contributed to this work.

Supporting Information Available Molecular models and validation, details of trajectory analysis, and numerical data are provided in the Supporting Information. This material is available free of charge via the Internet at http: //pubs.acs.org.

References 1. Sinturel, C.; Bates, F. S.; Hillmyer, M. A. High χ-Low N Block Polymer: How Far can We Go? ACS Macro Lett. 2015, 4, 1044–1050 2. Sanders, D. P. Advances in Patterning Materials for 193 nm Immersion Lithography. Chem. Rev. 2010, 110, 321–360 3. Bates, C. M.; Maher, M. J.; Janes, D. W.; Ellison, C. J.; Willson, C. G. Block Copolymer Lithography. Macromolecules 2013, 47, 2–12 4. Bates, F. S.; Fredrickson, G. H. Block Copolymers–Designer Soft Materials. Phys. Today 1999, 52, 32–38 5. Hillmyer, M. A.; Bates, F. S.; Almdal, K.; Mortensen, K.; Ryan, A. J.; Fairclough, J. P. A. Complex Phase Behavior in Solvent-Free Nonionic Surfactants. Science 1996, 271, 976 6. Kim, S. H.; Misner, M. J.; Xu, T.; Kimura, M.; Russell, T. P. Highly Oriented and Ordered Arrays from Block Copolymers via Solvent Evaporation. Adv. Mater. 2004, 16, 226–231 23

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

7. Jung, Y. S.; Ross, C. A. Orientation-Controlled Self-Assembled Nanolithography Using a Polystyrene–Polydimethylsiloxane Block Copolymer. Nano Lett. 2007, 7, 2046–2050 8. Bang, J.; Kim, S. H.; Drockenmuller, E.; Misner, M. J.; Russell, T. P.; Hawker, C. J. DefectFree Nanoporous Thin Films from ABC Triblock Copolymers. J. Am. Chem. Soc. 2006, 128, 7622–7629 9. Sinturel, C.; Vayer, M.; Morris, M.; Hillmyer, M. A. Solvent Vapor Annealing of Block Polymer Thin Films. Macromolecules 2013, 46, 5399–5415 10. Li, W.; Nealey, P. F.; de Pablo, J. J.; Müller, M. Defect Removal in the Course of Directed Self-Assembly is Facilitated in the Vicinity of the Order-Disorder Transition. Phys. Rev. Lett. 2014, 113, 168301 11. Blachut, G.; Sirard, S. M.; Maher, M. J.; Asano, Y.; Someya, Y.; Lane, A. P.; Durand, W. J.; Bates, C. M.; Dinhobl, A. M.; Gronheid, R. A Hybrid Chemo-/Grapho-Epitaxial Alignment Strategy for Defect Reduction in Sub-10 nm Directed Self-Assembly of Silicon-Containing Block Copolymers. Chem. Mater. 2016, 28, 8951–8961 12. Simão, C.; Khunsin, W.; Kehagias, N.; Salaun, M.; Zelsmann, M.; Morris, M. A.; Torres, C. M. S. Order Quantification of Hexagonal Periodic Arrays Fabricated by In Situ SolventAssisted Nanoimprint Lithography of Block Copolymers. Nanotechnology 2014, 25, 175703 13. Durand, W. J.; Carlson, M. C.; Maher, M. J.; Blachut, G.; Santos, L. J.; Tein, S.; Ganesan, V.; Ellison, C. J.; Willson, C. G. Experimental and Modeling Study of Domain Orientation in Confined Block Copolymer Thin Films. Macromolecules 2015, 49, 308–316 14. Wang, Q.; Yan, Q.; Nealey, P. F.; de Pablo, J. J. Monte Carlo Simulations of Diblock Copolymer Thin Films Confined between Two Homogeneous Surfaces. J. Chem. Phys. 2000, 112, 450–464

24

ACS Paragon Plus Environment

Page 24 of 34

Page 25 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

15. Bates, C. M.; Seshimo, T.; Maher, M. J.; Durand, W. J.; Cushen, J. D.; Dean, L. M.; Blachut, G.; Ellison, C. J.; Willson, C. G. Polarity-Switching Top Coats Enable Orientation of Sub–10-nm Block Copolymer Domains. Science 2012, 338, 775–779 16. Herr, D. J. C. Directed Block Copolymer Self-Assembly for Nanoelectronics Fabrication. J. Mater. Res. 2011, 26, 122–139 17. Werber, J. R.; Osuji, C. O.; Elimelech, M. Materials for Next-Generation Desalination and Water Purification Membranes. Nat. Rev. Mater. 2016, 1, 16018 18. Xu, J.; Howard, M. J.; Mittal, V.; Bates, F. S. Block Copolymer Micelle Toughened Isotactic Polypropylene. Macromolecules 2017, 50, 6421–6432 19. Li, T.; Zhang, J.; Schneiderman, D. K.; Francis, L. F.; Bates, F. S. Toughening Glassy Poly(lactide) with Block Copolymer Micelles. ACS Macro Lett. 2016, 5, 359–364 20. Bates, C. M.; Bates, F. S. 50th Anniversary Perspective: Block Polymers–Pure Potential. Macromolecules 2016, 50, 3–22 21. Xu, T.; Kim, H. C.; DeRouchey, J.; Seney, C.; Levesque, C.; Martin, P.; Stafford, C. M.; Russell, T. P. The Influence of Molecular Weight on Nanoporous Polymer Films. Polymer 2001, 42, 9091–9095 22. Jung, Y. S.; Chang, J. B.; Verploegen, E.; Berggren, K. K.; Ross, C. A. A Path to Ultranarrow Patterns using Self-Assembled Lithography. Nano Lett. 2010, 10, 1000–1005 23. Jeong, J. W.; Park, W. I.; Kim, M.; Ross, C. A.; Jung, Y. S. Highly Tunable Self-Assembled Nanostructures from a Poly(2-vinylpyridine-b-dimethylsiloxane) Block Copolymer. Nano Lett. 2011, 11, 4095–4101 24. Pitet, L. M.; Wuister, S. F.; Peeters, E.; Kramer, E. J.; Hawker, C. J.; Meijer, E. W. Well-Organized Dense Arrays of Nanodomains in Thin Films of Poly(dimethylsiloxane)-bPoly(lactide) Diblock Copolymers. Macromolecules 2013, 46, 8289–8295 25

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

25. Cushen, J. D.; Otsuka, I.; Bates, C. M.; Halila, S.; Fort, S.; Rochas, C.; Easley, J. A.; Rausch, E. L.; Thio, A.; Borsali, R. Oligosaccharide/silicon-containing Block Copolymers with 5 nm Features for Lithographic Applications. ACS Nano 2012, 6, 3424–3433 26. Aissou, K.; Otsuka, I.; Rochas, C.; Fort, S.; Halila, S.; Borsali, R. Nano-Organization of Amylose-b-Polystyrene Block Copolymer Films Doped with Bipyridine. Langmuir 2011, 27, 4098–4103 27. Durand, W. J.; Blachut, G.; Maher, M. J.; Sirard, S.; Tein, S.; Carlson, M. C.; Asano, Y.; Zhou, S. X.; Lane, A. P.; Bates, C. M. Design of High-χ Block Copolymers for Lithography. J. Polym. Sci., Part A: Polym. Chem. 2015, 53, 344–352 28. Park, S.; Liang, X.; Harteneck, B. D.; Pick, T. E.; Hiroshiba, N.; Wu, Y.; Helms, B. A.; Olynick, D. L. Sub-10 nm Nanofabrication via Nanoimprint Directed Self-Assembly of Block Copolymers. ACS Nano 2011, 5, 8523–8531 29. Mayeda, M. K.; Hayat, J.; Epps, T. H.; Lauterbach, J. Metal Oxide Arrays from Block Copolymer Thin Film Templates. J. Mater. Chem. A 2015, 3, 7822–7829 30. Zha, R. H.; de Waal, B. F. M.; Lutz, M.; Teunissen, A. J. P.; Meijer, E. W. End Groups of Functionalized Siloxane Oligomers Direct Block-Copolymeric or Liquid-Crystalline SelfAssembly Behavior. J. Am. Chem. Soc. 2016, 138, 5693–5698 31. Van Genabeek, B.; de Waal, B. F. M.; Gosens, M. M. J.; Pitet, L. M.; Palmans, A. R. A.; Meijer, E. W. Synthesis and Self-Assembly of Discrete Dimethylsiloxane–Lactic Acid Diblock Co-Oligomers: The Dononacontamer and Its Shorter Homologues. J. Am. Chem. Soc. 2016, 138, 4210–4218 32. Zha, W.; Han, C. D.; Lee, D. H.; Han, S. H.; Kim, J. K.; Kang, J. H.; Park, C. Origin of the Difference in Order–Disorder Transition Temperature between Polystyrene-block-Poly(2vinylpyridine) and Polystyrene-block-Poly(4-vinylpyridine) Copolymers. Macromolecules 2007, 40, 2109–2119 26

ACS Paragon Plus Environment

Page 26 of 34

Page 27 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

33. Kennemur, J. G.; Hillmyer, M. A.; Bates, F. S. Synthesis, Thermodynamics, and Dynamics of Poly(4-tert-butylstyrene-b-methyl methacrylate). Macromolecules 2012, 45, 7228–7236 34. Sweat, D. P.; Kim, M.; Larson, S. R.; Choi, J. W.; Choo, Y.; Osuji, C. O.; Gopalan, P. Rational Design of a Block Copolymer with a High Interaction Parameter. Macromolecules 2014, 47, 6687–6696 35. Kennemur, J. G.; Yao, L.; Bates, F. S.; Hillmyer, M. A. Sub-5 nm Domains in Ordered Poly(cyclohexylethylene)-block-Poly(methyl methacrylate) Block Polymers for Lithography. Macromolecules 2014, 47, 1411–1418 36. Cochran, E. W.; Morse, D. C.; Bates, F. S. Design of ABC Triblock Copolymers Near the ODT with the Random Phase Approximation. Macromolecules 2003, 36, 782–792 37. Zalusky, A. S.; Olayo-Valles, R.; Wolf, J. H.; Hillmyer, M. A. Ordered Nanoporous Polymers from Polystyrene–Polylactide Block Copolymers. J. Am. Chem. Soc. 2002, 124, 12761–12773 38. Baruth, A.; Seo, M.; Lin, C. H.; Walster, K.; Shankar, A.; Hillmyer, M. A.; Leighton, C. Optimization of Long-Range Order in Solvent Vapor Annealed Poly(styrene)-block-Poly(lactide) Thin Films for Nanolithography. ACS Appl. Mater. Interfaces 2014, 6, 13770–13781 39. Otsuka, I.; Zhang, Y.; Isono, T.; Rochas, C.; Kakuchi, T.; Satoh, T.; Borsali, R. Sub-10 nm Scale Nanostructures in Self-Organized Linear Di-and Triblock Copolymers and Miktoarm Star Copolymers Consisting of Maltoheptaose and Polystyrene. Macromolecules 2015, 48, 1509–1517 40. Vanderlaan, M. E.; Hillmyer, M. A. “Uncontrolled” Preparation of Disperse Poly(lactide)block-Poly(styrene)-block-Poly(lactide) for Nanopatterning Applications. Macromolecules 2016, 49, 8031–8040 41. Carter, M. C. D.; Jennings, J.; Speetjens, F. W.; Lynn, D. M.; Mahanthappa, M. K. A Reactive

27

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Platform Approach for the Rapid Synthesis and Discovery of High χ/Low N Block Polymers. Macromolecules 2016, 49, 6268–6276 42. Nowak, S. R.; Hwang, W.; Sita, L. R. Dynamic Sub-10-nm Nanostructured Ultrathin Films of Sugar–Polyolefin Conjugates Thermoresponsive at Physiological Temperatures. J. Am. Chem. Soc. 2017, 139, 5281–5284 43. Jeong, G.; Yu, D. M.; Mapas, J. K. D.; Sun, Z.; Rzayev, J.; Russell, T. P. Realizing 5.4 nm Full Pitch Lamellar Microdomains by a Solid-State Transformation. Macromolecules 2017, 50, 7148–7154 44. Feng, X.; Tousley, M. E.; Cowan, M. G.; Wiesenauer, B. R.; Nejati, S.; Choo, Y.; Noble, R. D.; Elimelech, M.; Gin, D. L.; Osuji, C. O. Scalable Fabrication of Polymer Membranes with Vertically Aligned 1 nm Pores by Magnetic Field Directed Self-Assembly. ACS Nano 2014, 8, 11977–11986 45. Chandrasekhar, S. Liquid Crystals; Cambridge University, 1992 46. Tiddy, G. J. T. Surfactant-Water Liquid Crystal Phases. Phys. Rep. 1980, 57, 1–46 47. Ward, M. D.; Horner, M. J. Structure and Order in Soft Matter. CrystEngComm 2004, 6, 401– 407 48. Kikuchi, M.; Binder, K. Microphase Separation in Thin Films of the Symmetric DiblockCopolymer Melt. J. Chem. Phys. 1994, 101, 3367–3377 49. Wang, Q.; Nealey, R. F.; de Pablo, J. J. Monte Carlo Simulations of Asymmetric Diblock Copolymer Thin Films Confined between Two Homogeneous Surfaces. Macromolecules 2001, 34, 3458–3470 50. Larson, R. G. Self-Assembly of Surfactant Liquid Crystalline Phases by Monte Carlo Simulation. J. Chem. Phys. 1989, 91, 2479–2488

28

ACS Paragon Plus Environment

Page 28 of 34

Page 29 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

51. Peters, A. J.; Lawson, R. A.; Nation, B. D.; Ludovice, P. J.; Henderson, C. L. Coarse-Grained Molecular Dynamics Modeling of the Kinetics of Lamellar Block Copolymer Defect Annealing. J. Micro/Nanolithogr., MEMS, MOEMS 2016, 15, 013508–013508 52. Chremos, A.; Nikoubashman, A.; Panagiotopoulos, A. Z. Flory-Huggins Parameter χ, from Binary Mixtures of Lennard-Jones Particles to Block Copolymer Melts. J. Chem. Phys. 2014, 140, 054909 53. Srinivas, G.; Discher, D. E.; Klein, M. L. Self-Assembly and Properties of Diblock Copolymers by Coarse-Grain Molecular Dynamics. Nat. Mater. 2004, 3, 638–644 54. Lee, W. B.; Mezzenga, R.; Fredrickson, G. H. Anomalous Phase Sequences in Lyotropic Liquid Crystals. Phys. Rev. Lett. 2007, 99, 187801 55. Martin, M. G.; Siepmann, J. I. Transferable Potentials for Phase Equilibria. 1. United-Atom Description of n-alkanes. J. Phys. Chem. B 1998, 102, 2569–2577 56. Martin, M. G.; Siepmann, J. I. Novel Configurational-Bias Monte Carlo Method for Branched Molecules. Transferable Potentials for Phase Equilibria. 2. United-Atom Description of Branched Alkanes. J. Phys. Chem. B 1999, 103, 4508–4517 57. Chen, B.; Potoff, J. J.; Siepmann, J. I. Monte Carlo Calculations for Alcohols and Their Mixtures with Alkanes. Transferable Potentials for Phase Equilibria. 5. United-Atom Description of Primary, Secondary, and Tertiary Alcohols. J. Phys. Chem. B 2001, 105, 3093–3104 58. Stubbs, J. M.; Potoff, J. J.; Siepmann, J. I. Transferable Potentials for Phase Equilibria. 6. United-Atom Description for Ethers, Glycols, Ketones, and Aldehydes. J. Phys. Chem. B 2004, 108, 17596–17605 59. Panagiotopoulos, A. Z. Direct Determination of Phase Coexistence Properties of Fluids by Monte Carlo Simulation in a New Ensemble. Mol. Phys. 1987, 61, 813–826

29

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

60. Panagiotopoulos, A. Z.; Quirke, N.; Stapleton, M.; Tildesley, D. J. Phase Equilibria by Simulation in the Gibbs Ensemble: Alternative Derivation, Generalization and Application to Mixture and Membrane Equilibria. Mol. Phys. 1988, 63, 527–545 61. Jou, F. Y.; Otto, F. D.; Mather, A. E. The Solubility of Propane in 1,2-Ethanediol at Elevated Pressures. J. Chem. Thermo. 1993, 25, 37–40 62. Chen, B.; Siepmann, J. I. Microscopic Structure and Solvation in Dry and Wet Octanol. J. Phys. Chem. B 2006, 110, 3555–3563 63. Israelachvili, J. N.; Mitechell, D. J.; Ninham, B. W. Theory of Self-Assembly of Hydrocarbon Amphiphiles into Micelles and Bilayers. J. Chem. Soc. Faraday Trans. 2 1976, 72, 1525–1568 64. Grason, G. M.; Kamien, R. D. Interfaces in Diblocks: A Study of Miktoarm Star Copolymers. Macromolecules 2004, 37, 7371–7380 65. Grason, G. M.; DiDonna, B. A.; Kamien, R. D. Geometric Theory of Diblock Copolymer Phases. Phys. Rev. Lett. 2003, 91, 058304 66. Milner, S. T. Chain Architecture and Asymmetry in Copolymer Microphases. Macromolecules 1994, 27, 2333–2335 67. Minehara, H.; Pitet, L. M.; Kim, S.; Zha, R. H.; Meijer, E. W.; Hawker, C. J. Branched Block Copolymers for Tuning of Morphology and Feature Size in Thin Film Nanolithography. Macromolecules 2016, 49, 2318–2326 68. Dehghan, A.; Shi, A.-C. Modeling Hydrogen Bonding in Diblock Copolymer/homopolymer Blends. Macromolecules 2013, 46, 5796–5805 69. Chen, Q. P.; Chu, J. D.; DeJaco, R. F.; Lodge, T. P.; Siepmann, J. I. Molecular Simulation of Olefin Oligomer Blend Phase Behavior. Macromolecules 2016, 49, 3975–3985 70. Vavasour, J. D.; Whitmore, M. D. Self-Consistent Mean Field Theory of the Microphases of Diblock Copolymers. Macromolecules 1992, 25, 5477–5486 30

ACS Paragon Plus Environment

Page 30 of 34

Page 31 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

71. Bates, F. S.; Fredrickson, G. H. Block Copolymer Thermodynamics: Theory and Experiment. Annu. Rev. Phys. Chem. 1990, 41, 525–557 72. Khandpur, A. K.; Forster, S.; Bates, F. S.; Hamley, I. W.; Ryan, A. J.; Bras, W.; Almdal, K.; Mortensen, K. Polyisoprene-Polystyrene Diblock Copolymer Phase-Diagram Near the OrderDisorder Transition. Macromolecules 1995, 28, 8796–8806 73. Matsen, M. W.; Bates, F. S. Unifying Weak- and Strong-Segregation Block Copolymer Theories. Macromolecules 1996, 29, 1091–1098 74. Pryamitsyn, V.; Ganesan, V. Self-Assembly of Rod–Coil Block Copolymers. J. Chem. Phys. 2004, 120, 5824–5838 75. Tenneti, K. K.; Chen, X.; Li, C. Y.; Tu, Y.; Wan, X.; Zhou, Q. F.; Sics, I.; Hsiao, B. S. Perforated Layer Structures in Liquid Crystalline Rod–Coil Block Copolymers. J. Am. Chem. Soc. 2005, 127, 15481–15490 76. Matsen, M. W.; Barrett, C. Liquid-Crystalline Behavior of Rod-Coil Diblock Copolymers. J. Chem. Phys. 1998, 109, 4108–4118 77. Zhou, Y.; Nicolas, A.; Thomas, K. R.; Steiner, U. Interplay of Electrohydrodynamic Structure Formation and Microphase Alignment in Lamellar Block Copolymers. Soft Matter 2012, 14, 3841–3847 78. Cheng, J. Y.; Rettner, C. T.; Sanders, D. P.; Kim, H.; Hinsberg, W. D. Dense Self-Assembly on Sparse Chemical Patterns: Rectifying and Multiplying Lithographic Patterns using Block Copolymers. Adv. Mater. 2008, 20, 3155–3158 79. Stoykovich, M. P.; Daoulas, K. C.; Müller, M.; Kang, H.; de Pablo, J. J.; Nealey, P. F. Remediation of Line Edge Roughness in Chemical Nanopatterns by the Directed Assembly of Overlying Block Copolymer Films. Macromolecules 2010, 43, 2334–2342

31

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

80. Kim, J. M.; Kim, Y.; Park, W. I.; Hur, Y. H.; Jeong, J. W.; Sim, D. M.; Baek, K. M.; Lee, J. H.; Kim, M.; Jung, Y. S. Eliminating the Trade-Off between the Throughput and Pattern Quality of Sub-15 nm Directed Self-Assembly via Warm Solvent Annealing. Adv. Funct. Mater. 2015, 25, 306–315 81. Matsen, M. W.; Schick, M. Stable and Unstable Phases of a Diblock Copolymer Melt. Phys. Rev. Lett. 1994, 72, 2660 82. Matsen, M. W. Effect of Architecture on the Phase Behavior of AB-Type Block Copolymer Melts. Macromolecules 2012, 45, 2161–2165 83. Russell, T. P.; Hjelm Jr, R. P.; Seeger, P. A. Temperature Dependence of the Interaction Parameter of Polystyrene and Poly(methyl methacrylate). Macromolecules 1990, 23, 890–893 84. Qin, J.; Morse, D. C. Fluctuations in Symmetric Diblock Copolymers: Testing Theories Old and New. Phys. Rev. Lett. 2012, 108, 238301 85. Mao, S.; MacPherson, Q.; Spakowitz, A. J. Fluctuation Effects in Semiflexible Diblock Copolymers. ACS Macro Lett. 2018, 6, 59–64 86. Ladanyi, B. M.; Skaf, M. S. Computer Simulation of Hydrogen-Bonding Liquids. Annu. Rev. Phys. Chem. 1993, 44, 335–368 87. Olsen, B. D.; Segalman, R. A. Self-Assembly of Rod–Coil Block Copolymers. Mater. Sci. Eng., R 2008, 62, 37–66 88. Kim, K.; Schulze, M. W.; Arora, A.; Lewis, R. M.; Hillmyer, M. A.; Dorfman, K. D.; Bates, F. S. Thermal Processing of Diblock Copolymer Melts Mimics Metallurgy. Science 2017, 356, 520–523 89. Li, W.; Duan, C.; Shi, A.-C. Nonclassical Spherical Packing Phases Self-Assembled from AB-Type Block Copolymers. ACS Macro Lett. 2017, 6, 1257–1262

32

ACS Paragon Plus Environment

Page 32 of 34

Page 33 of 34 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

90. Wood, W. W.; Parker, F. R. Monte Carlo Equation of State of Molecules Interacting with the Lennard-Jones Potential. I. A Supercritical Isotherm at about Twice the Critical Temperature. J. Chem. Phys. 1957, 27, 720–733 91. Maitland, G. C.; Rigby, M.; Smith, E.; Wakeham, W.; Henderson, D. Intermolecular Forces: Their Origin and Determination; Oxford University Press, 1981 92. Jorgensen, W. L.; Maxwell, D. S.; Tirado-Rives, J. Development and Testing of the OPLS AllAtom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 1996, 118, 11225–11236 93. Abraham, M. J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J. C.; Hess, B.; Lindahl, E. GROMACS: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX 2015, 1, 19–25 94. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A. E.; Berendsen, H. J. GROMACS: Fast, Flexible, and Free. J. Comput. Chem. 2005, 26, 1701–1718 95. Nosé, S. A Unified Formulation of the Constant Temperature Molecular Dynamics Methods. J. Chem. Phys. 1984, 81, 511–519 96. Hoover, W. G. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A 1985, 31, 1695 97. Parrinello, M.; Rahman, A. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys. 1981, 52, 7182–7190 98. Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577–8593 99. Hess, B. P-LINCS: A Parallel Linear Constraint Solver for Molecular Simulation. J. Chem. Theory Comput. 2008, 4, 116–122

33

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

100. Siepmann, J. I.; Martin, M. G.; Chen, B.; Wick, C. D.; Stubbs, J. M.; Potoff, J. J.; Eggimann, B. L.; McGrath, M. J.; Zhao, X. S.; Anderson, K. E.; Rafferty, J. L.; Rai, N.; Maerzke, K. A.; Keasler, S. J.; Bai, P.; Fetisov, E. O.; Shah, M. S.; Chen, Q. P.; Xue, B.; DeJaco, R. F. et al. Monte Carlo for Complex Chemical Systems–Minnesota Version 17.1; Minneapolis, MN, 2017. 101. Chen, Q. P.; Xue, B.; Siepmann, J. I. Using the k-d Tree Data Structure to Accelerate Monte Carlo Simulations. J. Chem. Theory and Comput. 2017, 13, 1556–1565 102. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Molec. Graphics 1996, 14, 33–38 103. Stone, J. An Efficient Library for Parallel Ray Tracing and Animation; University of Missouri-Rolla: Rolla, MO, 1998 104. Brehm, M.; Kirchner, B. TRAVIS - A Free Analyzer and Visualizer for Monte Carlo and Molecular Dynamics Trajectories. J. Chem. Inf. Model. 2011, 51, 2007–2023 105. Arora, A.; Morse, D. C.; Bates, F. S.; Dorfman, K. D. Commensurability and Finite Size Effects in Lattice Simulations of Diblock Copolymers. Soft Matter 2015, 11, 4862–4867 106. Medapuram, P.; Glaser, J.; Morse, D. C. Universal Phenomenology of Symmetric Diblock Copolymers Near the Order–Disorder Transition. Macromolecules 2015, 48, 819–839

34

ACS Paragon Plus Environment

Page 34 of 34