This study reports a general scenario for the out-of-equilibrium features of collapsing polymeric architectures. We use molecular dynamics simulations to characterize the coarsening kinetics, in bad solvent, for several macromolecular systems with an increasing degree of structural complexity. In particular, we focus on: flexible and semiflexible polymer chains, star polymers with 3 and 12 arms, and microgels with both ordered and disordered networks. Starting from a powerful analogy with critical phenomena, we construct a density field representation that removes fast fluctuations and provides a consistent characterization of the domain growth. Our results indicate that the coarsening kinetics presents a scaling behaviour that is independent of the solvent quality parameter, in analogy to the time–temperature superposition principle. Interestingly, the domain growth in time follows a power-law behaviour that is approximately independent of the architecture for all the flexible systems; while it is steeper for the semiflexible chains. Nevertheless, the fractal nature of the dense regions emerging during the collapse exhibits the same scaling behaviour for all the macromolecules. This suggests that the faster growing length scale in the semiflexible chains originates just from a faster mass diffusion along the chain contour, induced by the local stiffness. The decay of the dynamic correlations displays scaling behavior with the growing length scale of the system, which is a characteristic signature in coarsening phenomena.

We perform simulations to compute the effective potential between the centers-of-mass of two polymers with reversible bonds. We investigate the influence of the topology on the potential by employing linear and ring backbones for the precursor (unbonded) polymer, finding that it leads to qualitatively different effective potentials. In the linear and ring cases the potentials can be described by Gaussians and generalized exponentials, respectively. The interactions are more repulsive for the ring topology, in analogy with known results in the absence of bonding. We also investigate the effect of the specific sequence of the reactive groups along the backbone (periodic or with different degrees of randomness), establishing that it has a significant impact on the effective potentials. When the reactive sites of both polymers are chemically orthogonal so that only intramolecular bonds are possible, the interactions become more repulsive the closer to periodic the sequence is. The opposite effect is found if both polymers have the same types of reactive sites and intermolecular bonds can be formed. We test the validity of the effective potentials in solution, in a broad range of concentrations from high dilution to far above the overlap concentration. For this purpose, we compare simulations of the effective fluid and test particle route calculations with simulations of the real all-monomer system. Very good agreement is found for the reversible linear polymers, indicating that unlike in their nonbonding counterparts many-body effects are minor even far above the overlap concentration. The agreement for the reversible rings is less satisfactory, and at high concentration the real system does not show the clustering behavior predicted by the effective potential. Results similar to the former ones are found for the partial self-correlations in ring/linear mixtures. Finally, we investigate the possibility of creating, at high concentrations, a gel of two interpenetrated reversible networks. For this purpose we simulate a 50/50 two-component mixture of reversible polymers with orthogonal chemistry for the reactive sites, so that intermolecular bonds are only formed between polymers of the same component. As predicted by both the theoretical phase diagram and the simulations of the effective fluid, the two networks in the all-monomer mixture do not interpenetrate, and phase separation (demixing) is observed instead.

No abstract

This study reports a general scenario for the out-of-equilibrium features of collapsing polymeric architectures. We use molecular dynamics simulations to characterize the coarsening kinetics, in bad solvent, for several macromolecular systems with an increasing degree of structural complexity. In particular, we focus on: flexible and semiflexible polymer chains, star polymers with 3 and 12 arms, and microgels with both ordered and disordered networks.Starting from a powerful analogy with critical phenomena, we construct a density field representation that removes fast fluctuations and provides a consistent characterization of the domain growth. Our results indicate that the coarsening kinetics presents a scaling behaviour that is independent of the solvent quality parameter, in analogy to time-temperature superposition principle. Interestingly, the domain growth in time follows a power-law behaviour that is approximately independent of the architecture for all the flexible systems; while it is steeper for the semiflexible chains. Nevertheless, the fractal nature of the dense regions emerging during the collapse exhibits the same scaling behaviour for all the macromolecules.This suggests that the faster growing length scale in the semiflexible chains originates just from a faster mass diffusion along the chain contour, induced by the local stiffness. The decay of the dynamic correlations displays scaling behavior with the growing length scale of the system, which is a characteristic signature in coarsening phenomena. I. INTRODUCTIONUnderstanding the collapse of, fully polymeric, topologically complex objects in a bad solvent is of broad importance, because of its relevance in the early stages of the protein folding [1,2], and its connections with arresting processes or aging phenomena [3,4]. The thermodynamic and structural properties of the macromolecular collapse in solution, occurring when the quality of the solvent is decreased below a critical value, have been exhaustively investigated over the years [5][6][7][8][9][10], and at equilibrium, the dynamic and static behavior are well known above and below the volume phase transition [11]. Comparatively, a general framework for the non-equilibrium aspects, as the kinetics of the collapse, is still lacking. The development of new experimental setups for small angle X-ray scattering or single-molecule fluorescence spectroscopy allows to monitor the collapse arXiv:1912.03551v1 [cond-mat.soft]