Optimization of trial wave functions for use in quantum Monte Carlo with application to LiH

Methods for both partial and full optimization of wavefunction parameters are explored, and these are applied to the LiH molecule. A partial optimization can be easily performed with little difficulty. But to perform a full optimization we must avoid a wrong minimum, and deal with linear-dependency, time step-dependency and ensemble-dependency problems. Five basis sets are examined. The optimized wavefunction with a 3-function set gives a variational energy of -7.998 + 0.005 a.u., which is comparable to that (-7.990 + 0.003) 1 of Reynold's unoptimized \fin ( a double-~ set of eight functions). The optimized wavefunction with a double~ plus 3dz2 set gives ari energy of -8.052 + 0.003 a.u., which is comparable with the fixed-node energy (-8.059 + 0.004)1 of the \fin. The optimized double-~ function itself gives an energy of -8.049 + 0.002 a.u. Each number above was obtained on a Bourrghs 7900 mainframe computer with 14 -15 hrs CPU time.

Systematic construction of models for the exchange hole of density functional theory.

Dans ce travail, nous étendons le nombre de conditions physiques actuellement con- nues du trou d’échange exact avec la dérivation de l’expansion de quatrième ordre du trou d’échange sphérique moyenne exacte. Nous comparons les expansions de deux- ième et de quatrième ordre avec le trou d’échange exact pour des systèmes atomiques et moléculaires. Nous avons constaté que, en général, l’expansion du quatrième ordre reproduit plus fidèlement le trou d’échange exact pour les petites valeurs de la distance interélectronique. Nous démontrons que les ensembles de base de type gaussiennes ont une influence significative sur les termes de cette nouvelle condition, en étudiant com- ment les oscillations causées par ces ensembles de bases affectent son premier terme. Aussi, nous proposons quatre modèles de trous d’échange analytiques auxquels nous imposons toutes les conditions actuellement connues du trou d’échange exact et la nou- velle présentée dans ce travail. Nous évaluons la performance des modèles en calculant des énergies d’échange et ses contributions à des énergies d’atomisation. On constate que les oscillations causeés par les bases de type gaussiennes peuvent compromettre la précision et la solution des modèles.

Ab initio study of MF2 (M=Mn, Fe, Co, Ni) rutile-type compounds using the periodic unrestricted Hartree-Fock approach

The ab initio periodic unrestricted Hartree-Fock method has been applied in the investigation of the ground-state structural, electronic, and magnetic properties of the rutile-type compounds MF2 (M=Mn, Fe, Co, and Ni). All electron Gaussian basis sets have been used. The systems turn out to be large band-gap antiferromagnetic insulators; the optimized geometrical parameters are in good agreement with experiment. The calculated most stable electronic state shows an antiferromagnetic order in agreement with that resulting from neutron scattering experiments. The magnetic coupling constants between nearest-neighbor magnetic ions along the [001], [111], and [100] (or [010]) directions have been calculated using several supercells. The resulting ab initio magnetic coupling constants are reasonably satisfactory when compared with available experimental data. The importance of the Jahn-Teller effect in FeF2 and CoF2 is also discussed.

Computation of relativistic symmetry orbitals for finite double point groups

A program is presented for the construction of relativistic symmetry-adapted molecular basis functions. It is applicable to 36 finite double point groups. The algorithm, based on the projection operator method, automatically generates linearly independent basis sets. Time reversal invariance is included in the program, leading to additional selection rules in the non-relativistic limit.

Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

Effective atomic orbitals for fuzzy atoms

The method of extracting effective atomic orbitals and effective minimal basis sets from molecular wave function characterizing the state of an atom in a molecule is developed in the framework of the "fuzzy" atoms. In all cases studied, there were as many effective orbitals that have considerable occupation numbers as orbitals in the classical minimal basis. That is considered to be of high conceptual importance

CASPT2 Study of the Potential Energy Surface of the HSO(2) System

The importance of the HSO(2) system in atmospheric and combustion chemistry has motivated several works dedicated to the study of associated structures and chemical reactions. Nevertheless controversy still exists in connection with the reaction SH + O(2) -> H + SO(2) and also related to the role of the HSOO isomers in the potential energy surface (PES). Here we report high-level ab initio calculation for the electronic ground state of the HSO(2) system. Energetic, geometric, and frequency properties for the major stationary states of the PES are reported at the same level of calculations:,CASPT2/aug-cc-pV(T+d)Z. This study introduces three new stationary points (two saddle points and one minimum). These structures allow the connection of the skewed HSOOs and the HSO(2) minima defining new reaction paths for SH + O(2) -> H + SO(2) and SH + O(2) -> OH + SO. In addition, the location of the HSOO isomers in the reaction pathways have been clarified.

Electronic properties of a methane-water solution

Electronic properties of a methane-water solution were investigated by a sequential quantum mechanical/molecular dynamics approach. Upon hydration methane acquires an induced dipole moment of similar to 0.5 +/- 0.2 D. This is related to polarisation effects and to weak methane-water hydrogen bond interactions. From gas phase to solution, the first vertical excitation and ionisation energies of methane are red-shifted by 0.45 +/- 0.25 and 0.87 +/- 0.40 eV, respectively. We also report results for the dynamic polarisability of methane in water. In comparison with water, no difference was found for the average monomeric dipole moment of water molecules in close interaction with methane. (c) 2011 Elsevier B.V. All rights reserved.

Large Changes of Static Electric Properties Induced by Hydrogen Bonding: An ab Initio Study of Linear HCN Oligomers

We report the partitioning of the interaction-induced static electronic dipole (hyper)polarizabilities for linear hydrogen cyanide complexes into contributions arising from various interaction energy terms. We analyzed the nonadditivities of the studied properties and used these data to predict the electric properties of an infinite chain. The interaction-induced static electric dipole properties and their nonadditivities were analyzed using an approach based on numerical differentiation of the interaction energy components estimated in an external electric field. These were obtained using the hybrid variational-perturbational interaction energy decomposition scheme, augmented with coupled-cluster calculations, with singles, doubles, and noniterative triples. Our results indicate that the interaction-induced dipole moments and polarizabilities are primarily electrostatic in nature; however, the composition of the interaction hyperpolarizabilities is much more complex. The overlap effects substantially quench the contributions due to electrostatic interactions, and therefore, the major components are due to the induction and exchange induction terms, as well as the intramolecular electron-correlation corrections. A particularly intriguing observation is that the interaction first hyperpolarizability in the studied systems not only is much larger than the corresponding sum of monomer properties, but also has the opposite sign. We show that this effect can be viewed as a direct consequence of hydrogen-bonding interactions that lead to a decrease of the hyperpolarizability of the proton acceptor and an increase of the hyperpolarizability of the proton donor. In the case of the first hyperpolarizability, we also observed the largest nonadditivity of interaction properties (nearly 17%) which further enhances the effects of pairwise interactions.

Spectroscopic and structural studies of bis[isopropoxy(thiocarbonyl)]sulfide, [(CH(3))(2)CHOC(S)](2)S

Structural and conformational properties of the molecule bis[isopropoxy(thiocarbonyl)]sulfide, [(CH(3))(2)CHOC(S)](2)S, have been studied by vibrational spectroscopy (IR and Raman) and quantum chemical calculations (HF and B3LYP with 6-31+G* basis sets). The crystal and molecular structure of the title compound was determined by X-ray diffraction methods. It crystallizes in the monoclinic C2/c space group with a = 8.4007(4), b = 13.5936(5), c = 10.3648(5) angstrom, beta = 106.024(4)degrees and Z = 4 molecules per unit cell. The molecules are sited on a crystallographic twofold axis passing through the sulphide atom and arranged in layers perpendicular to the b-axis. The solid state IR and Raman spectra of the compound give no sign of any other rotamer. (C) 2009 Elsevier B.V. All rights reserved.

Vibrational properties, crystal X-ray diffraction structure and quantum chemical calculations on a divalent sulfur substituted phthalimide: 1H-Isoindole-1,3(2H)-dione, 2-[(methoxycarbonyl)thio]

Structural and conformational properties of 1H-Isoindole-1,3(2H)-dione, 2-[(methoxycarbonyl)thio] (S-phthalimido O-methyl thiocarbonate) are analyzed using a combined approach including X-ray diffraction, vibrational spectra and theoretical calculation methods. The vibrational properties have been studied by infrared and Raman spectroscopies along with quantum chemical calculations (B3LYP and B3PW91 functional in connection with the 6-311++G** and aug-cc-pVDZ basis sets). The crystal structure was determined by X-ray diffraction methods. The substance crystallizes in the monoclinic P2(1)/c space group with a = 6.795(1), b = 5.109(1), c = 30.011(3) angstrom, beta = 90.310(3)degrees and Z = 4 molecules per unit cell. The conformation adopted by the N-S-C=O group is syn (C=O double bond in synperiplanar orientation with respect to the N-S single bond). The experimental molecular structure is well reproduced by the MP2/aug-cc-pVDZ method. (C) 2010 Elsevier B.V. All rights reserved.

Cyclocreatine, an anticancer and neuroprotective agent. Spectroscopic, structural and theoretical study

The structural, spectroscopic and theoretical study of cyclocreatine (1-carboxymethyl-2-iminoimidazolidine, CyCre) has been performed prompted by the biological relevance of the molecule and its potential role as a ligand in biometalic compounds. The crystal structure of CyCre has been determined by X-ray diffraction methods. The compound crystallizes as a zwitterion in the monoclinic system, space group P2(1)/c. The crystal is further stabilized by a network of N-H center dot center dot center dot O bonds. Infrared and Raman spectra of the solid, electronic spectra of aqueous solutions at different pH values and (1)H and (13)C NMR spectra have been recorded and analyzed. Band assignments were accomplished with the help of theoretical calculations. Optimized molecular geometries, harmonic vibrational frequencies and molecular electrostatic potentials were calculated using methods based on the density functional theory. (C) 2010 Elsevier B.V. All rights reserved.

Electronic Structure of the Ground and Low-Lying Electronic States of MoB and MoB(+)

Multiconfigurational second-order perturbation theory (CASSCF//CASPT2) and quadruple-zeta ANO-RCC basis sets were employed to investigate the ground and low-lying electronic states of MoB and MoB(+). Spectroscopic constants, potential energy curves, wavefunctions, Mulliken population analyses, and ionization energies are given. The ground state of MoB is of X(6)Pi symmetry (R(e) = 1.968 angstrom, omega(e) = 664 cm(-1), and mu = 2.7 D), giving rise to a Omega = 7/2 ground state after including spin-orbit coupling. For MoB(+), the ground state is computed to be of X(7)Sigma(+) symmetry (R(e) = 2.224 angstrom, omega(e) = 141 cm(-1), and mu = 1.2 D), with an adiabatic ionization energy of 7.19 eV and a vertical one of 7.53 eV. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem 111: 3362-3370, 2011

Predicting and characterizing new molecular species at a high-level theoretical approach: HSeCl, HClSe, and SeCl

New molecular species HSeCl, HClSe, and SeCl were investigated at a high level of theory, CCSD(T), with a series of correlation consistent basis sets with extrapolation to the CBS limit. Account has been taken for valence-only and core-valence correlation effects, and of anharmonic effects on the vibrational frequencies. HSeCl is 43.25 kcal mol (1) more stable than HClSe. A barrier (Delta G(#)) of 47.20 kcal mol (1) separates these species. Internuclear distances are generally overestimated by 0.008 angstrom in the valence-only correlation calculations. Inclusion of anharmonicity leads to much improved vibrational frequencies. For SeCl, we estimate Delta H(f) (0 K) = 23.96 and Delta H(f) (298.15 K) = 24.64 kcal mol (1); for HSeCl, we had 4.20 and 4.97 kcal mol (1), respectively. (C) 2010 Elsevier B.V. All rights reserved.

Solvatochromic Shift of the pi -> pi* Transition in All-Trans, Cis-13, Cis-11, Cis-9, and Cis-7 Retinal Isomers Induced by Water and Methanol

The solvatochromic shift of the lowest singlet it pi -> pi* electronic transition in the all-trans, cis-13, cis-11, cis-9, and cis-7 retinal isomers were computed under the influence of water, methanol, and benzene solvents. Excitation energies were calculated in gas phase and in solution. The calculations in solution were performed considering the sequential Monte Carlo (MC) /Quantum Mechanical approach. The MC simulations were performed considering the full retinal isomer molecules and 900 water molecules, 900 methanol, or 400 benzene ones. The OPLS/AA parametrization was chosen for retinal, methanol, and benzene molecules and the SPC model was used for water one. From the MC calculations 100 independent configurations were selected, with 100 solvent molecules in thermodynamical equilibrium at T = 298.15 K. Average point-charges were obtained from those independent configurations for water, methanol, and benzene solvent. TDDFT and CASSCF//CASPT2 methodologies were used to compute the vertical excitation energy of the retinal isomers in different environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2076-2087, 2010