Local-orbital-based correlated ab initio band structure calculations in insulating solids: LiF

 A local-orbital-based ab initio approach to calculate correlation effects on quasi-particle energies in insulating solids is presented. The use of localized Wannier-type Hartree–Fock orbitals allows correlation effects to be efficiently assessed. First a Green's function approach based on exact diagonalization is introduced and this is combined with an incremental scheme, while subsequently different levels of perturbative approximations are derived from the general procedure. With these methods the band structure of LiF is calculated and good agreement with experiment is found. By comparing the different approximations proposed, including the exact diagonalization procedure, their relative quality is established. Received: 25 June 2001 / Accepted: 31 August 2001 / Published online: 19 December 2001     

Role of twisting and sliding on the solvation of a stacked cytosine dimer: an ab initio study

Ab initio calculations with inclusion of correlation effects at the MP2/6-31G* level have been used to predict the interaction energy of stacked cytosine dimer (C/C) as a function of twisting and sliding in the gas phase. Systematic calculations have also been carried out on the solvation free energies of various rotated and translated C/C dimers using a polarized continuum model approach at the HF/6-31G* level with a view to probe the role of various degrees of freedom on the free energy of solvation of the C/C dimer. The interaction energy of the C/C dimer decreases upon changing from a parallel to an antiparallel conformation in the gas phase. The 180°-rotated conformation has been found to be the most stable arrangement when compared to other rotated positions. The rotated and translated dimers exhibit lower solvation free energy than the parallel conformation. The decrease in the dipole moment upon rotation from the parallel to the antiparallel conformation indicates the cancellation of charge distribution upon rotation in the z direction of one cytosine base with respect to the other. The calculation reveals that the present approach could not yield association energy, ΔΔG _Asso, in a solvent medium. This may be due to the fact that in the case of floppy molecules the contribution from translational, rotational and vibrational free energies plays a significant role in the calculation of ΔΔG _Asso. Received: 13 December 2001 / Accepted: 25 March 2002 / Published online: 13 June 2002     

Direct ab initio dynamics calculations on the rate constants for the hydrogen-abstraction reaction of C2H5F with O (3P)

The hydrogen-abstraction reaction C₂H₅F+O → C₂H₄F+OH has been studied by a dual-level direct dynamics method. For the reaction, three reaction channels, one for α-abstraction and two for β-abstraction, have been identified. The potential-energy surface information is obtained at the MP2(full)/6-311G(d,p) and PMP2(full)/6-311G(3df,3pd) (single-point) levels. By canonical variational transition-state theory, rate constants for each reaction channel are calculated with a small-curvature tunneling correction. The total rate constant is calculated from the sum of the individual rate constants and the temperature dependence of the branching ratios is obtained over a wide range of temperatures from 300 to 5,000 K. The agreement of the rate constants with experiment is good in the experimental temperature range from 1,000 to 1,250 K. The calculated results indicate that at low temperatures α-abstraction is most likely to be the major reaction channel, while β-abstraction channels will significantly contribute to the whole reaction rate as the temperature increases. Received: 23 January 2002 / Accepted: 23 June 2002 / Published online: 20 September 2002     

Ab initio study of bonding trends for f0 actinide oxyfluoride species

 Fully relativistic, four-component Dirac–Fock calculations and quasirelativistic pseudopotential calculations at different ab initio levels are used to study the bonding trends among the naked, triatomic [OAnO] ^q+ groups or the oxyfluorides [AnO _n F _m ] ^q with f ⁰ configurations. The triatomic f ⁰ series is suggested to range from the bent ThO₂ via the linear OPaO⁺ to at least NpO₂ ³⁺, a possible new gas-phase species. The neutral oxyfluoride molecules include the experimentally unknown NpO₂F₃ and PuO₂F₄. The latter is a candidate for the so far unknown oxidation state Pu(VIII), which is found to lie considerably above Pu(VI), but to be locally stable. Their all-oxygen isoelectronic analogues are NpO₅ ³⁻, known in the solid state, and the unknown PuO₆ ⁴⁻. Further possible candidates for Pu(VIII) are PuO₄(D _4h ) and the cube-shaped PuF₈(O _h ). Isoelectronic UF₈ ²⁻ is calculated to be D _4d , in agreement with experiment. Received: 18 May 2001 / Accepted: 21 June 2001 / Published online: 11 October 2001     

Action of HCl on 3-hydroxypyrazolo(iso)quinolines to give 1-chloropyrazoles: evidence for an addition–elimination mechanism by ab initio calculations in gas phase and water

 Addition–elimination reactions involving a nucleophile and a remote leaving group [S^H _N(AE)^tele] are well-known under basic conditions, especially amongst electron-poor six-membered heterocycles, but are less commonly encountered for five-membered heterocycles and are rare under acidic conditions. Concentrated HCl converts 1-hydroxy-1H-pyrazolo[3,4-c] isoquinoline and 1-hydroxy-1H-pyrazolo[3,4-c]quinoline into 3-chloro-1H-pyrazolo[3,4-c]isoquinoline and 3-chloro-1H-pyrazolo[3,4-c]quinoline, respectively. However, apparently neither the isomeric 1-hydroxy-1H-pyrazolo[4,3-c](iso)-quinolines nor the parent 1-hydroxypyrazole undergo this reaction. Additionally, all these systems are refractory under basic conditions. We present a plausible mechanism for the reaction, involving the 3-addition of Cl^- to the diprotonated heterocycle, followed by the elimination of water. Calculations of the initial transition states and intermediates, using optimisation at B3LYP/6-311+G(d,p), including thermochemistry [HF/6-31+G(d)], and single-point Poisson–Boltzmann self-consistent reaction field determination of the free energy of solvation (Jaguar Poisson–Boltzmann self-consistent reaction field), support this mechanism and reproduce the observed order of reactivity, the addition step being 2–4 kcal less favourable for the isomeric 1-hydroxy-1H-pyrazolo[4,3-c](iso)quinolines and provide a rationalisation for the role of strong acid. Received: 27 June 2002 / Accepted: 6 September 2002 / Published online: 14 February 2003     

Ab initio molecular orbital study of Fe(CO) n (n = 1–3)

Geometry optimization was performed for the ground states of FeCO, Fe(CO)₂, and Fe(CO)₃ at various levels of ab initio calculations, and the bond lengths and dissociation energies obtained were in reasonable agreement with experimental results. The nature of bonding was studied for these molecules using a complete-active-space self-consistent-field method. From the Mulliken population analysis, it was found that the traditional donation and back donation mechanism is valid for these molecules, including Fe(CO)₃, which has a pyramidal structure. Received: 27 September 1999 / Accepted: 13 January 2000 / Published online: 19 April 2000     

Internal rotation in squaramide and related compounds. A theoretical ab initio study

The structural and energetic changes associated with C–N bond rotation in a squaric acid derivative as well as in formamide, 3-aminoacrolein and vinylamine have been studied theoretically using ab initio molecular orbital methods. Geometry optimizations at the MP2(full)/6-31+G* level confirmed an increase in the C–N bond length and a smaller decrease in the C=O length on going from the equilibrium geometry to the twisted transition state. Other geometrical changes are also discussed. Energies calculated at the QCISD(T)/6-311+G** level, including zero-point-energy correction, show barrier heights decreasing in the order formamide, squaric acid derivative, 3-aminoacrolein and vinylamine. The origin of the barriers were examined using the atoms-in-molecules approach of Bader and the natural bond orbital population analysis. The calculations agree with Pauling's resonance model, and the main contributing factor of the barrier is assigned to the loss of conjugation on rotating the C–N bond. Finally, molecular interaction potential calculations were used to study the changes in the nucleophilicity of N and O (carbonyl) atoms upon C–N rotation, and to obtain a picture of the abilities of the molecules to act in nonbonded interactions, in particular hydrogen bonds. The molecular interaction potential results confirm the suitability of squaramide units for acting as binding units in host–guest chemistry. Received: 13 March 2002 / Accepted: 23 June 2002 / Published online: 21 August 2002     

QM/QM study of the coverage effects on the adsorption of amino-cyclopentene at the Si(100) surface

In this work, we have tested 30 different adsorption situations in several coverage scenarios for the 1-amino-3-cyclopentene (ACP) molecule on the Si(100) surface. We have used a five-spot testing zone inserted in the high-level part of a quantum-mechanical/quantum-mechanical study performed in a big cluster. By defining several different scenarios, each one with a typical adsorption energy, we were able to understand in detail the process of surface functionalization. We are able to justify why the functionalization of this silicon surface achieves only a coverage of approximately 0.5 ML (half monolayer) and why the completely covered surface should be thermodynamically impossible to obtain.     

Evaluation of an ab initio quantum mechanical/molecular mechanical hybrid-potential link-atom method

 Hybrid potentials have become a common tool in the study of many condensed-phase processes and are the subject of much active research. An important aspect of the formulation of a hybrid potential concerns how to handle covalent bonds between atoms that are described with different potentials and, most notably, those at the interface of the quantum mechanical (QM) and molecular mechanical (MM) regions. Several methods have been proposed to deal with this problem, ranging from the simple link-atom method to more sophisticated hybrid-orbital techniques. Although it has been heavily criticized, the link-atom method has probably been the most widely used in applications, especially with hybrid potentials that use semiempirical QM methods. Our aim in this paper has been to evaluate the link-atom method for ab initio QM/MM hybrid potentials and to compare the results it gives with those of previously published studies. Given its simplicity and robustness, we find that the link-atom method can produce results of comparable accuracy to other methods as long as the charge distribution on the MM atoms at the interface is treated appropriately. Received: 27 September 2002 / Accepted: 21 October 2002 / Published online: 8 January 2003 Correspondence to: M. J. Field e-mail: mjfield@ibs.fr Acknowledgements. The authors thank the Institut de Biologie Structurale – Jean-Pierre Ebel, the Commissariat à l'Energie Atomique and the Centre National de la Recherche Scientifique for support of this work.     

Heteroelectrocyclic reactions of 3-acyl-4-azido heterocycles: ab initio and density functional calculations on 3-azido-propenal as a model system

 Concerted as well as stepwise reaction pathways for cyclization of 3-azido-propenal to isoxazole have been investigated by density functional (B3LYP) and ab initio methods up to CCSD(T)/cc-pVQZ methods. These calculations clearly establish the pathway with concerted albeit asynchronous nitrogen extrusion and ring closure as the most feasible mechanism. Barriers for cyclization increase in the order Hartree–Fock<B3LYP<ACPF<CCSD(T). According to the geometrical parameters and the electronic structure of the TS as evidenced by natural bond order analysis this cyclization can be interpreted as a pseudopericyclic (heteroelectrocyclic) reaction. Received: 26 May 2002 / Accepted: 18 June 2002 / Published online: 14 February 2003     

A genetic algorithm for the structural optimization of Morse clusters

This article describes the application of a genetic algorithm for the structural optimization of 19–50-atom clusters bound by medium-range and short-range Morse pair potentials. The GA is found to be efficient and reliable for finding the geometries corresponding to the previously published global minima [Doye JPK, Wales DJ (1997) J Chem Soc Faraday Trans 93: 4233]. Using the genetic algorithm, only a relatively small number of energy evaluations and minimizations are required to find the global minima. By contrast, a simple random search algorithm often cannot find the global minima of the larger clusters, even after many thousands of searches. Received: 27 October 1999 / Accepted: 7 December 1999 / Published online: 19 April 2000     

Computation of the influence of chemical substitution on the pK a of pyridine using semiempirical and ab initio methods

The physical properties of chemicals are strongly influenced by their protonation state, affecting, for example, solubility or hydrogen-bonding characteristics. The ability to accurately calculate protonation states in the form of pK _as is, therefore, desirable. Calculations of pK _a changes in a series of substituted pyridines are presented. Computations were performed using both ab initio and semiempirical approaches, including free energies of solvation via reaction-field models. The selected methods are readily accessible with respect to both software and computational feasibility. Comparison of calculated and experimental pK _as shows the experimental trends to be reasonably reproduced by the computations with root-mean-square differences ranging from 1.22 to 4.14 pK _a units. Of the theoretical methods applied the best agreement occurred using the second-order M?ller–Plesset/6-31G(d)/isodensity surface polarized continuum solvation model, while the more computationally accessible Austin model 1/Solvent model 2 (SM2) approach yielded results similar to the ab initio methods. Analysis of component contributions to the calculated pK _as indicates the largest source of error to be associated with the free energies of solvation of the protonated species followed by the gas-phase protonation energies; while the latter may be improved via the use of higher levels of theory, enhancements in the former require improvements in the solvation models. The inclusion of alternate minimum in the computation of pK _as is also indicated to contribute to differences between experimental and calculated pK _a values. Received: 27 April 1999 / Accepted: 27 July 1999 / Published online: 2 November 1999     

Infrared Spectroscopic Analysis of Plasma-Treated Si(100)-Surfaces

 Infrared reflection spectroscopy (specular reflection, attenuated total reflection) has been applied in combination with spectroscopic ellipsometry and electron microscopy to analyze the surface structure of plasma-treated Si(100) surfaces. It is shown that plasma treatments in oxygen and fluorine or chlorine-containing gases cause the formation of a thin surface layer having thicknesses of a few nanometers. The layer was identified to consist of SiO₂ for treatments in an oxygen plasma. Analyses of layers formed by treatments in a fluorine-containing plasma do not confirm the generally assumed model. Different Si-F vibration modes were identified in the surface layer caused by a SF₆ plasma. They correlate, however, with SiF and SiF₂ molecules. There are no indications of the existence of the generally assumed SiF₄. Neither has SiOF₂ been proven in layers produced by etching in a SF₆/O₂ plasma.     

Morphological evolution of Ge islands on the Si(100) surface: from huts to pits

The morphological evolution on the size and shape of Ge island on the Si(100) surface by depositing and subsequent annealing processes is studied in situ by using scanning tunneling microscopy at ultrahigh vacuum environment. A slower growth rate is maintained when the islands grow to larger sizes beyond the wetting layers. While at room temperature, the epitaxial strain is relieved by the formation of three‐dimensional islands (so‐called ‘hut’ clusters). When the sample is annealed at 200 °C, the strain is relieved by forming pits, having the circular cone shape but with their apex pointing down, with Ge clusters formed at the rim of pits. Copyright © 2016 John Wiley & Sons, Ltd.     

Theoretical study of the potential-energy surface of C2NP

 The B3LYP/6-311G(d) and CCSD(T)/6-311G(2df) (single-point) methods have been used to investigate the singlet potential energy surface of C₂NP, in which seven stationary isomers and seventeen interconversion transition states are involved. At the final CCSD(T)/6-311G(2df)//B3LYP6-311G(d) level with zero-point vibrational energy correction the lowest-lying isomer is a linear NCCP followed by two linear CNCP isomers at 23.9  and CCNP at 65.8 kcal mol⁻¹, respectively. The three isomers are kinetically very stable towards both isomerization and dissociation, and CCNP is even more kinetically stable than CNCP – by 14.3 kcal mol⁻¹ despite its high energy. Further comparative calculations were performed at the QCISD and QCISD(T) levels with the 6-311G(d) and 6-311G(2d) basis sets to obtain more reliable structures and spectroscopy for the three isomers. The calculated bond lengths, rotational constant, and dipole moment for NCCP were in close agreement with the experimentally determined values. Finally, similarities and discrepancies between the potential energy surface of C₂NP and those of the analogous species C₂N₂ and C₂P₂ were compared. The results presented in this paper might be helpful for future identification of the two still unknown yet kinetically very stable isomers CNCP and CCNP, both in the laboratory and in interstellar space. Received: 3 January 2001 / Accepted: 6 June 2001 / Published online: 30 October 2001     

The electron affinities of transition metal atoms at the CCSD(T) and density functional levels of theory

 The electron affinities of Ti, V, Cr, Fe, Co, Ni, and Cu are computed using the density function theory and CCSD(T) approaches. Overall the CCSD(T) approach yields the best results. For this property, the B3LYP, BLYP, and BP86 functionals perform better than the BPW91, PBEPBE, and PBE1PBE ones. The accuracy of all the methods is higher if the number of 3delectrons is the same in the neutral atom and the anion. This is especially true for the density functional theory methods. Received: 23 January 2002 / Accepted: 1 April 2002 / Published online: 24 June 2002     

Molecular dynamics study of the tautomeric equilibrium in the 4-nitro- and 2,4,6-trichloro derivatives of 2-(N,N-dialkyloaminomethyl)phenol

 Molecular dynamics thermodynamic integration (MDTI) method and quantum chemical calculations at the density functional theory B3LYP 6-31+(d,p) level, which included the Tomasi model of the solvent reaction field, were applied to study the tautomeric equilibrium of Mannich base in methanol solution. The values obtained for the free-energy difference are in good agreement with experimental data. However, the results from quantum mechanical calculations were not as good as the results of MDTI simulations owing to inappropriate treatment of intermolecular hydrogen bonds between the solute molecule and the first shell of solvent molecules in the Tomasi model of the solvent reaction field. The radial distribution functions between solute atoms and solvent atoms confirmed the formation of hydrogen bonds between the solute molecule and surrounding methanol molecules and indicated that the zwitterionic form is associated more with an organized solvent structure at the level of the first solvation shell than is the molecular form. Received: 26 April 2002 / Accepted: 9 September 2002 / Published online: 31 March 2003     

Bonding structure of phenylacetylene on hydrogen-terminated Si(111) and Si(100): surface photoelectron spectroscopy analysis and ab initio calculations

Interfaces between phenylacetylene (PA) monolayers and two silicon surfaces, Si(111) and Si(100), are probed by X-ray photoelectron spectroscopy (XPS), ultraviolet photoelectron spectroscopy (UPS), and near-edge X-ray absorption fine structure (NEXAFS) spectroscopy, and the results are analyzed using ab initio molecular orbital calculations. The monolayer systems are prepared via the surface hydrosilylation reaction between PA and hydrogen-terminated silicon surfaces. The following spectral features are obtained for both of the PA-Si(111) and PA-Si(100) systems: a broad π-π* shakeup peak at 292 eV (XPS), a broad first ionization peak at 3.8 eV (UPS), and a low-energy C 1s → π* resonance peak at 284.3 eV (NEXAFS). These findings are ascribed to a styrene-like π-conjugated molecular structure at the PA-Si interface by comparing the experimental data with theoretical analysis results. A conclusion is drawn that the vinyl group can keep its π-conjugation character on the hydrogen-terminated Si(100) [H:Si(100)] surface composed of the dihydride (SiH(2)) groups as well as on hydrogen-terminated Si(111) having the monohydride (SiH) group. The formation mechanism of the PA-Si(100) interface is investigated within cluster ab initio calculations, and the possible structure of the H:Si(100) surface is discussed based on available data.