Molecular dynamics and neutron scattering study of the dependence of polyelectrolyte dendrimer conformation on counterion behavior

Atomistic molecular dynamics (MD) simulations and contrast variation small angle neutron scattering (SANS) have been combined to investigate the Generation-5 polyelectrolyte polyamidoamine starburst dendrimer. This work reveals the dendrimer conformational dependence on counterion association at different levels of molecular charge. The accuracy of the simulations is verified through satisfactory comparison between modeled results, such as excess intra-dendrimer scattering length density distribution and hydration level, and their experimental counterparts. While the counterion distributions are not directly measureable with SANS, the spatial distribution of the counterions and their dendrimer association are extracted from the validated MD equilibrium trajectories. It is found that the conformation of the charged dendrimer is strongly dependent on the counterion association. Sensitivity of the distribution of counterions around charged amines to the counterion valency is qualitatively explained by adopting Langmuir adsorption theory. Moreover, via extending the concept of electrical double layer for compact charged colloids, we define an effective radius of a charged dendrimer including the spatial distribution of counterions in its vicinity. Within the same framework, the correlation between the strength of intra-dendrimer electrostatic repulsion and the counterion valency and dynamics is also addressed.     

Anomalous mixing behavior of polyisobutylene/polypropylene blends: molecular dynamics simulation study

The unusual mixing behavior of polyisobutylene (PIB) with head-to-head (hhPP) and head-to-tail polypropylene (PP) is studied using large-scale molecular dynamics (MD). The heats of mixing and Flory chi parameters were computed from MD simulations of both blends using a united atom model. The chi parameters from the simulations were estimated from the structure factors using the random phase approximation in analogy with neutron scattering (SANS) experiments. MD simulations for syndiotactic hhPP/PIB predicted a lower critical solution temperature with a chi parameter in very good agreement with SANS experiments on the atactic hhPP/PIB blend. MD simulations also predicted that the isotactic PP/PIB blend was immiscible at high molecular weight in qualitative agreement with cloud point measurements on atactic PP/PIB.     

Conformational effect on small angle neutron scattering behavior of interacting polyelectrolyte solutions: a perspective of integral equation theory

We present small angle neutron scattering (SANS) measurements of deuterium oxide (D(2)O) solutions of linear and star sodium poly(styrene sulfonate) (NaPSS) as a function of polyelectrolyte concentration. Emphasis is on understanding the dependence of their SANS coherent scattering cross section I(Q) on the molecular architecture of single polyelectrolyte. The key finding is that for a given concentration, star polyelectrolytes exhibit more pronounced characteristic peaks in I(Q), and the position of the first peak occurs at a smaller Q compared to their linear counterparts. Based on a model of integral equation theory, we first compare the SANS experimental I(Q) of salt-free polyelectrolyte solutions with that predicted theoretically. Having seen their satisfactory qualitative agreement, the dependence of counterion association behavior on polyelectrolyte geometry and concentration is further explored. Our predictions reveal that the ionic environment of polyelectrolyte exhibits a strong dependence on polyelectrolyte geometry at lower polyelectrolyte concentration. However, when both linear and star polyelectrolytes exceed their overlap concentrations, the spatial distribution of counterion is found to be essentially insensitive to polyelectrolyte geometry due to the steric effect.     

Structure of aqueous proline via parallel tempering molecular dynamics and neutron diffraction

The structure of aqueous L-proline amino acid has been the subject of much debate centering on the validity of various proposed models, differing widely in the extent to which local and long-range correlations are present. Here, aqueous proline is investigated by atomistic, replica exchange molecular dynamics simulations, and the results are compared to neutron diffraction and small angle neutron scattering (SANS) data, which have been reported recently (McLain, S.; Soper, A.; Terry, A.; Watts, A. J. Phys. Chem. B 2007, 111, 4568). Comparisons between neutron experiments and simulation are made via the static structure factor S(Q) which is measured and computed from several systems with different H/D isotopic compositions at a concentration of 1:20 molar ratio. Several different empirical water models (TIP3P, TIP4P, and SPC/E) in conjunction with the CHARMM22 force field are investigated. Agreement between experiment and simulation is reasonably good across the entire Q range although there are significant model-dependent variations in some cases. In general, agreement is improved slightly upon application of approximate quantum corrections obtained from gas-phase path integral simulations. Dimers and short oligomeric chains formed by hydrogen bonds (frequently bifurcated) coexist with apolar (hydrophobic) contacts. These emerge as the dominant local motifs in the mixture. Evidence for long-range association is more equivocal: No long-range structures form spontaneously in the MD simulations, and no obvious low-Q signature is seen in the SANS data. Moreover, associations introduced artificially to replicate a long-standing proposed mesoscale structure for proline correlations as an initial condition are annealed out by parallel tempering MD simulations. However, some small residual aggregates do remain, implying a greater degree of long-range order than is apparent in the SANS data.     

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X‐ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra‐molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. © 2012 Wiley Periodicals, Inc.     

Structural investigation of PAMAM dendrimers in aqueous solutions using small-angle neutron scattering: effect of generation

We investigate a series of poly(amidoamine) starburst dendrimers (PAMAM) of different generations in acidic, aqueous solutions using small-angle neutron scattering (SANS). While the overall molecular size is found to be practically unaffected by a pD change, a strong generational dependence of counterion association is revealed. Upon increasing the dendrimer generation, the effective charge obtained from our SANS experiments only shows a small increase in contrast to the nearly exponential increase predicted by a recent atomic simulation. We also find that with the same degree of molecular protonation the specific counterion association, which is defined as the ratio of bound chloride anions to positively charged amines in solutions, is larger for higher-generation PAMAM dendrimer. The associated counterion density also increases upon increasing generation number.     

Structure and rheology of hyperbranched and dendritic polymers. I. Modification and characterization of poly(propyleneimine) dendrimers with acetyl groups

Commercially available fourth and fifth generation poly(propyleneimine) (PPI) dendrimers were functionalized with acetyl chloride and deuterated acetyl chloride. Their solution properties in water and D₂O were measured with dilution viscometry, densitometry, rheology, and small‐angle neutron scattering (SANS) and compared to molecular modeling. Both the acetylated and PPI dendrimers exhibited Newtonian rheology in solution at all concentrations, but the functionalized dendrimers were less viscous than the nonacetylated dendrimers at an equal weight fraction (50 wt %). The acetylated dendrimers exhibited a pronounced structure peak in SANS, however, that was not evident for PPI in solution and a greatly enhanced solubility. This structure peak, evident at concentrations as low as 0.2 wt %, was evidence for long‐range electrostatic interdendrimer forces, which were screened by added salt. A quantitative agreement was obtained between the dilute‐limiting absolute scattering spectra of both the nonacetylated and acetylated dendrimers in solution with model calculations via a homogeneous spherical model and input parameters independently obtained from dilution viscometry or direct calculation. The combined measurements verified significant solvent penetration for both dendrimer types. The form factors measured in this manner were also in good quantitative agreement with the results of molecular dynamics simulations, which pointed to significant backfolding of the terminal groups. SANS and rheology measurements at higher concentrations suggested dendrimer clustering and interpenetration with increasing concentration, leading to less structure and lower viscosity than would be predicted from the dilute‐limiting behavior. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 857–873, 2000     

Similarities between polyelectrolyte gels and biopolymer solutions

Small angle neutron scattering (SANS) measurements and osmotic swelling pressure measurements are reported for polyelectrolyte gels and solutions under nearly physiological conditions. A synthetic polymer (sodium-polyacrylate) and three biopolymers (DNA, hyaluronic acid, and polyaspartic acid) are studied. The neutron scattering response of these anionic polyelectrolytes is closely similar, indicating that at larger length scales the organization of the polymer molecules is not significantly affected by the fine details of the molecular architecture (e.g., size and chemical structure of the monomer unit, type of polymer backbone). The results suggest that specific interactions between the polyelectrolyte chains and the surrounding monovalent cations are negligible. It is found that the osmotic compression modulus of these biopolymer solutions determined from the analysis of the SANS response decreases with increasing chain persistence length. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 3679–3686, 2006     

N,N-dimethylformamide-induced phase separation of hexafluoroisopropanol-water mixtures

We found that addition of N,N-dimethylformamide (DMF) induces phase separation of 1,1,1,3,3,3-hexafluoro-2-propanol (HFIP)-water mixtures. The phase diagram of a DMF-HFIP-water ternary system at 298 K showed that phase separation occurs in a closed-loop area in the water-rich mole fraction range of x(H(2)O) > ～0.4. To clarify the mechanism of DMF-induced phase separation of DMF-HFIP-water mixtures at the molecular level, small-angle neutron scattering (SANS) and (1)H and (13)C NMR measurements were conducted on the mixtures with varying DMF concentrations along a volume ratio of HFIP to water of 1?:?1 (x(S)(HFIP) = 0.147). Additionally, the solvation structure of DMF in water and HFIP-water mixtures was elucidated by molecular dynamics (MD) simulations. The SANS results revealed that the inherent heterogeneity of HFIP-water mixtures is increased with increasing DMF concentration toward the lower phase separation concentration, but decreased when the DMF concentration further increased beyond the upper phase separation one. (1)H and (13)C NMR measurements and MD simulations suggested that preferential solvation of the hydrophobic moiety of DMF by HFIP is the main driver of the phase behaviour of the DMF-HFIP-water system.     

Experimental evidences for molecular origin of low-Q peak in neutron/x-ray scattering of 1-alkyl-3-methylimidazolium bis(trifluoromethanesulfonyl)amide ionic liquids

Short- and long-range liquid structures of [C(n)mIm(+)][TFSA(-)] with n = 2, 4, 6, 8, 10, and 12 have been studied by high-energy x-ray diffraction (HEXRD) and small-angle neutron scattering (SANS) experiments with the aid of MD simulations. Observed x-ray structure factor, S(Q), for the ionic liquids with the alkyl-chain length n > 6 exhibited a characteristic peak in the low-Q range of 0.2-0.4 A?(-1), indicating the heterogeneity of their ionic liquids. SANS profiles I(H)(Q) and I(D)(Q) for the normal and the alkyl group deuterated ionic liquids, respectively, showed significant peaks for n = 10 and 12 without no form factor component for large spherical or spheroidal aggregates like micelles in solution. The peaks for n = 10 and 12 evidently disappeared in the difference SANS profiles ΔI(Q) [=I(D)(Q) - I(H)(Q)], although that for n = 12 slightly remained. This suggests that the long-range correlations originated from the alkyl groups hardly contribute to the low-Q peak intensity in SANS. To reveal molecular origin of the low-Q peak, we introduce here a new function; x-ray structure factor intensity at a given Q as a function of r, S(Q) (peak)(r). The S(Q) (peak)(r) function suggests that the observed low-Q peak intensity depending on n is originated from liquid structures at two r-region of 5-8 and 8-15 A? for all ionic liquids examined except for n = 12. Atomistic MD simulations are consistent with the HEXRD and SANS experiments, and then we discussed the relationship between both variations of low-Q peak and real-space structure with lengthening the alkyl group of the C(n)mIm.     

The structures of complexes between polyethylene imine and sodium dodecyl sulfate in D2O: a scattering study

The association between a highly branched polyelectrolyte with ionizable groups, polyethylene imine (PEI), and an anionic surfactant, sodium dodecyl sulfate (SDS), has been investigated at two pH values, using small-angle neutron and light scattering. The scattering data allow us to obtain a detailed picture of the association structures formed. Small-angle neutron scattering (SANS) measurements in solutions containing highly charged PEI at low pH and low SDS concentrations indicate the presence of disklike aggregates. The aggregates change to a more complex three-dimensional structure with increasing surfactant concentration. One pronounced feature in the scattering curves is the presence of a Bragg-like peak at high q-values observed at a surfactant concentration of 4.2 mM and above. This scattering feature is attributed to the formation of a common well-ordered PEI/SDS structure, in analogue to what has been reported for other polyelectrolyte-surfactant systems. Precipitation occurred at the charge neutralization point, and X-ray diffraction measurements on the precipitate confirmed the existence of an ordered structure within the PEI/SDS aggregates, which was identified as a lamellar internal organization. Polyethylene imine has a low charge density in alkaline solutions. At pH 10.1 and under conditions where the surfactant was contrast matched, the SANS scattering curves showed only small changes with increasing surfactant concentration. This suggests that the polymer acts as a template onto which the surfactant molecules aggregate. Data from both static light scattering and SANS recorded under conditions where SDS and to a lower degree PEI contribute to the scattering were found to be consistent with a structure of stacked elliptic bilayers. These structures increased in size and became more compact as the surfactant concentration was increased up to the charge neutralization point.     

Solvent quality changes the structure of G8 PAMAM dendrimer, a disagreement with some experimental interpretations

We have performed approximately 20-40 ns of molecular dynamics (MD) simulations for the generation 8 PAMAM dendrimer in explicit water under varying pH conditions to study the structure of the dendrimer (approximately 156,738 atoms at low pH). This is the first report of such a long MD simulation of a larger generation PAMAM dendrimer including the effect of salt and counterions with explicit water molecules. We find that changing the pH from a high value (approximately 12) to a low value (approximately 3) changes the radius of gyration from Rg = 37.8 to 43.1 A (increasing by 13%). We also find significant back-folding of the primary amines and a large amount of water penetration inside the polymer. The increase in size with decrease in pH is consistent with our earlier studies on G3-G6 and agrees with the Monte Carlo theory by Welch and Muthukumar of G8 (Macromolecules, 1998, 31, 5892) and the experiments on G5 and G8 PAMAM dendrimer by Topp et al. (Macromolecules, 1999, 32, 7232). However, these results disagree dramatically with the interpretations of SANS experiments of G8 PAMAM dendrimers by Nisato et al. (Macromolecules, 2000, 33, 4172) who observe no change in the size of the dendrimer with variations of solution pH and ionic strength. We assume that the disagreement might arise from neglecting nonspherical shape, penetration of water and ions into the core, and aggregation, all of which might depend on pH.     

Structure and properties of the dendron-encapsulated polyoxometalate (C(52)H(60)NO(12))(12)[(Mn(H(2)O))(3)(SbW(9)O(33))(2)], a first generation dendrizyme

Combining analytical and theoretical methods, we present a detailed study of a heteropolytungstate cluster encapsulated in a shell of dendritically branching surfactants, namely (C(52)H(60)NO(12))(12)[(Mn(H(2)O))(3)(SbW(9)O(33))(2)], 3. This novel surfactant-encapsulated cluster (SEC) self-assembles spontaneously from polyoxometalate-containing solutions treated with a stoichiometric amount of dendrons. Compound 3 exhibits a discrete supramolecular architecture in which a single polyoxometalate anion resides in a compact shell of dendrons. Our approach attempts to combine the catalytic activity of polyoxometalates with the steric properties of tailored dendritic surfactants into size-selective catalytic systems. The structural characterization of the SEC is based on analytical ultracentrifugation (AUC) and small-angle neutron scattering (SANS). The packing arrangement of dendrons at the cluster surface is gleaned from molecular dynamics (MD) simulations, which suggests a highly porous shell structure due to the dynamic formation of internal clefts and cavities. From analysis of the MD trajectory of 3, a theoretical neutron-scattering function is derived that is in good agreement with experimental SANS data. Force field parameters used in MD simulations are partially derived from a quantum mechanical geometry optimization of [(Zn(H(2)O))(3)(SbW(9)O(33))(2)](12)(-), 2b, at the density functional theory (DFT) level. DFT calculations are corroborated by X-ray structure analysis of Na(6)K(6)[(Zn(H(2)O))(3)(SbW(9)O(33))(2)].23H(2)O, which is isostructural with the catalytically active Mn derivative 2a. The combined use of theoretical and analytical methods aims at rapidly prototyping smart catalysts ("dendrizymes"), which are structurally related to naturally occurring metalloproteins.     

A hydrogen-sensitive analysis of the structure of deuterated liquid formamide

Structural investigations of deuterated liquid formamide were performed by using neutron scattering, ab-initio calculations and classical Molecular Dynamics (MD) simulations. The recorded neutron data are analysed to yield the total structure factor S_M (Q), the molecular form factor F (Q), the distinct pair correlation function g_L(r) and particularly the deuterium-oxygen signature of H-bond interactions. Neutron scattering data, as well as recent x-ray studies, clearly show that the local order of the liquid is largely described by one dimer, two trimers and one tetramer. Molecular Dynamics simulations show that neutron scattering data can be reproduced by three different force fields.     

Self-aggregation of mixtures of oppositely charged polyelectrolytes and surfactants studied by rheology, dynamic light scattering and small-angle neutron scattering

In this study, the phase behavior, structure and properties of systems composed of the cationic, cellulose-based polycation JR 400 and the anionic surfactants sodium dodecylbenzenesulfonate (SDBS) or sodium dodecylethoxysulfate (SDES), mainly in the semidilute regime, were examined. This system shows the interesting feature of a very large viscosity increase by nearly 4 orders of magnitude as compared to the pure polymer solution already at very low concentrations of 1 wt%. By using rheology, dynamic light scattering (DLS), and small-angle neutron scattering (SANS), we are able to deduce systematic correlations between the molecular composition of the systems (characterized by the charge ratio Z=[+(polymer)]/[?(surfactant)]), their structural organization and the resulting macroscopic flow behavior. Mixtures in the semidilute regime with an excess of polycation charge form highly viscous network structures containing rodlike aggregates composed of surfactant and polyelectrolyte that are interconnected by the long JR 400 chains. Viscosity and storage modulus follow scaling laws as a function of surfactant concentration (η~c(s)(4); G(0)~c(s)(1.5)) and the very pronounced viscosity increase mainly arises from the strongly enhanced structural relaxation time of the systems. In contrast, mixtures with excess surfactant charges form solutions with viscosities even below those of the pure polymer solution. The combination of SANS, DLS, and rheology shows that the structural, dynamical, and rheological properties of these oppositely charged polyelectrolyte/surfactant systems can be controlled in a systematic fashion by appropriately choosing the systems composition.     

Side chain length affects backbone dynamics in poly(3‐alkylthiophene)s

Charge transport in conjugated polymers may be governed not only by the static microstructure but also fluctuations of backbone segments. Using molecular dynamics simulations, we predict the role of side chains in the backbone dynamics for regiorandom poly(3‐alkylthiophene‐2,5‐diyl)s (P3ATs). We show that the backbone of poly(3‐dodecylthiophene‐2‐5‐diyl) (P3DDT) moves faster than that of poly(3‐hexylthiophene‐2,5‐diyl) (P3HT) as a result of the faster motion of the longer side chains. To verify our predictions, we investigated the structures and dynamics of regiorandom P3ATs with neutron scattering and solid state NMR. Measurements of spin‐lattice relaxations (T₁) using NMR support our prediction of faster motion for side chain atoms that are farther away from the backbone. Using small‐angle neutron scattering (SANS), we confirmed that regiorandom P3ATs are amorphous at about 300 K, although microphase separation between the side chains and backbones is apparent. Furthermore, quasi‐elastic neutron scattering (QENS) reveals that thiophene backbone motion is enhanced as the side chain length increases from hexyl to dodecyl. The faster motion of longer side chains leads to faster backbone dynamics, which in turn may affect charge transport for conjugated polymers. © 2018 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2018 , 56, 1193–1202     

Light,Small Angle Neutron and X-Ray Scattering from Gels

This paper gives two examples of experiments that demonstrate the power of small angle scattering techniques in the study of swollen polymer networks. First, it is shown how the partly ergodic character of these systems is directly detected by neutron spin echo experiments. The observed total field correlation function of the intensity scattered from a neutral gel allows the ergodic contribution to be directly distinguished from the non ergodic part, at values of transfer wave vector q that lie well beyond the range of dynamic light scattering. The results can be compared with those obtained at much lower q from visible light scattering. Second, a recent application of small angle X-ray (SAXS) and neutron (SANS) scattering is described for a polyelectrolyte molecule, DNA, in semi-dilute solutions under near-physiological conditions. For these observations, the divalent ion normally present, calcium, is replaced by an equivalent ion, strontium. The comparison between SANS and SAXS yields a quantitative picture of the cloud of divalent counter-ions around the central DNA core. At physiological conditions, the cloud is thinner than that predicted on the basis of the Debye screening length but thicker than if the counter-ions were condensed on the DNA chain.     

Cationic polythiophene-surfactant self-assembly complexes: phase transitions, optical response, and sensing

The absorption and photoluminescence spectra of the cationic conjugated polyelectrolyte poly[3-(6-trimethylammoniumhexyl)thiophene] (P3TMAHT) were observed to be dramatically altered in the presence of anionic surfactants due to self-assembly through ionic complex formation. Small-angle neutron scattering (SANS), UV/vis, and photoluminescence spectroscopy were used to probe the relationship between the supramolecular complex organization and the photophysical response of P3TMAHT in the presence of industrially important anionic surfactants. Subtle differences in the surfactant mole fraction and chemical structure (e.g., chain length, headgroup charge density, perfluorination) result in marked variations in the range and type of complexes formed, which can be directly correlated to a unique colorimetric and fluorimetric fingerprint. Our results show that P3TMAHT has potential as an optical sensor for anionic surfactants capable of selectively identifying distinct structural subgroups through dual mode detection.     

A powerful computational crystallography method to study ice polymorphism

Classical molecular dynamics (MD) simulations are employed as a tool to investigate structural properties of ice crystals under several temperature and pressure conditions. All ice crystal phases are analyzed by means of a computational protocol based on a clustering approach following standard MD simulations. The MD simulations are performed by using a recently published classical interaction potential for oxygen and hydrogen in bulk water, derived from neutron scattering data, able to successfully describe complex phenomena such as proton hopping and bond formation/breaking. The present study demonstrates the ability of the interaction potential model to well describe most ice structures found in the phase diagram of water and to estimate the relative stability of 16 known phases through a cluster analysis of simulated powder diagrams of polymorphs obtained from MD simulations. The proposed computational protocol is suited for automated crystal structure identification.