00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029 #ifndef _chemistry_qc_dft_hsosks_h
00030 #define _chemistry_qc_dft_hsosks_h
00031
00032 #ifdef __GNUC__
00033 #pragma interface
00034 #endif
00035
00036 #include <chemistry/qc/scf/hsosscf.h>
00037 #include <chemistry/qc/dft/integrator.h>
00038 #include <chemistry/qc/dft/functional.h>
00039
00040 namespace sc {
00041
00042
00043
00044 class HSOSKS: public HSOSSCF {
00045 protected:
00046 Ref<DenIntegrator> integrator_;
00047 Ref<DenFunctional> functional_;
00048 RefSymmSCMatrix vxc_a_;
00049 RefSymmSCMatrix vxc_b_;
00050
00051 RefSymmSCMatrix cl_vxc();
00052 RefSymmSCMatrix op_vxc();
00053
00054 public:
00055 HSOSKS(StateIn&);
00056 HSOSKS(const Ref<KeyVal>&);
00057 ~HSOSKS();
00058
00059 void save_data_state(StateOut&);
00060
00061 void print(std::ostream&o=ExEnv::out0()) const;
00062
00063 void two_body_energy(double &ec, double &ex);
00064
00065 int value_implemented() const;
00066 int gradient_implemented() const;
00067
00068 protected:
00069 void ao_fock(double accuracy);
00070 double exc_;
00071 double scf_energy();
00072 Ref<SCExtrapData> extrap_data();
00073 RefSymmSCMatrix effective_fock();
00074
00075 void init_vector();
00076 void done_vector();
00077
00078 void two_body_deriv(double*);
00079
00080 RefSymmSCMatrix lagrangian();
00081 };
00082
00083 }
00084
00085 #endif
00086
00087
00088
00089
00090