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
#ifndef _chemistry_qc_scf_tcscf_h
00029
#define _chemistry_qc_scf_tcscf_h
00030
00031
#ifdef __GNUC__
00032
#pragma interface
00033
#endif
00034
00035
#include <chemistry/qc/scf/scf.h>
00036
00037
namespace sc {
00038
00039
00040
00041
class TCSCF:
public SCF {
00042
protected:
00043
int user_occupations_;
00044
int tndocc_;
00045
int nirrep_;
00046
int *ndocc_;
00047
int osa_;
00048
int osb_;
00049
00050
double occa_;
00051
double occb_;
00052
00053
double ci1_;
00054
double ci2_;
00055
00056 ResultRefSymmSCMatrix focka_;
00057 ResultRefSymmSCMatrix fockb_;
00058 ResultRefSymmSCMatrix ka_;
00059 ResultRefSymmSCMatrix kb_;
00060
00061
public:
00062 TCSCF(StateIn&);
00063 TCSCF(
const Ref<KeyVal>&);
00064 ~TCSCF();
00065
00066
void save_data_state(StateOut&);
00067
00068
void print(std::ostream&o=ExEnv::out0())
const;
00069
00070
double occupation(
int ir,
int vectornum);
00071
double alpha_occupation(
int irrep,
int vectornum);
00072
double beta_occupation(
int irrep,
int vectornum);
00073
00074
int n_fock_matrices() const;
00075 RefSymmSCMatrix fock(
int);
00076 RefSymmSCMatrix effective_fock();
00077 RefSymmSCMatrix density();
00078 RefSymmSCMatrix alpha_density();
00079 RefSymmSCMatrix beta_density();
00080
00081
void symmetry_changed();
00082
00083
int spin_polarized();
00084
00085 protected:
00086
00087 RefSymmSCMatrix cl_dens_;
00088 RefSymmSCMatrix cl_dens_diff_;
00089 RefSymmSCMatrix op_densa_;
00090 RefSymmSCMatrix op_densa_diff_;
00091 RefSymmSCMatrix op_densb_;
00092 RefSymmSCMatrix op_densb_diff_;
00093
00094 RefSymmSCMatrix ao_gmata_;
00095 RefSymmSCMatrix ao_gmatb_;
00096 RefSymmSCMatrix ao_ka_;
00097 RefSymmSCMatrix ao_kb_;
00098
00099 RefSymmSCMatrix cl_hcore_;
00100
00101
void set_occupations(const RefDiagSCMatrix& evals);
00102
00103
00104
void init_vector();
00105
void done_vector();
00106
void reset_density();
00107
double new_density();
00108
double scf_energy();
00109
00110 Ref<SCExtrapData> extrap_data();
00111
00112
00113
void init_gradient();
00114
void done_gradient();
00115
00116 RefSymmSCMatrix lagrangian();
00117 RefSymmSCMatrix gradient_density();
00118
00119
00120
void init_hessian();
00121
void done_hessian();
00122 };
00123
00124 }
00125
00126 #endif
00127
00128
00129
00130
00131