Date of Award
12-20-2024
Date Published
January 2023
Degree Type
Dissertation
Degree Name
Doctor of Philosophy (PhD)
Department
Chemistry
Advisor(s)
Radhakrishna Sureshkumar
Second Advisor
Shalabh Maroo
Keywords
copolymer;molecular dynamics;nanovesicle;self assembly;shear-induced structure;Vesicle Morphogenesis
Subject Categories
Chemical Engineering | Engineering
Abstract
Block copolymers (BCPs) have been a focal point of research for several decades. They have been engineered to produce various morphologies by altering the combinations of chemically distinct polymer segments that interact selectively with their environment, such as in aqueous or organic solvents, or within a polymer matrix. The ability to tailor morphological features by manipulating polymer chain length, composition, and monomer chemistry has led to applications in advanced material manufacturing, catalysis, emulsification, environmental remediation, targeted drug delivery, gene therapy, and medical diagnostics. Consequently, understanding the relationship between copolymer architecture and the self-assembled morphologies at equilibrium, the pathways of structure evolution and the stability of equilibrium structures to perturbations in environmental variables such as flow shear and temperature, is crucial for designing polymers for diverse applications. This study employs CGMD simulations, focusing on prototypical amphiphilic copolymers composed of Poly(Butadiene) (PB) and Poly(Ethylene Oxide) (PEO). The behavior of single polymer chains in aqueous solutions is investigated. The dependence of the radius of gyration and configurational relaxation time of the copolymers on chain length L and PB/PEO ratio are studied. The study then examines the self-assembly of PB-PEO diblock copolymers. Across varying copolymer compositions (L and PB/PEO ratio), a diverse phase portrait of emergent morphologies is observed, including simple structures such as spherical micelles, vesicles (polymersomes), bilayer lamellae, linear wormlike micelles, and tori, as well as complex forms like branched micellar networks and composite aggregates. The topological and compositional details of these morphologies are quantitatively characterized. Further, the energetic and entropic metrics underlying morphology selection are analyzed. The simulation results align well with experimental observations, offering insights into the mechanisms that render morphological diversity in amphiphilic copolymer systems. The study also investigates how structural changes occur with variations in concentration and temperature. Further, a mechanism of vesicle formation by water diffusion into a spherical micelle is discovered. The study further focuses on the mechanism of vesicle formation by the self-assembly of AB and BAB type BCPs, where A and B represent the hydrophilic (PEO) and hydrophobic (PB) segments, respectively. For AB BCPs, evolution follows the order of spherical and rod-like micelles, wormlike structures, lamellae, cavities, and vesicles. BAB BCPs first form interconnected copolymer networks, then deform to lamellar cages, ultimately transitioning into stable vesicles. Molecular reorganization at constant aggregation number to reduce solvent accessible surface area (SASA) and consequently the unfavorable hydrophobic interactions is a common motif of vesiculation. Subsequently, this work investigates the behavior of stable BAB vesicles under shear flow. Non-equilibrium molecular dynamics simulations reveal how unilamellar triblock copolymer vesicles respond to external hydrodynamic forces. Their deformation mode is affected by Weissenberg number Wi, defined as the ratio of the time scale of vesicle shape fluctuations to the inverse shear rate. For Wi < 10, a spherical vesicle deforms into a flow-aligned ellipsoidal bilayer executing tank-treading motion. For Wi > 10, pronounced variations in bilayer thickness and polymer extension manifest along the contour of the elongated vesicle, which breaks up into lamellar fragments. Below a critical strain, the deformed vesicle upon flow cessation returns to initial spherical morphology. However, for larger strains, structure reorganization after flow stoppage results in the formation of a Novel Equilibrated Shear-Induced Structure (NoESIS), in which two vesicles are connected by a dynamic molecular bridge which can accommodate additional layers of copolymers leading to a reduction in the polymer-solvent interface area. At even larger strains, the deformed vesicle breaks apart after flow cessation and eventually reorganizes into a smaller vesicle and a composite vesicular structure.
Access
Open Access
Recommended Citation
Liu, Senyuan, "THERMODYNAMICS AND SHEAR FLOW DYNAMICS OF BLOCK COPOLYMER MORPHOLOGIES IN SOLUTION" (2024). Dissertations - ALL. 2046.
https://surface.syr.edu/etd/2046
Video S4.1. Evolution of vesicle (C = 10.1 mol%, N = 11, xPEO=0.55 )
LiuS2025AppendixB.gif (7164 kB)
Video S4.2. Evolution of spherical micelle (C = 11.7 mol% , N = 13 , xPEO=0.15
LiuS2025AppendixC.gif (6967 kB)
Video S4.3. Evolution of irregular shape and wormlike micelle with Y branch (C = 10.1 mol%, N = 11, xPEO = 0.73 )
LiuS2025AppendixD.gif (9495 kB)
Video S4.4. Evolution of flexible rodlike micelle (C = 11.7 mol% , N=13 , xPEO=0.62 )
LiuS2025AppendixE.GIF (45285 kB)
Video S5.1. The vesiculation process
LiuS2025AppendixF.GIF (3342 kB)
Video S5.2. 3D visualization of wormlike micelle intermediate
LiuS2025AppendixG.GIF (3271 kB)
Video S5.3. 3D visualization of rectangular lamella intermediate
LiuS2025AppendixH.GIF (4180 kB)
Video S5.4. 3D visualization of disk lamella intermediate
LiuS2025AppendixI.GIF (5756 kB)
Video S5.5. 3D visualization of cavity micelle intermediate
LiuS2025AppendixJ.GIF (6278 kB)
Video S5.6. 3D visualization of vesicle
LiuS2025AppendixK.mp4 (2140 kB)
Video S7.1. Flow startup and cessation, γmax = 15 for Wi = 12
LiuS2025AppendixL.mp4 (4790 kB)
Video S7.2. Flow startup and cessation, γmax = 20 for Wi = 12
LiuS2025AppendixM.mp4 (1909 kB)
Video S7.3. Flow startup and cessation, γmax = 25 for Wi = 12
LiuS2025AppendixN.mp4 (7282 kB)
Video S7.4. Copolymer chain motion within the NoESIS bilayer