Microsolvation of glycine: a DFT study

PBE1PBE/6-311+G(d,p) computations exploring the microsolvation of neutral and zwitterionic glycine are reported. A broad configuration search was performed to identify the lowest energy clusters of glycine with one to seven water molecules. The structures of the clusters are analyzed on the basis of the hydrogen-bonding network established between the water molecules and between water and glycine. Neutral glycine is favored when associated with zero to six water molecules, but with seven water molecules the two structures are isoenergetic.     

DFT Calculation of Glycine with Methanols Clusters

Density functional theory （DFT） calculations are reported for the structures of neutral and zwitterionic glycine-（CHaOH）n where n=1-6. Initial geometries of the clusters of neutral and zwitterionic glycine with 1-6 methanol molecules are fully optimized at B3LYP/6-31＋G^＊ level of theory. The lowest energy configurations are located and their hydrogen bond structures are analyzed. Theoretical prediction reveals that the methanols prefer to locate near the carboxylic acid group for the small clusters （n_〈3） with the neutral form while the configurations with the methanols bridging the acid and the amino group are favorite in the zwitterionic form clusters. When the number of the methanol molecules in the clusters reaches five and six, the two forms tend to be isoenergetic.     

Density functional theory study of 1:1 glycine-water complexes in the gas phase and in solution

We present a systematic study of 1:1 glycine-water complexes involving all possible glycine conformers. The complex geometries are fully optimized for the first time both in the gas phase and in solution using three DFT methods (B3LYP, PBE1PBE, X3LYP) and the MP2 method. We calculate the G3 energies and use them as the reference data to gauge hydrogen bond strength in the gas phase. The solvent effects are treated via the integral equation formalism-polarizable continuum model (IEF-PCM). Altogether, we locate fifty-two unique nonionized (N) structures and six zwitterionic (Z) structures in the gas phase, and fifty-five N structures and thirteen Z structures in solution. Both correlation and solvation are shown to be important in geometry determination. We found that in the gas phase, a water molecule binds more strongly to the carboxylic acid group of glycine than to its amine group, whereas in solution phase the reverse is true. The most stable Z structure is isoenergetic with the most stable N structure.     

Performance of six functionals (LDA, PBE, PBESOL, B3LYP, PBE0, and WC1LYP) in the simulation of vibrational and dielectric properties of crystalline compounds. The case of forsterite Mg2SiO4

The performance of six different density functionals (LDA, PBE, PBESOL, B3LYP, PBE0, and WC1LYP) in describing the infrared spectrum of forsterite, a crystalline periodic system with orthorhombic unit cell (28 atoms in the primitive cell, Pbmn space group), is investigated by using the periodic ab initio CRYSTAL09 code and an all‐electron Gaussian‐type basis set. The transverse optical (TO) branches of the 35 IR active modes are evaluated at the equilibrium geometry together with the oscillator strengths and the high‐frequency dielectric tensor ?_∞. These quantities are essential to compute the dielectric function ?(ν), and then the reflectance spectrum R(ν), which is compared with experiment. It turns out that hybrid functionals perform better than LDA and GGA, in general; that B3LYP overperforms WC1LYP and, in turn, PBE0; that PBESOL is better than PBE; that LDA is the worst performing functional among the six under study. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

A density functional theory study of the structure and vibrational spectra of beta-carotene, capsanthin, and capsorubin

A theoretical study of the structure and the vibrational spectra of the beta-carotene molecule and its derivatives capsanthin and capsorubin is carried out. We first investigate systematically the theoretical method which provides the best results for beta-carotene by performing ab initio calculations at the HF/6-31G(d), SVWN/6-31G(d), PBE0/6-31G(d), BLYP/6-31G(d), B3LYP/6-31G(d), B3LYP/6-31G(d,p), B3LYP/6-311G(d), and B3LYP/6-311G(d,p) levels and by using previous theoretical results available in the literature obtained at the AM1 and BPW91/6-31G(d) levels. The influence of both the level of calculation and the size of the basis set used in the geometry optimization and in the determination of the IR and Raman spectra of this molecule is thus analyzed. It is confirmed that the hybrid functional B3LYP with the basis 6-31G(d) is the method that gives the best results as a whole. By use of this level of calculation, we next optimize the molecular geometries of related molecules of capsanthin and capsorubin, which to the best of our knowledge have only been studied at the semiempirical AM1 level. In addition we calculate the IR and Raman spectra of these molecules at the B3LYP/6-31G(d) level of theory. The results obtained for capsanthin show on the one hand that the double bond of the beta-ionone ring is outside the polyene chain plane, due to the repulsion between the hydrogen atoms of the ring methyl groups and the hydrogen atoms of the polyene chain, and on the other hand that the carbonyl double bond in the other headgroup is very close to planarity with the polyene chain, since in this case such a repulsion does not exist. For the molecule of capsorubin the two carbonyl groups also take the same coplanar orientation relative to the polyene chain. The IR and Raman spectra theoretically computed for these two molecules are finally compared with their experimental spectra and the vibrational normal modes of the main signals are interpreted.     

Theoretical 13C chemical shift, 14N, and 2H quadrupole coupling- constant studies of hydrogen bonding in L-alanylglycine dipeptide

(13)C chemical shieldings and (14)N and (2)H electric field gradient (EFG) tensors of L-alanylglycine (L-alagly) dipeptide were calculated at RHF/6-31 + + G** and B3LYP/6-31 + + G** levels of theory respectively. For these calculations a crystal structure of this dipeptide obtained from X-ray crystallography was used. Atomic coordinates of different clusters containing several L-alagly molecules were used as input files for calculations. These clusters consist of central and surrounding L-alagly molecules, the latter forming short, strong, hydrogen bonds with the central molecule. Since the calculations did not converge for these clusters, the surrounding L-alagly molecules were replaced by glycine molecules. In order to improve the accuracy of calculated chemical shifts and nuclear quadrupole coupling constants (NQCCs), different geometry-optimization strategies were applied for hydrogen nuclei. Agreement between calculated and experimental data confirms that our optimized coordinates for hydrogen nuclei are more accurate than those obtained by X-ray diffraction.     

Computational study of static first hyperpolarizability of donor-acceptor substituted (E)-benzaldehyde phenylhydrazone

Various hybrid functionals (B3LYP, B97-2, PBE0, BMK, BH&HLYP, CAM-B3LYP, and LC-ωPBE) implemented in density functional theory were applied to give estimate of static first hyperpolarizabilty (β(0)) of (E)-benzaldehyde phenylhydrazone designated as (E)-BPH. Against those of MP2 computations as a function of the underlying density functional, good agreement was obtained with the BH&HLYP and CAM-B3LYP functionals. The LC-ωPBE functional and the B3LYP, PBE0, B97-2, and BMK functionals underestimated and overestimated β(0), respectively. The basis set effect on the calculated β(0) was also investigated. It turned out that the 6-311+G(2d,p) basis set provided excellent converged value of β(0). On the basis of the calculated results, we investigated the substituent effect on β(0) of donor-acceptor (D-A) substituted (E)-BPH systematically by using the BH&HLYP and CAM-B3LYP computations with the 6-311+G(2d,p) basis set. We proposed a Zwitterion structure to explain the calculated trend in the substituent effect and the enhanced hyperpolarizability of type II compounds (A-(E)-BPH-D) than type I compounds (D-(E)-BPH-A). Natural bonding orbital analysis carried out at BH&HLYP/6-311+G(2d,p)//B3LYP/6-31G(2df,p) level of theory substantiated the claim.     

Structures and stability of medium-sized silicon clusters. III. Reexamination of motif transition in growth pattern from Si15 to Si20

It has been established from experiments that stable medium-sized ionic clusters Si15-Si20 are prolate in shape. Density-functional theories (DFTs) also predict that nearly all low-lying neutral clusters in this size range are prolate in shape. Moreover, most of them are built onto two generic structural motifs, either the tricapped-trigonal-prism (TTP) Si9 motif or the six/six Si6Si6 (sixfold-puckered hexagonal ring Si6 plus six-atom tetragonal bipyramid Si6) motif. However, it appears that the exact location of the TTP-to-six/six motif transition is dependent on the functional (e.g., PBE or BLYP) used in the DFT calculations. Here, we present total-energy calculations for two series of clusters (one series containing six/six motif and the other containing the TTP motif) in the size range of Si16-Si20. The calculations were based on all-electron DFT methods with a medium [6-311G (2d)] and a large (cc-pVTZ) basis sets, as well as coupled-cluster single and double substitutions (including triple excitations) [CCSD(T)] method with a modest (cc-pVDZ) basis set. In the DFT calculations, two popular hybrid density functionals, the B3LYP and PBE1PBE, were selected. It is found that the B3LYP total-energy calculations slightly favor the six/six motif, whereas the PBE1PBE calculations slightly favor the TTP motif. The CCSD(T) total-energy calculations, however, show that isomers based on the six/six motif are energetically slightly favorable in the size range of Si16-Si20. Hence, the TTP-to-six/six motif transition is more likely to occur at Si16.     

Remarks on GIAO-DFT predictions of 13C chemical shifts

Predicting (13)C chemical shifts by GIAO-DFT calculations appears to be more accurate than frequently expected provided that: (a) the comparison between experimental and theoretical data is performed using the linear regression method, (b) a sufficiently high level of theory [e.g. B3LYP/6-311 + + G(2d,p)//B3LYP/6-311 + + G(2d,p) or PBE1PBE/6-311 + G(2df,p)//B3LYP/6-311 + + G(2d,p)] is used, (c) the experimental data originate from the measurements performed in one solvent whose influence is taken into account at the molecular geometry optimization step and, first of all, during the shielding calculation, (d) the experimental data are free of heavy atom effects or such effects are appropriately treated in calculations, and finally (e) the conformational compositions of the investigated objects are known.     

Comparison of model chemistry and density functional theory thermochemical predictions with experiment for formation of ionic clusters of the ammonium cation complexed with water and ammonia; atmospheric implications

The G2, G3, CBS-QB3, and CBS-APNO model chemistry methods and the B3LYP, B3P86, mPW1PW, and PBE1PBE density functional theory (DFT) methods have been used to calculate deltaH(o) and deltaG(o) values for ionic clusters of the ammonium ion complexed with water and ammonia. Results for the clusters NH4(+) (NH3)n and NH4(+) (H2O)n, where n = 1-4, are reported in this paper and compared against experimental values. Agreement with the experimental values for deltaH(o) and deltaG(o) for formation of NH4(+) (NH3)n clusters is excellent. Comparison between experiment and theory for formation of the NH4(+) (H2O)n clusters is quite good considering the uncertainty in the experimental values. The four DFT methods yield excellent agreement with experiment and the model chemistry methods when the aug-cc-pVTZ basis set is used for energetic calculations and the 6-31G* basis set is used for geometries and frequencies. On the basis of these results, we predict that all ions in the lower troposphere will be saturated with at least one complete first hydration shell of water molecules.     

DFT和热力学研究氢键协同效应及对关联1H NMR的影响

用DFT方法在B3LYP／6—311 G(d,p)水平下研究了甲醇线性和环状分子簇．对于不同大小的分子簇之间定义了协同因子,计算得到的协同因子可以验证氢键的强协同效应,环状分子簇之间的协同效应远远大于线性分子簇,做为理论验证和比较,热力学模型分别采用含氢键缔合的格子流体状态方程(LFHB),以及含氢键协同效应的LFHB,关联醇一惰性体系的^1H核磁共振化学位移．考虑协同效应的关联结果优于原始的LFHB,比较量子化学计算的和热力学模型中采用的协同因子,认为甲醇和乙醇在溶液中更可能大部分以线性缔合形式存在。     

用密度泛函方法研究质子化苯基丙酮-水团簇

The structure and growth trend of the protonated acetophenone-water clusters have been investigated using the DFT-B3LYP method combined with the standard 6-31+G(d,p) basis set. In order to obtain more accurate single-point energy the B3LYP/6-311++G(3df,2p) method was adapted. The results show that the formation of H+C8H8O-H2O is a barrierless reaction process and the equilibrium distance between the proton and the O atom in C8H8O molecule is 1.015 A. For H+C8H8O-(H2O)n(n=1,2,3) clusters, the proton lies between the acetophenone molecule C8H8O and the water molecule H2O. The distance between the proton and the O atom of the C8H8O molecule increased from n=1 to n=3; C8H8O-H+-H2O can be regarded as an solvation shell. For H+C8H8O (H2O)n (n=4,5,6,7,8) clusters, the proton lies between the two H2O molecules forming a H5O2+ structure, C8H8O-H5O2+ is an important structure, which the other H2O molecules will attack from different sides.     

Theoretical study of hydrogen bonded clusters of water and fulminic acid

Ab initio and density functional calculations are used to analyze the interaction between a molecule of fulminic acid with 1, 2, 3, and 4 molecules of water along with a 2:2 complex at B3LYP/6-31+G(d,p) and MP2/6-311++G(d,p) computational levels. Cooperative effect (CE) in terms of stabilization energy of clusters are calculated and discussed as well. CE is increased with increasing cluster size of studied clusters. Red shifts of H–C stretching frequency for complexes involving HCNO as H-donor are predicted.     

DFT study of the electronic and charge transfer properties of perfluoroarene–thiophene oligomers

The ground state structures of 5,5″-diperfluorophenyl-2,2′:5′,2″:5″,2‴-quaterthiophene (1), 5,5′-bis{1-[4-(thien-2-yl)perfluorophenyl]}-2,2′-dithiophene (2), 4,4′-bis[5-(2,2′-dithiophenyl)]-perfluorobiphenyl (3), 5,5″-diperfluorophenyl-2,2′:5′,2″-tertthiophene (4), 5,5′-diperfluorophenyl-2,2′-dihiophene (5), and 5,5-diperfluorophenylthiophene (6) have been optimized at the B3LYP/6-31G(d), B3LYP/6-31G(d,p), PBE0/6-31G(d), and PBE0/6-31G(d,p) level of theories. The B3LYP/6-31+G(d) and PBE0/6-31+G(d) level of theories have been applied to investigate the absorption spectra. The PBE0 functional is good to predict the C–S bond lengths while the C–F bond lengths are good envisaged with B3LYP functional. The increment of thiophene rings between two perfluoroarene rings leads to red shift in absorption spectra. The electron affinities are energetically destabilized while energetic stabilization of the radical-cation increases by decreasing the thiophene rings from four to one. The perfluoroarene rings leads to enhance the electron injection.     

Absorption and Fluorescence Emission Spectra of Molecules Containing Azo and/or Oxadiazole Units

The ground state structures of a series of organic molecules containing azo and/or oxadiazole units were obtained by means of density functional theory B3LYP/6-31G(d) method. The first singlet excited state structures were optimized by virtue of singlet-excitation configuration interaction CIS/6-31G(d) method. The absorption and fluorescence emission spectra were then evaluated via the time-dependent density functional theory B3LYP and PBE1PBE methods with 6-311++G(3df,2p) basis set. The calculation results show that compared with those of their parent molecules A-H, B-H, C-H, D-H, the absorption and emission wavelength values of all the derivatives show red shifts. The derivatives containing both the oxadiazole and methoxyl units are good candidates for longer absorption wavelength materials. The effects of azo, oxadiazole, and methoxyl units on the absorption and emission wavelength were discussed.     

7-Ethyl-2,3,5,6,8-pentahydroxy-1,4-naphthoquinone (echinochrome A): a DFT study of the mechanism of the antioxidant action. 4. Migration of protons and the Na<Superscript>+</Superscript> cation in echinochrome a monosodium salts

According to B3LYP/6-31G(d,p) and B3LYP/6-31G(d) calculations, intramolecular migrations of H⁺ and Na⁺ cations in the molecule of 7-ethyl-2,3,5,6,8-pentahydroxy-1,4-naphthoquinone (echinochrome A) monosodium salt in the gas phase and in clusters containing five water molecules can proceed asynchronously through the chain of O-H...O hydrogen bonds by the jump-over mechanism.     

Doubly hybrid density functional xDH-PBE0 from a parameter-free global hybrid model PBE0

Following the XYG3 model which uses orbitals and density from B3LYP, an empirical doubly hybrid (DH) functional is developed by using inputs from PBE0. This new functional, named xDH-PBE0, has been tested on a number of different molecular properties, including atomization energies, bond dissociation enthalpies, reaction barrier heights, and nonbonded interactions. From the results obtained, xDH-PBE0 not only displays a significant improvement with respect to the parent PBE0, but also shows a performance that is comparable to XYG3. Arguably, while PBE0 is a parameter-free global hybrid (GH) functional, the B3LYP GH functional contains eight fit parameters. From a more general point of view, the present work points out that reliable and general-purpose DHs can be obtained with a limited number of fit parameters.     

Absolute thermodynamic measurements of alkali metal cation interactions with a simple dipeptide and tripeptide

Absolute bond dissociation energies (BDEs) of glycylglycine (GG) and glycylglycylglycine (GGG) to sodium and potassium cations and sequential bond energies of glycine (G) in Na+G2 were determined experimentally by threshold collision-induced dissociation (TCID) in a guided ion beam tandem mass spectrometer. Experimental results showed that the binding energies follow the order of Na+ > K+ and M+GGG > M+GG > M+G. Theoretical calculations at the B3LYP/6-311+G(d) level show that all complexes had charge-solvated structures (nonzwitterionic) with either [CO,CO] bidentate or [N,CO,CO] tridentate coordination for M+GG complexes, [CO,CO,CO] tridentate or [N,CO,CO,CO] tetradentate coordination for M+GGG complexes, and [N,CO,N,CO] tetradentate coordination for Na+G2. Ab initio calculations at three different levels of theory (B3LYP, B3P86, and MP2(full) using the 6-311+G(2d,2p) basis set with geometries and zero-point energies calculated at the B3LYP/6-311+G(d) level) show good agreement with the experimental bond energies. This study demonstrates for the first time that TCID measurements of absolute BDEs can be successfully extended to biological molecules as complex as a tripeptide.