Test Suite and Sample Inputs¶
PSI4 is distributed with an extensive test suite, which can
be found in psi4/tests. After building the source code, these
can automatically be run by running ctest
in the compilation
directory. More info on ctest
options can be found
here. Sample input files
can be found in the psi4/samples subdirectory of the toplevel Psi
directory. The samples and a brief description are provided below.
Sample inputs accessible through interfaced executables are bulleted below.
Sample inputs for PSI4 as distributed are below.
Input File 
Description 

Accesses basis sets, databases, plugins, and executables in noninstall locations 

ROHFCCSD/ccpVDZ \(^{3}B_1\) CH2 geometry optimization via analytic gradients 

BH single points, checking that program can run multiple instances of DETCI in a single input, without an intervening clean() call 

OMP3 ccpCVDZ energy with ROHF initial guess for the NO radical 

SCF STO3G finitedifference tests 

CC2(RHF)/ccpVDZ energy of H2O. 

Test of the superposition of atomic densities (SAD) guess, using a highly distorted water geometry with a ccpVDZ basis set. This is just a test of the code and the user need only specify guess=sad to the SCF module’s (or global) options in order to use a SAD guess. The test is first performed in C2v symmetry, and then in C1. 

SCF STO3G finitedifferences frequencies from gradients for H2O 

Triple and Singlet Oxygen energy SOSCF, also tests nonsymmetric density matrices 

MP2.5 ccpVDZ gradient for the H2O molecule. 

RHF CCSD(T) ccpVDZ frozencore energy of C4NH4 Anion 

tdwb97x excitation energies of singlet states of h2o, wfn passing 

MkMRCCSD(T) single point. \(^1A_1\) O$_3` state described using the Ms = 0 component of the singlet. Uses TCSCF orbitals. 

SCF with various combinations of pk/densityfitting, castup/nocastup, and spherical/cartesian settings. Demonstrates that puream setting is getting set by orbital basis for all df/castup parts of calc. Demonstrates that answer doesn’t depend on presence/absence of castup. Demonstrates (by comparison to castup3) that output file doesn’t depend on options (scf_type) being set global or local. This input uses global. 

CC3/ccpVDZ H2O \(R_e\) geom from Olsen et al., JCP 104, 8007 (1996) 

Compute the dipole polarizability for water with custom basis set. 

DC06 calculation for the O2 molecule (triplet ground state). This performs geometry optimization using twostep and simultaneous solution of the response equations for the analytic gradient. 

MP2/augccpvDZ many body energies of an arbitrary Helium complex, addressing 4body formulas 

density fitted REMP/ccpVDZ energies for the CH3 radical 

OMP2 ccpVDZ energy for the H2O molecule. 

ZAPT(n)/631G NH2 Energy Point, with n=225 

Test initial SCF guesses on FH and FH+ in ccpVTZ basis 

Test case for Binding Energy of C4H5N (Pyrrole) with CO2 using MP2/def2TZVPP 

This test case shows an example of running and analyzing a standard FSAPT0/junccpvdz procedure for phenol dimer from the S22 database. 

ROHFEOMCCSD/DZ on the lowest two states of each irrep in \(^{3}B_1\) CH2. 

This test case shows an example of running and analyzing a difference FSAPT0/junccpvdz procedure for phenol dimer from the S22 database. 

tdcamb3lyp with DiskDF and method/basis specification 

This checks that all energy methods can run with a minimal input and set symmetry. 

OMP2.5 ccpVDZ energy for the H2O molecule. 

analog of fsaptextabc with molecule and external potentials in Bohr 

DFCCSDL ccpVDZ energy for the H2O molecule. 

CASSCF/631G** energy point 

631G** H2O+ Test CISD Energy Point 

SCF level shift on an ROHF computation 

CONV SCF 631G analytical vs finitedifference tests Tests UHF hessian code for Ca != Cb 

Test parsed and exotic calls to energy() like zapt4, mp2.5, and cisd are working 

UHF Dipole Polarizability Test 

Sample HF/ccpVDZ H2O computation all derivatives 

Test that Python Molecule class processes geometry like psi4 Molecule class. 

SAPT0(ROHF) openshell computation of CN  Ne interaction energy First with junccpVDZ and density fitted integrals with ROHF Then with ccpVDZ and direct integrals, except for dispersion that is computed with ccpVDZri density fitting with ROHF. 

DFCCSD(AT) ccpVDZ energy for the H2O molecule. 

Lithium test for coverage 

tduhf test on triplet states of methylene (tda), wfn passing 

ROHFCCSD(T) ccpVDZ energy for the \(^2\Sigma^+\) state of the CN radical, with Zmatrix input. 

mtd/basis syntax examples 

CASSCF/631G** energy point 

DF SCF 631G UHFl vs RHF test Tests DF UHF hessian code for Ca = Cb 

Tests all grid pruning options available and screening of small weights. Check against grid size. 

DFT (hybrids) test of implementations in: hybrid_superfuncs.py 

HF/ccpVDZ many body energies of an arbitrary noble gas trimer complex Size vs cost tradeoff is rough here 

This is a shorter version if isapt1  does not do cube plots. See isapt1 for full details 

optimization with method defined via cbs 

External potential calculation involving a TIP3P water and a QM water for DFMP2. Finite different test of the gradient is performed to validate forces. 

Test of SFX2C1e on water uncontracted ccpVDZDK The reference numbers are from Lan Cheng’s implementation in Cfour 

631G H2O Test FCI Energy Point 

Extrapolated water energies  densityfitted version 

Convergence of manybody gradients of different BSSE schemes 

631G** H2O CCSD optimization by energies, with ZMatrix input 

Optimize H2O HF/ccpVDZ 

tdwb97x singlet excitation energies of methylene (tda) 

Multilevel computation of water trimer energy (geometry from J. Chem. Theory Comput. 11, 21262136 (2015)) 

Intercalls among python wrappers database, cbs, optimize, energy, etc. Though each call below functions individually, running them all in sequence or mixing up the sequence is aspirational at present. Also aspirational is using the intended types of gradients. 

tdcamb3lyp with DiskDF and method/basis specification 

DFCCSD ccpVDZ gradient for the NH molecule. 

RHF orbitals and density for water. 

DFT Functional Smoke Test 

Vibrational and thermo analysis of water trimer (geometry from J. Chem. Theory Comput. 11, 21262136 (2015)) 

WaterArgon complex with ECP present; check of energies and forces. 

Patch of a glycine with a methyl group, to make alanine, then DFSCF energy calculation with the ccpVDZ basis set 

DFBP86D2 ccpVDZ frozen core gradient of S22 HCN updated ref gradient due to new BraggSlater radii 

631G H2O Test FCI Energy Point 

External potential calculation involving a TIP3P water and a QM water. Finite different test of the gradient is performed to validate forces. 

SCF level shift on a UHF computation 

RKS Linear Exchange Algorithm test for benzene 

Single point gradient of 12B1 state of H2O+ with EOMCCSD 

Compute the IRC for HCN <> NCH interconversion at the RHF/DZP level of theory. 

A rangeseperated gradient for SO2 to test disk algorithms by explicitly setting low memory 

Single point gradient of 12B2 state of H2O+ with EOMCCSD 

RHFCCSD(T) ccpVQZ frozencore energy of the BH molecule, with Cartesian input. After the computation, the checkpoint file is renamed, using the PSIO handler. 

sapt example with orbital freezing with alkali metal and dMP2 

631G(d) optimization of SF4 starting from linear bond angle that is not linear in the optimized structure but is in a symmetry plane of the molecule. 

Benzene vertical singlettriplet energy difference computation, using the PubChem database to obtain the initial geometry, which is optimized at the HF/STO3G level, before computing single point energies at the RHF, UHF and ROHF levels of theory. 

Patch of a glycine with a methyl group, to make alanine, then DFSCF energy calculation with the ccpVDZ basis set 

DSDPBEP86 S22 Ammonia test 

RHFCC2LR/ccpVDZ optical rotation of H2O2. gauge = both, omega = (589 355 nm) 

MP(n)/augccpVDZ BH Energy Point, with n=219. Compare against M. L. Leininger et al., J. Chem. Phys. 112, 9213 (2000) 

OMP3 ccpVDZ gradient for the H2O molecule. 

Extrapolated energies with delta correction 

OMP2 ccpVDZ energy for the H2O molecule. 

RHF interaction energies using nbody and cbs parts of the driver Ne dimer with mp2/v[dt]z + d:ccsd(t)/vdz 

ROHFCCSD ccpVDZ energy for the \(^2\Sigma^+\) state of the CN radical 

DFMP2 frequency by difference of energies for H2O 

RHFCCSD 631G** allelectron optimization of the H2O molecule 

OMP3 ccpVDZ gradient for the NO radical 

OMP2 ccpVDZ energy for the NO molecule. 

SCF DZ finite difference frequencies by energies for C4NH4 

Restricted DFDCT ODC12 energies with linearly dependent basis functions 

SCF/ccpVDZ optimization example with frozen cartesian 

RHF ccpVQZ energy for the BH molecule, with Cartesian input. 

check SP basis Fortran exponent parsing 

DCT calculation for the HF+ using DC06 functional. This performs both twostep and simultaneous update of the orbitals and cumulant using DIIS extrapolation. Fourvirtual integrals are first handled in the MO Basis for the first two energy computations. In the next two the ao_basis=disk algorithm is used, where the transformation of integrals for fourvirtual case is avoided. The computation is then repeated using the DC12 functional with the same algorithms. 

631G H2O Test for coverage 

SCF with various combinations of pk/densityfitting, castup/nocastup, and spherical/cartesian settings. Demonstrates that puream setting is getting set by orbital basis for all df/castup parts of calc. Demonstrates that answer doesn’t depend on presence/absence of castup. Demonstrates (by comparison to castup2) that output file doesn’t depend on options (scf_type) being set global or local. This input uses local. 

Check that C++ Molecule class and qcdb molecule class are reading molecule input strings identically 

OLCCD ccpVDZ gradient for the NO radical 

usapt example with empty beta due to frozen core 

DFT Functional Test for RangeSeperated Hybrids and Ghost atoms 

check that CC is returning the same values btwn CC*, FNOCC, and DFOCC modules 

EOMCC3(ROHF) on CH radical with userspecified basis and properties for particular root 

Computation of VMFCcorrected water trimer gradient (geometry from J. Chem. Theory Comput. 11, 21262136 (2015)) 

SOSOMP3 ccpVDZ geometry optimization for the H2O molecule. 

MkMRPT2 single point. \(^1A_1\) F2 state described using the Ms = 0 component of the singlet. Uses TCSCF singlet orbitals. 

SCF STO3G finitedifference frequencies from energies for H2O 

Extrapolated water energies 

DFOMP3 ccpVDZ energy for the H2O+ cation 

density fitted REMP/ccpVDZ energies for the CO2 molecule. 

631G H2O Test FCI Energy Point 

testing aligner on enantiomers based on Table 1 of 10.1021/ci100219f aka J Chem Inf Model 2010 50(12) 21292140 

Compute the IRC for HOOH torsional rotation at the RHF/DZP level of theory. 

CCSD/ccpVDZ dipole polarizability at two frequencies 

Various gradients for a strained helium dimer and water molecule 

RHF Linear Exchange Algorithm test for water 

Test LDA stability analysis against QChem. 

Vibrational and thermo analysis of several water isotopologs. Demonstrates Hessian reuse for different temperatures, pressures, and isotopologs 

631G** H2O Test CISD Energy Point with subspace collapse 

Triple and Singlet Oxygen energy SOSCF, also tests nonsymmetric density matrices 

External potential calculation with one Ghost atom and one point charge at the same position. 

This test case shows an example of running and analyzing an FISAPT0/junccpvdz computation for 2,4pentanediol (targeting the intramolecular hydrogen bond between the two hydroxyl groups) 

DC06, DC12, ODC06 and ODC12 calculation for the He dimer. This performs a simultaneous update of the orbitals and cumulant, using DIIS extrapolation. Fourvirtual integrals are handled in the MO Basis. 

RHFCCSDLR/ccpVDZ static polarizability of HOF 

Density fitted MP2 energy of H2, using density fitted reference and automatic looping over ccpVDZ and ccpVTZ basis sets. Results are tabulated using the built in table functions by using the default options and by specifiying the format. 

Extrapolated water energies  conventional integrals version 

Test of SFX2C1e on Water uncontracted ccpVDZ The reference numbers are from Lan Cheng’s implementation in Cfour 

Computation of NoCPcorrected water trimer gradient (geometry from J. Chem. Theory Comput. 11, 21262136 (2015)) 

Compute the dipole, quadrupole, and traceless quadrupoles for water. 

631G* C2 Test RASCI Energy Point, testing two different ways of specifying the active space, either with the ACTIVE keyword, or with RAS1, RAS2, RESTRICTED_DOCC, and RESTRICTED_UOCC 

Matches Table II aCCSD(T)/ccpVDZ H2O @ 2.5 * Re value from Crawford and Stanton, IJQC 98, 601611 (1998). 

SCF STO3G geometry optimzation, with Zmatrix input, by finitedifferences 

RHF STO3G (Cartesian) and ccpVDZ (spherical) water Hessian test, against Psi3 reference values. 

OLCCD ccpVDZ gradient for the H2O molecule. 

RHFCCSD(T) ccpVQZ frozencore energy of the BH molecule, with Cartesian input. This version tests the FROZEN_DOCC option explicitly 

EOMCCSD/631g excited state transition data for water with two excited states per irrep 

EOMCC2/ccpVDZ on H2O2 with two excited states in each irrep 

DFT Functional Test 

Tests RHF/ROHF/UHF SCF gradients 

Test if the the guess read in the same basis converges. 

SCF DZ allene geometry optimization, with Cartesian input, first in c2v symmetry, then in Cs symmetry from a starting point with a nonlinear central bond angle. 

tduhf test on triplet states of methylene (rpa) 

MBIS calculation on NaCl 

Example potential energy surface scan and CPcorrection for Ne2 

Tests to determine full point group symmetry. Currently, these only matter for the rotational symmetry number in thermodynamic computations. 

Test Gibbs free energies at 298 K of N2, H2O, and CH4. 

SCF DZ allene geometry optimzation, with Cartesian input 

631G** H2O Test RASSCF Energy Point will default to only singles and doubles in the active space 

SAPT(DFT) augccpVDZ interaction energy between Ne and Ar atoms. 

MOM excitation from LUMO HOMO+4 

SAPT(DFT) augccpVDZ interaction energy between Ne and Ar atoms. 

Check flavors of B3LYP (b3lyp3/b3lyp5) against other programs 

RHF STO3G dipole moment computation, performed by applying a finite electric field and numerical differentiation. 

BHH2+ FCI/ccpVDZ Transition Dipole Moment 

An example of using BLAS and LAPACK calls directly from the Psi input file, demonstrating 

OMP3 ccpVDZ energy for the H2O molecule 

Cholesky decomposed REMP/ccpVDZ energies for the CH3 radical 

DFOMP2.5 ccpVDZ gradients for the H2O molecule. 

DFOMP2.5 ccpVDZ energy for the H2O molecule. 

SCF 631G(d) optimization of TS for HCN to HNC Performs finite difference hessian calculation. Then optimizes using previous orbitals for scf guess, in subsequent calculations. The last two displacements of the hessian break the plane of symemtry, This test confirms that only the reference geometry, with the correct symmetry, writes orbitals to disk. SCF will fail (ValidationError) otherwise. 

ROHF stability analysis check for CN with ccpVDZ. This test corresponds to the rohfstab test from Psi3. 

Computation of VMFCcorrected water trimer Hessian (geometry from J. Chem. Theory Comput. 11, 21262136 (2015)) 

This test case shows an example of running and analyzing a standard FSAPT0/junccpvdz procedure for HSG18dimer from the HSG database. 

TCSCF ccpVDZ energy of asymmetrically displaced ozone, with Zmatrix input. 

External potential calculation involving a TIP3P water and a QM water. Gradient on the external charges is compared to gradient on the QM atoms to validate the gradient on the charges. 

TDHF test variable access 

CC3(ROHF)/ccpVDZ H2O \(R_e\) geom from Olsen et al., JCP 104, 8007 (1996) 

MOM excitation from LUMO HOMO+3 

apply linear fragmentation algorithm to a water cluster 

DFMP2 ccpVDZ gradient for the NO molecule. 

OLCCD ccpVDZ energy with ROHF initial guess for the NO radical 

SAPT0 augccpVDZ computation of the benzenemethane interaction energy, using the augpVDZJKFIT DF basis for SCF, the augccpVDZRI DF basis for SAPT0 induction and dispersion, and the augpVDZJKFIT DF basis for SAPT0 electrostatics and induction. This example uses frozen core as well as asyncronous I/O while forming the DF integrals and CPHF coefficients. 

FSAPT with external charge on trimer 

Example of stateaveraged CASSCF for the C2 molecule see C. D. Sherrill and P. Piecuch, J. Chem. Phys. 122, 124104 (2005) 

ROHF 631G** energy of the \(^{3}B_1\) state of CH2, with Zmatrix input. The occupations are specified explicitly. 

H2 with tiny basis set, to test basis set parser’s handling of integers 

test roundtripness of dict repr for psi4.core.Molecule and qcdb.Molecule 

check that methods can act on single atom 

SAPT calculation on bimolecular complex where monomers are unspecified so driver autofragments it. Basis set and auxiliary basis sets are assigned by atom type. 

Allelectron MP2 631G** geometry optimization of water 

Test of the superposition of atomic densities (SAD) guess, using a highly distorted water geometry with a ccpVDZ basis set. This is just a test of the code and the user need only specify guess=sad to the SCF module’s (or global) options in order to use a SAD guess. The test is first performed in C2v symmetry, and then in C1. 

Spinrestricted DC06 counterpart of dct1. 

CASSCF/631G** energy point 

Maximum Overlap Method (MOM) Test. MOM is designed to stabilize SCF convergence and to target excited Slater determinants directly. 

Single point energies of multiple excited states with EOMCCSD 

DFCCSD ccpVDZ gradients for the H2O molecule. 

DFOMP3 ccpVDZ energy for the H2O molecule. 

Frozencore CCSD(T)/ccpVDZ on C4H4N anion with disk ao algorithm 

SAPT(DFT) augccpVDZ interaction energy between Ne and Ar atoms. 

CCSD dipole with userspecified basis set 

LibXC density screening test. Tests empty, Conly, Xonly and XC superfunctionals. ‘super_mix’ showcases how to use different screening values for X and C parts. SCF will fail or crash (nans) without screening! 

Test fnocc with linear dependencies 

UHFCCSD(T) ccpVDZ frozencore energy for the \(^2\Sigma^+\) state of the CN radical, with Zmatrix input. 

SCF ccpVDZ geometry optimzation, with Zmatrix input 

Tests SCF gradient in the presence of a dipole field 

MkMRCCSD single point. \(^3 \Sigma ^\) O2 state described using the Ms = 0 component of the triplet. Uses ROHF triplet orbitals. 

OLCCD ccpVDZ energy for the H2O molecule. 

HF and DFT variants singlepoints on zmat methane, mostly to test that PSI variables are set and computed correctly. Now also testing that CSX harvesting PSI variables correctly update ref_dft_2e/xc due to new BraggSlater radii 

A demonstration of mixed Cartesian/ZMatrix geometry specification, using variables, for the benzenehydronium complex. Atoms can be placed using ZMatrix coordinates, whether they belong to the same fragment or not. Note that the Cartesian specification must come before the ZMatrix entries because the former define absolute positions, while the latter are relative. 

Various constrained energy minimizations of HOOH with ccpvdz RHF Internalcoordinate constraints in internalcoordinate optimizations. 

A general test of the MintsHelper function 

Example of stateaveraged CASSCF for the C2 molecule 

CI/MCSCF ccpvDZ properties for Potassium nitrate (rocket fuel!) 

test scf castup with custom basis sets 

The multiple guesses for DCT amplitudes for ODC12. 

SAPT0 augccpVTZ computation of the charge transfer energy of the water dimer. 

RHFCCSD/ccpVDZ energy of H2O partitioned into pair energy contributions. 

DFBP86D2 ccpVDZ frozen core gradient of S22 HCN update ref gradient due to new BraggSlater radii 

sapt0 of charged system in ECP basis set 

RHFBCCD(T)/631G** H2O singlepoint energy (fzc, MObasis \(\langle abcd \rangle\)) 

Analytic vs. finite difference DFSCF frequency test for water. 

Extrapolated water energies 

SAPT2+3(CCD) augccpVDZ+midbond computation of the water dimer interaction energy, using the augccpVDZJKFIT DF basis for SCF and augccpVDZRI for SAPT. 

He2+ FCI/ccpVDZ Transition Dipole Moment 

Example SAPT computation for ethene*ethine (i.e., ethylene*acetylene), test case 16 from the S22 database 

DFMP2 ccpVDZ frozen core gradient of benzene, computed at the DFSCF ccpVDZ geometry 

DCT calculation for the NH3+ radical using the ODC12 and ODC13 functionals. This performs both simultaneous and QC update of the orbitals and cumulant using DIIS extrapolation. Fourvirtual integrals are first handled in the MO Basis for the first two energy computations. In the next computation ao_basis=disk algorithm is used, where the transformation of integrals for fourvirtual case is avoided. 

Benzene Dimer OutofCore HF/ccpVDZ 

CASSCF/631G** energy point. Check energy with frozen core/virtual orbs. after semicanonicalization. 

incremental Cholesky filtered SCF 

EOMCCSD/ccpVDZ on H2O2 with two excited states in each irrep 

631G H2O Test FCI Energy Point 

ROHFEOMCCSD/DZ analytic gradient lowest \(^{2}A_1\) excited state of H2O+ (B1 excitation) 

631G** H2O Test CISD Energy Point 

Cholesky decomposed REMP/ccpVDZ energies for the CO2 molecule. 

MP2 ccpvDZ properties for Nitrogen oxide 

CCSD/sto3g optical rotation calculation (length gauge only) at two frequencies on methyloxirane 

UHF STO3G (Cartesian) and ccpVDZ (spherical) water Hessian test, against Psi3 reference values. This test should match RHF values exactly 

EOMCCSD/631g excited state transition data for water cation 

DFOMP2.5 ccpVDZ gradients for the H2O+ cation. 

Test case for some of the PSI4 outofcore codes. The code is given only 2.0 MB of memory, which is insufficient to hold either the A1 or B2 blocks of an ovvv quantity incore, but is sufficient to hold at least two copies of an oovv quantity incore. 

CCSD/sto3g optical rotation calculation (both gauges) at two frequencies on methyloxirane 

OMP2 ccpVDZ energy for the H2O molecule. 

MP2/augccpv[DT]Z many body energies of an arbitrary Helium complex Size vs cost tradeoff is rough here 

MkMRCCSD(T) single point. \(^1A_1\) CH2 state described using the Ms = 0 component of the singlet. Uses RHF singlet orbitals. 

OMP2 ccpVDZ gradient for the H2O molecule. 

Test FNOQCISD(T) computation 

ROHFCCSD ccpVDZ frozencore energy for the \(^2\Sigma^+\) state of the CN radical, with Cartesian input. 

MP2 ccpVDZ gradient for the H2O molecule. 

Carbon/UHF FractionallyOccupied SCF Test Case 

MP3 ccpVDZ gradient for the NO radical 

EDIIS test case from 10.1063/1.1470195 

WaterArgon complex with ECP present; check of UHF Hessian 

An example of using BLAS and LAPACK calls directly from the Psi input file, demonstrating matrix multiplication, eigendecomposition, Cholesky decomposition and LU decomposition. These operations are performed on vectors and matrices provided from the Psi library. 

MkMRCCSD(T) single point. \(^1A_1\) CH2 state described using the Ms = 0 component of the singlet. Uses RHF singlet orbitals. 

DFMP2 ccpVDZ gradients for the H2O molecule. 

SAPT0 ccpVDZ computation of the etheneethyne interaction energy, using the ccpVDZJKFIT RI basis for SCF and ccpVDZRI for SAPT. Monomer geometries are specified using Cartesian coordinates. 

DFCCSD ccpVDZ gradients for the H2O molecule. 

DFT Functional Test 

Various DCT analytic gradients for the O2 molecule with 631G basis set 

SCF DZ finite difference frequencies by gradients for C4NH4 

Computation of CPcorrected water trimer gradient (geometry from J. Chem. Theory Comput. 11, 21262136 (2015)) 

SCF level shift on an RKS computation 

updated dldf reference to new BraggSlater radii Dispersionless density functional (dlDF+D) internal match to Psi4 Extensive testing has been done to match supplemental info of Szalewicz et. al., Phys. Rev. Lett., 103, 263201 (2009) and Szalewicz et. al., J. Phys. Chem. Lett., 1, 550555 (2010) 

usapt example with empty beta 

wB97XD ccpVDZ gradient of S22 HCN update df/pk_ref values due to new BraggSlater radii 

Kr–Kr nocp energies with allelectron basis set to check frozen core 

He Dimer VV10 functional test. notes: DFT_VV10_B/C overwrites the NL_DISPERSION_PARAMETERS tuple updated ‘bench’ reference values for new BraggSlater radii. 

RHFCC2LR/ccpVDZ dynamic polarizabilities of HOF molecule. 

Compute three IP and 2 EA’s for the PH3 molecule 

Database calculation, so no molecule section in input file. Portions of the full databases, restricted by subset keyword, are computed by sapt0 and dfmp2 methods. 

SCF ccpVTZ geometry optimzation, with Zmatrix input 

LCCD ccpVDZ gradient for the NO radical 

comparison of DFMP2 and DLPNOMP2 with a cartesian basis set 

DFT Functional Test all values update for new BraggSlater radii 

RHF orbitals and density for water. 

Scan fractional occupation of electrons updated values due to new BraggSlater radii 

DFT custom functional test 

Tests RHF CCSD(T)gradients 

Test SAD SCF guesses on noble gas atom 

Test of the superposition of atomic densities (SAD) guess, using a highly distorted water geometry with a ccpVDZ basis set. This is just a test of the code and the user need only specify guess=sad to the SCF module’s (or global) options in order to use a SAD guess. The test is first performed in C2v symmetry, and then in C1. 

EOMCC3/ccpVTZ on H2O 

This checks that all energy methods can run with a minimal input and set symmetry. 

SAPT2+(3) augccpVDZ computation of the formamide dimer interaction energy, using the augccpVDZJKFIT DF basis for SCF and augccpVDZRI for SAPT. This example uses frozen core as well as MP2 natural orbital approximations. 

density fitted OOREMP/ccpVDZ engrad single points for the H2O+ molecule. 

RHFCC2LR/STO3G optical rotation of (S)methyloxirane. gauge = both, omega = (589 355 nm) 

Similar to mints2, but using the BSE to specify the basis sets 

Test QCISD(T) for H2O/ccpvdz Energy 

Tests SAPT0D corrections, with a variety of damping functions/parameters 

CC2(UHF)/ccpVDZ energy of H2O+. 

Generation of NBO file 

SCSOMP2 ccpVDZ geometry optimization for the H2O molecule. 

DFCCSD(T) ccpVDZ gradient for the NH molecule. 

comparison of DFMP2 and DLPNOMP2 

Frozencore CCSD(ROHF)/ccpVDZ on CN radical with diskbased AO algorithm 

ROHF frontier orbitals of CH2(s) and CH2(t). 

OMP3 ccpCVDZ energy with B3LYP initial guess for the NO radical 

Vibrational and thermo analysis of several water isotopologs. Demonstrates Hessian reuse for different temperatures and pressures but not for different isotopologs. 

Singlepoint gradient, analytic and via finitedifferences of 21A1 state of H2O with EOMCCSD 

DCT calculation for the triplet O2 using DC06 and DC12. Only twostep algorithm is tested. 

Tests SAPT0D corrections, with a variety of damping functions/parameters 

RKS Density Matrix basedIntegral Screening Test for benzene 

CCSD/ccpVDZ optical rotation calculation (length gauge only) on Zmat H2O2 

Sample HF/ccpVDZ H2O computation 

631G** H2O+ Test CISD Energy Point 

ccpvdz H2O Test coupledpair CISD against DETCI CISD 

MBIS calculation on OH (Expanded Arrays) 

Compute three IP and 2 EA’s for the PH3 molecule 

Omega optimization for LRC functional wB97 on water 

Various basis set extrapolation tests 

CCSD Response for H2O2 

run some BLAS benchmarks 

MkMRCCSD frequencies. \(^1A_1\) O$_3` state described using the Ms = 0 component of the singlet. Uses TCSCF orbitals. 

Test of all different algorithms and reference types for SCF, on singlet and triplet O2, using the ccpVTZ basis set. 

check nonphysical masses possible 

ROHFCCSD ccpVDZ frozencore energy for the \(^2\Sigma^+\) state of the CN radical, with Cartesian input. 

CCSD/ccpVDZ optical rotation calculation (both gauges) on Cartesian H2O2 

DFT integral algorithms test, performing wB97 RKS and UKS computations on water and its cation, using all of the different integral algorithms. This tests both the ERI and ERF integrals. 

Single point gradient of 11B2 state of H2O with EOMCCSD 

Single point energies of multiple excited states with EOMCCSD 

Sample UHF/ccpVDZ H2O computation on a doublet cation, using RHF/ccpVDZ orbitals for the closedshell neutral as a guess 

UHFODC12 and RHFODC12 singlepoint energy for H2O. This performs a simultaneous update of orbitals and cumulants, using DIIS extrapolation. Fourvirtual integrals are handled in the AO basis, where integral transformation is avoided. In the next RHFODC12 computation, AO_BASIS=NONE is used, where fourvirtual integrals are transformed into MO basis. 

Check that basis sets can be input with explicit angular momentum format 

MP3 ccpVDZ gradient for the H2O molecule. 

B3LYP ccpVDZ geometry optimzation of phenylacetylene, starting from not quite linear structure updated reference due to new BraggSlater radii 

SCF/sto3g optimization with a hessian every step 

CASSCF/631G** energy point 

MkMRCCSD(T) single point. \(^1A_1\) CH2 state described using the Ms = 0 component of the singlet. Uses RHF singlet orbitals. 

DFCCSD(T) ccpVDZ energy for the NH molecule. 

MP2 with a PBE0 reference computation 

check all variety of options parsing 

Density fitted MP2 ccPVDZ/ccpVDZRI computation of formic acid dimer binding energy using explicit specification of ghost atoms. This is equivalent to the dfmp2_1 sample but uses both (equivalent) specifications of ghost atoms in a manual counterpoise correction. 

SCF STO3G geometry optimzation, with Zmatrix input 

ROHFEOMCCSD/DZ analytic gradient lowest \(^{2}B_1\) state of H2O+ (A1 excitation) 

Test FNODFCCSD(T) energy 

OMP2.5 ccpVDZ energy for the H2O molecule. 

wB97XD test for a large UKS molecule update ref gradient due to new BraggSlater radii 

DFMP2 ccpVDZ gradients for the H2O molecule. 

UHFCCSD(T) ccpVDZ frozencore energy for the \(^2\Sigma^+\) state of the CN radical, with Zmatrix input. 

manybody different levels of theory on each body of helium tetramer 

integral conventional OOREMP/ccpVDZ engrad single points for the H2O molecule. 

Tests DFMP2 gradient in the presence of a dipole field 

RHFEOMCC2/ccpVDZ lowest two states of each symmetry of H2O. 

force occupations in scf 

Tests OMP2 gradient in the presence of a dipole field 

RHF/ccpvdzdecontract HCl singlepoint energy Testing the in line decontract option for basis sets 

integral conventional OOREMP/ccpVDZ engrad single points for the H2O molecule. single point energies were independently checked using the original wavels code 

density fitted OOREMP/ccpVDZ engrad single points for the H2O+ molecule. 

Analytic SVWN frequencies, compared to finite difference values 

MBIS calculation on H2O 

ccpvdz H2O Test ACPF Energy/Properties 

reproduces dipole moments in J.F. Stanton’s “biorthogonal” JCP paper 

Superficial test of PubChem interface 

NeXe dimer MP2 energies with ECP, with electrons correlated then frozen. 

631G** H2O Test RASSCF Energy Point will default to only singles and doubles in the active space 

FSAPT0/junccpvdz procedure for methane dimer 

A test of the basis specification. A benzene atom is defined using a ZMatrix containing dummy atoms and various basis sets are assigned to different atoms. The symmetry of the molecule is automatically lowered to account for the different basis sets. 

DFOMP3 ccpVDZ gradients for the H2O molecule. 

integral conventional REMP/ccpVDZ energies for the H2O molecule. results were independently verified against the initial wavels implementation 

H2O CISD/631G** Optimize Geometry by Energies 

DFT (LDA/GGA) test of custom implementations in: gga_superfuncs.py 

A test of the basis specification. Various basis sets are specified outright and in blocks, both orbital and auxiliary. Constructs libmints BasisSet objects through the constructor that calls qcdb.BasisSet infrastructure. Checks that the resulting bases are of the right size and checks that symmetry of the Molecule observes the basis assignment to atoms. 

DFCCDL ccpVDZ energy for the H2O molecule. 

UHF>UHF stability analysis test for BH with ccpVDZ Test direct SCF with and without symmetry, test PK without symmetry 

check distributed driver is correctly passing function kwargs 

SAPT0 augccpVDZ computation of the waterwater interaction energy, using the three SAPT codes. 

Test individual integral objects for correctness. 

RHF 631G** energy of water, using the MCSCF module and Zmatrix input. 

SAPT0 openshell computation of H2OHO2 interaction energy First with ccpVDZ and density fitted integrals with UHF Then with 631g and direct integrals, except for dispersion that is computed with ccpVDZri density fitting with UHF. 

Tests CAM gradients with and without XC pieces to narrow grid error 

This checks that all energy methods can run with a minimal input and set symmetry. 

ccpvdz H2O Test CEPA(1) Energy 

CASSCF/631G** energy point 

Gradient regularized asymptotic correction (GRAC) test. 

MBIS calculation on ZnO 

Secondorder SCF convergnece: Benzene 

UFH and B3LYP ccpVQZ properties for the CH2 molecule. 

Tests RHF CCSD(T)gradients 

Test SFX2C1e with a static electric field on He augccpVTZ 

SCF ccpVDZ geometry optimzation of ketene, starting from bent structure 

631G** H2O Test CISD Energy Point 

A very quick correctness test of FSAPT (see fsapt1 for a real example) 

Test SCF dipole derivatives against old Psi3 reference values 

UHF gradient for a oneelectron system (no beta electrons). 

apply linear fragmentation algorithm to a water cluster 

Doublehybrid density functional B2PYLP. Reproduces portion of Table I in S. Grimme’s J. Chem. Phys 124 034108 (2006) paper defining the functional. 

Test method/basis with disk_df 

This test case shows an example of running the ISAPT0/augccpVDZ computation for a positively charged system, illustrating the cationpi interaction. The SIAO1 link partitioning algorithm is used. The system is taken from http://dx.doi.org/10.1016/j.comptc.2014.02.008 

UHFCCSD/ccpVDZ \(^{3}B_1\) CH2 geometry optimization via analytic gradients 

Frequencies for H2O B3LYP/631G* at optimized geometry 

DFSCF ccpVDZ of benzenehydronium ion, scanning the dissociation coordinate with Python’s builtin loop mechanism. The geometry is specified by a Zmatrix with dummy atoms, fixed parameters, updated parameters, and separate charge/multiplicity specifiers for each monomer. Oneelectron properties computed for dimer and one monomer. 

comparison of DFMP2 and DLPNOMP2 with a CBS extrapolation 

routing check on lccd, lccsd, cepa(0). 

DFCCSD(T) ccpVDZ gradients for the H2O molecule. 

OLCCD ccpVDZ energy with B3LYP initial guess for the NO radical 

CASSCF/631G** energy point 

631G** H2O+ Test CISD Energy Point 

DFMP2 frequency by difference of energies for H2O 

RHF ccpVDZ energy for water, automatically scanning the symmetric stretch and bending coordinates using Python’s builtin loop mechanisms. The geometry is specified using a Zmatrix with variables that are updated during the potential energy surface scan, and then the same procedure is performed using polar coordinates, converted to Cartesian coordinates. 

CC3(UHF)/ccpVDZ H2O \(R_e\) geom from Olsen et al., JCP 104, 8007 (1996) 

RHFODC12 analytic gradient computations for H2O use AO_BASIS=DISK and AO_BASIS=NONE, respectively. RHFODC06 analytic gradient computations for H2O use AO_BASIS=DISK and AO_BASIS=NONE, respectively. 

LCCD ccpVDZ gradient for the H2O molecule. 

SAPT(DFT) augccpVDZ interaction energy between Ne and Ar atoms. 

DC06 calculation for the He dimer. This performs a twostep update of the orbitals and cumulant, using DIIS extrapolation. Fourvirtual integrals are handled in the MO Basis. 

DC06 calculation for the He dimer. This performs a simultaneous update of the orbitals and cumulant, using DIIS extrapolation. Fourvirtual integrals are handled in the AO Basis, using integrals stored on disk. 

DFCCSD(T) ccpVDZ energy for the H2O molecule. 

ROHF and UHFBCCD(T)/ccpVDZ \(^{3}B_1\) CH2 singlepoint energy (fzc, MObasis \(\langle abcd \rangle\) ) 

Fractional occupation with symmetry 

WaterArgon complex with ECP present; check of RHF Hessian 

ROHFCCSD(T) ccpVDZ frozencore energy for the \(^2\Sigma^+\) state of the CN radical, with Cartesian input. 

RISCF ccpVTZ energy of water, with Zmatrix input and ccpVTZRI auxilliary basis. 

Advanced python example sets different sets of scf/postscf conv crit and check to be sure computation has actually converged to the expected accuracy. 

DFOMP3 ccpVDZ gradients for the H2O+ cation. 

MBIS calculation on H2O 

OMP2 ccpVDZ energy for the NO radical 

Test of SAD/Castup (mainly not dying due to file weirdness) 

MP2.5 ccpVDZ gradient for the NO radical 

Various gradients for a strained helium dimer and water molecule 

Spectroscopic constants of H2, and the full ci ccpVTZ level of theory 

Cholesky decomposed OOREMP/ccpVDZ energy for the H2O molecule. 

RHF Density Matrix basedIntegral Screening Test for water 

Numpy interface testing 

Symmetry tests for a range of molecules. This doesn’t actually compute any energies, but serves as an example of the many ways to specify geometries in Psi4. 

OMP2 ccpVDZ energy for the NO molecule. 

UHF and ROHF Linear Exchange Algorithm test for benzyl cation 

Tests analytic CC2 gradients 

EOMCC3(UHF) on CH radical with userspecified basis and properties for particular root 

UHFCCSD(T)/ccpVDZ \(^{3}B_1\) CH2 geometry optimization via analytic gradients 

File retention, docc, socc, and bond distances specified explicitly. 

Test computing values of basis functions (puream and nonpuream) at points 

FSAPT0/junccpvdz procedure for methane dimer 

SCSOMP3 ccpVDZ geometry optimization for the H2O molecule. 

DF SCF 631G analytical vs finitedifference tests Tests DF UHF hessian code for Ca != Cb 

This test case shows an example of running the ISAPT0/junccpVDZ computation for 2,4pentanediol (targeting the intramolecular hydrogen bond between the two hydroxyl groups) The SIAO1 link partitioning algorithm is used. An FSAPT partitioning follows ISAPT. 

RHFCC2LR/STO3G optical rotation of (S)methyloxirane. gauge = length, omega = (589 355 nm) 

Transitionstate optimizations of HOOH to both torsional transition states. 

OMP2 ccpVDZ energy for the NO molecule. 

Benzene Dimer DFHF/ccpVDZ 

density fitted OOREMP/ccpVDZ engrad single points for the H2O molecule. 

Sample UHF/631G** CH2 computation 

ADIIS test case, from 10.1063/1.3304922 

cc3: RHFCCSD/631G** H2O geometry optimization and vibrational frequency analysis by finitedifferences of gradients 

Test if the the guess read in the same basis converges. 

Optimization followed by frequencies H2O HF/ccpVDZ 

metaGGA gradients of water and ssh molecules reference gradients updated due to new BraggSlater radii 

Quick test of external potential in FSAPT (see fsapt1 for a real example) 

DFOMP2.5 ccpVDZ energy for the H2O+ cation 

SCF level shift on a CUHF computation 

Ne atom RASCI/ccpVQZ Example of splitvirtual CISD[TQ] from Sherrill and Schaefer, J. Phys. Chem. XXX This uses a “primary” virtual space 3s3p (RAS 2), a “secondary” virtual space 3d4s4p4d4f (RAS 3), and a “tertiary” virtual space consisting of the remaining virtuals. First, an initial CISD computation is run to get the natural orbitals; this allows a meaningful partitioning of the virtual orbitals into groups of different importance. Next, the RASCI is run. The splitvirtual CISD[TQ] takes all singles and doubles, and all triples and quadruples with no more than 2 electrons in the secondary virtual subspace (RAS 3). If any electrons are present in the tertiary virtual subspace (RAS 4), then that excitation is only allowed if it is a single or double. 

OMP2 ccpVDZ energy with ROHF initial guess orbitals for the NO radical 

conventional and densityfitting mp2 test of mp2 itself and setting scsmp2 

DFCCSD ccpVDZ energy for the H2O molecule. 

SOSOMP2 ccpVDZ geometry optimization for the H2O molecule. 

RASCI/631G** H2O Energy Point 

check mixing ECP and nonECP orbital/fitting basis sets in a session 

Various extrapolated optimization methods for the H2 molecule 

Unrestricted DFDCT ODC12 gradient for O2 with ccpVTZ/ccpVTZRI standard/auxiliary basis set 

DFOMP2 ccpVDZ gradients for the H2O molecule. 

OMP2 ccpVDZ energy for the NO molecule. 

MBIS calculation on OH radical 

DFCCD ccpVDZ energy for the H2O molecule. 

Electrostatic potential and electric field evaluated on a grid around water. 

Test FNODFCCSD(T) energy 

OMP2 ccpVDZ gradient for the NO radical 

OMP2 ccpVDZ energy for the H2O molecule. 

Various constrained energy minimizations of HOOH with ccpvdz RHF. For “fixed” coordinates, the final value is provided by the user. 

Test omega is setable updated wb97x_20,wb97x_03 to account for new BraggSlater radii 

MP2 ccpVDZ gradient for the NO radical 

FSAPT with external charge on dimer 

SAPT(DFT) augccpVDZ computation for the water dimer interaction energy. 

DFMP2 ccpVDZ gradient for the NO molecule. 

Analytic UKS SVWN frequencies, compared to finite difference values 

SAPT2+3 with S^inf exchind30 Geometries taken from the S66x10 database, the shortestrange point (R = 0.7 R_e) 

Tests CCENERGY’s CCSD gradient in the presence of a dipole field 

SAPT0 with S^inf exchdisp20 

OMP2.5 ccpVDZ gradient for the NO radical 

Computation of VMFCcorrected HF dimer Hessian 

This test case shows an example of running the ISAPT0/junccpVDZ computation for 2,4pentanediol (targeting the intramolecular hydrogen bond between the two hydroxyl groups) The SIAO1 link partitioning algorithm is used. 

Tests the Psi4 SFSAPT code 

DFT JK ondisk test 

OMP2.5 ccpVDZ gradient for the H2O molecule. 

Cholesky filter a complete basis 

RHFCC2LR/ccpVDZ static polarizabilities of HOF molecule. 

MkMRCCSD single point. \(^3 \Sigma ^\) O2 state described using the Ms = 0 component of the triplet. Uses ROHF triplet orbitals. 

631G** UHF CH2 3B1 optimization. Uses a ZMatrix with dummy atoms, just for demo and testing purposes. 

Various constrained energy minimizations of HOOH with ccpvdz RHF. Cartesiancoordinate constrained optimizations of HOOH in internals. 

RHFCC2LR/ccpVDZ optical rotation of H2O2. gauge = length, omega= (589 355 nm) 

Test G2 method for H2O 

RHF augccpVQZ energy for the BH molecule, with Cartesian input. Various gradients for a strained helium dimer and water molecule 

Various constrained energy minimizations of HOOH with ccpvdz RHF. Cartesiancoordinate constrained optimizations of HOOH in Cartesians. 

integral conventional unrestricted REMP/ccpVDZ energies for the H2O+ molecule. results were independently verified against the initial wavels implementation 

checks that all SAPT physical components (elst, exch, indc, disp) and total IE are being computed correctly for SAPT2+3(CCD)dMP2/augccpvdz and all lesser methods thereof. 

test FCIDUMP functionality for rhf/uhf 

Test of SFX2C1e on Water ccpVDZDK. In this test the Dirac equation is solved in the uncontracted ccpVDZDK basis. The reference numbers are from Lan Cheng’s implementation in Cfour 

DFSCF ccpVDZ multipole moments of benzene, up to 7th order and electrostatic potentials evaluated at the nuclear coordinates 

OLCCD ccpVDZ freqs for C2H2 

Density fitted MP2 ccPVDZ/ccpVDZRI computation of formic acid dimer binding energy using automatic counterpoise correction. Monomers are specified using Cartesian coordinates. 

DFACCSD(T) ccpVDZ energy for the NH molecule. 

Restricted DFDCT ODC12 gradient for ethylene with ccpVDZ/ccpVDZRI standard/auxiliary basis set 

631G(d) optimization of SF4 starting from linear bond angle that is not linear in the optimized structure but is in a symmetry plane of the molecule. 

DCT calculation for the triplet O2 using ODC06 and ODC12 functionals. Only simultaneous algorithm is tested. 