Coupled-Cluster study of ‘no-pair’ bonding in the tetrahedral Cu4 cluster

Abstract Ab initio Coupled-Cluster calculations with single and double excitations and perturbative correction to the triple, CCSD(T), have been carried out for the high-spin electronic state, ( 5 A 2 ) , of the copper cluster Cu 4 in its tetrahedral arrangement. Like alkali metals clusters, tetrahedral Cu 4 presents a bound quintet state, i.e., a situation where all the valence electrons are unpaired. This rather exotic wavefunction, also known as no-pair bonding state, is examined in detail. The influence of the basis set is also analyzed, as well as the importance of the core correlation and the effect of the basis-set superposition errors.

Low-lying Rydberg states of HCl.

Vertical excitation energies belonging to some different Rydberg series of hydrogen chloride have been determined with a coupled-cluster theoretical approach. These excitation energies have allowed us to calculate electric dipole transition intensities in HCl and allow additional assessment of the calculation approach presently used to provide an adequate description of the valence and Rydberg states of HCl. The molecular quantum defect orbital has been applied to the calculation of oscillator strengths. In particular, new insight is given on the assignment of states as the R1Pi, the 1Delta(4dpi and 5ppi), the 1Sigma+(4dpi), and the nddelta(1Pi, 1Phi) and 4f states.

Electronic structure of the ground and excited states of beta-carboline.

Coupled-cluster calculations are used to compute the energy of conversion between the neutral and the zwitterionic forms of beta-carboline. The stability of the different species is discussed in terms of charge separation and aromatic character, which is related to magnetic criteria. By means of a linear response formalism the vertical excitation energies and oscillator strengths of the lowest singlet states of both structures as well as of the cationic species are determined. General agreement of the relative position and intensity of the different peaks with experimental data is achieved, but the overall spectra are slightly displaced because of solvent effects.

Elementary presentation of self‐consistent intermediate Hamiltonians and proposal of two totally dressed singles and doubles configuration interaction methods

Intermediate Hamiltonians are effective Hamiltonians which are defined on an N‐dimensional model space but which only provide n<N exact eigenvalues and the projections of the corresponding eigenvectors onto the model space. For a single root research, the intermediate Hamiltonian may be obtained from the restriction of the Hamiltonian to the model space by an appropriate, uniquely defined dressing of the diagonal energies or of the first column. Approximate self‐consistent dressings may be proposed. The simplest perturbative form gives the same result as the original 2nd order intermediate Hamiltonian or the ‘‘shifted Bk’’ technique but it is of easier implementation. Self‐consistent inclus…

Grand Canonical Monte Carlo Simulations to Determine the Optimal Interlayer Distance of a Graphene Slit-Shaped Pore for Adsorption of Methane, Hydrogen and their Equimolar Mixture.

The adsorption—for separation, storage and transportation—of methane, hydrogen and their mixture is important for a sustainable energy consumption in present-day society. Graphene derivatives have proven to be very promising for such an application, yet for a good design a better understanding of the optimal pore size is needed. In this work, grand canonical Monte Carlo simulations, employing Improved Lennard–Jones potentials, are performed to determine the ideal interlayer distance for a slit-shaped graphene pore in a large pressure range. A detailed study of the adsorption behavior of methane, hydrogen and their equimolar mixture in different sizes of graphene pores is obtained through ca…

Improved AMYR program: an algorithm for the theoretical simulation of molecular associations, including geometrical and topological characterization of the dimers.

Program AMYR, originally written by S. Fraga (University of Alberta, Canada), allows for the calculation of molecular associations using a pair-wise atom-atom potential. The interaction energy is evaluated through a 1/R expansion. Our improved version includes a dispersion energy term in the potential corrected by damping functions, the possibility of carrying out energy minimizations through variable metric methods, as well as the optional calculation of geometrical and topological indices. Program AMYR has been adapted also for high-performance computing and vectorization. An interactive version of the program carries out real-time molecular graphics showing simultaneously the energy prof…

Full configuration interaction calculation of the low lying valence and Rydberg states of BeH

The all-electron full configuration interaction (FCI) vertical excitation energies for some low lying valence and Rydberg excited states of BeH are presented in this article. A basis set of valence atomic natural orbitals has been augmented with a series of Rydberg orbitals that have been generated as centered onto the Be atom. The resulting basis set can be described as 4s2p1d/2s1p (Be/H) + 4s4p3d. It allows to calculate Rydberg states up to n= {3,4,5} of the s, p, and d series of Rydberg states. The FCI vertical ionization potential for the same basis set and geometry amounts to 8.298 eV. Other properties such as FCI electric dipole and quadrupole moments and FCI transition dipole and qua…

Vertical spectrum of ethene: uncontracted versus contracted correlation methods and the role of the adapted molecular orbitals

Abstract Ethylene vertical excitation energies (VEE) below 9.5 eV have been calculated with the complete active-space singles and doubles configuration interaction self-consistent size-consistent dressing (SC) 2 CAS–SDCI. The mean of the absolute deviation from experiment for the calculated VEE are 0.14, 0.05 and 0.03 eV for three different molecular orbital (MO) sets tried. The results show that an uncontracted method and a realistic MO set avoid the difficulties present in the calculation of the strongly mixed Rydberg-valence 1B 1u states.

Application of the open-shell size-consistent self-consistent singles and doubles configuration interaction method to multielectronic transitions in photoelectron spectra

The size-consistent self-consistent matrix dressing method has been applied on an open-shell single-configuration reference state. Once the reference state is converged, several low-lying roots can be obtained for the dressed configuration interaction (CI) matrices of appropriate symmetry. The CI matrices were built with a complete-active-space singles and doubles CI method in order to deal properly with multiconfiguration excited states. The vertical ionization and ionization–excitation transitions are obtained from the difference to the closed shell ground-state energy of the neutral molecule. The method has been applied to NH+3 and N+2 using atomic natural orbital basis sets and state-av…

Pair potential calculation of molecular associations: a vectorized version

Abstract The program AMYRVF is a vectorized and largely modified version of a previous program called AMYR for calculating molecular associations by means of Fraga's pair-wise atom-atom potential. Three new minimization procedures have been implemented as well as other improvements such as the inclusion of new pair-wise dispersion energy terms with damping functions, and the calculation of topological indices. Benchmark tests have been carried out on an IBM 3090 150E VF; the timings for the new vector algorithms and for the standard scalar computations, as well as the dependence of the overall performance gain on the size of interacting systems are reported for the VS Fortran 2.4 compiler e…

Ab initio study of rotational isomerism and electronic structure of isomeric bipyrroles

Abstract Ab initio calculations using STO-3G and 4-31G basis sets have been performed on the internal rotation barriers and conformational stabilities for 2,3′- and 3,3′-bipyrrole. The twofold rotation potential predicted for both isomers at minimal basis level becomes a more involved fourfold potential when the split-valence basis set is employed, because it takes into account more properly the nonbonded interannular interactions. A transoid-gauche minimum is predicted to have the minimal absolute conformational energy in both isomers. The electronic structure of the highest occupied MOs of 2,2′-, 2,3′- and 3,3′-bipyrrole are analyzed in terms of the single pyrrole MO pattern and a similar…

Totally dressed SDCI calculations: An application to HF and F2

A previously proposed procedure including the linked and unlinked contributions due to Triple and Quadruple excitations into a size-consistent SDCI-like model has been applied to HF and F2 single-bond systems. The procedure is a non-iterative approximation to the more general total dressing model, which is based on the intermediate Hamiltonians theory. Three basis sets have been employed: the correlation consistent cc-pVTZ basis, a similar one including 3d1f polarization functions, and another including one set of g polarization functions. Excellent agreement with experiment and high-quality calculations is obtained for both equilibrium distances and spectroscopic constants. The possibiliti…

Ab Initio Study on the Mechanism of Tropospheric Reactions of the Nitrate Radical with Haloalkenes:  Chloroethene

Ab Initio Study on the Mechanism of the Reactions of the Nitrate Radical with Haloalkenes:  1,2-Dichloroethene, 1,1-Dichloroethene, Trichloroethene, and Tetrachloroethene

A general mechanism for the reactions of the NO3 radical with 1,2-dichloroethene, 1,1-dichloroethene, trichloroethene, and tetrachloroethene is proposed from ab initio DFT calculations. The calculated mechanism shows three main parallel reaction pathways. For the systems where the two carbon atoms are differently substituted, the study includes both the attacks with Markownikoff and contra-Markownikoff orientation. The first reaction pathway leads to the formation of an epoxide along with the NO2 radical, the second one to the formation of carbonyl compounds, and the third one leads, through the cleavage of the C−C bond, to the formation of carbonyl compounds with a lower number of carbon a…

Ab initio calculations of the transfer parameters and coulombic repulsion and estimation of their effects on the electron delocalization and magnetic coupling in mixed-valence Keggin polyoxotungstates

International audience; In this work, we present ab initio calculations on embedded fragments that permit to extract the value of the effective electron transfer integral and coulombic repulsion between W nearest neighbour atoms in a mixed-valence αPW12O40 Keggin polyoxoanion. This allows us to perform a quantitative study of the influence of these two parameters on the magnetic properties of Keggin polyoxoanions reduced by two electrons. We surprisingly find that the electron transfer between edge-sharing and corner-sharing WO6 octahedra have very close values, and show that the punctual charges estimation of coulombic repulsion may not be accurate enough to study the electronic distributi…

A computational study of some electric and magnetic properties of gaseous BF3 and BCl3

We present the results of an extended computational study of the electric and magnetic properties connected to Cotton-Mouton birefringences, on the trifluoro- and trichloroborides in the gas phase. The electric dipole polarizabilities, magnetizabilities, quadrupole moments, and higher-order hypersusceptibilities—expressed as quadratic and cubic frequency-dependent response functions—are computed within Hartree-Fock, density-functional, and coupled-cluster response theories employing singly and doubly augmented correlation-consistent basis sets and London orbitals in the magnetic property calculations. The results, which illustrate the capability of time-dependent density-functional theory f…

Size-consistent single-reference methods for electronic correlation: a unified formulation through intermediate hamiltonian theory

Using the intermediate hamiltonian theory as a unique conceptual frame and the technique of CI matrix dressing, a wide series of single-reference methods for the treatment of the ground state correlation are reviewed, compared, and sometimes improved. These methods range from independent excitation approximation (the very next step beyond MP2) to coupled cluster, going through the so-called electron pair approximations and the (SC)2CI formalism. A hierarchy of these methods can be established according to two criteria: These formulations in terms of diagonalizations of dressed CI matrices avoid convergence problems, but their main advantage is their flexibility, since they apply to multi-re…

Size-consistent self-consistent configuration interaction from a complete active space : Excited states

The self-consistent size consistent on a complete active space singly and doubly configuration interaction (SC)2CAS-SDCI method is applied to excited states. The (SC)2 correction is performed on a closed shell state, and the excited states are obtained by diagonalization of the dressed matrix. A theoretical justification of the transferability of the improvement concerning the dressing state to all roots of the matrix is presented. The method is tested by three tests on the spectrum of small molecules. sanchezm@uv.es ; nebot@uv.es

Interacting induced dipoles polarization model for molecular polarizabilities. Reference molecules, amino acids and model peptides

Abstract We outline a method for the calculation of molecular dipole ( μ ) and quadrupole ( Θ = ) moments and dipole–dipole polarizabilities ( α = ) which we have successfully applied to a series of reference molecules, amino acids and model peptides. The results for μ are in line with CPHF reference calculations. In particular, the calculated positive value of CO is in agreement with both experimental and CI calculations. The computation of ( α = ) has been performed by the interacting induced dipoles polarization model that calculates tensor effective anisotropic point polarizabilities (method of Applequist et al.). The POLAR program cannot be used as a black box. Some tests should be per…

Potential models for the simulation of methane adsorption on graphene: development and CCSD(T) benchmarks

Different force fields for the graphene–CH4 system are proposed including pseudo-atom and full atomistic models. Furthermore, different charge schemes are tested to evaluate the electrostatic interaction for the CH4 dimer. The interaction parameters are optimized by fitting to interaction energies at the DFT level, which were themselves benchmarked against CCSD(T) calculations. The potentials obtained with both the pseudo-atom and full atomistic approaches describe accurately enough the average interaction in the methane dimer as well as in the graphene–methane system. Moreover, the atom–atom potentials also correctly provide the energies associated with different orientations of the molecu…

Theoretical study of the low‐lying states of trans‐1,3‐butadiene

We present extensive ab initio calculations on the low‐lying electronic states of trans‐1,3‐butadiene within the multireference configuration interaction (MRCI) framework by selecting the configurations with a perturbative criterion. The X 1Ag ground state and 1 3Bu, 1 3Ag, 2 1Ag, and 1 1Bu valence excited states have been calculated at a fixed geometry. The results obtained are in good agreement with previous experimental and calculated values, and could help to understand polyene spectroscopy, photochemistry, and photophysics. The advantages of a MRCI method where the most important contributions to the total MRCI wave function, perturbatively selected, are treated variationally, and the …

Are most of the stationary points in a molecular association minima? Application of Fraga's potential to benzene-benzene

The importance of characterizing the stationary points of the intermolecular potential by means of Hessian eigenvalues is illustrated for the calculation of the benzene–benzene interaction using an atom-to-atom pair potential proposed by Fraga (FAAP). Two models, the standard one-center-per atom and another using three-centers-per atom due to Hunter and Sanders, are used to evaluate the electrostatic contributions and the results are compared. It is found in both cases that although using low-gradient thresholds allows optimization procedures to avoid many stationary points that are not true minima computing time considerations makes the usual procedure of using high-gradient thresholds [sa…

Full configuration interaction calculation of BeH adiabatic states.

An all-electron full configuration interaction (FCI) calculation of the adiabatic potential energy curves of some of the lower states of BeH molecule is presented. A moderately large ANO basis set of atomic natural orbitals (ANO) augmented with Rydberg functions has been used in order to describe the valence and Rydberg states and their interactions. The Rydberg set of ANOs has been placed on the Be at all bond distances. So, the basis set can be described as 4s3p2d1f3s2p1d(BeH)+4s4p2d(Be). The dipole moments of several states and transition dipole strengths from the ground state are also reported as a function of the R(Be-H) distance. The position and the number of states involved in sever…

Assessment for the mean value total dressing method: Comparison with coupled cluster including triples methods for BF, NO+, CN+, C2, BeO, NH3, CH2, H2O, BH, HF, SiH2, Li2, LiNa, LiBe+, NeH+, and O3

Limited previous experience with the mean value total dressing (MVTD) method had shown that MVTD energies for closed shell systems are generally better than CCSD(T) ones compared to FCI. The method, previously published as total dressing 2′(td-2′), is based on the single reference intermediate Hamiltonian theory. It is not a CC method but deals in a great part with the same physical effects that CC methods that incorporate amplitudes of triples such as CCSDT or its CCSDT-1n approaches. A number of test calculations comparing to diverse CC methods, as well as FCI and experiment when available, have been performed. The tests concern equilibrium energies in NH3 and CH2, equilibrium energies an…

Vertical spectrum of the C2H 2+ system. An open shell (SC)2-CAS-SDCI study

The open shell (SC)(2)-CAS-SDCI method along with a basis set of atomic natural orbitals (ANO) has been applied for calculating the main ionization potentials of acetylene, as well as the manifold of excited states of the different symmetries up to 32 eV. In this method, the single and double excitations of a CAS space are generated and the corresponding CI matrix is corrected by means of the (SC)(2) procedure that cancels the size-extensivity error and adds some high order contributions. The mean absolute error for the outer-valence X (2)Pi(u)(1pi(u) (-1)), A (2)Sigma(g) (+)(3sigma(g) (-1)), and B (2)Sigma(u) (+)(2sigma(u) (-1)) states, and the inner-valence C (2)Sigma(g) (+)(2sigma(g) (-1…

Interacting induced dipoles polarization model for molecular polarizabilities: application to benzothiazole (A)-benzobisthiazole (B) oligomers AB13A

Abstract We have outlined a method for the calculation of atomic net charges and molecular dipole-dipole,  α , dipole-quadrupole,  A , and quadrupole-quadrupole,  C , polarizabilities which we have successfully applied to benzothiazole (A)-benzobisthiazole (B) oligomers in the series AB1A to AB13A, and the results have been extrapolated to give predictions for polybenzobisthiazole. The calculation of  α ,  A and  C has been carried out by the interacting induced dipoles polarization model that calculates tensor effective anisotropic point polarizabilities (method of Applequist). The mean effective atomic αeffpolarizabilities for the monomeric unit are of the same order of magnitud…

Excitation energies and photoabsorption oscillator strengths of the Rydberg series in CF3Cl. A linear response and quantum defect study.

Vertical excitation energies of the CF(3)Cl molecule have been obtained from a response function approach with a CC reference function to determine absolute photoabsorption oscillator strengths in the molecular-adapted quantum defect orbital formalism (MQDO). The present work covers more highly excited Rydberg states than have been experimentally reported. Assessing of the reliability of the present calculations is provided through a comparative analysis between the results of the molecule and the Cl atom. This can be used to allow for predictions of the same type of properties in other analogous systems.

Ab Initio Study on the Mechanism of Tropospheric Reactions of the Nitrate Radical with Alkenes:  Propene

DFT Theoretical Study on the Reaction Mechanism of the Nitrate Radical with Alkenes:  2-Butene, Isobutene, 2-Methyl-2-butene, and 2,3-Dimethyl-2-butene

A general mechanism for the reactions of the NO3 radical with 2-butene, isobutene, 2-methyl-2-butene, and 2,3-dimethyl-2-butene is proposed on the basis of density functional theory (DFT) calculati...

Size-consistent ab initio calculation of the electric quadrupole moment of Cl2

Abstract The molecular electric quadrupole moment ( Θ ) of Cl 2 has been calculated using SDCI, and (SC) 2 -SDCI wave functions as well as CCSD, CCSD(T), and CC3 methods. All these correlation methods are single reference. All of them, but SDCI, are free of the size-extensivity error. The variation of Θ from the separated atoms to the equilibrium region is reported. The present results leads to an estimated value of 2.3520 a.u. (10.55 × 10 −40 Cm 2 ) corresponding to a CC(3) calculation at the CBS approach and including the ro-vibrational and thermal averaging corrections. This value is compatible with two experimental values and points to one of them as slightly more reliable.

Characterizing cavities in model inclusion molecules: a comparative study

We have selected fullerene-60 and -70 cavities as model systems in order to test several methods for characterizing inclusion molecules. The methods are based on different technical foundations such as a square and triangular tessellation of the molecule taken as a unitary sphere, spherical tessellation of the molecular surface, numerical integration of the atomic volumes and surfaces, triangular tessellation of the molecular surface, and a cubic lattice approach to a molecular space. Accurate measures of the molecular volume and surface area have been performed with the pseudo-random Monte Carlo (MCVS) and uniform Monte Carlo (UMCVS) methods. These calculations serve as a reference for the…

Full Configuration-Interaction Study on the Tetrahedral Li4 Cluster

International audience; The Li4 cluster low lying electronic states were studied. In particular we investigated the tetrahedral geometry at full CI and coupled cluster level, with basis sets of increasing quality. The 5A2 electronic state, characterized by having all the valence electrons unpaired, forming a quite stable no-pair bonding state, was studied in greater detail. In order to compare the energies we also studied the Li4 rhombus singlet ground state. The ability of coupled cluster with perturbative triples to correctly reproduce energy levels in a quasi-degenerate system was validated with respect to the full CI.

Electrically Conductive Phthalocyanine Assemblies. Structural and Non-Integer Oxidation Number Considerations

Aggregation is a well-known phenomenon in phthalocyanine chemistry. Interactions can occur between adjacent phthalocyanine rings, both in organic and aqueous phases, resulting in coupling between the electronic states of two, or more, phthalocyanine units [1].

Torsional effects on the molecular polarizabilities of the benzothiazole (A)-benzobisthiazole (B) oligomer A-B13-A

Abstract We outline a method for the calculation of multipole moments and molecular dipole-dipole ( ), dipole-quadrupole ( ), and quadrupole-quadrupole ( ) polarizabilities, which we have successfully applied to benzothiazole (A)-benzobisthiazole (B) oligomer A-B13-A. Three model rotational isomers have been characterized: (1) the fully planar (000) rotational isomer; (2) a conformation with each unit rotated 10° in the alternate direction (+−+), and (3) a rotational isomer with each unit rotated 10° in the same direction (+++). The dipole moment, μ , is smaller for isomers 000 and +−+ than for isomer +++. The calculation of , , and has been performed by use of the interacting induced dipol…

Modeling the Interaction of Carbon Monoxide with Flexible Graphene: From Coupled Cluster Calculations to Molecular-Dynamics Simulations

The interaction of CO with graphene was studied at different theoretical levels. Quantum-mechanical calculations on finite graphene models with the use of coronene for coupled cluster calculations and circumcoronene for B97D calculations showed that there was no preferential site for adsorption and that the most important factor was the orientation of CO relative to graphene. The parallel orientation was preferred, with binding energies around 9 kJ mol-1 at the CCSD(T) and B97D levels, which was in good agreement with experimental findings. From a large number of CO-circumcoronene and CO-CO interactions, computed at different distances and randomly generated orientations, parameters were fi…

Lower Rydberg series of methane: a combined coupled cluster linear response and molecular quantum defect orbital calculation.

Vertical excitation energies as well as related absolute photoabsorption oscillator strength data are very scarce in the literature for methane. In this study, we have characterized the three existing series of low-lying Rydberg states of CH4 by computing coupled cluster linear response (CCLR) vertical excitation energies together with oscillator strengths in the molecular-adapted quantum defect orbital formalism from a distorted Cs geometry selected on the basis of outer valence green function calculations. The present work provides a wide range of data of excitation energies and absolute oscillator strengths which correspond to the Rydberg series converging to the three lower ionization p…

Electron delocalization in mixed-valence Keggin polyoxometalates. Ab initio calculation of the local effective transfer integrals and its consequences on the spin coupling.

We present a quantitative evaluation of the influence of the electron transfer on the magnetic properties of mixed-valence polyoxometalates reduced by two electrons. For that purpose, we extract from valence-spectroscopy ab initio calculations on embedded fragments the value of the transfer integrals between W nearest-neighbor atoms in a mixed-valence alphaPW(12)O(40) polyoxowolframate Keggin anion. In contradiction with what is usually assumed, we show that the electron transfer between edge-sharing and corner-sharing WO(6) octahedra have very close values. Considering fragments of various ranges, we analyze the accuracy of calculations on fragments based on only two WO(5) pyramids which s…

Reducing CAS-SDCI space. Using selected spaces in configuration interaction calculations in an efficient way

A new method is presented, which allows an important reduction of the size of some Configuration Interaction (CI) matrices. Starting from a Complete Active Space (CAS), the numerous configurations that have a small weight in the CAS wave function are eliminated. When excited configurations (e.g., singly and doubly excited) are added to the reference space, the resulting MR-SDCI space is reduced in the same proportion as compared with the full CAS–SDCI. A set of active orbitals is chosen, but some selection of the most relevant excitations is performed because not all the possible excitations act as SDCI generators. Thanks to a new addressing technique, the computational time is drastically …

Study of a medium-size biological molecular association by means of a pair potential semiempirical approach: β-carboline-lumiflavin

The molecular association between some substituted β-carbolines and a model flavin, lumiflavin, has been studied by means of a semiempirical approach proposed by Fraga based on a 1/R expansion of atom-atom pair potentials. Only stacked minima have been considered because of their possible biological interest. The calculations characterize 15 different minimum stacked conformations, most of them occurring in the complexes considered. The structural and energetic effects of 1–Me, 6–OH and 7–OH substituents of β-carboline are discussed. An equilibrium conformation involving geometric overlap between the pteridinic portion of flavin and the indole group of β-carboline is predicted to be the mos…

Multi-scale theoretical investigation of molecular hydrogen adsorption over graphene: coronene as a case study

The physisorption of molecular hydrogen onto coronene is studied using a multi-scale theoretical approach with Density Functional Theory (DFT) calculations and Molecular Dynamics (MD) simulations. We consider two different kinds of model conformation for the approach of hydrogen towards the coronene i.e., systematic and random. For the systematic attack of hydrogen over coronene, the resulting potential energy profiles from DFT analysis are further found to resemble the Morse potential, and even the highly flexible Murrell–Sorbie (M–S) potential. The resulting M–S fitting also shows a zero-point energy correction of ∼16–17%. On the other hand, the potential energies from the random approach…

A theoretical study of the rotational structure of the ϵ(0,0) band of NO

This study has been focused on the ϵ(0,0) band of the nitric oxide molecule, associated with the absorption electronic transition D2Σ+X 2Π, in the energetic vacuum ultraviolet region. A temperature of 295 K has been considered. The Molecular Quantum Defect Orbital (MQDO) methodology, with which reliable spectroscopic data have been reached in the γ(0,0), δ(0,0), and δ(1,0) bands of the same molecule, has also been used for these calculations. We hope that the present results might be of straightforward use in atmospheric and interstellar chemistry. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010

Nitrogen Gas on Graphene: Pairwise Interaction Potentials

We investigate different types of potential parameters for the graphene-nitrogen interaction. Interaction energies calculated at DFT level are fitted with the semi-emperical Improved Lennard-Jones potential. Both a pseudo-atom potential and a full atomistic potential are considered. Furthermore, we consider the influence of the electrostatic part on the parameters using different charge schemes found in the literature as well as optimizing the charges ourselves. We have obtained parameters for both the nitrogen dimer and the graphene-nitrogen system. For the former, the four-charges Cracknell scheme reproduces with high precision the CCSD(T) interaction energy as well as the experimental di…

Full configuration interaction calculation of Be3.

The full configuration interaction (FCI) study of the ground state of the neutral beryllium trimer has been performed using an atomic natural orbitals [3s2p1d] basis set. Both triangular and linear structures have been considered for the Be(3) cluster. The optimal geometry for the equilateral triangle has been calculated. The potential energy cut sections along the normal a(1)(') mode and one of the components of the e(') mode have then been studied. The FCI symmetric atomization potential of the linear cluster is also reported. It shows a secondary van der Waals minimum at a long bond distance. All singular points in the potential energy curves are characterized. Other properties, like dis…

Molecular electric quadrupole moments calculated with matrix dressed SDCI

Abstract We have calculated the molecular electric quadrupole moment (MEQM) for the set of molecules N 2 , C 2 H 2 , CO, CO 2 , CS 2 , HF, and BH. We have used SR-SDCI and (SC) 2 -SR-SDCI methods and we have compared our results with high-level theoretical ones, including FCI values for HF and BH, and with experimental values. The calculated MEQM provides a test of the effect that the energy converged (SC) 2 dressing method brings to the SDCI wavefunctions. The results suggest that the (SC) 2 -SR-SDCI method can be a cost-effective and quite accurate method for the calculation of post-SCF effects on electric quadrupole moments.

Code Interoperability and Standard Data Formats in Quantum Chemistry and Quantum Dynamics: The Q5/Q5cost Data Model

Code interoperability and the search for domain-specific standard data formats represent critical issues in many areas of computational science. The advent of novel computing infrastructures such as computational grids and clouds make these issues even more urgent. The design and implementation of a common data format for quantum chemistry (QC) and quantum dynamics (QD) computer programs is discussed with reference to the research performed in the course of two Collaboration in Science and Technology Actions. The specific data models adopted, Q5Cost and D5Cost, are shown to work for a number of interoperating codes, regardless of the type and amount of information (small or large datasets) …

High-spin states in tetrahedral X4 clusters (X = H, Li, Na, K)

The high-spin electronic states for lithium, sodium, and potassium four-atom clusters were studied. In particular, we performed coupled cluster geometry optimization of the quintet state in tetrahedral geometry. The quintet state of these systems is characterized by having all the valence electron Unpaired, giving rise to the so-called no-pair bonding. Single-point full configuration interaction computations on the equilibrium geometries for the various Clusters are also presented. The analysis of the valence orbitals in a localized representation confirms the importance of the p atomic orbitals to explain this unusual type of bond. (C) 2009 Wiley Periodicals, Inc. Int J Quantum Chem 110: 8…

Rydberg, valence and mixed states in the vertical spectrum of HF

Abstract Size-consistent self-consistent multireference singles and doubles configuration interaction (SC)2–MR–SDCI, and linear response coupled cluster of singles and doubles including triples LR–CCSDR(3) methods, with a basis set of ANO’s augmented with a single series of molecular Rydberg functions, have been applied to the calculation of vertical excitation energies of HF in the 10–16 eV region. Special care has been put in the description of the valence, Rydberg and mixed states. Some advantage has been taken from the different physical contents of the methods to discuss some assignments. The (4,5) f Rydberg states of HF are predicted at 15.2 and 15.5 eV.

Theoretical study of the discrete and continuum spectrum of BeH

The transition intensities supplied in this Letter are directly connected with a description of the discrete and continuum spectrum of BeH. An attempt to meet our goal requires the calculation of the absorption oscillator strengths of several transitions to Rydberg states of BeH, together with differential oscillator strengths which give rise to different dipole-allowed photoionization channels from the molecular ground state. The calculations have been performed with the molecular-adapted quantum defect orbital (MQDO) approach. Predictions of new spectroscopic data on BeH at energies where high Rydberg transitions can take place, including the continuum region of the spectrum have been mad…

Self‐consistent intermediate Hamiltonians : A coupled cluster type formulation of the singles and doubles configuration interaction matrix dressing

This paper presents a new self‐consistent dressing of a singles and doubles configuration interaction matrix which insures size‐consistency, separability into closed‐shell subsystems if localized molecular orbitals (MOs) are used, and which includes all fourth order corrections. This method yields, among several schemes, a reformulation of the coupled cluster method, including fully the cluster operators of single and double excitations, and partially those of the triples (Bartlett’s algorithm named CCSDT‐1a). Further improvement can be easily included by adding exclusion principle violating corrections. Since it leads to a matrix diagonalization, the method behaves correctly in case of nea…

Vectorized ‘‘TOPO’’ Program for the Theoretical Simulation of Molecular Shape

The TOPO program for theoretical simulation of molecular shape is presented here. Molecular shape is characterized by a set of electrostatic and geometrical descriptors and topological indices including the fractal dimension of the accessible surface. An atom-atom analysis of all descriptors has been implemented. The program has been adapted for the realization of high-performance computing. Results are reported for porphine and phthalocyanine molecules.

Theoretical Study on the Structures and Electronic Spectra of TCNE2−

Investigations into the charge-separated states and electron-transfer transitions in tetracyanoethylene (TCNE) compiles have recently generated much interest. In this work we present theoretical calculations showing that the most stable structure of the dianion TCNE 2- has D 2d symmetry in vacuum as well as in the solvents dichloromethane and cicetonitrile. By means of the coupled cluster linear response, we compute the vertical electronic spectrum in both the gas phase and solution. The theoreitcal results, are compared to the experimental data and good agreement is achieved.

Theoretical ab initio study of the vinylcarbene-CuO complex: application to last step of the propylene partial oxidation mechanism on Cu2O

Abstract Using non-empirical pseudopotentials, we have performed a theoretical study on the last step of the propylene partial oxidation mechanism, the acrolein desorption from the catalyst surface. To represent the catalyst we have used a simplified model of a catalyst active site. We have study the C 3 H 4 CuO complex, analysing the dissociation process of the oxygenated complex from a metallic centre. The results indicate that this complex dissociates into cis - and trans -acrolein and a metal site, in a process without energetic barrier.

Universal model for the calculation of all organic solvent–water partition coefficients

Abstract We present the basis for building a universal organic solvation model to calculate solubility in any organic solvent and in water, as well as the organic solvent–water partition coefficient ( P ). Log P values are of the same order of magnitude as reference calculations but for a few cases which are discussed. Normalized log P contributions are sensitive to the rest of the atoms. When comparing porphin with phthalocyanine, the latter results in an amphipathic molecule. For C 70 , the contribution of a–e carbons to log P correlates with the distances from the nearest pentagon. The method has been also applied to benzobisthiazole oligomers and phenyl alcohols.

Conformational behaviour of 2,2?-bipyrrole

The rotational potential around the interannular bond in 2,2′-bipyrrole has been calculated making use of standard minimal STO-3G and split valence 4-31G basis sets. Geometrical optimization concerning the most significant interannular internal parameters has been performed with both basis sets. The trans conformer is predicted to be more stable than the cis. The minimal basis set predicts the existence of a cisoid-gauche minimum which after limited optimization becomes very shallow and it seems to be an artifact of the rigid rotor approximation. At 4-31G level, both the trans and cis conformers represent maxima in the potential curve and two gauche minima appear at θ=46.0° and θ=147.6°, th…

Flexibility in the Graphene Sheet: The Influence on Gas Adsorption from Molecular Dynamics Studies

Despite being considered completely rigid in most studies, graphene is really flexible leading to out-of-plane movements. In this work, the influence of such flexibility on the adsorption of methane and nitrogen on graphene is studied using molecular dynamics. Indeed, we have used intramolecular force fields for graphene with in-plane and out-of-plane components that allow for describing the movements and deformations of the graphene sheets and providing a more realistic description of the adsorbent. In addition, intermolecular force fields validated at the CCSD(T) level are used. We show that considering the movement of graphene in the adsorption study significantly improves the performanc…

The vertical spectrum of H2CO revisited: (SC)2-CI and CC calculations

The vertical electronic spectrum of formaldehyde has been studied by means of (SC)2-MR-SDCI and CCLR methods. Two basis sets of atomic natural orbitals (ANOs) complemented with a one-centre series of Rydberg orbitals were used. The first was taken from the CASPT2 study by Merchan, M., and Roos, B. O., 1995, Theoret. Chim. Acta, 92, 221, and may be described as C,O[4s3pld]/H[2slp] with a lslpld Rydberg series centred in the charge centroid of the 2B2 state of the cation. The second was a larger basis set that may be described as C,O[6s5p3d2f]/H[4s3p2d] + 3s3p3d in the same centre. The (SC)2 dressing may be applied efficiently to an MR-SDCI method and comparison with the dressed CAS-SDCI is s…

The problem of interoperability: A common data format for quantum chemistry codes

A common format for quantum chemistry (QC), enhancing code interoperability and communication between different programs, has been designed and implemented. An XML-based format, QC-ML, is presented for representing quantities such as geometry, basis set, and so on, while an HDF5-based format is presented for the storage of large binary data files. Some preliminary applications that use the format have been implemented and are also described. This activity was carried out within the COST in Chemistry D23 project “MetaChem,” in the Working Group “A meta-laboratory for code integration in ab initio methods.” © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007

Full configuration interaction calculation of singlet excited states of Be3

The full configuration interaction (FCI) study of the singlets vertical spectrum of the neutral beryllium trimer has been performed using atomic natural orbitals [3s2p1d] basis set. The FCI triangular equilibrium structure of the ground state has been used to calculate the FCI vertical excitation energies up to 4.8 eV. The FCI vertical ionization potential for the same geometry and basis set amounts to 7.6292 eV. The FCI dipole and quadrupole transition moments from the ground state are reported as well. The FCI electric quadrupole moment of the X (3)A(1) (') ground state has been also calculated with the same basis set (Theta(zz)=-2.6461 a.u., Theta(xx)=Theta(yy)=-1/2Theta(zz)). Twelve of …

Multistate active spaces from local CAS-SCF molecular orbitals: the photodissociation of HFCO as an example.

A recently developed algorithm to generate localized molecular orbitals (LMO) is applied to the study of excited states along a photodissociation process. The LMOs allow for the selection of a consistent complete active space (CAS) for the simultaneous description of all the electronic states involved in a multistate process on the basis of simple chemical criteria. The local nature of the orbitals is used to label them in a unique way that does not depend on the molecular geometry. The selection of the electronic configurations of interest for the set of target states on only the basis of the dominant excitations required by the simplest configuration interaction (CI) descriptions for both…

AMYR 2: A new version of a computer program for pair potential calculation of molecular associations

AMYR is a computer program for the calculation of molecular associations using Fraga's pairwise atom-atom potential. The interaction energy is evaluated through a 1R expansion. The electrostatic energy is calculated through either the one-centre-per atom or the three-centres-per atom model by Hunter and Sanders. A pairwise dispersion energy term is included in the potential and corrected by a damping function. The program carries out energy minimizations through variable metric methods. The new version allows for the stationary point analysis of the intermolecular potential by means of the Hessian eigenvalues. Although using low-gradient thresholds optimization procedures to avoid many stat…

Bending Carbon Nanoforms for Supramolecular Recognition: A Topological Study on Hemifullerene-Based Aggregates

International audience; Buckybowls have risen as appealing fullerene fragment derivatives. Their intrinsic curvature has been exploited in the generation of host–guest supramolecular assemblies, not only through concave–convex complementarity but also through less-known concave–concave staggered arrangements. Whereas the stabilization of bowl-in-bowl dispositions has been ascribed to efficient π–π forces together with favorable dipole–dipole interactions, a detailed analysis on the forces guiding the formation of the staggered arrangements is missing so far. Herein, we present a thorough theoretical characterization of bowl-in-bowl vs staggered hemifullerene-based homodimers and heterodimer…

A study of coronene?coronene association using atom?atom pair potentials

A study of the coronene—coronene association using different interaction potentials based on an atom-atom pair potential proposed by Fraga has been performed. The interaction potentials employed differ in the way the electrostatic and/or dispersion contributions are computed. The influence of both contributions on the geometries predicted for the coronene dimer is discussed in order to analyze the effectiveness of the different interaction potentials. The stationary points found in each interaction energy hypersurface are characterized by calculating the Hessian eigenvalues. Results are discussed in the light of those previously reported for the benzene dimer. Stacked-displaced structures a…