Computational modeling of polyoxotungstates by relativistic DFT calculations of (183)W NMR chemical shifts

The (183)W nuclear shielding in a variety of tungsten polyoxometalates (POM) (Lindqvist, Anderson, decatungstates, Keggin) of different shapes and charges has been modeled by DFT calculations that take into account relativistic effects, by means of the zero-order regular approximation (ZORA), and solvent effects, by the conductor-like screening model (COSMO) continuum method. The charge/surface area ratio (q/A) is proposed as an indicator of the charge density to which the solvation energies of all POMs are correlated in a satisfactory way. Among the various theoretical levels tested (ZORA scalar or spin-orbit, frozen-core or all-electron basis set, geometry optimization in the gas phase or in the continuum solvent, etc.), the best results are obtained when both geometry optimization in solvent and spin-orbit shielding are included (mean absolute error of delta=35 ppm). The quality of the computed chemical shifts depends systematically on the charge density as expressed by q/A; thus, POMs with low q/A ratios display the best agreement with experimental data. The performance of the method is such that computed values can aid the assignment of the (183)W NMR spectra of polyoxotungstates, as shown by the case of alpha-[PW(11)TiO(40)](5-), whose six signals are ranked computationally so as to almost reproduce the experimental ordering even though the signals are spaced by as little as 5 ppm.     

All-electron density functional calculation on insulin with quasi-canonical localized orbitals

An all-electron density functional (DF) calculation on insulin was performed by the Gaussian-based DF program, ProteinDF. Quasi-canonical localized orbitals (QCLOs) were used to improve the initial guess for the self-consistent field (SCF) calculation. All calculations were carried out by parallel computing on eight processors of an Itanium2 cluster (SGI Altix3700) with a theoretical peak performance of 41.6 GFlops. It took 35 h for the whole calculation. Insulin is a protein hormone consisting of two peptide chains linked by three disulfide bonds. The numbers of residues, atoms, electrons, orbitals, and auxiliary functions are 51, 790, 3078, 4439, and 8060, respectively. An all-electron DF calculation on insulin was successfully carried out, starting from connected QCLOs. Regardless of a large molecule with complicated topology, the differences in the total energy and the Mulliken atomic charge between initial and converged wavefunctions were very small. The calculation proceeded smoothly without any trial and error, suggesting that this is a promising method to obtain SCF convergence on large molecules such as proteins.     

Elastic bag model for molecular dynamics simulations of solvated systems: application to liquid argon

A new approach is developed to study the dynamics of the localized process in solutions and other condensed phase systems. The approach employs a fluctuating elastic boundary (FEB) model which encloses the simulated system in an elastic bag that mimics the effects of the bulk solvent. This alleviates the need for periodic boundary conditions and allows for a reduction in the number of solvent molecules that need to be included in the simulation. The boundary bag is modeled as a mesh of quasi-particles connected by elastic bonds. The FEB model allows for volume and density fluctuations characteristic of the bulk system, and the shape of the boundary fluctuates during the course of the simulation to adapt to the configuration fluctuations of the explicit solute-solvent system inside. The method is applied to the simulation of a Lennard-Jones model of liquid argon. Various structural and dynamical quantities are computed and compared with those obtained from conventional periodic boundary simulations. The agreement between the two is excellent in most cases, thus validating the viability of the FEB method.     

Simulation of IR and Raman spectra of p-hydroxyanisole and p-nitroanisole based on scaled DFT force fields and their vibrational assignments

This work deals with the vibrational spectroscopy of p-hydroxyanisole (PHA) and p-nitroanisole (PNA) by means of quantum chemical calculations. The mid and far FT-IR and FT-Raman spectra were recorded in the condensed state. The fundamental vibrational frequencies and intensity of vibrational bands were evaluated using density functional theory (DFT) with the standard B3LYP/6-31G* method and basis set combination and were scaled using various scale factors which yield a good agreement between observed and calculated frequencies. The vibrational spectra were interpreted with the aid of normal coordinate analysis based on scaled density functional force field. The results of the calculations were applied to simulate infrared and Raman spectra of the title compounds, which showed excellent agreement with the observed spectra.     

Observing the 3D Chemical Bond and its Energy Distribution in a Projected Space

Our curiosity-driven desire to “see” chemical bonds dates back at least one-hundred years, perhaps to antiquity. Sweeping improvements in the accuracy of measured and predicted electron charge densities, alongside our largely bondcentric understanding of molecules and materials, heighten this desire with means and significance. Here we present a method for analyzing chemical bonds and their energy distributions in a two-dimensional projected space called the condensed charge density. Bond “silhouettes” in the condensed charge density can be reverse-projected to reveal precise three-dimensional bonding regions we call bond bundles. We show that delocalized metallic bonds and organic covalent bonds alike can be objectively analyzed, the formation of bonds observed, and that the crystallographic structure of simple metals can be rationalized in terms of bond bundle structure. Our method also reproduces the expected results of organic chemistry, enabling the recontextualization of existing bond models from a charge density perspective.     

Conformational and vibrational study of platinum(II) anticancer drugs: cis-diamminedichloro-platinum(II) as a case study

A conformational and vibrational analysis of cisplatin [cis-diamminedichloro-platinum(II)] (cDDP) is reported. Several theory methods (from Hartree-Fock to Moller-Plesset and density functional theory) combined with different all-electron basis sets are evaluated, in view of determining the best suited strategy for accurately representing this molecule. This choice is based on the best compromise between accuracy and computational requirements. Different scaling models of the cDDP vibrational modes were tested for obtaining the best scaling factors to be used in this type of inorganic systems. The structural parameters and vibrational results predicted by the calculations are compared with the corresponding experimental data, namely, x-ray structure, and Raman and inelastic neutron scattering spectra. Finally, a complete assignment of the cDDP vibrational spectra is presented.     

Dimetalloendofullerene U(2)@C(60) has a U-U multiple bond consisting of sixfold one-electron-two-center bonds

Endohedral metallofullerenes (EMFs) have been extensively studied since their discovery in 1985. Metal-metal bonds, nevertheless, have never been explicitly observed in EMFs synthesized so far. In this contribution, we show by means of all-electron relativistic density functional computations that the dimetalloendofullerene, U(2)@C(60), has an unprecedented U-U multiple bond consisting solely of sixfold ferromagnetically coupled one-electron-two-center bonds with the electronic configuration (5fpi(u))(2)(5fsigma(g))(1)(5fdelta(g))(2)(5fphi(u))(1), which are dominated by the uranium 5f atomic orbitals. This bonding scheme is completely distinct from the metal-metal bonds discovered thus far in the d- and f-block polynuclear metal complexes. This finding initiates a connection of the metal-metal multiple bonding chemistry and the fullerene chemistry.     

Theoretical binding affinities and spectra of complexes formed by a cyclic β-peptoid with amino acids

Binding affinities of a cyclic β-peptoid to amino acids were studied using the density functional theory (DFT) at the B3LYP/6-311G(d,p) level after the basis set superior error (BSSE). The host molecule possesses binding ability to amino acids since the binding energies of the complexes formed are negative. The complexes were stabilized via hydrogen bonds between the host and the guest molecules. Based on the B3LYP/6-31G(d) optimized geometries, electronic spectra of the complexes were calculated using the INDO/CIS method. ¹³C NMR spectra and nucleus-independent chemical shift (NICS) values of the complexes were computed at the B3LYP/6-31G(d) level. Carbon atoms in the carboxyl groups of the complexes are shifted downfield relative to those of the host. Some complexes exhibit aromaticity although the host shows anti-aromaticity. Formation of hydrogen bonds leads to cyclic current formation in these complexes.     

A first principles theory of nuclear magnetic resonance J-coupling in solid-state systems

A method to calculate NMR J-coupling constants from first principles in extended systems is presented. It is based on density functional theory and is formulated within a planewave-pseudopotential framework. The all-electron properties are recovered using the projector augmented wave approach. The method is validated by comparison with existing quantum chemical calculations of solution-state systems and with experimental data. The approach has also been applied to the silicophosphate, Si(5)O(PO(4))(6), giving (31)P-(29)Si-couplings which are in excellent agreement with experiment.     

The chemical bond between Au(I) and the noble gases. Comparative study of NgAuF and NgAu+ (Ng = Ar, Kr, Xe) by density functional and coupled cluster methods

The nature of the chemical bond between gold and the noble gases in the simplest prototype of Au(I) complexes (NgAuF and NgAu+, where Ng = Ar, Kr, Xe), has been theoretically investigated by state of art all-electron fully relativistic DC-CCSD(T) and DFT calculations with extended basis sets. The main properties of the molecules, including dipole moments and polarizabilities, have been computed and a detailed study of the electron density changes upon formation of the Ng-Au bond has been made. The Ar-Au dissociation energy is found to be nearly the same in both Argon compounds. It almost doubles along the NgAuF series and nearly triples in the corresponding NgAu+ series. The formation of the Ng-Au(I) bonds is accompanied by a large and very complex charge redistribution pattern which not only affects the outer valence region but reaches deep into the core-electron region. The charge transfer from the noble gas to Au taking place in the NgAu+ systems is largely reduced in the fluorides but the Ng-Au chemical bond in the latter systems is found to be tighter near the equilibrium distance. The density difference analysis shows, for all three noble gases, a qualitatively identical nature of the Ng-Au bond, characterized by the pronounced charge accumulation in the middle of the Ng-Au internuclear region which is typical of a covalent bond. This bonding density accumulation is more pronounced in the fluorides, where the Au-F bond is found to become more ionic, while the overall density deformation is more evident and less localized in the NgAu+ systems. Accurate density difference maps and charge-transfer curves help explain very subtle features of the chemistry of Au(I), including its peculiar preference for tight linear bicordination.     

苝二酸酐与嘧啶衍生物的氢键组装

Semi-empirical AM1 method was used to study 1:1 and 1:2 hydrogen bond complexes formed with perylene dianhydride and pyridine derivatives. The weak interaction energy become bigger as the number of hydrogen bonds increases. The donor groups on the host and electron-withdrawing groups on the guest molecules favor hydrogen bonding interactions, and the formation of hydrogen bonding leads to electron density flow from the host to the guest molecules. Electronic spectra of these complexes were computed using INDO/SCI method. Blue-shift of the electronic absorption spectra for the complexes, comparing that of the host,takes place, and the first peaks for different complexes changed slightly. These are in agreement with the experimental results. The cause of blue-shift was discussed, and the electronic transitions were assigned based on theoretical calculations. The potential curve of double proton transfer in the complex was calculated, and the transition state and activated energy relative to the N-H bond were obtained.     

All-electron time-dependent density functional theory with finite elements: time-propagation approach

We present an all-electron method for time-dependent density functional theory which employs hierarchical nonuniform finite-element bases and the time-propagation approach. The method is capable of treating linear and nonlinear response of valence and core electrons to an external field. We also introduce (i) a preconditioner for the propagation equation, (ii) a stable way to implement absorbing boundary conditions, and (iii) a new kind of absorbing boundary condition inspired by perfectly matched layers.     

Simulating infrared spectra and hydrogen bonding in cellulose Iβ at elevated temperatures

We have modeled the transformation of cellulose Iβ to a high temperature (550 K) structure, which is considered to be the first step in cellulose pyrolysis. We have performed molecular dynamics simulations at constant pressure using the GROMOS 45a4 united atom forcefield. To test the forcefield, we computed the density, thermal expansion coefficient, total dipole moment, and dielectric constant of cellulose Iβ, finding broad agreement with experimental results. We computed infrared (IR) spectra of cellulose Iβ over the range 300-550 K as a probe of hydrogen bonding. Computed IR spectra were found to agree semi-quantitatively with experiment, especially in the O-H stretching region. We assigned O-H stretches using a novel synthesis of normal mode analysis and power spectrum methods. Simulated IR spectra at elevated temperatures suggest a structural transformation above 450 K, a result in agreement with experimental IR results. The low-temperature (300-400 K) structure of cellulose Iβ is dominated by intrachain hydrogen bonds, whereas in the high-temperature structure (450-550 K), many of these transform to longer, weaker interchain hydrogen bonds. A three-dimensional hydrogen bonding network emerges at high temperatures due to formation of new interchain hydrogen bonds, which may explain the stability of the cellulose structure at such high temperatures.     

Accurate ab initio structure, dissociation energy, and vibrational spectroscopy of the F(-)-CH4 anion complex

Accurate equilibrium structure, dissociation energy, global potential energy surface (PES), dipole moment surface (DMS), and the infrared vibrational spectrum in the 0-3000 cm(-1) range of the F(-)-CH4 anion complex have been obtained. The equilibrium electronic structure calculations employed second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster (CC) method up to single, double, triple, and perturbative quadruple excitations using the aug-cc-p(C)VXZ [X = 2(D), 3(T), 4(Q), and 5] correlation-consistent basis sets. The best equilibrium geometry has been obtained at the all-electron CCSD(T)/aug-cc-pCVQZ level of theory. The dissociation energy has been determined based on basis set extrapolation techniques within the focal-point analysis (FPA) approach considering (a) electron correlation beyond the all-electron CCSD(T) level, (b) relativistic effects, (c) diagonal Born-Oppenheimer corrections (DBOC), and (d) variationally computed zero-point vibrational energies. The final D(e) and D0 values are 2398 +/- 12 and 2280 +/- 20 cm(-1), respectively. The global PES and DMS have been computed at the frozen-core CCSD(T)/aug-cc-pVTZ and MP2/aug-cc-pVTZ levels of theory, respectively. Variational vibrational calculations have been performed for CH4 and F(-)-CH4 employing the vibrational configuration interaction (VCI) method as implemented in Multimode.     

Tuning the interaction energy of hydrogen bonds: the effect of the substituent

The effect of the substituent R in the hydrogen bonding properties of FH···FR (R = H, Al, Li, Cl and CCH) complexes has been studied by theoretical calculations. The dependency of the interaction energy with the hydrogen bond distance and R is explained in terms of the topologies of the electron density and the electrostatic potential. A simple model of the hydrogen bond interaction energy, which can be assimilated to an interaction potential, is defined in terms of a stabilizing mutual polarization of the monomers and an overall destabilizing contribution associated with the electron density reorganization when the overlap of the closed shells is large enough. This model shows an excellent agreement with the ab initio interaction energies and is common for all the analyzed complexes. The substituent effect is represented in the model by a single parameter that can be calculated from the electron distribution in the acceptor atom region. The perturbation in the hydrogen bonding interaction induced by the change of R presents a close similarity with that produced by an external electric field of the same order of magnitude than those found in crystalline solids, indicating that both perturbations should play a significant and similar role on the properties of hydrogen bonds in condensed matter.