Hydration of OCS with One to Four Water Molecules in Atmospheric and Laboratory Conditions

Carbonyl sulfide is the most abundant sulfur gas in the atmosphere. We have used MP2 and CCSD(T) theory to study the structures and thermochemistries of carbonyl sulfide interacting with one to four water molecules. We have completed an extensive search for clusters of OCS(H2O)n, where n = 1−4. We located three dimers, two trimers, five tetramers, and four pentamers with the MP2/aug-cc-pVDZ method. In each of the complexes with two or more waters, OCS preferentially interacts with low-energy water clusters. Our results match current theoretical and experimental literature, showing correlation with available geometries and frequencies for the OCS(H2O) species. The CCSD(T)/aug-cc-pVTZ thermochemical values combined with the average amount of OCS and the saturated concentration of H2O in the troposphere, lead to the prediction of 106 OCS(H2O) clusters·cm−3 and 102 OCS(H2O)2 clusters·cm−3 at 298 K. We predict the structures of OCS(H2O)n, n = 1−4 that should predominate in a low-temperature molecular beam and identify specific infrared vibrations that can be used to identify these different clusters.

Accurate Predictions of Water Cluster Formation, (H2O)n=2−10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Structure and thermodynamics of H3O+(H2O)8 clusters: A combined molecular dynamics and quantum mechanics approach

We have studied the structure and stability of (H3O+)(H2O)8 clusters using a combination of molecular dynamics sampling and high-level ab initio calculations. 20 distinct oxygen frameworks are found within 2 kcal/mol of the electronic or standard Gibbs free energy minimum. The impact of quantum zero-point vibrational corrections on the relative stability of these isomers is quite significant. The box-like isomers are favored in terms of electronic energy, but with the inclusion of zero-point vibrational corrections and entropic effects tree-like isomers are favored at higher temperatures. Under conditions from 0 to 298.15 K, the global minimum is predicted to be a tree-like structure with one dangling singly coordinated water molecule. Above 298.15 K, higher entropy tree-like isomers with two or more singly coordinated water molecules are favored. These assignments are generally consistent with experimental IR spectra of (H3O+)(H2O)8 obtained at 150 K.

Broadband Fourier transform rotational spectroscopy for structure determination: The water heptamer

Over the recent years chirped-pulse, Fourier-transform microwave (CP-FTMW) spectrometers have chan- ged the scope of rotational spectroscopy. The broad frequency and large dynamic range make possible structural determinations in molecular systems of increasingly larger size from measurements of heavy atom (13C, 15N, 18O) isotopes recorded in natural abundance in the same spectrum as that of the parent isotopic species. The design of a broadband spectrometer operating in the 2–8 GHz frequency range with further improvements in sensitivity is presented. The current CP-FTMW spectrometer performance is benchmarked in the analyses of the rotational spectrum of the water heptamer, (H2O)7, in both 2– 8 GHz and 6–18 GHz frequency ranges. Two isomers of the water heptamer have been observed in a pulsed supersonic molecular expansion. High level ab initio structural searches were performed to pro- vide plausible low-energy candidates which were directly compared with accurate structures provided from broadband rotational spectra. The full substitution structure of the most stable species has been obtained through the analysis of all possible singly-substituted isotopologues (H218O and HDO), and a least-squares rm(1) geometry of the oxygen framework determined from 16 different isotopic species compares with the calculated O–O equilibrium distances at the 0.01 Å level.

Fluorite and Mixed-Metal Kagome-Related Topologies in Metal-Organic Framework Compounds: Synthesis, Structure, and Properties

Two new three-dimensional metal-organic frameworks (MOFs) [Mn-2(mu(3)-OH)(H2O)(2)(BTC)]-2 H2O, I, and [NaMn(BTC)], II (BTC=1,2,4-benzenetricarboxylate = trimellitate) were synthesized and their structures determined by single-crystal X-ray diffraction (XRD). In I, the Mn-4 cluster, [Mn-4(mu(5)-OH)(2)(H2O)(4)O-12], is connected with eight trimellitate anions and each trimellitate anion connects to four different Mn-4 clusters, resulting in a fluorite-like structure. In II, the Mn2O8 dimer is connected with two Na+ ions through carboxylate oxygen to form mixed-metal distorted Kagome-related two-dimensional -M-O-M- layers, which are pillared by the trimellitate anions forming the three-dimensional structure. The extra-framework water molecules in I are reversibly adsorbed and are also corroborated by powder XRD studies. The formation of octameric water clusters involving free and coordinated water molecules appears to be new. Interesting magnetic behavior has been observed for both compounds. Electron spin resonance (ESR) studies indicate a broadening of the signal below the ordering temperature and appear to support the findings of the magnetic studies.

Photodissociation and reaction dynamics of chlorine-containing species important in stratospheric ozone chemistry

Chlorine oxide species have received considerable attention in recent years due to their central role in the balance of stratospheric ozone. Many questions pertaining to the behavior of such species still remain unanswered and plague the ability of researchers to develop accurate chemical models of the stratosphere. Presented in this thesis are three experiments that study various properties of some specific chlorine oxide species.

In the first chapter, the reaction between ClONO_2 and protonated water clusters is investigated to elucidate a possible reaction mechanism for the heterogeneous reaction of chlorine nitrate on ice. The ionic products were various forms of protonated nitric acid, NO_2 +(H_20)_m, m = 0, 1, 2. These products are analogous to products previously reported in the literature for the neutral reaction occurring on ice surfaces. Our results support the hypothesis that the heterogeneous reaction is acid-catalyzed.

In the second chapter, the photochemistry of ClONO_2 was investigated at two wavelengths, 193 and 248 nm, using the technique of photofragmentation translational spectroscopy. At both wavelengths, the predominant dissociation pathways were Cl + NO_3 and ClO + NO_2. Channel assignments were confirmed by momentum matching the counterfragments from each channel. A one-dimensional stratospheric model using the new 248 nm branching ratio determined how our results would affect the predicted Cl_x and NO_x partitioning in the stratosphere.

Chapter three explores the photodissociation dynamics of Cl_2O at 193, 248 and 308 nm. At 193 nm, we found evidence for the concerted reaction channel, Cl_2 + O. The ClO + Cl channel was also accessed, however, the majority of the ClO fragments were formed with sufficient internal energies for spontaneous secondary dissociation to occur. At 248 and 308 nm, we only observed only the ClO + Cl channel. . Some of the ClO formed at 248 nm was formed internally hot and spontaneously dissociated. Bimodal translational energy distributions of the ClO and Cl products indicate two pathways leading to the same product exist.

Appendix A, B and C discuss the details of data analysis techniques used in Chapters 1 and 2. The development of a molecular beam source of ClO dimer is presented in Appendix D.

Machine learning for predictive condensed-phase simulation

We show how machine learning techniques based on Bayesian inference can be used to reach new levels of realism in the computer simulation of molecular materials, focusing here on water. We train our machine-learning algorithm using accurate, correlated quantum chemistry, and predict energies and forces in molecular aggregates ranging from clusters to solid and liquid phases. The widely used electronic-structure methods based on density-functional theory (DFT) give poor accuracy for molecular materials like water, and we show how our techniques can be used to generate systematically improvable corrections to DFT. The resulting corrected DFT scheme gives remarkably accurate predictions for the relative energies of small water clusters and of different ice structures, and greatly improves the description of the structure and dynamics of liquid water.

Syntheses, structures, and characterizations of four new silver(I) sulfonate coordination polymers with neutral ligands

In this paper, four novel silver(I) sulfonate coordination polymers containing neutral ligands, namely, [Ag(2)Ll (biim)(2)]center dot 2H(2)O (1). AgL2(biim) (2), [Ag(HL3)(Pic)(2)]center dot H2O (3), and [Ag-3(L3)(HL3)(4,4'-bipy)(3)(H2O)(2)]center dot 4H(2)O (4), have been synthesized [L1 = 3-carboxy-4-hydroxybenzenesulfonate, L2 = p-aminobenzenesulfonate, H(2)L3 = p-hydroxybenzenesulfonic acid, biim = 1,1'-(1.4-butanediyl)-bis(imidazole), Pic = beta-picoline, 4,4'-bipy = 4,4'-bipyridine]. For compounds 1 and 2, Ag(I) cations are bridged by biim ligands to form a one-dimensional (1D) "zigzag" chain, and L1 and L2 sulfonate ligands are not coordinated to the silver cation. Compound 3 has a dimeric structure in which two silver cations are bridged by two HL3 ligands. For compound 4, L3 ligand coordinates to a silver cation as a monodentate ligand, and Ag(l) cations are bridged by 4,4'-bipy ligands to form a ID chain. Compound 1 contains water dimers, while compound 4 contains water trimers. Compounds 1-3 display room-temperature photoluminescence.

Simulations of Biased Agonists in the β2 Adrenergic Receptor with Accelerated Molecular Dynamics

The biased agonism of the G protein-coupled receptors (GPCRs), where in addition to a traditional G protein-signalling pathway a GPCR promotes intracellular signals though ß-arrestin, is a novel paradigm in pharmacology. Biochemical and biophysical studies have suggested that a GPCR forms a distinct ensemble of conformations signalling through the G protein and ß-arrestin. Here we report on the dynamics of the ß2 adrenergic receptor bound to the ß-arrestin and G protein biased agonists and the empty receptor to further characterize the receptor conformational changes caused by biased agonists. We use conventional and accelerated molecular dynamics (aMD) simulations to explore the conformational transitions of the GPCR from the active state to the inactive state. We found that aMD simulations enable monitoring the transition within the nanosecond timescale while capturing the known microscopic characteristics of the inactive states, such as the ionic lock, the inward position of F6.44, and water clusters. Distinct conformational states are shown to be stabilized by each biased agonist. In particular, in simulations of the receptor with the ß-arrestin biased agonist, N-cyclopentylbutanepherine we observe a different pattern of motions in helix 7 when compared to simulations with the G protein biased agonist, Salbutamol that involves perturbations of the network of interactions within the NPxxY motif. Understanding the network of interactions induced by biased ligands and the subsequent receptor conformational shifts will lead to development of more efficient drugs. © 2013 American Chemical Society

Organic synthesis in the interstellar medium by low-energy carbon irradiation

We present a first principles molecular dynamics (FPMD) study of the interaction of low-energy neutral carbon projectiles with amorphous solid water clusters at 30 K. Reactions involving the carbon atom at an initial energy of 11 and 1.7 eV with 30-molecule clusters have been investigated. Simulations indicate that the formation of hydroxymethylene, an intermediate in formaldehyde production, dominates at the higher energy. The reaction proceeds by fragmenting a water molecule, binding the carbon to the OH radical, and saturating the C valence with a hydrogen atom that can arise from the originally dissociated water molecule, or through a chain of proton transfer events. We identified several possible pathways for the formation of HCOH. When the initial collision occurs at the periphery of the cluster, we observe the formation of CO and the evaporation of water molecules. At the lower energy water fragmentation is not favorable, thus leading to the formation of weakly bound carbon-water complexes. © 2013 American Chemical Society.

Structure and dynamics of aqueous 2-propanol: A THz-TDS, NMR and neutron diffraction study

Aqueous liquid mixtures, in particular, those involving amphiphilic species, play an important role in many physical, chemical and biological processes. Of particular interest are alcohol/water mixtures; however, the structural dynamics of such systems are still not fully understood. Herein, a combination of terahertz time-domain spectroscopy (THz-TDS) and NMR relaxation time analysis has been applied to investigate 2-propanol/water mixtures across the entire composition range; while neutron diffraction studies have been carried out at two specific concentrations. Excellent agreement is seen between the techniques with a maximum in both the relative absorption coefficient and the activation energy to molecular motion occurring at ∼90 mol% H₂O. Furthermore, this is the same value at which well-established excess thermodynamic functions exhibit a maximum/minimum. Additionally, both neutron diffraction and THz-TDS have been used to provide estimates of the size of the hydration shell around 2-propanol in solution. Both methods determine that between 4 and 5 H₂O molecules per 2-propanol are found in the 2-propanol/water clusters at 90 mol% H₂O. Based on the acquired data, a description of the structure of 2-propanol/water across the composition range is presented.

Interacções intermoleculares fracas em materiais funcionais avançados: um estudo espectroscópico e computacional

A utilização combinada de espectroscopia vibracional e de cálculos envolvendo a teoria do funcional de densidade (DFT) possibilita o estudo de ligações de hidrogénio em fase condensada, assim como a análise da estrutura molecular dos sistemas em estudo. Por um lado, a espectroscopia vibracional permite a detecção de associações moleculares, enquanto os métodos computacionais auxiliam na obtenção de informação referente aos mecanismos de associação, nomeadamente no que diz respeito à possível estrutura de dímeros e compostos de inclusão em ciclodextrinas e às energias de interacção e de inclusão. O estudo que originou a presente dissertação pretende contribuir para o reforço da aplicação de estudos espectroscópicos e computacionais na elucidação de diversos fenómenos químicos, com especial destaque para o papel desempenhado por interacções intermoleculares fracas na estrutura e propriedades de materiais moleculares. No âmbito desta tese foram investigados os seguintes tópicos: polimorfismo e pseudopolimorfismo em sólidos farmacêuticos, transições de fase em misturas binárias de ácidos gordos, inclusão em ciclodextrinas, interacção de compostos farmacêuticos com superfícies metálicas e formação de agregados de água em materiais híbridos orgânicos-inorgânicos. Os sistemas foram analisados utilizando a espectroscopia vibracional – particularmente a espectroscopia de difusão de Raman – como técnica fundamental. Para uma melhor caracterização de processos envolvendo transições de fase, foram efectuados estudos com variação de temperatura, variação de humidade relativa e substituição isotópica. O estudo da interacção com superfícies metálicas foi realizado por espectroscopia de Raman intensificada à superfície. Dada a complexidade dos sistemas em estudo, a informação obtida por espectroscopia vibracional foi complementada por resultados de cálculos mecânico-quânticos. Em particular, os cálculos DFT foram utilizados para a optimização de geometrias e previsão de frequências vibracionais de moléculas e associações moleculares, permitindo assim a análise e interpretação de espectros vibracionais e a caracterização da estrutura de materiais.

Cobalt complexes bearing scorpionate ligands: Synthesis, characterization, cytotoxicity and DNA cleavage

A number of novel, water-stable redox-active cobalt complexes of the C-functionalized tripodal ligands tris(pyrazolyl)methane XC(pz)(3) (X = HOCH2, CH2OCH2Py or CH2OSO2Me) are reported along with their effects on DNA. The compounds were isolated as air-stable solids and fully characterized by IR and FIR spectroscopies, ESI-MS(+/-), cyclic voltammetry, controlled potential electrolysis, elemental analysis and, in a number of cases, also by single-crystal X-ray diffraction. They showed moderate cytotoxicity in vitro towards HCT116 colorectal carcinoma and HepG2 hepatocellular carcinoma human cancer cell lines. This viability loss is correlated with an increase of tumour cell lines apoptosis. Reactivity studies with biomolecules, such as reducing agents, H2O2, plasmid DNA and UV-visible titrations were also performed to provide tentative insights into the mode of action of the complexes. Incubation of Co(II) complexes with pDNA induced double strand breaks, without requiring the presence of any activator. This pDNA cleavage appears to be mediated by O-centred radical species.

Molecular model for hydrated biological tissues

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Mejillonesite, a new acid sodium, magnesium phosphate mineral, from Mejillones, Antofagasta, Chile

Mejillonesite, ideally NaMg(2)(PO(3)OH)(PO(4))(OH)center dot H(5)O(2), is a new mineral approved by the CNMNC (IMA 2010-068). It occurs as isolated crystal aggregates in thin zones in fine-grained opal-zeolite aggregate on the north slope of Cerro Mejillones, Antofagasta, Chile. Closely associated minerals are bobierrite, opal, clinoptilolite-Na, clinoptilolite-K, and gypsum. Mejillonesite forms orthorhombic, prismatic, and elongated thick tabular crystals up to 6 mm long, usually intergrown in radiating aggregates. The dominant form is pinacoid {100}. Prisms {hk0}, {h0l}, and {0kl} are also observed. The crystals are colorless, their streak is white, and the luster is vitreous. The mineral is transparent. It is non-fluorescent under ultraviolet light. Mohs' hardness is 4, tenacity is brittle. Cleavage is perfect on {100}, good on {010} and {001}, and fracture is stepped. The measured density is 2.36(1) g/cm(3); the calculated density is 2.367 g/cm(3). Mejillonesite is biaxial (-), alpha= 1.507(2), beta= 1.531(2), gamma= 1.531(2), 2V(meas) = 15(10)degrees, 2V(calc) = 0 degrees (589 nm). Orientation is X= a, Z= elongation direction. The mineral is non-pleochroic. Dispersion is r> v, medium. The IR spectrum contains characteristic bands of the Zundel cation (H(5)O(2)(+), or H(+)center dot 2H(2)O) and the groups P-OH and OH(-). The chemical composition is (by EDS, H(2)O by the Alimarin method, wt%): Na(2)O 9.19, MgO 26.82, P(2)O(5) 46.87, H(2)O 19, total 101.88. The empirical formula, based on 11 oxygen atoms, is Na(0.93)Mg(2.08)(PO(3)OH)(1.00) (PO(4)) (OH)(0.86) .0.95H(5)O(2) The strongest eight X-ray powder-diffraction lines [d in angstrom(I)(hkl)] are: 8.095(100)(200), 6.846(9) (210), 6.470(8)(111), 3.317(5)(302), 2.959(5)(132), 2.706(12)(113), 2.157(19)(333), and 2.153(9) (622). The crystal structure was solved on a single crystal (R = 0.055) and gave the following data: orthorhombic, Pbca, a = 16.295(1), b = 13.009(2), c = 8.434(1) angstrom, V= 1787.9(4) angstrom(3), Z = 8. The crystal structure of mejillonesite is based on a sheet (parallel to the b-c plane) formed by two types of MgO(6) octahedra, isolated tetrahedra PO(4) and PO(3)OH whose apical vertices have different orientation with respect to the sheet. The sheets are connected by interlayer, 5-coordinated sodium ions, proton hydration complexes, and hydroxyl groups. The structure of mejillonesite is related to that of angarfite, NaFe(5)(3+)(PO(4))(4)(OH)(4).4H(2)O and bakhchisaraitsevite, Na(2)Mg(5)(PO(4))(4)center dot 7H(2)O.