Aqueous solvation dynamics at metal oxide surfaces

Broadband transient absorption (TA) spectroscopy, three-pulse photon echo peak shift (3PEPS), and anisotropy decay measurements were used to study the solvation dynamics in bulk water and interfacial water at ZrO(2) surfaces, using Eosin Y as a probe. The 3PEPS results show a multiexponential behavior with two subpicosecond components that are similar in bulk and interfacial water, while a third component of several picoseconds is significantly lengthened at the interface. The bandwidth correlation function from TA spectra exhibits the same behavior, and the TA spectra are well reproduced using the doorway-window picture with the time constants from PEPS. Our results suggest that interfacial water is restricted to a thickness of less than 5 A. Also the high-frequency collective dynamics of water does not seem to be affected by the interface. On the other hand, the increase of the third component may point to a slowing down of diffusional motion at the interface, although other effects, may play a role, which are discussed.     

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.     

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.     

On the (im)possibility of evaluating correct individual rate constants of chain propagation kp and chain termination kt by combining kp2/kt and kp/kt data for chain‐length dependent termination

On the basis of simulated data two ways of evaluating individual rate constants by combining k_p²/k_t and k_p /k_t (k_p , k_t = rate constants of chain propagation and termination, respectively) were checked considering the chain‐length dependence of k_t. The first way tried to make use of the fact that pseudostationary polymerization yields data for k_p²/k_t as well as for k_p /k_t referring to the very same experiment, in the second way k_p²/k_t (from steady state experiments) and k_p/k_t data referring to the same mean length of the terminating radical chains were compared. In the first case no meaningful data at all could be obtained because different averages of k_t are operative in the expressions for k_p /k_t and k_p²/k_t. In spite of the comparatively small difference between these two averages (≈15% only) this makes the method collapse. The second way, which can be regarded as an intelligent modification of the “classical” method of determining individual rate constants, at least succeeded in reproducing the correct order of magnitude of the individual rate constants. However, although stationary and pseudostationary experiments independently could be shown to return the same k_t for the same average chain‐length of terminating radicals within extremely narrow limits no reasonable chain‐length dependence of k_t could be derived in this way. The reason is an extreme sensitivity of the pair of equations for k_p/k_t and k_p²/k_t towards small errors and inconsistencies which renders the method unsuccessful even for the high quality simulation data and most probably makes it even collapse for real data. This casts a characteristic light on the unsatisfactory situation with respect to individual rate constants determined in the classical way, regardless of a chain‐length dependence of termination. As a consequence, all efforts of establishing the chain‐length dependence of k_t are recommended to avoid this way and should rather resort to methods based on inserting a directly determined k_p into the equations characteristic of k_p²/k_t or k_p/k_t, properly considering the chain‐length dependent character of k_t.     

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.     

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.