henniggroup / VASPsol

Solvation model for the plane wave DFT code VASP.
Other
138 stars 85 forks source link

Debye length setting when electrolyte containing bivalent ion #38

Closed Julyjing closed 3 years ago

Julyjing commented 5 years ago

According to the Debye length formula, Debye length is related to concentration and valence of ion. May I ask if 1M H2SO4 electrolyte should be expressed as 2M monovalent cation and 1M bivalent anion?

rhennig commented 5 years ago

Dear Julyjing,

For multivalent electrolyte solutions, the Debye length is given by

\kappa\mathrm{D}^{-1} = \sqrt \frac{\epsilon k\mathrm{B} T}{e^2\sum_i{Z_i n_i^b}}

Where Z_i is the valence of the ions in the electrolyte and n_i^b there bulk concentrations. Based on this square root dependence you can determine how the Debye length changes with ionic concentration and their valence.

Best wishes, Richard

— Richard G. Hennig Professor, Materials Science and Engineering, University of Florida Associate Director, Quantum Theory Project, University of Florida Adjunct Professor, Materials Science and Engineering, Cornell University Email: rhennig@ufl.edumailto:rhennig@ufl.edu Web: http://hennig.mse.ufl.eduhttp://hennig.mse.ufl.edu/

On April 14, 2019 at 11:28:23 PM, Julyjing (notifications@github.commailto:notifications@github.com) wrote:

According to the Debye length formula, Debye length is related to concentration and valence of ion. May I ask if 1M H2SO4 electrolyte should be expressed as 2M monovalent cation and 1M bivalent anion?

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_henniggroup_VASPsol_issues_38&d=DwMCaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=bSw49KO0hWxPjfCyM4jz2g&m=EvrtVOaN9f_DTzD0IZJKAprKXO6CVBxnYJAeXkIbdnc&s=VKwKinZN4eWyOU5euGnQRxmfeMaI7tlcmLYRwOStytE&e=, or mute the threadhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_AITp6vKIzEDLTWFBbPQWKpSjayRxp4-5Fyks5vg-2Dq1gaJpZM4cvC2Z&d=DwMCaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=bSw49KO0hWxPjfCyM4jz2g&m=EvrtVOaN9f_DTzD0IZJKAprKXO6CVBxnYJAeXkIbdnc&s=ZiyIR0GpQl11hg4gU234qWzk3ss4tdQCDmTODwf0evw&e=.

xnmao commented 2 years ago

Dear @rhennig,

I have a similar work. I tried to get the the Debye length 3.04 Å of 1M electrolyte following $ \kappa\mathrm{D}^{-1} = \sqrt \frac{\epsilon k\mathrm{B} T}{e^2\sum_i{Z_i n_i^b}} $

epsilon = 78.4
kB = 1.38e-23
T = 298.15
e = 1.60e-19
sumi_Zi_nib = 1*1 + 1*1

kappaD = 1 / (epsilon*kB*T/(e**2*sumi_Zi_nib))**0.5 # 3.9840061403095744e-10, not 3.04

What is wrong with that? Can you fix it? Thank you!

ericfonseca95 commented 2 years ago

Hello Xnmao,

Here is the calculation worked out:

charge_electron = e = 1.6022 E-19 Coulombs T = 298.15 Kelvin K_b = 1.3802 E-23Joules / Kelvin epsilon_0 = 8.85 E-12 Farads/Meter epsilon_r = 78.4 epsilon = epsilon_0 epsilon_r = 6.94 E -10 Debye_Length (Gouy-Chapman theory) = [ (epsilon k_b T ) / ( 2 e^2 z^2 c_i0) ] ^ (1/2)

For a monovalent salt in water, z_i = 1, and c_i0 (bulk concentration of salt) = 1 Mol / Liter = 6.022 E+26 particles/m^3. Plugging in:

Debye_Length = [ (6.94 E -10 1.38 E-23 298 ) / (2 1.602 E-19 ^ 2 1.0^2 * 6.022 E +26)] ^ 1/2 Debye_Length = 3.04 E -10 Meters = 3.04 Angstroms

Hope this helps,

Eric Fonseca

xnmao commented 2 years ago

Cool. It works. Thank you for your help! @ericfonseca95