MOLECULAR DYNAMICS STUDIES OF SIMPLE MODEL FLUIDS AND WATER CONFINED IN CARBON NANOTUBE

Molecular Dynamics (MD) simulation is one of the most important computational techniques with broad applications in physics, chemistry, chemical engineering, materials design and biological science. Traditional computational chemistry refers to quantum calculations based on solving Schrodinger equations. Later developed Density Functional Theory (DFT) based on solving Kohn-Sham equations became the more popular ab initio calculation technique which could deal with ~1000 atoms by explicitly considering electron interactions. In contrast, MD simulation based on solving classical mechanics equations of motion is a totally different technique in the field of computational chemistry. Electron interactions were implicitly included in the empirical atom-based potential functions and the system size to be investigated can be extended to ~106 atoms. The thermodynamic properties of model fluids are mainly determined by macroscopic quantities, like temperature, pressure, density. The quantum effects on thermodynamic properties like melting point, surface tension are not dominant. In this work, we mainly investigated the melting point, surface tension (liquid-vapor and liquid-solid) of model fluids including Lennard-Jones model, Stockmayer model and a couple of water models (TIP4P/Ew, TIP5P/Ew) by means of MD simulation. In addition, some new structures of water confined in carbon nanotube were discovered and transport behaviors of water and ions through nano-channels were also revealed.

Computation of (3)J(HH) coupling constants with a combination of density functional theory and semiempirical calculations. Application to complex molecules

This work evaluates the efficiency of economic levels of theory for the prediction of (3)J(HH) spin-spin coupling constants, to be used when robust electronic structure methods are prohibitive. To that purpose, DFT methods like mPW1PW91. B3LYP and PBEPBE were used to obtain coupling constants for a test set whose coupling constants are well known. Satisfactory results were obtained in most of cases, with the mPW1PW91/6-31G(d,p)//B3LYP/6-31G(d,p) leading the set. In a second step. B3LYP was replaced by the semiempirical methods PM6 and RM1 in the geometry optimizations. Coupling constants calculated with these latter structures were at least as good as the ones obtained by pure DFT methods. This is a promising result, because some of the main objectives of computational chemistry - low computational cost and time, allied to high performance and precision - were attained together. (C) 2012 Elsevier B.V. All rights reserved.

delta-Lactam derivative from thermophilic soil fungus exhibits in vitro anti-allergic activity

From cultures of thermophilic soil fungus Humicola grisea var thermoidea, a delta-lactam derivative (3-(2-(4-hydroxyphenyl)-2-oxoethyl)-5,6-dihydropyridin-2( 1H)-one) that displayed anti-allergic activity was isolated, which was predicted by in silico computational chemistry approaches. The in vitro anti-allergic activity was investigated by beta-hexosaminidase release assay in rat basophilic leukaemia RBL-2H3 cells. The delta-lactam derivative exhibited similar anti-allergic activity (IC50 = 18.7 +/- 6.7 mu M) in comparison with ketotifen fumarate (IC50 = 15.0 +/- 1.3 mu M) and stronger anti-allergic activity than azelastine (IC50 = 32.0 mu M). Also, the MTT cytotoxicity assay with RBL-2H3 cells showed that delta-lactam does not display cytotoxicity at concentrations lower than 50 mu M. This study suggests that the delta-lactam derivative has the potential to be used as a lead compound in the development of anti-allergic drugs for clinical use in humans.

Systematic UPLC-ESI-MS/MS Study on the Occurrence of Staurosporine and Derivatives in Associated Marine Microorganisms from Eudistoma vannamei

Former bioactivity-guided analysis of the marine invertebrate Eudistoma vannamei led to the isolation of staurosporine derivatives, which revealed strong cytotoxic activity against several human cancer cell lines. The occurrence of such alkaloids in E. vannamei may be correlated to the presence of associated biota, such as Streptomyces bacteria. In agreement to this hypothesis, marine microorganisms associated with E. vannamei were recovered and cultured, leading to a total of 84 isolated bacterial strains. Gas phase fragmentation reactions of staurosporine and derivatives were systematically studied and the analyzed results further supported by computational chemistry studies. The resulting fragment patterns were used to search for the presence of different derivatives in extracts of isolated microorganisms, thereby using LC-MS/MS analysis in MRM mode. These results evidenced that one isolated Streptomyces sp. was able to generate staurosporine, while none of the hydroxy-7-oxo derivatives were detected. Finally, significant cytotoxic activity against human cancer lines was observed for one of the extracts.

Systematic UPLC-ESI-MS/MS study on the occurrence of staurosporine and derivatives in associated marine microorganisms from Eudistoma vannamei

Former bioactivity-guided analysis of the marine invertebrate Eudistoma vannamei led to the isolation of staurosporine derivatives, which revealed strong cytotoxic activity against several human cancer cell lines. The occurrence of such alkaloids in E. vannamei may be correlated to the presence of associated biota, such as Streptomyces bacteria. In agreement to this hypothesis, marine microorganisms associated with E. vannamei were recovered and cultured, leading to a total of 84 isolated bacterial strains. Gas phase fragmentation reactions of staurosporine and derivatives were systematically studied and the analyzed results further supported by computational chemistry studies. The resulting fragment patterns were used to search for the presence of different derivatives in extracts of isolated microorganisms, thereby using LC-MS/MS analysis in MRM mode. These results evidenced that one isolated Streptomyces sp. was able to generate staurosporine, while none of the hydroxy-7-oxo derivatives were detected. Finally, significant cytotoxic activity against human cancer lines was observed for one of the extracts.

Large-scale coupled-cluster calculations

Coupled-cluster theory provides one of the most successful concepts in electronic-structure theory. This work covers the parallelization of coupled-cluster energies, gradients, and second derivatives and its application to selected large-scale chemical problems, beside the more practical aspects such as the publication and support of the quantum-chemistry package ACES II MAB and the design and development of a computational environment optimized for coupled-cluster calculations. The main objective of this thesis was to extend the range of applicability of coupled-cluster models to larger molecular systems and their properties and therefore to bring large-scale coupled-cluster calculations into day-to-day routine of computational chemistry. A straightforward strategy for the parallelization of CCSD and CCSD(T) energies, gradients, and second derivatives has been outlined and implemented for closed-shell and open-shell references. Starting from the highly efficient serial implementation of the ACES II MAB computer code an adaptation for affordable workstation clusters has been obtained by parallelizing the most time-consuming steps of the algorithms. Benchmark calculations for systems with up to 1300 basis functions and the presented applications show that the resulting algorithm for energies, gradients and second derivatives at the CCSD and CCSD(T) level of theory exhibits good scaling with the number of processors and substantially extends the range of applicability. Within the framework of the ’High accuracy Extrapolated Ab initio Thermochemistry’ (HEAT) protocols effects of increased basis-set size and higher excitations in the coupled- cluster expansion were investigated. The HEAT scheme was generalized for molecules containing second-row atoms in the case of vinyl chloride. This allowed the different experimental reported values to be discriminated. In the case of the benzene molecule it was shown that even for molecules of this size chemical accuracy can be achieved. Near-quantitative agreement with experiment (about 2 ppm deviation) for the prediction of fluorine-19 nuclear magnetic shielding constants can be achieved by employing the CCSD(T) model together with large basis sets at accurate equilibrium geometries if vibrational averaging and temperature corrections via second-order vibrational perturbation theory are considered. Applying a very similar level of theory for the calculation of the carbon-13 NMR chemical shifts of benzene resulted in quantitative agreement with experimental gas-phase data. The NMR chemical shift study for the bridgehead 1-adamantyl cation at the CCSD(T) level resolved earlier discrepancies of lower-level theoretical treatment. The equilibrium structure of diacetylene has been determined based on the combination of experimental rotational constants of thirteen isotopic species and zero-point vibrational corrections calculated at various quantum-chemical levels. These empirical equilibrium structures agree to within 0.1 pm irrespective of the theoretical level employed. High-level quantum-chemical calculations on the hyperfine structure parameters of the cyanopolyynes were found to be in excellent agreement with experiment. Finally, the theoretically most accurate determination of the molecular equilibrium structure of ferrocene to date is presented.

Reactivity of activated electrophiles and nucleophiles: labile intermediates and properties of the reaction products

The main topic of my Ph.D. thesis is the study of nucleophilic and electrophilic aromatic substitution reaction, in particular from a mechanistic point of view. The research was mainly focused on the reactivity of superactivated aromatic systems. In spite of their high reactivity (hence the high reaction’s rate), we were able to identify and in some case to isolate -complexes until now only hypothesized. For example, interesting results comes from the study of the protonation of the supernucleophiles tris(dialkylamino)benzenes. However, the best result obtained in this field was the isolation and structural characterization of the first stables zwitterionic Wheland-Meisenheimer complexes by using 2,4-dipyrrolidine-1,3-thiazole as supernucleophile and 4,6-dinitrobenzofuroxan or 4,6-dinitrotetrazolepyridine as superelectrophile. These reactions were also studied by means of computational chemistry, which allowed us to better investigate on the energetic and properties of the reactions and reactants studied. We also discovered, in some case fortuitously, some relevant properties and application of the compounds we synthesized, such as fluorescence in solid state and nanoparticles, or textile dyeing. We decided to investigate all these findings also by collaborating with other research groups. During a period in the “Laboratoire de Structure et Réactivité des Systèmes Moléculaires Complexes-SRSMC, Université de Lorraine et CNRS, France, I carried out computational studies on new iron complexes for the use as dyes in Dye Sensitized Solar Cells (DSSC). Furthermore, thanks to this new expertise, I was involved in a collaboration for the study of the ligands’ interaction in biological systems. A collaboration with University of Urbino allowed us to investigate on the reactivity of 1,2-diaza-1,3-dienes toward nucleophiles such as amino and phosphine derivatives, which led to the synthesis of new products some of which are 6 or 7 member heterocycles containing both phosphorus and nitrogen atoms.

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.

Multiscale constitutive modeling of polymer materials

Materials are inherently multi-scale in nature consisting of distinct characteristics at various length scales from atoms to bulk material. There are no widely accepted predictive multi-scale modeling techniques that span from atomic level to bulk relating the effects of the structure at the nanometer (10-9 meter) on macro-scale properties. Traditional engineering deals with treating matter as continuous with no internal structure. In contrast to engineers, physicists have dealt with matter in its discrete structure at small length scales to understand fundamental behavior of materials. Multiscale modeling is of great scientific and technical importance as it can aid in designing novel materials that will enable us to tailor properties specific to an application like multi-functional materials. Polymer nanocomposite materials have the potential to provide significant increases in mechanical properties relative to current polymers used for structural applications. The nanoscale reinforcements have the potential to increase the effective interface between the reinforcement and the matrix by orders of magnitude for a given reinforcement volume fraction as relative to traditional micro- or macro-scale reinforcements. To facilitate the development of polymer nanocomposite materials, constitutive relationships must be established that predict the bulk mechanical properties of the materials as a function of the molecular structure. A computational hierarchical multiscale modeling technique is developed to study the bulk-level constitutive behavior of polymeric materials as a function of its molecular chemistry. Various parameters and modeling techniques from computational chemistry to continuum mechanics are utilized for the current modeling method. The cause and effect relationship of the parameters are studied to establish an efficient modeling framework. The proposed methodology is applied to three different polymers and validated using experimental data available in literature.

Cyclobutadiene automerization and rotation of ethylene: Energetics of the barriers by using Spin-polarized wave functions

Spin-projected spin polarized Møller–Plesset and spin polarized coupled clusters calculations have been made to estimate the cyclobutadiene automerization, the ethylene torsion barriers in their ground state, and the gap between the singlet and triplet states of ethylene. The results have been obtained optimizing the geometries at MP4 and/or CCSD levels, by an extensive Gaussian basis set. A comparative analysis with more complex calculations, up to MP5 and CCSDTQP, together with others from the literature, have also been made, showing the efficacy of using spin-polarized wave functions as a reference wave function for Møller–Plesset and coupled clusters calculations, in such problems.

Towards Accurate and Efficient Description of Excited States

Thesis (Ph.D.)--University of Washington, 2016-06

On the microscopic mechanism of carbon gasification: A theoretical study

In this paper we report the results of ab initio calculations on the energetics and kinetics of oxygen-driven carbon gasification reactions using a small model cluster, with full characterisation of the stationary points on the reaction paths. We show that previously unconsidered pathways present significantly reduced barriers to reaction and must be considered as alternative viable paths. At least two electronic spin states of the model cluster must be considered for a complete description. (C) 2004 Elsevier Ltd. All rights reserved.

Wurst: a protein threading server with a structural scoring function, sequence profiles and optimized substitution matrices

Wurst is a protein threading program with an emphasis on high quality sequence to structure alignments (http://www.zbh.uni-hamburg.de/wurst). Submitted sequences are aligned to each of about 3000 templates with a conventional dynamic programming algorithm, but using a score function with sophisticated structure and sequence terms. The structure terms are a log-odds probability of sequence to structure fragment compatibility, obtained from a Bayesian classification procedure. A simplex optimization was used to optimize the sequence-based terms for the goal of alignment and model quality and to balance the sequence and structural contributions against each other. Both sequence and structural terms operate with sequence profiles.