Can the Gibbs free energy of adsorption be predicted efficiently and accurately: an M05-2X DFT study

This study presents new insight into the prediction of partitioning of organic compounds between a carbon surface (soot) and water, and it also sheds light on the sluggish desorption of interacting molecules from activated and nonactivated carbon surfaces. This paper provides details about the structure and interactions of benzene, polycyclic aromatic hydrocarbons, and aromatic nitrocompounds with a carbon surface modeled by coronene using a density functional theory approach along with the M05-2X functional. The adsorption was studied in vacuum and from water solution. The molecules studied are physisorbed on the carbon surface. While the intermolecular interactions of benzene and hydrocarbons are governed by dispersion forces, nitrocompounds are adsorbed also due to quite strong electrostatic interactions with all types of carbon surfaces. On the basis of these results, we conclude that the method of prediction presented in this study allows one to approach the experimental level of accuracy in predicting thermodynamic parameters of adsorption on a carbon surface from the gas phase. The empirical modification of the polarized continuum model leads also to a quantitative agreement with the experimental data for the Gibbs free energy values of the adsorption from water solution.     

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.     

Interactions of Aluminum(III) with Phosphates

In order to obtain information about aluminum(III)-phosphate interactions, potentiometric measurements were carried out to characterize the complex forming properties of Al(III) with organic phosphates, phosphonates, and nucleoside-5'-monophosphates. The aluminum(III)-orthophosphate system is difficult to study due to AlPO(4) precipitation. To overcome this problem, the stability constant logarithms of the 1:1 Al(III) complexes of ligands with the same donor groups (log K(1:1)) were plotted against the basicities of the ligands (log K(PO)3(H)). The resulting linear free energy relation (LFER) indicates that organic phosphates, phosphonates, and uridine-, thymidine-, and guanosine 5'-monophosphates similarly bind Al(III). Adenosine and cytidine 5'-monophosphate fall above the LFER owing to the presence of a second microform with the nucleic base protonated and a hydroxide bound to the Al(III). From the LFER the log stability constant for Al(III) binding to HPO(4)(2-) is estimated as 6.13 +/- 0.05. From the weakness of any soluble orthophosphate complexes of Al(III) we confirm the importance of citrate as the main small molecule Al(3+) binder in the blood serum. The study includes investigation of Al(III) binding to di- and triphosphates, which bind metal ion differently than monophosphates. Structures of the complexes were supported by (31)P NMR measurements.     

Accurate ab initio intermolecular potential energy surface for the quintet state of the O2(3Sigma(g)-)-O2(3Sigma(g)-) dimer

A new potential energy surface (PES) for the quintet state of rigid O(2)((3)Sigma(g)(-)) + O(2)((3)Sigma(g)(-)) has been obtained using restricted coupled-cluster theory with singles, doubles, and perturbative triple excitations [RCCSD(T)]. A large number of relative orientations of the monomers (65) and intermolecular distances (17) have been considered. A spherical harmonic expansion of the interaction potential has been built from the ab initio data. It involves 29 terms, as a consequence of the large anisotropy of the interaction. The spherically averaged term agrees quite well with the one obtained from analysis of total integral cross sections. The absolute minimum of the PES corresponds to the crossed (D(2d)) structure (X shape) with an intermolecular distance of 6.224 bohrs and a well depth of 16.27 meV. Interestingly, the PES presents another (local) minimum close in energy (15.66 meV) at 6.50 bohrs and within a planar skewed geometry (S shape). We find that the origin of this second structure is due to the orientational dependence of the spin-exchange interactions which break the spin degeneracy and leads to three distinct intermolecular PESs with singlet, triplet, and quintet multiplicities. The lowest vibrational bound states of the O(2)-O(2) dimer have been obtained and it is found that they reflect the above mentioned topological features of the PES: The first allowed bound state for the (16)O isotope has an X structure but the next state is just 0.12 meV higher in energy and exhibits an S shape.     

On the mechanism of nanopore filling of SAPO-5 molecular sieve by nitrogen molecules

SAPO-5 molecular sieve was synthesized according to patent literature and characterized with X-ray diffraction, electron microscopy (SEM, TEM), and solid state magic angle spinning NMR spectroscopy. The material of particles in the micrometer region was found to consist of ca. 20 nm microcrystallites packed in mostly parallel orientation to ca. 200 nm sized agglomerates. The nitrogen adsorption isotherm was measured at 77.6 K over ca. 7 decades of pressure up to pore saturation. The course of the isotherm is interpreted to consist of filling of the nanopores (diameter, 0.73 nm) up to 2N(2)/unit cell, subsequent multilayer adsorption on the outer surface of the agglomerates, and, finally, pore condensation in the interparticle adsorption space. The nanopore adsorption can be quantitatively reproduced with the statistical mechanical model of a quasi one-dimensional lattice gas taking intermolecular interactions into account. The evaluated energy parameters are of physically reasonable magnitude and agree with literature data. The multilayer part of the adsorption isotherm can be well represented by the Brunauer-Emmett-Teller model yielding a specific outer space area (63 m(2) g(-)(1)) which is consistent with estimated geometrical and pore size analysis data.     

Interactions between GaO4Al12(OH)24(H2O)12(7+) and cellulosic materials

The adsorption qualities of GaO(4)Al(12)(OH)(24)(H(2)O)(12)(7+), a polycation with ε-Keggin structure, and its stability in contact with anionic cellulosic materials, was investigated under different concentration and ionic strength conditions. The cellulosic materials employed were two different fully bleached fibre materials, carboxyl methyl cellulose (CMC), and a spin-coated cellulose model surface. As analytical techniques, pH-measurements, potentiometric titrations, ICP-OES, QCM-D, equilibrium calculations and Extended X-ray Absorption Fine Structure (EXAFS) were used. The adsorption is substantial and the addition of GaO(4)Al(12)(OH)(24)(H(2)O)(12)(7+) to a fibre suspension results in a rapid decrease in pH, followed by a small and slow increase in pH. This behaviour can be explained as due to a rapid and strong (log β>2) equilibrium adsorption of intact GaO(4)Al(12)(OH)(24)(H(2)O)(12)(7+) ions, followed by a slow, and minor, 3-8%, decomposition into different monomers. Alternative layer by layer adsorption of this ion, and CMC, on a spin-coated cellulose model surface constitutes further evidence for the strong interactions between the anionic cellulose materials and GaO(4)Al(12)(OH)(24)(H(2)O)(12)(7+). It is shown that the adsorption observed could not be described as due to an unspecific Donnan adsorption behaviour, neither of GaO(4)Al(12)(OH)(24)(H(2)O)(12)(7+) nor Ga and Al monomers, and specific surface complex formation is therefore discussed and applied. The (≡COO)(7)GaO(4)Al(12)(OH)(24)(H(2)O)(12) species found to explain the pH- and metal adsorption data should be considered strictly as a stoichiometric entity.     

Chromatographic models for the sorption of neutral organic compounds by soil from water and air.

The solvation parameter model is used to construct models for the estimation of the soil-water and soil-air distribution constants and to characterize the contribution of fundamental intermolecular interactions to the underlying sorption processes. Wet soil is shown to be quite cohesive and polar but relatively non-selective for dipole-type, lone-pair electron and hydrogen-bond interactions. Using a comparison of system constant ratios chromatographic systems employing reversed-phase liquid chromatography on polar bonded phases are shown to provide suitable models for estimating soil-water distribution constants. No suitable gas chromatographic models were found for the soil-air distribution constant but the requirements for such a system are indicated. Models are also provided for adsorption at the air-water interface. Estimation methods based on either the solvation parameter model or chromatographic model reproduce experimental distribution constants for a wide variety of compounds with a similar error (0.2-0.3 log units) to that expected in the experimental data.     

Intermolecular interactions in (arene)chromium carbonyl compounds: prediction of chiral crystal packing from racemate structure

Six X-ray crystal structures are reported, all containing substituted triphenylmethanol derivative 4 either alone or as its mono or bis(chromium tricarbonyl) complexes. All four chromium complexes studied crystallize with two independent molecules in the crystallographic asymmetric unit. It is demonstrated that from the X-ray crystal structure of the acentric racemic (+/-)-(1pR,1' 'R)(1pS,1' 'S)-[Cr(CO)(3)(eta(6)-t-BuC(6)H(3)(CMeOMe)CPh(2)OH)], (+/-)-3, it is possible to deduce the 4-fold helical structure of the chiral (-)-(1pR,1' 'R) isomer, (-)-3. The bimetallic derivatives demonstrate the ability to control intermolecular interactions by the positioning of relative stereochemistry.     

Experimental observations on the scaling of adsorption isotherms for nonionic surfactants at a hydrophobic solid-water interface

The self-assembly of nonionic surfactants in bulk solution and on hydrophobic surfaces is driven by the same intermolecular interactions, yet their relationship is not clear. While there are abundant experimental and theoretical studies for self-assembly in bulk solution and at the air-water interface, there are only few systematic studies for hydrophobic solid-water interfaces. In this work, we have used optical reflectometry to measure adsorption isotherms of seven different nonionic alkyl polyethoxylate surfactants (CH3(CH2)I-1(OCH2CH2)JOH, referred to as CIEJ surfactants, with I = 10-14 and J = 3-8), on hydrophobic, chemically homogeneous self-assembled monolayers of octadecyltrichlorosilane. Systematic changes in the adsorption isotherms are observed for variations in the surfactant molecular structure. The maximum surface excess concentration decreases (and minimum area/molecule increases) with the square root of the number of ethoxylate units in the surfactant (J). The adsorption isotherms of all surfactants collapse onto the same curve when the bulk and surface excess concentrations are rescaled by the bulk critical aggregation concentration (CAC) and the maximum surface excess concentration. In an accompanying paper we compare these experimental results with the predictions of a unified model developed for self-assembly of nonionic surfactants in bulk solution and on interfaces.     

Diffuse reflectance infrared Fourier transform spectroscopy as a tool to characterise water in adsorption/confinement situations

We present experimental data acquired by diffuse reflectance infrared spectroscopy in the mid-IR (4000-400 cm(-1)), on micrometric-sized mineral grain powders. The spectral evolution of the OH-stretching band is followed when the adsorbed water film is thinned under dry conditions, from high to low hydration states. The IR bands are found to be characteristic of the degree of adsorption/confinement of the liquid water. The OH-stretching band is shifted toward shorter wavenumbers than in bulk water, showing that a significant portion of adsorbed water has a higher intermolecular bonding energy. Complementary treatment of the kinetics of water desorption, varying with the surface forces in the water film, confirms the relationships of these bands with the constrained water state. We distinguish different water types obeying liquid-liquid interactions (free and capillary water) or dominated by solid-water interactions (confined and adsorbed water). Part of this study is devoted to mesoporous silica MCM-41, of interest due to the restricted geometries of its mesopores (4.7 nm) favouring the confined water state. The methodology allows us to distinguish bulk and adsorbed/confined water, using spectral analysis coupled with an understanding of the dynamic behaviour of the desorption process.     

Molecular Structure and Dynamics of Water at the Water–Air Interface Studied with Surface‐Specific Vibrational Spectroscopy

Water interfaces provide the platform for many important biological, chemical, and physical processes. The water–air interface is the most common and simple aqueous interface and serves as a model system for water at a hydrophobic surface. Unveiling the microscopic (<1 nm) structure and dynamics of interfacial water at the water–vapor interface is essential for understanding the processes occurring on the water surface. At the water interface the network of very strong intermolecular interactions, hydrogen‐bonds, is interrupted and the density of water is reduced. A central question regarding water at interfaces is the extent to which the structure and dynamics of water molecules are influenced by the interruption of the hydrogen‐bonded network and thus differ from those of bulk water. Herein, we discuss recent advances in the study of interfacial water at the water–air interface using laser‐based surface‐specific vibrational spectroscopy.     

Insights into electrocatalysis

When standard reversible potentials for bulk solution reactions, U(0), are known, the reversible potentials when the reactant and product are adsorbed on an electrocatalyst surface, U(surf)(rev), are given in terms of these potentials and the adsorption Gibbs energy bond strengths: U(surf)(rev) = U0 + D(ads)G (Ox)/F-Δ(ads)G (R)/F (i). When the Δ(ads)G (Ox) and Δ(ads)G (Red) values are known at potential U, this equation is exact. When the overpotential for a multi-electron transfer reaction is minimal, each electron transfer takes place at the standard reversible potential for the overall reaction. In the case of O(2) reduction to water via the intermediate step OOH(aq) → O(aq) + OH(aq), or via O(2)(g) → 2O(aq), the respective endergonic O-O dissociation Gibbs energies are shown to be 2.52 eV and 4.76 eV. When the oxygen product and water reactant adsorb weakly, as on platinum, the adsorption Gibbs energies, Δ(ads)G, for O, OH, and OOH intermediates can be uniquely predicted using these data. All of the above depend exclusively on experimentally determined data. Reversible potentials have been calculated for oxygen reduction steps on the platinum electrocatalyst surface using Interface 1.0, a comprehensive computational code for the potential dependence of the electrochemical interface. Using these results as benchmarks, is found to be accurate to around 0.1 V when the Δ(ads)G are values calculated for the potentials of zero charge, instead of 1.229 V, which is a significant simplification. The variation in Δ(ads)G values between the calculated potentials of zero charge and 1.229 V are found to be 0.2 eV V(-1) or less. Prior work, using internal adsorption energies calculated at the potential of zero charge in place of Gibbs energies in was found to be accurate to within about 0.2 V. On platinum Δ(ads)G of the reaction OOH(ads) → O(ads) + OH(ads) is calculated at the potential of zero charge for the reactant and product to be about 1.2 eV exergonic under Langmuir conditions, and this Gibbs energy loss reduces the 1.229 V four-electron reversible potential on the platinum surface to an effective reversible potential of about 0.93 V for this mechanism on platinum. The effective reversible potential is a consequence of efficiency loss, not kinetics. Based on these values, the onset potential for four-electron oxygen reduction will be less than or equal to the effective reversible potential and on pure Pt(111) it appears to be equal to it.     

Chemistry of NO2 on oxide surfaces: formation of NO3 on TiO2(110) and NO2<-->O vacancy interactions

Synchrotron-based high-resolution photoemission, X-ray absorption near-edge spectroscopy, and first-principles density functional (DF) slab calculations were used to study the interaction of NO(2) with a TiO(2)(110) single crystal and powders of titania. The main product of the adsorption of NO(2) on TiO(2)(110) is surface nitrate with a small amount of chemisorbed NO(2). A similar result is obtained after the reaction of NO(2) with polycrystalline powders of TiO(2) or other oxide powders. This trend, however, does not imply that the metal centers of the oxides are unreactive toward NO(2). An unexpected mechanism is seen for the formation of NO(3). Photoemission data and DF calculations indicate that the surface nitrate forms through the disproportionation of NO(2) on Ti sites (2NO(2,ads) --> NO(3,ads) + NO(gas)) rather than direct adsorption of NO(2) on O centers of titania. Complex interactions take place between NO(2) and O vacancies of TiO(2)(110). Electronic states associated with O vacancies play a predominant role in the bonding and surface chemistry of NO(2). The adsorbed NO(2), on its part, affects the thermochemical stability of O vacancies, facilitating their migration from the bulk to the surface of titania. The behavior of the NO(2)/titania system illustrates the importance of surface and subsurface defects when using an oxide for trapping or destroying NO(x)() species in the prevention of environmental pollution (DeNOx operations).     

Palladium-catalyzed alpha-arylation of esters

A new and simple one-pot procedure for the palladium-catalyzed intermolecular alpha-arylation of esters is described. A number of esters can be functionalized with a wide range of aryl bromides using Pd(OAc)(2) or Pd(2)(dba)(3) and bulky electron-rich o-biphenyl phosphines 1-3. Under the reaction conditions, using LiHMDS as base, alpha-arylation proceeds at room temperature or at 80 degrees C with very good yields and high selectivities for monoarylation. Important nonsteroidal antiinflammatory drug derivatives such as (+/-)-naproxen tert-butyl ester and (+/-)-flurbiprofen tert-butyl ester can be prepared in 79% and 86% yield, respectively. The catalyst system based on the di-tert-butylphosphine (2) is also active for the alpha-arylation of esters using aryl chlorides. Furthermore, using (3) the alpha-arylation of trisubstituted ester enolates can be accomplished to provide compounds that have quaternary centers.     

Mechanistic aspects of nitrogen-heterocyclic compound adsorption on bamboo charcoal

Adsorption of three nitrogen-heterocyclic compounds (NHCs), pyridine, indole, and quinoline, on a new porous carbonaceous adsorbent, bamboo charcoal (BC), is investigated. Different structures and surface properties of BC are created by HNO(3) and NaOH treatment as well as by microwave (MW) radiation. HNO(3) treatment decreases surface area, micropore volume, and surface hydrophobicity, whereas NaOH treatment increases surface area and micropore volume. MW treatment dramatically increases hydrophobicity and fraction of aromatic structure. Adsorption isotherms of NHCs are nonlinear and better fitted by Freundlich model (FM) compared with Langmuir model (LM) and Polanyi-Manes model (PMM). The maximum adsorption capacities for pyridine, indole, and quinoline reach 42.92, 93.24, and 91.74mg g(-1), respectively, at an initial concentration of 200mg L(-1). Surface area, hydrophobicicty, and electrostatic and π-π electron-donor-acceptor (EDA) interactions are accountable for NHC adsorption. A model relating NHC adsorption (log K) and adsorbent-adsorbate physicochemical properties is developed to measure the relative contribution of these interactions giving the sequence of surface area>hydrophobic interaction>electrostatic interaction>π-π EDA interaction.     

Surface-active properties of new cationic gemini surfactants with cyclic spacer

Three cationic gemini surface active compounds of the type (1r,4r)-1,4-dialkyl-1,4-dimethy-l-piperazine-1,4-diium bromide (Ia, Ib, and Ic), were synthesized. They were characterized using elemental analysis and ¹H-NMR spectra. Their surface-active properties were measured in aqueous solutions with different concentrations at different temperatures (25, 40, and 55°C). Various surface measurements of these gemini surfactants, (compared to the conventional one, 1-Dodecyl-1-methylpiperidinium bromide (a)) were estimated, specifically critical micelle concentration (CMC), effectiveness (π_CMC), efficiency (P_C20) as well as maximum surface excess (Γ_max) and minimum surface area (A_min). The measurements of the gemini compounds gave low CMC, high efficiency in reducing the surface tension, and intense adsorption at air/water interface. These surfactants have lower Krafft points and thus better solubility. Thermodynamic data, free energy, entropy, and enthalpy changes (ΔG°, ΔS°, and ΔH°) for micellization at the air/water interface and also for adsorption in the bulk of surface-active solutions were calculated.     

The effect of solvent polarity on the balance between charge transfer and non-charge transfer pathways in the sensitization of singlet oxygen by pipi triplet states

A large set of literature kinetic data on triplet (T(1)) sensitization of singlet oxygen by two series of biphenyl and naphthalene sensitizers in solvents of strongly different polarity has been analyzed. The rate constants and the efficiencies of singlet oxygen formation are quantitatively reproduced by a model that assumes the competition of a non-charge transfer (nCT) and a CT deactivation channel. nCT deactivation occurs from a fully established spin-statistical equilibrium of (1)(T(1)(3)Sigma) and (3)(T(1)(3)Sigma) encounter complexes by internal conversion (IC) to lower excited complexes that dissociate to yield O(2)((1)Sigma(g)(+)), O(2)((1)Delta(g)), and O(2)((3)Sigma(g)(-)). IC of (1,3)(T(1)(3)Sigma) encounter complexes is controlled by an energy gap law that is generally valid for the transfer of electronic energy to and from O(2). (1,3)(T(1)(3)Sigma) nCT complexes form in competition to IC (1)(T(1)(3)Sigma) and (3)(T(1)(3)Sigma) exciplexes if CT interactions between T(1) and O(2) are important. The rate constants of exciplex formation depend via a Marcus type parabolic model on the corresponding free energy change DeltaG(CT), which varies with sensitizer triplet energy, oxidation potential, and solvent polarity. O(2)((1)Sigma(g)(+)), O(2)((1)Delta(g)), and O(2)((3)Sigma(g)(-)) are formed in the product ratio (1/6):(1/12):(3/4) in the CT deactivation channel. The balance between nCT and CT deactivation is described by the relative contribution p(CT) of CT induced deactivation calculated for a sensitizer of known triplet energy from its quenching rate constant. It is shown how the change of p(CT) influences the quenching rate constant and the efficiency of singlet oxygen formation in both series of sensitizers. p(CT) is sensitive to differences of solvent polarity and varies for the biphenyls and the naphthalenes as sigmoidal with DeltaG(CT). This quantitative model represents a realistic and general mechanism for the quenching of pipi triplet states by O(2), surpassing previous advanced models.     

Design of surface active soluble peptide molecules at the air/water interface

Surface active molecules collect at interfaces and have the potential to be used for water evaporation reduction. The objective of this work is to design surface active soluble peptides that collect at the air/water interface using molecular simulations. Rotational isomeric state Monte Carlo (RISMC) sampling together with a solvation model that we recently invented, the AAD solvation model [Gu, C.; Lustig, S.; Trout, B. J. Phys. Chem. B 2006, 110 (3), 1476-1484] was applied to calculate the adsorption free energy of the peptide molecule at the air/water interface. The results were validated by both molecular dynamics simulations with an explicit solvent model and surface tension measurements on synthesized peptides. It was demonstrated that this approach is able to give a reasonable prediction of surface activity with an approximately 50% hit rate in terms of designed surface active molecules actually being surface active. The relationship between the chemical composition and the surface morphology is also discussed.     

Surface rheology and foaming properties of sodium oleate and C12(EO)6 aqueous solutions

The dynamic surface tension (DST) and the surface viscoelastic modulus of sodium oleate aqueous solutions at different concentrations were measured using an image analysis tensiometer based on the oscillating bubble technique. The diffusion coefficient of oleate moieties was calculated from DST measurements and the surface viscoelastic modulus using the Langmuir-Szyszkowski and the diffusion-controlled adsorption models. The viscoelastic moduli obtained from model calculations were compared with the corresponding experimental values. The diffusion coefficient of C(12)(EO)(6) in water and the parameters of the Langmuir-Szyszkowski adsorption isotherm were taken from the literature and used to calculate the surface viscoelastic modulus of its aqueous solutions at different concentrations. The foaming properties of both C(12)(EO)(6) and sodium oleate solutions, viz., the foam conductance and the water volume fraction in the foam, were measured using a commercial Foamscan device. Foaming experiments with C(12)(EO)(6) and sodium oleate solutions were carried out either under static conditions; i.e., the foam conductance and the water volume fraction were measured as a function of time after the generation of a fixed volume of foam, or under dynamic conditions; i.e., the foam conductance and the water volume fraction were measured during foam formation. The variations in the foam permeability as a function of surfactant concentration were related to the viscoelastic properties of the air/water interface and to the presence of micelles in the foam films. With foams in which the water volume fraction was higher than 0.05, the foam electrical conduction could be described using a simple parallel resistor model and their conductance measurements were related to the foam water volume fraction. The results related to water drainage under static conditions were used to interpret water drainage under dynamic conditions. Preliminary conjectures on the influence of foam permeability and water volume fraction on the yield of the flotation deinking process were drawn from these results.