mirror of
https://github.com/mveril/qp-demo
synced 2024-11-03 20:54:01 +01:00
50 lines
1.0 KiB
Plaintext
Executable File
50 lines
1.0 KiB
Plaintext
Executable File
#!/usr/bin/env qpsh
|
|
|
|
# In this third tutorial, we will compute the ground state and the lowest
|
|
# singly excited of N2.
|
|
|
|
# 1. Create the EZFIO database, with the def2-svp basis set:
|
|
|
|
qp create_ezfio -b def2-svp n2.zmt
|
|
|
|
# 2. Run the Hartree-Fock calculation:
|
|
|
|
qp run scf | tee scf.out
|
|
|
|
# 3. Set the core MOs as frozen:
|
|
|
|
qp set_frozen_core
|
|
|
|
# 4. Request for two states
|
|
|
|
qp set determinants n_states 2
|
|
|
|
# 5. Run a CIS
|
|
|
|
qp run cis
|
|
|
|
# 6. Truncate the wave function to the 2 most important determinants : the
|
|
# Hartree-Fock and one a singly excited determinant (spin symmetry will
|
|
# be automatically restored in the next calculation)
|
|
|
|
qp edit -n 2
|
|
|
|
# 7. Set the CIS wave function as the starting trial wave function
|
|
|
|
qp set determinants read_wf True
|
|
|
|
# 8. To avoid collapsing to the ground state and a doubly excited state,
|
|
# force to use a state-following algorithm
|
|
|
|
qp set davidson state_following True
|
|
|
|
# 9. Terminate when we have 50000 determinants
|
|
|
|
qp set determinants n_det_max 50*10**3
|
|
|
|
# 10. Run the selected FCI calculation
|
|
|
|
qp run fci | tee fci.out
|
|
|
|
|