Applying many-body expansion aiming to make ab in initio molecular dynamics more efficient

In this present work we present a methodology that aims to apply the many-body expansion to decrease the computational cost of ab initio molecular dynamics, keeping acceptable accuracy on the results. We implemented this methodology in a program which we called ManBo. In the many-body expansion approach, we partitioned the total energy E of the system in contributions of one body, two bodies, three bodies, etc., until the contribution of the Nth body [1-3]: E = E1 + E2 + E3 + …EN. The E1 term is the sum of the internal energy of the molecules; the term E2 is the energy due to interaction between all pairs of molecules; E3 is the energy due to interaction between all trios of molecules; and so on. In Manbo we chose to truncate the expansion in the contribution of two or three bodies, both for the calculation of the energy and for the calculation of the atomic forces. In order to partially include the many-body interactions neglected when we truncate the expansion, we can include an electrostatic embedding in the electronic structure calculations, instead of considering the monomers, pairs and trios as isolated molecules in space. In simulations we made we chose to simulate water molecules, and use the Gaussian 09 as external program to calculate the atomic forces and energy of the system, as well as reference program for analyzing the accuracy of the results obtained with the ManBo. The results show that the use of the many-body expansion seems to be an interesting approach for reducing the still prohibitive computational cost of ab initio molecular dynamics. The errors introduced on atomic forces in applying such methodology are very small. The inclusion of an embedding electrostatic seems to be a good solution for improving the results with only a small increase in simulation time. As we increase the level of calculation, the simulation time of ManBo tends to largely decrease in relation to a conventional BOMD simulation of Gaussian, due to better scalability of the methodology presented. References [1] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 3, 46 (2007). [2] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 4, 1 (2008). [3] R. Rivelino, P. Chaudhuri and S. Canuto; J. Chem. Phys., 118, 10593 (2003).

ESTIMATING REACTION CONSTANTS BY AB INITIO MOLECULAR MODELING: A STUDY ON THE OXIDATION OF PHENOL TO CATECHOL AND HYDROQUINONE IN ADVANCED OXIDATION PROCESSES

Molecular modeling is growing as a research tool in Chemical Engineering studies, as can be seen by a simple research on the latest publications in the field. Molecular investigations retrieve information on properties often accessible only by expensive and time-consuming experimental techniques, such as those involved in the study of radical-based chain reactions. In this work, different quantum chemical techniques were used to study phenol oxidation by hydroxyl radicals in Advanced Oxidation Processes used for wastewater treatment. The results obtained by applying a DFT-based model showed good agreement with experimental values available, as well as qualitative insights into the mechanism of the overall reaction chain. Solvation models were also tried, but were found to be limited for this reaction system within the considered theoretical level without further parameterization.

Estimating reaction constants by ab initio molecular modeling: a study on the oxidation of phenol to catechol and hydroquinone in advanced oxidation processes

Molecular modeling is growing as a research tool in Chemical Engineering studies, as can be seen by a simple research on the latest publications in the field. Molecular investigations retrieve information on properties often accessible only by expensive and time-consuming experimental techniques, such as those involved in the study of radical-based chain reactions. In this work, different quantum chemical techniques were used to study phenol oxidation by hydroxyl radicals in Advanced Oxidation Processes used for wastewater treatment. The results obtained by applying a DFT-based model showed good agreement with experimental values available, as well as qualitative insights into the mechanism of the overall reaction chain. Solvation models were also tried, but were found to be limited for this reaction system within the considered theoretical level without further parameterization.

Accurate ab initio potential energy surfaces for the (3)A '' and (3)A ' electronic states of the O(P-3) plus HBr system

In this work, we report the construction of potential energy surfaces for the (3)A '' and (3)A' states of the system O(P-3) + HBr. These surfaces are based on extensive ab initio calculations employing the MRCI+Q/CBS+SO level of theory. The complete basis set energies were estimated from extrapolation of MRCI+Q/aug-cc-VnZ(-PP) (n = Q, 5) results and corrections due to spin-orbit effects obtained at the CASSCF/aug-cc-pVTZ(-PP) level of theory. These energies, calculated over a region of the configuration space relevant to the study of the reaction O(P-3) + HBr -> OH + Br, were used to generate functions based on the many-body expansion. The three-body potentials were interpolated using the reproducing kernel Hilbert space method. The resulting surface for the (3)A '' electronic state contains van der Waals minima on the entrance and exit channels and a transition state 6.55 kcal/mol higher than the reactants. This barrier height was then scaled to reproduce the value of 5.01 kcal/mol, which was estimated from coupled cluster benchmark calculations performed to include high-order and core-valence correlation, as well as scalar relativistic effects. The (3)A' surface was also scaled, based on the fact that in the collinear saddle point geometry these two electronic states are degenerate. The vibrationally adiabatic barrier heights are 3.44 kcal/mol for the (3)A '' and 4.16 kcal/mol for the (3)A' state. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4705428]

New insights on reaction pathway selectivity promoted by crown ether phase-transfer catalysis: Model ab initio calculations of nucleophilic fluorination

Crown ethers have the ability of solubilizing inorganic salts in apolar solvents and to promote chemical reactions by phase-transfer catalysis. However, details on how crown ethers catalyze ionic S(N)2 reactions and control selectivity are not well understood. In this work, we have used high level theoretical calculations to shed light on the details of phase-transfer catalysis mechanism of KF reaction with alkyl halides promoted by 18-crown-6. A complete analysis of the of the model reaction between KF(18-crown-6) and ethyl bromide reveals that the calculations can accurately predict the product ratio and the overall kinetics. Our results point out the importance of the K* ion and of the crown ether ring in determining product selectivity. While the K* ion favors the S(N)2 over the E2 anti pathway, the crown ether ring favors the S(N)2 over E2 syn route. The combination effects lead to a predicted 94% for the S(N)2 pathway in excellent agreement with the experimental value of 92%. A detailed analysis of the overall mechanism of the reaction under phase-transfer conditions also reveals that the KBr product generated in the nucleophilic fluorination acts as an inhibitor of the 18-crown-6 catalyst and it is responsible for the observed slow reaction rate. (C) 2012 Elsevier B.V. All rights reserved.

Ab initio studies of pristine and oxidized Cu3Au(100) and (111) surfaces

We review the previous literature and our recent work on first-principles studies of Cu3Au(100) and (111) surfaces, with focus on the segregation of atomic species to the surface at pristine conditions and in the presence of oxygen. In particular, the combined use of experimental and theoretical tools to achieve chemical identification at an atomic level of the surface species is emphasized and discussed.

A novel set of ss-N-biaryl ether sulfonamide hydroxamates as potential MMPs inhibitors: Molecular dynamics simulations and molecular properties evaluation

Matrix metalloproteinases (MMPs) constitute a family of zinc-dependent proteases involved in the extracellular matrix degradation. MMP-2 and MMP9 are overexpressed in several human cancer types, including melanoma, thus the development of new compounds to inhibit MMPs' activity is desirable. Molecular dynamic simulation and molecular properties calculations were performed on a set of novel beta-N-biaryl ether sulfonamide-based hydroxamates, reported as MMP-2 and MMP-9 inhibitors, for providing data to develop an exploratory analysis. Thermodynamic, electronic, and steric descriptors have significantly discriminated highly active from moderately and less active inhibitors of MMP-2 whereas apparent partition coefficient at pH 1.5 was also significant for the MMP-9 data set. Compound 47 was considered an outlier in all analysis, indicating the presence of a bulky substituent group in R3 is crucial to this set of inhibitors for the establishment of molecular interactions with the S1 subsite of both enzymes, but there is a limit. (C) 2012 Wiley Periodicals, Inc.

The Structural and Electronic Properties of Tin Oxide Nanowires: An Ab Initio Investigation

We performed an ab initio investigation on the properties of rutile tin oxide (SnOx) nanowires. We computed the wire properties determining the equilibrium geometries, binding energies, and electronic band structures for several wire dimensions and surface facet configurations. The results allowed us to establish scaling laws for the structural properties, in terms of the nanowire perimeters. The results also showed that the surface states control most of the electronic properties of the nanowires. Oxygen incorporation in the nanowire surfaces passivated the surface-related electronic states, and the resulting quantum properties and scaling laws were fully consistent with electrons confined inside the nanowire. Additionally, oxygen incorporation in the wire surfaces generated an unbalanced concentration of spin up and down electrons, leading to magnetic states for the nanowires.

Ab-initio calculations for a realistic sensor: A study of CO sensors based on nitrogen-rich carbon nanotubes

The use of nanoscale low-dimensional systems could boost the sensitivity of gas sensors. In this work we simulate a nanoscopic sensor based on carbon nanotubes with a large number of binding sites using ab initio density functional electronic structure calculations coupled to the Non-Equilibrium Green's Function formalism. We present a recipe where the adsorption process is studied followed by conductance calculations of a single defect system and of more realistic disordered system considering different coverages of molecules as one would expect experimentally. We found that the sensitivity of the disordered system is enhanced by a factor of 5 when compared to the single defect one. Finally, our results from the atomistic electronic transport are used as input to a simple model that connects them to experimental parameters such as temperature and partial gas pressure, providing a procedure for simulating a realistic nanoscopic gas sensor. Using this methodology we show that nitrogen-rich carbon nanotubes could work at room temperature with extremely high sensitivity. Copyright 2012 Author(s). This article is distributed under a Creative Commons Attribution 3.0 Unported License. [http://dx.doi.org/10.1063/1.4739280]

Photodynamic Efficiency of Cationic meso-Porphyrins at Lipid Bilayers: Insights from Molecular Dynamics Simulations

Porphyrin derivatives have applications as photoactive drugs in photodynamic therapy. However, little is known about their interactions with phospholipid membranes at the molecular level. We employed molecular dynamics simulations to model the binding between a series of cationic meso-(N-methyl-4-pyridinium)phenylporphyrins and anionic phosphatidylglycerol lipid bilayers. This was done in the presence of molecular oxygen within the membrane. The ability of various porphyrins to cause photodamage was quantified in terms of their immersion depth and degree of exposition to a higher oxygen concentration inside the membrane. Simulations showed that the photodynamic efficiency could be improved as the number of hydrophobic phenyl substituents attached to the porphyrinic ring increased. In the specific case of porphyrins containing two hydrophobic and two charged substituents, the cis isomer was significantly more efficient than the trans. These results correlate well with previous experimental observations. They highlight the importance of both the total charge and amphiphilicity of the photosensitizer for its performance in photodynamic therapy.

Influence of Molecular Dynamics on the Dielectric Properties of Poly(9,9-di-n-octylfluorene-alt-benzothiadiazole) -Based Devices

This paper uses Nuclear Magnetic Resonance (NMR) and Differential Scanning Calorimetry (DSC) techniques to study the molecular relaxations and phase transitions in poly(9,9-di-n-octylfluorene-alt-benzothiadiazole) (F8BT), which has been extensively studied as the active thin film in organic devices. Besides the identification of the glass transition, beta relaxation and crystal-to-crystal phase transition, we correlate such phenomena with dielectric and transport mechanisms in diodes with F8BT as the active layer. The beta relaxation has been assigned to a transition at about 210 K measured by H-1 and C-13 solid state NMR, and can be attributed to local motions in the side chains. The glass transition has been detected by DSC and H-1 NMR. Dielectric spectroscopy (DS) carried out at low frequencies on diodes made from F8BT show two peaks which are coincident with the above transitions. This allowed us to correlate the electrical changes in the film with the onset of specific molecular motions. In addition, DS indicates a third peak related with a crystal-to-crystal phase transition. Finally, these transitions were correlated with changes in the carrier mobility recorded in thin films and published recently.

Effect of the thermostat in the molecular dynamics simulation on the folding of the model protein chignolin

Molecular dynamics simulations of the model protein chignolin with explicit solvent were carried out, in order to analyze the influence of the Berendsen thermostat on the evolution and folding of the peptide. The dependence of the peptide behavior on temperature was tested with the commonly employed thermostat scheme consisting of one thermostat for the protein and another for the solvent. The thermostat coupling time of the protein was increased to infinity, when the protein is not in direct contact with the thermal bath, a situation known as minimally invasive thermostat. In agreement with other works, it was observed that only in the last situation the instantaneous temperature of the model protein obeys a canonical distribution. As for the folding studies, it was shown that, in the applications of the commonly utilized thermostat schemes, the systems are trapped in local minima regions from which it has difficulty escaping. With the minimally invasive thermostat the time that the protein needs to fold was reduced by two to three times. These results show that the obstacles to the evolution of the extended peptide to the folded structure can be overcome when the temperature of the peptide is not directly controlled.

Molecular Dynamics Investigations of PRODAN in a DLPC Bilayer

Molecular dynamics computer simulations have been performed to identify preferred positions of the fluorescent probe PRODAN in a fully hydrated DLPC bilayer in the fluid phase. In addition to the intramolecular charge-transfer first vertical excited state, we considered different charge distributions for the electronic ground state of the PRODAN molecule by distinct atomic charge models corresponding to the probe molecule in vacuum as well as polarized in a weak and a strong dielectric solvent (cyclohexane and water). Independent on the charge distribution model of PRODAN, we observed a preferential orientation of this molecule in the bilayer with the dimethylamino group pointing toward the membrane's center and the carbonyl oxygen toward the membrane's interface. However, changing the charge distribution model of PRODAN, independent of its initial position in the equilibrated DLPC membrane, we observed different preferential positions. For the ground state representation without polarization and the in-cyclohexane polarization, the probe maintains its position close to the membrane's center. Considering the in-water polarization model, the probe approaches more of the polar headgroup region of the bilayer, with a strong structural correlation with the choline group, exposing its oxygen atom to water molecules. PRODAN's representation of the first vertical excited state with the in-water polarization also approaches the polar region of the membrane with the oxygen atom exposed to the bilayer's hydration shell. However, this model presents a stronger structural correlation with the phosphate groups than the ground state. Therefore, we conclude that the orientation of the PRODAN molecule inside the DLPC membrane is well-defined, but its position is very sensitive to the effect of the medium polarization included here by different models for the atomic charge distribution of the probe.

The T-2 phase in the Nb-Si-B system studied by ab initio calculations and synchrotron X-ray diffraction

The solid solution based on Nb5Si3 (Cr5B3 structure type, D8(l), tl32, 14/mcm, No140, a=6.5767 angstrom, c=11.8967 angstrom) in the Nb-Si-B system was studied from the structural and thermodynamic point of view both experimentally and by ab initio calculations. Rietveld refinement of powder X-ray synchrotron data allowed to determine the boron to silicon substitution mechanism and the structural parameters. Ab initio calculations of different ordered compounds and selected disordered alloys allowed to obtain in addition to the enthalpy of formation of the solution, substitution mechanism and structural parameters which are in excellent agreement with the experimental data. The stability of the phase is discussed. (C) 2012 Elsevier Inc. All rights reserved.

Selectivity and Mechanisms Driven by Reaction Dynamics: The Case of the Gas-Phase OH-+CH3ONO2 Reaction

Well-established statistical approaches such as transition-state theory based on high-level calculated potential energy profiles are unable to account for the selectivity observed in the gas-phase OH- + CH3ONO2 reaction. This reaction can undergo bimolecular nucleophilic displacement at either the carbon center (S(N)2@C) or the nitrogen center (S(N)2@N) as well as a proton abstraction followed by dissociation (E(CO)2) pathway. Direct dynamics simulations yield an S(N)2:E(CO)2 product ratio in close agreement with experiment and show that the lack of reactivity at the nitrogen atom is due to the highly negative electrostatic potential generated by the oxygen atoms in the ONO2 group that scatters the incoming OH-. In addition to these dynamical effects, the nonstatistical behavior of these reactions is attributed to the absence of equilibrated reactant complexes and to the large number of recrossings, which might be present in several ion-molecule gas-phase reactions.