Questions tagged [density-functional-theory]
For questions about simulations using Density Functional Theory.
1,134
questions
4
votes
1
answer
66
views
How can I solve the error in TRANSIESTA calculations? [duplicate]
Recently, I've been trying to calculate the transmission coefficients for the molecular junction (Au lattice - molecule - Au lattice) using SIESTA, TRANSIESTA, and TBtrans.
As you know, the procedure ...
15
votes
3
answers
488
views
What is the best program to manipulate plane-wave DFT wavefunctions to calculate custom matrix elements?
For pedagogical reasons, I am looking for ways to calculate quantum-mechanical quantities such as $\langle m | \dot{m}\rangle, \langle m | \dot{n}\rangle, \langle m | \ddot{n}\rangle$ using ...
3
votes
0
answers
43
views
How to calculate band gap alignment for different layers used in double perovskite solar cell?
I want to compare the electronic bands (VBM and CBM) of a double perovskite material with other layer using like transports both holes and electrons .I tried to calculate them, but the values I ...
3
votes
0
answers
32
views
Relax calculation not converging
I am using a DFT+U+V approach and using a pseudohybrid
density functional for extended Hubbard Coulomb interactions along with GBRV ultrasoft pseudopotentials (PBEsol). This enables me to open a gap ...
3
votes
0
answers
9
views
How to build input files to obtain vibrational and thermodynamic measurements using Quntum Espresso?
I'm working with a SiPOCa structure and I've already done the SCF, NSCF, DOS, BANDS, CHARGE DENSITY, FTIR and RAMAN calculations. Now I would like to calculate the vibrational and thermodynamic ...
2
votes
0
answers
13
views
Error in the PH file for RAMAN and FTIR calculation
I am doing the Raman and FTIR calculations using quantum espresso, I ran the following codes
...
4
votes
0
answers
41
views
vc-relax calculation not converging with some pseudopotentials
I am new to Quantum ESPRESSO and it's my first calculations with $\mathrm{CsPbBr}_3$ perovskite structure. I am trying to perform vc-relax calculations. The ...
4
votes
1
answer
56
views
Is it neccessary to apply dipole corrections in TMD bilayers?
I am at the beginning of my PhD and I just learned about the dipole correction in ab-initio calculations. I have seen that it is used mainly when an electric field is applied to a system. But I have ...
3
votes
0
answers
20
views
How to construct the input files to obtain the structural measurements using Quntum Espresso?
I'm working with a SiPOCa structure and I've already done the SCF, NSCF, DOS, BANDS, CHARGE DENSITY, FTIR and RAMAN calculations. Now I would like to calculate the structural measurements of this ...
4
votes
1
answer
64
views
Supercell Elastic Constant
Suppose, that one wants to study a bulk structure which is doped by several impurities. More precisely, I want to calculate the elastic constants using a DFT code (VASP) by fitting stress to strain (...
6
votes
1
answer
258
views
Formation Energy Calculations
How may I calculate each term of the formation energy equation using DFT?
I am assuming that for every software (i.e. Quantum ESPRESSO; Wien2k), the calculation procedure must be the same.
4
votes
1
answer
123
views
Is there an Optimized Norm-Conserving Vanderbilt (ONCV) pseudopotential of Eu atom?
I'm currently searching for an Optimized Norm-Conserving Vanderbilt (ONCV) pseudopotential of Eu, which has a localized f-orbital. But no matter how much I search online, I can't find it.
There is an ...
4
votes
1
answer
40
views
Simple Quantum Espresso-Phonon
I am completely new to material modelling. I am trying to calculate phonon frequency at X point for Aluminium fcc.
I picked an fcc lattice and displaced the atoms at x=0 by (+0.1) and those in x=0.5 ...
4
votes
1
answer
77
views
Help with Geometry Optimization in CASTEP
I am a beginner in CASTEP simulation and I need some guidance on the general steps for geometry optimization. I have tried the following three approaches:
(1) Build unit cell > geometry ...
5
votes
1
answer
170
views
Calculation of material properties
I am in process of learning some simulation codes. There some properties for a material which are input parameters. I want to know if DFT or other related codes can be used to calculate them for any ...
1
vote
0
answers
93
views
Disappearing vacuum layer when relaxing bilayered hBN with VASP
I'm doing geometry optimization of bilayer h-BN, a 72-atom supercell (6x6x1) with vacuum thickness of 20 Å, and I set the parameters in INCAR as below:
...
10
votes
1
answer
2k
views
How to carry out HSE06+SOC band structure calculation using VASP?
I could not find any good source for this type of calculation. I'd be very appreciated if someone who is experienced with that to show it step by step
3
votes
1
answer
49
views
The Dynamic Constant Matrix for the Green Function in the Thermal Conductance Calculation
When computing the thermal transmission coefficient, the green function needs to be computed first; for example, $G^{R}_{scattering\:region}=\left[\left(\omega+i0^{+}\right)^{2}I-D\right]^{-1}$, where ...
2
votes
0
answers
57
views
Prerelaxation before AIMD with VASP
I am investigating a bulk structure by ab initio molecular dynamics (AIMD) implemented in VASP.
I was wondering, how to choose the initial geometry (= POSCAR file) to run the AIMD simulation. In the ...
9
votes
1
answer
435
views
What should be the criteria of convergence over ENCUT?
For example, should a computation be considered to converge when the differences in energy per atom are smaller than a certain value?
2
votes
0
answers
31
views
Dynamical Matrix for the Thermal Transmission Coefficient Calculation
The dynamical matrix for the thermal transmission coefficient calculation is written as follows:
$$D_{I\alpha,J\beta}(\overrightarrow{q})=\frac{1}{\sqrt{M_{I}M_{J}}}\sum_{b}K_{I\alpha,J\beta}(a,b)e^{i\...
4
votes
0
answers
52
views
Struggling with Mixed Basis Set Inclusion for IRC Calculations
I'm trying to compute the IRC for a W complex, and I'm using Def2TZVP for W and 6-311G** for the rest of the atoms. However, when I include the mixed basis set, Gaussian 16 generates a syntax error.
I ...
2
votes
0
answers
60
views
Difference of the Transmission Coefficient between Thermal and Charge Conductance by Nonequilibirum Green Function Method
The equation 57 in the reference [Jian-Sheng Wang, Jian Wang and J. T. Lu, Quantum thermal transport in nanostructures, Eur. Phys. J. B 62, 381 (2008)] explains the the transmission coefficient for ...
2
votes
1
answer
92
views
Errors in CASTEP supercell geometry optimization
These are the errors when I run my supercell on material studio using CASTEP module.
I have tried decreasing k points and also increased the cutoff energy. Furthermore I deleted all the files related ...
4
votes
2
answers
86
views
Confused About Convergence Values in Quantum ESPRESSO: Seeking Help with Unit Conversion
I'm using Quantum ESPRESSO in my calculations. And I'm a bit confused about the convergence for energy between two consecutive steps and the force allowed on each atom.
In my input file, I retained ...
4
votes
1
answer
60
views
Can the Electron Charge Density (ECD) have negative values and why?
In certain charge densities I am getting negative values, in the place where the atomic cores would be. I do not have access to the vasp.xml or OUTCAR files, so there is no way to determine whether ...
4
votes
1
answer
161
views
Is the one-electron reduced density matrix enough in DFT?
There is a conceptual mismatch in my personal understanding between the role of reduced density matrices in DFT and post-HF methods (possibly due to different discussions in textbooks), and I can not ...
4
votes
1
answer
369
views
How to choose the best U value for studying the electronic properties?
I'm new to DFT+U and trying to understand how FeTiO3 works by looking at its electronic properties. But there's a big problem: when I use normal DFT, I can't get the w
right band-gap value.
How can I ...
5
votes
1
answer
95
views
How may I analyze the oxidation state of any atom by analyzing the DOS plot?
I was doing DFT calculations using VASP of some transition metal compounds and plotted the density of states (DOS) along with the projected density of states (PDOS).
I am wondering how can someone ...
4
votes
0
answers
53
views
Source for Hammett Parameters useable in QSAR model training?
I tried to create a QSAR (Quantitative Structure Activity Relationship) model for my thesis that tries to predict the melting points of certain ZIFs (Zeolitic Imidazolate Frameworks) based on relative ...
6
votes
1
answer
92
views
Is it essential to include solvent cage effect in radical reactions?
My question is straightforward and is outlined in the post title: Is it essential to include the solvent cage effect in radical reactions during the computational estimation of Gibbs free energies and ...
5
votes
0
answers
67
views
What do Negative Atomic Orbital Contributions in a Mulliken Population Analysis mean?
I am performing a Mulliken population analysis (MPA) for a molecular system of interest (calculated with B3LYP/Def2-TZVP). I am analyzing the LCAO-MO coefficients/contributions associated with each &...
3
votes
0
answers
29
views
How to choose the correct ibrav value and Identifying the 3-fold Axis
I'm trying to perform relaxation calculations for FeTiO3 with the ilmenite structure in Quantum Espresso. The documentation mentions two possible ibrav values for triagonal systems: 5 and -5. However, ...
2
votes
0
answers
60
views
Temperature effect on electronic properties
I'm particularly interested in studying how temperature affects electronic properties. However, I've been working with Quantum Espresso and its DFT + U functionality, and I'm struggling to understand ...
5
votes
1
answer
255
views
SCF calculations for new materials with unknown properties (e.g. charge, spin, magnetization) using Quantum Espresso
Performing SCF calculations for crystal structures using Quantum Espresso sometimes gives me this error.
...
3
votes
1
answer
77
views
vdW Corrections for Comparing the Effects on Small Molecule Adsorption
Is it necessary to perform a new geometric relaxation with vdW corrections enabled for each implementation (DFT-D2, DFT-D3, etc.) when comparing their effects on the adsorption of small molecules on a ...
6
votes
1
answer
544
views
pyscf: Why should the choice of pseudopotential determine the total energy, and not the physics?
I have limited experience with DFT, but as an exercise I have written a DFT program of an atom (from scratch using nothing more than python built in functions and numpy) by following Richard Martin's ...
2
votes
0
answers
33
views
Help with understanding the method of MBD & MBD-NL
I am having a hard time understanding the steps of Tkatchenko's method of many-body dispersion (MBD) and nonlocal many-body dispersion (MBD-NL).
I would greatly appreciate a detailed & clear ...
2
votes
0
answers
53
views
Frequency calculation using Quantum Espresso
I want to calculate the frequency using quantum espresso. How can I get the input file of frequency calculation using QE?
...
3
votes
1
answer
93
views
How may I verify that my A matrix for TDA (Tamm-Dancoff Approximation of TD-DFT) from PySCF is correct?
This is a follow-up to a series of questions that I have asked here:
Is a D-Wave quantum computer able to do TD-DFT calculations? (answer was "yes, but it might not be efficient")
Inquiry ...
4
votes
0
answers
104
views
Electronic Convergence issues for DFT+U+SOC in VASP for magnetic materials
While performing DFT+U+SOC calculation in VASP for the material $\mathrm{Mn_3Sb}$ my calculation is not at all converging. Here I am attaching the necessary input files. I have come across some ...
7
votes
0
answers
70
views
Degenerate local-orbital radial functions
I am trying to do a simple groundstate calculation with the exciting code for a perovskite structure with 4d-electron elements. The code uses a muffin-tin approach and expands the wavefunctions inside ...
3
votes
2
answers
282
views
Synergy of Exchange and Correlation, Why are they Combined in the Term "Exchange-Correlation" Energy
The exchange interaction is a quantum mechanical effect that only occurs between identical particles. The effect is due to the wave function of indistinguishable particles being subject to exchange ...
2
votes
0
answers
16
views
Can I define molecular units on NBO3.1 or NBO6.0?
I'm studying donor-acceptor charge transfer (CT) process between tetrakis(m2-acetato)-di-Wolframium (acceptor) with an acetonitrile (donor) at the B3LYP-D3BJ/def2-TZVPD level of theory.
The W···N ...
6
votes
2
answers
475
views
What are some pros/cons of Quantum ESPRESSO and FLEUR?
I want to know the pros/cons of plane wave (PW) pseudopotential (PP) based DFT codes and all-electron full-potential linearised augmented-plane wave (FLAPW) DFT codes. In the title, I mentioned ...
6
votes
1
answer
175
views
Modeling techniques for simulating the electronic and structural properties of Transition Metal Dichalcogenide nanotubes
Dear Matter Modeling Stack Exchange community,
I am currently working on calculating the electronic and structural properties of Transition Metal Dichalcogenide (TMD) nanotubes using DFT, with a ...
3
votes
1
answer
86
views
Error in TranSIESTA: "The electrode down-folding region is 0 in the device region"
When I run a tranSIESTA calculation, the SCF converges, and the TBtrans output shows the following error:
...
7
votes
2
answers
117
views
Self-interaction in density functional approximations for many-electron systems
In their famous paper Self-interaction correction to density-functional approximations for many-electron systems Perdew and Zunger said
The exact density functional for the ground-state energy is ...
6
votes
0
answers
81
views
Predicting the band gap or some other properties of materials before running any calculations [closed]
Can we somehow predict what sort of modulation doping will do to a material without actually running any calculations? For instance, could we use the already known density of states (DOS) of Zn to ...
3
votes
1
answer
216
views
About Van der Waals (vdW) correction in VASP
I would like to know when we need to apply Van der Waals (vdW) correction in our system. If I am dealing with 2d systems (monolayer or bilayer), then should I need to add IVDW (vdW (dispersion) ...