Lattice model of equilibrium polymerization. VII. Understanding the role of "cooperativity" in self-assembly

Stringlet excitation model of the boson peak

The boson peak (BP), a low-energy excess in the vibrational density of states over the Debye contribution, is often identified as a characteristic of amorphous solid materials. Despite decades of efforts, its microscopic origin still remains a mystery. Recently, it has been proposed, and corroborated with simulations, that the BP might stem from intrinsic localized modes involving one-dimensional (1D) string-like excitations ("stringlets"). We build on a theory originally proposed by Lund that describes the localized modes as 1D vibrating strings, but we specify the stringlet size distribution to be exponential, as observed in simulations. We provide an analytical prediction for the BP frequency ωBP in the temperature regime well below the observed glass transition temperature Tg. The prediction involves no free parameters and accords quantitatively with prior simulation observations in 2D and 3D model glasses based on inverse power law potentials. The comparison of the string model to observations is more uncertain when compared to simulations of an Al-Sm metallic glass material at temperatures well above Tg. Nonetheless, our stringlet model of the BP naturally reproduces the softening of the BP frequency upon heating and offers an analytical explanation for the experimentally observed scaling with the shear modulus in the glass state and changes in this scaling in simulations of glass-forming liquids. Finally, the theoretical analysis highlights the existence of a strong damping for the stringlet modes above Tg, which leads to a large low-frequency contribution to the 3D vibrational density of states, observed in both experiments and simulations.

Thermoreversible Gelation with Supramolecularly Polymerized Cross-Link Junctions

Structure and reversibility of cross-link junctions play pivotal roles in determining the nature of thermoreversible gelation and dynamic mechanical properties of the produced polymer networks. We attempt to theoretically explore new types of sol-gel transitions with mechanical sharpness by allowing cross-links to grow without upper bound. We consider thermoreversible gelation of the primary molecules R{Af} carrying the number f of low molecular weight functional groups (gelators) A. Gelators A are assumed to form supramolecular assemblies. Some examples are: telechelic polymers (f=2) carrying π-π stacking benzene derivatives at their both ends, and trifunctional star molecules (f=3) bearing multiple hydrogen-bonding gelators. The sol-gel transition of the primary molecules becomes sharper with the cooperativity parameter of the stepwise linear growth of the cross-links. There is a polymerization transition (crossover without singularity) of the junctions in the postgel region after the gel point is passed. If the gelator A tends to form supramolecular rings competitively with linear chains, there is another phase transition in the deep postgel region where the average molecular weight of the rings becomes infinite (Bose-Einstein condensation of rings). As a typical example of binary cross-links where gelators A and B form mixed junctions, we specifically consider metal-coordinated binding of ligands A by metal ions B. Two types of multi-nuclear supramolecular complexes are studied: (i) linear stacking (ladder) of the sandwich A2B units, and (ii) linear train of egg-box A4B units. To find the strategy towards experimental realization of supramolecular cross-links, the average molecular weight, the gel fraction, the average length of the cross-link junctions are numerically calculated for all of these models as functions of the functionality f, the concentration of the solute molecules, and the temperature. Potential candidates for the realization of these new types of thermoreversible gelation are discussed.

Approach to hyperuniformity in a metallic glass-forming material exhibiting a fragile to strong glass transition

We investigate a metallic glass-forming (GF) material (Al₉₀Sm₁₀) exhibiting a fragile-strong (FS) glass-formation by molecular dynamics simulation to better understand this highly distinctive pattern of glass-formation in which many of the usual phenomenological relations describing relaxation times and diffusion of ordinary GF liquids no longer apply, and where instead genuine thermodynamic features are observed in response functions and little thermodynamic signature is exhibited at the glass transition temperature, T_g. Given the many unexpected similarities between the thermodynamics and dynamics of this metallic GF material with water, we first focus on the anomalous static scattering in this liquid, following recent studies on water, silicon and other FS GF liquids. We quantify the "hyperuniformity index" H of our liquid, which provides a quantitative measure of molecular "jamming". To gain insight into the T-dependence and magnitude of H, we also estimate another more familiar measure of particle localization, the Debye-Waller parameter 〈u²〉 describing the mean-square particle displacement on a timescale on the order of the fast relaxation time, and we also calculate H and 〈u²〉 for heated crystalline Cu. This comparative analysis between H and 〈u²〉 for crystalline and metallic glass materials allows us to understand the critical value of H on the order of 10^-3 as being analogous to the Lindemann criterion for both the melting of crystals and the "softening" of glasses. We further interpret the emergence of FS GF and liquid-liquid phase separation in this class of liquids to arise from a cooperative self-assembly process in the GF liquid.

Extended Wertheim theory predicts the anomalous chain length distributions of divalent patchy particles under extreme confinement

Colloidal patchy particles with divalent attractive interaction can self-assemble into linear polymer chains. Their equilibrium properties in 2D and 3D are well described by Wertheim's thermodynamic perturbation theory which predicts a well-defined exponentially decaying equilibrium chain length distribution. In experi- mental realizations, due to gravity, particles sediment to the bottom of the suspension forming a monolayer of particles with a gravitational height smaller than the particle diameter. In accordance with experiments, an anomalously high monomer concentration is observed in simulations which is not well understood. To account for this observation, we interpret the polymerization as taking place in a highly confined quasi-2D plane and extend the Wertheim thermodynamic perturbation theory by defining addition reactions constants as functions of the chain length. We derive the theory, test it on simple square well potentials, and apply it to the experimental case of synthetic colloidal patchy particles immersed in a binary liquid mixture that are described by an accurate effective critical Casimir patchy particle potential. The important interaction parameters entering the theory are explicitly computed using the integral method in combination with Monte Carlo sampling. Without any adjustable parameter, the predictions of the chain length distribution are in excellent agreement with explicit simulations of self-assembling particles. We discuss generality of the approach, and its application range.

Comparative experimental and computational study of synthetic and natural bottlebrush polyelectrolyte solutions

We systematically investigate model synthetic and natural bottlebrush polyelectrolyte solutions through an array of experimental techniques (osmometry and neutron and dynamic light scattering) along with molecular dynamics simulations to characterize and contrast their structures over a wide range of spatial and time scales. In particular, we perform measurements on solutions of aggrecan and the synthetic bottlebrush polymer, poly(sodium acrylate), and simulations of solutions of highly coarse-grained charged bottlebrush molecules having different degrees of side-branch density and inclusion of an explicit solvent and ion hydration effects. While both systems exhibit a general tendency toward supramolecular organization in solution, bottlebrush poly(sodium acrylate) solutions exhibit a distinctive "polyelectrolyte peak" in their structure factor, but no such peak is observed in aggrecan solutions. This qualitative difference in scattering properties, and thus polyelectrolyte solution organization, is attributed to a concerted effect of the bottlebrush polymer topology and the solvation of the polymer backbone and counterions. The coupling of the polyelectrolyte topological structure with the counterion distribution about the charged polymer molecules along with direct polymer segmental hydration makes their solution organization and properties "tunable," a phenomenon that has significant ramifications for biological function and disease as well as for numerous materials applications.

Dynamic heterogeneity, cooperative motion, and Johari-Goldstein [Formula: see text]-relaxation in a metallic glass-forming material exhibiting a fragile-to-strong transition

We investigate the Johari-Goldstein (JG) [Formula: see text]-relaxation process in a model metallic glass-forming (GF) material ([Formula: see text]), previously studied extensively by both frequency-dependent mechanical measurements and simulation studies devoted to equilibrium properties, by molecular dynamics simulations based on validated and optimized interatomic potentials with the primary aim of better understanding the nature of this universal relaxation process from a dynamic heterogeneity (DH) perspective. The present relatively low temperature and long-time simulations reveal a direct correspondence between the JG [Formula: see text]-relaxation time [Formula: see text] and the lifetime of the mobile particle clusters [Formula: see text], defined as in previous DH studies, a relationship dual to the corresponding previously observed relationship between the [Formula: see text]-relaxation time [Formula: see text] and the lifetime of immobile particle clusters [Formula: see text]. Moreover, we find that the average diffusion coefficient D nearly coincides with [Formula: see text] of the smaller atomic species (Al) and that the 'hopping time' associated with D coincides with [Formula: see text] to within numerical uncertainty, both trends being in accord with experimental studies. This indicates that the JG [Formula: see text]-relaxation is dominated by the smaller atomic species and the observation of a direct relation between this relaxation process and rate of molecular diffusion in GF materials at low temperatures where the JG [Formula: see text]-relaxation becomes the prevalent mode of structural relaxation. As an unanticipated aspect of our study, we find that [Formula: see text] exhibits fragile-to-strong (FS) glass formation, as found in many other metallic GF liquids, but this fact does not greatly alter the geometrical nature of DH in this material and the relation of DH to dynamical properties. On the other hand, the temperature dependence of the DH and dynamical properties, such as the structural relaxation time, can be significantly altered from 'ordinary' GF liquids.

Mass-Balance Models for Scrutinizing Supramolecular (Co)polymerizations in Thermodynamic Equilibrium

Recent years have witnessed increasing attention on supramolecular polymerization, i.e., the formation of one-dimensional aggregates in which the monomeric units bind together via reversible and usually highly directional non-covalent interactions. Because of the presence of these reversible interactions, such as hydrogen bonding, π-π interactions, or metal coordination, supramolecular polymers exhibit numerous desirable properties ranging from high thermoresponsiveness to self-healing and great capacity for processability and recycling. These properties relate to intriguing experimentally observed nonlinear effects such as the monomer-dependent presence of a critical temperature for aggregation and a solvent- and temperature-tunable aggregate morphology. For coassemblies this is complemented with monomer-ratio- and monomer-compatibility-dependent internal order as well as majority-rules-type chiral amplification. However, the dynamic nature of the (co)polymers and the intricate interplay of many interactions make these effects difficult to rationalize without theoretical models. This Account presents recent advances in the development and use of equilibrium models for supramolecular copolymerization based on mass balances, mainly developed by our group. The basic idea of these models is that we describe a supramolecular (co)polymerization by a set of independent equilibrium reactions, like monomer associations and dissociations, and that in thermodynamic equilibrium the concentrations of the reactants and products in each reaction are coupled via the equilibrium constant of that reaction. Recursion then allows the concentration of each possible aggregate to be written as a function of the free monomer concentrations. Because a monomer should be present either as a free monomer or in one of the aggregates, a set of n equations can be formed with the n free monomer concentrations as the only unknowns. This set of mass-balance equations can then be solved numerically, yielding the free monomer concentrations, from which the complete system can be reconstituted. By a step-by-step extension of the model for the aggregation of a single monomer type to include the formation of multiple aggregate types and the coassembly of multiple monomer types, we can capture increasingly complex supramolecular (co)polymerizations. In each step we illustrate how the extended model explains in detail another of the experimentally observed nonlinear effects, with the common denominator that small differences in association energies are intricately amplified at the supramolecular level. We finally arrive at our latest and most general approach to modeling (cooperative) supramolecular (co)polymerization, which encompasses all of our earlier models and shows great promise to help rationalize also future systems featuring ever-increasing complexity.

Lattice theory for binding of linear polymers to a solid substrate from polymer melts: I. Influence of chain connectivity on molecular binding and adsorption

Most theories of the binding of molecules to surfaces or for the association between molecules treat the binding species as structureless entities and neglect their rigidity and the changes in their stiffness induced by the binding process. The binding species are also taken to be "ideal," meaning that the existence of van der Waals interactions and changes in these interactions upon molecular binding are also neglected. An understanding of the thermodynamics of these multifunctional molecular binding processes has recently come into focus in the context of the molecular binding of complex molecules, such as dendrimers and DNA grafted nanoparticles, to surfaces where the degree of binding cooperativity and selectivity, as well as the location of the binding transition, are found to be sensitive to the number of binding units constrained to a larger scale polymeric scaffold. We address the fundamental problem of molecular binding by extending classical Langmuir theory to describe the particular example of the reversible binding of semiflexible polymer chains to a solid substrate under melt conditions. The polymer chains are assumed to have a variable number N of binding units (segments) and to exhibit variable bending energies and van der Waals interactions in the bulk and on the surface, in addition to strong directional interactions with the surface. The resulting generalized Langmuir theory is applied to the examination of the influence of the chain connectivity of ideal polymers on the surface coverage Θ, transition binding temperature T_1/2 at which Θ = 1/2, and on the derivative |dΘ/dT|_T=T1/2 and the constant volume specific heat of binding, C_v ^bind, measures of the cooperativity and "sharpness" of the binding transition, respectively. Paper II is devoted to the impact of the van der Waals attractive interactions and chain stiffness on the reversible binding of nonideal polymer chains to a solid surface, including the enthalpy-entropy compensation phenomenon observed experimentally in many molecular and particle binding processes.

Equilibrium Model for Supramolecular Copolymerizations

The coassembly of different building blocks into supramolecular copolymers provides a promising avenue to control their properties and to thereby expand the potential of supramolecular polymers in applications. However, contrary to covalent copolymerization which nowadays can be well controlled, the control over sequence, polymer length, and morphology in supramolecular copolymers is to date less developed, and their structures are more determined by the delicate balance in binding free energies between the distinct building blocks than by kinetics. Consequently, to rationalize the structures of supramolecular copolymers, a thorough understanding of their thermodynamic behavior is needed. Though this is well established for single-component assemblies and over the past years several models have been proposed for specific copolymerization cases, a generally applicable model for supramolecular cooperative copolymers is still lacking. Here, we provide a generalization of our earlier mass-balance models for supramolecular copolymerizations that encompasses all our earlier models. In this model, the binding free energies of each pair of monomer types in each aggregate type can be set independently. We provide scripts to solve the model numerically for any (co)polymerization of one or two types of monomer into an arbitrary number of distinct aggregate types. We illustrate the applicability of the model on data from literature as well as on new experimental data of triarylamine triamide-based copolymers in three distinct solvents. We show that apart from common properties such as the degree of polymerization and length distributions, our approach also allows us to investigate properties such as the copolymer microstructure, that is, the internal ordering of monomers within the copolymers. Moreover, we show that in some cases, also intriguing analytical approximations can be derived from the mass balances.

Lattice theory of competitive binding: Influence of van der Waals interactions on molecular binding and adsorption to a solid substrate from binary liquid mixtures

The reversible binding of molecules to surfaces is one of the most fundamental processes in condensed fluids, with obvious applications in the molecular separation of materials, chromatographic characterization, and material processing. Motivated in particular by the ubiquitous occurrence of binding processes in molecular biology and self-assembly, we have developed a lattice type theory of competitive molecular binding to solid substrates from binary mixtures of two small molecule liquids that interact between themselves by van der Waals forces in addition to exhibiting binding interactions with the solid surface. The derived theory, in contrast to previously existing theoretical frameworks, enables us to investigate the influence of van der Waals interactions on interfacial binding and selective molecular adsorption. For reference, the classic Langmuir theory of adsorption is recovered when all van der Waals interaction energies between the molecules in the bulk liquid phase and those on the surface are formally set to zero. Illustrative calculations are performed for the binding of molecules to a solid surface from pure liquids and from their binary mixtures. The properties analyzed include the surface coverage θ, the binding transition temperature T_bind, the individual surface coverages, θ_A and θ_C, and the relative surface coverages, σ_AC≡θ_A/θ_C or σ_CA≡θ_C/θ_A. The latter two quantities coincide with the degrees of adsorption directly determined from experimental adsorption measurements. The Langmuir theory is shown to apply formally under a wide range of conditions where the original enthalpies (Δh or Δh_A and Δh_C) and entropies (Δs or Δs_A and Δs_C) of the binding reactions are simply replaced by their respective "effective" counterparts (Δh_eff or Δh_A^eff and Δh_C^eff and Δs_eff or Δs_A^eff and Δs_C^eff), whose values depend on the strength of der Waals interactions and of the "bare" free energy parameters (Δh or Δh_A and Δh_C, and Δs or Δs_A and Δs_C). Numerous instances of entropy-enthalpy compensation between these effective free energy parameters follow from our calculations, confirming previous reports on this phenomenon obtained from experimental studies of molecular binding processes in solution.

Macrofibril Formation

Macrofibrils are the main structural component of the hair cortex, and are a composite material in which trichokeratin intermediate filaments (IFs) are arranged as organised arrays embedded in a matrix composed of keratin-associated proteins (KAPs) and keratin head groups. Various architecture of macrofibrils is possible, with many having a central core around which IFs are helically arranged, an organisation most accurately described as a double-twist arrangement. In this chapter we describe the architecture of macrofibrils and then cover their formation, with most of the material focusing on the theory that the initial stages of macrofibril formation are as liquid crystals.

Mixtures of two self- and mutually-associating liquids: Phase behavior, second virial coefficients, and entropy-enthalpy compensation in the free energy of mixing

The theoretical description of the phase behavior of polymers dissolved in binary mixtures of water and other miscible solvents is greatly complicated by the self- and mutual-association of the solvent molecules. As a first step in treating these complex and widely encountered solutions, we have developed an extension of Flory-Huggins theory to describe mixtures of two self- and mutually-associating fluids comprised of small molecules. Analytic expressions are derived here for basic thermodynamic properties of these fluid mixtures, including the spinodal phase boundaries, the second osmotic virial coefficients, and the enthalpy and entropy of mixing these associating solvents. Mixtures of this kind are found to exhibit characteristic closed loop phase boundaries and entropy-enthalpy compensation for the free energy of mixing in the low temperature regime where the liquid components are miscible. As discussed by Widom et al. [Phys. Chem. Chem. Phys. 5, 3085 (2003)], these basic miscibility trends, quite distinct from those observed in non-associating solvents, are defining phenomenological characteristics of the "hydrophobic effect." We find that our theory of mixtures of associating fluids captures at least some of the thermodynamic features of real aqueous mixtures.

Supramolecular Copolymers: Structure and Composition Revealed by Theoretical Modeling

Supramolecular copolymers, non-covalent analogues of synthetic copolymers, constitute a new and promising class of polymers. In contrast to their covalent counterparts, the details of their mechanism of formation, as well as the factors determining their composition and length, are still poorly understood. Here, the supramolecular copolymerization between two slightly structurally different benzene-1,3,5-tricarboxamide (BTA) monomers functionalized with either oligodimethylsiloxane (oDMSi) or alkyl side chains is unraveled by combining experimental and theoretical approaches. By applying the "sergeant-and-soldiers" approach using circular dichroism (CD) experiments, we are able to obtain detailed insights into the structure and composition of these supramolecular copolymers. Moreover, we observe an unexpected chiral induction upon mixing two independently CD-silent solutions of the achiral (soldier) and chiral (sergeant) monomers. We find that the subtle differences in the chemical structure of the two monomers impact their homopolymerization mechanism: whereas alkyl-BTAs cooperatively self-assemble, oDMSi-BTAs self-assemble in an isodesmic manner. The effect of these mechanistic differences in the supramolecular copolymerization process is investigated as a function of the composition of the two monomers and explicitly rationalized by mathematical modeling. The results show that, at low fractions of oDMSi-BTA sergeants (<10 mol%), the polymerization process is cooperative and the supramolecular helicity is biased toward the helical preference of the sergeant. However, at higher fractions of oDMSi-BTA sergeant (>25 mol%), the isodesmic assembly of the increasing amounts of sergeant becomes more dominant, and different species start to coexist in the copolymerization process. The analysis of the experimental data with a newly developed theoretical model allows us to quantify the thermodynamic parameters, the distribution of different species, and the compositions and stack lengths of the formed supramolecular copolymers existing at various feed ratios of the two monomers.

Fragmentation and Coagulation in Supramolecular (Co)polymerization Kinetics

The self-assembly of molecular building blocks into one-dimensional supramolecular architectures has opened up new frontiers in materials science. Due to the noncovalent interactions between the monomeric units, these architectures are intrinsically dynamic, and understanding their kinetic driving forces is key to rationally programming their morphology and function. To understand the self-assembly dynamics of supramolecular polymerizations (SP), kinetic models based on aggregate growth by sequential monomer association and dissociation have been analyzed. However, fragmentation and coagulation events can also play a role, as evident from studies on peptide self-assembly and the fact that aggregations can be sensitive to mechanical agitations. Here, we analyze how fragmentation and coagulation events influence SP kinetics by theoretical analysis of self-assembling systems of increasing complexity. Our analysis starts with single-component systems in which aggregates are able to grow via an isodesmic or cooperative nucleation-elongation mechanism. Subsequently, equilibration dynamics in cooperative two-component supramolecular copolymerizations are investigated. In the final part, we reveal how aggregate growth in the presence of competing, kinetically controlled pathways is influenced by fragmentation and coagulation reactions and reveal how seed-induced growth can give rise to block copolymers. Our analysis shows how fragmentation and coagulation reactions are able to modulate SP kinetics in ways that are highly system dependent.

Coupling of isotropic and directional interactions and its effect on phase separation and self-assembly

The interactions of molecules and particles in solution often involve an interplay between isotropic and highly directional interactions that lead to a mutual coupling of phase separation and self-assembly. This situation arises, for example, in proteins interacting through hydrophobic and charged patch regions on their surface and in nanoparticles with grafted polymer chains, such as DNA. As a minimal model of complex fluids exhibiting this interaction coupling, we investigate spherical particles having an isotropic interaction and a constellation of five attractive patches on the particle's surface. Monte Carlo simulations and mean-field calculations of the phase boundaries of this model depend strongly on the relative strength of the isotropic and patch potentials, where we surprisingly find that analytic mean-field predictions become increasingly accurate as the directional interactions become increasingly predominant. We quantitatively account for this effect by noting that the effective interaction range increases with increasing relative directional to isotropic interaction strength. We also identify thermodynamic transition lines associated with self-assembly, extract the entropy and energy of association, and characterize the resulting cluster properties obtained from simulations using percolation scaling theory and Flory-Stockmayer mean-field theory. We find that the fractal dimension and cluster size distribution are consistent with those of lattice animals, i.e., randomly branched polymers swollen by excluded volume interactions. We also identify a universal functional form for the average molecular weight and a nearly universal functional form for a scaling parameter characterizing the cluster size distribution. Since the formation of branched clusters at equilibrium is a common phenomenon in nature, we detail how our analysis can be used in experimental characterization of such associating fluids.

String model for the dynamics of glass-forming liquids

We test the applicability of a living polymerization theory to describe cooperative string-like particle rearrangement clusters (strings) observed in simulations of a coarse-grained polymer melt. The theory quantitatively describes the interrelation between the average string length L, configurational entropy Sconf, and the order parameter for string assembly Φ without free parameters. Combining this theory with the Adam-Gibbs model allows us to predict the relaxation time τ in a lower temperature T range than accessible by current simulations. In particular, the combined theories suggest a return to Arrhenius behavior near Tg and a low T residual entropy, thus avoiding a Kauzmann "entropy crisis."

Cooperative polymerization of one-patch colloids

We numerically investigate cooperative polymerization in an off-lattice model based on a pairwise additive potential using particles with a single attractive patch that covers 30% of the colloid surface. Upon cooling, these particles self-assemble into small clusters which, below a density-dependent temperature, spontaneously reorganize into long straight tubes. We evaluate the partition functions of clusters of all sizes to provide an accurate description of the chemical reaction constants governing this process. Our calculations show that, for intermediate sizes, the partition functions retain contributions from two different structures, differing in both energy and entropy. We illustrate the microscopic mechanism behind the complex polymerization process in this system and provide a detailed evaluation of its thermodynamics.

Phase field method for nonequilibrium dynamics of reversible self-assembly systems

Phase field methods are extended to describe the nonequilibrium dynamics of reversible self-assembly systems, an extension that is complicated by the mutual coupling of many non-conserved order parameters into a set of highly nonlinear partial differential equations. Further complications arise because the sum of all non-conserved order parameters equals a conserved order parameter. The theory is developed for the simplest model of reversible self-assembly in which no additional constraints are imposed on the self-assembly process since the extension to treat more complex self-assembly models is straightforward. Specific calculations focus on the time evolution of the cluster size distribution for a free association system that is rapidly dropped from one ordered state to a more ordered state within the one-phase region. The dynamics proceed as expected, thereby providing validation of the theory which is also capable of treating systems with spatial inhomogeneities.

Lattice cluster theory of associating telechelic polymers. III. Order parameter and average degree of self-assembly, transition temperature, and specific heat

The lattice cluster theory of strongly interacting, structured polymer fluids is applied to determine the thermodynamic properties of solutions of telechelic polymers that may associate through bifunctional end groups. Hence, this model represents a significant albeit natural extension of a diverse array of prior popular equilibrium polymerization models in which structureless "bead" monomers associate into chain-like clusters under equilibrium conditions. In particular, the thermodynamic description of the self-assembly of linear telechelic chains in small molecule solvents (initiated in Paper II) is systematically extended through calculations of the order parameter Φ and average degree <N> of self-assembly, the self-assembly transition temperature T(p), and the specific heat C(V) of solutions of telechelic molecules. Special focus is placed on examining how molecular and thermodynamic parameters, such as the solution composition φ, temperature T, microscopic interaction energies (ε(s) and ε), and length M of individual telechelic chains, influence the computed thermodynamic quantities that are commonly used to characterize self-assembling systems.

An equilibrium model for chiral amplification in supramolecular polymers

We describe a model that rationalizes amplification of chirality in cooperative supramolecular copolymerization. The model extends nucleation-elongation based equilibrium models for growth of supramolecular homopolymers to the case of two monomer and aggregate types. Using the principle of mass-balance for the two monomer types, we derive a set of two nonlinear equations, describing the thermodynamic equilibrium state of the system. These equations can be solved by numerical methods, but also analytical approximations are derived. The equilibrium model allows two-sided growth of the aggregates and can be applied to symmetric supramolecular copolymerizations, corresponding to the situation in which the monomers are enantiomerically related, as well as to the more general case of nonsymmetric supramolecular copolymerizations. In detail, so-called majority-rules phenomena in supramolecular systems with isodesmic as well as cooperative growth are analyzed. Comparison of model predictions with experimental data shows that the model gives a very good description of both titration and melting curves. When the system shows cooperative growth, the model leads to a phase diagram in which the presence of the various aggregate types is given as a function of composition and temperature.

Stretched Exponential Stress Relaxation in a Thermally Reversible, Physically Associating Block Copolymer Solution

The shear stress relaxation of a thermally reversible, physically associating solution formed from a triblock copolymer in solvent selective for the mid-block was found to be well described over a broad temperature range by a stretched exponential function with a temperature independent ‘stretching exponent’, β ≈ 1/3. This same exponent value has been suggested to have particular significance in describing structural relaxation in a wide range of disordered viscoelastic materials ranging from associating polymer materials (‘gels’) to glass-forming liquids. We quantify the temperature dependence of the high frequency, or short time, shear modulus as function of temperature and find that this property also follows a variation often observed in gels and glass-forming materials.

Theoretical models of nonlinear effects in two-component cooperative supramolecular copolymerizations

The understanding of multi-component mixtures of self-assembling molecules under thermodynamic equilibrium can only be advanced by a combined experimental and theoretical approach. In such systems, small differences in association energy between the various components can be significantly amplified at the supramolecular level via intricate nonlinear effects. Here we report a theoretical investigation of two-component, self-assembling systems in order to rationalize chiral amplification in cooperative supramolecular copolymerizations. Unlike previous models based on theories developed for covalent polymers, the models presented here take into account the equilibrium between the monomer pool and supramolecular polymers, and the cooperative growth of the latter. Using two distinct methodologies, that is, solving mass-balance equations and stochastic simulation, we show that monomer exchange accounts for numerous unexplained observations in chiral amplification in supramolecular copolymerization. In analogy with asymmetric catalysis, amplification of chirality in supramolecular polymers results in an asymmetric depletion of the enantiomerically related monomer pool.

Elasticity of networks with permanent and thermoreversible cross-links

Simplified models of flexible chain and stiff fiber networks are introduced to address how the network elasticity becomes modified when the cross-linking is thermoreversible in nature and changes in the stability of the network with deformation. These idealized models apparently able to capture many aspects of the elastic properties of real networks.

Atomic force microscopy nanomanipulation of shape persistent, spherical, self-assembled aggregates of gold nanoparticles

Gold (Au) nanoparticles have been synthesized that are stabilized by an organic ligand bearing a dithiolane functional group for binding to Au, an oligo(p-phenylene vinylene) (OPV) chromophoric group to drive self-assembly via π-π interactions, and a hydroxy functionality for interparticle hydrogen bonding. The OPV-Au particles reversibly self-assemble in n-heptane solution, forming shape persistent, spherical, nanometer-sized aggregates that do not collapse on a substrate. Optical absorption and transmission electron microscopy tomography studies show that the size and shape persistency can be tuned by modification of the ligands, adjustment of the core size, and variation of the concentration. The spherical assemblies can be manipulated with the tip of an atomic force microscope: an aggregate can be pushed over the surface for at least 20 times with nanometer precision and without substantial loss of material.

Self-consistent field lattice model study on the phase behavior of physically associating polymer solutions

The phase behavior of physically associating polymer solutions, where the polymer chain contains a small fraction of "stickers" regularly placed along the backbone, is studied using self-consistent field lattice model. Two inhomogenous morphologies are observed. One is a microfluctuation homogenous (MFH) morphology, where the mean-field values of the local average concentrations of polymers phi(P)(r) and stickers phi(st)(r) slightly fluctuate around their respective bulk average values phi(P) and phi(st) and regularly from site to site. The other is a randomly close-packed micelle (RCPM) morphology. The structure of the micelle in RCPM morphology is similar to that of the "flower micelle" in the telechelic associative polymer system, where stickers are located in the core of the micelle and nonsticky groups in the corona. When phi(P) approximately or > 0.08, if homogenous associating polymer solutions are cooled, MFH morphology appears, and the system entirely changes from homogenous solutions (HS) to MFH morphology; If the solutions are cooled further, RCPM morphology appears. When phi(P) < 0.08, however, RCPM morphology appears immediately. If phi(P) < 0.53, a macroscopic phase separation, where the polymer rich phase is RCPM morphology, occurs. If phi(P) approximately or > 0.53, only RCPM morphology is found in the system. A peak appears in the temperature-dependent specific-heat curve C(V)(chi) at each transition point. For the HS-MFH transition, C(V)(chi) has an abrupt increase and a slow decrease, whereas for the MFH-RCPM transition, both the increase and the decrease in C(V)(chi) are slow. Furthermore, the system with only MFH morphology may be trapped in one of the two energy basins in a experimental time scale. However, the appearance of RCPM morphology means that the system is trapped in one of a series of "deeper" energy basins, and it is very difficult to jump off this deep basin into the one of MFH morphology or one of the other RCPM morphologies through thermal fluctuations.

Development of minimal models of the elastic properties of flexible and stiff polymer networks with permanent and thermoreversible cross-links

We review the elasticity of flexible and stiff polymer networks with permanent cross-links and synthesize these results into a unifying polymer chain network model. This framework is then used to address how the network elasticity becomes modified when the network cross-linking is thermoreversible in nature, changes in the stability of the network with deformation, and the effect of a variable rate of network deformation on the non-linear elastic response. Comparisons are made between this class of simplified network models with elasticity measurements performed on flexible chain and stiff fiber networks, both with permanent and associative cross-links. Although these network models are highly idealized, they are apparently able to capture many aspects of the elastic properties of diverse real networks.

Crowding induced self-assembly and enthalpy-entropy compensation

We develop a general virial expansion to describe the influence of molecular additives on the equilibrium self-assembly of proteins or other supermolecularly assembling species M in solution. When specialized to high molar mass polymer additives, the cross-virial coefficient between the polymer and M, which dominates this effect, is found to vanish at a particular temperature T_{Theta} corresponding to an enthalpy-entropy compensation condition. Specifically, the increased stability of the assembled form of M, due to the modification of the entropy of the assembly by repulsive polymer-protein interactions, is progressively compensated by attractive interactions that alter the enthalpy of assembly.

Theoretical issues relating to thermally reversible gelation by supermolecular fiber formation

Existing models of the thermodynamics and dynamics of self-assembly are summarized to provide a context for discussing the difficulties that arise in modeling supermolecular fiber assembly and the formation of thermally reversible gels through fiber growth and branching. Challenging problems in this field, such as the physical origin of fibers of uniform diameter and fiber twisting, the kinetics of fiber growth, the hierarchical bundling of fibers into "superfibers", fiber branching, gelation through fiber impingement and the associated phenomenon of fractal fiber network and spherulite formation, and the origin and control of structural polymorphism in the fiber and superfiber geometry, are discussed from a personal perspective. Suggestions are made for integrating current research efforts into a more coherent multiscale description of fiber formation and gelation on molecular, mesoscopic, and macroscopic scales.

Equilibrium polymerization models of re-entrant self-assembly

As is well known, liquid-liquid phase separation can occur either upon heating or cooling, corresponding to lower and upper critical solution phase boundaries, respectively. Likewise, self-assembly transitions from a monomeric state to an organized polymeric state can proceed either upon increasing or decreasing temperature, and the concentration dependent ordering temperature is correspondingly called the "floor" or "ceiling" temperature. Motivated by the fact that some phase separating systems exhibit closed loop phase boundaries with two critical points, the present paper analyzes self-assembly analogs of re-entrant phase separation, i.e., re-entrant self-assembly. In particular, re-entrant self-assembly transitions are demonstrated to arise in thermally activated equilibrium self-assembling systems, when thermal activation is more favorable than chain propagation, and in equilibrium self-assembly near an adsorbing boundary where strong competition exists between adsorption and self-assembly. Apparently, the competition between interactions or equilibria generally underlies re-entrant behavior in both liquid-liquid phase separation and self-assembly transitions.