Jump to content
Geochemist's Workbench Support Forum

can't retrieve H+ molal in python script


Helge

Recommended Posts

Hello everybody,

I can't retrieve the molal H+ concentration. With

#######

    species = myGWB.results("species")
    conc = myGWB.results("concentration aqueous", "molal")
    print "\nThere are" , len(species) , "aqueous species.\n"
    for i in range(len(species)):
        print "%-4s = %10.4g mol/kg" % (species, conc)
    print "concentration H+ = ", myGWB.results('concentration H+', 'molal')[0]

#######

I get

Y:\Rechnungen\gwb>0001_test.py

There are 9 aqueous species.

H+   =  3.027e-07 mol/kg
HSO4- =  3.231e-06 mol/kg
K+   =     0.7233 mol/kg
KMg(SO4)+ =  6.463e-05 mol/kg
Mg(OH)+ =  3.447e-06 mol/kg
Mg++ =      3.114 mol/kg
Mg3(OH)4++ =  8.826e-15 mol/kg
OH-  =  8.746e-08 mol/kg
SO4-- =      3.475 mol/kg
concentration H+ =  -3.47366172108e-10

It seems that with myGWB.results('concentration H+', 'molal')[0] I retrieve something else (charge balance?) but not the H+ molality. Of course I could work around the problem and get the H+ molality from the aqueous species list,  but this wouldn't be elegant.

Helge

 

 

Link to comment
Share on other sites

Guest Melika Sharifi

Hi Helge,

Concentration H+, molal will give you the concentration of component H+ which refers to the sum of all species with H+, such as HSO4-, OH-, Mg(OH)3, etc. If you want to get the concentration of a particular species, such as H+ which is called free H+, you will need to add "aqueous" to your command (i.e. ('concentration aqueous H+', 'molal').

Hope this helps.

 

Bests,

Melika Sharifi

Aqueous Solutions LLC

Link to comment
Share on other sites

Means: if I want to retrieve the total molality of Tc(+IV) in solution, and TcO(OH)2 is the secondary master for the formation of all Tc(+IV)-species, I use

myGWB.results('concentration TcO(OH)2', 'molal')[0]

?

 

Helge

Link to comment
Share on other sites

Guest Melika Sharifi

Hi Helge,

It will give you the total concentration of Tc in your system (all oxidation states), unless you have decoupled Tc(IV) from other states.

Bests,

Melika

Link to comment
Share on other sites

Taken the case I look at the transition of Tc(+VII) to Tc(+IV) by inreasing pH: is there any means to retrieve the proportion of the two different oxidation states of Tc in solution?

Helge

Link to comment
Share on other sites

Guest Melika Sharifi

Hi Helge, 

I'm not sure which dataset you are using, so I am not able to give you an exact response. It looks like in your script, you have not decoupled Tc(VII) from Tc(IV). If you use "aqueous" in your command, it will give you the free concentration of that specific species in solution. Consider thermo.tdat dataset when Fe(II) is not decoupled from Fe(III),  myGWB.results('concentration aqueous Fe++', 'molal')[0], will give you the concentration of Fe++ aqueous species, and myGWB.results('concentration aqueous Fe+++', 'molal')[0] will give you the concentration of Fe+++ aqueous species. While myGWB.results('concentration Fe++', 'molal')[0], will give you the total concentration of component Fe(II), and myGWB.results('concentration Fe+++', 'molal')[0] will give you the total concentration of component Fe(III) (i.e. all species with Fe(III) such as Fe(OH)3, Fe+++, Fe2O3, etc).

Hope this helps.

Bests,

Melika Sharifi

Aqueous Solutions LLC

Link to comment
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.

Guest
Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

Loading...
×
×
  • Create New...