Next Article in Journal
Equilibrium and Dynamic Surface Tension Behavior in Colloidal Unimolecular Polymers (CUP)
Next Article in Special Issue
Highly Self-Healable Polymeric Coating Materials with Enhanced Mechanical Properties Based on the Charge Transfer Complex
Previous Article in Journal
Characterization of Immunogenicity Associated with the Biocompatibility of Type I Collagen from Tilapia Fish Skin
Previous Article in Special Issue
Stimulus-Responsiveness of Thermo-Sensitive Polymer Hybridized with N-Doped Carbon Quantum Dots and Its Applications in Solvent Recognition and Fe3+ Ion Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Loop and Bridge Conformations of ABA Triblock Comb Copolymers: A Conformational Assessment for Molecular Composites

1
Department of Chemical and Biological Engineering, Korea University, Seoul 02841, Korea
2
College of Medicine, Seoul National University, Seoul 03080, Korea
3
Department of Life Sciences, Korea University, Seoul 02841, Korea
*
Authors to whom correspondence should be addressed.
Polymers 2022, 14(11), 2301; https://doi.org/10.3390/polym14112301
Submission received: 15 May 2022 / Revised: 2 June 2022 / Accepted: 2 June 2022 / Published: 6 June 2022
(This article belongs to the Special Issue Polymer Blends and Synthesis for the Fabrication of Nanocomposites)

Abstract

:
We computationally investigate the conformational behavior, “bridging” chain, between different the phase-separated domains vs “looping” chain on the same domain, for two chain architectures of ABA triblock copolymers, one with a linear architecture (L-TBC) and the other with comb architecture (C-TBC) at various segregation regimes using dissipative particle dynamics (DPD) simulations. The power-law relation between the bridge fraction ( Φ ) and the interaction parameter ( χ ) for C-TBC is found to be Φ χ 1.6 in the vicinity of the order-disorder transition ( χ O D T ), indicating a drastic conversion from the bridge to the loop conformation. When χ further increases, the bridge-loop conversions slow down to have the power law, Φ χ 0.18 , approaching the theoretical power law Φ χ 1 / 9 predicted in the strong segregation limit. The conformational assessment conducted in the present study can provide a strategy of designing optimal material and processing conditions for triblock copolymer either with linear or comb architecture to be used for thermoplastic elastomer or molecular nanocomposites.

Graphical Abstract

1. Introduction

Triblock copolymer (TBC), comprised of three polymer blocks linearly linked together in either ABA or ABC form, is an industrially important polymer used in a wide range of fields such as thermoplastic elastomer [1,2,3,4,5,6] and molecular composites, refs. [7,8,9] where the conformational behavior of TBC is often a key factor for determining the resultant performance of such applications. For instance, well-designed TBC consisting of two terminal blocks with hard segments and a middle block with rubbery segments can form hard nanodomains embedded in rubbery matrix, where chains bridging between two different hard domains, which can efficiently persist the mechanical deformation, play a critical role for the overall mechanical properties of thermoplastic elastomer or molecular composite systems [10,11,12]. For this reason, there have been steady interest in the chain conformation of TBC with respect to the preferred conformation, “bridging” between different domains vs “looping” on the same domain (Figure 1). Previous theoretical approach based on the self-consistent theory predicted the fraction of bridge conformation ( Φ ) scale as Φ ( χ N ) 1 / 9 [13] in the limit of strong segregation where χ is the Flory–Huggins interaction parameter between unlike segments and N is the degree of polymerization of TBC. The assessment of bridging fraction have been reported by several groups. Earlier experimental works reported the loop/bridge conformation ratios of TBC determined from viscoelastic, dielectric or mechanical behavior [12,14,15,16]. In the theoretical side, some computational works based on self-consistent field theory or coarse-grained simulations were also carried out to investigate the bridging fraction of the ordered TBC [13,15,17,18,19,20,21,22,23].
Comb copolymers, where two or more dissimilar types of side chains as macromers are grafted to a linear polymer backbone with a certain macromer sequence, have attracted much interest recently owing to their intriguing self-assembly behaviors [24,25,26,27,28,29]. Analogous to TBC having a triblock sequence based on monomeric units, comb copolymers can also have triblock sequence based on macromer units (Figure 1), which envisions its use for a novel type of molecular composites or thermoplastic elastomer [30,31,32] with a number of additional advantages such as more variety of options for functional design of macromers as a reinforcing component. However, despite growing interest in comb copolymers as a novel type of molecular composite materials, their conformational behavior, in particular, bridge vs loop conformations have not yet been investigated. In this brief report, as the first investigation of bridge conformation of comb copolymer, we report a simulation assessment of a bridge fraction of comb copolymer consisting of A and B macromers with an ABA triblock sequence (hereafter referred to as C-TBC) in a molten state focusing on its dependence on χ in the various segregation regime.

2. Simulation Methods

All simulations were carried out by a dissipative particle dynamics (DPD) [33,34] using the HOOMD package [35]. In the DPD model, polymers are represented by bead-spring chains, [36,37] where each bead representing a statistical monomer interacts with each other via a pairwise additive force. The force f i acting on bead i of mass m i at a position vector r i is given as:
f i = m i r ¨ i = j i F i j ( C ) + F i j ( D ) + F i j ( R ) + F i j ( S ) ,
where F i j ( C ) , F i j ( D ) , F i j ( R ) , and F i j ( S ) are a conservative force, a drag force, a random force, and spring force between bead i and j, respectively. The conservative force F i j ( C ) is modeled as a soft core repulsion,
F i j ( C ) = a i j 1 r i j R c r ^ i j for r i j < R c 0 otherwise ,
where a i j is a maximum repulsion ( a i j > 0 ) between beads i and j, r i j is the distance between bead i and j, r ^ i j is a unit vector along the direction from bead i to bead j, and R c is the cutoff distance. The drag force F i j ( D ) and the random force F i j ( R ) are given as:
F i j ( D ) = γ w ( r i j ) 2 r ^ i j · r ˙ i j r ^ i j
F i j ( R ) = ζ i j ( t ) w ( r i j ) 6 k B T γ δ t r ^ i j ,
where γ is the friction coefficient, w ( r i j ) is a weight function related to r i j , ζ i j is a random number uniformly distributed in the range of [−1,1] generated independently for each pair of bead i and j at each time step, k B T is thermal energy, and δ t is the time step size. Equations (3) and (4) ensure the consistency between kinetic energy and thermal energy via the amplitude of random noise ( 6 k B T γ / δ t ), refs. [33,38] with the weight function w chosen to have the following form:
w ( r ) = 1 r R c for r < R c 0 otherwise .
The bonding between bead i and j, responsible for chain connectivity, is taken into account by a spring force, F i j ( S ) :
F i j ( S ) = K ( r i j r o ) r ^ i j ,
where K is the spring constant and r o is the equilibrium bond length. The equation of motions (Equation (1)) for beads in the system were time-integrated using velocity-Verlet algorithm [39]:
r i ( t + δ t ) = r i ( t ) + r ˙ i ( t ) δ t + 1 2 r ¨ i ( t ) δ t 2
r ˙ i ( t + δ t ) = r ˙ i ( t ) + 1 2 r ¨ i ( t ) δ t + 1 2 r ¨ i ( t + δ t ) δ t .
The basic units for length, mass, energy, and time in the simulation are set to be R c = 1 , m = 1 , k B T = 1 , and t = R c m / k B T = 1 , respectively, and the time step δ t is set to be δ t = 0.01 which is specified from the unit thermal energy, k B T = 1 , for the consistency between thermal and kinetic energy. All DPD parameters introduced Equations (11)–(15) are rescaled according to these basic units, which are listed in Table 1.
Using the bead-spring chain model by DPD, architecturally monodisperse C-TBCs were generated in a 40 R c × 40 R c × 40 R c simulation box with a number density of beads ρ = 3 R c 3 chosen for the molten state. An ABA-type C-TBC chain consists of a backbone of M beads where A-macromers each with N 1 beads are grafted to the backbone in the two terminal backbone region each with M / 4 beads, respectively, and B macromers each with N 1 beads are grafted to the backbone in the middle backbone region with M / 2 beads. The periodic boundary conditions were applied in all axes of the simulation box. The unfavorable interaction between the beads of type A and type B was modeled using the maximum repulsion a i j introduced in Equation (2), where its value is given from the Flory interaction parameter between A- and B-bead, χ , using the relation, c χ = Δ a / k B T , where Δ a = a A B ( a A A + a B B ) / 2 and the density-dependent parameter c is given as c = 3.27 for the present choice of bead density [34]. The maximum repulsion between the same kind of beads is set to be a A A = a B B = 25 k B T / R c and that between A and B beads, a A B , is given according to a desired χ . The ordered state of each system was obtained by stepwise-increasing Δ a / k B T from an athermal state ( Δ a / k B T = 0.0 ) to a desired Δ a / k B T with the simulated annealing scheme [40] where the thermal profile ( τ k B T / Δ a ) with time t is given by:
τ ( t + Δ t ) = τ ( t ) k ( τ ( t ) ) .
Here, k ( τ ( t ) ) is the cooling schedule, i.e., the change of τ in the time interval Δ t = 4 × 10 4 δ t , which has the form of:
k ( τ ( t ) ) = k o for 0.8 τ O D T < τ ( t ) < 2.5 τ O D T 10 k o otherwise ,
where τ O D T is the τ at the ODT and k o is the cooling rate in the region nearby ODT. The value of k o in this study is set to be in the range of 0.003–0.006 depending on the statistical inefficiency analyzed for the given system [41].
Having obtained the ordered structures at desired Δ a / k B T , the systems were further equilibrated at each of Δ a / k B T for 2 × 10 6 δ t followed by the production step for 2 × 10 6 δ t to produce configuration samples for thermodynamic average.

3. Results

We consider a symmetric ABA-type C-TBC chain architecture comprised of two terminal blocks each with M / 4 A-macromers and a middle block with M / 2 B-macromers where each macromer consists of N 1 beads (Figure 1c).
It is noted that the the C-TBC chain with N = 1 reduces to the ABA triblock chain with linear architecture (hereafter referred to as L-TBC) whose conformational behavior is to be compared to that of comb architecture, i.e., C-TBC with N 1 . To investigate the conformational behavior of C-TBC chain in the disordered and ordered states, we first located the order–disorder transition (ODT) of molten C-TBC from the DPD simulations. To do this, the density fluctuation of A-bead and B-bead in the systems simulated at a desired χ was analyzed using scattering function, computed by:
S ( q ) = 1 V i < j e i q · ( r i r j ) Ψ i Ψ j ,
where q is the wave vector, V is the volume of the system, r i is the coordinates of the bead i, Ψ i is the occupation variable having values of −1 or 1 if the bead i is an A bead or a B bead, respectively, and the bracket   indicates a thermodynamic average. The order parameter was then used for the determination of ODT, which can be computed by the second order Legendre polynomial using the scattering function, P 2 = 3 2 q ( q ^ · q ^ 1 ) 2 S ( q ) / q S ( q ) 1 where q ^ and q ^ 1 are the unit vectors in the direction of a wave vector q and in the direction of the dominant wave vector q 1 , respectively. In the disordered state where χ < χ O D T , the density fluctuations of A-beads, described by wave vector q , are broadly distributed and therefore the order parameter P 2 fluctuates around zero. When the ordered phase is formed such that χ > χ O D T , a certain wave vector becomes dominant (i.e., S ( q ) is sharply peaked at the dominant wave vector q 1 ), which leads to the increase in P 2 . Figure 2 shows an example of the determination of ODT by the order parameter, plotted against Δ a for molten C-TBC with { M = 24 , N = 4 } . Table 2 summarizes the value of χ at ODT, χ O D T and the domain spacing, L, for the L-TBC and C-TBC systems simulated in this study. It should be noted that the values of N and M for L-TBC ( N = 1 ) and C-TBC ( N 1 ) samples were chosen such that the domain spacing in all copolymer samples are nearly the same to exclude the effect of domain size on the conformational behavior about the bridge and loop.
For measuring the fraction of bridge and loop conformation, we followed the method used in our previous work for L-TBC [20] where the bridge and loop fractions were estimated using the angle θ between two vectors,
cos θ = r B A · r B A | r B A | | r B A | ,
where r B A and r B A are the two vectors from the middle position B pointing to two terminal positions, A and A′, in the backbone, respectively. The fraction of bridge conformation, Φ , was then measured by the following formula,
Φ = 1 0 P ( u ) d u ,
where u cos θ and P ( u ) is the probability density function of u. The fraction of loop conformation is then given as 1 Φ .
Figure 3a shows the behavior of P ( cos θ ) for L-TBC (i.e., N = 1 ) with M = 48 at athermal state ( χ = 0 ), χ < χ O D T , χ = χ O D T , and χ χ O D T . For chains satisfying random walk so that the two vectors pointing termini r B A and r B A are randomly oriented to each other, P ( cos ϕ ) must follow a uniform distribution. However, because of excluded volume interactions, the distribution shows a slight slope even at athermal condition ( χ = 0 ) and rapidly decays at cos ϕ 1 where the two terminal A-blocks overlap with their radii of gyrations. The higher value of P ( cos θ ) at cos θ = 0.5 1.0 at χ = 0.75 when compared to that at χ = 0 is due to the formation of longer-lived transient domain at χ = 0.75 . Thus, the system at χ = 0.75 creates more bridge conformations, which is reflected by the higher value of P ( cos θ ) at cos θ = 0.5 1.0 . As the system becomes ordered ( χ > χ O D T = 0.83 ), the shape of the distribution becomes parabolic, which indicates that the L-TBC chains adopt either bridge or loop conformation as the phase-separated domains are formed. Similar χ -evolution of the distribution but more pronounced shape (i.e., steeper slopes at χ = 0 and more convex curves in the region of cos θ = 0.5 1.0 at higher χ ) were observed for C-TBC (Figure 3b), which reflect the more significant effect of excluded volume interactions owing to bulky macromers.
In Figure 4, the fractions of bridge conformation, Φ , measured by Equations (12) and (13), are plotted as a function of χ for various L- and C-TBC systems. It is noted from Figure 4 that, as the system enters the ordered regime, the bridge fraction for L-TBC ( { N = 1 , M = 48 } ) and for C-TBC with relatively shorter side chains ( { N = 2 , M = 32 } , { N = 4 , M = 24 } ) scales similarly as χ 0.08 , showing slightly smaller exponent than the strong segregation theory result, ν = 1 / 9 0.11 for χ ν , predicted for L-TBC [13]. An interesting result is found for C-TBC with relatively longer side chains ( { N = 8 , M = 16 } ) where the bridge fraction decreases rapidly as χ increases in the vicinity to χ O D T but still in the disordered region, following the χ 0.55 power law. In this disordered region near ODT, long-lived but still transient A-domains are responsible for such uneven conformational behavior. The exponent becomes even larger when this C-TBC system enters the ordered regime showing a very rapid decay of Φ χ 1.6 followed by an another regime of slow conversion from loop to bridge, Φ χ 0.18 . It is also of interest that the bridge conformation for C-TBC with a longer side chain is preferred over the loop conformation in the very low- χ regime ( Φ 0.65 at χ χ O D T ) and becomes rapidly unfavored as χ passes χ O D T ( Φ 0.35 at χ χ O D T ), which is in contrast with the behavior of L-TBC showing Φ 0.52 at χ χ O D T and Φ 0.5 at χ χ O D T . We interpret this conformational behavior of C-TBC as follows. In the very low χ regime, the loop conformation of C-TBC is less favored due to the steric repulsion between two terminal blocks consisting of A-macromers, resulting in a larger fraction of bridge conformation. As χ increases so that the transient domains at χ < χ O D T or the phase-separated domains at χ > χ O D T are formed, the bridge conformation starts to be less preferred, because the number of possible ways for backbone paths, which is obliged to be across the midplane of B-domain, is smaller than that of the loop conformation which does not have to pass through the midplane (See schematic in Figure 1a,b) [31,42]. This different restriction on bridge and loop conformation imposed by the phase-separated domains gives rise to the preference of the loop conformation over the bridge conformation as χ increases and this propensity becomes even more pronounced for comb architecture by its side chain which can save the stretching energy of looping C-TBC. In general, the domain spacing of the block copolymer system is a result of the balance between the minimal interfacial energy and the minimal stretching energy: the former wants to expand the domain spacing while the latter wants to compress the spacing. It can be noticed that the difference in backbone stretching energy between the loop conformation and the bridge conformation for C-TBC is larger than that for L-TBC because the stretching energy of backbone of a looping C-TBC can be saved more than that of a looping L-TBC due to the presence of side chains filling in the region nearby the midplane.

4. Conclusions

In conclusion, the computational assessments of bridge and loop conformation for two chain architectures of ABA triblock copolymers, one with a linear architecture (L-TBC) and the other with comb architecture (C-TBC), were performed at various segregation regimes using DPD simulations. It is found from DPD simulation that the bridge conformation of C-TBC favored in the very low χ regime is converted to the loop conformation much more rapidly than the bridge-loop conversion in L-TBC. The power law relation between the bridge fraction ( Φ ) and the interaction parameter ( χ ) for C-TBC is found to be Φ χ 1.6 in the ordered state but in the vicinity of χ O D T indicating a drastic conversion from the bridge to the loop as the system enters the regime of the ordered state. When the segregation power becomes stronger such that χ χ O D T , the bridge–loop conversions slow down to the power law, Φ χ 0.18 , presumably approaching the theoretical power law Φ χ 1 / 9 in the strong segregation limit. Assuming that a large amount of bridge conformation in the TBC system is a key requirement for the reliable mechanical performance of TBC materials to be used for thermoplastic elastomers or molecular composites, the finding for C-TBC in the present study is somewhat discouraging, suggesting that the prolonged equilibration of the C-TBC system could cause an undesirable TBC system containing the loop conformation dominant over bridge conformation. However, noting in C-TBC that the bridge conformation is dominant in the low χ regime, finding an efficient thermal profile for manipulating the C-TBC system can be a solution to detour the rapid bridge–loop conversion, which will be investigated in our future work.

Author Contributions

Conceptualization, J.B. and J.H.; methodology, J.P., J.-Y.J. and J.H.; software, J.P., J.-Y.J. and J.H.; validation, H.-W.S., J.-W.P., J.B. and J.H.; formal analysis, J.P., J.-Y.J. and J.H.; investigation, J.P., J.-Y.J. and J.H.; resources, J.P., J.-Y.J. and J.H.; data curation, J.P., J.-Y.J. and J.H.; writing—original draft preparation, J.H.; writing—review and editing, J.H.; visualization, J.P. and J.H.; supervision, J.B. and J.H.; project administration, J.H.; funding acquisition, H.-W.S., J.-W.P. and J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Research Foundation of Korea via grant number 2020R1A4A2002903, 2020R1F1A1065951 and by the Creative Materials Discovery Program funded by the Ministry of Science and ICT via grant number NRF-2018M3D1A1058536.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Holden, G.L.; Legge, N.R.; Quirk, R.P.; Schroeder, H.E. Thermoplastic Elastomers; Hanser: Munich, Germany, 1996. [Google Scholar]
  2. Ding, W.; Wang, S.; Yao, K.; Ganewatta, M.S.; Tang, C.; Robertson, M.L. Physical Behavior of Triblock Copolymer Thermoplastic Elastomers Containing Sustainable Rosin-Derived Polymethacrylate End Blocks. ACS Sustain. Chem. Eng. 2017, 5, 11470–11480. [Google Scholar] [CrossRef]
  3. Gregory, G.L.; Sulley, G.S.; Carrodeguas, L.P.; Chen, T.P.D.; Santmarti, A.; Terrill, N.J.; Lee, K.Y.; Williams, C.K. Triblock Polyester Thermoplastic Elastomers with Semi-Aromatic Polymer End Blocks by Ring-Opening Copolymerization. Chem. Sci. 2020, 11, 6567–6581. [Google Scholar] [CrossRef]
  4. Jia, M.; Zhang, D.; de Kort, G.W.; Wilsens, C.H.R.M.; Rastogi, S.; Hadjichristidis, N.; Gnanou, Y.; Feng, X. All-Polycarbonate Thermoplastic Elastomers Based on Triblock Copolymers Derived from Triethylborane-Mediated Sequential Copolymerization of CO2 with Various Epoxides. Macromolecules 2020, 53, 5297–5307. [Google Scholar] [CrossRef] [PubMed]
  5. Kalita, U.; Samanta, S.; Banerjee, S.L.; Das, N.C.; Singha, N.K. Biobased Thermoplastic Elastomer Based on an SMS Triblock Copolymer Prepared via RAFT Polymerization in Aqueous Medium. Macromolecules 2021, 54, 1478–1488. [Google Scholar] [CrossRef]
  6. Kawarazaki, I.; Hayashi, M. Enhancement of Mechanical Properties of ABA Triblock CopolymerBased Elastomers by Incorporating Partial Cross-Links on the Soft Bridge Chains. ACS Appl. Polym. Mater. 2021, 3, 1271–1275. [Google Scholar] [CrossRef]
  7. Wu, H.; Thakur, V.K.; Kessler, M.R. Novel Low-Cost Hybrid Composites from Asphaltene/SBS Triblock Copolymer with Improved Thermal and Mechanical Properties. J. Mater. Sci. 2016, 51, 2394–2403. [Google Scholar] [CrossRef]
  8. Thakur, V.K.; Grewell, D.; Thunga, M.; Kessler, M.R. Novel Composites from Eco-Friendly SoyFlour/SBS Triblock Copolymer. Macromol. Mater. Eng. 2014, 299, 953–958. [Google Scholar] [CrossRef]
  9. Zhang, Y.; Chen, S.; Zhang, H.; Ding, X.; Fu, P.; Du, F. Fabrication of Conjugated Triblock Copolymer/Single-Walled Carbon Nanotubes Composite Films with Enhanced Thermoelectric Performance. Compos. Commun. 2021, 27, 10083. [Google Scholar] [CrossRef]
  10. Watanabe, H.; Sato, T.; Osaki, K.; Yao, M.L.; Yamagishi, A. Rheological and Dielectric Behavior of a Styrene-Isoprene-Styrene Triblock Copolymer in Selective Solvents. 2. Contribution of Loop-Type Middle Blocks to Elasticity and Plasticity. Macromolecules 1997, 30, 5877–5892. [Google Scholar] [CrossRef]
  11. Honeker, C.C.; Thomas, E.L. Impact of Morphological Orientation in Determining Mechanical Properties in Triblock Copolymer Systems. Chem. Mater. 1996, 8, 1702–1714. [Google Scholar] [CrossRef]
  12. Takano, A.; Kamaya, I.; Takahashi, Y.; Matsushita, Y. Effect of Loop/Bridge Conformation Ratio on Elastic Properties of the Sphere-Forming ABA Triblock Copolymers: Preparation of Samples and Determination of Loop/Bridge Ratio. Macromolecules 2005, 38, 9718–9723. [Google Scholar] [CrossRef]
  13. Milner, S.T.; Witten, T.A. Bridging Attraction by Telechelic Polymers. Macromolecules 1992, 25, 5495–5503. [Google Scholar] [CrossRef]
  14. Watanabe, H. Slow Dielectric Relaxation of a Styrene—Isoprene—Styrene Triblock Copolymer with Dipole Inversion in the Middle Block: A Challenge to a Loop/Bridge Problem. Macromolecules 1995, 28, 5006–5011. [Google Scholar] [CrossRef]
  15. Karatasos, K.; Anastasiadis, S.H.; Pakula, T.; Wayanabe, H. On the Loops-to-Bridges Ratio in Ordered Triblock Copolymers: An Investigation by Dielectric Relaxation Spectroscopy and Computer Simulations. Macromolecules 2000, 33, 523–541. [Google Scholar] [CrossRef]
  16. Mori, Y.; Lim, L.S.; Bates, F.S. Consequences of Molecular Bridging in Lamellae-Forming Triblock/Pentablock Copolymer Blends. Macromolecules 2003, 36, 9879–9898. [Google Scholar] [CrossRef]
  17. Zhulina, E.B.; Halperin, A. Lamellar Mesogels and Mesophases: A Self-Consistent-Field Theory. Macromolecules 1992, 25, 5730–5741. [Google Scholar] [CrossRef]
  18. Avalas, J.B.; Johner, A.; Joanny, J.F. Bridging by Adsorbed Polymers Between Two Surfaces. J. Chem. Phys. 1994, 101, 9181–9194. [Google Scholar] [CrossRef]
  19. Matsen, M.W.; Schick, M. Lamellar Phase of a Symmetric Triblock Copolymer. Macromolecules 1994, 27, 187–192. [Google Scholar] [CrossRef]
  20. Huh, J.; Jo, W.H.; Brinke, G.T. Conformational Analysis in ABA Triblock Melts by Monte Carlo Simulation. Macromolecules 2002, 35, 2413–2416. [Google Scholar] [CrossRef] [Green Version]
  21. Aoyagi, T.; Honda, T.; Doi, M. Microstructural Study of Mechanical Properties of the ABA Triblock Copolymer Using Self-Consistent Field and Molecular Dynamics. J. Chem. Phys. 2002, 117, 8153–8161. [Google Scholar] [CrossRef]
  22. Abu-Sharkh, B.; AlSunaidi, A. Morphology and Conformation Analysis of Self-Assembled Triblock Copolymer Melts. Macromol. B Theory Simul. 2006, 15, 507–515. [Google Scholar] [CrossRef]
  23. Morita, H.; Miyamoto, A.; Kotani, M. Recoverably and Destructively Deformed Domain Structures in Elongation Process of Thermoplastic Elastomer Analyzed by Graph Theory. Polymer 2020, 188, 122098. [Google Scholar] [CrossRef]
  24. Gu, W.; Huh, J.; Hong, S.W.; Sveinbjornsson, B.R.; Park, C.; Grubbs, R.H.; Russell, T.P. Self-Assembly of Symmetric Brush Diblock Copolymers. ACS Nano 2013, 7, 2551–2558. [Google Scholar] [CrossRef] [Green Version]
  25. Rzayev, J. Synthesis of Polystyrene-Polylactide Bottlebrush Block Copolymers and Their Melt Self-Assembly into Large Domain Nanostructures. Macromolecules 2009, 42, 2135–2141. [Google Scholar] [CrossRef]
  26. Bolton, J.; Bailey, T.S.; Rzayev, J. Large Pore Size Nanoporous Materials from the Self-Assembly of Asymmetric Bottlebrush Block Copolymers. Nano Lett. 2011, 11, 998–1001. [Google Scholar] [CrossRef]
  27. Fenyves, R.; Schmutz, M.; Horner, I.J.; Bright, F.V.; Rzayev, J. Aqueous Self-Assembly of Giant Bottlebrush Block Copolymer Surfactants as Shape-Tunable Building Blocks. J. Am. Chem. Soc. 2014, 136, 7762–7770. [Google Scholar] [CrossRef]
  28. Song, D.P.; Zhao, T.H.; Guidetti, G.; Vignolini, S.; Parker, R.M. Hierarchical Photonic Pigments via the Confined Self-Assembly of Bottlebrush Block Copolymers. ACS Nano 2019, 13, 1764–1771. [Google Scholar] [CrossRef] [Green Version]
  29. Kim, E.J.; Shin, J.J.; Do, T.; Lee, G.S.; Park, J.; Thapar, V.; Choi, J.; Bang, J.; Yi, G.R.; Hur, S.M.; et al. Molecular Weight Dependent Morphological Transitions of Bottlebrush Block Copolymer Particles: Experiments and Simulations. ACS Nano 2021, 15, 5513–5522. [Google Scholar] [CrossRef]
  30. Yang, J.X.; Long, Y.Y.; Pan, L.; Men, Y.F.; Li, Y.S. Spontaneously Healable Thermoplastic Elastomers Achieved through One-Pot Living Ring-Opening Metathesis Copolymerization of Well-Designed Bulky Monomers. ACS Appl. Mater. Interfaces 2016, 8, 12445–12455. [Google Scholar] [CrossRef]
  31. Sunday, D.F.; Chang, A.B.; Liman, C.D.; Gann, E.; Delongchamp, D.M.; Thomsen, L.; Matsen, M.W.; Grubbs, R.H.; Soles, C.L. Self-Assembly of ABC Bottlebrush Triblock Terpolymers with Evidence for Looped Backbone Conformations. Macromolecules 2018, 51, 7178–7185. [Google Scholar] [CrossRef]
  32. Nian, S.; Lian, H.; Gong, Z.; Zhernenkov, M.; Qin, J.; Cai, L.H. Molecular Architecture Directs Linear-Bottlebrush-Linear Triblock Copolymers to Self-Assemble to Soft Reprocessable Elastomers. ACS Macro Lett. 2019, 8, 1528–1534. [Google Scholar] [CrossRef]
  33. Español, P.; Warren, P.B. Statistical Mechanics of Dissipative Particle Dynamics. Europhys. Lett. 1995, 30, 191–196. [Google Scholar] [CrossRef] [Green Version]
  34. Groot, R.D.; Warren, P.B. Dissipative Particle Dynamics: Bridging the Gap between Atomistic and Mesoscopic Simulation. J. Chem. Phys. 1997, 107, 4423–4435. [Google Scholar] [CrossRef]
  35. Andersona, J.A.; Glaser, J.; Glotzer, S.C. HOOMD-blue: A Python Package for High-Performance Molecular Dynamics and Hard Particle Monte Carlo Simulations. Comput. Mater. Sci. 2020, 173, 109363. [Google Scholar] [CrossRef] [Green Version]
  36. Groot, R.D.; Madden, T.J. Dynamic Simulation of Diblock Copolymer Microphase Separation. J. Chem. Phys. 1998, 108, 8713–8724. [Google Scholar] [CrossRef]
  37. Spenley, N.A. Scaling Laws for Polymers in Dissipative Particle Dynamics. Europhys. Lett. 2000, 49, 534–540. [Google Scholar] [CrossRef]
  38. Phillips, C.L.; Anderson, J.A.; Glotzer, S.C. Pseudo-Random Number Generation for Brownian Dynamics and Dissipative Particle Dynamics Simulations on GPU Devices. J. Comput. Phys. 2011, 230, 7191–7201. [Google Scholar] [CrossRef]
  39. Allen, M.P.; Tildesley, D.J. Computer Simulation of Liquids, 2nd ed.; Clarendon Press: Oxford, UK, 2017. [Google Scholar]
  40. Karabin, M.; Stuartb, S.J. Simulated Annealing with Adaptive Cooling Rates. J. Chem. Phys. 2020, 153, 114103. [Google Scholar] [CrossRef]
  41. Bishop, M.; Frinks, S. Error Analysis in Computer Simulations. J. Chem. Phys. 2020, 87, 3675–3676. [Google Scholar] [CrossRef]
  42. Mikhailov, I.V.; Leermakers, F.A.; Darinskii, A.A.; Zhulina, E.B.; Borisov, O.V. Theory of microphase segregation in ABA triblock comb-shaped copolymers: Lamellar mesophase. Macromolecules 2021, 54, 4747–4759. [Google Scholar] [CrossRef]
Figure 1. Illustration of (a) bridge and loop conformation for ABA linear triblock copolymer (L-TBC), (b) those of ABA comb triblock copolymer (C-TBC) in the phase-separated domains, and (c) the chain architecture of C-TBC investigated in this study. In (a,b), the red-colored regions and the blue-colored region are the A- and the B-domain, respectively, and the dashed line in the B-domain represents the midplane.
Figure 1. Illustration of (a) bridge and loop conformation for ABA linear triblock copolymer (L-TBC), (b) those of ABA comb triblock copolymer (C-TBC) in the phase-separated domains, and (c) the chain architecture of C-TBC investigated in this study. In (a,b), the red-colored regions and the blue-colored region are the A- and the B-domain, respectively, and the dashed line in the B-domain represents the midplane.
Polymers 14 02301 g001
Figure 2. The order parameter versus the interaction parameter for the C-TBC with { M = 24 , N = 4 } . The open circles and filled circles represent the points where the disordered and the ordered phases are stable, respectively, and the blue solid line is fit to a three-parameter sigmoidal function. The inset images show the two example structures simulated at the disordered and ordered region.
Figure 2. The order parameter versus the interaction parameter for the C-TBC with { M = 24 , N = 4 } . The open circles and filled circles represent the points where the disordered and the ordered phases are stable, respectively, and the blue solid line is fit to a three-parameter sigmoidal function. The inset images show the two example structures simulated at the disordered and ordered region.
Polymers 14 02301 g002
Figure 3. P ( cos θ ) for (a) L-TBC (i.e., N = 1 ) with M = 48 and for (b) C-TBC with M = 16 and M = 8 at various χ values. The χ O D T for (a,b) were χ = 1.08 and χ = 0.72 , respectively.
Figure 3. P ( cos θ ) for (a) L-TBC (i.e., N = 1 ) with M = 48 and for (b) C-TBC with M = 16 and M = 8 at various χ values. The χ O D T for (a,b) were χ = 1.08 and χ = 0.72 , respectively.
Polymers 14 02301 g003
Figure 4. The fractions of bridge conformation, Φ , versus χ for the simulated L-TBC and C-TBC samples. The vertical dashed lines represent the χ O D T for the corresponding copolymer samples.
Figure 4. The fractions of bridge conformation, Φ , versus χ for the simulated L-TBC and C-TBC samples. The vertical dashed lines represent the χ O D T for the corresponding copolymer samples.
Polymers 14 02301 g004
Table 1. The list of the DPD parameters used in the present study.
Table 1. The list of the DPD parameters used in the present study.
ParameterValueUnit 1Equations
a i i ( a A A , a B B ) 25.0 k B T / R c (2)
γ 4.5 m k B T / R c (3) and (4)
K100.0 k B T / R c 2 (6)
r o 0.7 R c (6)
1 The basic units for length, mass, and energy are set to be Rc = 1, m = 1, and kBT = 1, respectively.
Table 2. The chain architectures of triblock copolymers simulated in this study and the list of the χ values at ODT, χ O D T and the domain spacing, L, for these triblock copolymers in molten state.
Table 2. The chain architectures of triblock copolymers simulated in this study and the list of the χ values at ODT, χ O D T and the domain spacing, L, for these triblock copolymers in molten state.
ArchitecureMN χ ODT L1
L-TBC4811.088.26 ± 0.19
C-TBC3220.998.18 ± 0.26
C-TBC2440.778.17 ± 0.16
C-TBC1680.728.26 ± 0.18
1 measured at χ MN = 150.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Park, J.; Jung, J.-Y.; Shin, H.-W.; Park, J.-W.; Bang, J.; Huh, J. Loop and Bridge Conformations of ABA Triblock Comb Copolymers: A Conformational Assessment for Molecular Composites. Polymers 2022, 14, 2301. https://doi.org/10.3390/polym14112301

AMA Style

Park J, Jung J-Y, Shin H-W, Park J-W, Bang J, Huh J. Loop and Bridge Conformations of ABA Triblock Comb Copolymers: A Conformational Assessment for Molecular Composites. Polymers. 2022; 14(11):2301. https://doi.org/10.3390/polym14112301

Chicago/Turabian Style

Park, Jihoon, Je-Yeon Jung, Hyun-Woo Shin, Jong-Wan Park, Joona Bang, and June Huh. 2022. "Loop and Bridge Conformations of ABA Triblock Comb Copolymers: A Conformational Assessment for Molecular Composites" Polymers 14, no. 11: 2301. https://doi.org/10.3390/polym14112301

APA Style

Park, J., Jung, J. -Y., Shin, H. -W., Park, J. -W., Bang, J., & Huh, J. (2022). Loop and Bridge Conformations of ABA Triblock Comb Copolymers: A Conformational Assessment for Molecular Composites. Polymers, 14(11), 2301. https://doi.org/10.3390/polym14112301

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop