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)

Filter by
Sorted by
Tagged with
4
votes
1answer
66 views

Finding mathematically the ground state density in DFT

To find the ground state density in DFT, you set the following Lagrangian: $$L = E[\rho(\vec r)] - \mu\left(\int \rho(\vec r) \mathrm{d}\vec r - N\right)$$ While minimising with respect to the ...
8
votes
1answer
204 views

DFT vs. MP2 for stacked dimer

Recently, I looked at these two papers analyzing the excited-state properties of modified DNA bases (2-aminopurine and 8-vinyl-A) and how they are influenced by stacking with natural nucleobases: ...
0
votes
0answers
32 views

How can I obtain the wave function of the dihydrogen system, by hand using Self Consistent Field DFT calculations? [duplicate]

After doing some energy calculations for the $\ce{H2}$ system in Gaussian, using ωB97X-D/STO-3G level of theory, and generating a .wfn file, I wanted to understand better how the computer does ...
2
votes
0answers
46 views

How to account for solvation of a molecule explicitly in quantum calculation?

I'm looking for a way to account for hydrogen bonds formed with my molecule in aqueous solution explicitly in DFT calculation. I'll do it in order to compare stability of two structures in water. As ...
9
votes
1answer
133 views

How to calculate the amount charge transfer from a natural population analysis?

I am trying to reproduce the quantity referred to as the degree of charge transfer reported in this publication by Zhu et al. [Ref. 1] for the F4-TCNQ molecule. I am interested in section 3.2 in the ...
0
votes
2answers
60 views

How can I find the smaller symmetric structure from big crystal unit cell?

I have Pyrope (Mg3Al2(SiO4)3) crystal structure, downloaded from materialsproject: https://materialsproject.org/materials/mp-6073/# It appeared as triclinic cell, but if you download CIF file and ...
0
votes
0answers
36 views

What is meaning of LAV2P basis set?

I am working on energy optimization of molecular structure using DFT with Schrödinger software and I found that job is running on LAV2P basis set. I want to know what is the meaning of this basis set ...
4
votes
1answer
190 views

Can I compare DFT calculations with different grids?

When doing DFT calculations, some integrations are commonly done numerically on grids. [In fact, more than a single grid may be used at the same time for different integrals, e.g. approximations such ...
1
vote
0answers
23 views

How one can calculate Si-O-Ca/Mg valence angle energy?

I want to perform the relaxed scan (or relaxed surface scan) for the valence angle energy of Si–O–Ca, Si–O–Mg, Al–O–Ca, and Al–O–Mg. I will use DFT (ORCA), probably B3LYP and triple zeta basis set (...
2
votes
0answers
27 views

Guide to DFT functionals [duplicate]

Having recently begun studying Quantum and computational chemistry in depth... I am looking for a manual, text, or guide of any form that could explain the basic formulations and potential of the ...
1
vote
1answer
27 views

Limitations of extremes for the exchange-correlation functional - DFT

In DFT the exchange-correlation function is critical for the correct description of the system. I have a question about the applications and shortcoming of some functionals. LDA delocalizes electrons ...
0
votes
0answers
45 views

Getting radial Kohn-Sham orbital from the output of an electronic structure code using gaussian basis sets

As we know that in solving the Schrodinger equation (Kohn-Sham equation), we expand the wavefunction in terms of some basis sets and them optimize the energy (the expectation value of Hamiltonian) by ...
1
vote
0answers
28 views

Is there an easy explanation of MP2 calculation [closed]

for my master thesis I calculated some ring strains with MP2 calculation. It was not hard to calculate these things but I'm struggling to understand what the theory is about. Now I'm defending my work ...
3
votes
0answers
95 views

DFT Exchange-Correlation Two Electron Integrals

Are there any electronic structure programs that print the exchange-correlation (XC) integrals for a given DFT functional? These integrals are of the form $$\langle pq|w|rs\rangle$$ where $w$ is an XC ...
0
votes
0answers
16 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
2answers
83 views

What's the most complete basis set for calculating chemical shifts with 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 ...
4
votes
3answers
98 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 ...
3
votes
0answers
53 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 ...
5
votes
0answers
85 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 ...
10
votes
2answers
149 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,...
1
vote
0answers
101 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, ...
4
votes
1answer
65 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 ...
6
votes
0answers
323 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
48 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
153 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 ...
1
vote
0answers
20 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. ...
0
votes
1answer
62 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 ...
2
votes
0answers
71 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 ...
6
votes
1answer
355 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 ...
2
votes
0answers
24 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 ...
1
vote
0answers
32 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 ...
2
votes
1answer
52 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 ...
6
votes
1answer
143 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...
4
votes
1answer
1k 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 ...
2
votes
1answer
74 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 ...
3
votes
1answer
168 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 ...
2
votes
1answer
151 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 ...
2
votes
0answers
32 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 ...
10
votes
1answer
140 views

Spin in meta-GGAs

In density functional theory (DFT) we know that the local density approximation (LDA) leads to an exchange-correlation energy kernel of the form: $$E_\mathrm{xc}^\mathrm{LDA} =\int e_\mathrm{xc}\left(...
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
110 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
810 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 ...
2
votes
1answer
74 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 ...
2
votes
1answer
487 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
221 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
196 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
55 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
293 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 ...
2
votes
0answers
73 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. ...