Density functional crystal vs. cluster models as applied to zeolites

In order to compare solid and cluster models of zeolites, we have studied the substitution Si⁴⁺→Al³⁺+H⁺ on the T₁ site of mordenite in the dilute limit using a self-consistent, full potential, local density functional (LDF) approach. Clusters size ranged from 9 to 105 atoms. Two crystal models with different Al concentrations were used. The first contained one substitution site per primitive cell of 72 atoms, the other one per conventional cell, containing 144 atoms. The unrelaxed substitution energies as computed with cluster and crystal models correspond well if the cluster results are extrapolated to infinite radius. Size effects are much smaller in crystal models. In addition, a structure relaxation (with fixed unit cell) was carried out for pure-silica offretite, a zeolite with 54 atoms per unit cell, and pure-silica mordenite, with 144 atoms per unit cell, starting from the low aluminum content X-ray crystallographic structure. In the offretite and mordenite optimizations full use was made of the D_3h¹–P $ \bar 6 $m2 and the nonsymmorphic D_2h¹⁷–Cmcm space group symmetry, respectively. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 68: 135–144, 1998     

Method for evaluation of density functional integrals in molecular calculations

Numerical methods for computing variationally optimized molecular orbitals within the Hartree–Fock approximation are augmented to include correlation functionals of the density in the energy and the numerical methods for carrying this out are described. The approach is applied explicitly to the Colle–Salvetti correlation energy functional. It is found that the gradient terms in the Colle–Salvetti functional present numerical problems associated with the low-density behavior, but also that they make a relatively small contribution to the correlation energy. In the three cases considered, HF, H₂O and N₂, it is found that the Colle–Salvetti correction considerably underestimates the correlation energies obtained in coupled-cluster theory.     

A theoretical study on the cyclopropane adsorption onto the copper surfaces by density functional theory and quantum chemical molecular dynamics methods

We report a theoretical study on the cyclopropane adsorption onto Cu(1 1 1) surfaces by density functional theory (DFT) and quantum chemical molecular dynamics methods. The equilibrium geometry of the physisorbed species was obtained using both periodic and cluster models by DFT methods that employ Cambridge serial total energy package (CASTEP), DMol ab initio quantum chemistry software of Accelrys’ materials studio (DMol), and Amsterdam density functional (ADF) program. It was found that the adsorbate molecule was tilted towards the metal surface with one C---C bond (upwards) parallel to the surface and that the physisorption occurred via a third carbon atom pointing (downwards) towards the surface. The electronic distribution and geometrical structure of physisorbed cyclopropane were slightly deviated from its gas phase molecule. The calculated vibrational frequencies and adsorption energies are close to experimental data, confirming the reliability of our DFT results. The adsorption process was simulated using our novel tight-binding quantum chemical molecular dynamics program, ‘Colors’. The calculation results indicated that both the adsorption and desorption processes of cyclopropane took place molecularly. The electron transfer and structural properties of equilibrium position obtained by ‘Colors’ are consistent with those by the first principles DFT methods.     

Approximation of the molecular electrostatic potential in a gaussian density functional method

The recently developed Asymptotic Density Model (ADM) [6, 9] is here implemented in the density functional framework using the program deMon-KS [13]. While the original implementation divided the atoms into a core shell and a valence shell, the present version allows for an arbitrary number of shells making it therefore more flexible and, as shown with benzene, potentially more accurate. Moreover, since this method is derived through Poisson's equation, an expression for the electronic charge density is also obtained. However, the present discussion will restrict itself to the electrostatic potential. Finally, even though this method requires parametrization, it is shown that the parameters obtained for homonuclear diatomic species, and used as is in molecular calculations, yield satisfactory results. Indeed, the ADM reproduces almost all basic features of the MEP for all molecules presented here, (water, ammonia, ethylene, acetylene, hydrogen cyanide, carbon monoxide, benzene, nitrous acid). Received: 5 July 1996 / Accepted: 12 November, 1996     

Determination of extremely localized molecular orbitals in the framework of density functional theory

Extremely localized molecular orbitals are rigorously localized on only a preselected set of atoms and do not have any tails outside the localization region. The importance of these orbitals lies in their ability to be transferred from one molecule to another one. A new algorithm to determine extremely localized molecular orbitals in the framework of the density functional theory method is presented. This could also be a valuable tool in the quantum mechanics/molecular mechanics methodology where localized molecular orbitals are used to describe covalent bonds across the frontier region. The present approach is used to build up the electron density of thymopentin, a polypeptide constituted by five residues, starting from extremely localized molecular orbitals determined on a set of model molecules. The results obtained confirm good transferability properties for these orbitals.Proceedings of the 11th International Congress of Quantum Chemistry satellite meeting in honor of Jean-Louis Rivail     

Local and nonlocal density functional calculations of the molecular structure of isomeric thiadiazoles

The chemistry of thiadiazoles and their derivatives is of considerable interest in chemistry owing to their pharmacological and potential industrial applications. In this context, a detailed study of isomeric thiadiazole molecules has been done using local (SVWN; Slater, and Vosko, Wilk and Nusair) and nonlocal (BLYP; Becke, and Lee, Yang and Parr) density functionals and optimizing the molecular geometries by means of the gradient technique. A charge sensitivity analysis of the studied molecule has been performed by resorting to density functional theory, obtaining several sensitivity coefficients such as the molecular energy, net atomic charges, global and local hardness, global and local softness and Fukui functions. With these results and the analysis of the dipole moments, the molecular electrostatic potentials and the total electron density maps, several conclusions have been inferred about the preferred sites of chemical reaction of the studied compounds. The condensed Fukui functions are shown to be one of the best criteria for predicting chemical reactivity.     

Ab initio and density functional theory calculations of molecular structure and vibrational spectrum of ethyl azidoacetate

Here we report ab initio and density functional results for molecular properties of ethyl azidoacetate (N₃CH₂COOC₂H₅) and for the corresponding singly ionized structure (N₃CH₂COOC₂H₅⁺). Ab initio ionization energies based on Koopmans’ theorem are in excellent agreement with the experimental data from ultraviolet photoelectron spectroscopy. DFT adiabatic energy differences between neutral and ionized structures are very sensitive to electronic correlation effects and are not in very good agreement with experiment. The results for the structure and vibrational frequencies are compared with the experimental data of related molecular structures.     

Study on electronic density topology of various cluster models of Mg/Al hydrotalcite by density functional theory

The geometry and electronic topology properties of Mg/Al hydrotalcite cluster models were comparatively investigated by means of density functional theory at GGA/DND levels. The results suggested that cluster model containing seven octahedral cations was the smallest size to be employed to simulate other properties. The fact that the n+ charge of cluster models containing n aluminum atoms can reflect electronic properties of anionic clay layer sheet. The bond lengths of clusters can be modified by terminating with or without OH~-/H_2O groups in terms of principle of bond order conservation.     

Chemisorption of OCN on Cu (100) surface: a density functional study

The chemisorption of cyanato radial (OCN) on Cu (100) surface is studied by using density functional theory (DFT) and the cluster model method. Cu₁₄ cluster is used to simulate the surface. Vertical bonding geometries with the nitrogen or oxygen atom down, and a parallel bonding geometry are considered, respectively. The present calculations show that cyanato-N species absorbed on the surface is more favorable than the other configurations. It indicates that OCN species is linearly bonded to the Cu (100) surface via the nitrogen atom, and is in good agreement with the experimental result. The cyanato-N species at the bridge site is most stable. For cyanato-N, the calculated symmetric and asymmetric OCN stretch frequencies are all blue-shifted compared with the calculated gaseous values, which is consistent with the experiment result. The charge transfer from the surface to OCN causes a work function increase on the surface. Bonding of OCN to the metal surface is largely ionic.     

Structure and ionization energies of some analogues of iron-only hydrogenases studied by density functional theory methods

Density functional theory calculations using both the B3LYP and BP86 functional in conjunction with a medium and large size basis set have been used to predict the structures and ionization energies of 12 models of iron-only hydrogenases. Although the structural predictions do not allow a clear discrimination between the different computational models, these models do yield significantly different adiabatic and vertical ionization energies. The closest agreement with experiment is given by the BP86 functional and the large all-electron basis. At this level of theory the adiabatic ionization energies are very close to experiment, but the vertical values are uniformly too small, leading to an underestimation of the reorganization energies. The calculations also suggest that measured ionization energies may help in identifying both the bridge-head group and whether CO bridging takes place upon ionization.     

Energy landscapes of nucleophilic substitution reactions: a comparison of density functional theory and coupled cluster methods

We have carried out a detailed evaluation of the performance of all classes of density functional theory (DFT) for describing the potential energy surface (PES) of a wide range of nucleophilic substitution (SN2) reactions involving, amongst others, nucleophilic attack at carbon, nitrogen, silicon, and sulfur. In particular, we investigate the ability of the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA as well as hybrid DFT to reproduce high-level coupled cluster (CCSD(T)) benchmarks that are close to the basis set limit. The most accurate GGA, meta-GGA, and hybrid functionals yield mean absolute deviations of about 2 kcal/mol relative to the coupled cluster data, for reactant complexation, central barriers, overall barriers as well as reaction energies. For the three nonlocal DFT classes, the best functionals are found to be OPBE (GGA), OLAP3 (meta-GGA), and mPBE0KCIS (hybrid DFT). The popular B3LYP functional is not bad but performs significantly worse than the best GGA functionals. Furthermore, we have compared the geometries from several density functionals with the reference CCSD(T) data. The same GGA functionals that perform best for the energies (OPBE, OLYP), also perform best for the geometries with average absolute deviations in bond lengths of 0.06 A and 0.6 degrees, even better than the best meta-GGA and hybrid functionals. In view of the reduced computational effort of GGAs with respect to meta-GGAs and hybrid functionals, let alone coupled cluster, we recommend the use of accurate GGAs such as OPBE or OLYP for the study of SN2 reactions.     

Selective adsorption of bacterial cells onto zeolites

Zeolites adsorb microbial cells on their surfaces and selective adsorption for specific microorganisms was seen with certain zeolites. Tests for the adsorption ability of zeolites were conducted using various established microbial cell lines. Specific cell lines were shown to selectively absorb to certain zeolites, species to species.

In order to understand the selectivity of adsorption, we tested adsorption under various pH conditions and determined the zeta-potentials of zeolites and cells. The adsorption of some cell lines depended on the pH, and some microorganisms were preferentially adsorbed at acidic pH. The values of zeta-potentials were used for calculating the electric double layer interaction energy between zeolites and microbial cells. There was a correlation between the experimental adsorption results and the interaction energy. Moreover, we evaluated the surface hydrophobicity of bacterial cells by using the microbial adherence to hydrocarbon (MATH) assay. In addition, we also applied this method for zeolites to quantify relative surface hydrophobicity. As a result, we found a correlation between the adsorption results and the hydrophobicity of bacterial cells and zeolites. These results suggested that adsorption could be explained mainly by electric double layer interactions and hydrophobic interactions.

Finally, by using the zeolites Na-BEA and H-Y, we succeeded in clearly separating three representative microbes from a mixture of Escherichia coli, Bacillus subtilis and Staphylococcus aureus. Zeolites could adsorb each of the bacterial cell species with high selectivity even from a mixed suspension. Zeolites can therefore be used as effective carrier materials to provide an easy, rapid and accurate method for cell separation.     

The evaluation of bond dissociation energies for simple sulphur containing molecules using ab initio and density functional methods

The S–H and C–S bond dissociation energies for simple alkylthiols and dialkylsulphides, along with the S–S bond dissociation energy for dimethyl disulphide, compounds which have been used in the metal–organic chemical vapour deposition (MOCVD) growth of wide band gap II–VI (12–16) Zn- and Cd-based compound semiconductors, have been computed using the ab initio (ROHF and MP2) and density functional theory (DFT) methods (BHandH, BHandHLYP, B3LYP, B3P86, B3PW91, BLYP and BP86) with the 6-311+G(2d,p) basis set along with high accuracy complete basis set, CBS-4 and CBS-Q energy computations. The computed energies are compared with experimental results and the suitability of the DFT methods, for the computational study of these systems, is discussed.     

Excited state geometries within time-dependent and restricted open-shell density functional theories

Singlet excited state geometries of a set of medium sized molecules with different characteristic lowest excitations are studied. Geometry optimizations of excited states are performed with two closely related restricted open-shell Kohn–Sham methods and within linear response to time-dependent density functional theory. The results are compared to wave-function based methods. Excitation energies (vertical and adiabatic) calculated from the open-shell methods show systematic errors depending on the type of excitation. However, for all states accessible by the restricted methods a good agreement for the geometries with time-dependent density functional theory and wave-function based methods is found. An analysis of the energy with respect to the mixing angle for the singly occupied orbitals reveals that some states (mostly [n→π^*]) are stable when symmetry constraints are relaxed and others (mostly [π→π^*]) are instable. This has major implications on the applicability of the restricted open-shell methods in molecular dynamics simulations.     

A short-range correlation energy density functional with multi-determinantal reference

We introduce a short-range correlation density functional defined with respect to a multi-determinantal reference which is meant to be used in a multi-determinantal extension of the Kohn–Sham scheme of density functional theory based on a long-range/short-range decomposition of the Coulomb electron–electron interaction. We construct the local density approximation for this functional and discuss its performance on the He atom.     

Stable isomers for trifluoroacetic acid (TFA) pentahydrates obtained from density functional calculations

Stable isomers of trifluoroacetic acid (TFA) pentahydrate clusters, TFA-(H₂O)₅, have been explored by using density functional theory calculations. As done for TFA-(H₂O)₄ (Ito, 2013), structure optimization and vibrational calculations were performed for 70 isomeric structures (68 for neutral and 2 for ion-pair species, respectively) at the B971/6–311++G(3df,3pd) level. We found that the edge-sharing bicyclic isomer is at the global minimum and that three other isomers lie energetically within 100 cm⁻¹. Two types of ion-pair species were found to be unstable by 1100 cm⁻¹ in comparison with the global minimum. The results were compared with infrared spectra observed in nitrogen matrix.     

The Gaussian and augmented-plane-wave density functional method for ab initio molecular dynamics simulations

A new algorithm for density-functional-theory-based ab initio molecular dynamics simulations is presented. The Kohn–Sham orbitals are expanded in Gaussian-type functions and an augmented-plane-wave-type approach is used to represent the electronic density. This extends previous work of ours where the density was expanded only in plane waves. We describe the total density in a smooth extended part which we represent in plane waves as in our previous work and parts localised close to the nuclei which are expanded in Gaussians. Using this representation of the charge we show how the localised and extended part can be treated separately, achieving a computational cost for the calculation of the Kohn–Sham matrix that scales with the system size N as O(NlogN). Furthermore, we are able to reduce drastically the size of the plane-wave basis. In addition, we introduce a multiple-cutoff method that improves considerably the performance of this approach. Finally, we demonstrate with a series of numerical examples the accuracy and efficiency of the new algorithm, both for electronic structure calculations and for ab initio molecular dynamics simulations. Received: 15 December 1998 /Accepted: 18 February 1999 /Published online: 14 July 1999