Questions tagged [basis-sets]
For questions about the usage and construction of electronic structure basis sets.
83
questions
36
votes
3
answers
2k
views
What are the positives and negatives of periodic DFT codes that don't use plane-wave basis sets?
Most periodic density functional theory (DFT) codes use plane-wave basis sets in conjunction with three-dimensional periodic boundary conditions. In contrast, for molecular systems of finite size, ...
- 7,401
31
votes
1
answer
2k
views
Approximate equivalence table between Pople's basis sets and Jensen's DFT optimized polarization consistent basis sets (pcseg-n)
At first I posted this question in chemistry.stackexchange.com, before materials.stackexchange.com was launched in beta. As it went unaswered, I'm giving a try here.
I have some calculations results ...
- 4,128
23
votes
5
answers
793
views
What functions other than Gaussians are used for orbital basis sets?
The Gaussian function $\propto\exp((x-a)^2/b)$ with $b>0$ is one of the most common functions used in molecular modelling (e.g. Gaussian type orbitals).
What are some examples of applications of ...
- 1,627
23
votes
1
answer
223
views
Which basis set families were optimized to give the lowest possible variational energy for a given number of orbitals?
Dunning's basis set families (such as cc-pVXZ and aug-cc-pCVXZ) were not optimized to give the lowest variational energy for a given number of orbitals. For example the S-type orbitals are optimized, ...
21
votes
2
answers
302
views
Why were no correlation consistent basis sets for potassium contracted for a non-relativistic Hamiltonian?
On basis set exchange, the only correlation consistent (cc) basis sets for potassium, are relativistic (labelled by X2C). Very recently (I consider 2017 quite recent, because most Dunning basis sets ...
20
votes
1
answer
454
views
Are plane-wave basis sets reliable for modeling adsorption processes?
Plane-wave basis sets are useful in the calculation of periodic systems and can be used in combination with pseudopotentials. However, when dealing with surfaces the system is no longer periodic in, ...
- 2,731
20
votes
2
answers
393
views
When using the "even tempered method" for augmenting basis sets, what justification is there for the numerical factor which is used?
Diffuse functions are often added to basis sets using "even-tempered" exponents, or sometimes I have heard the phrase "even tempering".
In this paper by Jacek Koput, the author says "the customary ...
19
votes
3
answers
674
views
How trustworthy is the Tran-Blaha modified Becke-Johnson potential and how does it perform across basis sets?
When modeling solid-state materials and particularly semiconductors, one must go beyond LDA and GGA. One alternative is to use hybrid functionals or the $GW$ method. However, this can be very ...
- 4,889
17
votes
3
answers
996
views
How is "basis set projection" done?
I recently came across something called the basis set projection in Q-chem's manual. What it seems to do is to take the converged SCF in a smaller basis set, and then "projects" that onto a ...
- 7,161
17
votes
2
answers
512
views
How does a complete basis set remove Pulay forces/stress?
When finding nuclear forces for the true electronic wavefunction, thanks to the Hellmann-Feynman theorem we only need to consider explicit derivatives with respect to the nuclear coordinates $\mathbf{...
- 15.7k
17
votes
1
answer
3k
views
Gaussian-type orbitals vs Slater-type orbitals in terms of computational efficiency
When creating molecular orbitals via the linear combination of atomic orbitals method I am told that Slater orbitals provide a more accurate representation whilst Gaussian orbitals are computationally ...
- 851
17
votes
2
answers
1k
views
Understanding the difference between Dunning and Karlsruhe basis sets
I have heard that TURBOMOLE works extremely fast (even compared to MOLPRO) with def2 basis sets, but is slow if Dunning basis sets are used.
What is it about def2 basis sets that allows calculations ...
- 407
16
votes
2
answers
595
views
What's the best way to compare two DFT codes?
What's the most efficient way to compare the two DFT codes in terms of computation time, basis sets (for example, one plane wave vs one atomic orbitals), functional performance, etc.?
- 1,051
16
votes
1
answer
324
views
What methods can be used to calculate excited-states in molecular crystals?
In periodic systems plane-wave basis functions are often used. In a previous question on this site it was questioned if plane-wave basis-sets are reliable for modeling adsorption-processes
The ...
- 8,015
15
votes
3
answers
824
views
How should I compare different basis sets for post-Hartree-Fock methods?
HF or DFT methods are variational so I can confidently say that any basis set that gives the lowest energy is the best for that system. However, I have learnt that post-HF methods such as MP2, MP3, ...
- 7,161
15
votes
1
answer
354
views
What are good basis sets and functionals for Ni, Cd and Pb?
I will setup a single point energy calculation using Gaussian for a system containing atoms of B-N-Ni, B-N-Pb and B-N-Cd.
The system are slabs of boron-nitride [1,2] interacting with those metals.
...
- 23.9k
15
votes
1
answer
124
views
Functionals/basis sets optimized with respect to properties?
In general (at least for molecular calculations) basis sets and DFT functionals are fit to some high level calculation or experimental energy. It is speculated that an accurate energy will result in ...
- 15.7k
14
votes
2
answers
431
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,...
- 741
14
votes
2
answers
930
views
When are atomic-orbital-basis (rather than plane-wave) methods appropriate in periodic DFT?
For example, GPAW supports both plane-wave and atomic-orbital basis methods. I know that atomic-orbital basis methods can have difficulty with electrons occupying vacancies for example, but what ...
- 11.4k
14
votes
1
answer
333
views
Is basis set superposition error reduced when using the GAPW method?
CP2K implements the Gaussian and Augmented Planewaves (GAPW) approach for all-electron calculations. My understanding is that the GAPW method involves using atom-centered Gaussian type orbitals to ...
- 950
14
votes
1
answer
148
views
What are the physical consequences of adding a constant to the diagonal of the effective Hamiltonian of monolayer materials?
Effective Hamiltonians modeling many-layered materials are often tuned using some sort of bias voltage. For instance, in a $4\times 4$ Hamiltonian matrix to describe biased bilayer graphene using some ...
- 2,381
13
votes
2
answers
1k
views
What are the scientific justifications of the binding energy equation?
There are many problems where we want to calculate the binding energy between two systems.
Normally, we have a system A, adsorbent (a surface, a nanotube, a protein, etc.), that interact with a system ...
- 23.9k
12
votes
2
answers
160
views
Avoiding charge-migration in electronic structure calculations with solvent-like environments
When performing electronic structure calculations on a molecule, one can surround it with point charges to mimic a solvent environment, and lend polarization effects. There is a bit of a caveat ...
- 4,406
11
votes
2
answers
1k
views
Critical difference and roles of Slater- and Gaussian-type orbitals
I was very interested in the following question and answer about the Slater- and Gaussian-type orbitals (STO and GTO) used in quantum chemistry and a practical calculation.
Gaussian-type orbitals vs ...
- 1,789
11
votes
1
answer
651
views
What is the minimum basis set one should use?
If I'm interested in getting an answer that is minimally correct what is the minimum basis set level I should use? (i.e. anything less would be much less reliable)
Obviously this depends on the ...
- 8,015
11
votes
1
answer
277
views
Predicting linear dependencies in a basis set in advance of calculating any integrals
I've been doing some calculations for a water molecule, with a big uncontracted basis set: For oxygen I am using aug-cc-pV9Z supplemented with the "tight" functions from cc-pCV7Z. The s-type exponents ...
11
votes
1
answer
143
views
Is it possible to make segmented contracted basis sets for correlated calculations?
The segmented contracted basis sets are usually handled well by the QM programs than generally contracted basis sets, as the primitives are not repeated. The correlation consistent basis sets (cc-PVnZ)...
- 7,161
11
votes
0
answers
118
views
Should augmented basis be used for Hirshfeld charges? [closed]
I'm calculating condensed Fukui functions to explain the reactivity of my molecules.
The manual for Multiwfn suggests using atomic dipole moment corrected Hirshfeld charges for it. However I get ...
- 2,383
10
votes
4
answers
2k
views
How can I find the parameters in a basis set?
I have a question about the basis functions of Gaussian. When the basis function is set (e.g., 6-31G), I would like to know the specific value of Gaussian's parameters for each atom.
In particular, I ...
- 1,789
10
votes
2
answers
150
views
How to create a differential spectrum from two simulated spectra, calculated by two different basis sets?
In my bachelor thesis, I compare simulated IR-spectra computed by different basis sets. The simulated spectra from sets like cc-pVTZ and cc-pVQZ are too similar to compare them by plotting them next ...
- 589
10
votes
1
answer
860
views
Molecular orbitals and active space
During the same minute as asking this question, I also asked this at Quantum Computing SE.
In Qiskit, each qubit corresponds to one spin orbital. For example, the $\ce{N2}$ molecule would have 10 ...
- 457
10
votes
1
answer
418
views
How to interpret molden AO/MO coefficients?
A simple calculation with Molpro 2012 on H2 molecule in cc-pvtz basis produces a molden file with 28 molecular orbitals. Each MO is defined with 30 AO/MO coefficient, but the number of symmetric AOs ...
- 103
9
votes
3
answers
1k
views
Which software is good with generally contracted basis sets?
Recently, I learnt that basis sets can be contracted as segmented (like def2 or 6-31G) or general (like cc-PVTZ, or ANO). The general contraction allows all the primitives to appear in all shells, ...
- 7,161
9
votes
3
answers
446
views
The number of orbitals in Gaussian/GAMESS log file
After calculating the energy of a $\ce{CH4}$ molecule with a 6-31G(d,p) basis set, I obtained the Gaussian/GAMESS log file and describe a part of the log file as follows.
...
- 1,789
9
votes
1
answer
384
views
Basis sets optimized for CCSDT(Q) calculations
CCSDT(Q) and CCSDTQ calculations are extraordinarily expensive, both in terms of FLOP count and RAM capacity requirements, eg. CCSDTQ requires a computational effort that scales as ~$N^{10}$, where N ...
- 405
9
votes
2
answers
1k
views
What are the advantages of plane wave basis? [duplicate]
For systems with periodic boundary condition , plane wave basis is usually adopted. While for molecular systems, gaussian basis set is normally adopted. For periodic systems, the gaussian basis is ...
- 3,973
9
votes
1
answer
694
views
How to build cartesian representations of spherical Gaussian basis functions?
I'm trying to code an integral library for my own QC software from scratch. I need to use the spherical basis functions like in any other popular program. So, a D basis function is supposed to make ...
- 133
9
votes
1
answer
316
views
How are molecular orbitals constructed from the coefficients and basis set?
When an SCF calculation is done in a QM program, the MO coefficients of the solution are printed into the output. Now, graphical softwares can plot 3D molecular orbitals from these coefficients, ...
- 7,161
9
votes
2
answers
445
views
Quantifying relativistic effects in DFT calculations
I am wondering what is the best approach to quantify the contribution of relativistic corrections to a result of an electronic structure calculation. As far as I understand, from a user perspective ...
- 624
8
votes
2
answers
580
views
Basis Set to analyze a molecule that involves silver
I am a bit new to computational chemistry, but I've been looking into a reaction of a molecule that contains C, H, and F reacting with AgF in a polar solvent. I am looking into one pathway that ...
- 183
8
votes
2
answers
3k
views
Is the number of basis functions same in the different atoms in the 6-31G basis set?
I'm now studying the quantum chemistry calculation and Gaussian basis set.
For example, when using the 6-31G basis set,
the C atom has 6 basis functions in 1s orbital:
3+1=4 in $2s$,
3+1=4 in $2p_x$,
...
- 1,789
8
votes
1
answer
894
views
Which basis set should be used for a system that contains both anions and cations?
I am trying to model a system that has both anions and cations interacting with each other. (In this case it is calcium ion interacting with two small carboxylic acid molecules). Now, looking through ...
- 7,161
8
votes
1
answer
288
views
Difference between original STO-3G and BSE STO-3G
I am trying to work with the STO-3G basis set, but I noticed that the STO-3G in Basis set exchange is different from the STO-3G originally published, and not only that, the BSE version seem to be ...
- 615
8
votes
1
answer
74
views
PES of Tetramethylthiuramdisulfid / Problem with Disulfidbridge
I like to perform a DFT Study on Tetramethylthiuramdisulfid. So I made a structure optimization for my start structur with PW6B95D3/aug-cc-pvtz to get a good structure to begin my dihedral scan to get ...
- 866
7
votes
2
answers
751
views
p-orbitals in STO-3G basis set
I look at the STO-3G basis set for Carbon atom (downloaded from the basissetexchange.org) and I see one set of exponents and two sets of coefficients for the p-orbitals.
How do these exponents and ...
- 321
7
votes
2
answers
257
views
How does pyscf's internal basis format distinguish between 1s and 2s orbitals?
For example: gto.basis.load("6-31g", "Li") will show:
...
- 231
7
votes
2
answers
197
views
What quantity to perform CBS extrapolation on for adsorption/binding energies?
When calculating complete basis set extrapolations for quantities, such as adsorption energies or binding energies, which depend on the difference between three quantities (i.e. $E_\textrm{prop} = E_\...
- 418
7
votes
1
answer
69
views
Are there properties that can only be calculated with a given set-up?
If you do a bibliography search about the codes used in matter modeling, you will find a large list.
I played with a few of them (both free and commercial) and one of my findings is that the ...
- 23.9k
7
votes
1
answer
54
views
Any molecular codes to calculate non-local susceptibility?
I am interested in calculating the non-interacting non-local susceptibility for simple molecules.
In principle, if you have a single Slater determinant, then you should just take a sum over the ...
user4626
7
votes
1
answer
64
views
Why am I getting wrong information about basis functions on running gaussian calculation using the 6-31G basis set for the krypton atom?
This is the information about basis functions which I am getting on doing Gaussian calculation. But according to the best of my knowledge Kr has six primitive Gaussians for 1s, 2s,2p, 3s, and 3p ...
user139977