Stack Exchange Network

Stack Exchange network consists of 174 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)

1
vote
0answers
9 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, ...
-1
votes
0answers
19 views

TD-DFT analysis using Gaussian

For the Pentacene molecule, there is an absorption chance in the range of 400nm-650nm that I can’t view by using TD-DFT analysis [![using Gaussian][1]][1] software. I tried different number of states ...
6
votes
0answers
53 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
32 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
72 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
15 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
39 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
36 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
205 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
16 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
21 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
45 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 ...
5
votes
1answer
125 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
217 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
64 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
107 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
89 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
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 ...
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
59 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
555 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
62 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
308 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
177 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
114 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
47 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
2answers
270 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
45 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
857 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
370 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
811 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 ...
0
votes
0answers
19 views

Scan type of DFT calculations adding Freq=Raman keyword [duplicate]

I´m doing a scan calculation by adding Freq=Raman keyword. My objective is to obtain the Raman spectrum of each conformation in my compound. However, when i write the Freq=Raman Keyword, only the ...
10
votes
1answer
181 views

How does one actually get the energy from the Kohn–Sham equations?

This might be a silly question, but how does one acquire the energy of the system from the Kohn–Sham equations? $\left[ -\frac{1}{2}\nabla^2 + V_{eN}(\vec{r}) + V_{ee}(\vec{r}) + V_{xc}(\vec{r}) \...
4
votes
0answers
235 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-...
1
vote
1answer
82 views

How to reconcile conflict between assumed oxidation state and the most stable multiplicity from electronic structure calculations?

Consider either of the small, neutral molecules containing two Zn atoms below. If I were to determine the oxidation state of Zn from this molecule, I would formally assume each Zn to be Zn(I) in order ...
14
votes
1answer
473 views

Determining Kohn-Sham and Hartree Fock virtual orbitals: The underlying field

In Frank Jensen Introduction to Computational Chemistry from 2007, I stumbled upon this paragraph about whether to assign meaning to the KS orbitals or not Another difference is that the unoccupied ...
10
votes
2answers
3k 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 ...
4
votes
1answer
334 views

How to optimize graphene geometry with Gaussian?

I am new to Gaussian and I want to optimize graphene using DFT. After the the geometry optimization it has been severely changed, and some bonds were set as double and triple. Why am I unable to get ...
1
vote
0answers
53 views

How do I get a value of the Gibbs energy of formation at 0 K and 0 bar?

The standard Gibbs free energy of formation for copper(II) oxide, $\ce{CuO}$, is $\Delta_\mathrm fG (\pu{298.15 K},\pu{1 bar}) = \pu{-129.7 kJ mol-1}$. How can I get $\Delta_\mathrm fG (\pu{0 K},\pu{...
2
votes
0answers
61 views

DFT: Can the calculated enthalpy of two systems that aren't isoelectronic be compared?

As a simple example, say I wish to calculate the product of a reaction between methyl iodide and triethylamine, then compare it with the product of benzyl iodide and triethylamine. My goal is to ...
6
votes
0answers
181 views

How to define the initial guess for the electron density?

This is a technical question. If I have a one-dimensional box of length L, and an electron density existing in the box and integrating to N electrons. And if I want to determine the electron density ...
11
votes
1answer
510 views

When is it necessary to check wavefunction stability in density functional theory?

In Gaussian, there is a stable keyword that checks the stability of the wavefunction. Using stable=opt reoptimizes the ...
6
votes
1answer
132 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 ...
4
votes
1answer
579 views

Is there a relation between transition density and density differences?

When I have an excited state, is there a relation between the transition density associated with this electronic transition and the density difference between the excited state density and the ground ...
3
votes
0answers
295 views

Units on rate constants from transition state theory and DFT

I am using density functional theory (DFT) to compute the energy profile and, by extension, kinetic parameters for a catalytic process. I know that the rate constant from transition state theory (TST) ...
4
votes
1answer
142 views

What is the “curse of dimensionality” in molecular dynamics?

The 'curse of dimensionality' or the 'bottleneck problem' in molecular dynamics is explained in page 5 and 6 of Ab Initio Molecular Dynamics: Basic Theory and Advanced Methods by Dominic Marx and Jurg ...
5
votes
2answers
817 views

Partial Atomic Charge in GAMESS

In GAMESS, a calculation of optimisation gives data about partial atomic charges (like Mulliken or Lowdin charges). Can you please let me know how to calculate the other types of partial atomic ...
4
votes
1answer
85 views

Do dipole dipole interactions need to be added explicitly to vdW forces while doing DFT simulations for organic molecules?

I am modelling the adsorption of organic molecules on metal surfaces using plane-wave dft. While there is literature on how important van der Waals forces are for organic systems, and there are also ...
1
vote
1answer
55 views

Building a solid-liquid interface

I am trying to build a solid-liquid interface to be simulate using molecular dynamics. Since, I would be using ab initio molecular dynamics (cpmd or cp2k), I need the minimum cell dimensions. Please ...
4
votes
1answer
272 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 ...