Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Au@Hg nanoalloy formation through direct amalgamation: Structural, spectroscopic, and computational evidence for slow nanoscale diffusion

Dynamic core-shell nanoparticles have received increasing attention in recent years. This paper presents a detailed study of Au-Hg nanoalloys, whose composing elements show a large difference in cohesive energy. A simple method to prepare Au@Hg particles with precise control over the composition up to 15 atom% mercury is introduced, based on reacting a citrate stabilized gold sol with elemental mercury. Transmission electron microscopy shows an increase of particle size with increasing mercury content and, together with X-ray powder diffraction, points towards the presence of a core-shell structure with a gold core surrounded by an Au-Hg solid solution layer. The amalgamation process is described by pseudo-zero-order reaction kinetics, which indicates slow dissolution of mercury in water as the rate determining step, followed by fast scavenging by nanoparticles in solution. Once adsorbed at the surface, slow diffusion of Hg into the particle lattice occurs, to a depth of ca. 3 nm, independent of Hg concentration. Discrete dipole approximation calculations relate the UV-vis spectra to the microscopic details of the nanoalloy structure. Segregation energies and metal distribution in the nanoalloys were modeled by density functional theory calculations. The results indicate slow metal interdiffusion at the nanoscale, which has important implications for synthetic methods aimed at core-shell particles.

Conformationally Controlled Electron Delocalization in n-Type Rods: Synthesis, Structure, and Optical, Electrochemical, and Spectroelectrochemical Properties of Dicyanocyclophanes

A series of dicyanobiphenyl-cyclophanes 1-6 with various pi-backbone conformations and characteristic n-type semiconductor properties is presented. Their synthesis, optical, structural, electrochemical, spectroelectrochemical, and packing properties are investigated. The X-ray crystal structures of all n-type rods allow the systematic correlation of structural features with physical properties. In addition, the results are supported by quantum mechanical calculations based on density functional theory. A two-step reduction process is observed for all n-type rods, in which the first step is reversible. The potential gap between the reduction processes depends linearly on the cos(2) value of the torsion angle phi between the pi-systems. Similarly, optical absorption spectroscopy shows that the vertical excitation energy of the conjugation band correlates with the cos(2) value of the torsion angle phi. These correlations demonstrate that the fixed intramolecular torsion angle phi is the dominant factor determining the extent of electron delocalization in these model compounds, and that the angle phi measured in the solid-state structure is a good proxy for the molecular conformation in solution. Spectroelectrochemical investigations demonstrate that conformational rigidity is maintained even in the radical anion form. In particular, the absorption bands corresponding to the SOMO-LUMO+i transitions are shifted bathochromically, whereas the absorption bands corresponding to the HOMO-SOMO transition are shifted hypsochromically with increasing torsion angle phi.

Intermolecular Clamping by Hydrogen Bonds: 2-Pyridone center dot NH3

A combined spectroscopic and ab initio theoretical study of the doubly hydrogen-bonded complex of 2-pyridone (2PY) with NH3 has been performed. The S-1 <- S-0 spectrum extends up to approximate to 1200 cm(-1) above the 0(0)(0) band, close to twice the range observed for 2PY. The S-1 state nonradiative decay for vibrations above approximate to 300 cm(-1) in the NH3 complex is dramatically slowed down relative to bare 2PY. Also, the Delta v=2,4,... overtone bands of the v(1)' and v(2)' out-of-plane vibrations that dominate the low-energy spectral region of 2PY are much weaker or missing for 2PY center dot NH3, which implies that the bridging (2PY)NH center dot center dot center dot NH3 and H2NH center dot center dot center dot O=C H-bonds clamp the 2PY at a planar geometry in the S-1 state. The mass-resolved UV vibronic spectra of jet-cooled 2PY center dot NH3 and its H/D mixed isotopomers are measured using two-color resonant two-photon ionization spectroscopy. The S-0 and S-1 equilibrium structures and normal-mode frequencies are calculated by density functional (B3LYP) and correlated ab initio methods (MP2 and approximate second-order coupled-cluster, CC2). The S-1 <- S-0 vibronic assignments are based on configuration interaction singles (CIS) and CC2 calculations. A doubly H-bonded bridged structure of C-S symmetry is predicted, in agreement with that of Held and Pratt [J. Am. Chem. Soc. 1993, 115, 9718]. While the B3LYP and MP2 calculated rotational constants are in very good agreement with experiment, the calculated H2NH center dot center dot center dot O=C H-bond distance is approximate to 0.7 angstrom shorter than that derived by Held and Pratt. On the other hand, this underlines their observation that ammonia can act as a strong H-bond donor when built into an H-bonded bridge. The CC2 calculations predict the H2NH center dot center dot center dot O distance to increase by 0.2 angstrom upon S-1 <- S-0 electronic excitation, while the (2PY)NH center dot center dot center dot NH3 H-bond remains nearly unchanged. Thus, the expansion of the doubly H-bonded bridge in the excited state is asymmetric and almost wholly due to the weakening of the interaction of ammonia with the keto acceptor group.

A tetrathiafulvalene-functionalized naphthalene diimide: synthesis, electrochemical and photophysical properties

A tetrathiafulvalene donor has been attached to the naphthalene diimide core via a rigid bridge affording a new planar molecular dyad. Its electronic properties have been studied experimentally by the combination of electrochemistry and UV-vis-NIR spectroscopy. Various electronic excited charge-transfer states are generated in different oxidation states, leading to almost full absorption in the visible to near-IR region with high extinction coefficients. The observed electronic properties are explained on the basis of density-functional-theory. In particular, the oxidized radical species show a strong tendency to undergo aggregation, in which the long-distance attractive interactions overcome the electrostatic repulsions. (C) 2011 Elsevier Ltd. All rights reserved.

A tetrathiafulvalene-functionalized schiff base macrocycle: synthesis, electrochemical, and photophysical properties

Highly selective formation of 2+2 macrocycle 1 from 2,5-bis(3-formyl-2-hydroxyphenyl)-1,3,4-oxadiazole and a diamine-functionalized tetrathiafulvalene derivative is reported. Its electronic properties have been studied experimentally by the combination of electrochemistry and UV-vis-NIR spectroscopy. Particularly, its largely extended pi-conjugation renders this novel macrocycle simultaneously a good multielectron donor and a strong chromophore, which is rationalized on the basis of density functional theory. (C) 2011 Elsevier Ltd. All rights reserved.

Excited state properties of formamide in water solution: An ab initio study

We present ab initio quantum calculation of the optical properties of formamide in vapor phase and in water solution. We employ time dependent density functional theory for the isolated molecule and many-body perturbation theory methods for the system in solution. An average over several molecular dynamics snapshots is performed to take into account the disorder of the liquid. We find that the excited stateproperties of the gas-phase formamide are strongly modified by the presence of the water solvent: the geometry of the molecule is distorted and the electronic and optical properties are severely modified. The important interaction among the formamide and the water molecules forces us to use fully quantum methods for the calculation of the excited stateproperties of this system. The excitonic wave function is localized both on the solute and on part of the solvent.

Structural Analysis of alpha-Fetoprotein (AFP)-like Peptides with Anti-Breast-Cancer Properties

The abundance of alpha-fetoprotein (AFP), a natural protein produced by the fetal yolk sac during pregnancy, correlates with lower incidence of estrogen receptor positive (ER+) breast cancer. The pharmacophore region of AFP has been narrowed down to a four amino acid (AA) region in the third domain of the 591 AA peptide. Our computational study focuses on a 4-mer segment consisting of the amino acids threonine-proline-valine-asparagine (TPVN). We have run replica exchange molecular dynamics (REMD) simulations and used 120 configurational snapshots from the total trajectory as starting configurations for quantum chemical calculations. We optimized structures using semiempirical (PM3, PM6, PM6-D2, PM6-H2, PM6-DH+, PM6-DH2) and density functional methods (TPSS, PBE0, M06-2X). By comparing the accuracy of these methods against RI-MP2 benchmarks, we devised a protocol for calculating the lowest energy conformers of these peptides accurately and efficiently. This protocol screens out high-energy conformers using lower levels of theory and outlines a general method for predicting small peptide structures.

A Computational Analysis of the Collision-Induced Dissociation Mechanisms of the Dipeptide Glycine-Serine Under Mass Spectrometry

Collision-induced dissociation (CID) of peptides using tandem mass spectrometry (MS) has been used to determine the identity of peptides and other large biological molecules. Mass spectrometry (MS) is a useful tool for determining the identity of molecules based on their interaction with electromagnetic fields. If coupled with another method like infrared (IR) vibrational spectroscopy, MS can provide structural information, but in its own right, MS can only provide the mass-to-charge (m/z) ratio of the fragments produced, which may not be enough information to determine the mechanism of the collision-induced dissociation (CID) of the molecule. In this case, theoretical calculations provide a useful companion for MS data and yield clues about the energetics of the dissociation. In this study, negative ion electrospray tandem MS was used to study the CID of the deprotonated dipeptide glycine-serine (Gly-Ser). Though negative ion MS is not as popular a choice as positive ion MS, studies by Bowie et al. show that it yields unique clues about molecular structure which complement positive ion spectroscopy, such as characteristic fragmentations like the loss of formaldehyde from the serine residue.2 The increase in the collision energy in the mass spectrometer alters the flexibility of the dipeptide backbone, enabling isomerizations (reactions not resulting in a fragment loss) and dissociations to take place. The mechanism of the CID of Gly-Ser was studied using two computational methods, B3LYP/6-311+G* and M06-2X/6-311++G**. The main pathway for molecular dissociation was analyzed in 5 conformers in an attempt to verify the initial mechanism proposed by Dr. James Swan after examination of the MS data. The results suggest that the loss of formaldehyde from serine, which Bowie et al. indicates is a characteristic of the presence of serine in a protein residue, is an endothermic reaction that is made possible by the conversion of the translational energy of the ion into internal energy as the ion collides with the inert collision gas. It has also been determined that the M06-2X functional¿s improved description of medium and long-range correlation makes it more effective than the B3LYP functional at finding elusive transition states. M06-2X also more accurately predicts the energy of those transition states than does B3LYP. A second CID mechanism, which passes through intermediates with the same m/z ratio as the main pathway for molecular dissociation, but different structures, including a diketopiperazine intermediate, was also studied. This pathway for molecular dissociation was analyzed with 3 conformers and the M06-2X functional, due to its previously determined effectiveness. The results suggest that the latter pathway, which meets the same intermediate masses as the first mechanism, is lower in overall energy and therefore a more likely pathway of dissociation than the first mechanism.

Effect of external electric field on hydrogen adsorption over activated carbon separated by dielectric materials

Energy crisis and worldwide environmental problem make hydrogen a prospective energy carrier. However, storage and transportation of hydrogen in large quantities at small volume is currently not practical. Lots of materials and devices have been developed for storage hydrogen, but to today none is able to meet the DOE targets. Activated carbon has been found to be a good hydrogen adsorbent due to its high surface area. However, the weak van der Waals force between hydrogen and the adsorbent has limited the adsorption capacity. Previous studies have found that enhanced adsorption can be obtained with applied electric field. Stronger interaction between the polarized hydrogen and the charged sorbents under high voltage is considered as the reason. This study was initiated to investigate if the adsorption can be further enhanced when the activated carbon particles are separated with a dielectric coating. Dielectric TiO2 nanoparticles were first utilized. Hydrogen adsorption measurements on the TiO2-coated carbon materials, with or without an external electric field, were made. The results showed that the adsorption capacity enhancement increased with the increasing amount of TiO2 nanoparticles with an applied electric field. Since the hydrogen adsorption capacity on TiO2 particles is very low and there is no hydrogen adsorption enhancement on TiO2 particles alone when electric field is applied, the effect of dielectric coating is demonstrated. Another set of experiments investigated the behavior of hydrogen adsorption over TiO2-coated activated carbon under various electric potentials. The results revealed that the hydrogen adsorption first increased and then decreased with the increase of electric field. The improved storage was due to a stronger interaction between charged carbon surface and polarized hydrogen molecule caused by field induced polarization of TiO2 coating. When the electric field was sufficient to cause considerable ionization of hydrogen, the decrease of hydrogen adsorption occurred. The current leak detected at 3000 V was a sign of ionization of hydrogen. Experiments were also carried out to examine the hydrogen adsorption performances over activated carbon separated by other dielectric materials, MgO, ZnO and BaTiO3, respectively. For the samples partitioned with MgO and ZnO, the measurements with and without an electric field indicated negligible differences. Electric field enhanced adsorption has been observed on the activated carbon separated with BaTiO3, a material with unusually high dielectric constant. Corresponding computational calculations using Density Functional Theory have been performed on hydrogen interaction with charged TiO2 molecule as well as TiO2 molecule, coronene and TiO2-doped coronene in the presence of an electric field. The simulated results were consistent with the observations from experiments, further confirming the proposed hypotheses.

Electronic properties of zig-zag carbon nanotubes : a first-principles study

Carbon nanotube (CNT) is a one dimensional (1-D) nanostructured material, which has been the focal point of research over the past decade for intriguing applications ranging from nanoelectronics to chemical and biological sensors. Using a first-principles gradient corrected density functional approach, we present a comprehensive study of the geometry and energy band gap in zig-zag semi-conducting (n,0) carbon nanotubes (CNT) to resolve some of the conflicting findings. Our calculations confirm that the single wall (n,0) CNTs fall into two distinct classes depending upon n mod 3 equal to 1 (smaller band gaps) or 2 (larger gaps). The effect of longitudinal strain on the band gap further confirms the existence of two distinct classes: for n mod 3 = 1 or 2, changing Eg by ~ ±110 meV for 1% strain in each case. We also present our findings for the origin of metallicity in multiwall CNTs.

First-principles studies of boron nanostructures

Boron is an 'electron deficient' element which has a rather fascinating chemical versatility. In the solid state, the elemental boron has neither a pure covalent nor a pure metallic character. As a result, its vast structural dimensionally and peculiar bonding features hold a unique place among other elements in the periodic table. In order to understand and properly describe these unusual bonding features, a detailed and systematic theoretical study is needed. In this work, I will show that some of the qualitative features of boron nanostructures, including clusters, sheets and nanotubes can easily be extracted from the results of first principles calculations based on density functional theory. Specifically, the size-dependent evolution of topological structures and bonding characteristics of boron clusters, Bn will be discussed. Based on the scenario observed in the boron clusters, the unique properties of boron sheets and boron nanotubes will be described. Moreover, the ballistic electron transport in single-walled carbon nanotubes will be considered. It is expected that the theoretical results obtained in the present thesis will initiate further studies on boron nanostructures, which will be helpful in understanding, designing and realizing boron-based nanoscale devices.

Development of a high performance parallel computing platform and its use in the study of nanostructures : clusters, sheets and tubes .

Small clusters of gallium oxide, technologically important high temperature ceramic, together with interaction of nucleic acid bases with graphene and small-diameter carbon nanotube are focus of first principles calculations in this work. A high performance parallel computing platform is also developed to perform these calculations at Michigan Tech. First principles calculations are based on density functional theory employing either local density or gradient-corrected approximation together with plane wave and gaussian basis sets. The bulk Ga2O3 is known to be a very good candidate for fabricating electronic devices that operate at high temperatures. To explore the properties of Ga2O3 at nonoscale, we have performed a systematic theoretical study on the small polyatomic gallium oxide clusters. The calculated results find that all lowest energy isomers of GamOn clusters are dominated by the Ga-O bonds over the metal-metal or the oxygen-oxygen bonds. Analysis of atomic charges suggest the clusters to be highly ionic similar to the case of bulk Ga2O3. In the study of sequential oxidation of these slusters starting from Ga2O, it is found that the most stable isomers display up to four different backbones of constituent atoms. Furthermore, the predicted configuration of the ground state of Ga2O is recently confirmed by the experimental result of Neumark's group. Guided by the results of calculations the study of gallium oxide clusters, performance related challenge of computational simulations, of producing high performance computers/platforms, has been addressed. Several engineering aspects were thoroughly studied during the design, development and implementation of the high performance parallel computing platform, rama, at Michigan Tech. In an attempt to stay true to the principles of Beowulf revolutioni, the rama cluster was extensively customized to make it easy to understand, and use - for administrators as well as end-users. Following the results of benchmark calculations and to keep up with the complexity of systems under study, rama has been expanded to a total of sixty four processors. Interest in the non-covalent intereaction of DNA with carbon nanotubes has steadily increased during past several years. This hybrid system, at the junction of the biological regime and the nanomaterials world, possesses features which make it very attractive for a wide range of applicatioins. Using the in-house computational power available, we have studied details of the interaction between nucleic acid bases with graphene sheet as well as high-curvature small-diameter carbon nanotube. The calculated trend in the binding energies strongly suggests that the polarizability of the base molecules determines the interaction strength of the nucleic acid bases with graphene. When comparing the results obtained here for physisorption on the small diameter nanotube considered with those from the study on graphene, it is observed that the interaction strength of nucleic acid bases is smaller for the tube. Thus, these results show that the effect of introducing curvature is to reduce the binding energy. The binding energies for the two extreme cases of negligible curvature (i.e. flat graphene sheet) and of very high curvature (i.e. small diameter nanotube) may be considered as upper and lower bounds. This finding represents an important step towards a better understanding of experimentally observed sequence-dependent interaction of DNA with Carbon nanotubes.

Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.