/usr/share/psi/samples/rasci-ne/input.dat is in psi4-data 1:0.3-5.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 | #! Ne atom RASCI/cc-pVQZ
#! Example of split-virtual 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 split-virtual 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.
#
memory 250 mb
# First, get the natural orbitals
molecule ne {
Ne
}
set globals {
wfn detci
basis cc-pVQZ
d_convergence 10
# docc [2, 0, 0, 0, 0, 1, 1, 1]
frozen_docc [1, 0, 0, 0, 0, 0, 0, 0]
ex_level 2
}
set detci {
nat_orbs_write true
}
set transqt2 {
delete_tei false
}
psioh = psi4.IOManager.shared_object()
psioh.set_specific_retention(32, True)
psioh.set_specific_retention(33, True)
psioh.set_specific_retention(35, True)
thisenergy = energy('detci')
clean()
# !Warning: DGAS intentionally broke the second half here.
## Now, run the RASCI
#
#set globals {
# wfn detci
# d_convergence 10
# docc [2, 0, 0, 0, 0, 1, 1, 1]
# frozen_docc [1, 0, 0, 0, 0, 0, 0, 0]
# ras1 [1, 0, 0, 0, 0, 1, 1, 1]
# ras2 [1, 0, 0, 0, 0, 1, 1, 1]
# ras3 [5, 2, 2, 2, 1, 3, 3, 3]
# ex_level 2
# val_ex_level 2
#}
#
#set detci {
# maxiter 25
# nat_orbs_write false
#
# # This keyword is at the heart of the "split-virtual" CISD[TQ].
# # If electrons are in RAS 4, do not allow the excitation to go beyond
# # that specified by ex_level (here, only singles & doubles allowed if
# # RAS 4 is populated)
# mixed4 false
#
# # this will just reduce the number of strings generated in the
# # computation; because of the mixed4=false specification above, not
# # all of the strings normally allowed by the RAS will actually be needed
# # in this computation
# r4s true
#}
#
#thisenergy = energy('detci',bypass_scf='True')
#
|