MPQC 2.3.1
tcscf.h
1//
2// tcscf.h --- definition of the two-configuration SCF class
3//
4// Copyright (C) 1996 Limit Point Systems, Inc.
5//
6// Author: Edward Seidl <seidl@janed.com>
7// Maintainer: LPS
8//
9// This file is part of the SC Toolkit.
10//
11// The SC Toolkit is free software; you can redistribute it and/or modify
12// it under the terms of the GNU Library General Public License as published by
13// the Free Software Foundation; either version 2, or (at your option)
14// any later version.
15//
16// The SC Toolkit is distributed in the hope that it will be useful,
17// but WITHOUT ANY WARRANTY; without even the implied warranty of
18// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19// GNU Library General Public License for more details.
20//
21// You should have received a copy of the GNU Library General Public License
22// along with the SC Toolkit; see the file COPYING.LIB. If not, write to
23// the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
24//
25// The U.S. Government is granted a limited license as per AL 91-7.
26//
27
28#ifndef _chemistry_qc_scf_tcscf_h
29#define _chemistry_qc_scf_tcscf_h
30
31#ifdef __GNUC__
32#pragma interface
33#endif
34
35#include <chemistry/qc/scf/scf.h>
36
37namespace sc {
38
39// //////////////////////////////////////////////////////////////////////////
40
41class TCSCF: public SCF {
42 protected:
43 int user_occupations_;
44 int tndocc_;
45 int nirrep_;
46 int *ndocc_;
47 int osa_;
48 int osb_;
49
50 double occa_;
51 double occb_;
52
53 double ci1_;
54 double ci2_;
55
60
61 public:
62 TCSCF(StateIn&);
63 TCSCF(const Ref<KeyVal>&);
64 ~TCSCF();
65
67
68 void print(std::ostream&o=ExEnv::out0()) const;
69
70 double occupation(int ir, int vectornum);
71 double alpha_occupation(int irrep, int vectornum);
72 double beta_occupation(int irrep, int vectornum);
73
74 int n_fock_matrices() const;
75 RefSymmSCMatrix fock(int);
76 RefSymmSCMatrix effective_fock();
80
82
84
85 protected:
86 // these are temporary data, so they should not be checkpointed
87 RefSymmSCMatrix cl_dens_;
88 RefSymmSCMatrix cl_dens_diff_;
89 RefSymmSCMatrix op_densa_;
90 RefSymmSCMatrix op_densa_diff_;
91 RefSymmSCMatrix op_densb_;
92 RefSymmSCMatrix op_densb_diff_;
93
94 RefSymmSCMatrix ao_gmata_;
95 RefSymmSCMatrix ao_gmatb_;
96 RefSymmSCMatrix ao_ka_;
97 RefSymmSCMatrix ao_kb_;
98
99 RefSymmSCMatrix cl_hcore_;
100
101 void set_occupations(const RefDiagSCMatrix& evals);
102
103 // scf things
104 void init_vector();
105 void done_vector();
106 void reset_density();
107 double new_density();
108 double scf_energy();
109
110 Ref<SCExtrapData> extrap_data();
111
112 // gradient things
113 void init_gradient();
114 void done_gradient();
115
116 RefSymmSCMatrix lagrangian();
117 RefSymmSCMatrix gradient_density();
118
119 // hessian things
120 void init_hessian();
121 void done_hessian();
122};
123
124}
125
126#endif
127
128// Local Variables:
129// mode: c++
130// c-file-style: "ETS"
131// End:
static std::ostream & out0()
Return an ostream that writes from node 0.
The RefDiagSCMatrix class is a smart pointer to an DiagSCMatrix specialization.
Definition: matrix.h:380
The RefSymmSCMatrix class is a smart pointer to an SCSymmSCMatrix specialization.
Definition: matrix.h:261
A template class that maintains references counts.
Definition: ref.h:332
The SCF class is the base for all classes that use a self-consistent field procedure to solve an effe...
Definition: scf.h:49
Restores objects that derive from SavableState.
Definition: statein.h:70
Serializes objects that derive from SavableState.
Definition: stateout.h:61
Definition: tcscf.h:41
int spin_polarized()
Return 1 if the alpha density is not equal to the beta density.
RefSymmSCMatrix alpha_density()
Return alpha electron densities in the SO basis.
void print(std::ostream &o=ExEnv::out0()) const
Print information about the object.
RefSymmSCMatrix beta_density()
Return beta electron densities in the SO basis.
RefSymmSCMatrix density()
Returns the SO density.
double occupation(int ir, int vectornum)
Returns the occupation.
double beta_occupation(int irrep, int vectornum)
Returns the beta occupation.
void save_data_state(StateOut &)
Save the base classes (with save_data_state) and the members in the same order that the StateIn CTOR ...
double alpha_occupation(int irrep, int vectornum)
Returns the alpha occupation.
void symmetry_changed()
Call this if you have changed the molecular symmetry of the molecule contained by this MolecularEnerg...

Generated at Thu Jan 19 2023 00:00:00 for MPQC 2.3.1 using the documentation package Doxygen 1.9.6.