9
1
mirror of https://github.com/QuantumPackage/qp2.git synced 2024-06-16 01:25:18 +02:00
qp2/src/mo_optimization
Anthony Scemama 4f296efb66
Some checks failed
continuous-integration/drone/push Build is failing
Fixed qp_export_as_tgz
2023-12-05 17:19:47 +01:00
..
org Fixed qp_export_as_tgz 2023-12-05 17:19:47 +01:00
83.mo_optimization.bats add mo optimization 2023-04-18 13:56:30 +02:00
cipsi_orb_opt.irp.f working on casscf_cipsi 2023-10-06 11:28:20 +02:00
class.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
constants.h add mo optimization 2023-04-18 13:56:30 +02:00
debug_gradient_list_opt.irp.f removed UGLY NON ASCII CHARACTERS 2023-06-29 18:31:48 +02:00
debug_gradient_opt.irp.f removed UGLY NON ASCII CHARACTERS 2023-06-29 18:31:48 +02:00
debug_hessian_list_opt.irp.f removed UGLY NON ASCII CHARACTERS 2023-06-29 18:31:48 +02:00
debug_hessian_opt.irp.f removed UGLY NON ASCII CHARACTERS 2023-06-29 18:31:48 +02:00
diagonal_hessian_list_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
diagonal_hessian_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
diagonalization_hessian.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
EZFIO.cfg Change defaults in mo_optimization 2023-08-22 11:27:36 +02:00
first_diagonal_hessian_list_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
first_diagonal_hessian_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
first_gradient_list_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
first_gradient_opt.irp.f Changed formats E to ES 2023-07-04 22:17:31 +02:00
first_hessian_list_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
first_hessian_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
gradient_list_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
gradient_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
hessian_list_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
hessian_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
NEED add mo optimization 2023-04-18 13:56:30 +02:00
orb_opt.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
README.md add mo optimization 2023-04-18 13:56:30 +02:00
run_orb_opt_trust_v2.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
save_energy.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
state_average_energy.irp.f OpenMP in MO optimization 2023-08-22 10:49:59 +02:00
state_weight_normalization.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
update_parameters.irp.f add mo optimization 2023-04-18 13:56:30 +02:00
update_st_av_ci_energy.irp.f add mo optimization 2023-04-18 13:56:30 +02:00

Orbital optimization

Methods

Different methods are available:
- full hessian

qp set orbital_optimization optimization_method full  
  • diagonal hessian
qp set orbital_optimization optimization_method diag  
  • identity matrix
qp set orbital_optimization optimization_method none  

After the optimization the ezfio contains the optimized orbitals

For a fixed number of determinants

To optimize the MOs for the actual determinants:

qp run orb_opt

For a complete optimization, i.e, with a larger and larger wave function

To optimize the MOs with a larger and larger wave function:

qp run optimization  

The results are stored in the EZFIO in “mo_optimization/result_opt”, with the following format:
(1) (2) (3) (4)
1: Number of determinants in the wf,
2: Cispi energy before the optimization,
3: Cipsi energy after the optimization,
4: Energy difference between (2) and (3).

The optimization process if the following: - we do a first cipsi step to obtain a small number of determinants in the wf - we run an orbital optimization for this wf - we do a new cipsi step to double the number of determinants in the wf - we run an orbital optimization for this wf - … - we do that until the energy difference between (2) and (3) is
smaller than the targeted accuracy for the cispi (targeted_accuracy_cipsi in qp edit) or the wf is larger than a given size (n_det_max_opt in qp_edit) - after that you can reset your determinants (qp reset -d) and run a clean Cispi calculation

End of the optimization

You can choos the number of determinants after what the optimization will stop:

qp set orbital_optimization n_det_max_opt 1e5 # or any number

Weight of the states

You can change the weights of the differents states directly in qp edit.
It will affect ths weights used in the orbital optimization.

Tests

To run the tests:

qp test

Org files

The org files are stored in the directory org in order to avoid overwriting on user changes. The org files can be modified, to export the change to the source code, run

./TANGLE_org_mode.sh
mv *.irp.f ../.