ERROR: Numerical method failed on all combinations of convergence parameters

[andyql]

hello, I'm trying to run a template with phreeqc, but I'm getting this error: ERROR: Numerical method failed on all combinations of convergence parameters, cell/soln/mix 1.
not sure what is the issue, here is the error output:

SOLUTION_RAW                 1
  -temp                      100
  -pressure                  1
  -potential                 0
  -total_h                   111.0199097274
  -total_o                   55.515392259472
  -cb                        0.0001610658469348
  -density                   0.95837940139191
  -totals
    C(4)                     6.8424205425676e-05
    Cl                       0.00073601424034852
    Na                       0.00081721534164311
    O(0)                     0.0052658559723307
  -pH                        4.1093763816428
  -pe                        14
  -mu                        0.00081740311342841
  -ah2o                      0.99992226406847
  -mass_water                1
  -soln_vol                  1.0435673277967
  -total_alkalinity          6.9382773428519e-05
  -activities
    Alkalinity               -11.791345675074
    C(4)                     -11.791345675074
    Cl                       -3.1606827598409
    E                        -14
    H(0)                     -39.44663977085
    Na                       -3.1040926089278
    O(0)                     -2.5795610192223
  -gammas
    CO3-2                    -0.06586838301126
    Cl-                      -0.016416548492049
    ClO4-                    -0.0164312783792
    H+                       -0.016362057084756
    H2O-0.01                 -1.6450751448853e-06
    HCO3-                    -0.016450751448853
    Na(CO3)-                 -0.016450751448853
    Na+                      -0.016429043763174
    OH-                      -0.01641806284093
GAS_PHASE_RAW                2
  # GAS_PHASE_MODIFY candidate identifiers #
  -type                      1
  -total_p                   1
  -volume                    1
  -component                 CO2(g)
    # GAS_PHASE_MODIFY candidate identifiers #
    -moles                   0.43408664046987
    # GAS_PHASE_MODIFY candidate identifiers with new_def=true #
    -p_read                  0
    # GasComp workspace variables #
    -initial_moles           0
    -p                       0
    -phi                     0
    -f                       0
  -component                 H2O(g)
    # GAS_PHASE_MODIFY candidate identifiers #
    -moles                   5.6943618235537
    # GAS_PHASE_MODIFY candidate identifiers with new_def=true #
    -p_read                  0
    # GasComp workspace variables #
    -initial_moles           0
    -p                       0
    -phi                     0
    -f                       0
  -component                 H2S(g)
    # GAS_PHASE_MODIFY candidate identifiers #
    -moles                   0
    # GAS_PHASE_MODIFY candidate identifiers with new_def=true #
    -p_read                  0
    # GasComp workspace variables #
    -initial_moles           0
    -p                       0
    -phi                     0
    -f                       0
  -component                 N2(g)
    # GAS_PHASE_MODIFY candidate identifiers #
    -moles                   106.05420004798
    # GAS_PHASE_MODIFY candidate identifiers with new_def=true #
    -p_read                  0
    # GasComp workspace variables #
    -initial_moles           0
    -p                       0
    -phi                     0
    -f                       0
  -component                 NH3(g)
    # GAS_PHASE_MODIFY candidate identifiers #
    -moles                   0
    # GAS_PHASE_MODIFY candidate identifiers with new_def=true #
    -p_read                  0
    # GasComp workspace variables #
    -initial_moles           0
    -p                       0
    -phi                     0
    -f                       0
  -component                 O2(g)
    # GAS_PHASE_MODIFY candidate identifiers #
    -moles                   28.516552441627
    # GAS_PHASE_MODIFY candidate identifiers with new_def=true #
    -p_read                  0
    # GasComp workspace variables #
    -initial_moles           0
    -p                       0
    -phi                     0
    -f                       0
  # GAS_PHASE_MODIFY candidate identifiers with new_def=true #
  -new_def                   0
  -solution_equilibria       0
  -n_solution                -999
  -temperature               303.15
  # GasPhase workspace variables #
  -total_moles               0
  -v_m                       0
  -pr_in                     1
  -totals                   

thanks very much!

[dlparkhurst]

Looks like you are trying to put 100 moles of N2(g) in 1 liter of head space. I increased the head space to 10 and equilibrated with the solution to yield a pressure of 400 atm. So, I think 1 L would be something like 4000 atm, which is outside the range of PHREEQC calculations.

[andyql]

thanks very much! it works now!