Eﬀect of sticker clustering on the dynamics of associative networks † Soft Matter

Recent experimental and theoretical work has shown that sticker clustering can be used to enhance properties such as toughness and creep resistance of polymer networks. While it is clear that the changes in properties are related to a change in network topology, the mechanistic relationship is still not well understood. In this work, the eﬀect of sticker clustering was investigated by comparing the dynamics of random copolymers with those where the stickers are clustered at the ends of the chain in the unentangled regime using both linear mechanics and diﬀusion measurements. Copolymers of N , N dimethyl acrylamide (DMA) and pendant histidine groups were synthesized using reversible addition– fragmentation chain transfer (RAFT) polymerization. The clustered polymers were synthesized using a bifunctional RAFT agent, such that the midblock consisted of PDMA and the two end blocks were random copolymers of DMA and the histidine-functionalized monomer. Upon addition of Ni ions, transient metal-coordinate crosslinks are formed as histidine–Ni complexes. Combined studies of rheology, neutron scattering and self-diffusion measurements using forced Rayleigh scattering revealed changes to the network topology and stress relaxation modes. The network topology is proposed to consist of aggregates of the histidine–Ni complexes bridged by the non-associative midblock. Therefore, stress relaxation requires the cooperative dissociation of multiple bonds, resulting in increased relaxation times. The increased relaxation times, however, were accompanied by faster diffusion. This is attributed to the presence of defects such as elastically inactive chain loops. This study demonstrates that the effects of cooperative sticker dissociation can be observed even in the presence of a significant fraction of loop defects which are known to alter the nonlinear properties of conventional telechelic polymers.


Introduction
Associative networks are ubiquitous both in natural and synthetic materials, and the dynamics within these networks dictate many of their desirable properties such as self-healing and stress relaxation. 1,2 The versatility of these materials is the result of the wide array of design options, ranging from different chain architectures (linear vs. star) 3,4 to binding chemistry 5,6 to solvent environments 7 to junction functionality. 8,9 One strategy that has garnered interest of late is the clustering of stickers such that they are concentrated at the chain ends. [10][11][12] This results in the chain ends participating in multivalent interactions due to the proximity of the stickers to one another. In biological systems, the presence of multiple binding sites in the form of multivalent ligands can enhance the strength and specificity of interactions compared to weak binding affinities of monovalent ligands. 13,14 For example, protein-carbohydrate recognition events mediate processes such as pathogen-cell adhesion and inflammatory response, and carbohydrate epitopes are often present as multivalent arrays at the cell membrane to serve as highly efficient ligands. In synthetic systems, sticker clustering has been demonstrated to be a promising approach for designing networks with improved toughness, dissipation and creep resistance. 10,15 However, the molecular mechanism for these improvements is not fully understood.
As an accessible model of clustered systems, several studies have now been performed on networks formed by chains with a triblock architecture of a soluble midblock with two end blocks that are random copolymers of the stickers and diluent monomer. In comparison to random copolymers of the same molecular weight and number of stickers per chain, a delay in the terminal relaxation time was reported, accompanied by a higher activation energy for the onset of flow for the clustered copolymers. 10,11,16 This delay has been attributed to the need for cooperative dissociation of multiple bonds before the chains are able to relax. 15,17 The higher activation energy indicates that these networks have a stronger temperature dependence, and this has been suggested to originate from the difference in the formation of bonds as the temperature is reduced. 17 While lowering the temperature should drive the system to favor bond association in both random and clustered sticker configurations, the proximity of the stickers in the clustered polymers enhances this effect. In some systems, the triblock chain architecture leads to microphase separation and results in the formation of structures such as cylinders 16 and lamellae. 18 However, changes in the viscoelastic properties were observed even for networks that did not undergo microphase separation, 10 which indicates that the physics behind these changes in mechanical properties are not solely the result of microphase separation.
Self-diffusion measurements using forced Rayleigh scattering (FRS) provide an orthogonal probe to mechanical property measurements that can provide insight into the dynamics of associative polymer networks. Studies of several unentangled networks 3,5,6 have revealed further details of the dynamics of these gels over length scales of several time the radius of gyration, R g of the polymer. The previously investigated systems include a protein gel with pentavalent coiled-coiled associations, 5 a four-arm star polymer end-functionalized with terpyridine and crosslinked with Zn 2+ in DMF 6 and linear random copolymers with pendant histidine groups crosslinked with Ni 2+ in water. 3 In all these systems, superdiffusive scaling was observed at the smaller range of length scales that is experimentally accessible, prior to transitioning to a regime with Fickian scaling. The observation of superdiffusive scaling in these networks was attributed to the presence of two diffusive modes with distinct diffusivities, which are walking and hopping, in the molecular model developed by Ramirez et al. 19 Walking refers to diffusive modes where motion of the chains require sequential dissociations and reassociations of the stickers while hopping refers to diffusive modes where the chains dissociate all of its stickers to diffuse over several times the R g . These dynamics were not detectable from studies of the polymers' viscoelastic properties and provide further insights into the dynamics of associative networks.
In this work, the dynamics and mechanics of a model associative network was investigated using rheology, smallangle neutron scattering (SANS) and FRS to provide a more detailed picture of the effect of sticker clustering. The model polymer system consists of random copolymers of N,Ndimethyl acrylamide (DMA) and a histidine-functionalized monomer, synthesized using RAFT polymerization. The clustered polymer was synthesized using a bifunctional RAFT agent, such that the midblock consists of only PDMA while the end blocks are random copolymers of DMA and the histidine-functionalized monomer ( Fig. 1(A)). On addition of Ni 2+ , the crosslinks are formed as histidine-Ni complexes. The histidine-Ni complex was chosen because its kinetics have been thoroughly characterized, 20 making it particularly suited to elucidate the effect of sticker clustering. This publication compares this new data on clustered polymers to those with randomly distributed stickers, published in a previous work. 3 The effect of sticker clustering on the network structure was probed through SANS experiments, while network stiffness and terminal relaxation were characterized using rheology. To complement these measurements, self-diffusion within the network was probed using FRS, a technique which has not been applied in earlier studies investigating the effect of sticker clustering. 10,11 Comparisons of the results to existing theories provide insights into the molecular mechanism behind the observed behaviors.

Materials
Fluorescein-5-maleimide was purchased from ThermoFisher Scientific. N-(3-N-Boc-Nim-trityl-N-3-methacrylamidopropyl-Lhistidinamide (HisMA) 20 and 2-(ethylthio-carbonothioylthio)-2-methylpropionic acid (EMP) 21 were synthesized following published procedures. N,N-Dimethyl-acrylamide (DMA) was purified through a basic alumina column to remove inhibitor before polymerization. All other chemical reagents were purchased from Sigma-Aldrich or VWR and used as received. Polymers PDHM5 and PDHM10 were synthesized for a previous publication, 3 and the same polymers were used in this study. All data for PDHM10 was taken from the previous study, while some additional data was collected on PDHM5 for this study.

Characterizations
Gel permeation chromatography (GPC) measurements were performed on an Agilent 1260 LC system with two ResiPore columns (300 Â 7.5 mm, Agilent Technologies, Santa Clara, CA) in series at a flow rate of 1 mL min À1 at 70 1C, where DMF with 0.02 M LiBr was used as the mobile phase. The molecular weights were determined using a Wyatt miniDAWN TREOS multiangle light scattering detector and a Wyatt Optilab T-rEX differential refractive index detector. Liquid chromatrographymass spectrometry (LC-MS) analysis was performed using an Agilent 1260 Infinity LC system coupled with a 6130 quadrupole mass spectrometer. A mixture of 0.1% formic acid in water and MeCN was used as the mobile phase. NMR spectra were recorded on a Varian 400 MHz spectrometer. The residual undeuterated solvent peaks were used as references (7.27 ppm for CDCl 3 and 4.79 ppm for D 2 O).

Synthesis of EMP dimer
In a Schlenk flask, EMP (0.763 g, 3.4 mmol) was dissolved in 10 mL of anhydrous DMF. Hexafluorophosphate azabenzotriazole tetramethyl uronium (HATU) (1.29 g, 3.4 mmol) and 4-methylmorpholine (0.5 g, 5.0 mmol) were added, and the solution was stirred, under nitrogen atmosphere, for 30 minutes. Ethylene diamine (0.11 mL, 1.7 mmol) was added slowly, and the reaction was stirred overnight (Scheme 1). Then, the solvent was removed under vacuum, and the crude product was dissolved in 50 mL of DCM. The organic solution was washed with deionized water (3 Â 50 mL), washed with brine (3 Â 50 mL), and dried over Na 2 SO 4 , The product was then filtered and concentrated under vacuum. The crude product was purified via column chromatography with silica gel as the stationary phase and hexanes and ethyl acetate (70 : 30 to 40 : 60) as eluents. The product was obtained as a yellow solid in 89% yield. 1  Synthesis of PDMA polymers with clustered pendant histidine side groups (PDHMc8). Clustered copolymers from DMA and HisMA were synthesized by reversible addition-fragmentation chain transfer (RAFT) polymerization (Scheme 2). The total monomer concentration in polymerization was 2.0 M, and the ratio of DMA/HisMA/EMP dimer/AIBN was 333 : 8 : 1 : 0.2. This monomer pair was chosen because the reactivity ratios for acrylamide and methacrylamide have been shown to be close to unity. 22,23 Therefore, the two monomers are expected to copolymerize in a nearly statistical manner, with HisMA distributed evenly in the end blocks. The polymerization was performed in MeCN at 60 1C for 7 h. In the first stage of the synthesis, DMA was polymerized until a molar mass of 15.9 kg mol À1 was achieved, as determined by DMF GPC. Then HisMA, dissolved in MeCN, was cannulated into the reaction vial. Once the desired conversion was achieved, the reaction was quenched by exposure to air and cooling to room temperature. The polymer was purified by precipitation into diethyl ether once and dried under vacuum. The mole fraction of HisMA in the polymer was determined to be 2.9 mol% by 1 H NMR (Fig. S1, ESI †), close to the feed composition of 2.7 mol%. The molar mass of polymer was characterized by DMF GPC prior to the deprotection step and was determined to be 29.5 kg mol À1 (Fig. S2, ESI †). To remove the Boc and Trt protecting groups, the resulting polymers (700 mg) were dissolved in DCM (11.7 mL). Water (291.7 mL), triisopropylsilane (TIPS, 291.7 mL), and trifluoracetic acid (TFA, 11.7 mL) were sequentially added to the solution. The mixture was stirred at room temperature for 2 h. The volatiles were then removed under vacuum, and the residue was dissolved in MeOH. The polymers were recovered by precipitation into diethyl ether twice. The polymers were then dissolved in water, transferred to a centrifugal filter (3 kDa MWCO), and spun at 4000 Â g for 1 h. More water was then added, and the filtration was repeated four times. The polymers were then filtered through a 0.45 mm filter and lyophilized to yield 550.6 mg of product. Complete removal of the Boc and Trt groups was evidenced by 1 H NMR (Fig. S3, ESI †).

Synthesis of fluorescein-labeled PDHMc8 polymers
The deprotected polymers (25 mg, 0.9 mmol) were first dissolved in 1.25 mL DMF, and then hexylamine (4.79 mL, 36 mmol) was added. The reaction was stirred under a nitrogen atmosphere overnight to minimize undesirable cysteine oxidation and to ensure complete aminolysis. Then tris(2-carboxyethyl)phosphine hydrochloride (TCEPÁHCl, 5.2 mg, 18 mmol) and maleimide-functionalized fluorescein (7.74 mg in 180 mL of DMSO, 18 mmol) were added to the reaction mixture. After the reaction was stirred overnight in the dark, the solution was diluted with 30 mL 5% DMSO in water. The mixture was transferred to a centrifugal filter (3 kDa MWCO), spun at 4000 Â g for 1 h at 4 1C, and further diluted 30 mL with 5% DMSO in water. This process was repeated several times until the spin-through fraction was colorless. A final spin with water then removed the DMSO. Polymers were filtered through a 0.45 mm filter and lyophilized to yield 13.9 mg of final product.

Gel preparation
The gels of the PDMA polymers with pendant histidine groups were prepared following previously published procedures. 3,20 The polymers were first dissolved in a Bis-Tris buffer (100 mM, pH 7.0), and complete dissolution was confirmed when the solution appeared clear. The appropriate volume of a stock solution containing 200 mM NiCl 2 and 100 mM Bis-Tris was then added, and the mixture was vortexed for 15 s. The volume of NiCl 2 stock solution required to prepare the gels at 2 : 1 of histidine : Ni was determined by calculating the concentration of histidine in the polymer using 1 H NMR (Fig. S1, ESI †). Once mixed, the appropriate volume of a stock solution of 1 M NaOH with 100 mM Bis-Tris buffer was then added to adjust the pH to 7.0. The volume of NaOH stock solution required to adjust the pH to 7.0 was determined by titration experiments in dilute solution ( Fig. S4, ESI †). The gels were then mixed with a micro spatula until a macroscopically homogenous gel was obtained, and the gels were centrifuged at 21 100 Â g to remove air bubbles introduced during mixing.

Small-angle neutron scattering (SANS)
SANS experiments were conducted at the Oak Ridge National Laboratory Spallation Neutron Source (ORNL SNS) Extended Q-Range SANS Diffractometer. The beam aperture was 8 mm. To facilitate neutron scattering contrast in SANS experiments, Bis-Tris buffered D 2 O (75 mM Bis-Tris and 25 mM HCl in D 2 O) was prepared by gravimetric measurements of the buffer components to reach pD 7.0. Additionally, the 200 mM stock solution of NiCl 2 was prepared by dissolving 252.9 mg of NiCl 2 in 10 mL of the Bis-Tris buffered D 2 O with pD 7.0. The 1 M NaOD stock solution was prepared by diluting 30% (w/w) solution of NaOD in D 2 O with the Bis-Tris buffer with pD 7.0. All stock solutions were filtered using 0.2 mm Acrodisc syringe filters (PALL Corporation). The gels were prepared as previously described and pressed between two quartz disks with a Teflon spacer (1 mm thickness, 13 mm inner diameter, 17 mm outer diamater). The quartz sample sandwich was then loaded in a titanium cell. Scattering patterns were measured using sample-todetector distances of 2.5, 4.0 and 9.0 m, with neutron wavelength bands of 2.0-6.0 Å, 4.0-7.5 Å and 15.0-17.5 Å, respectively. This corresponds to a Q-range of 0.02-10 nm À1 . All experiments were performed at 25 1C. The raw scattering intensity was reduced using the Mantid reduction package 24 and corrected for the background by subtracting the scattering from an empty sample cell and a sample cell containing buffered D 2 O. The absolute intensity was calibrated using a porous silica standard sample. The reduced SANS curves were fit using non-linear least squares regression to a correlation length model. [25][26][27][28] Rheology Frequency sweep experiments were performed on an Anton Paar 301 Physica rheometer, using a stainless-steel cone-plate upper geometry (25 mm in diameter, 11 angle). Inertial calibration and motor adjustment were performed before each measurement. All hydrogel samples were centrifuged at 21 100 Â g for 10 min at 4 1C to remove bubbles before loading onto the rheometer. Dehydration was minimized by adding mineral oil to the sample edge. Experiments were performed at four temperatures: 5, 15, 25, and 35 1C. The temperature was controlled by a Peltier plate. Time-temperature superposition was used to construct master curves, and the procedure is described in Section D (p. 8) of the ESI. † Experiments were performed at 1% strain, which was within the linear viscoelastic (LVE) region as determined by strain sweep experiments.

Forced Rayleigh scattering
For self-diffusion measurements, 80 mM of the fluoresceinlabelled polymers was added before the addition of NiCl 2 stock solution during the gel preparation. This allows for thorough mixing of the fluorescein-labelled polymers into the solution before the formation of the gel. All samples were sealed between two quartz disks (17 mm in diameter) separated by a 0.2 mm thick Teflon spacer. To eliminate shear history from loading, all samples were left overnight at room temperature. Samples were equilibrated at the desired temperature for 1 h before further experiments were performed. The self-diffusion measurements were performed using forced Rayleigh scattering (FRS), as previously described. 5,29,30 Briefly, a 100 mW continuous wave laser with l = 488 nm was split into two beams which were individually refocused and crossed at an angle of y onto the sample. This generated a holographic grating of characteristic spacing d that is defined by the following equation (1) On exposure of the sample for 100-500 ms, the photochromic fluorescein dye conjugated to the tracer molecules was irreversibly isomerized in the volumes of constructive interference, producing an amplitude grating of dye concentration. Diffusion of the dye results in a sinusoidal concentration profile by diffusion, which was monitored by diffraction of a single reading beam at the same wavelength. The intensity of the reading beam was attenuated by 10 À4 so that it was low enough to ensure the change of the profile was only due to diffusion. The time constant, t, can be extracted from fitting a stretched exponential function to the signal: where I is the intensity, b is the stretched exponent ranging from 0 to 1, and B is the incoherent background. The average decay time constant was calculated as the first moment of the stretched exponential: where G is the gamma function. In FRS experiments that measure simple Fickian diffusion, the diffusion coefficient is given by:

Results and discussion
Effect of HisMA distribution on critical concentration regimes When categorizing the dynamic regime of associative polymers, the important concentration regimes to consider are the chain overlap, f overlap , strand (between stickers) overlap, f s , and entanglement concentration, f e (see ESI † for calculation). 31,32 Since the degree of polymerization, N is approximately B250 for all the polymers investigated, f overlap and f e (which depends on N) are similar for the three polymers (Table 1). In contrast, f s is proportional to the spacing between the stickers, l. Since l = N/S, f s is inversely proportional to the average number of stickers per chain, S, for the random copolymers.
For the clustered polymer, there are two spacings that must be considered, which are the midblock separating the two chains ends, l mid = 160, and the average spacing between stickers in the two end blocks, l end = 14, such that two f s can be calculated ( Table 1). The gels investigated in this work were prepared at 25% and 30% (w/v), which is well above f overlap and well below f e , such that an unentangled, percolated network is formed. Additionally, f s defines the limit above which most stickers should be in interchain bonds, such that gels prepared below this concentration are predicted to have a significant fraction of intrachain bonds. For the random copolymers, f = 25% (w/v) is above f s . For the clustered polymer, it is above f s for the midblock but below f s for the end blocks. Thus, for the clustered polymer, the chains are overlapping enough to form interchain bonds between different chains, but these bonds likely exist as aggregates bridged by the midblocks. Since the same trends were observed for the gels at 30% (w/v) (see Fig. S6-S8 in Section E of ESI †) the following analysis will focus on the gels prepared at 25% (w/v).

Effect of clustering on linear viscoelastic properties and network topology
The linear viscoelastic response of the gels indicates that the network topology and mode of stress relaxation is altered by sticker clustering. The frequency sweeps for all the gels at a concentration of 25% (w/v) show a single plateau, G p , in the storage modulus, G 0 at high angular frequencies which crosses over with the loss modulus, G 00 with crossover angular frequency, o c (Fig. 2(A-C)). While the relaxation times, t = 2p/ o c , show an increase with the number of stickers per chain (t 5 o t c8 o t 10 ), the modulus does not show this trend. Rather, the plateau modulus of the clustered sticker polymer is lower than either of the two polymers with randomly distributed stickers: G p,c8 o G p,5 o G p, 10 . The plateau modulus, G p is a measure of the concentration of elastically active strands, and the trend observed for G p can be explained by considering the relation between G p and the average spacing between crosslinks, l under the affine network assumption. For gels prepared at a concentration of f (volume fraction), where a = 1.3 nm 33,34 is the monomer size and l is the average spacing between crosslinks. 32 Fig. 3 plots the experimentally observed G p vs. the curve of predicted G p vs. l from eqn (4). When l for the clustered sticker polymer is taken to be the midblock length, there is close agreement between the theory and the experiment. This suggests that the lower G p for PDHMc8 is consistent with the midblocks of the polymers acting as the only elastically active strands in the network, which is in support of the above proposed network topology wherein midblocks serve as a bridge between end block cluster bond aggregates ( Fig. 1(B)). Further evidence of the proposed changes in network topology is provided by small angle neutron scattering (SANS) experiments. The appearance of an upturn at low q (onset shown by black single arrows in Fig. 4) is a feature often observed in disordered hydrogels, 25 as expected for hydrogels formed from polymers with associative groups along their backbone. This upturn begins at a larger q for the PDHMc8 gel which indicates the presence of inhomogeneity at smaller length scales than expected for networks made with random copolymers with the same average composition. To quantify this effect, the SANS data was fit to an empirical correlation length model developed by Hammouda et al., 35 that has been used to analyze scattering from other hydrogels. [25][26][27][28] The scattering intensity in the correlation length model is given by where I(q) is the scattering intensity, q is the scattering vector, and B is the incoherent background (see ESI † for further details and results of the fits). The first term, A/q n captures the low-q scattering and describes Porod scattering from the network, while the second term, C 1 þ qx ð Þ m captures the high-q scattering and is the Lorentzian function which characterizes local network structure. The correlation length x represents a weighted average of the polymer blob size in the network, and given that the gel concentration was kept constant, it is expected to be similar for all three gels. 25 The Porod and Lorentzian scale (A and C) and the Porod and Lorentzian exponents (n and m) along with x were obtained by a nonlinear least-squares fit of the data. The m exponent of all the gels are approximately 2 (Table S1, ESI †), which indicates that the polymers are behaving as though in a good solvent consistent with the assumption of good solvents conditions used to calculate the overlap concentration.
The main difference between the three polymers is captured by the first term, A/q n which defines the clustering strength from a large network structure and has been used as a method to evaluate the clustering strength of random polymer networks. 25,27,35,36 Note that while no quantitative relation can be inferred from this factor, a high clustering strength is associated with networks while low clustering strength corresponds to dissolved chains. 25,27,35,36 The clustering strength is significantly higher for PDHMc8 compared to PDHM5 and PDHM10 (Fig. 5). The value of q was chosen to be 0.04 nm À1   because it is low enough to be well within the Porod scattering regime. 25,27,35,36 As shown in Fig. 5, the clustering strength increases with decreasing average spacing between the stickers, l, for the random copolymers. For the clustered stickers gel, the high clustering strength is consistent with the scattering response being dominated by a network of average sticker spacing corresponding to l end . Thus, this result provides further evidence that the end blocks in the clustered polymer network have formed aggregates of the histidine-Ni complexes. Within these aggregates, the network structure shows similarity to the structure formed by random copolymers with average sticker spacing of l end . From Fig. 4 (and the fit parameters in Table S1, ESI †) the second term from fitting to the correlation length model is very similar for all three gels, which confirms that the gels are a disordered one-phase system. 25 As such, the SANS data indicates that the sticker clustering in the PDHMc8 polymers leads to the formation of aggregates of the histidine-Ni complexes, without inducing phase separation in the gels. As indicated by the lower plateau modulus (Fig. 3), these aggregates are connected by the midblock, which acts as the elastically active chains under linear deformation.

Effect of sticker clustering on mode of stress relaxation
Along with changes in the network topology, sticker clustering alters the mode of stress relaxation. In the PDHMc8 gels, network stress relaxation requires cooperative dissociation of multiple bonds, while the random copolymer gel networks relax stress when single bonds dissociate and bind to a new partner. This difference in stress relaxation mechanism is reflected by the increased temperature dependence of the relaxation time for the PDHMc8 gel compared to the random copolymers (Fig. 6). As shown in Fig. 6, the relaxation times of the gels increase in the order of t 5 o t c8 o t 10 . Additionally, the relaxation times of the gels are all longer than t d , which is the histidine-Ni 2+ complex dissociation time as measured in dilute solution reported by Tang et al. 20 These observations can be explained by considering the molecular mechanism for stress relaxation in the random and clustered copolymers as summarized in Fig. 7.
For the random copolymers, the stress relaxation time scale order of t d o t 5 o t 10 is consistent with the concept of bond renormalization put forth in the sticky Rouse model for linear polymers with stickers distributed evenly along the chain. 32 In the sticky Rouse model, the stress relaxation times measured in frequency sweeps correspond to the bond exchange times. For networks where the equilibrium constant, K eq c 1 (K eq = k a /k d where k a and k d are the rate constants for association and dissociation respectively), as is the case for the gels studied in this work most of the stickers are in the associated state. 20 As such, once a sticker dissociates there are very few exchange partners that are available. Thus, the newly dissociated stickers would have to explore the surrounding volume to find a new partner that is in the dissociated state. In the presence of   neighboring stickers, the volume a sticker can explore is reduced such that the stickers have a lower probability of finding a new partner that is in a dissociated state 32 (compare panel B and C in Fig. 7). As a result, the dissociated sticker will have to return to the same partner multiple times before successfully exchanging partners. The sticky Rouse model additionally predicts that the need for multiple bond dissociations results in the apparent activation energy for chain stress relaxation to be approximately 1.3 times higher than for single bond dissociation, but that activation energy should be independent of the number of stickers along the chain (S). 32 These predictions are consistent with our data as the activation energies for the random copolymers obtained from the Arrhenius fits (Fig. 6) are relatively similar (E a,5 = 69 AE 1 kJ mol À1 and E a,10 = 78 AE 6 kJ mol À1 ), but significantly higher than E a,d = 56 AE 4 kJ mol À1 (reported in ref. 20).
For the clustered copolymer, the increased relaxation time compared to the bond dissociation time (t d o t c8 ) has a different origin, as indicated by the higher activation energy of E a,c8 = 84 AE 1 kJ mol À1 for the PDHMc8 gel. The higher activation energy can be explained by the model proposed by Sing et al. 15 for telechelic polymers with multipart stickers at the chain ends. In this model, stress relaxation requires the cooperative dissociation of the stickers for pull-out of the chain ends. The model predicts that as the number of stickers at each chain end increases, the relaxation time, t, and activation energy, E a , will both increase because multiple stickers must be released prior to stress relaxation. The number of stickers that must be cooperatively dissociated for stress relaxation to occur can be estimated as so that x B 1.5. Since the average number of stickers per chain for the PDHMc8 polymer was 7.57 (Table 1), the average number of stickers per chain end is 3.8 under the assumption that the probability of the histidine-functionalized monomer being added to either end of the chain is equal. Thus, the estimate of x B 1.5 is lower than the calculated average number of stickers per chain end of 3.8 which indicates that the remaining stickers are either in a dissociated state or in intrachain bonds which do not hinder chain pull-out from occurring. The formation of just one intrachain bond within a chain end would reduce the number of interchain bonds from 3.8 to 1.8, which is close to the estimated value of x B 1.5. As demonstrated by several authors, the close proximity of the stickers within the chain ends can increase the propensity for formation of intrachain bonds. 31,37 Effect of sticker clustering on self-diffusion Self-diffusion measurements of the gels showed that diffusion was faster for the PDHMc8 gels compared to PDHM5 and PDHM10 gels despite the fact that PDHMc8 has more stickers than PDHM5 (Fig. 8). As previously reported for the PDHM5 and PDHM10 gels, 3 phenomenological superdiffusive scaling was observed for the PDHMc8 gel at smaller length scales prior to transitioning to Fickian scaling at d 2 B 200 mm 2 . Therefore, while sticker clustering has increased the diffusivity of the polymers compared to the random copolymers, it did not affect the length scales over which the apparent superdiffusive scaling occurs. New self-diffusion measurements for PDHM5 and PDHMc8 at lower temperatures of 25, 20 and 15 1C, showed that the same trends were observed across all the temperatures. Note that self-diffusion measurements for PDHM10 at temperatures below 35 1C were not experimentally accessible. Superdiffusive scaling that transitions to Fickian diffusion at length scales larger than the R g has been previously reported in other unentangled associative networks as well. 3,5,6,38 A previously proposed two-state model 5 demonstrated that the presence of two diffusive modes, with distinct diffusivities can lead to the appearance of superdiffusive scaling at length All gels prepared at a concentration of 25% (w/v). The dashed lines are fits to the two-state model 5 . Error bars represent one standard deviation of measurements performed in triplicate. Note that the data shown for PDHM5 and PDHM10 at 35 1C (filled symbols) was reported in an earlier publication, 3 while remaining data for PDHM5 was newly measured on the same polymers that were previously synthesized in ref. 3. scales larger than R g . This has been confirmed through simulations performed by Ramirez et al. 19 The two-state model 5 was able to capture self-diffusion data for such studies of unentangled associative networks showing superdiffusive scaling and is likewise able to capture the self-diffusion data for gels in this work (dashed lines in Fig. 8). As discussed in detail in an earlier publication, 3 the two modes of diffusion in the random copolymers PDHM5 and PDHM10 were proposed to be walking and hopping, based on the molecular model of Ramirez et al. 19 However, the molecular model does not consider the effect of sticker clustering. Thus, the analysis in this work will focus on fits to the two-state model which quantitively fits the data, but without assigning a molecular mechanism to the two diffusion modes. 5 The two-state model 5 hypothesizes that the polymers in an associative network exist in two states, the associated and mobile states, with distinct diffusivities, D A and D M (units: The polymers can interconvert between the two states with interconversion rates, k on and k off (units: s À1 ), with pseudo-first order kinetics. Since the physical details of the two diffusive states are not specified in the model, the model can be applied for the PDHMc8 gels without modification; however, k on and k off should not be taken as physical rate constants. While the individual model parameters D A , D M , k on and k off cannot be independently determined, as discussed by Tang et al., 5 the parameters of interest are the effective diffusivity in the large length-scale Fickian regime, given by D M,eff = D M /(1 + K eq ) and the anomaly index, gK eq = D A /D M Ák on / k off . Note that gK eq is inversely proportional to the extent of the superdiffusive regime and can take any value between 0 and 1.
The observation of faster diffusion in the PDHMc8 gels compared to the PDHM5 and PDHM10 gels across all temperatures was not expected based on the trend observed with the stress relaxation times (t 5 o t c8 o t 10 in Fig. 6). The faster diffusion is seen in Fig. 9(A) as a higher effective diffusivity in the large-length-scale Fickian regime, D M,eff . Following the approach of de Gennes, the diffusivity is related to the relaxation time through the relation D E R g 2 /t, 39 such that D B t À1 . Based on this relation, the effective diffusivity in the Fickian regime, D M,eff is expected to show the inverse of the trend with the relaxation times, such that D M,eff,5 4 D M,eff,c8 4 D M,eff,10 . Thus, while the trends observed for the random copolymers are consistent with the predictions of the sticky Rouse model, as discussed in ref. 3, the diffusing species measured for the PDHMc8 gels are not governed by the same time scales for mechanical relaxation as measured in the frequency sweeps. This discordant result can be further understood by considering the temperature dependence of the D M,eff along with the temperature dependence of the relaxation times. From the Arrhenius fits in Fig. 9(A), the activation energies for diffusion are E D,5 = 100 AE 10 kJ mol À1 and E D,c8 = 44 AE 20 kJ mol À1 . For PDHM5, the higher activation energy for diffusion compared to the activation energy for stress relaxation (E a,5 = 69 AE 1 kJ mol À1 ) indicates that more interchain bonds must be dissociated for the chain to diffuse several times its radius of gyration, R g . For the PDHMc8, not only is the activation energy for diffusion lower than the activation energy for stress relaxation (E a,c8 = 84 AE 1 kJ mol À1 ), the average value is lower than the activation energy for bond dissociation that was measured in dilute solution, E a,d = 56 AE 4 kJ mol À1 (reported in ref. 20). This suggests that self-diffusion in the PDHMc8 gels is mostly governed by single bond dissociations, which contrasts with the need for cooperative dissociation of multiple bonds for stress relaxation. The need for cooperative dissociation indicates that the elastically active chains are bound to the network through multiple interchain bonds. Thus, dissociation of a single bond in the elastically active chains will mean that several other interchain bonds are still in the associated state, and the chains will be unable to diffuse over length scales spanning several times its R g . This result implies that selfdiffusion measurement for the clustered polymer is dominated by defects such as chain loops or even ''superloops'' where just a few bond dissociations can result in a cluster of multiple chains diffusing a significant distance (Fig. 10(A)). 40 Similar results have been reported for diffusion measurements using fluorescence recovery after photobleaching (FRAP) on telechelic hydrophobically modified ethoxylated urethane (HEUR), where the defects were found to dominate self-diffusion measurements. 41 It should be noted that due to the statistical nature of the random copolymerization used to prepare the polymers in this work, the number of stickers per chain end will show a distribution. The distribution of the number of stickers per chain end can be approximated by a Poisson distribution 42 and as shown in Fig. 10(B), a small fraction of the chains (estimated to be 0.018) will exist as dangling chains (Fig. 10(A)). Defects such as chain loops and dangling chains can contribute to self-diffusion measurements but are elastically inactive since they cannot bridge two aggregates. This demonstrates that cooperative effects as indicated in the stress relaxation measurements can be seen even in the presence of a significant fraction of loop defects and dangling chains in the network. While the loop defects are elastically inactive, their presence in the telechelic hydrophobically modified ethoxylated urethane (HEUR) networks have been associated with the observation of shear thickening under nonlinear deformation. 43 This indicates that sticker clustering likely affects the nonlinear deformation behavior of these networks as well.
Sticker clustering and temperature have minimal effects on the extent of the superdiffusive scaling as seen in the very similar values of gK eq in Fig. 9(B). Since gK eq can be recast as gK eq = D A /D M,eff , gK eq can be interpreted as the ratio of apparent mobilities of molecules in the associative and mobile states. 5 Thus, the slightly larger values of gK eq for PDHMc8 suggests that changes in the mobility of molecules upon association are more pronounced for the clustered polymers. This is consistent with the results presented in this work, where the presence of the stickers in close proximity at the chain ends appears to drive the formation of multiple bonds for each chain end. While gK eq is larger for PDHMc8 compared to the random copolymers across all the temperatures investigated, this difference is small especially in comparison to the other unentangled associative networks previously investigated, which showed gK eq in the range of 0.06-0.001. 5,6 These other studies were performed on very different model systems, including hydrogels formed by linear proteins with four associating coiled-coil domains 5 and four-arm star-shaped polymers endfunctionalized with terpyridine moieties that are complexed with Zn 2+ in DMF. 6 Thus, these results suggest that gK eq is strongly influenced by the features of the gels that were kept constant between the random and clustered copolymers, including the binding chemistry and molecular weight of the polymers, compared to the relatively weak effect of sticker clustering and temperature.

Effect of temperature on the network topology
The opposite trends observed with the temperature dependence of the stress relaxation times and self-diffusion measurements indicates that temperature does not alter the network structure significantly. Decreasing the temperature should drive the system to favor bond association in both random and clustered sticker configurations. 15,31,32 In the clustered polymer, the close proximity of the stickers to one another enhances this effect, resulting in the stronger temperature dependence in the stress relaxation measurements. 17 However, the weaker dependence shown by the self-diffusion measurements indicates that the formation of more associations occurs within the loops and superloops themselves, such that the distribution of the interchain bonds is not altered between the elastically active and inactive parts of the network. This result was also reported by Feldman et al. who compared random and clustered copolymers using hydrogen-bonding polymer melts. 10 Their study reported that small angle X-ray scattering (SAXS) experiments performed over varying temperatures showed minimal change, indicating that temperature does not alter the network topology. Comparison to other studies Beyond the effect of sticker clustering, the bond chemistry is also an important factor in the final properties of the network. In addition to the study reported by Feldman et al. which compared random and clustered copolymers of hydrogenbonding polymer melts, 10 a second study was reported by Wu et al. which compared random and clustered copolymers of ionomers. 11 Both studies found that the stress relaxation times were increased with clustering; however, this effect was much stronger in the ionomers where the terminal relaxation times were not experimentally accessible. This is likely due to the propensity for the ionic groups to form large aggregates with high junction functionality 11 compared to the binary associations of the histidine-Ni complexes studies in this work. While the results reported by Wu et al. showed that the plateau in the storage modulus is lowered by clustering, as seen in this work, Feldman et al. found that clustering had no effect on the plateau in the storage modulus. The origin of this discrepancy is unknown, but it highlights the effect of the binding chemistry. It should be noted that the two studies focused on linear rheology experiments combined with SAXS and did not include self-diffusion measurements.

Conclusions
In this work, the effect of sticker clustering was investigated by comparing the properties of the model associative network to random copolymers with the same chemical composition. Sticker clustering was found to alter the network topology and stress relaxation mechanism, as indicated by frequency sweeps, small angle neutron scattering and self-diffusion studies using forced Rayleigh scattering. The network prepared from the clustered polymers consists of aggregates of the histidine-Ni complexes that are bridged by the nonassociative midblocks. The presence of multiple stickers at the chain ends results in an increase in the stress relaxation times due to the need for cooperative dissociation of multiple bonds for chain pull-out to occur. The weak effect of temperature on the self-diffusion measurements for the clustered polymer further revealed that the diffusion measurements were dominated by defects, such as superloops, that have been reported for other telechelic polymers (with a single associative group at the chain ends). This weak temperature dependence was not observed for the random copolymer, indicating that the clustering of the stickers drives the formation of these loop defects which are known to affect nonlinear rheology properties of associative networks. Additionally, the observation of phenomenological superdiffusive scaling here shows strong similarities to those observed in the random copolymer, which provide further insights for the development of the molecular model for diffusion.

Conflicts of interest
There are no conflicts to declare.