Binding mechanisms in dendrimer-surfactant complexes.

Molecular dynamics simulations were employed to investigate the impact of interactions between dendritic polyeclectrolytes and amphiphilic surfactants on the supramolecular complex formation. We recognize two crucial parameters that govern association of surfactants within dendrimers: surfactant hydrophobicity, ε^{*}, and dendrimer generation, G. We find that depending on the values of ε^{*} and G encapsulation of surfactants by dendrimers is either noncooperative or cooperative. The noncooperative binding is characterized by absorption of surfactants as unimers, whereas in cooperative binding absorption of unimers is followed by aggregate formation through hydrophobic attractions between the surfactant tails. Our results provide guidelines for controlled encapsulation of guest molecules in dendrimer-based guest-host complexes.

Spatial segregation of mixed-sized counterions in dendritic polyelectrolytes.

Langevin dynamics simulations are utilized to study the structure of a dendritic polyelectrolyte embedded in two component mixtures comprised of conventional (small) and bulky counterions. We vary two parameters that trigger conformational properties of the dendrimer: the reduced Bjerrum length, [Formula: see text], which controls the strength of electrostatic interactions and the number fraction of the bulky counterions, [Formula: see text], which impacts on their steric repulsion. We find that the interplay between the electrostatic and the counterion excluded volume interactions affects the swelling behavior of the molecule. As compared to its neutral counterpart, for weak electrostatic couplings the charged dendrimer exists in swollen conformations whose size remains unaffected by [Formula: see text]. For intermediate couplings, the absorption of counterions into the pervaded volume of the dendrimer starts to influence its conformation. Here, the swelling factor exhibits a maximum which can be shifted by increasing [Formula: see text]. For strong electrostatic couplings the dendrimer deswells correspondingly to [Formula: see text]. In this regime a spatial separation of the counterions into core-shell microstructures is observed. The core of the dendrimer cage is preferentially occupied by the conventional ions, whereas its periphery contains the bulky counterions.

Charged Dendrimers with Finite-Size Counterions.

We report on the structure of dendritic polyelectrolytes accompanied by counterions in a good, salt-free, implicit solvent using Langevin dynamics simulations and a Flory-type approach. Our focus is on the modification of charged dendrimer conformations via the strength of electrostatic interactions and the counterion excluded volume. We study the effects caused by charges by varying the reduced Bjerrum length, λB*, between the extremes of weak and strong electrostatic interactions. The counterion excluded volume was controlled by the size of ions. We investigate counterions ranging from conventional ones, with the size comparable to the monomer size, to bulky ions. Our results indicate that, as compared to neutral dendrimers, dendritic polyelectrolytes exist in swollen conformations, and the degree of swelling changes non-monotonically with increasing λB*. For weak electrostatic couplings, counterion density within dendrimers is minor and their radius of gyration subtly exceeds the size of neutral dendrimers. For intermediate electrostatic couplings, Coulomb attraction between opposite charges promotes absorption of ions into dendrimers' pervaded volume and counterion condensation on charged monomers. As a result, counterion density inside dendrimers abruptly increases and the ionic size starts to play a crucial role. In this regime, we observe that swelling of dendrimers reaches its maximum and is more pronounced for bulky counterions. For strong electrostatic couplings, complete condensation of conventional counterions proceeds, whereas for bulky ions condensation remains partial. In this regime, dendrimers deswell. In particular, in the presence of conventional ions, dendrimers collapse into globules, while, for bulky counterions, deswelling is suppressed.

Detachment of semiflexible polymer chains from a substrate: a molecular dynamics investigation.

Using Molecular Dynamics simulations, we study the force-induced detachment of a coarse-grained model polymer chain from an adhesive substrate. One of the chain ends is thereby pulled at constant speed off the attractive substrate and the resulting saw-tooth profile of the measured mean force ⟨f⟩ vs height D of the end-segment over the plane is analyzed for a broad variety of parameters. It is shown that the observed characteristic oscillations in the ⟨f⟩-D profile depend on the bending and not on the torsional stiffness of the detached chains. Allowing for the presence of hydrodynamic interactions (HI) in a setup with explicit solvent and dissipative particle dynamics-thermostat, rather than the case of Langevin thermostat, one finds that HI have little effect on the ⟨f⟩-D profile. Also the change of substrate affinity with respect to the solvent from solvophilic to solvophobic is found to play negligible role in the desorption process. In contrast, a changing ratio ε(s)(B)/ε(s)(A) of the binding energies of A- and B-segments in the detachment of an AB-copolymer from adhesive surface strongly changes the ⟨f⟩-D profile whereby the B-spikes vanish when ε(s)(B)/ε(s)(A)<0.15. Eventually, performing an atomistic simulation of (bio)-polymers, we demonstrate that the simulation results, derived from our coarse-grained model, comply favorably with those from the all-atom simulation.

Thermal decomposition of a honeycomb-network sheet: a molecular dynamics simulation study.

The thermal degradation of a graphene-like two-dimensional honeycomb membrane with bonds undergoing temperature-induced scission is studied by means of Molecular Dynamics simulation using Langevin thermostat. We demonstrate that at lower temperature the probability distribution of breaking bonds is highly peaked at the rim of the membrane sheet whereas at higher temperature bonds break at random everywhere in the hexagonal flake. The mean breakage time τ is found to decrease with the total number of network nodes N by a power law τ ∝ N(-0.5) and reveals an Arrhenian dependence on temperature T. Scission times are themselves exponentially distributed. The fragmentation kinetics of the average number of clusters can be described by first-order chemical reactions between network nodes n(i) of different coordination. The distribution of fragments sizes evolves with time elapsed from initially a δ-function through a bimodal one into a single-peaked again at late times. Our simulation results are complemented by a set of 1st-order kinetic differential equations for n(i) which can be solved exactly and compared to data derived from the computer experiment, providing deeper insight into the thermolysis mechanism.

Force-induced breakdown of flexible polymerized membrane.

We consider the fracture of a free-standing two-dimensional (2D) elastic-brittle network to be used as protective coating subject to constant tensile stress applied on its rim. Using a molecular-dynamics simulation with a Langevin thermostat, we investigate the scission and recombination of bonds, and the formation of cracks in the 2D graphenelike hexagonal sheet for different pulling force f and temperature T. We find that bond rupture occurs almost always at the sheet periphery, and the first mean breakage time <τ> of bonds decays with membrane size as <τ> ∝N(-ß), where ß≈0.50±0.03 and N denotes the number of atoms in the membrane. The probability distribution of bond scission times t is given by a Poisson function W(t)∝t(1/3)exp(-t/<τ>). The mean failure time <τ(r)> necessary to rip off the sheet declines with growing size N as a power law <τ(r)>∝N(-φ(f)). We also find <τ(r)>∝exp(ΔU(0)/k(B)T), where the nucleation barrier for crack formation ΔU(0)∝f(-2), in agreement with Griffith's theory. <τ(r)> displays an Arrhenian dependence of <τ(r)> on temperature T. Our results indicate a rapid increase in crack spreading velocity with growing external tension f.

Thermal degradation of unstrained single polymer chain: non-linear effects at work.

We examine the thermally induced fracture of an unstrained polymer chain of discrete segments coupled by an anharmonic potential by means of molecular dynamics simulation with a Langevin thermostat. Cases of both under- and over-damped dynamics are investigated, and a comparison with recent studies of bond scission in model polymers with harmonic interactions is performed. We find that the polymer degradation changes qualitatively between the inertial regime and that of heavily damped dynamics. The role of bond healing (recombination) is also studied and probability distributions for the recombination times and overstretched bond lengths are obtained. Our extensive simulations reveal many properties of the scission dynamics in agreement with the notion of random breakdown of independent bonds, e.g., the mean time of chain rupture, <τ> follows an Arrhenian behavior with temperature T, and depends on the number of bonds N in the polymer as <τ> ∝ N(-1). In contrast, the rupture rates of the individual bonds along the polymer backbone indicate clearly the presence of self-induced inhomogeneity resulting from the interplay of thermal noise and nonlinearity. Eventually we examine the fragmentation kinetics during thermolysis. We demonstrate that both the probability distribution function of fragment sizes as well as the mean length of fragments at subsequent times t characterize degradation as predominantly a first order reaction.