Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Some Failed ACI-DSRG Combinations #228

Open
lcyyork opened this issue Feb 19, 2021 · 0 comments
Open

Some Failed ACI-DSRG Combinations #228

lcyyork opened this issue Feb 19, 2021 · 0 comments

Comments

@lcyyork
Copy link
Contributor

lcyyork commented Feb 19, 2021

I modified aci-14 test case and found the following sci_excited_algorithm may not work with state-averaged DSRG at the moment.

The input:

molecule C2{
C  0.00000 0.00000 0.0
C  0.00000 0.00000 1.24253
}

set {
  basis 6-31g*
  scf_type pk
  e_convergence 13
  d_convergence 12
  r_convergence 12
  docc [2,0,0,0,0,2,1,1]
  guess gwh
}

set forte {
  active_space_solver aci
  frozen_docc [2,0,0,0,0,0,0,0]
  restricted_docc [0,0,0,0,0,0,0,0]
  active [1,0,1,1,0,3,2,2]
  sci_excited_algorithm ALGORITHM
  sigma 0.001
#  nroot 2
#  aci_n_average 2
  charge 0
  root_sym 0
  sci_preiterations 2
  dl_maxiter 100
  SCI_SAVE_FINAL_WFN true
  correlation_solver dsrg-mrpt2
  calc_type sa
  avg_state [[0,1,2]]
}

Escf, wfn = energy('scf', return_wfn=True)
energy('forte', ref_wfn=wfn)

ROOT_COMBINE: seg fault due to empty final_wfn_.
ROOT_ORTHOGONALIZE: this runs but the energy is wrong: the DSRG reference energy is not equal to the averaged energy from ACI. May need to set the eigenvectors and set of determinants from old_roots_.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant