Comments (14)
http://www.quantum-simulation.org/potentials/sg15_oncv/ seems to be the most used. The XML format is quite straightforward, so might be the simplest to parse. They give the local potential, the projectors and the dij, so it should be pretty easy to incorporate into DFTK.
from dftk.jl.
Partially done in #741 .
from dftk.jl.
I'm going to call this done for now. We still miss features, but the key things are there.
from dftk.jl.
Notably, fully-relativistic pseudos and pseudos where the local potential angular momentum is defined (i.e. not -1) are not yet supported. Most, if not all, other features of ONCVPSP-generated pseudos (local and non-local potentials, non-linear core correction, pseudo-atomic densities, and pseudo-atomic wavefunctions) are supported.
from dftk.jl.
pseudos where the local potential angular momentum is defined (i.e. not -1)
Is that common? Shouldn't be too hard
fully-relativistic pseudos
This means Dirac equation? Do people really do plane waves with that?
from dftk.jl.
pseudos where the local potential angular momentum is defined (i.e. not -1)
Is that common? Shouldn't be too hard
Not super common, and it shouldn't be very hard to handle.
fully-relativistic pseudos
This means Dirac equation? Do people really do plane waves with that?
Yes, they're from a solution of the Dirac equation. AFAIK, they're needed for spin-orbit coupling.
from dftk.jl.
Yes, they're from a solution of the Dirac equation. AFAIK, they're needed for spin-orbit coupling
OK I see, the PSP code solves the Dirac equation, but the PW code solves for a C^2 valued wavefunction, which has the same Hamiltonian as usual plus a sigma . (nabla V x p) term. I think it's just an info on how the PSP was generated and is meant to be used, but there's nothing more to be done at the PSP level in this case - at least I don't see any relevant field in the UPF format (just PP_SEMILOCAL but pseudos are usually given in nonlocal form anyway).
from dftk.jl.
but there's nothing more to be done at the PSP level in this case
I don't think this is exactly true; the data for the non-local projectors and the pseudo-waves are different, and the way fully-relativistic pseudos are treated in UPF and ABINIT's PSP8 is different enough that they're incompatible without some major finagling.
from dftk.jl.
The data is different but the pseudo-pw code (ie the pseudopotential part of the hamiltonian) is the same? Or what do you mean? If you have references please add them here #526, @ClementineBarat might be interested in implementing this at some point
from dftk.jl.
Just had a chat with a colleague, and it seems that indeed there will be some work to be done for spin-orbit. The relativistic beta and chi projectors need to be used in this case, which can amount to an adjustment to the non-local coupling Dij
or to how the projection vectors are constructed from the radial functions beta
. Seems like there's some Clebsch-Gordan wizardry to be done. (see QE's upflib/spinor.f90
, upflib/upf_spinorb.f90
, Modules/wavefunctions.f90
)
from dftk.jl.
Thanks! I need to read the papers...
from dftk.jl.
At least the ONCVPSP paper uses scalar relativistic pseudos, meaning the pseudos can just be used as is. So let's not worry about it yet...
from dftk.jl.
Ah, possibly not: see discussion eg here https://orbi.uliege.be/bitstream/2268/25694/1/prb_78_045119_verstraete_Pb_phonons_PREPRINT.pdf
from dftk.jl.
Looks like what's going on is that the pseudopotential term itself stays the same, but the spin-orbit coupling term needs the force nabla V. This requires a bit more algebra, but nothing more is required from the pseudos other than the betas and Dij.
from dftk.jl.
Related Issues (20)
- Precompilation causes warning HOT 1
- Rounding issues cause bad wannier interface HOT 4
- unit and equation of calculated force HOT 2
- how to customize/define an xc functional HOT 2
- Multicomponent HOT 4
- Minus sign in Fourier transform of projectors HOT 6
- Unexpected error when using Crystallographic Information Framework (CIF) File Input HOT 5
- Use TestItemRunner.jl for managing tests
- Some characters in `https://docs.dftk.org/dev/developer/symmetries/#Crystal-symmetries` are not rendered HOT 11
- Elastic constants HOT 7
- ExtXYZ segfault HOT 1
- Stresses need an unreasonable amount of memory and runtime
- Space group numbers and symmetries
- Fix phase when possible HOT 2
- Make sure guess_density returns something with correct total charge
- Include `Manifest.toml` in the repo? HOT 2
- Use BenchmarkCI to automatically detect major runtime regressions HOT 3
- Gpu CI Failure HOT 2
- Maybe plot non-converged bands differently HOT 8
- Exact exchange status
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from dftk.jl.