Questions tagged [hartree-fock]
For any questions about (or related to) the Hartree-Fock method.
44
questions
8
votes
1
answer
328
views
How to detect an "open shell" molecule?
I feel like this is a really stupid question, and I've been ashamed to ask for a while...
When I'm about to run some electronic structure calculation on a molecule, how do I know if I need to use an ...
- 183
5
votes
1
answer
207
views
Hartree Fock Algorithm too high values for Molecules
I am working on a high school graduation paper on Hartree Fock and have gotten pretty far already. With help from this community, I managed to generate all the necessary integrals and implement the ...
- 125
5
votes
1
answer
106
views
Implementation of one-elecron integrals in a Hartree-Fock program
Absolute newbie in quantum chemistry here. I was recently looking through some Fortran libraries from my university to see if I understand all of the mathematics. In a Hartree-Fock library, I found ...
- 53
2
votes
0
answers
67
views
How to read two electron integrals from pyscf
This is a follow-up question to this one I asked a few days ago. I in the meantime talked to my teacher and generate all the necessary integrals in PySCF. For the one electron integrals that isn't an ...
- 125
5
votes
1
answer
91
views
What is the best algorithm for ERIs of contracted gaussian atomic orbitals
I am currently working on python script that does Hartree-Fock-Calculations. I try to avoid packages as good as possible. I've figured out how to calculate the one electron integrals when it comes to ...
- 125
9
votes
2
answers
156
views
What is the physical meaning of the term "non-local interaction"?
This term pops up through out Hartree-Fock and DFT. Does this term have any physical significance to actual(real) system or it is a mathematical by-product?
- 121
5
votes
1
answer
220
views
Continuously-varying Hartree-Fock transition matrix?
Suppose we have a parametrized molecular Hamiltonian $H(R)$, where $R$ are the nuclear coordinates. If we let $B(R)$ be some chosen atomic basis set (like STO-3G) centred at $R$, then using $B(R)$, we ...
- 399
4
votes
1
answer
97
views
Molecular orbital values on grid points in PySCF?
I am looking for a way to easily evaluate individual molecular orbitals on a grid(assuming a single determinant method for now) in PySCF.
I am aware on how to generate efficient grids for real-space ...
- 178
4
votes
0
answers
76
views
Is it possible to obtain the current-voltage characteristic for a crystal using the Hartree–Fock (or DFT) calculation?
In continuation of the question of obtaining an electrical conductivity tensor (ECT) from the “first characteristics”, will my following suggestion be correct? Since a crystal with an applied-to-it ...
- 371
7
votes
2
answers
123
views
Simple question about Hartree-Fock algorithm
Classically in Hartree-Fock implementations we define the wave-function in term of a basis set:
$$\phi_i = \sum_{\alpha} C_{\alpha i}\chi_\alpha$$
If I understand correctly the problem is reduced to ...
- 971
3
votes
0
answers
73
views
Can the electrical conductivity tensor of a crystal be obtained from a Hartree-Fock calculation?
Let’s assume that a crystal has been calculated by Hartree–Fock (HF) method (or DFT one), i.d., so to speak, the "first characteristics", which are eigenvalues and eigenvectors, as well as ...
- 371
4
votes
0
answers
87
views
In which part of the Hartree–Fock algorithm (or DFT one) could a neural network be most effectively used?
I understand that the task of implementing machine learning in DFT and Hartree–Fock (HF) algorithm has already been solved, perhaps to some extent, but it is interesting to think about how to ...
- 371
3
votes
0
answers
41
views
Rys Quadrature method for integrating over Gaussian basis functions: 1D integral recursion relation
I'm trying to implement the Rys Quadrature method for integrations of Gaussians, from scratch. I'm using this paper as a reference.
So, how does the overall recursion relation for evaluating 1D ...
- 123
4
votes
2
answers
551
views
Which formalism is better suited to model a crystal with defects: Hartree–Fock, Kohn–Sham, or something else?
Let’s suppose that
(1) it’s required to calculate a crystal with certain parameters of a crystal lattice;
(2) this crystal has some defects (let it be Frenkel defects for simplicity);
(3) these ...
- 371
7
votes
1
answer
145
views
Hartree-Fock Method: Projection of density matrix onto a larger basis set
I am implementing a Hartree-Fock SCF program with an educational purpose. Although the code works fine using a Core Hamiltonian guess for the SCF iterative process, I want to implement a more accurate ...
- 73
5
votes
1
answer
82
views
About the order of indices in basis set expansion
In many books, e.g., Szabo and Ostlund, Helgaker Jorgensen and Olsen, etc etc, the expansion of Hartree-Fock orbital by basis set is written as a kind of
$$
\phi_i = \sum_{\mu} C_{\mu i} \chi_{\mu}
$$
...
- 172
17
votes
4
answers
799
views
What are the physical reasons if the SCF doesn't converge?
I've been mostly using ab initio methods as a routine calculation tool, so even though I know some things from the theoretical side, I still have mostly hands-on experience. One of the things that I ...
user4626
7
votes
0
answers
144
views
Evaluating Coulomb and exchange integrals in practice using cartesian coordinates
Many software packages provide means to calculate the Coulomb and exchange integrals, for example, the exchange
$$\int \phi^*(1) \, \psi^*(2) \frac{1}{|\mathbf{r_1} - \mathbf{r_2}|} \, \psi(1) \, \phi(...
- 593
7
votes
1
answer
142
views
Exploiting 8-fold symmetry of ERI tensor for building Coulomb and Exchange matrices
I'm trying to write a restricted Hartree Fock code in Fortran that reads in a file of zeroth-iteration 1 and 2 electron integrals (FCIDUMP format) and uses them to do the SCF procedure from an initial ...
- 71
8
votes
1
answer
349
views
What do size-extensivity and size-consistency mean?
I have heard both terms in various lectures and books on quantum chemistry, however, I have not found a proper explanation of them. As I understand now, size consistency of a method means for example, ...
- 5,834
6
votes
1
answer
338
views
How to compute the overlap matrix in Python
I want to calculate the overlap integral (S), I made the code, but it is only worked when I use a basis function that describe 1s and when I include 2s orbital, I start to get error.
The First Problem ...
7
votes
0
answers
66
views
Is the formula of a Dyson orbital independent of whether there is degeneracy or not in the HOMO of the neutral?
A Dyson orbital is defined as:
$$
\textrm{Dy}(x) = \int dx_1 \ldots dx_{N-1} \, (\Psi^+(x_1,\ldots,x_{N-1}))^* \, \Psi^0(x_1,\ldots,x_{N-1}, x),\tag{1}
$$
where $\Psi^0$ and $\Psi^+$ are the total ...
- 595
7
votes
1
answer
95
views
How to read data from an input and form two-electrons integral?
I am doing a code that can calculate the hartree-fock energy. I want to form the two-electrons integral by reading the data from the input like this one below. The first four number are the indices ...
6
votes
0
answers
123
views
Differences between DFTB3 energy and HF energy
I was performing DFTB3 calculation with DFTB+ on water monomer. This is my input.
...
17
votes
3
answers
1k
views
Why is uncertainty not a big problem in computational chemistry?
The molecular Hamiltonian (or, for simplicity, the Fock operator) contains coulomb potentials as well as momentum operators. To evaluate the coulomb potential, we need to know where the electrons are. ...
- 447
12
votes
2
answers
476
views
For DFT and Hartree-Fock, how can we know that we have a true minimum? Is there an equivalent to the "frequency analysis" for geometry optimization?
In a geometry optimisation, one can check if they have reached a minimum by checking the phonons. And to a certain extent, the global minimum can be confirmed by relaxing from different starting ...
- 1,829
8
votes
1
answer
695
views
Convert reduced density matrix from M.O. to A.O. basis
I am interested in converting the one-particle reduced matrix (rdm1) from the molecular orbital (M.O.) basis to the atomic orbital (A.O.) basis. Is the following method correct (in an identical ...
- 1,829
9
votes
2
answers
658
views
Visualize electron density using pyscf
I found the following example of code which uses density functional theory to compute the electron density $\rho$:
...
9
votes
1
answer
308
views
'Exchange' in Hartree-Fock and Kohn-Sham DFT
I am a bit confused about the notions of exchange in the Hartree-Fock and Kohn-Sham Density Functional Theory schemes.
In Hartree-Fock, one just writes down the many-electron wavefunction as a hartree-...
- 2,336
8
votes
3
answers
490
views
Does Kohn-Sham DFT use Slater determinants?
In the Hartree method, it is known that the wavefunction of the system does not obey the antisymmetry principle of fermions - that is when you swap two particles, they don't up a negative sign. ...
- 2,336
16
votes
1
answer
495
views
Why is CPHF/CPKS necessary for calculating second derivatives?
This question is coming from an answer to one of my previous questions. During optimizations, QM programs usually compute the gradient(first derivative) analytically, and take a guess of the hessian (...
- 5,834
14
votes
2
answers
279
views
Optimization of Gaussian basis sets within the Hartree-Fock Method
I am revisiting some exercises in Thijssen's Computational Physics book, particularly chapter 4 on the Hartree-Fock method. I am interested in the method of nonlinear optimisation for its own purposes,...
- 731
6
votes
1
answer
481
views
Why is my quantum chemistry software not printing as many 2e- integrals as I expect?
Background:
In Hartree-Fock theory, the two-electron integral is given by:
$[ij|kl]$ = $\int dx_{1}dx_{2}(\chi_{i}(x_{1})\chi_{j}(x_{1})\frac{1}{r_{12}}\chi_{k}^{*}(x_2)\chi_{l}(x_{2})$
I am ...
- 1,829
9
votes
2
answers
336
views
Normalization constant and Roothaan Equations
Roothaan Hall Equations:
The Hartree-Fock equations are a set of modified Schrodinger equations:
$f_{i}\psi_{m}=\epsilon_{m}\psi_{m}$
where:
The Fock operator ($f_{i}$) is given by (restricted case):
...
- 1,829
8
votes
1
answer
250
views
Semi-canonicalisation vs canonicalisation of the Fock matrix and orbitals
I have seen the terms semi-canonicalized and canonicalized used in relation to the Fock matrix, density matrices, and orbitals; however, I am unsure what these terms actually describe.
For instance:
...
- 1,829
9
votes
1
answer
122
views
Muonic Modeling
I have seen an article[1] that describes a HF approach to working with exotic muonic atoms. I am curious about how far DFT has been extended for modeling muonic chemistry. Are there psuedopotentials ...
- 10.2k
18
votes
5
answers
797
views
What are the types of SCF?
Many of us know the most common types of SCF
(though we can do better than Wikipedia at explaining them):
RHF (Restricted Hartree-Fock), RKS (Restricted Kohn-Sham) [link to answer]
UHF (Unrestricted ...
- 29.1k
12
votes
1
answer
215
views
Analog computing in matter modeling today: Any applications?
When we think about computers and computing today, almost invariably we think about digital computers, the ever-present objects today. But reading this question in the Chemistry stackexchange, How are ...
- 3,728
11
votes
1
answer
688
views
How to calculate the Fock matrix in the molecular orbital basis PySCF?
I am interested in calculating the Fock matrix in the molecular orbital basis with PySCF, though I am not clear on the methodology behind this task.
In my attempt, I use the following script (for the ...
- 1,829
12
votes
2
answers
818
views
What do the rows and columns of a Fock matrix represent?
I am using the pyscf code, where the Fock matrix can be obtained by:
...
- 1,829
15
votes
2
answers
644
views
Hartree-Fock density vs Kohn-Sham density
Hartree-Fock density is free of self-interaction but lacks electron correlation effects, while the density from KS-DFT (using an xc functional or potential, both which are explicitly density-dependent ...
- 936
12
votes
2
answers
356
views
One-center two-electron integrals between 1s STO
As per @user1271772 's suggestion I am asking this question here again.
I am trying to understand SCF cycle by trying to code up solved example from Quantum Chemistry by Levine (page 443, 5th edition)...
- 305
12
votes
2
answers
192
views
What is the mathematical condition that ensure that the self-consistent field (SCF) procedure must converge?
We know that there is convergence issue of DFT method on study real molecular system. It would be conducive if people could pre-determine if the SCF procedure is converge analytically before starting ...
- 3,609
16
votes
2
answers
697
views
How does the Hartree-Fock method improve on Hartree after considering the wavefunction that takes the anti-symmetry property into account?
I am not able to understand the math of this.
After taking the Slater determinant as a wavefunction in Hartree-Fock, does the procedure to find wavefunctions remain the same as in the Hartree method ...
- 656