Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [density-functional-theory]

Density functional theory (DFT) is a quantum mechanical model used to estimate the electronic structure of molecules and condensed matter. In broad terms, DFT works by treating all the electrons in the system as a single electron density, and computing physical quantities of interest as functionals of that electron density. (from Physics.SE)

-2
votes
0answers
44 views

Convergence failure [on hold]

I trying to run a DFT calculation with the following job input: aui2s2pc20h18 0 1 # PBEPBE/LANL2DZ Opt Freq NoSymm aui2s2pc20h18 0 1 Au 2.9788534146 ...
0
votes
0answers
14 views

Is the exchange interaction always less in value than coulomb interaction? (Absolute values of both.)

I wonder for 2 electron interactions whether exchange interaction (K) is always less than the coulomb interaction(J) value in all cases (HF or DFT) or there is no such necessity? Does it depend on the ...
3
votes
0answers
39 views

What's the most complete basis set for DFT?

I'm trying to benchmark some NMR chemical shifts. The literature suggests that the pcS-4 basis set will be the most accurate. I was wondering if there were any basis sets that are more complete than ...
5
votes
1answer
280 views

Hartree Fock and DFT - Exchange term

this might be a stupid question but I am really confused at the moment. I'm currently preparing for an exam in computational chemistry. In the lecture I noted always down that the exchange term (Pauli-...
2
votes
1answer
130 views

How to know about radical character of a diradical molecule?

I understand that the radical character is quantified by the occupation number of the lowest unoccupied natural orbital obtained from UHF calculations or from CASSCF calculations. My issue is that I ...
14
votes
3answers
525 views

How to define custom force in quantum chemistry packages?

I don't know if this is really applicable but I want to add a custom force (say, $F(x,y,z)$) to some atoms in my geometry optimization. Is there any way to add such forces to the DFT or molecular ...
4
votes
3answers
76 views

Reference states for molecular orbital energies in DFT calculations

When doing a population analysis using density functional theory (specifically as implemented in Gaussian), what are the orbital energies (i.e. alpha and beta eigenvalues) with respect to? If I have ...
11
votes
1answer
147 views

Estimating stability and mean life of anomalous tautomers

A pair of bases in the DNA, say, A-T, have a tautomerized form A*-T* (resulting from switching the sides of both protons along the hydrogen bridges). I have studied how, by means of DFT computations, ...
10
votes
3answers
4k views

PBE vs. PBEPBE functional

I've found good explanation what the PBE functional is. However, sometimes I see PBEPBE and I cannot find any good information about what this means and how it's different from PBE. The Gaussian ...
19
votes
2answers
432 views

Do modern dispersion-corrected DFT methods give more accurate molecular geometries?

I was discussing with a colleague the use of modern dispersion-corrected density functionals. I take it almost as a given that the methods generally produce "better" (for some definition) geometries, ...
3
votes
0answers
47 views

Why can we know a materials properties by knowing its ground state at 0K - DFT

I have had this doubt ever since I was introduced to DFT. For what I gathered whatever results you obtained with DFT (by using VASP, quantum expresso or any other software) you can say those ...
9
votes
2answers
134 views

Strange “far out” orbitals in HF and DFT calculations. Why do they occur?

I calculated some orbitals and I get quite strange results with orbital geometries I wouldn't expect from MO theory. A more experienced colleague just said: "That's always the case with big basis sets,...
5
votes
0answers
62 views

Ab initio MD simulations in VASP

I'm not sure if this is an appropriate question for Chemistry.SE, but I'm taking a stab. I'm running ab initio MD calculations in VASP and I'm trying to optimize different input tags. On the VASP ...
2
votes
1answer
43 views

Correlation energy in quantum chemistry

In density functional theory (DFT) there is classic Coulomb energy and non-classic correlation energy. But what is difference between them? In books these energies are explained very similarly. Both ...
1
vote
0answers
47 views

Natural orbital population analysis in Gaussian

Within (restricted) density functional theory and Hartree-Fock, respectively the Kohn-Sham orbitals and spin orbitals are the natural orbitals, with occupation numbers either 2 or 0. As I understood, ...
6
votes
0answers
120 views

Calculation of rate constant for bimolecular reaction from Gaussian DFT calculations

I am working with a bimolecular reaction for which I have computed the energy profile. I have output files for TS and reactants which I intend to use for rate constant calculations. From TST we have: ...
1
vote
0answers
38 views

How is occupation number of lowest unoccupied natural orbital defined?

I have run some calculations to analyse the occupation of natural orbitals. Since it is needed to run in triplet state, I get two orbitals with occupation number 1. Now, which orbital occupation ...
3
votes
1answer
98 views

What is density in density functional theory?

This may be dumb question. Please bear with me. Density functional theory (DFT) is a successful theory for electronic structure calculations of materials. In DFT, electron density is the fundamental ...
0
votes
1answer
54 views

Finding the weighted parameters of the vertex adjacency matrix for a derivative of azulene

In my last question ( Numbering in a derivative of azulene ), I asked about the numbered of a molecule in order to find the Vertex Adjacency Matrix (VAM) of the same molecule. To explain my question ...
1
vote
0answers
17 views

Relationship between the 2nd theorem of Hohenberg and Kohn and the Variational Principle

In Density Functional Theory, the 2nd theorem of Hohenberg and Kohn states: The electron density of the system that minimises the ground-state energy is the exact electron density of the system. ...
2
votes
0answers
59 views

Q-chem equivalent solvent settings for Gaussian NMR ab inito methods

I'm attempting to use Q-Chem 5.0 to replicate CHESHIRE-esque http://cheshirenmr.info/ ab initio calculations of chemical shift values. Unfortunately Gaussian seems to be the predominant tool in this ...
7
votes
1answer
279 views

Hartree product and the Slater determinant

I am very new to theoretical chemistry so my question may sound very basic to most. It is known that the Hartree product does not satisfy the antisymmetry principle for electrons while the Slater ...
25
votes
2answers
854 views

Equivalent of Szabo and Ostlund book for DFT

Szabo and Ostlund book Modern quantum chemistry [1] is extremely useful to understand Hartree-Fock and post-Hartree-Fock methods. Not only it explains the theory behind such methods, but it is also ...
2
votes
0answers
21 views

How to calculate the magnetic exchange coupling in a charged radical complex?

I have a neutral complex with two magnetic centers. For this case it is easy to calculate the exchange coupling as the broken symmetry state is converged properly. But when I am making the complex ...
2
votes
1answer
68 views

Why density functional theory favors a uniform electron density?

Why density functional theory (DFT) favors a uniform electron density when we use it to calculate ground state electronic structure self-consistently? For example, when we use DFT for calculating ...
5
votes
1answer
133 views

Why are the calculated NMR values on the DFT & GIAO level of cyclopentane-1,3-dione so far off?

I'm trying to simulate the 13C NMR spectrum of cyclopentane-1,3-dione (PubChem CID: 77466; CAS 3859-41-4; ChemSpider ID: 69875; SDBS No: 15258; SMILES: C1CC(=O)CC1=O...
2
votes
1answer
50 views

Magnetic couplings using DFT in a spin frustrated system

I’m working with a trimeric copper system that have the classic triangular shape and experiences spin frustration. This is my first time dealing with such system and, so far, from what I’m reading, I ...
1
vote
0answers
29 views

How does MP2 level of theory affect the natural orbital picture of the system?

I have performed CASSCF calculations for a series of molecules with varying input orbitals. I have used UNO, UHF, ROHF and MP2 based natural orbitals to perform the CASSCF calculations. For my system ...
4
votes
1answer
313 views

Defining the Turbomole named BP86 functional in Gaussian [closed]

According to Turbomole, the BP86 density functional has the following definition: $$ \mathrm{Exchange = LDA + B88\\ Correlation = LDA(VWN(V)) + P86} $$ What is the equivalent definition in Gaussian ...
4
votes
1answer
552 views

How to use dispersion correction with TPSSh in Gaussian 16?

A colleague of mine suggested to use TPSSh as a functional. Since I am applying Grimme's D3 dispersion correction by default, I was surprised that it wasn't implemented in the standard density ...
8
votes
1answer
203 views

Ab initio effective potentials

I have been reading a journal article on ab-initio pseudo-potentials, and I need some help understanding it. The article is C. F. Melius and W. A. Goddard, III. Phys. Rev. A 1974, 10, 1528. A summary ...
3
votes
1answer
143 views

How are imaginary frequencies and force constant calculated by DFT?

Computational calculations allow us to simulate the frequencies of molecules. They can even tell us if the optimized structure is a minimum, a saddle point or a maximum according to the number of ...
28
votes
1answer
482 views

d-orbital splittings in WS2 monolayer

The formerly degenerate $d$ orbitals of the tungsten atoms in the $\ce{WS2}$ monolayer are split into three groups: (1) $\mathrm{d}_{z^2}$, (2) $\mathrm{d}_{x^2-y^2}, \mathrm{d}_{xy}$ and (3) $\mathrm{...
2
votes
0answers
31 views

Determining global minimum structures for ion-water clusters

What is the best way to find the global minimum structure(or get close to it) for an ion and N waters? (i.e. ion cluster in vacuum) For example, consider Flouride ion cluster with N=5 waters. I want ...
2
votes
1answer
65 views

DFT calculation of solids with different periods

I would like to simulate a layer of graphene on a solid substrate, say Cu or quartz. I'm interested in their binding energy and other physical or chemical properties. However, their periods do not ...
1
vote
1answer
46 views

Is it possible to define the coordinates of an atom using other coordinates while optimizing in gaussian 16?

I know it is possible to fix the value of composite expressions using generalized internal coordinates (gic) in gaussian 16. But is it possible to assign that value to the coordinate of another atom. ...
3
votes
0answers
84 views

calculate SPE using dummy/ghost atoms in NWChem

I have been making some failed attempts to calculate interaction energy of 2 monomers using a dummy/ghost atom in NWChem6.5. The code I was using to calculate the interaction energy is: ...
13
votes
3answers
676 views

What to do with (large) imaginary frequencies for constrained minimum structures?

I am performing DFT calculations using ORCA 4.0.1 on an enzyme active site model. The model contains 89 atoms including the substrate (see Animation 1), five of which are fixed in space (the spherical ...
26
votes
1answer
2k views

Suggest methods and basis sets for a variety of systems

Please help me with any/all of the cases below. In the following cases, the named method and basis set are not suitable for the chemical systems. Why aren't they? Could you suggest a suitable method/...
2
votes
1answer
400 views

Periodic boundary conditions in DFT for cylindrical system

As I understand it, a DFT calculation uses periodic boundary conditions, whenever you want to simulate a bulk material. This is convenient because all information can then be extracted from the ...
1
vote
0answers
200 views

Understanding the basics of DFT [closed]

I have recently started doing DFT calculations as part of my PhD. The general approach to DFT in my department appears to be that it is regarded as a tool for doing calculations, which need not be ...
3
votes
1answer
166 views

Possible to Calculate Electronic Charge Distributions/Densities for Organic Molecules?

I was wondering if there are any programs that can calculate the electronic charge distributions around an organic molecule, possibly in the form of a volumetric charge volume density? I am reading up ...
1
vote
1answer
52 views

Small but modest displacements at negligible forces from DFT

When looking at a Gaussian log file today, I noticed the following information in the convergence criterion section after an analytical frequency calculation was performed: ...
6
votes
1answer
157 views

Ground state oxidation potential and Ionization potential from density functional theory (DFT)

Can someone please tell me how to calculate the ground state oxidation potential (GSOP)? I understand that according to Koopmans' theorem, the negative of HOMO can be considered as the GSOP. Is there ...
6
votes
2answers
281 views

Is there a reliable way to determine if a chemical system is more suited to be calculated using a multireference method instead of DFT?

Some time ago I recall hearing a professor say that looking at a DFT result (I don't remember exactly what, a TD-DFT, spin contamination etc...), he could tell that the system was multireferenced. My ...
5
votes
1answer
80 views

Atom-condensed softness matrix

In the ACKS2 polarizable force field paper, I found a thing called the atom-condensed softness matrix. In another paper, I found this expression for it: $$ \chi_{kl} = 2 \sum_{i}^{\text{occ MOs}} \...
2
votes
0answers
55 views

What exactly does setting the net charge in a DFT calculation do?

In the implementations of DFT using Gaussian basis sets, it is common to set the net charge. What exactly does this do? The Gaussian (program) manual, for instance, says that it introduces a ...
9
votes
1answer
1k views

Use of basis set in DFT (Density Functional Theory)

Basis sets are used to guess the electronic wave functions for Hartree Fock or similar methods, which are quite legitimate since these methods deal with the wave function of each and every electron. ...
3
votes
1answer
459 views

How to calculate the energy of a single proton? [closed]

Using e.g. the VWN5 functional, what is a simple way to calculate the energy of a single proton? Of course, it is not possible to perform a calculation of a proton in isolation, as it bears no ...
7
votes
1answer
889 views

Calculating dipole moments from atomic partial charges

How can I calculate dipole moment of a complex molecule where the partial charges are available from density functional theory? I have the coordinates of each atom from atomistic simulations of an ...