Experimental investigation into the effect of stress on dissolution and growth of very soluble brittle salts in aqueous solution

Das Studium der Auflösungs- und Wachstumsprozesse an Feststoff-Flüssigkeits-Grenzflächen unter nicht-hydrostatischen Beanspruchungen ist wesentlich für das Verständnis von Defor-mationsprozessen, die in der Erde ablaufen. Unter diesen genannten Prozessen gehört die Drucklösung zu den wichtigsten duktilen Deformationsprozessen, von der Diagenese bishin zur niedrig- bis mittelgradigen metamorphen Bedingungen. Bisher ist allerdings wenig darüber bekannt, welche mechanischen, physikalischen oder chemischen Potentialenergie-Gradienten die Drucklösung steuern. I.a. wird angenommen, daß die Drucklösung durch Un-terschiede kristallplastischer Verformungsenergien oder aber durch Unterschiede der Normal-beanspruchung an Korngrenzen gesteuert wird. Unterschiede der elastischen Verformungs-energien werden dabei allerdings als zu gering erachtet, um einen signifikanten Beitrag zu leisten. Aus diesem Grund werden sie als mögliche treibende Kräfte für die Drucklösung vernachlässigt. Andererseits haben neue experimentelle und theoretische Untersuchungen gezeigt, daß die elastische Verformung in der Tat einen starken Einfluß auf Lösungs- und Wachstumsmechanismen von Kristallen in einer Lösung haben kann. Da die in der Erdkruste vorherrschenden Deformationsmechanismen überwiegend im elastischen Verformungsbereich der Gesteine ablaufen, ist es sehr wichtig, das Verständnis für die Effekte, die die elastische Verformung verursacht, zu erweitern, und ihre Rolle während der Deformation durch Drucklösung zu definieren. Die vorliegende Arbeit beschäftigt sich mit Experimenten, bei denen der Effekt der mechanisch kompressiven Beanspruchung auf Lösungs- und Wachstumsprozesse von Einzelkristallen unterschiedlicher, sehr gut löslicher, elastisch/spröder Salze untersucht wurde. Diese Salze wurden als Analoga gesteinsbildender Minerale wie Quarz und Calcit ausgewählt. Der Einfluß von Stress auf die Ausbildung der Oberflächenmikrostrukturen in einer untersättigten Lösung wurde an Kaliumalaun untersucht.Lösungsrillen (20 – 40 µm breit, 10 – 40 µm tief und 20 – 80 µm Abstand) entwickelten sich in den Bereichen, in denen die Beanspruchung im Kristall am größten war. Sie verschwanden wieder, sobald der Kristall entlastet wurde. Diese Rillen entwickelten sich parallel zu niedrig indizierten kristallographischen Richtungen und sub-perpendikular zu den Trajektorien, die der maximalen, lokalen kompressiven Beanspruchung entsprachen. Die Größe der Lösungsrillen hing von der lokalen Oberflächenbeanspruchung, der Oberflächenenergie und dem Untersättigungsgrad der wässrigen Lösung ab. Die mikrostrukturelle Entwicklung der Kristalloberflächen stimmte gut mit den theoretischen Vorhersagen überein, die auf den Modellen von Heidug & Leroy (1994) und Leroy & Heidug (1994) basieren. Der Einfluß der Beanspruchung auf die Auflösungsrate wurde an Natriumchlorat-Einzelkristallen untersucht. Dabei wurde herausgefunden, daß sich gestresste Kristalle schneller lösen als Kristalle, auf die keine Beanspruchung einwirkt. Der experimentell beobachtete Anstieg der Auflösungsrate der gestressten Kristalle war ein bis zwei Größenordnungen höher als theoretisch erwartet. Die Auflösungsrate stieg linear mit dem Stress an, und der Anstieg war um so größer, je stärker die Lösung untersättigt war. Außerdem wurde der Effekt der Bean-spruchung auf das Kristallwachstum an Kaliumalaun- und Kaliumdihydrogenphosphat-Ein-zelkristallen untersucht. Die Wachstumsrate der Flächen {100} und {110} von Kalium-alaun war bei Beanspruchung stark reduziert. Für all diese Ergebnisse spielte die Oberflächenrauhigkeit der Kristalle eine Schlüsselrolle, indem sie eine nicht-homogene Stressverteilung auf der Kristalloberfläche verursachte. Die Resultate zeigen, daß die elastische Verformung eine signifikante Rolle während der Drucklösung spielen kann, und eine signifikante Deformation in der oberen Kruste verursachen kann, bei Beanspruchungen, die geringer sind, als gemeinhin angenommen wird. Somit folgt, daß die elastische Bean-spruchung berücksichtigt werden muß, wenn mikrophysikalische Deformationsmodelle entwickelt werden sollen.

Weathering of Fe-bearing minerals under extraterrestrial conditions, investigated by Mössbauer spectroscopy

The weathering of Fe-bearing minerals under extraterrestrial conditions was investigated by Mössbauer (MB) spectroscopy to gain insights into the role of water on the planet Mars. The NASA Mars Exploration Rovers Spirit and Opportunity each carry a miniaturized Mössbauer spectrometer MIMOS II for the in situ investigation of Martian soils and rocks as part of their payload. The MER flight instruments had to be modified in order to work over the Martian diurnal temperature range (180 K – 290 K) and within the unique electronic environment of the rovers. The modification required special calibration procedures. The integration time necessary to obtain a good quality Mössbauer spectrum with the MIMOS II flight instruments was reduced by 30 % through the design of a new collimator. The in situ investigation of rocks along the rover Spirit's traverse in Gusev crater revealed weakly altered olivine basalt on the plains and pervasively altered basalt in the Columbia Hills. Correlation plots of primary Fe-bearing minerals identified by MB spectroscopy such as olivine versus secondary Fe-bearing phases such as nanophase Fe oxides showed that olivine is the mineral which is primarily involved in weathering reactions. This argues for a reduced availability of water. Identification of the Fe-oxyhydroxide goethite in the Columbia Hills is unequivocal evidence for aqueous weathering processes in the Columbia Hills. Experiments in which mineral powders were exposed to components of the Martian atmosphere showed that interaction with the atmosphere alone, in the absence of liquid water, is sufficient to oxidize Martian surface materials. The fine-grained dust suspended in the Martian atmosphere may have been altered solely by gas-solid reactions. Fresh and altered specimens of Martian meteorites were investigated with MIMOS II. The study of Martian meteorites in the lab helped to identify in Bounce Rock the first rock on Mars which is similar in composition to basaltic shergottites, a subgroup of the Martian meteorites. The field of astrobiology includes the study of the origin, evolution and distribution of life in the universe. Water is a prerequisite for life. The MER Mössbauer spectrometers identified aqueous minerals such as jarosite and goethite. The identification of jarosite was crucial to evaluate the habitability of Opportunity's landing site at Meridiani Planum during the formation of the sedimentary outcrop rocks, because jarosite puts strong constrains on pH levels. The identification of olivine in rocks and soils on the Gusev crater plains provide evidence for the sparsity of water under current conditions on Mars. Ratios of Fe2+/Fe3+ were obtained with Mössbauer spectroscopy from basaltic glass samples which were exposed at a deep sea hydrothermal vent. The ratios were used as a measure of potential energy for use by a microbial community. Samples from Mars analogue field sites on Earth exhibiting morphological biosignatures were also investigated.

Global Search for Minimum Energy (H2O)n Clusters, n = 3−5

The Gaussian-3 (G3) model chemistry method has been used to calculate the relative ΔG° values for all possible conformers of neutral clusters of water, (H2O)n, where n = 3−5. A complete 12-fold conformational search around each hydrogen bond produced 144, 1728, and 20 736 initial starting structures of the water trimer, tetramer, and pentamer. These structures were optimized with PM3, followed by HF/6-31G* optimization, and then with the G3 model chemistry. Only two trimers are present on the G3 potential energy hypersurface. We identified 5 tetramers and 10 pentamers on the potential energy and free-energy hypersurfaces at 298 K. None of these 17 structures were linear; all linear starting models folded into cyclic or three-dimensional structures. The cyclic pentamer is the most stable isomer at 298 K. On the basis of this and previous studies, we expect the cyclic tetramers and pentamers to be the most significant cyclic water clusters in the atmosphere.

Modeling of temperature and turbidity in a natural lake and a reservoir connected by pumped-storage operations

Pumped-storage (PS) systems are used to store electric energy as potential energy for release during peak demand. We investigate the impacts of a planned 1000 MW PS scheme connecting Lago Bianco with Lago di Poschiavo (Switzerland) on temperature and particle mass concentration in both basins. The upper (turbid) basin is a reservoir receiving large amounts of fine particles from the partially glaciated watershed, while the lower basin is a much clearer natural lake. Stratification, temperature and particle concentrations in the two basins were simulated with and without PS for four different hydrological conditions and 27 years of meteorological forcing using the software CE-QUAL-W2. The simulations showed that the PS operations lead to an increase in temperature in both basins during most of the year. The increase is most pronounced (up to 4°C) in the upper hypolimnion of the natural lake toward the end of summer stratification and is partially due to frictional losses in the penstocks, pumps and turbines. The remainder of the warming is from intense coupling to the atmosphere while water resides in the shallower upper reservoir. These impacts are most pronounced during warm and dry years, when the upper reservoir is strongly heated and the effects are least concealed by floods. The exchange of water between the two basins relocates particles from the upper reservoir to the lower lake, where they accumulate during summer in the upper hypolimnion (10 to 20 mg L−1) but also to some extent decrease light availability in the trophic surface layer.

The elusive S2 state, the S1/S2 splitting, and the excimer states of the benzene dimer

We observe the weak S 0 → S 2 transitions of the T-shaped benzene dimers (Bz)2 and (Bz-d 6)2 about 250 cm−1 and 220 cm−1 above their respective S 0 → S 1 electronic origins using two-color resonant two-photon ionization spectroscopy. Spin-component scaled (SCS) second-order approximate coupled-cluster (CC2) calculations predict that for the tipped T-shaped geometry, the S 0 → S 2 electronic oscillator strength f el (S 2) is ∼10 times smaller than f el (S 1) and the S 2 state lies ∼240 cm−1 above S 1, in excellent agreement with experiment. The S 0 → S 1 (ππ ∗) transition is mainly localized on the “stem” benzene, with a minor stem → cap charge-transfer contribution; the S 0 → S 2 transition is mainly localized on the “cap” benzene. The orbitals, electronic oscillator strengths f el (S 1) and f el (S 2), and transition frequencies depend strongly on the tipping angle ω between the two Bz moieties. The SCS-CC2 calculated S 1 and S 2 excitation energies at different T-shaped, stacked-parallel and parallel-displaced stationary points of the (Bz)2 ground-state surface allow to construct approximate S 1 and S 2 potential energy surfaces and reveal their relation to the “excimer” states at the stacked-parallel geometry. The f el (S 1) and f el (S 2) transition dipole moments at the C 2v -symmetric T-shape, parallel-displaced and stacked-parallel geometries are either zero or ∼10 times smaller than at the tipped T-shaped geometry. This unusual property of the S 0 → S 1 and S 0 → S 2 transition-dipole moment surfaces of (Bz)2 restricts its observation by electronic spectroscopy to the tipped and tilted T-shaped geometries; the other ground-state geometries are impossible or extremely difficult to observe. The S 0 → S 1/S 2 spectra of (Bz)2 are compared to those of imidazole ⋅ (Bz)2, which has a rigid triangular structure with a tilted (Bz)2 subunit. The S 0 → S 1/ S 2 transitions of imidazole-(benzene)2 lie at similar energies as those of (Bz)2, confirming our assignment of the (Bz)2 S 0 → S 2 transition.

Structure-based conformational preferences of amino acids

Proteins can be very tolerant to amino acid substitution, even within their core. Understanding the factors responsible for this behavior is of critical importance for protein engineering and design. Mutations in proteins have been quantified in terms of the changes in stability they induce. For example, guest residues in specific secondary structures have been used as probes of conformational preferences of amino acids, yielding propensity scales. Predicting these amino acid propensities would be a good test of any new potential energy functions used to mimic protein stability. We have recently developed a protein design procedure that optimizes whole sequences for a given target conformation based on the knowledge of the template backbone and on a semiempirical potential energy function. This energy function is purely physical, including steric interactions based on a Lennard-Jones potential, electrostatics based on a Coulomb potential, and hydrophobicity in the form of an environment free energy based on accessible surface area and interatomic contact areas. Sequences designed by this procedure for 10 different proteins were analyzed to extract conformational preferences for amino acids. The resulting structure-based propensity scales show significant agreements with experimental propensity scale values, both for α-helices and β-sheets. These results indicate that amino acid conformational preferences are a natural consequence of the potential energy we use. This confirms the accuracy of our potential and indicates that such preferences should not be added as a design criterion.

Peptide interfacial adsorption is kinetically limited by the thermodynamic stability of self association

We present a study of the adsorption of two peptides at the octane–water interface. The first peptide, Lac21, exists in mixed monomer–tetramer equilibrium in bulk solution with an appreciable monomer concentration. The second peptide, Lac28, exists as a tetramer in solution, with minimal exposed hydrophobic surface. A kinetic limitation to interfacial adsorption exists for Lac28 at moderate to high surface coverage that is not observed for Lac21. We estimate the potential energy barrier for Lac28 adsorption to be 42 kJ/mol and show that this is comparable to the expected free energy barrier for tetramer dissociation. This finding suggests that, at moderate to high surface coverage, adsorption is kinetically limited by the availability of interfacially active monomeric “domains” in the subinterfacial region. We also show how the commonly used empirical equation for protein adsorption dynamics can be used to estimate the potential energy barrier for adsorption. Such an approach is shown to be consistent with a formal description of diffusion–adsorption, provided a large potential energy barrier exists. This work demonstrates that the dynamics of interfacial adsorption depend on protein thermodynamic stability, and hence structure, in a quantifiable way.

Sorption and desorption behaviour of acetochlor in surface, subsurface and size-fractionated soil

Pesticides leaching through a soil profile will be exposed to changing environmental sorption and desorption conditions as different horizons with distinct physical and chemical properties are encountered. Soil cores were taken from a clay soil profile and samples taken from 0.0 to 0.3 m (surface), 1.0-1.3 m (mid) and 2.7-3.0 m (deep) and treated with the chloroacetanilide herbicide, acetochlor. Freundlich isotherms revealed that sorption and desorption behaviour varied with each depth sampled. As soil depth increased, the extent and strength of sorption decreased, indicating that the potential for leaching was increased in the subsoils compared with the surface soil. Hysteresis was evident at each of the three depths sampled, although no significant correlations between soil properties and the hysteresis coefficients were evident. Desorption studies using soil fractions with diameters of > 2000, 250-2000, 53-250, 20-53, 2-20, 0-2 and 0-1 mum separated from each of the three soil depths showed that differential desorption kinetics occurred and that the retention of acetochlor significantly correlated (R-2 = 0.998) with organic matter content. A greater understanding of the influence of soil components on the overall sorption and desorption potential of surface and subsurface soils is required to allow accurate prediction of acetochlor retention in the soil. In addition, it is likely that the proportion of each size fraction in a soil horizon would influence acetochlor bioavailability and movement to groundwater.

Comparisons of diffusive and viscous contributions to transport coefficients of light gases in single-walled carbon nanotubes

We examine here the relative importance of different contributions to transport of light gases in single walled carbon nanotubes, using methane and hydrogen as examples. Transport coefficients at 298 K are determined using molecular dynamics simulation with atomistic models of the nanotube wall, from which the diffusive and viscous contributions are resolved using a recent approach that provides an explicit expression for the latter. We also exploit an exact theory for the transport of Lennard-Jones fluids at low density considering diffuse reflection at the tube wall, thereby permitting the estimation of Maxwell coefficients for the wall reflection. It is found that reflection from the carbon nanotube wall is nearly specular, as a result of which slip flow dominates, and the viscous contribution is small in comparison, even for a tube as large as 8.1 nm in diameter. The reflection coefficient for hydrogen is 3-6 times as large as that for methane in tubes of 1.36 nm diameter, indicating less specular reflection for hydrogen and greater sensitivity to atomic detail of the surface. This reconciles results showing that transport coefficients for hydrogen and methane, obtained in simulation, are comparable in tubes of this size. With increase in adsorbate density, the reflection coefficient increases, suggesting that adsorbate interactions near the wall serve to roughen the local potential energy landscape perceived by fluid molecules.

Diffusion of n-decane in mesoporous MCM-41 silicas

We investigate here the diffusion of n-decane in nanoporous MCM-41 silicas with pore diameters between 3.0 and 4.3 nm, and at various temperatures and purge flow rates, by the Zero Length Column method. A complete-time-range analysis of desorption curves is used to derive the diffusion coefficient, and the effect of pore size, purge flow rate and temperature on the diffusion character is systematically studied. The results show that the calculated low-coverage diffusivity values are strongly dependent on temperature but only weakly dependent on pore size. The study reveals that transport is controlled by intracrystalline diffusion and dominated by sorbate-sorbent interaction, with the experimental isosteric heat matching the potential energy of flat-lying n-decane molecules on the surface, determined using a united atom model. The diffusion activation energy and adsorption isosteric heat at zero loading for the different pore size MCM-41 samples vary in a narrow range respectively, and their ratio is essentially constant over the pore size range studied. The study shows that the ZLC method is an effective tool to investigate the diffusion kinetics of hydrocarbons in mesoporous MCM-41 materials. (c) 2005 Elsevier Inc. All rights reserved.

Energy profiles for ketene cyclizations. interconversion of 1,3-oxazin-6-ones, mesoionic 1,3-oxazinium olates and acylketenes, imidoylketenes, oxoketenimines, and cyclization products

The energy surface connecting oxazinium olates 9, several possible conformers of ketenes 10 and 11, and the final cyclization products 12, 13 and 14, as well as the isomeric 1,3-oxazine-6-ones 15, ring opening of the latter to N-acylimidoylketenes 16, and subsequent rearrangement of 16 to oxoketenimines 17, azetinones 18, and the cyclization products 19 and 20 are evaluated computationally at the B3LYP/6-31G* and B3LYP/6-311+G*//B3LYP/6-31G* levels. The cyclizations of ketenes to oxazinium olates 9 and oxazines 15 have the characteristics of pseudopericyclic reactions. Plots of the energy vs internal reaction coordinate for the cyclization of transoid acylketenes such as 10 to 9 (via TS1) and 16 to 15 (via TS7) feature two inflection points and indicate that the part of the energy surface above the lower inflection points describe internal rotation of the acyl function in the ketene moiety, and the part below this point describes the cyclization of the cisoid ketene to the planar mesoionic oxazinium olate 9 or oxazinone 15. The 1,3-shifts of the OR group that interconvert ketenes 16 and ketenimines 17 via four-membered cyclic transition states TS8 behave similarly, the first portion (from the ketenimine side) of the activation barrier being due largely to internal rotation of substituents, and the top part being due to the 1,3-shift proper.

Adsorption in slit pores and pore-size distribution: A molecular layer structure theory

A new approach is developed to analyze the thermodynamic properties of a sub-critical fluid adsorbed in a slit pore of activated carbon. The approach is based on a representation that an adsorbed fluid forms an ordered structure close to a smoothed solid surface. This ordered structure is modelled as a collection of parallel molecular layers. Such a structure allows us to express the Helmholtz free energy of a molecular layer as the sum of the intrinsic Helmholtz free energy specific to that layer and the potential energy of interaction of that layer with all other layers and the solid surface. The intrinsic Helmholtz free energy of a molecular layer is a function (at given temperature) of its two-dimensional density and it can be readily obtained from bulk-phase properties, while the interlayer potential energy interaction is determined by using the 10-4 Lennard-Jones potential. The positions of all layers close to the graphite surface or in a slit pore are considered to correspond to the minimum of the potential energy of the system. This model has led to accurate predictions of nitrogen and argon adsorption on carbon black at their normal boiling points. In the case of adsorption in slit pores, local isotherms are determined from the minimization of the grand potential. The model provides a reasonable description of the 0-1 monolayer transition, phase transition and packing effect. The adsorption of nitrogen at 77.35 K and argon at 87.29 K on activated carbons is analyzed to illustrate the potential of this theory, and the derived pore-size distribution is compared favourably with that obtained by the Density Functional Theory (DFT). The model is less time-consuming than methods such as the DFT and Monte-Carlo simulation, and most importantly it can be readily extended to the adsorption of mixtures and capillary condensation phenomena.

Influence of adsorbate interaction on transport in confined spaces

This article provides a review of the recent theory of transport in nanopores developed in the author's laboratory. In particular the influence of fluid-solid interactions on the transport coefficient is examined, showing that such interactions reduce the value of the coefficient by almost an order of magnitude in comparison to the Knudsen theory for non-interacting systems. The activation energy and potential energy barriers for diffusion in smooth pores with a one-dimensional potential energy profile are also discussed, indicating the inadequacy of the commonly used assumption of proportionality between the activation energy and heat of adsorption or the minimum pore potential energy. A further feature affected by fluid-solid interactions is the nature of the reflection of fluid molecules colliding with a pore wall surface, varying from being nearly specular - such as in carbon nanotubes - to nearly diffuse for amorphous solids. Diffuse reflection leads to momentum loss and reduced transport coefficients. However, fluid-solid interactions do not affect the transport coefficient in the single-file diffusion regime when the surface reflection is diffuse, and the transport coefficient in this case is largely independent of the adsorbed density.

Steric stabilization and dissolution characteristics of aqueous suspensions of pharmaceutical interest

The adsorption of two qroups of nonionic surface active agents and a series of hiqh molecular weiqht hydrophilic polymer fractions onto a polystyrene latex and a drug substance diloxanide furoate B.P. has been investigated. The presence of pores within the drug surface has been demonstrated and this is shown to increase the adsorption of low molecular weight polymer species. Differences in the maximum amount of polymer adsorbed at both solid-solution interfaces have been ascribed to the different hydrophobicities of the surface as determined by contact angle measurements. Adsorbed layer thicknesses of polymer on polystyrene latex have been determined by three techniques: microelectrophoresis, intensity fluctuation spectroscopy and by viscometric means. These results, in combination with adsorption data, were used to interpret the configuration of the adsorbed polymer molecules at the interface. The type of druq suspension produced on adsorbing the different polymers in the absence of electrostatic stabilization was correlated with theoretical prediuctions of suspension characteristics deduced from potential energy diagrams, The agreement was good for the adsorption of short chain length surfactants, but for the polyvinylalcohols, discrepancies were found between experiment and theory. This was attributed to the inappropriate use of a mean segment density approximation within the adsorbed layer to calculate attractive potentials between particles. A maximum in the redispersibility values for suspensions coated with adsorbed nonylphenylethoxylates was attributed to "partial static stabilization" of the particles in conjunction with the attractive forces operating in the sediment between bare surface patches on neighbouring particles. No significant change in the dissolution of the drug was observed when nonylphenylethoxylates were adsorbed due to desorption upon contact with the dissolution medium. Pluronic F68 and all the polyvinylalcohol fractions caused a reduction in the dissolution rate which is explained by the decreased diffusion of drug' through the adsorbed polymer layer.

The effect of certain hydrophilic polymers on suspension stability

The adsorption of nonionic surface active agents of polyoxyethylene glycol monoethers of n hexadecanols on polystyrene latex and nonionic cellulose polymers of hydroxyethyl cellulose, hydroxypropyl cellulose and hydroxypropyl methylcellulose on polystyrene latex and ibuprofen drug particles have been studied. The adsorbed layer thicknesses were determined by means of microelectrophoretic and viscometric methods. The conformation of the adsorbed molecules at the solid-liquid interface was deduced from the molecular areas and the adsorbed layer thicknesses. Comparison of the adsorption results obtained from polystyrene latex and ibuprofen particles was made to explain the conformation difference between these two adsorbates. Sedimentation volumes and redispersibility values were the main criteria used to evaluate suspension stability. At low concentrations of surface active agents, hard caked suspensions were found, probably due to the attraction between the uncoated areas or, the mutual adsorption of the adsorbed molecules on the bare surface of the particles in the sediment. At high concentrations of hydroxypropyl cellulose and hydroxypropyl methylcellulose, heavily caked sediments were attributed to network structure formation by the adsorbed molecules. An attempt was made to relate the characteristics of the suspensions to the potential energy of interaction curves. Generally, the agreement between theory and experiment was good, but for hydroxyethyl cellulose-ibuprofen systems discrepancies were found. Experimental studies showed that hydroxyethyl cellulose flocculated polystyrene latex over a rather wide range of concentrations; similarly, hydroxyethyl cellulose-ibuprofen suspensions were also flocculated. Therefore, it ls suggested that a term to account for flocculation energy of the polymer should be added to the total energy of interaction. A rheometric method was employed to study the flocculation energy of the polymer.