Scaling and Interactions of Linear and Ring Polymer Brushes via DPD Simulations.

Single and double layers of polymer coated surfaces are investigated by means of Dissipative Particle Dynamics (DPD), focusing on the difference between grafted ring and linear chains. Several different surface coverages σ , as well as chain lengths N and brush separations D, are analyzed for athermal, i.e., good solvent, conditions. The size in the form of the radius of gyration R g , the shape as asphericity δ ∗ , and orientation ß ∗ , as well as density profiles as functions of distance from grafting plane ρ ( z ) , are studied. The effect of an added bond repulsion potential to suppress bond crossing in DPD is analyzed. Scaling laws of R g and its components R g ⊥ and R g ∥ are investigated. We find R g ∝ N ν , ν = 0.588 for surface coverages below the overlap surface concentration σ ∗ . For σ > σ ∗ we find R g ⊥ ∝ N ν ⊥ , ν ⊥ ≅ 1 and R g ∥ ∝ N ν ∥ , ν ∥ = 1 / 2 of ring brushes with the standard DPD model and ν ∥ ≅ 2 / 5 with added bond repulsion. The σ dependence of the radius of gyration was found to be R g ∝ σ µ with µ = 1 / 3 for surface coverages grater than σ ∗ . The perpendicular component R g ⊥ scales independent of the bond repulsion potential as R g ⊥ ∝ σ µ âŠ¥ , µ ⊥ = 1 / 3 , whereas the scaling of the parallel component exhibits a topological repulsion dependence R g ∥ ∝ σ µ âˆ¥ , µ ∥ = - 1 / 12 for standard DPD and µ ∥ = - 1 / 6 for bond repulsion.

Simulation of high-density water: its glass transition for various water models.

High-density amorphous water is simulated by use of isothermal-isobaric molecular dynamics at a pressure of 0.3 GPa making use of several water models (SPC/E, TIP3P, TIP4P variants, and TIP5P). Heating/cooling cycles are performed in the temperature range 80-280 K and quantities like density, total energy, and mobility are analysed. Raw data as well as the glass transition temperatures Tg observed in our studies depend on the water model used as well as on the treatment of intramolecular bonds and angles. However, a clear-cut evidence for the occurrence of a glass-to-liquid transition is found in all cases. Thus, all models indicate that high-density amorphous ice found experimentally may be a low-temperature proxy of an ultraviscous high-density liquid.

Shielding effects in polymer-polymer reactions. V. Concentration dependence of contact formation between star-branched and linear chains.

By use of the Dissipative Particle Dynamics (DPD) simulation technique mixtures of star-branched (arm number F = 4) and linear chains in athermal (good) solvent are analyzed regarding probabilities for intermolecular contacts of various reactive sites within different polymer coils. The accompanying sterical hindrances are described in the framework of shielding factors in order to investigate reactions and side reactions in radical polymerization and other techniques that involve polymer-polymer coupling. The shielding factors are studied as a function of total concentration from high dilution up to the bulk for different chain lengths of star-shaped and linear chains. Results indicate that their concentration dependence can be described by a power law for systems above the overlap concentration, whereas the chain length dependence vanishes when extrapolating to infinite chain lengths in that concentration range. Also the influence of the ratio of star chains and linear chains is studied for various concentrations.

Monte Carlo simulation studies of ring polymers at athermal and theta conditions.

By use of an intramolecular criterion, i.e., the direct proportionality between mean square dimension and chain length, theta conditions for linear chains and ring shaped polymers are evaluated for several types of cubic lattice chains (simple cubic, body centered cubic, and face centered cubic). The properties of the rings are evaluated for the same thermodynamic conditions under which they are prepared thus allowing for a natural amount of knots which have been identified by use of Alexander polynomials. For the limit of infinite chain lengths the same theta parameter is found for linear chains and rings. On the contrary, a significant theta point depression occurs due to an additional excluded volume effect if unknots are exclusively regarded. Parameters characteristic of the shape of rings and chains under theta conditions extrapolated to infinite chain length fairly well coincide with respective data for random walks. Mean square dimensions (characteristic of the size) of theta systems are slightly in excess as compared to nonreversal random walks due to the necessity of avoiding overlaps on a local scale. Furthermore athermal systems are studied as well for comparison; mean square dimensions are described by use of scaling relations with proper short chain corrections, shape parameters are given in the limit of infinite chain length.

Shielding effects in polymer-polymer reactions. IV. Intermolecular contact formation between star-branched molecules.

Pairs of star-branched molecules--taken from ensembles of athermal five-way cubic lattice chains prepared by Monte Carlo simulation--are analyzed for the relative probability of mutual contact formation between particularly specified segments i and j belonging to different chains within these pairs. These contact probabilities--termed shielding factors K(ij)--are calculated by means of exact enumeration as a function of chain length ranging from n = 8 to 256 bonds per arm, as a function of functionality (i.e., the number of arms) ranging from F = 2 to 6 arms, and as a function of segment position within the arms ranging from central to terminal segments. In addition, changes of properties that characterize the size and shape of the involved molecules while approaching and penetration are evaluated as a function of chain separation.

Molecular dynamics simulations of beta-cyclodextrin-aziadamantane complexes in water.

Force-field-based atomistic simulations of host-guest supramolecular complexes between beta-cyclodextrin and several aziadamantane derivatives have been analyzed with respect to relative orientation and interaction energies, explicitly considering solvent (water) molecules. For each case, the calculations revealed two stable orientations of the guest within the host that are different in interaction energy. Fluctuation of and correlation between characteristic properties were analyzed. Among other things, it turned out that orientation angle and inclusion depth are clearly correlated. In addition, for the unsubstituted aziadamantane, the enthalpy of complex formation was calculated and compared to experimental results.

Shielding effects in polymer-polymer reactions. II. Reactions between linear and star-branched chains with up to six arms.

The shielding effect of surrounding arms and chains on the encounter probability of reactive sites located both at the end of a linear chain and at several positions along the arms of star-branched chains with up to six arms is calculated by means of exact enumeration of samples prepared by Monte Carlo simulation. The changes of parameters that characterize the size and the shape of chain configurations during the approach of reactive centers located at the end of the linear chain and at the center of the star are evaluated. In addition to this specific case, which represents the central reaction step in reversible addition-fragmentation chain transfer star polymerization following the Z-group approach, a general discussion is given on the chain-length dependence of shielding factors associated with distinct segment positions.

Lactonisation--a degradation pathway for active pharmaceutical compounds: an in silico study in amorphous trehalose.

The lactonisation of a CCR1 inhibitor (CC chemokine receptor 1, involved in autoimmune diseases) featuring a hydroxyl group in a gamma-position (gamma-OH) with respect to an amide group has been investigated in silico. The two key steps of the lactonisation reaction are (i) rearrangement to an optimal conformation and (ii) the formation of the lactone (ring closure) and expulsion of NH3. Quantum chemical calculations in the gas phase were employed to identify conformers of the molecule with favorable starting geometries for a lactonisation reaction. In total, calculations of 1296 conformers revealed that it is energetically feasible for an inhibitor molecule to adopt a conformation where the carbon atom of the amide group (C(amide)) is suitably close to the oxygen atom of the gamma-OH (O(gamma)) to facilitate a successful lactonisation reaction. Additionally, molecular dynamics methods were used to show that rearrangement to a suitable conformer for lactonisation to occur happens to a lesser extent when the CCR1 inhibitor was embedded in an amorphous trehalose matrix (a model carbohydrate excipient). The mechanism of the actual lactonisation was investigated using the complete Linear Synchronous Transit/Quadratic Synchronous Transit (LST/QST) method. This was performed in both the gas phase and in water and was found to be a concerted reaction.

The glass transition temperatures of amorphous trehalose-water mixtures and the mobility of water: an experimental and in silico study.

Isothermal-isobaric molecular dynamics simulations are used to calculate the specific volume of models of trehalose and three amorphous trehalose-water mixtures (2.9%, 4.5% and 5.3% (w/w) water, respectively) as a function of temperature. Plots of specific volume versus temperature exhibit a characteristic change in slope when the amorphous systems change from the glassy to the rubbery state and the intersection of the two regression lines provides an estimate of the glass transition temperature T(g). A comparison of the calculated and experimental T(g) values, as obtained from differential scanning calorimetry, shows that despite the predicted values being systematically higher (about 21-26K), the trend and the incremental differences between the T(g) values have been computed correctly: T(g)(5.3%(w/w))

A computational investigation of the different intermediates during organoalkoxysilane hydrolysis.

Using a combination of atomistic molecular dynamics (MD) simulations and density functional theory (DFT) calculations, the four steps of hydrolysis of aminopropyl-, thiolpropyl-, and butyltrimethoxysilane have been studied. Large box MD simulations at constant pressure and temperature yield appropriate pair distribution functions--which allows us to quantify the number of surrounding water molecules--as well as the density of the systems. These densities serve as input for small box DFT calculations, which allow further geometry optimization and calculation of the electronic structure of the systems. The periodic DFT calculations are compared with gas-phase simulations. In all cases, the first step of hydrolysis is exothermic with the extent depending on the type of silane as well as on the number of hydrogen bonds in the initial stage.

Glass transition temperature of glucose, sucrose, and trehalose: an experimental and in silico study.

Isothermal-isobaric molecular dynamics simulations are used to calculate the specific volume of models of different amorphous carbohydrates (glucose, sucrose, and trehalose) as a function of temperature. Plots of specific volume vs temperature exhibit a characteristic change in slope when the amorphous systems change from the glassy to the rubbery state. The intersection of the regression lines of data below (glassy state) and above (rubbery state) the change in slope provides the glass transition temperature (T(g)). These predicted glass transition temperatures are compared to experimental T(g) values as obtained from differential scanning calorimetry measurements. As expected, the predicted values are systematically higher than the experimental ones (about 12-34 K) as the cooling rates of the modeling methods are about a factor of 10(12) faster. Nevertheless, the calculated trend of T(g) values agrees exactly with the experimental trend: T(g)(glucose) < T(g)(sucrose) < T(g)(trehalose). Furthermore, the relative differences between the glass transition temperatures were also computed precisely, implying that atomistic molecular dynamics simulations can reproduce trends of T(g) values in amorphous carbohydrates with high quality.

Adsorption of organosilanes at a Zn-terminated ZnO (0001) surface: molecular dynamics study.

Four different organosilanes (octyltrihydroxysilane, butyltrihydroxysilane, aminopropyltrihydroxysilane, and thiolpropyltrihydroxysilane) adsorbed at a reconstructed Zn-terminated polar ZnO (0001) surface are studied via constant temperature (298 K) molecular dynamics simulations. Both single adsorbed silane molecules as well as adsorbed silane layers are modeled, and the energy, distance, orientation, and alignment of these adsorbates are analyzed. The adsorbed silane molecules exhibit behavior depending on the chemical nature of their tail (nonpolar or polar) as well as on the silane concentration at the solid surface (single adsorption or silane layer). In contrast to the O-terminated ZnO surface studied previously, now adsorption can only occur at the vacancies of this reconstructed crystal surface, thus leading to an arched structure of the liquid phase near the crystal surface. Nevertheless, both nonpolar and polar single adsorbed silanes show a similar orientation and alignment at the surface (orthogonal in the former, parallel in the latter case) as for the O-terminated ZnO surface, although the interaction energy with the surface is considerably increased for nonpolar silanes while it is nearly unaffected for the polar ones. For adsorbed silanes within silane layers, the difference to single adsorbed silanes depends on the polarity of the tail: nonpolar silanes again show an orthogonal alignment, while polar silanes exhibit two different orientations at the solid surface-a head and a tail down configuration. This leads to two completely different but nevertheless stable orientations of these silanes at the Zn-terminated ZnO surface.

Atomistic molecular-dynamics simulations of the size and shape of polyethylene in hexane at infinite dilution.

Parameters characteristic of size and shape of single polyethylene chains consisting of 15-60 monomer units dissolved in hexane are calculated by use of molecular-dynamics simulations based on a fully atomistic representation of the system. Results are compared with corresponding calculations in vacuum as well as Monte Carlo simulations of coarse-grained chains. The major concern of the study is a careful check of actual limits and possibilities of atomistic simulations of global properties of polymers. As expected such simulations are still restricted to rather small chain lengths but are already large enough to obey the characteristics of polymer coils.

Multilayer adsorption of water at a rutile TiO2)(110) surface: towards a realistic modeling by molecular dynamics.

We present a model combining ab initio concepts and molecular dynamics simulations for a more realistic treatment of complex adsorption processes. The energy, distance, and orientation of water molecules adsorbed on stoichiometric and reduced rutile TiO(2)(110) surfaces at 140 K are studied via constant temperature molecular dynamics simulations. From ab initio calculations relaxed atomic geometries for the surface and the most probable adsorption sites were derived. The study comprises (i) large two-dimensional surface supercells, providing a realistically low concentration of surface oxygen defects, and (ii) a water coverage sufficiently large to model the onset of the growth of a bulk phase of water on the surface. By our combined approach the influence of both, the metal oxide surface, below, and the bulk water phase, above, on the water molecules forming the interface between the TiO(2) surface and the water bulk layer is taken into account. The good agreement of calculated adsorption energies with experimental temperature programmed desorption spectra demonstrates the validity of our model.

Monte Carlo investigations of dense copolymer systems. III. Properties of triblock copolymers in good and theta solvent.

The present article gives an analysis of XYX triblock copolymers in a good solvent and in a theta solvent, the segments of type X and type Y being repulsive for each other. The results are compared to homopolymers as well as to copolymers in a selective solvent that is a good one for the outer blocks and a theta solvent for the inner one and vice versa, the strength of repulsion between blocks being the same as in the present types of copolymers. A lattice model is used for the investigations and the concentration ranges from a volume fraction phi = 0 up to phi = 0.8. In the limit phi --> 0 the triblocks in good solvent are slightly more expanded than homopolymers and in theta solvent mean square dimensions of triblocks are considerably increased compared to homopolymers due to the repulsion between blocks. With increasing concentration the dimensions decrease but then they increase again and for large concentrations they become similar for all types of copolymers studied, as the effect of the solvent levels off making the repulsive interaction between blocks the dominant interaction. This leads to an orientation effect and as a consequence to microphase separation which is demonstrated by the concentration dependence of various quantities as well as by visualization of snapshots.

The initial stages of aminosilanol polymerisation.

For polysiloxanes to be used as a protective coating it is important that proton transfer, a trigger to polymerisation, is a facile process. Here we investigate the initial stages of polycondensation and compare different silanol tail groups and the effect of solvent (isopropanol). In the case of (3-aminopropyl)trihydroxysilane we see the potential for self catalysis as the tail group is a proton acceptor, while thiolpropyltrihydroxysilane and isopropanol do not promote proton transfer.