Questions tagged [density-functional-theory]
For questions about simulations using Density Functional Theory.
1,045
questions
3
votes
0
answers
16
views
Carbon Capture Research - Force Constants & Basis Sets
Greetings Matter Modeling,
I'm researching carbon capture dynamics in the context of sorbent regeneration.
The UN IPCC Net-Zero goal requires 75 Mt-CO₂/year captured by 2030. According to the IEA, we ...
0
votes
0
answers
18
views
What (executable) command to calculate charge density using quantum espresso?
I'm trying to calculate the charge density of a structure using quantum espresso, I've already used the DENS.X and DEN.X commands, but in both tests the executable not found error appears. Can someone ...
1
vote
0
answers
13
views
Free energy calculation
I want to plot the relationship between free energy and temperature. I want to consider two plots, when phonons are considered and when phonons are not considered how the free energy changes with ...
2
votes
0
answers
42
views
Wrong lattice parameter of CeO2 after structure relaxation with VASP
I am calculating the lattice parameter of CeO2 with VASP, but the optimized lattice parameter is wrong according to many references (Physical Review B 75, 035109 (2007); The Journal of Chemical ...
2
votes
0
answers
16
views
Strange change in total energy when lattice relaxing of EuTiO3 material and then atomic relaxing in a specific direction using Quantum Espresso
These days, I'm having a trouble in atomic relaxing of "EuTiO3" with tetragonal symmetry using Quantum espresso 7.1 for a period of more that a week. This material have a transverse optical ...
1
vote
0
answers
15
views
Finding optical properties (epsilon) by Quantum ESPRESSO
I am calculating the dielectric ( Epsilon ) for ZnO by quantum espresso, when I run the eps.in file, I always get this error. ' Reading namelist energy_grid'
How do I solve this problem.
4
votes
0
answers
18
views
QE - transition state SCF not converging
Background
Using Quantum Espresso examining a porous material, that has a ligand with an aromatic group. The aim is to rotate the aromatic group and determine an energy profile for the rotation. The ...
2
votes
0
answers
18
views
Creating nanocomposite with material studio
I need to create ZnO/ graphene sheet nanocomposite with material studio and i dont know how. Please, help me or guide me to some one who can help me. I'm blocked.
5
votes
1
answer
71
views
Bader charge calculation in Quantum ESPRESSO DFT
I am searching for calculating Bader charge using Quantum ESPRESSO. I am using pseudopotential based DFT calculation. And I have charge density files (cube and XSF formats).
My objective is to ...
3
votes
0
answers
27
views
Constrained geometry optmization: Molecule interacts with its surroundings
I have a water cluster(suppose) generated from MD run. Now, I want to optimize the geometry of only one water molecule of interest using QC methods (e.g. B3LYP/6-31+G(d,p)). I intend that the water ...
3
votes
0
answers
24
views
How to perform atomistic spin dynamics (ASD) Ab Initio Molecular Dynamics (AIMD) using VASP?
I have tried normal AIMD simulation using VASP for Non magnetic materials. Now, I am interested to perform the molecular dynamics simulations for the finite temperature for Magnetic materials to ...
4
votes
1
answer
96
views
Spin-orbit coupling and lattice constant
Cross-posted on Chem.SE.
Is there any evidence that spin-orbit coupling (SOC) leads to a crystal lattice distortion? Using Quantum Espresso and DFT, I am getting different relaxed lattice constants ...
4
votes
2
answers
230
views
Can the the relaxed energy of a material system (eV/atom) computed by DFT using VASP be positive?
If it is positive, why is it so? Is it due to specific convergence criteria or choice of pseudopotentials? I cannot check the convergence of the DFT calculations as they have been downloaded from a ...
4
votes
1
answer
46
views
Is it possible to find H-bond enrgy with O-H distance less than 0.95A?
I am trying to calculate the energy of a hydrogen bond between RCOO- and Tyrosine (figure below). I freeze the O-O bond at 2.8A and calculate the energy by changing the distance between O-H. But when ...
4
votes
0
answers
56
views
DFT+U: localizing electrons on combination of atomic orbitals
TLDR: DFT+U is a cheap scheme to help localize electrons on atomic orbitals by penalizing fractional occupancies of such orbitals. However, when one tries to localize an electron/hole on an atomic ...
2
votes
0
answers
20
views
Calculation of formation energy of a nanotube using DFT
Suppose that we have a bulk structure of a two-dimensional (2D) material,.e.g graphite, and we want to construct a single-wall carbon nanotube tube starting from this bulk structure.
If we already ...
3
votes
1
answer
46
views
Adsorption of Li on GeC (optimization problem)
Actually I am trying to study the adsorption of Li or K on the sheet starting from the hexagon centre the system optimises with Li or (K) adsorbed over one of the Ge atoms and not remain at the centre....
4
votes
1
answer
79
views
Determining the Nature of a Material's Band Structure: Metallic or Semi-Metallic?
I have obtained a band structure plot (see image below) for a crystalline material. The Fermi level is set at the zero energy level. Based on the band structure, I am uncertain whether the material ...
2
votes
0
answers
46
views
Utilizing DFT-Computed Band Structures in Solid Physics to Investigate Superconductivity: Clues and Indicators in Band Structures
How can band structures and Fermi surfaces calculated with Density Functional Theory (DFT) be effectively employed in the study of superconductivity?
What are the key features or indications related ...
4
votes
1
answer
572
views
Fermi Energy in Band Structure
I am new to this field and after going through some tutorials I was making band structure diagram of Silicon as shown in figure below, and the Fermi energy I got from SCF calculation is 6.208 eV.
...
3
votes
1
answer
47
views
How to read UPF pseudopotential for ultrasoft pseudopotential simple in a fortran program
I'm currently working on calculating semiempirical local and nonlocal matrix elements for Ga. To proceed with my calculations, I need to extract specific parameters from the UPF (Unified ...
3
votes
0
answers
25
views
Discrepancy between Force Convergence in Geometrical Optimization and SCF Calculation in VASP: Reasons and Troubleshooting Suggestions
During a geometrical optimization of ions and lattice parameters in VASP, the forces acting on atoms were minimized to less than 0.01 eV/Å. However, when performing the subsequent self-consistent ...
5
votes
1
answer
87
views
How to calculate the charge density using Quantum Espresso?
How do I modify the input file below to calculate the charge density (DENS) using Quantum Espresso?
Input file:
...
5
votes
0
answers
68
views
Can the electron charge density of two material systems be similar?
Broadly given two electron charge density (ECD) fields, is it possible to distinguish between the elements present in both the systems?
Since the ECD does not explicitly present any information on the ...
3
votes
1
answer
73
views
How to set the LMAXMIX tag in VASP for the material BaLaMn2O6?
I am trying to do DFT+U calculation for the material BaLaMn2O6 . As per the vaspwiki it is recommended to set LMAXMIX = 4 (d-elements) and ...
2
votes
1
answer
52
views
Optimal choice between ISIF = 3 and ISIF = 4 for lattice parameter and atomic position relaxation in VASP calculations
In my VASP calculations, I am interested in simultaneously optimizing both the lattice parameters and atomic positions of my system.
I understand that using ISIF = 3 allows for relaxation of both ...
3
votes
0
answers
29
views
How do I calculate XPS spectra for molecules to match experimental spectra?
I am trying to calculate the XPS spectra using the ICORELEVEL method in VASP, both using initial and final state approaches, for C1s and F1s electrons. My systems are single molecules.
However, my ...
5
votes
1
answer
102
views
What do we get by integrating over the charge density data of a material system calculated with DFT?
I am confused as to whether integrating over the charge density gets us the number of valence electrons in the material system or the total number of valence electrons in the material system.
Let's ...
4
votes
1
answer
146
views
How may I calculate atomic orbital energies using VASP?
I am trying to calculate the d orbital energy of the Cu single atom, like in this page, I have tried to use the Eigenband energies minus the vacuum energy level(about 0.007eV).
...
3
votes
0
answers
44
views
My INCAR file, POSCAR, and KPOINTS are same for GGA, however after the PBE calculation, I found that the bandgap is Changing. Why is that happened?
INCAR, POSCAR, and KPOINTS are the same but for the same calculation, the bandgap is changing.
First, I relaxed the structure then did the self-consistent calculation, and bandstructure calculation. ...
5
votes
1
answer
122
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
48
views
Does VASP include energy penalty contribution in OUTCAR energy of constrained moment calculation?
In VASP, one can use I_CONSTRAINED_M to constrain local magnetic moments to a certain direction in SOC calculations. However, this comes at the cost of adding an energy penalty to the total energy. ...
3
votes
0
answers
20
views
How to calculate kslice_b2 and kslice_b1 in wannierization?
How is reduced coordinates calculated in wannier90? examples 17- 18 , 29
3
votes
0
answers
52
views
Are then any material property solvers (from Charge Density) available in python?
I am looking to predict properties like the Bulk Modulus, total energy, etc., all from the charge density I have as input (CHGCAR file)
7
votes
1
answer
88
views
What is the difference between non self consisent calculation and single-shot calculation?
I am new to density functional theory. My goal is to calculate magnetic anisotropy energy using the force theorem. For that, I have been told that one can run a self-consistent calculation with some ...
2
votes
0
answers
56
views
Clarifying local magnetic moment axes in VASP: cartesian vs crystal (for MAGMOM/M_CONSTR)
In VASP, one can set the local magnetic moments of each ion in a non-collinear SOC calculation along a certain direction using the MAGMOM/M_CONSTR tags. The manual linked to states that these local ...
5
votes
0
answers
38
views
How to calculate anisotropic exchange J $J_{xy}, J_{xz}, J_{yz}, J_{yx}, ...$ using 4-state mapping and DFT?
In this review, the authors present several schemes to calculate various quantities like single-ion anisotropy, DMI, and anisotropic exchange J using DFT. For instance, the DMI term in the y-direction ...
4
votes
1
answer
73
views
computing cholesky errer in nscf calculations?
after running nscf calculations I faced this problem and I don't know what could be the cause of it.
...
6
votes
0
answers
64
views
SCF calculations using HSE hybrid functional dying/stopping? [closed]
I am trying to perform an SCF calculation using HSE hybrid functional to measure the band gap of a 2D semiconductor, but my calculation reached the level of calculating the exact exchange and ...
4
votes
0
answers
34
views
supercell nscf kpoints for optical properties
I am calculating the optical properties of 200 atoms supercell with quantum espresso. What is the approximate number of the nscf kpoints (nk nk nk)should I start with, to avoid lengthy kpoints ...
4
votes
0
answers
58
views
Explicit calculation of Poisson's ratio using energies from a potential energy surface [closed]
I would like to calculate the above quantities for a given unit cell (or supercell) with arbitrary symmetry. Specifically, I am trying to use a machine-learned interatomic potential (MLIP), and since ...
3
votes
0
answers
25
views
Convex hull calculations with WIEN2k
How does one do the calculation of a convex hull with the WIEN2k code for the formation energy of alloys?
5
votes
1
answer
112
views
Contents of a CHGCAR file and how to read it appropriately in Python? [duplicate]
Rephrasing the above, what is the content of the data in the CHGCAR file?
Is the data part shown below representative of the electron-charge density or is it the actual charge in the system and I have ...
4
votes
0
answers
105
views
Is DFT+U an ab initio based method? [duplicate]
As the title already states is DFT+U considered an ab initio method?
6
votes
1
answer
159
views
Fcc lattice with two atom basis equivalent to tetragonal lattice?
I am trying to follow this tutorial from the FLEUR DFT code's official website. The tutorial goes over calculating the ground state energies for the ferromagnetic and antiferromagnetic configurations ...
5
votes
0
answers
62
views
VESTA using command line
Can anyone please show me how to make a supercell with VESTA
using a single command line without opening a GUI.
4
votes
1
answer
61
views
How to find the Fermi energy with highest occupied and lowest unoccupied level from a single nscf caculation using quantum espresso?
Is it possible to find the Fermi energy with highest occupied and lowest unoccupied level from a single nscf caculation using quantum espresso?
5
votes
1
answer
84
views
Is it necessary to use DFT+U when calculating binding energy?
Recently I am trying to calculate the binding energy between CO2/ C2H2 and MOF (containing Ni, C, H, O, N).
I want to know that do I have to use DFT+U for Ni element in this binding energy calculation?...
4
votes
0
answers
50
views
How to Install GIPAW on Ubuntu with QE installation from Ubuntu Repositories
I've just installed Quantum Espresso on Ubuntu using the command:
sudo apt install quantum-espresso
everything went well and I can run calculations on QE normally. ...
2
votes
0
answers
58
views
When calculating DOS and band structures, do I still need to provide a regular mesh for MetaGGA functionals as in the case with hybrid functionals?
In VASP manual it is mentioned that the deorbitalized versions of SCAN and R2SCAN MetaGGA functionals do not depend on the kinetic energy density but instead on the Laplacian of the charge density.
...