THE CRYSCOR PROJECT : STATUS AND PROSPECTS
|
|
- Kimberly Martin
- 5 years ago
- Views:
Transcription
1 THE CRYSCOR PROJECT : STATUS AND PROSPECTS ( # ) Cesare PISANI, Massimo BUSSO, Gabriella CAPECCHI, Silvia CASASSA, Roberto DOVESI, Lorenzo MASCHIO (*) Claudio ZICOVICH-WILSON ( ) Vic R. SAUNDERS ( & ) Martin SCHÜTZ ( # ) Università di Torino (Italy) (*) Universidad de Cuernavaca (Mexico) ( ) CRLC Laboratories, Daresbury (GB) ( & ) Universität Regensburg (Germany) Cesare Pisani - LCC2004 1
2 THE STRATEGY Combining two well assessed, robust, compatible technologies Geometrical and structural analysis of periodic system Accurate HF solution in terms of AOs Local representation of occupied manifold (WF) CRYSTAL On Local Correlation (LC) theory and techniques at various levels of sophistication Meyer-Pulay-Werner ( MOLPRO ) Cesare Pisani - LCC2004 2
3 Local-correlation approach for crystals MEYER, PULAY, SAEBØ,,WERNER, KNOWLES, HETZER, MANBY, SCHÜTZ, (Arkansas, Stuttgart, Birmingham, Bristol,...) Local, orthonormal representation of occupied HF manifold Local Molecular Orbitals Symmetry Adapted Wannier Functions (SAWF) Local, non-orthonormal, redundant representation of virtual manifold Projected Atomic Orbitals (PAO) (PAO) Reformulation of standard approaches ( MPn, CCSD(T), ) TRUNCATION STRATEGY Dynamical correlation effects are short-ranged: Ignore excitations from very distant pairs of SAWFs Exploit translational (and point) symmetry of SAWFs and PAOs N--SCALING n--scaling (n is the size of the irreducible part of the crystalline cell) Explore the role of computational parameters Cesare Pisani - LCC2004 3
4 THE OBJECTIVE To produce (in a reasonably short time) a public domain ab initio code for estimating electron correlation effects in non-conducting crystals, characterized by: generality simplicity robustness acceptable efficiency, and open to improvements As a first step, an MP2 code has been prepared and is here presented. Cesare Pisani - LCC2004 4
5 Is the objective worth the effort? Is the strategy appropriate? Cesare Pisani - LCC2004 5
6 Is the objective worth the effort? In our opinion it is. A number of ab-initio periodic correlation codes are ready, in preparation or in project (Stoll, Bartlett, Scuseria, Birkenheuer, Saunders-Orlando, Malrieu-Evangelisti, Monkhorst.). However, a conceptually simple, reliable, easily accessible code, could serve the purpose of providing reference data, exploring the role of computational parameters, checking the usefulness of alternative techniques, etc. Cesare Pisani - LCC2004 6
7 Is the strategy appropriate? With respect to other approaches in this area, CRYSCOR is characterized by its being strictly founded on the HF program CRYSTAL (language, basic techniques, etc.), so as to be fully compatible with it in a sense, CRYSCOR is the post-hf option of CRYSTAL. Disadvantages Advantages Cesare Pisani - LCC2004 7
8 Disadvantages CRYSTAL s computational technology is efficient but somewhat rigid and outdated. Inserting new parts of code may require a lot of work and attention. Cesare Pisani - LCC2004 8
9 Advantages Full space group symmetry of system and basis functions are fully and efficiently provided by CRYSTAL. Quasi-HF periodic solutions are obtained, and their quality can be easily assessed. Efficient and well tested techniques for generating localized and symmetry adapted WFs are available. Fourier transformation techniques are available (for instance, for calculating the inverse of quasi-diagonal translationally invariant matrices). Cesare Pisani - LCC2004 9
10 CRYSTAL 2003 : Main Features VR Saunders, R Dovesi, C Roetti, R Orlando, CM Zicovich-Wilson, NM Harrison, K Doll, B Civalleri, IJ Bush, Ph D Arco, M Llunell Distributed starting 9/03. Info : (Rev. Computational Chemistry, in press) The code FORTRAN90, fully parallelized, dynamic memory allocation The periodic model Consistent treatment of periodicity : 3D, 2D, 1D, 0D Ewald techniques for lattice sums (specific for 1D, 2D, 3D) Full exploitation of point symmetry in direct and reciprocal space Basis set Bloch functions as Linear Combination of Atomic Orbitals (contractions of Hermite Gaussian Functions) All-Electron or Valence-only-plus-Pseudopotential basis set Hamiltonians RHF, UHF Kohn-Sham techniques with Local and Gradient-corrected exchange and correlation functionals Hybrid DFT-HF exchange functionals Energy derivatives Automated geometry optimisation based on analytical gradient Wave function analysis and manipulation Band structure, PDOS, Charge, spin, electron momentum density, Structure factors, Compton profiles Elastic, dielectric constant, piezoelectric, hyperfine and nuclear quadrupole coupling tensors, [ Vibrational frequencies, based on analytical gradient ] Localized Wannier Functions [Symmetry adapted] Cesare Pisani - LCC
11 A pre-requisite for Local correlation methods in crystals: Efficient generation of Wannier functions (WF) to span the occupied HF space {Ψ n (k)} occ FT + localization criterion { w sg } occ Edmiston-Ruedenberg (1965) : Maximum intra-lo repulsion (N 5 ) Boys (1966) : Maximum distance between LO centroids (N 3 ) Pipek-Mezey (1989) : Minimum number of atoms per LO (N 3 ) What does efficient mean? Computationally inexpensive Well localized WFs Strictly orthonormal WFs Zicovich, Dovesi & Saunders, J. Chem. Phys. 115, 9708 (2001) Symmetry adapted WFs Zicovich, Casassa (2004) Cesare Pisani - LCC
12 CRYSCOR work to date and scheme of presentation 1. Generalization of the LOCALI part of CRYSTAL03, to produce symmetry adapted Wannier Functions (see Casassa-Zicovich poster) 2. Reformulation of LC-MP2 equations, so as to exploit translational and point symmetry 3. CRYSCOR code preparation (from CRYSTAL output to final results) (see Casassa et al. poster) 4. Refinement work on the integral part of the code to obtain 2-el integrals (i a j b ) either exactly or in a multipolar approximation (see Capecchi-Maschio poster) 5. Test of computational parameters (molecular cases + Diamond, Silicon, SiC, BN, BeS) Cesare Pisani - LCC
13 LC Reformulation of LMP2 and use of Symmetry R ij ab = 0 = K ij ab + Σ cd [ f ac T ij cd S db + S ac Tij cd f dc ] + Σ cd [ S ac Σ k (f ik T kj cd + Tik cd f kj )S cb ] E 2 = Σ (ij) Σ ab (ij) (Kij ab + Rij ab ) (2 Tij ab Tij ba ) ( a sum of pair energies) Ψ (1) = Σ (ij) Σ ( 2 ab (ij) Tij ab Tij ba ) 2 Φ ij ab Φij ba K ij ab = ( i a j b) = Σ µρνσ µρνσ c WF iµ cpao aρ cwf jν cpao bσ ( µ ρ ν σ) Symmetry exploitation Translational : The first WF index (i) is always confined in the zero cell Rotational : Only irreducible WF pairs (i j) need to be considered Irreducible pairs may have a residual symmetry which can be exploited in the K integral evaluation (Capecchi-Maschio poster). Cesare Pisani - LCC
14 Use of Symmetry in the Update step R ij ab = Σ cd [ S ac Σ k (f ik T kj cd + Tik cd f kj )S cb ] = Σ cd [ S ac β ij cd S cb ] T ij ab = Rij ab / ε εij ab To update amplitudes of irreducible WF pairs, we need those of all pairs. For each irreducible pair ij do Update T ij For each symmetry operator of i j (n) do Obtain amplitude of rotated pair T injn For each irreducible pair kl do If (k i n If (l j n and l (close-to) j n ) then β kl = β kl + T inj n f jn l and k (close-to) i n ) then β kl = β kl + f k in T injn enddo enddo enddo Cesare Pisani - LCC
15 BLOCK DIAGRAM OF CRYSCOR CRYSCOR (on disk) CRYSTAL MP2MAIN ORIQAO SYMPAIR DOMAINS MULTIPOLES MP2INT MP2CORE MP2LOOP Crystal type and symmetry; AO, S AOg,F AOg, WF (coefficients & symmetry) Computational parameters from input cards Construct PAOs, Calculate F g PAO, S PAO PAO g PAO,F g WF Recognize irreducible ww pairs and wp-wp bi-pairs Domains (Boughton Boughton-Pulay Pulay), Pair domains and their classification: Calculate multipolar expansion of WF-PAO products Calculate 2-el el-integrals ( exact or multipolar ) (see Capecchi-Maschio poster) Solve LMP2 equations (next dia) Cesare Pisani - LCC
16 SOLUTION OF LMP2 EQUATIONS MP2MAIN MP2CORE For each irreducible WW pair, do LON BSETR FWPRIM IGUESS Construct local orthonormal orbitals (lon lon) Organize quantities for iterative loop (close close-by pairs,, etc.) Calculate local pseudo-canonical virtual orbitals Provide initial guess for amplitudes T and for coupling coupling β quantities Enddo MP2LOOP For each irreducible WW pair, do End YES NO conv? ECALC UPDATE Read K, T, β, Update T,Update, E2 Update β Enddo Cesare Pisani - LCC
17 Main computational parameters in LMP2 for crystals R ij ab = Kij ab + Σ cd [ f ac T ij cd S db + S ac Tij cd f dc ] Σ cd [ S ac β ij cd S cb ] E 2 = Σ (ij) Σ ab (ij) (Kij ab + Rij ab ) (2 Tij ab Tij ba ) K ij ab = ( i a j b) = Σ µρνσ µρνσ c WF iµ cpao aρ cwf jν cpao bσ ( µ ρ ν σ) Basis set {µ} (representation of WFs and PAOs) Truncation of WF and PAO tails ( c WF iµ > tow ; cpao aρ > toq) Exact / multipolar treatment of K integrals (Capecchi-Maschio ) Prescreening of exact K integrals (Schwarz+density screening) Size of WF domains (range of occ. to virt. excitations in Σ ) ab (ij) Maximum WF-WF distance (range of Σ (ij) ) Cesare Pisani - LCC
18 C Shell exponent s coefficient p coefficient d coefficient 6-21G* 1s sp sp d Cesare Pisani - LCC
19 Si 6-21G* Shell exponent s coefficient p coefficient d coefficient s sp sp sp d Cesare Pisani - LCC
20 Main computational parameters in LMP2 for crystals R ij ab = Kij ab + Σ cd [ f ac T ij cd S db + S ac Tij cd f dc ] Σ cd [ S ac β ij cd S cb ] E 2 = Σ (ij) Σ ab (ij) (Kij ab + Rij ab ) (2 Tij ab Tij ba ) K ij ab = ( i a j b) = Σ µρνσ µρνσ c WF iµ cpao aρ cwf jν cpao bσ ( µ ρ ν σ) Basis set {µ} (representation of WFs and PAOs) Truncation of WF and PAO tails ( c WF iµ > tow ; cpao aρ > toq) Exact / multipolar treatment of K integrals (Capecchi-Maschio ) Prescreening of exact K integrals (Schwarz+density screening) Size of WF domains (range of occ. to virt. excitations in Σ ) ab (ij) Maximum WF-WF distance (range of Σ (ij) ) Cesare Pisani - LCC
21 PAOs and WFs (SiC and BeS) Cesare Pisani - LCC
22 Main computational parameters in LMP2 for crystals R ij ab = Kij ab + Σ cd [ f ac T ij cd S db + S ac Tij cd f dc ] Σ cd [ S ac β ij cd S cb ] E 2 = Σ (ij) Σ ab (ij) (Kij ab + Rij ab ) (2 Tij ab Tij ba ) K ij ab = ( i a j b) = Σ µρνσ µρνσ c WF iµ cpao aρ cwf jν cpao bσ ( µ ρ ν σ) Basis set {µ} (representation of WFs and PAOs) Truncation of WF and PAO tails ( c WF iµ > tow ; cpao aρ > toq) Exact / multipolar treatment of K integrals (Capecchi-Maschio ) Prescreening of exact K integrals (Schwarz+density screening) Size of WF domains (range of occ. to virt. excitations in Σ ) ab (ij) Maximum WF-WF distance (range of Σ (ij) ) Cesare Pisani - LCC
23 Influence of Computational parameters on Diamond results Basis set: 6-21G* ; valence-only MP2 calculations Standard parameter settings are in bold letters (outlined in yellow in tables) Parameters WF Size: There is only one type of WF (C-C bond localized orbital). The corresponding domain comprises 2 (C C) or 8 ( C 3 C C C 3 ) atoms Tow : threshold of c WF iµ coefficients ( 0.01, 0.005, 0.001, , 0.001) Toq : threshold of c PAO aρ coefficients ( 0.01, 0.005, 0.001, , 0.001) Pairs: n. of irreducible pairs considered [exact/multipolar] (4,19 / 0,154,301,407) The threshold for disregarding 2-el integrals according to Schwarz criterion has been set to 10-8, except in the cases indicated with an asterisk, where it is set to Results [bb] : contribution to -E 2 from [C C C C] pair (microhartree) [bcb ] : contribution to -E 2 from [C C CC ] (microhartree) -E 2 : absolute value of MP2 energy per cell (microhartree) Time: CPU time/10 3 sec on AMD Athlon, 1.6 GHz Cesare Pisani - LCC
24 Influence of Computational parameters on Diamond results Domain size and number of pairs All energies in microhartrees, times in 10 3 sec on AMD Athlon 1.6 GHz Size Tow Toq Pairs [bb] [bcb ] -E2 Time / / (+2363) 19/ / /0 19/ (+3861) Cesare Pisani - LCC
25 Influence of Computational parameters on Diamond results Truncation of tails and number of pairs Size Tow Toq Pairs [bb] [bcb ] -E 2 Time / / / / / / / / / / * / * / * / * Cesare Pisani - LCC
26 Comparison between 6-21G* results for diamond-like molecules (MOLPRO) and Diamond (CRYSCOR, standard setting) [bb] [bcb ] [b b ] [b b ] System -E2 4 C-C Neopentane 12 C-H C 5 H C-C Adamantane 16 C-H C 10 H C 35 H Diamond m=4 m=24 m=24 m= C All energies in µhartree; the best fit estimate of diamond E2 energy per cell from present molecular results is µhartree Cesare Pisani - LCC
27 Cesare Pisani - LCC
28 Comparison between 6-21G* results for diamond-like molecules (MOLPRO) and Diamond (CRYSCOR, standard setting) [bb] [bcb ] [b b ] [b b ] System -E2 Neopentane C 5 H Adamantane C 10 H C-C 36 C-H. C 35 H Diamond m=4 m=24 m=24 m= C All energies in µhartree; the best fit estimate of diamond E2 energy per cell from present molecular results is µhartree Cesare Pisani - LCC
29 Cesare Pisani - LCC
30 Comparison between 6-21G* results for diamond-like molecules (MOLPRO) and Diamond (CRYSCOR, standard setting) [bb] [bcb ] [b b ] [b b ] System -E2 4 C-C Neopentane 12 C-H C 5 H C-C Adamantane 16 C-H C 10 H C-C 36 C-H. C 35 H C-C Diamond m=4 m=24 m=24 m= per cell C All energies in µhartree; the best fit estimate of diamond E2 energy per cell from present molecular results is µhartree Cesare Pisani - LCC
31 Summary of MP2 valence-only results (6-21 G* basis set) BE HF E MP2 -E at MP2 BE TOT BE INCR BE EXP Diamond Silicon * SiC * BN * BeS * * Still not at convergence Cesare Pisani - LCC
32 STATE OF PROJECT AND PROSPECTS Current work : Refinement and standardization of basic program Higher efficiency, generalization of WF symmetry treatment Tests Basis set quality, different systems (MgO, Polymers, Rare gases ) Open problems and lines of development (Schütz, Usvjat, ): Efficient calculation and transformation of K integrals (Density fitting, etc.) Auxiliary basis set (to complement the HF set) Extension to other local correlation schemes, either within the MOLPRO framework (CCSD, MP4,.. ), or different (IEPA + NO ; MP2-R12 (Kutzelnigg, Klopper, Manby) ;?) Cesare Pisani - LCC
A local MP2 periodic study of crystalline argon
Journal of Physics: Conference Series A local MP2 periodic study of crystalline argon To cite this article: S Casassa et al 2008 J. Phys.: Conf. Ser. 117 012007 Recent citations - Laplace transformed MP2
More informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for CrystEngComm. This journal is The Royal Society of Chemistry 2014 Electronic Supplementary Information Configurational and energetical study of the (100) and
More informationMP2 Description of Solids: The CRYSCOR Project
MP2 Description of Solids: The CRYSCOR Project Lorenzo Maschio lorenzo.maschio@unito.it Dipartimento di Chimica and NIS centre, University of Torino, Italy The Regensburg group C. Pisani S. Casassa M.
More informationTranslation Symmetry, Space Groups, Bloch functions, Fermi energy
Translation Symmetry, Space Groups, Bloch functions, Fermi energy Roberto Orlando and Silvia Casassa Università degli Studi di Torino July 20, 2015 School Ab initio Modelling of Solids (UniTo) Symmetry
More informationAb initio structure prediction for molecules and solids
Ab initio structure prediction for molecules and solids Klaus Doll Max-Planck-Institute for Solid State Research Stuttgart Chemnitz, June/July 2010 Contents structure prediction: 1) global search on potential
More informationQuantum-chemical approach to cohesive properties of metallic beryllium
Quantum-chemical approach to cohesive properties of metallic beryllium Elena Voloshina 1, Beate Paulus 2, and Hermann Stoll 3 1 Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Straße 38, 01187
More informationAb initio treatment of electron correlations in polymers: Lithium hydride
JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 10 8 MARCH 2000 Ab initio treatment of electron correlations in polymers: Lithium hydride chain and beryllium hydride polymer Ayjamal Abdurahman a) Max-Planck-Institut
More informationCalculations of band structures
Chemistry and Physics at Albany Planning for the Future Calculations of band structures using wave-function based correlation methods Elke Pahl Centre of Theoretical Chemistry and Physics Institute of
More informationWavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts.
Wavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts. Silvia Casassa Università degli Studi di Torino July 12, 2017 Minnesota Workshop on ab initio MSC Symmetry
More informationWavefunction-based Correlation Calculations for Hole and Electron Capture States in Solids and Polymers
Wavefunction-based orrelation alculations for ole and Electron apture States in Solids and Polymers Uwe Birkenheuer 1, hrista Willnauer, Malte von Arnim, Walter Alsheimer, Dmitry Izotov Wavefunction-based
More informationTowards gas-phase accuracy for condensed phase problems
Towards gas-phase accuracy for condensed phase problems Fred Manby Centre for Computational Chemistry, School of Chemistry University of Bristol STC 2006: Quantum Chemistry Methods and Applications Erkner,
More informationarxiv:physics/ v1 [physics.comp-ph] 19 Dec 2005
Analytical Hartree-Fock gradients with respect to the cell parameter: systems periodic in one and two dimensions K. Doll Institut für Mathematische Physik, TU Braunschweig, Mendelssohnstraße 3, D-38106
More informationLinear Scaling Local Coupled Cluster Theory with Density Fitting. I: 4-External Integrals
Linear Scaling Local Coupled Cluster Theory with Density Fitting. I: 4-External Integrals Martin Schütz Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
More informationVibrational frequencies in solids: tools and tricks
Vibrational frequencies in solids: tools and tricks Roberto Dovesi Gruppo di Chimica Teorica Università di Torino Torino, 4-9 September 2016 This morning 3 lectures: R. Dovesi Generalities on vibrations
More informationHigh Accuracy Local Correlation Methods: Computer Aided Implementation
High Accuracy Local Correlation Methods: Computer Aided Implementation Marcel Nooijen Alexander Auer Princeton University University of Waterloo USA Canada So Hirata, PNNL Supported by: NSF ITR (Information
More informationThe orbital-specific-virtual local coupled cluster singles and doubles method
The orbital-specific-virtual local coupled cluster singles and doubles method Jun Yang, Garnet Kin-Lic Chan, Frederick R. Manby, Martin Schütz, and Hans-Joachim Werner Citation: J. Chem. Phys. 136, 144105
More informationASSESSMENT OF DFT METHODS FOR SOLIDS
MSSC2009 - Ab Initio Modeling in Solid State Chemistry ASSESSMENT OF DFT METHODS FOR SOLIDS Raffaella Demichelis Università di Torino Dipartimento di Chimica IFM 1 MSSC2009 - September, 10 th 2009 Table
More information7/29/2014. Electronic Structure. Electrons in Momentum Space. Electron Density Matrices FKF FKF. Ulrich Wedig
Electron Density Matrices Density matrices Γ, an alternative to the wavefunction Ψ, for the description of a quantum system Electronic Structure The N-particle density matrix Electrons in Momentum Space
More informationGaussian Basis Sets for Solid-State Calculations
Gaussian Basis Sets for Solid-State Calculations K. Doll Molpro Quantum Chemistry Software Institute of Theoretical Chemistry, D-70569 Stuttgart, Germany MW-MSSC 2017, Minneapolis, July 10, 2017 Introduction
More informationBasis sets for SIESTA. Emilio Artacho. Nanogune, Ikerbasque & DIPC, San Sebastian, Spain Cavendish Laboratory, University of Cambridge
Basis sets for SIESTA Emilio Artacho Nanogune, Ikerbasque & DIPC, San Sebastian, Spain Cavendish Laboratory, University of Cambridge Solving: Basis set Expand in terms of a finite set of basis functions
More informationInternational Journal of Quantum Chemistry
International Journal of Quantum Chemistry First-principles calculation of second-order elastic constants and equations of state for Lithium Azide, LiN, and Lead Azide, Pb(N ) Journal: International Journal
More informationAN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University
AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,
More informationCRYSTAL in parallel: replicated and distributed (MPP) data. Why parallel?
CRYSTAL in parallel: replicated and distributed (MPP) data Roberto Orlando Dipartimento di Chimica Università di Torino Via Pietro Giuria 5, 10125 Torino (Italy) roberto.orlando@unito.it 1 Why parallel?
More informationOrbital Density Dependent Functionals
Orbital Density Dependent Functionals S. Kluepfel1, P. Kluepfel1, Hildur Guðmundsdóttir1 and Hannes Jónsson1,2 1. Univ. of Iceland; 2. Aalto University Outline: Problems with GGA approximation (PBE, RPBE,...)
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationFast and accurate Coulomb calculation with Gaussian functions
Fast and accurate Coulomb calculation with Gaussian functions László Füsti-Molnár and Jing Kong Q-CHEM Inc., Pittsburgh, Pennysylvania 15213 THE JOURNAL OF CHEMICAL PHYSICS 122, 074108 2005 Received 8
More informationFully Automated Implementation of the Incremental Scheme: Application. to CCSD Energies for Hydrocarbons and Transition Metal Compounds
Fully Automated Implementation of the Incremental Scheme: Application to CCSD Energies for Hydrocarbons and Transition Metal Compounds Joachim Friedrich 1, Michael Hanrath 1, Michael Dolg 1 1 Institute
More informationOne-Electron Properties of Solids
One-Electron Properties of Solids Alessandro Erba Università di Torino alessandro.erba@unito.it most slides are courtesy of R. Orlando and B. Civalleri Energy vs Wave-function Energy vs Wave-function Density
More informationSupplementary material for Electronic Structure of IrO 2 : the Role of the Metal D Orbitals
Supplementary material for Electronic Structure of IrO 2 : the Role of the Metal D Orbitals Yuan Ping 1, Giulia Galli 2 and William A. Goddard III 3 1 Joint Center for Artificial Photosynthesis, Lawrence
More informationUsing BLIS for tensor computations in Q-Chem
Using BLIS for tensor computations in Q-Chem Evgeny Epifanovsky Q-Chem BLIS Retreat, September 19 20, 2016 Q-Chem is an integrated software suite for modeling the properties of molecular systems from first
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More informationComputational Chemistry I
Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of
More informationChapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set
Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless
More informationPractical Guide to Density Functional Theory (DFT)
Practical Guide to Density Functional Theory (DFT) Brad Malone, Sadas Shankar Quick recap of where we left off last time BD Malone, S Shankar Therefore there is a direct one-to-one correspondence between
More informationAb-initio Electronic Structure Calculations β and γ KNO 3 Energetic Materials
ISSN 0974-9373 Vol. 15 No.3 (2011) Journal of International Academy of Physical Sciences pp. 337-344 Ab-initio Electronic Structure Calculations of α, β and γ KNO 3 Energetic Materials Pradeep Jain and
More informationTDDFT in Chemistry and Biochemistry III
TDDFT in Chemistry and Biochemistry III Dmitrij Rappoport Department of Chemistry and Chemical Biology Harvard University TDDFT Winter School Benasque, January 2010 Dmitrij Rappoport (Harvard U.) TDDFT
More informationMolecular Simulation I
Molecular Simulation I Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences
More informationCohesive properties of CeN and LaN from first principles
APS/23-QED Cohesive properties of CeN and LaN from first principles Elena Voloshina a and Beate Paulus b a Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Straße 38, 087 Dresden, Germany b
More informationElectronic energy optimisation in ONETEP
Electronic energy optimisation in ONETEP Chris-Kriton Skylaris cks@soton.ac.uk 1 Outline 1. Kohn-Sham calculations Direct energy minimisation versus density mixing 2. ONETEP scheme: optimise both the density
More informationElectron Correlation - Methods beyond Hartree-Fock
Electron Correlation - Methods beyond Hartree-Fock how to approach chemical accuracy Alexander A. Auer Max-Planck-Institute for Chemical Energy Conversion, Mülheim September 4, 2014 MMER Summerschool 2014
More informationElectronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch
Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation There are 10 20 atoms in 1 mm 3
More informationThe orbital-specific virtual local triples correction: OSV-L(T)
The orbital-specific virtual local triples correction: OSV-L(T) Martin Schütz, Jun Yang, Garnet Kin-Lic Chan, Frederick R. Manby, and Hans-Joachim Werner Citation: The Journal of Chemical Physics 138,
More informationLocal Approaches to the Simulation of Electron Correlation in complex systems
Local Approaches to the Simulation of Electron Correlation in complex systems Martin Schütz Institut für Physikalische und Theoretische Chemie, Universität Regensburg Universitätsstraße 31, D-93040 Regensburg
More informationCHEM6085: Density Functional Theory
Lecture 11 CHEM6085: Density Functional Theory DFT for periodic crystalline solids C.-K. Skylaris 1 Electron in a one-dimensional periodic box (in atomic units) Schrödinger equation Energy eigenvalues
More informationarxiv:cond-mat/ v1 10 May 1996
Cohesive energies of cubic III-V semiconductors Beate Paulus, Peter Fulde Max-Planck-Institut für Physik komplexer Systeme, Bayreuther Str. 40, 01187 Dresden, Germany arxiv:cond-mat/9605064v1 10 May 1996
More informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
More informationElectronic Structure Methodology 1
Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E
More informationHandbook of Computational Quantum Chemistry. DAVID B. COOK The Department of Chemistry, University of Sheffield
Handbook of Computational Quantum Chemistry DAVID B. COOK The Department of Chemistry, University of Sheffield Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1998 CONTENTS 1 Mechanics and molecules 1 1.1
More informationChemistry 4560/5560 Molecular Modeling Fall 2014
Final Exam Name:. User s guide: 1. Read questions carefully and make sure you understand them before answering (if not, ask). 2. Answer only the question that is asked, not a different question. 3. Unless
More information6: Plane waves, unit cells, k- points and all that
The Nuts and Bolts of First-Principles Simulation 6: Plane waves, unit cells, k- points and all that Durham, 6th- 13th December 2001 CASTEP Developers Group with support from the ESF ψ k Network Overview
More informationThe high-pressure phase transitions of silicon and gallium nitride: a comparative study of Hartree Fock and density functional calculations
J. Phys.: Condens. Matter 8 (1996) 3993 4000. Printed in the UK The high-pressure phase transitions of silicon and gallium nitride: a comparative study of Hartree Fock and density functional calculations
More informationValence-band structure of group-iv semiconductors by means of local increments
PHYSICAL REVIEW B VOLUME 55, NUMBER 20 15 MAY 1997-II Valence-band structure of group-iv semiconductors by means of local increments Jürgen Gräfenstein* Max-Planck-Institut für Physik komplexer Systeme
More informationSame idea for polyatomics, keep track of identical atom e.g. NH 3 consider only valence electrons F(2s,2p) H(1s)
XIII 63 Polyatomic bonding -09 -mod, Notes (13) Engel 16-17 Balance: nuclear repulsion, positive e-n attraction, neg. united atom AO ε i applies to all bonding, just more nuclei repulsion biggest at low
More informationAll electron optimized effective potential method for solids
All electron optimized effective potential method for solids Institut für Theoretische Physik Freie Universität Berlin, Germany and Fritz Haber Institute of the Max Planck Society, Berlin, Germany. 22
More informationThe Plane-Wave Pseudopotential Method
Hands-on Workshop on Density Functional Theory and Beyond: Computational Materials Science for Real Materials Trieste, August 6-15, 2013 The Plane-Wave Pseudopotential Method Ralph Gebauer ICTP, Trieste
More informationOn the selection of domains and orbital pairs in local correlation treatments
On the selection of domains and orbital pairs in local correlation treatments Hans-Joachim Werner and Klaus Pflüger Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More informationThe Interpretation of the Short Range Disorder in the Fluorene- TCNE Crystal Structure
Int. J. Mol. Sci. 2004, 5, 93-100 International Journal of Molecular Sciences ISSN 1422-0067 2004 by MDPI www.mdpi.net/ijms/ The Interpretation of the Short Range Disorder in the Fluorene- TCNE Crystal
More informationproperties Michele Catti Dipartimento di Scienza dei Materiali Università di Milano Bicocca, Italy
Elastic and piezoelectric tensorial properties Michele Catti Dipartimento di Scienza dei Materiali Università di Milano Bicocca, Italy (catti@mater.unimib.it) 1 Tensorial physical properties of crystals
More informationCoupled-Cluster Theory. Nuclear Structure
Coupled-Cluster Theory! for Nuclear Structure!!!! Sven Binder INSTITUT FÜR KERNPHYSIK! 1 Nuclear Interactions from Chiral EFT NN 3N 4N NLO LO N 2 LO +... N 3 LO +... +... +... 2 Nuclear Interactions from
More informationOrbital dependent correlation potentials in ab initio density functional theory
Orbital dependent correlation potentials in ab initio density functional theory noniterative - one step - calculations Ireneusz Grabowski Institute of Physics Nicolaus Copernicus University Toruń, Poland
More informationQUANTUM CHEMISTRY FOR TRANSITION METALS
QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationPseudopotential. Meaning and role
Pseudopotential. Meaning and role Jean-Pierre Flament jean-pierre.flament@univ-lille.fr Laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM) Université de Lille-Sciences et technologies MSSC2018
More informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationBand calculations: Theory and Applications
Band calculations: Theory and Applications Lecture 2: Different approximations for the exchange-correlation correlation functional in DFT Local density approximation () Generalized gradient approximation
More informationSupporting Information. for
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2016 Supporting Information for Isoreticular Zirconium-Based Metal-Organic Frameworks:
More informationSoftware implementation of correlated quantum chemistry methods. Exploiting advanced programming tools and new computer architectures
Software implementation of correlated quantum chemistry methods. Exploiting advanced programming tools and new computer architectures Evgeny Epifanovsky Q-Chem Septermber 29, 2015 Acknowledgments Many
More informationFundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering
Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering Outline PART 1: Fundamentals of Density functional theory (DFT)
More informationDensity Functional Theory - II part
Density Functional Theory - II part antonino.polimeno@unipd.it Overview From theory to practice Implementation Functionals Local functionals Gradient Others From theory to practice From now on, if not
More informationAB INITIO MODELING OF ALKALI METAL CHALCOGENIDES USING SOGGA THEORY
Int. J. Chem. Sci.: 13(4), 215, 163-1638 ISSN 972-768X www.sadgurupublications.com AB INITIO MODELING OF ALALI METAL CHALCOGENIDES USING SOGGA THEORY HITESH CHANDRA SWARNAR and GUNJAN ARORA a,* Department
More informationElectron Correlation Methods
Electron Correlation Methods HF method: electron-electron interaction is replaced by an average interaction E HF c = E 0 E HF E 0 exact ground state energy E HF HF energy for a given basis set HF E c
More informationGround State Projector QMC in the valence-bond basis
Quantum Monte Carlo Methods at Work for Novel Phases of Matter Trieste, Italy, Jan 23 - Feb 3, 2012 Ground State Projector QMC in the valence-bond basis Anders. Sandvik, Boston University Outline: The
More informationarxiv:cond-mat/ v1 17 May 1995
Projection of plane-wave calculations into atomic orbitals Daniel Sanchez-Portal, Emilio Artacho, and Jose M. Soler Instituto de Ciencia de Materiales Nicolás Cabrera and Departamento de Física de la Materia
More informationGeometry optimization of solids
The Minnesota Workshop on ab Initio Modeling in Solid State Chemistry with CRYSTAL Minneapolis, MN(U.S.A.) 9-14 July 2017 7 1 Geometry optimization of solids Bartolomeo Civalleri Dip. di Chimica IFM, Via
More informationWave function methods for the electronic Schrödinger equation
Wave function methods for the electronic Schrödinger equation Zürich 2008 DFG Reseach Center Matheon: Mathematics in Key Technologies A7: Numerical Discretization Methods in Quantum Chemistry DFG Priority
More informationNANOSTRUCTURED OXIDES: NEW MATERIALS FOR ENERGY AND ENVIRONMENT
NANOSTRUCTURED OXIDES: NEW MATERIALS FOR ENERGY AND ENVIRONMENT Quantum Chemistry Laboratory Dipartimento di Scienza dei Materiali Università Milano-Bicocca http://www.mater.unimib.it/utenti/pacchioni
More informationSpin densities and related quantities in paramagnetic defects
Spin densities and related quantities in paramagnetic defects Roberto Orlando Dipartimento di Scienze e Tecnologie Avanzate Università del Piemonte Orientale Via G. Bellini 25/G, Alessandria roberto.orlando@unipmn.it
More informationNWChem: Coupled Cluster Method (Tensor Contraction Engine)
NWChem: Coupled Cluster Method (ensor Contraction Engine) What we want to solve H Ψ = E Ψ Many Particle Systems Molecular/Atomic Physics, Quantum Chemistry (electronic Schrödinger equations) Solid State
More informationThe case of orbitals Trond Saue (Toulouse/Oslo) Trond Saue (LCPQ, Toulouse) The case of orbitals ESQC / 137
The case of orbitals Trond Saue (Toulouse/Oslo) Trond Saue (LCPQ, Toulouse) The case of orbitals ESQC 2017 1 / 137 Some facts The lanthanide series: Orbital energies Numerical 4-component atomic Hartree-Fock
More informationMD simulation: output
Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square
More informationThe calculation of the universal density functional by Lieb maximization
The calculation of the universal density functional by Lieb maximization Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry,
More informationHighly accurate quantum-chemical calculations
1 Highly accurate quantum-chemical calculations T. Helgaker Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, Norway A. C. Hennum and T. Ruden, University
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry Chemicum 4th floor vesa.hanninen@helsinki.fi September 10, 2013 Lecture 3. Electron correlation methods September
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationThe Schrödinger equation for many-electron systems
The Schrödinger equation for many-electron systems Ĥ!( x,, x ) = E!( x,, x ) 1 N 1 1 Z 1 Ĥ = " $ # " $ + $ 2 r 2 A j j A, j RAj i, j < i a linear differential equation in 4N variables (atomic units) (3
More informationBeyond Hartree-Fock: MP2 and Coupled-Cluster Methods for Large Systems
John von Neumann Institute for Computing Beyond Hartree-Fock: MP2 and Coupled-Cluster Methods for Large Systems Christof Hättig published in Computational Nanoscience: Do It Yourself, J. Grotendorst, S.
More informationLecture 9. Hartree Fock Method and Koopman s Theorem
Lecture 9 Hartree Fock Method and Koopman s Theorem Ψ(N) is approximated as a single slater determinant Φ of N orthogonal One electron spin-orbitals. One electron orbital φ i = φ i (r) χ i (σ) χ i (σ)
More informationHandbook of Computational Quantum Chemistry
Handbook of Computational Quantum Chemistry David B. Cook Dept. of Chemistry University of Sheffield DOVER PUBLICATIONS, INC. Mineola, New York F Contents 1 Mechanics and molecules 1 1.1 1.2 1.3 1.4 1.5
More informationA theoretical study of stability, electronic, and optical properties of GeC and SnC
JOURNAL OF APPLIED PHYSICS VOLUME 88, NUMBER 11 1 DECEMBER 2000 A theoretical study of stability, electronic, and optical properties of GeC and SnC Ravindra Pandey a) Department of Physics, Michigan Technological
More informationChemistry, Physics and the Born- Oppenheimer Approximation
Chemistry, Physics and the Born- Oppenheimer Approximation Hendrik J. Monkhorst Quantum Theory Project University of Florida Gainesville, FL 32611-8435 Outline 1. Structure of Matter 2. Shell Models 3.
More informationComputational chemistry with GAMESS: a very brief overview with examples
Computational chemistry with GAMESS: a very brief overview with examples PHY-6120 Molecular Physics (Spring 2015), UConn Phys. Dept. Feb 17 th 2015 H = ħ2 2μ i Intro: V(R) for diatomic molecules + k Z
More informationLecture 4: Hartree-Fock Theory
Lecture 4: Hartree-Fock Theory One determinant to rule them all, One determinant to find them, One determinant to bring them all and in the darkness bind them Second quantization rehearsal The formalism
More informationIntroduction to Computational Chemistry: Theory
Introduction to Computational Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC andrew.gilbert@anu.edu.au 3023 Course Lectures Introduction Hartree Fock Theory Basis Sets Lecture 1 1 Introduction
More information4 Post-Hartree Fock Methods: MPn and Configuration Interaction
4 Post-Hartree Fock Methods: MPn and Configuration Interaction In the limit of a complete basis, the Hartree-Fock (HF) energy in the complete basis set limit (ECBS HF ) yields an upper boundary to the
More informationHomologation of Boronic Esters with Organolithium Compounds: A Computational Assessment of Mechanism
Homologation of Boronic Esters with Organolithium Compounds: A Computational Assessment of Mechanism Stéphanie Essafi,*,1 Simone Tomasi, 2 Varinder K. Aggarwal, 1 Jeremy Harvey*,1 1 School of Chemistry,
More informationLarge Scale Electronic Structure Calculations
Large Scale Electronic Structure Calculations Jürg Hutter University of Zurich 8. September, 2008 / Speedup08 CP2K Program System GNU General Public License Community Developers Platform on "Berlios" (cp2k.berlios.de)
More informationNumerical construction of Wannier functions
July 12, 2017 Internship Tutor: A. Levitt School Tutor: É. Cancès 1/27 Introduction Context: Describe electrical properties of crystals (insulator, conductor, semi-conductor). Applications in electronics
More informationIntroduction to Computational Quantum Chemistry: Theory
Introduction to Computational Quantum Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC 3108 Course Lectures 2007 Introduction Hartree Fock Theory Configuration Interaction Lectures 1 Introduction
More information