

(6 intermediate revisions not shown) 
Line 1: 
Line 1: 
  =Comprehensive Suite of Scalable Capabilities=
 +  #REDIRECT [[Release66:Capabilities]] 
   +  
  NWChem provides many methods for computing the properties of molecular and periodic systems using standard quantum mechanical descriptions of the electronic wavefunction or density. Its classical molecular dynamics capabilities provide for the simulation of macromolecules and solutions, including the computation of free energies using a variety of force fields. These approaches may be combined to perform mixed quantummechanics and molecularmechanics simulations.
 +  
   +  
  The specific methods for determining molecular electronic structure, molecular dynamics, and pseudopotential planewave electronic structure and related attributes are listed in the following sections.
 +  
   +  
  ==Molecular Electronic Structure==
 +  
   +  
  Methods for determining energies and analytic first derivatives with respect to atomic coordinates include the following:
 +  
   +  
  • HartreeFock (RHF, UHF, highspin ROHF)
 +  
  • Gaussian orbitalbased density functional theory (DFT) using many local and nonlocal exchangecorrelation potentials (LDA, LSDA)
 +  
  • secondorder perturbation theory (MP2) with RHF and UHF references
 +  
  • complete active space selfconsistent field theory (CASSCF).
 +  
   +  
  Analytic second derivatives with respect to atomic coordinates are available for RHF and UHF, and closedshell DFT with all functionals.
 +  
   +  
  The following methods are available to compute energies only:
 +  
   +  
  • iterative CCSD, CCSDT, and CCSDTQ methods and their EOMCC counterparts for RHF, ROHF, and UHF references
 +  
  • activespace CCSDt and EOMCCSDt approaches
 +  
  • completely renormalized CRCCSD(T), and CREOMCCSD(T) correction to EOMCCSD excitation energies
 +  
  • locally renormalized CCSD(T) and CCSD(TQ) approaches
 +  
  • noniterative approaches based on similarity transformed Hamiltonian: the CCSD(2)T and CCSD(2) formalisms.
 +  
  • MP2 with RHF reference and resolution of the identity integral approximation MP2 (RIMP2) with RHF and UHF references
 +  
  • selected CI with secondorder perturbation correction.
 +  
   +  
  For all methods, the following may be performed:
 +  
   +  
  • single point energy calculations
 +  
  • geometry optimization with constraints (minimization and transition state)
 +  
  • molecular dynamics on the fully ab initio potential energy surface
 +  
  • automatic computation of numerical first and second derivatives
 +  
  • normal mode vibrational analysis in Cartesian coordinates
 +  
  • ONIOM hybrid calculations
 +  
  • ConductorLike Screening Model (COSMO) calculations
 +  
  • electrostatic potential from fit of atomic partial charges
 +  
  • spinfree oneelectron DouglasKroll calculations
 +  
  • electron transfer (ET)
 +  
  • vibrational SCF and DFT.
 +  
   +  
  At the SCF and DFT level of theory various (response) properties are available, including NMR shielding tensors and indirect spinspin coupling.
 +  
   +  
  ==Quantum Mechanics/Molecular Mechanics (QM/MM)==
 +  
   +  
  The QM/MM module in NWChem provides a comprehensive set of capabilities to study ground and excited state properties of largemolecular systems. The QM/MM module can be used with practically any quantum mechanical method available in NWChem. The following tasks are supported
 +  
   +  
  • single point energy and property calculations
 +  
  • excited states calculation
 +  
  • optimizations and transition state search
 +  
  • dynamics
 +  
  • free energy calculations.
 +  
   +  
  ==Pseudopotential PlaneWave Electronic Structure==
 +  
   +  
  The NWChem PlaneWave (NWPW) module uses pseudopotentials and planewave basis sets to perform DFT calculations. This method's efficiency and accuracy make it a desirable first principles method of simulation in the study of complex molecular, liquid, and solidstate systems. Applications for this first principles method include the calculation of free energies, search for global minima, explicit simulation of solvated molecules, and simulations of complex vibrational modes that cannot be described within the harmonic approximation.
 +  
   +  
  The NWPW module is a collection of three modules:
 +  
   +  
  • PSPW (PSeudopotential PlaneWave) A gamma point code for calculating molecules, liquids, crystals, and surfaces.
 +  
  • Band A band structure code for calculating crystals and surfaces with small band gaps (e.g. semiconductors and metals).
 +  
  • PAW (Projector Augmented Wave) a gamma point projector augmented planewave code for calculating molecules, crystals, and surfaces.
 +  
   +  
  These capabilities are available:
 +  
   +  
  • constant energy and constant temperature CarParrinello molecular dynamics (extended Lagrangian dynamics)
 +  
  • LDA, PBE96, and PBE0, exchangecorrelation potentials (restricted and unrestricted)
 +  
  • SIC, pertOEP, HartreeFock, and hybrid functionals (restricted and unrestricted)
 +  
  • Hamann, TroullierMartins, HartwigsenGoedeckerHutter normconserving pseudopotentials with semicore corrections
 +  
  • geometry/unit cell optimization, frequency, transitionstates
 +  
  • fractional occupation of molecular orbitals for metals
 +  
  • AIMD/MM capability in PSPW
 +  
  • constraints needed for potential of mean force (PMF) calculation
 +  
  • wavefunction, density, electrostatic, Wannier plotting
 +  
  • band structure and density of states generation
 +  
   +  
  ==Molecular Dynamics==
 +  
   +  
  The NWChem Molecular Dynamics (MD) module can perform classical simulations using the AMBER and CHARMM force fields, quantum dynamical simulations using any of the quantum mechanical methods capable of returning gradients, and mixed quantum mechanics molecular dynamics simulation and molecular mechanics energy minimization.
 +  
   +  
  Classical molecular simulation functionality includes the following methods:
 +  
   +  
  • single configuration energy evaluation
 +  
  • energy minimization
 +  
  • molecular dynamics simulation
 +  
  • free energy simulation (MCTI and MSTP with single or dual topologies, doublewide sampling, and separationshifted scaling).
 +  
   +  
  The classical force field includes the following elements:
 +  
   +  
  • effective pair potentials
 +  
  • firstorder polarization
 +  
  • selfconsistent polarization
 +  
  • smooth particle mesh Ewald
 +  
  • twinrange energy and force evaluation
 +  
  • periodic boundary conditions
 +  
  • SHAKE constraints
 +  
  • constant temperature and/or pressure ensembles
 +  
  • dynamic proton hopping using the QHOP methodology
 +  
  • advanced system setup capabilities for biomolecular membranes.
 +  