While the attempts currently in progress in several groups for the rigorous inclusion of dispersion interactions in density functional theory (DFT) calculations mature and evolve into practical methodology, we contribute to the debate on the applicability of current functionals to the calculation of weak interaction with a systematic investigation of a few, typical, weakly bound systems. We have used both pure DFT and a hybrid approach in which the total interaction energy is partitioned into two parts: (a) the dispersion energy which, in a first approximation is the contribution due to intermonomer correlations and (b) all other interactions. The first component is accurately obtained at all distances of interest by means of a well-known damped multipolar expansion of the dispersion energy while for the second component different approximations will be evaluated. The need to avoid double counting a fraction of the correlation energy when using the hybrid approach and the choice of the appropriate functional are also discussed. We consider four systems of increasing binding strength, namely the Ar2 and Kr2 dimers, the benzene dimer, the water dimer, and a few metal carbonyls. For pure DFT calculations we confirm the conclusion reached by others concerning (a) the strong dependence of the results on the choice of the GGA functional for dispersion-dominated interaction (noble gases and benzene) with the overall tendency to yield underbinding and (b) the relatively accurate, functional-independent, description for that DFT gives of water, which we attribute to the fact that this system is dominated by electrostatic interactions. For the carbonyls we find that DFT yields results which area again strongly dependent on the choice of the functional and show a tendency to give overbinding. Our hybrid method shows instead shortcomings only for the noble gases. The problem in this case is traceable to the well-known difficulties that all current functionals experience at medium–large intermonomer separations. The quality of the hybrid results improves markedly for benzene due to the large value of both dispersion and repulsive interactions at the equilibrium distance for this dimer, which makes the balance between the two, less delicate. Excellent results are also obtained for water (for the same reason as indicated above) and more significantly for the carbonyls where we find that dispersion contributes to the binding more than it could be guessed a priori. We do not claim to have found a general solution to this difficult problem, but we aim at providing a quantitative assessment to where the problems are pointing at directions from which a general solution may, eventually, emerge.

The origin of the superlattice present in the commensurate hexagonal structure of self-assembled monolayers of n-alkanethiols on gold and the question of whether the thiols are adsorbed onto the surface as dimers (disulfides) or monomers (thiolates) have been under debate for many years. Looking for a better understanding of the structural properties of these systems, we have performed a theoretical study of the molecular and dissociative adsorption of dimethyl disulfide on Au(111) as a function of coverage (0.25 ≤ ϑ ≤ 1), using gradient-corrected density functional (DFT) calculations with a slab geometry. For the dissociated state, our results indicate that the hcp hollow site is much less favorable than the fcc site. For the latter site, we find that, because of surface gold atom relaxation, the adsorption energy depends strongly on ϑ, changing from ∼18 kcal/mol at ϑ = 0.25 to ∼3 kcal/mol at ϑ = 1. For the bridge site, instead, the adsorption energy is a weak function of ϑ, and for all investigated coverages, this site is by far the most stable. According to our DFT approach, the adsorption of dimethyl disulfide is dissociative with a thermodynamic gain, at ϑ = 1, of ∼13 kcal/mol with respect to the adsorbed molecular state. We also find, however, that the energy of c(4 × 2) structures containing at least two inequivalent CH3S groups per unit cell (with a minimum S−S distance of ∼3.7 Å) is, within the accuracy of our approach, indistinguishable from the pure ( × ) hexagonal structure. Our results suggest that the full solution of this thorny problem will require, also for the shortest chains, an estimate of the energetic contribution of dispersion forces that are not included in the DFT calculations.

We fluidize a granular bed in a rectangular container by injecting energy through the lateral walls with high-frequency sinusoidal horizontal vibrations. In this way, the bed is brought to a steady state with no convection. We measured buoyancy forces on light spheres immersed in the bed and found that they obey Archimedes' principle. The buoyancy forces decrease when we reduce the injected energy. By measuring ascension velocities as a function of gamma, we can evaluate the frictional drag of the bed; its exponential dependence agrees very well with previous findings. Rising times of the intruders ascending through the bed were also measured, they increase monotonically as we increase the density.

A theoretical study of the hydration of arsenious acid is presented. This study included ab initio calculations and Monte Carlo simulations. The model potentials used for the simulations were ab initio derived and they include polarizability, nonadditivity, and molecular relaxation. It is shown that with these refined potentials it is possible to reproduce the available experimental evidence and therefore permit the study of clusters, as well as of the hydration process in solution. From the study of stepwise hydration and the Monte Carlo simulation of the condensed phase it is concluded that As(OH)(3) presents a hydration scheme similar to an amphipathic molecule. This phenomenon is explained as due to the existence of both a positive electrostatic potential and a localized lone pair in the vicinity of As. These results are used to rationalize the known passage of As(OH)(3) through aqua-glyceroporines.

We performed X-ray absorption spectroscopy measurements on the arsenic K-edge of As(III) in solution under acidic conditions. Extended X-ray absorption fine structure (EXAFS) and X-ray near edge structure (XANES) spectra were compared with theoretical calculations which use local atomic structure configurations, either derived from density functional theory (DFT) energy minimization (EM) calculations or based on classical Monte Carlo (MC) simulations, for a As(OH)3 cluster surrounded by water molecules. The nearest arsenic-oxygen distances obtained from the fit of the XAFS spectra are consistent with the distances present in configurations derived from Monte Carlo simulations but not with those obtained from DFT-EM calculations. Calculations of XANES using either DFT-EM or the average configuration obtained from MC simulations do not reproduce the XANES spectra in the vicinity of the absorption edge. However, specific local atomic structural configurations of the As(OH)3 and water molecules, obtained from MC simulations, which show some ordering of water molecules up to 5 Å from the arsenic, reproduce qualitatively the experimental spectra. These results highlight the capability of XANES to yield information about hydration of ions in solution.