Ab-initio simulation of liquid water by quantum Monte Carlo
|
|
- Mervyn Gardner
- 5 years ago
- Views:
Transcription
1 Ab-initio simulation of liquid water by quantum Monte Carlo Sandro Sorella G. Mazzola & Y. Luo SISSA, IOM DEMOCRITOS, Trieste A. Zen, L. Guidoni U. of L Aquila, L Aquila 28 July 2014, Mike Towler Institute, Apuan Alps
2 Outline/Motivations à The correlated wave function for realistic systems: from Hubbard to Hydrogen and Water à Few examples on accuracy and achievements on small molecules H 2 (H 2 O) 2 à Large number of electrons/long simulations with QMC now possible. à MD for realistic liquid à example on Hydrogen at high pressures à Can we do also liquid water?
3 We assume that HTc superconductivity shows up in a correlated systems due to strong correlation. The paradigm wave function is the Gutzwiller partially projected BCS (Mean Field) wavefunction: ψ VMC = exp( g n R n R )P N MF R where MF is the ground state of the BCS hamiltonian H = k,σ Variational Gutzwiller ansatz k [ 2(cosk x + cosk y ) µ VMC ]c + k,σ c kσ + Δ(cosk x cosk y ) % P N BCS = AGP = ' f k c k & k + + c k c + + k c k There are only 3 variational parameters for f k. ( * ) N /2 0 + h.c.
4 Important to optimize Jastrow and BCS toghether Hubbard Model : BCS order parameter H + + = t ci σ c jσ U ni ni < ij>, σ i J BCS> BCS> Hubbard model U/t= Energy minimization Iterations In mean field (BCS) no way to have BCS>0 for U>0 Theorem Lieb 90 Qualitative new features appear if Jastrow and BCS optimized toghether: RVB insulator or supercond.
5 Generalization of the wave function to continuous ψ T = J MF MF> may be a standard Slater determinant J is the so called Jastrow correlation term: " J = exp g( r i, % $ r j )' # $ i< j &' g is a generic function of two el. cooredinates The peculiarity of our approach (TurboRVB) is to fully optimize MF> and J in a localized basis of simple atomic orbitals (e.g. Gaussians 1s,2p ).
6 Generalization to realityà the Hamiltonian is: H = i Δ i 2 Z j r i R + ij j i< j 1 r i r + Z iz j j R i R j {R i } are atomic classical coordinates within the Born-Oppenheimer approximation r 1, r 2, r 3 r " N J SD = exp g( r i, % $ r j )' # $ i< j &' Det "# ψ i( r j )% & Given that, one can apply Variational Monte Carlo and compute all correlation functions by a statistical method No further approximation required fully ab-initio (no U, no double counting) But how to parametrize the function g?
7 The Gutzwiller for realistic systems u lr (r) = 1 2 g( r, r!) = u lr ( r r! )+ λ a,b i, j ψ a i ( r)ψ b j (! a,b,,i, j e.g. ψ a k ( r) = exp" Z k r R # a 2 $ a,b %, i.e. localized atomic orbitals, #λ kl # atoms variational parameters (say~1000) determined by: min a λkl, SD r ) r 1+ Br "a(b)" labels atom positions R a (R b ) The non-homogeneous part a=b is local like Gutzw. and useful to decrease # parameters (no 4-body a b) SD JHJ SD SD J 2 SD
8 Quantum Monte Carlo vs DFT, is it worth? In H 2 clear With a very small basis (2 gaussians/atom) one gets the essentially exact dispersion for H 2
9 The main question we want to address: What happens when we apply large pressure to a hydrogen molecular liquid? When the average distance between molecules is comparable with their bond length (~1.4 a.u.) we have a transition to a system where the molecule is no longer defined (atomic). According to band theory, from an insulator 2el/unit (H 2 ) to an half-filled band 1el./unit (H)à Metal, Wigner and Heterington prediction 35.
10 Example: DFT failure for hydrogen
11 As we have learned by Car and Parrinello (1998) Phase Diagram of realistic systems à Ab-Initio Molecular dynamics with Born-Oppenheimer approx. Evaluation of Forces are required within QMC Algorithmic differentiation helped much and thanks also to Tapenade (automatic diff.): SS & Luca Capriotti, JCP 133, (2010)
12 Cpu time referenced to simple VMC (only energy) for computing all 3M force components in water. (CPU VMC+forces)/(CPU VMC) AD Finite difference # Water molecules x 4 Use of pseudopotentials straightforward
13 Just to clarify a bit what we mean by AD Just a black box programming discipline allowing to compute all derivatives of e L ( r 1, r 2,, r N, R 1, R 2, R M ) With respect to all 3N electron coordinates And the 3M ionic coordinates at the same ~cost of computing the local energy e L And the same for the wave function.
14 Dynamics: 1. Efficient QMC forces Check: Newtonian dynamics of a H2 molecule. Verlet integrator.
15 Dynamics: 2. Generalized Langevin We have efficient but still noisy forces. We use Langevin dynamics in order to sample the canonical ensamble for the ions. Fluctuation - dissipation Covariance matrix of the forces Now the noise does not prevent the possibility of doing MD. Only renormalize the friction!
16 Dynamics: 2. Generalized Langevin We have efficient but still noisy forces. We use Langevin dynamics in order to sample the canonical ensamble for the ions.
17 At large N the first order is evident from the g(r) 6 g (r) r s = 1.28 r s = 1.24 r s = 1.23 r s = 1.22 r s = 1.21 r s = Jump T= 600 K r r s = 1.19
18 Liquid-liquid transition T= 2300 K
19 Our quantum Monte Carlo phase diagram ( ) for the first time with N=256 Hydrogen, is now much different from DFT!!!! Weir (molecular) Experiment 1996 We find that the molecular fluid is unexpectedly stable and the transition towards a fully atomic liquid occurs at much higher pressures.
20 Comparison small basis/good basis 54 atoms rs=1.44 T=1000K increased accuracyà more stable molecular!!!
21 Why is so different? QMC vs DFT PES for a single molecule inside the bulk (54H)
22 How to distinguish a metal from insulator? One can compute the density matrix: D(r,r')= ψ i i (r)ψ i (r') ψ i (r) are the optimized molecular orbitals Metalà Fermi surface à D(R,r) ~ R-r -2 Insulatorà Gap à D(R,r) ~exp(- R-r /ξ) Thus DM = 1 # atoms dr 3 D(R a, r) R a Metal DM à Insulator DM à Finite
23 There is some interesting crossover at T=2400K But we do not have enough large size for the MIT
24 And now few slides on liquid water Why water liquid simulation? It is fundamental in biological life, e.g. life with No waterà Non sense Phase diagram of immense difficulties, low energy and competing (e.g. Hydrogen bond and vdw long range) scales and still many things to understand by computer simulations. DFT problems, g(r) overstructured, eq. density large (20% off), supercooled liquid (melting at ~400K)
25 Reduction of number of parameters In QMC optimization the number of parameters is proportional to the dimension of the basis. It is useful to reduce them by hybrid contraction ϕ a,l,m,n (r) = lmn c a ψ GTO (r) lmn lmn,a This is not useful instead in chemistry as the molecular HF basis is used instead. See A. Zen, Y. Luo, SS and L.Guidoni JCTC 2013
26 Water dimer test: NB in dynamics we are interested in relative forces, i.e. derivative of binding energy 2.0 Water dimer bond energy Bond energy [kcal/mol] VMC-JSD small basis VMC-JSD-converged basis LRDMC Relevant region Distance oxygens [Ang]
27 Choice of Great freedom properly damping fast modes large time step reducing the correlation time Our choice γ estimated by QMC contains info of Hessian α α = α I QMC = Noise is useful!!! 0 cov( γ + Δ f ) 1 = α( R) 2T 0 α QMC
28 Water dimer integration test: Δ 0 = 8a.u. 300K DFT has not the covariance of forces and is much less efficient (smaller time steps) than QMC (fs)
29 Following the BO 150K fitting 1000K fitting up to V 2 up to V 4 Y. Luo, A. Zen and S. Sorella submitted to JCP
30 E pur si muove (and yet it moves) (Galileo Galilei) Second order Langevin dynamics of 32 water molecules at 300K with Variational Monte Carlo G. Mazzola, A. Zen, Y. Luo, L. Guidoni, and S. Sorella in preparation (2014)
31 Some info on this simulation 2048 nodes on BG/Q, 30 days simulation: 24milion core hours Time step = 1.54fs, Total time ~ 10ps Each step MDà 10 Optimization steps ~12000 Variational parameters ~ Sampling measurements/step Each sample after 1024 Metropolis step A huge computation, impossible without HPC
32 A different phylosophy is to use DMC to correct DFT (Blyp-2) Alfe et al JCP 2013.
33 But the quantum effects should play a role from J. A. Morrone and R.Car PRL, 2008 i) Peak positions are not changed by quantum ii) The radial distribution is substantially broaden
34 QMC water radial distribution function More than 4000 iterations~ 7ps Remaining differences: Ionic quantum effects? More accuracy? Size effects? A. Zen, G. Mazzola,, Y. Luo, L. Guidoni, and S. Sorella in preparation (2014)
35 64 waters possible with nodes (> cores!!!) 5 first peak in g OO (r) convergence r MAX region of fit VMC 32-wat g OO (r) VMC 64-wat r [Ang] / # steps Neutron diffraction - Soper, Chem. Phys. 258,121 (2000) Neutron diffraction - Soper, ISRN Phys. Chem.,2013 (2013) X-ray diffraction - LB Skinner et al., JCP 138, (2013) 32 waters VMC-based NVT MD simulation 32 waters DFT/BLYP-based NVT MD simulation 64 waters DFT/BLYP-based NVT MD simulation
36 The pressure problem: At the right density the average pressure is 0.47Gpa>> 1atm(=10-5 Gpa!!!) Pressure (Gpa) VMC LRDMC DMC reduces much (and change sign) this value a 2
37 Conclusions Realistic simulation of liquids are now possible also within fully many-body wave function based approach. 256H (64 H 2 O) are not so far from what is currently done within DFT ~500H (128 H 2 O) à Peak positions of the rdf are finally reproduced: no other ab-initio first principle simulation is able à Accuracy of VMC probably not enough, and also quantum effects should play a role. à Several applications are now possible, allowing to falsify or improve DFT predictions. Liquid water is currently under investigation.
38 TurboRVB Quantum Monte Carlo package main developers: Sandro Sorella Michele Casula Claudio Attaccalite A. Zen Ye Luo G. Mazzola E. Coccia Funding: MIUR-CNR, Riken-Tokio, PRACE Computational support CINECA Bologna, SISSA
Ab-initio molecular dynamics for High pressure Hydrogen
Ab-initio molecular dynamics for High pressure Hydrogen Claudio Attaccalite Institut d'electronique, Microélectronique et Nanotechnologie (IEMN), Lille Outline A brief introduction to Quantum Monte Carlo
More informationQMC dissociation energy of the water dimer: Time step errors and backflow calculations
QMC dissociation energy of the water dimer: Time step errors and backflow calculations Idoia G. de Gurtubay and Richard J. Needs TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay.
More informationMETALLIZATION AND DISSOCIATION IN HIGH PRESSURE LIQUID HYDROGEN BY AN EFFICIENT MOLECULAR DYNAMICS WITH QUANTUM MONTE CARLO
METALLIZATION AND DISSOCIATION IN HIGH PRESSURE LIQUID HYDROGEN BY AN EFFICIENT MOLECULAR DYNAMICS WITH QUANTUM MONTE CARLO Author: Guglielmo Mazzola Supervisor: Prof. Sandro Sorella A thesis submitted
More informationOptimization of quantum Monte Carlo wave functions by energy minimization
Optimization of quantum Monte Carlo wave functions by energy minimization Julien Toulouse, Roland Assaraf, Cyrus J. Umrigar Laboratoire de Chimie Théorique, Université Pierre et Marie Curie and CNRS, Paris,
More informationResonating Valence Bond wave function with molecular orbitals: application to diatomic molecules
Resonating Valence Bond wave function with molecular orbitals: application to diatomic molecules M. Marchi 1,2, S. Azadi 2, M. Casula 3, S. Sorella 1,2 1 DEMOCRITOS, National Simulation Center, 34014,
More informationFixed-Node quantum Monte Carlo for Chemistry
Fixed-Node quantum Monte Carlo for Chemistry Michel Caffarel Lab. Physique et Chimie Quantiques, CNRS-IRSAMC, Université de Toulouse e-mail : caffarel@irsamc.ups-tlse.fr. p.1/29 The N-body problem of Chemistry
More informationTowards ab-initio simulation of high temperature superconductors
Towards ab-initio simulation of high temperature superconductors Sandro Sorella SISSA, DEMOCRITOS, Trieste Santa Barbara 1 September 2009 Motivations Model systems, nice reduction schemes but how much
More informationComputational methods: Coupled Electron Ion Monte Carlo Path Integral Monte Carlo Examples: Electron gas Hydrogen
! Computational methods: Coupled Electron Ion Monte Carlo Path Integral Monte Carlo Examples: Electron gas Hydrogen WHO DID THE WORK? Miguel Morales, Livermore Carlo Pierleoni: L Aquila, Italy Jeff McMahon
More informationMetal-insulator transition with Gutzwiller-Jastrow wave functions
Metal-insulator transition with Gutzwiller-Jastrow wave functions Odenwald, July 20, 2010 Andrea Di Ciolo Institut für Theoretische Physik, Goethe Universität Frankfurt Metal-insulator transition withgutzwiller-jastrow
More informationSpin liquid phases in strongly correlated lattice models
Spin liquid phases in strongly correlated lattice models Sandro Sorella Wenjun Hu, F. Becca SISSA, IOM DEMOCRITOS, Trieste Seiji Yunoki, Y. Otsuka Riken, Kobe, Japan (K-computer) Williamsburg, 14 June
More informationQuantum Monte Carlo wave functions and their optimization for quantum chemistry
Quantum Monte Carlo wave functions and their optimization for quantum chemistry Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France CEA Saclay, SPhN Orme des Merisiers April 2015 Outline
More informationQMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters
QMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters Idoia G. de Gurtubay TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay. Quantum
More informationVariational Monte Carlo Optimization and Excited States
Variational Monte Carlo Optimization and Excited States Eric Neuscamman August 9, 2018 motivation charge transfer core spectroscopy double excitations the menu aperitif: number counting Jastrows main course:
More informationSize-extensive wave functions for QMC A linear-scaling GVB approach
Size-extensive wave functions for QMC A linear-scaling GVB approach Claudia Filippi, University of Twente, The Netherlands Francesco Fracchia, University of Pisa, Italy Claudio Amovilli, University of
More informationRecent advances in quantum Monte Carlo for quantum chemistry: optimization of wave functions and calculation of observables
Recent advances in quantum Monte Carlo for quantum chemistry: optimization of wave functions and calculation of observables Julien Toulouse 1, Cyrus J. Umrigar 2, Roland Assaraf 1 1 Laboratoire de Chimie
More informationElectronic structure simulations of water solid interfaces
Electronic structure simulations of water solid interfaces Angelos Michaelides London Centre for Nanotechnology & Department of Chemistry, University College London www.chem.ucl.ac.uk/ice Main co-workers:
More informationOptimization of quantum Monte Carlo (QMC) wave functions by energy minimization
Optimization of quantum Monte Carlo (QMC) wave functions by energy minimization Julien Toulouse, Cyrus Umrigar, Roland Assaraf Cornell Theory Center, Cornell University, Ithaca, New York, USA. Laboratoire
More informationQuantum Monte Carlo methods
Quantum Monte Carlo methods Lubos Mitas North Carolina State University Urbana, August 2006 Lubos_Mitas@ncsu.edu H= 1 2 i i 2 i, I Z I r ii i j 1 r ij E ion ion H r 1, r 2,... =E r 1, r 2,... - ground
More informationAuxiliary-field quantum Monte Carlo calculations of excited states and strongly correlated systems
Auxiliary-field quantum Monte Carlo calculations of excited states and strongly correlated systems Formally simple -- a framework for going beyond DFT? Random walks in non-orthogonal Slater determinant
More informationCorrelation in correlated materials (mostly transition metal oxides) Lucas K. Wagner University of Illinois at Urbana-Champaign
Correlation in correlated materials (mostly transition metal oxides) Lucas K. Wagner University of Illinois at Urbana-Champaign Understanding of correlated materials is mostly phenomenological FN- DMC
More informationSupplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water
Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water Thomas D. Kühne 1, 2 and Rustam Z. Khaliullin 1, 1 Institute of Physical
More informationQuantum Monte Carlo backflow calculations of benzene dimers
Quantum Monte Carlo backflow calculations of benzene dimers Kathleen Schwarz*, Cyrus Umrigar**, Richard Hennig*** *Cornell University Department of Chemistry, **Cornell University Department of Physics,
More informationarxiv: v3 [physics.comp-ph] 17 Jun 2014
The new Resonating Valence Bond Method for ab-initio Electronic Simulations Sandro Sorella Scuola Internazionale Superiore di Studi Avanzati (SISSA) and Democritos National Simulation Center, Istituto
More informationTime-dependent linear-response variational Monte Carlo.
Time-dependent linear-response variational Monte Carlo. Bastien Mussard bastien.mussard@colorado.edu https://mussard.github.io/ Julien Toulouse julien.toulouse@upmc.fr Sorbonne University, Paris (web)
More informationFermionic tensor networks
Fermionic tensor networks Philippe Corboz, Institute for Theoretical Physics, ETH Zurich Bosons vs Fermions P. Corboz and G. Vidal, Phys. Rev. B 80, 165129 (2009) : fermionic 2D MERA P. Corboz, R. Orus,
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationCCSD(T) benchmarks of non-equilibrium water clusters: the importance of monomer deformation
CCSD(T) benchmarks of non-equilibrium water clusters: the importance of monomer deformation Biswajit Santra 1, Angelos Michaelides 1,2, and Matthias Scheffler 1 1 Fritz-Haber-Institut der MPG, Berlin,
More informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationNew applications of Diffusion Quantum Monte Carlo
New applications of Diffusion Quantum Monte Carlo Paul R. C. Kent (ORNL) Graphite: P. Ganesh, J. Kim, C. Park, M. Yoon, F. A. Reboredo (ORNL) Platinum: W. Parker, A. Benali, N. Romero (ANL), J. Greeley
More informationAn Overview of Quantum Monte Carlo Methods. David M. Ceperley
An Overview of Quantum Monte Carlo Methods David M. Ceperley Department of Physics and National Center for Supercomputing Applications University of Illinois Urbana-Champaign Urbana, Illinois 61801 In
More informationarxiv: v1 [cond-mat.str-el] 18 Jul 2007
arxiv:0707.2704v1 [cond-mat.str-el] 18 Jul 2007 Determination of the Mott insulating transition by the multi-reference density functional theory 1. Introduction K. Kusakabe Graduate School of Engineering
More informationAb initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014
Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid
More informationAb initio molecular dynamics
Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field
More informationarxiv: v1 [physics.chem-ph] 2 Sep 2013
Molecular properties by Quantum Monte Carlo: an investigation on the role of the wave function ansatz and the basis set in the water molecule Andrea Zen Dipartimento di Fisica, La Sapienza - Università
More informationQuantum Monte Carlo tutorial. Lucas K. Wagner Dept. of Physics; University of Illinois at Urbana-Champaign
Quantum Monte Carlo tutorial Lucas K. Wagner Dept. of Physics; University of Illinois at Urbana-Champaign QMC: what is it good for? Theoretical gap (ev) 5 4 3 2 1 FN-DMC DFT(PBE) VO 2 (monoclinic) La 2
More informationProperties of reactive oxygen species by quantum Monte Carlo
Properties of reactive oxygen species by quantum Monte Carlo The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published
More informationAnion-π and π-π cooperative interactions
Chapter 5 Anion-π and π-π cooperative interactions 5.1 Introduction The design of selective receptors of anionic species is a very active area of research within supramolecular chemistry due to the potential
More informationPhonon wavefunctions and electron phonon interactions in semiconductors
Phonon wavefunctions and electron phonon interactions in semiconductors Bartomeu Monserrat bm418@cam.ac.uk University of Cambridge Quantum Monte Carlo in the Apuan Alps VII QMC in the Apuan Alps VII Bartomeu
More informationAb initio molecular dynamics
Ab initio molecular dynamics Molecular dynamics Why? allows realistic simulation of equilibrium and transport properties in Nature ensemble averages can be used for statistical mechanics time evolution
More informationAdvanced Electronic Structure Theory Density functional theory. Dr Fred Manby
Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)
More informationLecture 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationAb initio molecular dynamics: Basic Theory and Advanced Met
Ab initio molecular dynamics: Basic Theory and Advanced Methods Uni Mainz October 30, 2016 Bio-inspired catalyst for hydrogen production Ab-initio MD simulations are used to learn how the active site
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationQuantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid
Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid N. D. Drummond and R. J. Needs TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge
More informationAccelerated Quantum Molecular Dynamics
Accelerated Quantum Molecular Dynamics Enrique Martinez, Christian Negre, Marc J. Cawkwell, Danny Perez, Arthur F. Voter and Anders M. N. Niklasson Outline Quantum MD Current approaches Challenges Extended
More informationQuantum Monte Carlo Benchmarks Density Functionals: Si Defects
Quantum Monte Carlo Benchmarks Density Functionals: Si Defects K P Driver, W D Parker, R G Hennig, J W Wilkins (OSU) C J Umrigar (Cornell), R Martin, E Batista, B Uberuaga (LANL), J Heyd, G Scuseria (Rice)
More informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationElectronic structure quantum Monte Carlo methods and variable spins: beyond fixedphase/node
Electronic structure quantum Monte Carlo methods and variable spins: beyond fixedphase/node approximations Cody Melton, M. Chandler Bennett, L. Mitas, with A. Ambrosetti, F. Pederiva North Carolina State
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 informationthe abdus salam international centre for theoretical physics
united nations educational, scientific and cultural organization the abdus salam international centre for theoretical physics international atomic energy agency SMR 1595-16 Joint DEMOCRITOS - ICTP School
More informationThe Nature of the Interlayer Interaction in Bulk. and Few-Layer Phosphorus
Supporting Information for: The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus L. Shulenburger, A.D. Baczewski, Z. Zhu, J. Guan, and D. Tománek, Sandia National Laboratories, Albuquerque,
More informationHigh Temperature High Pressure Properties of Silica From Quantum Monte Carlo
High Temperature High Pressure Properties of Silica From Quantum Monte Carlo K.P. Driver, R.E. Cohen, Z. Wu, B. Militzer, P. Lopez Rios, M. Towler, R. Needs, and J.W. Wilkins Funding: NSF, DOE; Computation:
More informationMolecular Mechanics. I. Quantum mechanical treatment of molecular systems
Molecular Mechanics I. Quantum mechanical treatment of molecular systems The first principle approach for describing the properties of molecules, including proteins, involves quantum mechanics. For example,
More informationBinding energy of 2D materials using Quantum Monte Carlo
Quantum Monte Carlo in the Apuan Alps IX International Workshop, 26th July to 2nd August 2014 The Apuan Alps Centre for Physics @ TTI, Vallico Sotto, Tuscany, Italy Binding energy of 2D materials using
More informationAll-electron quantum Monte Carlo calculations for the noble gas atoms He to Xe
All-electron quantum Monte Carlo calculations for the noble gas atoms He to Xe A. Ma, N. D. Drummond, M. D. Towler, and R. J. Needs Theory of Condensed Matter Group, Cavendish Laboratory, University of
More informationNew Perspective on structure and bonding in water using XAS and XRS
New Perspective on structure and bonding in water using XAS and XRS Anders Nilsson Stanford Synchrotron Radiation Laboratory (SSRL) and Stockholm University, Sweden R. Ludwig Angew. Chem. 40, 1808 (2001)
More informationThe electronic structure of materials 1
Quantum mechanics 2 - Lecture 9 December 18, 2013 1 An overview 2 Literature Contents 1 An overview 2 Literature Electronic ground state Ground state cohesive energy equilibrium crystal structure phase
More informationQMC in the Apuan Alps IV Vallico Sotto, 28th July 2008
Van der Waals forces in graphitic nanostructures Norbert Nemec QMC in the Apuan Alps IV Vallico Sotto, 28th July 2008 1 History Johannes Diderik van der Waals (1837-1923): Postulation of intermolecular
More informationAb-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University
Ab-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University Hands-on Tutorial Workshop on Ab-Initio Molecular Simulations, Fritz- Haber-Institut, Berlin, July 12-21,
More informationSupporting information. Realizing Two-Dimensional Magnetic Semiconductors with. Enhanced Curie Temperature by Antiaromatic Ring Based
Supporting information Realizing Two-Dimensional Magnetic Semiconductors with Enhanced Curie Temperature by Antiaromatic Ring Based Organometallic Frameworks Xingxing Li and Jinlong Yang* Department of
More informationStrong Correlation Effects in Fullerene Molecules and Solids
Strong Correlation Effects in Fullerene Molecules and Solids Fei Lin Physics Department, Virginia Tech, Blacksburg, VA 2461 Fei Lin (Virginia Tech) Correlations in Fullerene SESAPS 211, Roanoke, VA 1 /
More informationComputational Approaches to Quantum Critical Phenomena ( ) ISSP. Fermion Simulations. July 31, Univ. Tokyo M. Imada.
Computational Approaches to Quantum Critical Phenomena (2006.7.17-8.11) ISSP Fermion Simulations July 31, 2006 ISSP, Kashiwa Univ. Tokyo M. Imada collaboration T. Kashima, Y. Noda, H. Morita, T. Mizusaki,
More informationXANES spectra of Bulk Water (Ph. Wernet, A. Nilsson, et al. Science 304, 995 (2004)) Pre-edge peak Post-edge
XANES spectra of Bulk Water (Ph. Wernet, A. Nilsson, et al. Science 304, 995 (2004)) Pre-edge peak Post-edge Transition State Potential Method Consider an electronic transition: Initial state Final state
More informationBEC-BCS Crossover in Cold Atoms
BEC-BCS Crossover in Cold Atoms (2 years later...) Andrew Morris Pablo López Ríos Richard Needs Theory of Condensed Matter Cavendish Laboratory University of Cambridge TTI 31 st July 2009 Outline Theory
More informationAb initio Molecular Dynamics Born Oppenheimer and beyond
Ab initio Molecular Dynamics Born Oppenheimer and beyond Reminder, reliability of MD MD trajectories are chaotic (exponential divergence with respect to initial conditions), BUT... With a good integrator
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 informationWhen (2 + 2) 4? Alexandre Tkatchenko Fritz Haber Institut der Max Planck Gesellschaft, Berlin, Germany
When (2 + 2) 4? Alexandre Tkatchenko Fritz Haber Institut der Max Planck Gesellschaft, Berlin, Germany QMC@TTI, Apuan Alps, Jul 29, 2013 When (2 + 2) 4 or Van der Waals Interactions in Complex (and Simple)
More informationFermionic tensor networks
Fermionic tensor networks Philippe Corboz, ETH Zurich / EPF Lausanne, Switzerland Collaborators: University of Queensland: Guifre Vidal, Roman Orus, Glen Evenbly, Jacob Jordan University of Vienna: Frank
More informationIs a system of fermions in the crossover BCS-BEC. BEC regime a new type of superfluid?
Is a system of fermions in the crossover BCS-BEC BEC regime a new type of superfluid? Finite temperature properties of a Fermi gas in the unitary regime Aurel Bulgac,, Joaquin E. Drut, Piotr Magierski
More informationOrbital-dependent backflow transformations in quantum Monte Carlo
transformations in quantum Monte Carlo P. Seth, P. López Ríos, and R. J. Needs TCM group, Cavendish Laboratory, University of Cambridge 5 December 2012 VMC and DMC Optimization Wave functions Variational
More informationAb initio molecular dynamics : BO
School on First Principles Simulations, JNCASR, 2010 Ab initio molecular dynamics : BO Vardha Srinivasan IISER Bhopal Why ab initio MD Free of parametrization and only fundamental constants required. Bond
More informationTheoretical Study of High Temperature Superconductivity
Theoretical Study of High Temperature Superconductivity T. Yanagisawa 1, M. Miyazaki 2, K. Yamaji 1 1 National Institute of Advanced Industrial Science and Technology (AIST) 2 Hakodate National College
More informationRustam Z. Khaliullin University of Zürich
Rustam Z. Khaliullin University of Zürich Molecular dynamics (MD) MD is a computational method for simulating time evolution of a collection of interacting atoms by numerically integrating Newton s equation
More informationAb initio molecular dynamics and nuclear quantum effects
Ab initio molecular dynamics and nuclear quantum effects Luca M. Ghiringhelli Fritz Haber Institute Hands on workshop density functional theory and beyond: First principles simulations of molecules and
More informationarxiv: v1 [cond-mat.str-el] 17 Aug 2016
arxiv:1608.04901v1 [cond-mat.str-el] 17 Aug 2016 Evolutionary algorithm based configuration interaction approach Rahul Chakraborty 1 1, a) and Debashree Ghosh Physical and Materials Chemistry Division,
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 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 informationIs the homogeneous electron gas homogeneous?
Is the homogeneous electron gas homogeneous? Electron gas (jellium): simplest way to view a metal homogeneous and normal Hartree-Fock: simplest method for many-electron systems a single Slater determinant
More informationAn Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method
Fakultät für Mathematik und Naturwissenschaften - Lehrstuhl für Physikalische Chemie I / Theoretische Chemie An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method Jan-Ole Joswig
More informationMichael W. Mahoney Department of Physics, Yale University, New Haven, Connecticut 06520
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 23 15 DECEMBER 2001 Quantum, intramolecular flexibility, and polarizability effects on the reproduction of the density anomaly of liquid water by simple potential
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Yesterday you (should have) learned this: http://upload.wikimedia.org/wikipedia/commons/e/ea/ Simple_Harmonic_Motion_Orbit.gif 1. deterministic
More information2D Bose and Non-Fermi Liquid Metals
2D Bose and Non-Fermi Liquid Metals MPA Fisher, with O. Motrunich, D. Sheng, E. Gull, S. Trebst, A. Feiguin KITP Cold Atoms Workshop 10/5/2010 Interest: A class of exotic gapless 2D Many-Body States a)
More informationChapter 2 Quantum chemistry using auxiliary field Monte Carlo
Chapter 2 Quantum chemistry using auxiliary field Monte Carlo 1. The Hubbard-Stratonovich Transformation 2. Neuhauser s shifted contour 3. Calculation of forces and PESs 4. Multireference AFMC 5. Examples
More informationComputational Chemistry. An Introduction to Molecular Dynamic Simulations
Computational Chemistry An Introduction to Molecular Dynamic Simulations Computational chemistry simulates chemical structures and reactions numerically, based in full or in part on the fundamental laws
More informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
More informationThe Overhauser Instability
The Overhauser Instability Zoltán Radnai and Richard Needs TCM Group ESDG Talk 14th February 2007 Typeset by FoilTEX Introduction Hartree-Fock theory and Homogeneous Electron Gas Noncollinear spins and
More informationExploring deep Earth minerals with accurate theory
Exploring deep Earth minerals with accurate theory K.P. Driver, R.E. Cohen, Z. Wu, B. Militzer, P. Lopez Rios, M. Towler, R. Needs, and J.W. Wilkins Funding: NSF, DOE; Computation: NCAR, TeraGrid, NCSA,
More informationComputational Physics. J. M. Thijssen
Computational Physics J. M. Thijssen Delft University of Technology CAMBRIDGE UNIVERSITY PRESS Contents Preface xi 1 Introduction 1 1.1 Physics and computational physics 1 1.2 Classical mechanics and statistical
More informationThe nature of superfluidity in the cold atomic unitary Fermi gas
The nature of superfluidity in the cold atomic unitary Fermi gas Introduction Yoram Alhassid (Yale University) Finite-temperature auxiliary-field Monte Carlo (AFMC) method The trapped unitary Fermi gas
More informationDFT EXERCISES. FELIPE CERVANTES SODI January 2006
DFT EXERCISES FELIPE CERVANTES SODI January 2006 http://www.csanyi.net/wiki/space/dftexercises Dr. Gábor Csányi 1 Hydrogen atom Place a single H atom in the middle of a largish unit cell (start with a
More informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
More informationAb Ini'o Molecular Dynamics (MD) Simula?ons
Ab Ini'o Molecular Dynamics (MD) Simula?ons Rick Remsing ICMS, CCDM, Temple University, Philadelphia, PA What are Molecular Dynamics (MD) Simulations? Technique to compute statistical and transport properties
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More informationCHE3935. Lecture 4 Quantum Mechanical Simulation Methods Continued
CHE3935 Lecture 4 Quantum Mechanical Simulation Methods Continued 1 OUTLINE Review Introduction to CPMD MD and ensembles The functionals of density functional theory Return to ab initio methods Binding
More informationCoupled Electron-Ion Monte Carlo of high pressure hydrogen
RPMBT 14; July 16-20, 2007, Barcelona Spain p. 1/6 Coupled Electron-Ion Monte Carlo of high pressure hydrogen Carlo Pierleoni Physics Department, University of L Aquila, Italy RPMBT 14; July 16-20, 2007,
More informationAnalysis of the ultrafast dynamics of the silver trimer upon photodetachment
J. Phys. B: At. Mol. Opt. Phys. 29 (1996) L545 L549. Printed in the UK LETTER TO THE EDITOR Analysis of the ultrafast dynamics of the silver trimer upon photodetachment H O Jeschke, M E Garcia and K H
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 informationJoint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012
2327-3 Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data 23-27 January 2012 Qunatum Methods for Plasma-Facing Materials Alain ALLOUCHE Univ.de Provence, Lab.de la Phys.
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 informationFrom real materials to model Hamiltonians with density matrix downfolding
From real materials to model Hamiltonians with density matrix downfolding (Converting a many-electron problem to a fewer-electron one) Hitesh J. Changlani Florida State University Based on HJC, H. Zheng,
More information