Albeiro Restrepo

Structural Studies of the Water Hexamer

In this paper we report the geometries and properties of 24 structural isomers located on the MP2/6-311++g** potential energy surface of the water hexamer. At least 15 structural patterns are located within 3 kcal/mol of the most stable conformation, leading to a very complex potential energy surface, several isomers having significant contributions. A quadratic correlation between the distance from the proton to the center of the hydrogen bond with the distance between oxygen atoms for all clusters is reported. MP2/6-311++g** and CCSD(T)/aug-cc-pvdz//MP2/6-311++g** predict different stabilization orderings but are in good agreement for binding energies. Compact structures are energetically favored by electronic energies with zero point energy corrections, while noncompact cyclic structures are preferred when temperature and entropy are accounted for.

Stochastic Search of the Quantum Conformational Space of Small Lithium and Bimetallic Lithium-Sodium Clusters

In this paper we report the results obtained by an implementation and application of the simulated annealing optimization procedure to the exploration of the conformational space of small neutral and charged lithium clusters (Li(n)(q), n = 5, 6, 7; q = 0, +/-1) and of the bimetallic lithium/sodium clusters (Li5Na) in their lowest spin states. Our methodology eliminates the structure guessing procedure in the process of generating cluster configurations. We evaluate the quantum energy, typically with the Hartree-Fock Hamiltonian, of randomly generated points in the conformational space and use a modified Metropolis test in the annealing algorithm to generate candidate structures for atomic clusters. The structures are further optimized by analytical methods (gradient following) at the Møller-Plesset second order perturbation theory level (MP2), in conjunction with basis sets including polarization functions with and without diffuse functions. High accuracy ab initio energies at the coupled clusters level, with single, double, and triple substitutions from the Hartree-Fock determinant (CCSD(T)), on the MP2 geometries were calculated and used to establish the relative stability of the isomers within each potential energy surface. Various cluster properties were computed and compared to existing values in order to validate our methods. Our results show excellent agreement with previous experimental and theoretical reports. Even at these small sizes, evidence for 10 new structures never reported before for the lithium clusters and four new structures for the bimetallic clusters is presented.

Structures, energies and bonding in neutral and charged Li microclusters

Structural and chemical properties of charged and neutral Lithium microclusters are investigated for

Structures, energies, and bonding in the water heptamer

{text{Li}}_n^q(n = {5} - {1}0,q = 0,pm {1})

Hydrophobic meddling in small water clusters

. A total of 18 quantum conformational spaces are randomly walked to produce candidate structures for local minima. Very rich potential energy surfaces are produced, with the largest structural complexity predicted for anionic clusters. Analysis of the electron charge distributions using the quantum theory of atoms in molecules (QTAIM) predicts major stabilizing roles of Non–nuclear attractors (NNAs) via NNA···Li interactions with virtually no direct Li···Li interactions, except in the least stable configurations. A transition in behavior for clusters containing more than seven nuclei is observed by using the recently introduced quantum topology to determine in a quantum mechanically consistent fashion the number of spatial dimensions each cluster has. We experiment with a novel scheme for extracting persistent structural motifs with increase in cluster size. The new structural motifs correlate well with the energetic stability, particularly in highlighting the least stable structures. Quantifying the degree of covalent character in Lithium bonding independently agrees with the observation in the transition in cluster behavior for lithium clusters containing more than seven nuclei. Good correlation with available experimental data is obtained for all properties reported in this work.

Electronic Structure Calculations on Helical Conducting Polymers

We present an exhaustive stochastic search of the quantum conformational spaces of the (CH(3)O)(2)PO(2)(-) + nH(2)O (n = 1,2,3) systems. We uncover structural, conformational and energetic features of the problem. As in the isolated species, clusters containing the gauche-gauche (gg) conformation of dimethylphosphate (DMP(-)) are energetically preferred, however, contributions from hydrated gauche-anti (ga) and anti-anti (aa) monomers cannot be neglected because such structures are quite common and because they are close in energy to those containing the gg monomer. At least seven distinct types of O∙∙∙H-O-H contacts lead to DMP(-) ↔ water interactions that are always stabilizing, but not strong enough to induce significant changes in the geometries of either DMP(-) or water units. Our results lead us to postulate DMP(-) to be a suitable model to study explicit and detailed aspects of microsolvation of cell membranes.

The rich and complex potential energy surface of the ethanol dimer

In this paper we report the geometries and properties of 38 distinct geometrical motifs located on the B3LYP/6-31+G(d), MP2/6-311++G(d, p) potential energy surfaces of the water heptamer. Binding energies of up to 45 kcal/mol are calculated. All motifs fall within 10 kcal/mol of the most stable conformation, with at least 13 structural patterns located no more than 3 kcal/mol above, leading to a very complex potential energy surface, populated by a multitude of motifs each one allowing large numbers of conformations. Cluster stability does not seem to be correlated with the number of hydrogen bonds. Compact structures are energetically favored by electronic energies with zero-point energy corrections, while more open structures are preferred when temperature and entropy are accounted for. The molecular interactions holding the clusters as discrete units lead to large binding energies but are not strong enough to cause significant changes in the geometries of the interacting monomers. Our results indicate that bonding in the water heptamers can be considered as largely non-shared interactions with contributions from intermediate character of increasing covalency.

A detailed look at the reaction mechanisms of substituted carbenes with water.

In this work, we report a detailed study of the microsolvation of anionic ibuprofen, Ibu(-). Stochastic explorations of the configurational spaces for the interactions of Ibu(-) with up to three water molecules at the DFT level lead to very rich and complex potential energy surfaces. Our results suggest that instead of only one preponderant structure, a collection of isomers with very similar energies would have significant contributions to the properties of the solvated drug. One of these properties is the shift on the vibrational frequencies of the asymmetric stretching band of the carboxylate group in hydrated Ibu(-) with respect to the anhydrous drug, whose experimental values are nicely reproduced using the weighted contribution of the structures. We found at least three types of stabilizing interactions, including conventional CO2(-)⋯H2O, H2O⋯H2O charge assisted hydrogen bonds (HBs), and less common H2O⋯H-C and H2O⋯π interactions. Biological water molecules, those in direct contact with Ibu(-), prefer to cluster around the carboxylate oxygen atoms via cyclic or bridged charge assisted hydrogen bonds. Many of those interactions are strongly affected by the formal carboxylate charge, resulting in enhanced HBs with increased strengths and degree of covalency. We found striking similarities between this case and the microsolvation of dymethylphosphate, which lead us to hypothesize that since microsolvation of phosphatidylcholine depends mainly on the formal charge of its ionic PO2(-) group in the polar head, then microsolvation of anionic ibuprofen and interactions of water molecules with eukaryotic cell membranes are governed by the same types of physical interactions.