Comments (7)
I have the same question.
What are the units of the numbers written to RHOION and RHOB?
It seems they are different from CHGCAR. Also it seems that the units in RHOION are different than in RHOB.
Many thanks for any ideas!
from vaspsol.
The units of the data in the RHOIOIN and RHOB files are electrons*grid_volume, same as the units of the CHGCAR file. Based on this, you can convert the units to moles/liter, and we made the plots for the VASPsol paper that way. Note that instead of having absolute ion charge densities or bound charge densities, the RHOION and RHOB files have Net ion/bound charge densities.
Hence, in CHGCAR, the sum(CHGCAR_data/GRIDPTS) == total number of valence electrons.
However, in RHOION, the sum(RHOION_data/GRIDPTS) will be close to zero (order of 1E-10) as the positive and negative ions charge cancel to most parts.
And, the sum(RHOB_data/GRIDPTS) sum would be even smaller (order of 1E-15).
from vaspsol.
Many thanks for the fast reply.
I just still don't understand why sum(RHOION_data/GRIDPTS) is so small.
If RHOION corresponds to the net ion charge density, then shouldn't sum(RHOION_data/GRIDPTS) be equal in magnitude to the net solute charge, so the Delta_NELECT?
As an example, in my calculation I have 308 electrons for the neutral system. Now I change that number to NELECT=308.1.
Then I get sum(CHGCAR_data/GRIDPTS)=308.10000 and sum(RHOION_data/GRIDPTS)=-0.00017924.
But I was expecting the latter to be equal to 0.1 to compensate in magnitude for the net electron number.
As a side note: My cell volume is 1287.7 A^3.
from vaspsol.
UPDATE:
I think I found the missing scaling factors:
The data in RHOION must be multiplied with CELL_VOLUME / LN(10) to obtain the same definition as CHGCAR.
The data in RHOB must be multiplied with 4 * PI to obtain the same definition as CHGCAR.
Is this possible?
from vaspsol.
I wonder if there's any follow up this from the developer @rhennig ? I am trying to visualize the total solvent charge distribution at a charged Cu(100) surface. I also see that, when the metal is charged to +0.1 e, the RHOION integral sum(RHOION/GRIDPTS) gives some -1e-5 e, far from the supposed -0.1 e. I wonder what's missing.
Also, would the total solvent charge simply be RHOB + RHOION? Thanks for your help!
from vaspsol.
Hi, I met a issue. When I opened the RHOION and LRHOB tags, RHOION and RHOB files weren't generated.
I was wondering if there is a solution for it?
I am using vasp 5.4.1 and vaspsol can work well, including lsol, ebk, lambda_d_k.
Thanks
from vaspsol.
Hello @veritas496,
I have not heard of this issue occurring before. Could you provide more detail if possible?
For now all I can recommend is conducting a fresh installation of VASPsol.
Let me know if I can be of further assistance,
Eric
from vaspsol.
Related Issues (20)
- vaspsol setting for the enhanced sampling HOT 3
- Patches for VASP 6.3 HOT 13
- 1. Free Energy or Energy calculated at T=0 K in the OUTCAR calculated using VASPsol. HOT 1
- Potential vs z calculated using VASPsol? HOT 3
- Confusion about DFT energy correction using FERMI_SHIFT HOT 3
- Unable to build VASP 5.4.4 with new Intel OneAPI compilers, mpi mkl HOT 1
- Water molecules at the intersection of implicit and explicit solvents behave abnormally HOT 3
- How to set the parameters in the isolated system of charged adsorbates? HOT 1
- Speed of solvation step with VASPsol HOT 1
- clear
- Fermi energy changes with the value of parameter NC_K
- VASPsol for solid-solid interfaces HOT 1
- How to get the profile of dielectric constant HOT 1
- add topic tags
- requesting pbz_patch_541 patch file
- VASPsol Patch File for VASP 6.4.3 HOT 1
- Will the uniform background charge be turned off after introducing VASPsol for charged slabs? HOT 2
- requesting pbz_patch_5.4.4 patch file
- VASP sol patch file for VASP 6.4 HOT 1
- Issue with optimization with vaspsol with EMIM (EB_K=8.81) HOT 2
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 vaspsol.