Example 10--Aragonite-Strontianite Solid Solution

PHREEQC has the capability to model multicomponent ideal and binary nonideal solid solutions. For ideal solid solutions, the activity of each end member solid is equal to its mole fraction. For nonideal solid solutions, the activity of each end member is the product of the mole fraction and an activity coefficient, which is determined from the mole fraction and Guggenheim excess free-energy parameters. Example 10 considers an aragonite (CaCO₃ )-strontianite (SrCO₃ ) solid solution and demonstrates how the composition of the solid solution and the aqueous phase change as strontium carbonate is added to an initially pure calcium carbonate system.

The example is derived from a diagram presented in Glynn and Parkhurst (1992). The equilibrium constants at 25 °C, and , and the Guggenheim parameters, and , are derived from Plummer and Busenberg (1987). The input file is shown in table 29. The PHASES data block defines the log K s for aragonite and strontianite and overrides any data for these minerals that might be present in the database file. The SOLID_SOLUTIONS data block defines the unitless Guggenheim excess free-energy parameters and the initial composition of the solid solution, which is zero moles of aragonite and strontianite. Initial solution 1 is defined to be a calcium bicarbonate solution. The solution is then equilibrated with aragonite at nearly 1 atm partial pressure of carbon dioxide and saved as the new composition of solution 1.

Table 29. Input file for example 10.
TITLE Example 10.--Solid solution of strontianite and aragonite.
PHASES
Strontianite
SrCO3 = CO3-2 + Sr+2
log_k -9.271
Aragonite
CaCO3 = CO3-2 + Ca+2
log_k -8.336
END
SOLID_SOLUTIONS 1
Ca(x)Sr(1-x)CO3
-comp1 Aragonite 0
-comp2 Strontianite 0
-Gugg_nondim 3.43 -1.82
END
SOLUTION 1
-units mmol/kgw
pH 5.93 charge
Ca 3.932
C 7.864
EQUILIBRIUM_PHASES 1
CO2(g) -0.01265 10
Aragonite
SAVE solution 1
END
#
# Total of 0.00001 to 0.005 moles of SrCO3 added
#
USE solution 1
USE solid_solution 1
REACTION 1
SrCO3 1.0
.005 in 500 steps
PRINT
-reset false
-echo true
-user_print true
USER_PRINT
-start
10 sum = (S_S("Strontianite") + S_S("Aragonite"))
20 if sum = 0 THEN GOTO 110
30 xb = S_S("Strontianite")/sum
40 xc = S_S("Aragonite")/sum
50 PRINT "Simulation number: ", SIM_NO
60 PRINT "Reaction step number: ", STEP_NO
70 PRINT "SrCO3 added: ", RXN
80 PRINT "Log Sigma pi: ", LOG10 (ACT("CO3-2") * (ACT("Ca+2") + ACT("Sr+2")))
90 PRINT "XAragonite: ", xc
100 PRINT "XStrontianite: ", xb
110 PRINT "XCa: ", TOT("Ca")/(TOT("Ca") + TOT("Sr"))
120 PRINT "XSr: ", TOT("Sr")/(TOT("Ca") + TOT("Sr"))
130 PRINT "Misc 1: ", MISC1("Ca(x)Sr(1-x)CO3")
140 PRINT "Misc 2: ", MISC2("Ca(x)Sr(1-x)CO3")
-end
SELECTED_OUTPUT
-file ex10.sel
-reset false
-reaction true
USER_PUNCH
-head lg_SigmaPi X_Arag X_Stront X_Ca_aq X_Sr_aq mol_Misc1 mol_Misc2 \
mol_Arag mol_Stront
-start
10 sum = (S_S("Strontianite") + S_S("Aragonite"))
20 if sum = 0 THEN GOTO 60
30 xb = S_S("Strontianite")/(S_S("Strontianite") + S_S("Aragonite"))
40 xc = S_S("Aragonite")/(S_S("Strontianite") + S_S("Aragonite"))
50 REM Sigma Pi
60 PUNCH LOG10(ACT("CO3-2") * (ACT("Ca+2") + ACT("Sr+2")))
70 PUNCH xc # Mole fraction aragonite
80 PUNCH xb # Mole fraction strontianite
90 PUNCH TOT("Ca")/(TOT("Ca") + TOT("Sr")) # Mole aqueous calcium
100 PUNCH TOT("Sr")/(TOT("Ca") + TOT("Sr")) # Mole aqueous strontium
110 x1 = MISC1("Ca(x)Sr(1-x)CO3")
120 x2 = MISC2("Ca(x)Sr(1-x)CO3")
130 if (xb < x1 OR xb > x2) THEN GOTO 250
140 nc = S_S("Aragonite")
150 nb = S_S("Strontianite")
160 mol2 = ((x1 - 1)/x1)*nb + nc
170 mol2 = mol2 / ( ((x1 -1)/x1)*x2 + (1 - x2))
180 mol1 = (nb - mol2*x2)/x1
190 REM # Moles of misc. end members if in gap
200 PUNCH mol1
210 PUNCH mol2
220 GOTO 300
250 REM # Moles of misc. end members if not in gap
260 PUNCH 1e-10
270 PUNCH 1e-10
300 PUNCH S_S("Aragonite") # Moles aragonite
310 PUNCH S_S("Strontianite") # Moles Strontianite
-end
USER_GRAPH Example 10
-headings x_Aragonite x_Srontianite
-chart_title "Aragonite-Strontianite Solid Solution"
-axis_titles "Log(SrCO3 added, in moles)" "Log(Mole fraction of component)"
-axis_scale x_axis -5 1 1 1
-axis_scale y_axis -5 0.1 1 1
-connect_simulations true
-start
10 sum = (S_S("Strontianite") + S_S("Aragonite"))
20 IF sum = 0 THEN GOTO 70
30 xb = S_S("Strontianite")/ sum
40 xc = S_S("Aragonite")/ sum
50 PLOT_XY LOG10(RXN), LOG10(xc), line_w = 2, symbol_size = 0
60 PLOT_XY LOG10(RXN), LOG10(xb), line_w = 2, symbol_size = 0
70 rem
-end
END
#
# Total of 0.005 to 0.1 moles of SrCO3 added
#
USE solution 1
USE solid_solution 1
REACTION 1
SrCO3 1.0
.1 in 20 steps
END
#
# Total of 0.1 to 10 moles of SrCO3 added
#
USE solution 1
USE solid_solution 1
REACTION 1
SrCO3 1.0
10.0 in 100 steps
END

In the next simulation, solution 1 is brought together with the solid solution (USE keywords) and 5 millimoles of strontium carbonate are added in 500 steps (REACTION data block). The PRINT keyword data block excludes all default printing to the output file and includes only the printing defined in the USER_PRINT data block. The USER_PRINT data block specifies that the following information about the solid solution be printed to the output file after each reaction step: the simulation number, reaction-step number, amount of strontium carbonate added, (log of the sum of the ion activity products), mole fractions of strontianite and aragonite, aqueous mole fractions of calcium and strontium, and the composition of the two solids that exist within the miscibility gap. The SELECTED_OUTPUT data block defines the selected-output file to be ex10.sel , cancels any default printing to the selected-output file ( -reset false), and requests that the amount of reaction added at each step (as defined in the REACTION data block) be written to the selected-output file ( -reaction true). The USER_PUNCH data block prints additional columns of information to the selected-output file, including all of the information needed to make figure 11. Two additional simulations add successively larger amounts of strontium carbonate to the system up to a total addition of 10 mol.

The excess free-energy parameters describe a nonideal solid solution that has a miscibility gap. For compositions that fall within the miscibility gap, the activities of calcium and strontium within the aqueous phase remain fixed and are in equilibrium with solids of two compositions, one solid with a strontium mole fraction of 0.0048 and one solid with a strontium mole fraction of 0.8579. For the simulations of example 10, each incremental addition of strontium carbonate increases the mole fraction of strontium carbonate in the solid until about 0.001 mol of strontium carbonate have been added (fig. 11A). That point is the beginning of the miscibility gap (fig. 11) and the composition of the solid is 0.0048 strontium mole fraction. The next increments of strontium carbonate (up to 0.005 mol strontium carbonate added) produce constant mole fractions of calcium and strontium in the solution (fig. 11B) and equilibrium with both the miscibility-gap end members. However, the amounts of calcium carbonate and strontium carbonate in the solid phases (fig. 11C) and the amounts of each of the miscibility gap end members (fig. 11D) vary with the amount of strontium carbonate added. Finally, the end of the miscibility gap is reached after about 0.005 mol of strontium carbonate have been added. At this point, the solution is in equilibrium with a single solid with a strontium mole fraction of 0.8579. Addition of more strontium carbonate increases the mole fractions of strontium in the aqueous phase and in the solid solution until both mole fractions are nearly 1.0 after the addition of 10 mol of strontium carbonate.

| Next || Previous || Top |