BioRisk 17: 389-406 (2022) Apeer-rev ‘iewed open-access journal 1] doi: 10.3897/biorisk.1 7.77487 RESEARCH ARTICLE & R te) Re IS https://biorisk.pensoft.net Development of accurate chemical thermodynamic database for geochemical storage of nuclear waste. Part Il: Models for predicting solution properties and solid-liquid equilibrium in binary nitrate systems Stanislav Donchev!, Tsvetan Tsenov!, Christomir Christov! | Department Chemistry, Faculty of Natural Sciences, Shumen University “Konstantin Preslavski”, Shu- men, Bulgaria Corresponding author: Christomir Christov (ch.christov@shu.bg) Academiceditor: Michaela Beltcheva | Received 2 November 2021 | Accepted 13 December2021 | Published21 April2022 Citation: Donchev S, Tsenov T, Christov C (2022) Development of accurate chemical thermodynamic database for geochemical storage of nuclear waste. Part II: Models for predicting solution properties and solid-liquid equilibrium in binary nitrate systems. In: Chankova S, Peneva V, Metcheva R, Beltcheva M, Vassilev K, Radeva G, Danova K (Eds) Current trends of ecology. BioRisk 17: 389-406. https://doi.org/10.3897/biorisk. 17.77487 Abstract The main purpose of this study is to develop new thermodynamic models for solution behavior and solid-liquid equilibrium in 10 nitrate binary systems of the type 2-1 (Mg(NO,),-H,O, Ca(NO,).- H,O, Ba(NO,),-H,O, Sr(NO,),-H,O, and UO,(NO,),-H,O), 3-1 (Cr(NO,),-H,O, Al(NO,),-H,O, La(NO,),-H,O, Lu(NO,),-H,O), and 4-1 (Th(NO,),-H,O) from low to very high concentration at 25 °C. To construct models, we used different versions of standard molality-based Pitzer approach. To param- eterize models, we used all available raw experimental osmotic coefficients data (p) for whole concentra- tion range of solutions, and up to supersaturation zone. The predictions of developed models are in excel- lent agreement with ¢-data, and with recommendations on activity coefficients (y,) in binary solutions from low to very high concentration. The Deliquescence Relative Humidity (DRH), and thermodynamic solubility product (as In K°,) of 12 nitrate solid phases, precipitating from saturated binary solutions have been calculated. The concentration-independent models for nitrate systems described in this study are of high importance for development of strategies and programs for nuclear waste geochemical storage. Keywords Nuclear waste sequestration, Chemical modelling,Pitzer approach, DRH and K°, of Mg(NO,),.6H,O\(), Ca(NO,),.4H,O(s), Ca(NO,),.3H,O(s), — Ba(NO,),(s),_- Sr(NO,),(s), 9 UO,(NO,),.6H,O(), AI(NO,),.9H,O(s), Cr(NO,),(s), La(NO,),.6H,O (s), La(NO,),(s), Lu(NO,),.5H,O(s) and Th(NO,),.6H,O(s) Copyright Stanislav Donchev et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 390 Stanislav Donchev et al. / BioRisk 17: 389-406 (2022) Introduction Computer models that predict solution behavior and solid-liquid-gas equilibria close to experimental accuracy have wide applicability. They can simulate the complex changes that occur in nature and can replicate conditions that are difficult or expensive to duplicate in the laboratory. Such models can be powerful predictive and interpretive tools to study the geochemistry of natural waters and mineral deposits, solve environmental problems and optimize industrial processes. However, development of comprehensive models for natural systems, with their complexity and sensitivity, is a very difficult, time consuming and challenging task. The specific interaction approach for describing electrolyte solutions to high concentration introduced by Pitzer (1973, 1991) represents a significant advance in physical chemistry that has facilitated the construction of accurate computer thermo- dynamic models. It was showed that this approach could be expanded to accurately cal- culate solubilities in complex brines, and to predict the behavior of natural and industrial fluids from very low to very high concentration at standard temperature of 25 °C (Harvie et al. 1984; Trendafelov et al. 1995a, 1995b; Christov 1996a, 1998, 1999, 2001, 2002a, 2002b, 2003a, 2003b, 2003c, 2004, 2005; Christov et al. 1998; Ojkova et al. 1999; Park et al. 2009; Kolev et al. 2013; Lach et al. 2018; Guignot et al. 2019; Donchev and Chris- tov 2020; Lassin et al. 2020; Donchev et al. 2021; Tsenov et al. 2021), and from 0 to 290 °C (Petrenko and Pitzer 1997; Christov and Moller 2004ab; Moller et al. 2006, 2007; Lassin et al. 2015; Christov 1995, 1996b, 2005, 2007, 2009, 2012, 2020). A long term safety assessment of a repository for radioactive waste requires evi- dence, that all relevant processes are known and understood, which might have a sig- nificant positive or negative impact on its safety. It has to be demonstrated, that the initiated chemical reactions don’t lead to an un-due release of radionuclides into the environmental geo-, hydro-, and bio-sphere. One key parameter to assess the propaga- tion of a radionuclide is its solubility in solutions interacting with the waste. Solubility estimations can either be based on experimental data determined at conditions close to those in the repository or on thermodynamic calculations. A so called “thermody- namic database” created from experimental data is the basis for thermodynamic model calculations. Since the disposal of radioactive waste is a task encompassing decades, the database is projected to operate on a long-term basis. Chemical models that predict equilibrium involving mineral, gas and aqueous phases over a broad range of solution compositions and temperatures are useful for studying the interactions between used nuclear fuel waste and its surroundings. ‘The reliability of such predictions depends largely on the thermodynamic database. Waters of high salinity are not a typical of many geochemical environments which may be chosen as future nuclear waste reposi- tory sites. This suggests that an accurate description of highly saline waters should be required for modeling of chemical interactions in and around nuclear repositories. Currently, the most accurate description of saline waters uses the Pitzer ion interac- tion model. Extensive thermodynamic databases, which are based on the Pitzer ion interaction model was developed within the Yucca Mountain Project (YMTDB: data0. ypf.r2) (Sandia National Laboratories 2007), and Thereda project (THermodynamic REference DAtabase, THEREDA-Final Report) (Altmaier et al. 2011). Unfortunately, Models for Nitrate Systems Seal many of introduced in YMTDB and in THEREDA databases Pitzer models are con- centration restricted and cannot describe correctly the solid-liquid equilibrium in geo- chemical and industrial systems of interest for nuclear waste programs. Nitrates are expected to play a significant role in the context of the underground ge- ochemical repository of nuclear waste (Lach et al. 2018; Guignot et al. 2019; Donchev and Christov 2020; Lassin et al. 2020). More precisely, long-lived, intermediate-level radioactive wastes that are planned to be stored in deep clay formations are composed of dried sludge from effluent treatments that contain significant quantities of nitrate amongst other elements. They are enclosed in specific containers which are placed in underground cavities dug in a very low-permeable argillite host rock. “The storage safe- ty analyses show that, despite the protection of the concrete or stainless steel-made ex- ternal layers of the containers, the formation water of the host rock is likely to migrate and reach the waste during the disposal period” (Donchev and Christov 2020; Lassin et al. 2020). This would result in the potential dissolution of large amounts of nitrate and other elements, resulting in a highly saline, corrosive and oxidative media with a high reactivity towards the containment materials and their surroundings, including the host rock. Several options for the management of radioactive waste involve their preliminary leaching using nitric acid or nitrate salts to recover U and Pu, followed by the incorporation of the leaching residues into concrete or metal packages, which are then stored underground. For safety analysis purpose, the scenarios envisaged for these options assume that formation water returns to the storage compartments sometime after the end of the operating period Dossier ANDRA (2005). With a natural pH value slightly above 7 in clayey formations, pore water has to flow through basic concrete ma- terials before being in contact with the acidic nuclear waste. A large range of pH lead- ing to various chemical behaviors can thus be expected in the vicinity of the waste. At near-neutral to basic pH values, reactions of hydrolysis, complexation, and formation of solid phases can take place and control the fate of radionuclides (Wang et al. 2006; Lassin et al. 2020). Therefore, this reactivity must be characterized by development of not concentration restricted thermodynamic models, which accurately describe not only solution behavior at low molality, but also low and high molality solid-liquid phase equilibrium in nitrate systems (Donchev and Christov 2020). The experimental data presented in Rard et al. (1977, 2004), Rard and Spedding (1981), Maliutin et al. (2020), El Guendouzi and Marouani (2003), and accurate models reported in previ- ous studies (Wang et al. 2006; Lach et al. 2018; Guignot et al. 2019; Donchev and Christov 2020; Lassin et al. 2020), and in the present work is a step towards this objec- tive. It should be noted that THEREDA (Altmaier et al. 2011) do not include models for nitrate solutions and solids. The models introduced in YMTDB (Sandia National Laboratories 2007), including these for NO,-systems are restricted up to 6 mol.kg". In our previous study (Donchev and Christov 2020) we reported very well vali- dated accurate thermodynamic models based on Pitzer ion interactions approach for 7 nitrate binary systems of the type 1-1 (HNO,-H,O, LiNO,-H,O, NaNO,-H,O, KNO,-H,O, RbNO,-H,O, CsNO,-H.O, and NH ,NO,-H,O) from low to very high concentration at 25 °C. In this study we developed very well validated not concentra- tion restricted thermodynamic models for solution behavior and solid-liquid equi- 392 Stanislav Donchev et al. / BioRisk 17: 389-406 (2022) librium in 10 nitrate binary systems of the type 2-1 (Mg(NO,),-H,O, Ca(NO,).- H,O, Ba(NO,),-H,O, Sr(NO,),-H,O, and UO,(NO,),-H,O), 3-1 (Cr(NO,),-H,O, AI(NO,),-H,O, La(NO,),-H,O, Lu(NO,),-H,O), and 4-1 (Th(NO,),-H,O) from low to very high concentration at 25 °C. Models are developed on the basis of Pitzer ion interactions approach. The models for nitrate systems described in this study are of high importance, especially in development of strategies and programs for nuclear waste geochemical storage. These models are also of interest for industrial application, such as production and purification of nitrate compounds. Methodology The models for nitrate binary systems have been developed on the basis of Pitzer’s semi- empirical equations (Pitzer 1973, 1991). Since the Pitzer’s representation of the aque- ous phase is based on the excess free energy, all the activity expressions are consistent, allowing different kinds of data (e.g., osmotic, emf, and solubility measurements) to be used in the parameter evaluations and other thermodynamic functions to be calcu- lated (Christov and Moller 2004a, Christov and Moller 2004b; Christov 2007, 2009, 2012). Pitzer approach has found extensive use in the modeling of the thermodynamic properties of aqueous electrolyte solutions. Several extensive parameter databases have been reported. These include: 25 °C database of Pitzer and Mayorga (1973, 1974) (summarized also in Pitzer 1991); of Kim and Frederick (1988); YMTDB (Sandia National Laboratories 2007), and THEREDA (Altmaier et al. 2011). However, some of the models in all of these databases are concentration restricted, and do not include all minerals precipitating from saturated and supersaturated binary and mixed systems. The most widely used are databases of Chemical Modelling Group at UCSD (Uni- versity California San Diego): at 25 °C (Harvie et al. 1984; Park et al. 2009), and T- variation (from 0 to 300 °C) (Christov and Moller 2004a, Christov and Moller 2004b; Moller et al. 2006, 2007; Christov 2009). Some of comprehensive minerals solubilitiy YMTDB (Sandia National Laboratories 2007), and THEREDA (Altmaier et al. 2011) databases also contain concentration restricted models for some low-, or high- concen- tration binary and mixed sub-systems with strong association reactions in unsaturated solutions. The concentration restricted sub-models are developed using experimental activity data in binary solutions, and solubility data in binary and high order systems up to maximum concentration (m(max)), which is much lower than concentration of saturated or supersaturated binary and mixed solutions (m(sat)). Such a restricted models predict minerals solubility, which is in pure agreement with experimental data. The Pitzer’s equations According to Pitzer theory electrolytes are completely dissociated and in the solution there are only ions interacting with one to another (Pitzer 1973; Pitzer and Mayorga 1973). Two kinds of interactions are observed: (i) specific Coulomb interaction be- Models for Nitrate Systems 393 tween distant ions of different signs, and (ii) nonspecific short-range interaction be- tween two and three ions. The first kind of interaction is described by an equation of the type of the Debye-Hueckel equations. Short-range interactions in a binary system (MX(aq)) are determined by Pitzer using the binary parameters of ionic interactions (8,8, C*, and B®). The Pitzer’s equations (1 to 4) are described and widely discussed in the literature (Harvie et al. 1984; Moller et al. 2006, 2007; Christov and Moller 2004a, Christov and Moller 2004b; Christov 2005). Here only the expression for the activity coefficient of the interaction of cation (M) with other solutes, Your is given: Inv, =2,,F +> m, (2By,(1)+ZCy,)+ > m, [29 LIne | LE MI + [Zu] >) >) m.11,Cog ve m, (24a )+ 2, Dial nal (1) Equation (1) is symmetric for anions. The subscripts c and a in eqn 1 refer to cations and anions, and m is their molality; z is the charge of the M* ion. B and ® represent measurable combinations of the second virial coefficients; C and ) represent measurable combinations of third virial coefficients. B and C are parameterized from single electrolyte data, and @ and ¢) are parameterized from mixed solution data. The function F is the sum of the Debye-Hueckel term, -A? [VI/ (1+ bV I) + (2/b) (n(1 + bv I) ], (2) and terms with the derivatives of the second virial coefficients with respect to ionic strength (see Harvie et al. 1984). In Eq. (2), b is a universal empirical constant assigned to be equal to 1.2. A? (Debye-Hiickel limiting law slope for the osmotic coefficient) is a function of temperature, density and the dielectric constant of water (Christov and Moller 2004b). For the interaction of any cation M and any anion X in a binary system MX-H,O, Pitzer assumes that in Eq. (1) B has the ionic strength dependent form: Bee = BOL er Bie o(a, VI) + (3) +8, g(a,VD, (3A) where g(x) = 2[1 - (1 + x)e*] / x’ with x = aV lor a,V I. x terms are function of elec- trolyte type and does not vary with concentration or temperature. In Eq. 1, the ® terms account for interactions between two ions i and j of like charges. In the expression for ®, ®, = 0, + =O (1), (4) 0,, is the only adjustable parameter. ‘The "0. (I) term accounts for electrostatic unsym- metric mixing effects that depend only on the charges of ions i and j and the total ionic 394 Stanislav Donchev et al. / BioRisk 17: 389-406 (2022) strength. ‘The ¢,,, parameters are used for each triple ion interaction where the ions are not all of the same sign. Their inclusion is generally important for describing solubili- ties in concentrated multicomponent systems. Therefore, according to the basic Pitzer equations, at constant temperature and pressure, the solution model parameters to be evaluated are: 1) pure electrolyte 8, 8“, and C® for each cation-anion pair; 2) mixing 0 for each unlike cation-cation or anion-anion pair; 3) mixing ) for each triple ion interaction where the ions are all not of the same sign. Fluids commonly encountered in natural systems include dissolved neutral species (such as carbon dioxide (CO; 3 SiO 544) and Al(OH),°(aq)). To account neutral specie interactions in aqueous solutions the UCSD Chemical Modelling Group included in their models additional terms to Pitzer equations, denoted as ese or hy Re and ¢ (Eq. (1)) (Harvie et al. 1984; Moller et al. 2006, 2007). N, A,X. The 8” parameter (Eqn. 3A) for 2—2 type of electrolytes Pitzer and Mayorga (1973) did not present analysis for any 2—2 (e.g. MgSO,-H,O) or higher {e.g. 3-2: Al,(SO,),-H,O} electrolytes. Indeed, they found that three 6, 6, and C? parameters approach (see Eqns. 1 and 3) could not accurately fit the activity data for these types of solutions. For these electrolytes mean activity (y,) and osmotic (y) coefficients drop very sharply in dilute solutions, while showing a very gradual increase, with a very wide minimum at intermediate concentration. Pitzer concluded that this behaviour is due to ion association reactions and that the standard approach with three evaluated solution parameters cannot reproduce this behaviour. This lead to a further (Pitzer and Mayorga 1974) modification to the original equations for the description of binary solutions: parameter 8*(M,X), and an associated a,V I term are added to the B,,, expression (see Eqn. (3A). Pitzer presented these parameterizations assuming that the form of the functions (i.e. 3 or 4 8 and C ® values, as well as the values of the « terms) vary with electrolyte type. For binary electrolyte solutions in which either the cationic or anionic species are univalent (e.g. NaCl, Na,SO,, or MgCl), the standard Pitzer approach use 3 parameters (i.e. omit the 8 term) and «, is equal to 2.0. For 2—2 type of electrolytes the model includes the 8 parameter and «, equals to 1.4 and «, equals to 12. This approach provides accurate models for many 2—2 binary sulfate (Pitzer and Mayorga 1974; Christov 1999, 2003a) and selenate (Christov 2003a; Christov et al. 1998) electrolytes, giving excellent representation of activity data covering the entire concentration range from low molality up to saturation and beyond. Inclusion of “standard Pitzer approach” 8” parameter into a models for I-I, 2-1, 3-1, 4-1, 1-2, 1-3, and 3—2 type of electrolytes Some authors found that there are some restrictions limited the potential of the model to describe correctly activity and solubility properties in some binary electrolyte sys- tems with minimum one univalent ion (see Petrenko and Pitzer 1997, 2012; Gruszk- iewicz and Simonson 2005; Lach et al. 2018; Guignot et al. 2019; Lassin et al. 2020), Models for Nitrate Systems Sie) and of 3—2 type (see Christov, 2002ab, 2003b) at very high molality using classical 3 parameters (8, 8", and C*) approach. According to discussion in Christov (2004, 2005, 2012) and in Lassin et al. (2015), there is one major factor which determined these restrictions: type of ¢ (osmotic coefficient) vs. m, or y, (activity coefficient) vs. m dependences at high concentration. For all these systems y vs. m, or y, vs. m curves have a wide maximum at molality approaching molality of saturation: “LiCl(aq) type”: see Lassin et al. (2015), “FeCl,(aq) and FeCl,(aq) type”: see Christov (2003c, 2004); “HNO, (aq) type”: see Donchev and Christov (2020);“Al (SO,),(aq), Cr,(SO,),(aq) type”: see Christov (2002a, 2002b, 2003b). To describe the high concentration solution behaviour of systems showing a “smooth” maximum on y, vs. m dependence, and to account strong association reac- tions at high molality, Christov (1996a, 1998, 1999, 2001, 2005) used a very simple modelling technology: introducing into a model a fourth ion interaction parameter from basic Pitzer theory {8 in Eqn. (3A)}, and varying the values of «, and «, terms in Eqs. (3 and 3A))). The author also found that by variation of the values of «, and a, terms it is possible to vary the concentration range of binary solutions at which as- sociation reactions become more important and should be account by introducing 8” parameter. According to Christov (2005), model which uses «, = 1.4 and a, = 12 ac- counts association only at low molality solutions (see also Christov and Moller (2004b) for Ca(OH),-H,O model). According to previous studies of one of the authors (Chris- tov) an approach with 4 ion interaction parameters (8,8, 8®,and C*), and accept- ing a, = 2, and varying in , values can be used for solutions for which ion association occurs in high molality region. This approach was used for binary electrolyte systems of different type: 1-1 type {such as HNO,-H,O, LiNO,-H,O (Donchev and Chris- tov 2020), CsF -H,O (Tsenov et al. 2021), and LiCl-H,O (Lassin et al. 2015)}, 2-1 {such as NiCL-H,O, CuCl,-H,O, MnCl -H,0O, CoCL,-H,O: (Christov 199G6a, 1999); FeCl,-H,O: (Christov 2004); Ca(NO,),-H,O: (Lach et al. 2018); UO,(NO,),-H20 (Lassin et al. 2020)}, 1-2 {such as Na,Cr,O,-H,O: (Christov 2001); K,Cr,O_-H,O: (Christov 1998)}, 3-1 {such as FeCl,-H,O: (Christov 2004); Ln(NO,),(aq): (Guignot et al. 2019)}, and 3-2 {such as Al (SO,),-H,O, Cr,(SO,),-H,O, and Fe,(SO,),-H,O: (Christov 2001, 2002a, 2003b, 2004, 2005)}. The resulting models reduce the sigma values of fit of experimental activity data, and extend the application range of models for binary systems to the highest molality, close or equal to molality of saturation {m(sat)}, and in case of data availability: up to supersaturation. For example, aque- ous complexes free 4 parameters model for LiCl-H,O system predicts LiCl.nH2O(s) solubilities from 0 to 200°C and up to 40 mol.kg" (Lassin et al. 2015). The resulting accurate 4 - parameters solution models are used directly to determine InK°, values of precipitated solid phases using solubility approach (Harvie et al. 1984; Christov 1995, Christov 1996a, Christov 1996b, Christov 2005, Christov 2012; Christov and Moller 2004a, Christov and Moller 2004b). Therefore, the developed not high molality re- stricted parameterization, were used without any changes for development of solid- liquid equilibrium models for high order systems. ‘Thus, models for Al,(SO,,),(aq) and Cr,(SO,),(aq) are used without additional adjustments to construct a model for multi- 396 Stanislav Donchev et al. / BioRisk 17: 389-406 (2022) component (Na+K+NH,+Mg+Al+Cr+SO ,+H,O) system (Christov, 2002ab, 2003b). Four parameters (8, 6, 8 and C*) models for NiCl,(aq), CuCl,(aq), MnCl (aq), and CoCl,(aq) are used for construction of Na-K-Rb-Cs-Ni-Co-Cu-Mn-Cl-H,O model (Christov 1996a, 1999). Four parameters models for FeCl,(aq)and FeCl, (aq) are directly used in development of high accuracy minerals solubility model for (Na +K+Mg+Fe(II)+Fe(II)+Cl+SO,+H,O) system (Christov 2004). A model for binary systems Na,Cr,O,(aq), K,Cr,O,(aq) was used without any changes to develop a com- prehensive model for: (Nat+K+Cl+SO,+Cr,O,+H,O) system (Christov 1998, 2001), and Ca(OH),(aq) model is used as s strong base for H+Na+K+Ca+OH+Cl+SO ,+H,O model from from 0 to 300°C (Christov and Moller, 2004b). Results and discussions Model parameterization and validation of models for binary 2—1, 3—I,and 4—| type nitrate systems In this study we developed new thermodynamic models for solution behavior and solid-liquid equilibrium in 10 nitrate binary systems of the type 2-1 (Mg(NO,).- H,O, Ca(NO,),-H,O, Ba(NO,),-H,O, Sr(NO,).-H,O, and UO,(NO,),-H,O), 3-1 (Cr(NO,),-H,O, AI(NO,),-H,O, La(NO,),-H,O, Lu(NO,),-H,O), and 4-1 (Th(NO,) ,-H,O) from low to very high concentration at 298.15 K. New sets of Pitzer ion interaction binary parameters are evaluated using available raw experimental osmot- ic coefhicients (p) data for whole molality range of solutions. Rard and co-authors (1977, 1981) reported an extensive experimental activity database for rare earth nitrate systems. Data of Rard and Spedding (1981) are used to parametrize the model for Lu(NO,),- HO system. The 9 vs. m data for remaining 9 nitrate solutions under study are given in Mikulin (1968), and Robinson and Stokes (1959). Reference » vs. m data sets of Miku- lin (1968), and Robinson and Stokes (1959) are in a good agreement. Data of Mikulin (1968) for La(NO,),-H,O are also in good agreement with the data of Rard (1987). However, the data of Mikulin (1968) cover the whole molality range of unsaturated and saturated solutions. In case of Ca(NO,),-H,O, UO,(NO,),-H,O, and Th(NO,),- HO systems, Mikulin also reported data for supersaturated solutions. In this study we parameterize the models using 1) all data of Mikulin (1968) for whole molalty range of unsaturated solutions from 0.1 m to m(max), 2) the data points at saturation (¢(sat)) (from Mikulin 1968), and 3) data for supersaturated Ca(NO,),-H,O, UO,(NO,).- HO, and Th(NO,),-H,O solutions (from Mikulin 1968), and 4) all experimental and recommended data of Rard and Spedding (1981) for Lu(NO,),-H,O system. In parameterization we used the value of Debye-Hiickel term (A®) equals to 0.39147 (Christov 2007, 2009, 2012). Following the parameterization scheme described in pre- vious paragraph the model for all 10 binary nitrate solutions is parameterized using two different approaches: (I) standard for N-1 electrolytes (N = 2, 3, or 4) approach with 3 ion interaction binary parameters (8, 6, and C’) and setting «, term equals Models for Nitrate Systems Oo 7. to 2, and «, = 0.0, and (II) an extended approach with four Pitzer ion interaction bi- nary parameters (8, 8", 8°,and C*) and varying in the values of «, and a, terms. As a first step in parameterization we used classical 3 parameters approach (I) and evaluate binary parameters using all available raw » data for whole molality range of solutions. As a next step, using the same ¢ data we re-parameterize the models on the basis of extended approach (II), and using three « — combinations: (Ia) «, = 2 and «, = 1, and (IIb) «, = 2 and «, = -1 (Christov 1996a, 1998, 1999, 2004, 2005) and (IIc) «, = 2 and a, = 0.3 (Guignot et al. 2019; Donchev and Christov 2020; Donchev et al. 2021). It was found that more combinations in “alfa” values do not improve the fit of data used in parameterization. The main criterion in the choice of established parameterization was the value of standard deviation (co) of fit of used ¢ data, i.e. parameterization with the lowest sigma value is accepted. For definition of sigma (co) see Christov and Moller (2004b), and Christov (2007, 2009, 2012). It was found that for 2 of studied systems Ba(NO,),-H,O, and Sr(NO,),-H,O, the approach (I) with 3 parameters (8, 6, C) give an acceptable agreement with the data. For these systems introducing into a model of fourth (8) parameter do not improve considerably the fit of data. For all other ni- trate systems under study we construct a model on the basis of extended approach (I), and using different combinations of “alfa’values: for Lu(NO,),-H,O system. 1) «, = 2 and a, = 0.3 (approach IIc), and 2) «, = 2 and a, = 1 (approach IIa), and 3) «, = 2 and a, = -1 (approach Ib). The resulting models fits the data up to supersaturation zone (m(max) = 14.77 m in Ca(NO,),-H,O) with sigma values, which is much less than the sigma values of models of Pitzer and Mayorga (1973), and of Kim and Frederick (1988). On next Figure 1 we present a comparison of osmotic coefhcients in nitrate bi- nary solutions 2-1 (Mg(NO,),-H,O, Ca(NO,),-H,O, Ba(NO,),-H,O, Sr(NO,).- H,O, and UO,(NO,),-H,O), 3-1 (Cr(NO,),-H,O, Al(NO,),-H,O, La(NO,),-H,O, Lu(NO,),-H,O), and 4-1 (Th(NO,) -H,O) calculated by the accepted models devel- oped here (heavy solid lines), and models developed by other authors (dashed lines and light solid lines: Pitzer and Mayorga (1973), Kim and Frederick (1988), Rard and Spedding (1981) (for Lu(NO,),-H,O system only), Rard et al. (2004) (for Mg(NO,).- H,O system only), and Wijesinghe and Rard (2005) (for Ca(NO,),-H,O system only). The recommended osmotic coefficients values given in literature at 25 °C are given on Fig. 1 by symbols. The vertical lines on the figures denote the molality of solutions saturated with corresponding nitrate solid phase (m(sat)), taken from Mikulin (1968). Excellent accepted new model (heavy solid line) — experiment (symbols) agreement has been obtained for all 10 systems and from low (see figures for Mg(NO,),-H,O, Ca(NO,),-H,O systems) and up to very high molality. As is shown Fig. 1, the new model for Ca(NO,),-H,O is in excellent agreement not only with the data at high mo- lality ((m(max) = 14.77 m), but contrary to the models of Kim and Frederick (1988) and Wijesinghe and Rard (2005) also in low molality range. It should be noted that all reference models presented on Fig. 1 by dashed, dashed-dotted, and light solid lines (Kim and Frederick (1988), Pitzer and Mayorga (1973), Rard et al. (2004), Wijesinghe and Rard (2005), and Rard and Spedding (1981)) have been constructed on the basis of standard Pitzer approach with 3 interaction parameters. Therefore, these models 398 Stanislav Donchev et al. / BioRisk 17: 389-406 (2022) ZO Mg(NOs HO 25°C:Mg(NOz)-H0 25°C Ca(NOz)-H20 Mg(NOz)2.6H20(cr 25 shapti eu YM=Rard et al. y, 15 Ca(NOs)2.4H2O(cr) a YM=Rard et al. (2004): a\=1.55; 4 # (2004): a,=1.55; ——,. light solid line; fe ’ light solid line: fp , wat ee He “ E This study: j ¢ & ad Ee This study: = solid line Z 2 eS solid line a 5 F Y 8 ; d gt wes 2 ye 2 an 2 a a] F, -s Kk S % P&M = fo Rard etal. 5 a 5 \+- - (2004): a,=2; & I Rard et al. fad “ess . \ ° “ dashed-dotted et (2004): a4=2: W&R(2005)=YM tN line; arr dashed-dotted line; xi ine: 05 ms , % ; 4 f 0.2 o4 O06 O8 - 0 2 4 6 8 m(Mg(NO,},)/mol.kg-" m(Mg(NOs),mol.kg” m(Ca(NO,),mol.kg" 1 ~ 1.1 25° Ca(NOs}zH20 2 25°C: Sr(NOs), +H,O 26°C Supersaturated Ca(NOz)}2-H20 Se Metast. Ca(NOs}.3H,O(cr) F Sr{NO-), (ct) 2 K&F(1988) E1, re new model, > : 3 This study = id li g @ solid line a See meee : , solid line > i ane ™. 3 P| = : 2 os 6 WAR(2005)=YM e | vg Eat ° 3 ‘ ra) Mikulin, empty omy 075 | Mit ants FE | triangles (exp) rs : Wa&R(2005)=¥M emal ae Stable Ca(NO3)z.4H20(cr) Pitzer+Mayorga: dashed-dotted line 07 1 1 1 4 : 0.5 L L in in i 4 n 0 01 02 03 04 05 7 g 11 13 15 005 115 2 25 3 35 4 m(Ca(NO,),)/mol.kg m(Ca(NO;);\/mol.kg"! $r(NO,),/mol.kg* - 7 2 | 25°C:Ba(NO3).+H,O 25°C:U0-(NO;);-H,0 8 (35°C, ONO, HO ol 18 e » Cr(NOy}s+Ho on : . 16 Mikulin: empty triangles; : Mikulin (1968) : = UO2(NOs)-6H20(cr) 13 Linesinew models 1,2, and » + empty triangles; = 14 212 oO a s 1 new model, e 42 2 1,1 8 solid line Hy new model, 8 = a Ba(NO,),(cr 3 solid line = E07 ae (NO3),(cr) é 08 20.9 i? O06 “a8 0 1 2 3 4 5 6 07 : : : key! tt) 05 1 1,5 04 m(UO,(NO,)\mol.kg Cr(NO,)'mol. ke" 0 02 OS O75 1 Ba(NO;)./mol.kg* Figure |. Comparison of model calculated (lines) osmotic coefficients (y) of Mg(NO,), Ca(NO,),, Ba(NO,),, Sr(NO,), UO,(NO,),, Cr(NO,) i AI(NO,) * La(NO,),, Lu(NO,),, and Th(NO,), in binary so- lutions 2-1 (Mg(NO,).-H,O, Ca(NO,),-H,O, Ba(NO,),-H,O, Sr(NO,).-H,O, and UO,(NO,),-H,O), 3-1 (Cr(NO,),-H,O, AI(NO,),-H,O, La(NO,),-H,O, Lu(NO,),-H,O), and 4-1 (Th(NO,),-H,O) against molality at T = 298.15 K, with recommendations in literature (symbols). For Mg(NO,),-H,O and Ca(NO,),-H,O systems an enlargement of the low molality corner is also given. Heavy solid lines represent the predictions of the developed in this study and accepted models. Dashed-dotted, dashed and light solid lines represent the predictions of the reference models of Kim and Frederick (1988), of Pitzer and Mayorga (1973), of Rard et al. 2004 (for Mg(NO,),.-H,O), of Wijesinghe and Rard (2005) (for Ca(NO,),-H,O), and of Rard and Spedding (1981) (for Lu(NO,),-H,O). On Figures YM denotes YMTDB (Sandia National Laboratories 2007). For Lu(NO,),-H,O the experimental data and recom- mended data are taken from Rard et al. (1977) (open squares), and Rard and Spedding (1981) (crosses), respectively. For all other systems the experimental data of Mikulin (1968) are used (open squares and open triangles). The molality of stable and metastable (for Ca(NO,),-H,O) crystallization of solid nitrate phases (m(sat)) is given on all figures by vertical lines (see Table 1). Models for Nitrate Systems 399 2,5 25°C; Al(NO3)s+H20 14 : a 20°C; La(NO3),+H30; not incl in YM 2 he sat point = 2.94 m Al(NO;);.9H,0 1.2 | This study: 4-parameters: B = a@;and a variation), solid line cI #2 15 a = = 1 ; : i | u ‘as we 1 open triangles: data of Mikulin 2 wo” 6 New Eva'ns: 2 (1968); 2 me standard approach | é solid line: model |; 508 “aa [ with 3-param's) 0.5 dashed line: mode! II _-°7" K&F(1988)- dashed line 0.6 4 i n i n 0 T T T 1 0 1 2 3 e a 2 7 4 m(La(NO3))imol.kg"? m(AI{NO,);)/mol.kg 2,5 , g |25°C: LU(NO2).+H,0 1.7 P5°C; Th(NOs),+H,O ; 1,5 + 217 r Th(NOz)46HsO(cr 12 Lu(NO3)..5H20 1,3 + ; = ' Ce. 27 311+ ett g = Pa 1,5 @ 0.9 € 9 P&M=YM § 13 2 OF a S ; 4—__ K&F git 505 | / E os Geo va 5°. 034 # Heavy solid line-New Eva'ns 07 4 \ Lf Ly 0.4 ni n 1 i 1 L 1 1 fl n 1 ee , 0,5 0 1 2 3 4 5 6 o 2 4 6 8 m(Th(NO,),)/mol.kg” m (Lu(NOs)3)/mol.kg* Figure |. Continue. cannot reproduce well the experimental data (Fig. 1). To illustrate this conclusion on Fig. 1 we give the predictions of two new models for La(NO,),-H,O system. As it is shown the 3 parameters model (light solid line) is in pure agreement with the data. The models for all nitrate binary systems under study are also validated by com- parison with recommendations given in literature (Rard and Spedding (1981) for Lu(NO,),-H,O); and Mikulin (1968)) (for all other 9 systems under study) on the mean activity coefficients (y,). These recommendations on y, are model-dependent. Therefore, they are not used in parameterization process, and only to validate the re- sulting models. The comparisons between predictions of new developed models and reference recommendations, which are not given here, show an excellent agreement from low to very high concentrations. Deliquescence relative humidity (DRH) calculations Deliquescence of single inorganic salt or their mixture is a process of spontaneous solid-liquid phase change. It is a process in which a soluble solid substance sorbs water vapor from the air to form a thermodynamically stable saturated aqueous solution on the surface of the particle. It is occurring when relative humidity (RH) in the gas-phase environment is at, or above deliquescence relative humidity (DRH) of the salt, or mutual deliquescence relative humidity (MDRH) of a salt mixture. Within the solid- 400 Stanislav Donchev et al. / BioRisk 17: 389-406 (2022) liquid equilibrium model, relative humidity is related to water activity(z,) (Clegg et al. 1998; Christov 2009, 2012; Donchev and Christov 2020) according to Eqn. (5): a,=P./P° = RH/100, (5) where P, and P° are the vapor pressure of the saturation solution and pure water, re- spectively, at given temperature. As a result, both DRH and MDRH of saturated surface solutions depend of temperature, the salt stoichiometry, and the solution composition. This process is of interest in many areas, such as heterogeneous chemistry of inorganic salts, corrosion of metals in wet atmosphere, in studies of chemistry of sea-type aero- sol atmospheric system (Kolev et al. 2013), and especially in development of strategies and programs for nuclear waste geochemical storage. Because of very high complicity of experiments, the relative humidity DRH experimental data are sparse. ‘Therefore, different sophisticated thermodynamic models have been proposed and developed to describe the deliquescence behavior of inorganic salts at wet conditions. In our previ- ous studies it was showed that calculations based on not high concentration restricted Pitzer models can be used for accurate determinations of both DRH and MDRH of saturated solutions in a wide range of temperatures, and compositions (Christov 2009, 2012; Donchev and Christov 2020). On the basis of evaluated binary parameters (8, 8, B,and C*) in this study we also determine water activity (4,) and Deliquescence Relative Humidity (DRH (%)) (eqn. 5) of 12 solid phases crystallizing from saturated binary nitrate solutions [Mg(NO,),.6H,O(s), Ca(NO,),.4H,O(s), Ca(NO,),.3H,O\(s), Ba(NO,),(s), Sr(NO,),(s), UO,(NO,),.6H,O(s), AI(NO,),.9H,O(s), Cr(NO,),(s), La(NO,),.6H,O(s), La(NO,),(s) Lu(NO,),.5H,O(s) and Th(NO,),.6H,O(s)]. Note that the widely used databases of Pitzer (1991), Pitzer and Mayorga (1973), Pitzer and Kim (1974) and Kim and Frederick (1988) do not consider solid phases. The results of calcu- lations are given in Table 1. The model DRH predictions are in excellent agreement with the experimental data determined using isopiestic method, and given in Mikulin (1968). According to model calculations the solid-liquid phase change of Ca(NO,),.3H,O(s), and Lu(NO,),.5H,O(s) occurs at lowest relative humidity of environment. It can be con- cluded that the solid-liquid phase change of solid nitrates of Lanthanide metals is more activated in the presence of calcium in the nuclear storage environment. Determination of thermodynamic solubility product (K°.,) of precipitates In this study we determine the thermodynamic solubility products (as K° ) of solid phases, precipitating from saturated nitrate binary solutions, s.a. anhydrous Ba(NO,),(s) and hydrate Ca(NO,),.3H,O(s), precipitating in Ba(NO,),-H,O and Ca(NO,),-H,O. The K°, have been determined on the basis of evaluated binary parameters and using experimental m(sat) solubility data, and using the following relationships (Christov 2005, 2007, 2009, 2012): K°sp (Ba(NO,),) = 4. Y (sat)? . m(sat)? K°sp(Ca(NO,),.3H,O) = 4. Y (sat) >, m(sat) °. a (sat) * (G) Models for Nitrate Systems 401 Table |. Comparison between model calculated and recommended values of the Deliquescence Relative Humidity [DRH (%) = a, (sat). 100; where a, (sat) is activity of water at saturation] and of the logarithm of the thermodynamic solubility product (as InK”_) of nitrate solid phases crystallizing from saturated binary solutions at T = 25° C. Salt composition m(sat) (exp) Ink iff DRH(%) (mol.kg’) This work Reference data This work Reference calculated calculated data* Mg(NO,),.6H,O(cr) 5.06 7.0098 7.02% 52.32 52.90 Ca(NO,) , 4H,O(cr) (stable solid) 8.41? 4.4362 4.53 49.07 49.10 Ca(NO,) , 3H,O(cr) (metastable solid) 14.77? 6.6449 5.34° (m/(sat) = 15.0 m) 22.52 - Ba(NO,),(cr) 0.39? -5.125 - 98.61 98.60 St(NO,), .4H,O (cr) 3.76% 0.0327 - 84.83 84.80 UO,(NO,), .6H,O(cr) 3.21 5.3022 5.251° 73.44 73.60 AI(NO,),.9H,O (cr) 3.16 4.3081 : 59.88 60.20 Cr(NO,), (cr) ¢ 1.4° 1.2097 - 86.38 - La(NO,), .6H,O (cr)¢ 4.6154 2.1599 297° 62.38 - La(NO,), (cr) * 2.94# 1.4704 < 77.73 77.60 Lu(NO,),5H,O (cr) 6.8154 10.7681 10.67 31.27 - Th(NO,),.6H,O(cr) 4.00 4.4886 4,71° (as Th(NO,),.5H,O) 54.46 55.0 Guignot et al. (2019); ‘Accepted m(sat) molality and stoichiometry of solid phase. As a next step, using the accepted new developed parameterizations, and experi- mentally determined molalities (m(sat) of the saturated binary solutions (Mikulin 1968; Guignot et al. 2019; Lassin et al. 2020)) we calculate the logarithm of the thermodynamic solubility product (In K°) of twelve nitrate solid phases crystal- lizing from saturated binary nitrate solutions at 25 °C (Eqn. (6)). The model cal- culations are given in Table 1. With only 2 exceptions (for Ca(NO,),.3H,O(s) and La(NO,),.6H,O(s)) a good agreement has been obtained with calculations of Lach et al. (2018), Lassin et al. (2020), and Guignot et al. (2019) for all nitrate solids. The In K°. | differences are mainly due on the 1) different m(sat) values used in calcula- tions (see Egn. (6)), and 2) different experimental data source with different m(max) values used in parameterization. Summary and conclusions In this study we developed new thermodynamic models for solution behavior and solid-liquid equilibrium in 10 nitrate binary systems of the type 2-1 (Mg(NO,).- H,O, Ca(NO,),-H,O, Ba(NO,),-H,O, Sr(NO,),-H,O, and UO,(NO,),-H,O), 3-1 (Cr(NO,),-H,O, Al(NO,),-H,O, La(NO,),-H,O, Lu(NO,),-H,O), and 4-1 (Th(NO,),-H,O) from low to very high concentration at 25 °C. To parameterize models for binary systems we used all available raw experimental osmotic coefficients data () for whole concentration range of solutions, and up to saturation point. Data for super- saturation zone, available for Ca(NO,),-H,O, UO,(NO,),-H,O, and Th(NO,),-H,O systems, are also included in parameterization. To construct models, we used different versions of standard molality-based Pitzer approach. It was established that with only 2 exceptions (Ba(NO,),-H,O, and UO,(NO,),-H,O) application of extended approach 402 Stanislav Donchev et al. / BioRisk 17: 389-406 (2022) with 4 parameters (8, 6", 8°,and C*) and variation of «, term in fundamental Pitzer equations leads to the lowest values of standard model-experiment deviation. The pre- dictions of new developed here models are in excellent agreement with experimental os- motic coeflicients data (see Fig. 1), and with recommendations on activity coefficients (not given here) in binary solutions from low to very high concentration: up to 14.77 mol. kg" in Ca(NO,),-H,O. The Deliquescence Relative Humidity (DRH), and thermody- namic salubiliey reduce (as In K°, ») of of 12 solid phases crystallizing from saturated bi- nary nitrate solutions [Mg(NO, ), Hae sO(s), Ca(NO,),.4H,O(s), Ca(NO,),.3H,O\(s), Ba(NO,),(s), Sr(NO,),(s), UO (NO, ),-6H,O(s), AL(NO, Jeet COKS) Cx(NO, ),(s), La(NO,),.6H,O(s), La(NO,),(s) Lu(NO,),.5H,0() and “Th(NO,),.6H,O(s)] fae been determined on the basis of evaluated binary parameters and using experimental m(sat) solubility data. Model predictions are in good agreement with available refer- ence data. The accurate solid-liquid equilibrium models for nitrate systems described in this study are of high importance for development of strategies and programs for nuclear waste geochemical storage. Acknowledgement We wish to thank the reviewers (Dr. Krasimir Kostov and anonymous reviewer) for their constructive suggestions and helpful comments. The manuscript was improved considerably through their comments. The work was supported by the European Re- gional Development Fund, Project BGO5M2OP001-1.001-0004, and by Shumen University Research Program, Project RD-08-131/04.02.2021 References Altmaier M, Brendler V, Bube C, Neck V, Marquardt C, Moog HC, Richter A, Scharge T, Voigt W, Wilhelm S, Wilms T, Wollmann G (2011) THEREDA-Thermodynamische Ref- erenzdatenbasis. Report GRS 265. Christov C (1995) Thermodynamic study of (b,LiBr + b,MgBr,)(aq), where b denotes mo- lality, at the temperature 348.15 K. The Journal of Chemical Thermodynamics 27(11): 1267-1273. https://doi.org/10.1006/jcht.1995.0133 Christov C (1996a) Thermodynamic study of the co-crystallization of 2RbCI.NiCl,.2H,O and 2RbCI.MnCL,.2H,O at the temperature 298.15 K. The Journal of Chemical Thermody- namics 28(7): 743-752. https://doi.org/10.1006/jcht.1996.0068 Christov C (1996b) Thermodynamics of the aqueous sodium and magnesium bromide sys- tem at the temperatures 273.15 K and 298.15 K. Calphad 20(4): 501-509. https://doi. org/10.1016/S0364-5916(97)00012-6 Christov C (1998) Thermodynamic study of the KCI-K,SO,-K,Cr,O,-H,O system at the temperature 298.15 K. Calphad 22(4): 449-457. https://doi.org/10.1016/S0364- 5916(99)00004-8 Models for Nitrate Systems 403 Christov C (1999) Study of (m,KCl + m,MeCL)(aq), and (m,K,SO, + m,MeSO,)(aq) where m denotes molality and Me denotes Cu or Ni, at the temperature 298.15 K. The Journal of Chemical Thermodynamics 31(1): 71-83. https://doi.org/10.1006/jcht.1998.0419 Christov C (2001) Thermodynamic study of the NaCl-Na,SO,-Na,Cr,O,-H,O system at the temperature 298.15 K. Calphad 25: 11-17.https://doi.org/10.1016/S0364- 5916(01)00025-6 Christov C (2002a) Thermodynamics of formation of ammonium, sodium, and potassium alums and chromium alums. Calphad 26(1): 85-94. https://doi.org/10.1016/S0364- 5916(02)00026-3 Christov C (2002b) Thermodynamic study of quaternary systems with participation of am- monium and sodium alums and chromium alums. Calphad 26(3): 341-352. https://doi. org/10.1016/S0364-5916(02)00049-4 Christov C (2003a) Thermodynamics of formation of double salts M,SO,.MeSO,.6H,O and M,SeO,,.MeSeO ,.6H,O where M denotes Rb, or Cs, and Me denotes Co, Ni or Zn. The Journal of Chemical Thermodynamics 35(11): 1775-1792. https://doi.org/10.1016/j. jct.2003.08.004 Christov C (2003b) Thermodynamic study of the co-crystallization of ammonium, sodium and potassium alums and chromium alums. Calphad 27(2): 153-160. https://doi.org/10.1016/ S0364-5916(03)00046-4 Christov C (2003c) Thermodynamic study of aqueous sodium, potassium and chromium chloride systems at the temperature 298.15 K. The Journal of Chemical Thermodynamics 35(6): 909-917. https://doi.org/10.1016/S0021-9614(03)00042-9 Christov C (2004) Pitzer ion-interaction parameters for Fe(II) and Fe(II) in the quinary {Na+K+Mg+Cl+SO,+H,O} system at T = 298.15 K. The Journal of Chemical Thermody- namics 36(3): 223-235. https://doi.org/10.1016/j.jct.2003.11.010 Christov C (2005) Thermodynamics of formation of double salts and solid solutions from aqueous solutions. The Journal of Chemical Thermodynamics 37(10): 1036-1060. htt- ps://doi.org/10.1016/j.jct.2005.01.008 Christov C (2007) An isopiestic study of aqueous NaBr and KBr at 50°C. Chemical Equilib- rium model of solution behavior and solubility in the NaBr-H,O, KBr-H,O and Na-K-Br- H,O systems to high concentration and temperature. Geochimica et Cosmochimica Acta 71(14): 3357-3369. https://doi.org/10.1016/j.gca.2007.05.007 Christov C (2009) Chemical equilibrium model of solution behavior and solubility in the Mg- CL,-H,O, and HCI-MgCL,-H,O systems to high concentration from 0°C to 100°C. Jour- nal of Chemical & Engineering Data 54: 2599-2608. https://doi.org/10.1021/je900135w Christov C (2012) Study of bromide salts solubility in the (m, KBr + m,CaBr,)(aq) system at T = 323.15 K. Thermodynamic model of solution behavior and solid-liquid equilibria in the ternary (m,KBr + m,CaBr,)(aq), and (m,MgBr, + m,CaBr,)(aq), and in quinary {Na+K+Mg+Ca+Br+H,O} systems to high concentration and temperature. The Journal of Chemical Thermodynamics 55: 7-22. https://doi.org/10.1016/j.jct.2012.06.006 Christov C (2020) Thermodynamic models for solid-liquid equilibrium of aluminum, and aluminum-silicate minerals in natural fluids. Current state and perspectives. Review of the Bulgarian Geological Society 81: 69-71. 404 Stanislav Donchev et al. / BioRisk 17: 389-406 (2022) Christov C, Moller N (2004a) Chemical equilibrium model of solution behavior and solu- bility in the H-Na-K-Cl-OH-HSO,-SO,-H,O system to high concentration and tem- perature. Geochimica et Cosmochimica Acta 68: 1309-1331. https://doi.org/10.1016/j. gca.2003.08.017 Christov C, Moller N (2004b) A chemical equilibrium model of solution behavior and solubil- ity in the H-Na-K-Ca-Cl-OH-HSO,-SO,-H,O system to high concentration and temper- ature. Geochimica et Cosmochimica Acta 68(18): 3717-3739. https://doi.org/10.1016/j. gca.2004.03.006 Christov C, Ojkova T, Mihov D (1998) Thermodynamic study of (m,Na,SeO, + m,NiSeO,) (aq), where m denotes molality, at the temperature 298.15 K. The Journal of Chemical Thermodynamics 30(1): 73-79. https://doi.org/10.1006/jcht.1997.0274 Christov C, Zhang M, Talman S, Reardon E, Yang T (2012) Review of issues associated with evaluation of Pitzer interaction parameters. Mineralogical Magazine 76: 1578. https:// goldschmidtabstracts.info/abstracts/abstract View?id=2012001520 Clegg S, Brimblecombe P, Wexler A (1998) Thermodynamic model of the system H*-NH,'- Na*- SO,7-NO,-CI-H,O at 298.15 K. The Journal of Physical Chemistry A 102(12): 2155-2171. https://doi.org/10.1021/jp973043; Donchev S, Christov C (2020) Development of Accurate Chemical Thermodynamic Database for Geochemical Storage of Nuclear Waste. Part I: Models for Predicting Solution Proper- ties and Solid-Liquid Equilibrium in Binary Nitrate Systems of the Type 1-1. Ecologia Balkanica, Special Edition 3: 195-210. http://eb.bio.uni-plovdiv.bg Donchev S, Tsenov T, Christov C (2021) Chemical and geochemical modeling. Thermody- namic models for binary fluoride systems from low to very high concentration (> 35 m) at 298.15 K. Acta Scientifica Naturalis 8(2): 1-15. https://doi.org/10.2478/asn-2021-0014 Dossier ANDRA (2005) Argile. Evaluation de La Faisabilité Du Stockage Géologique En For- mation Argileuse. https://www.andra.fr/sites/default/files/2017-12/266.pdf El Guendouzi M, Marouani M (2003) Water activities and osmotic and activity coefficients of aqueous solutions of nitrates at 25 °C by the hygrometric method. Journal of Solution Chemistry 32(6): 535-546. https://doi.org/10.1023/A:1025365900350 Gruszkiewicz M, Simonson J (2005) Vapor pressures and isopiestic molalities of concentrated CaCl,(aq), and CaBr,(aq), and NaCl(aq) to T = 523.15 K. The Journal of Chemical Ther- modynamics 37(9): 906-930. https://doi.org/10.1016/j.jct.2004.12.009 Guignot S, Lassin A, Christov C, Lach A, André L, Henocg P (2019) Modelling the osmotic and activity coefficients of lanthanide nitrate aqueous solutions at 298.15 K from low molalities to (super) saturation. Journal of Chemical & Engineering Data 64(1): 345-359. https://doi.org/10.1021/acs.jced.8b00859 Harvie C, Moller N, Weare J (1984) The prediction of mineral solubilities in natural wa- ters: The Na-K-Mg-Ca-H-Cl-SO,-OH-HCO,-CO,-CO2-H,0 system from zero to high concentration at 25°C. Geochimica et Cosmochimica Acta 48(4): 723-751. https://doi. org/10.1016/0016-7037(84)90098-X Kim H-T, Frederick Jr W (1988) Evaluation of Pitzer ion interaction parameters of aqueous electrolytes at 25°C. 1. Single salt parameters. Journal of Chemical & Engineering Data 33(2): 177-184. https://doi.org/10.1021/je00052a035 Models for Nitrate Systems 405 Kolev H, Tyuliev G, Christov C, Kostov K (2013) Experimental study of the surface chemical composition of sea salt crystallized at evaporation of sea water under natural conditions. Bulgarian Chemical Communications 45: 584-591. http://www.bcc.bas.bg/BCC_ Lach A, André L, Guignot S, Christov C, Henocq P, Lassin A (2018) A Pitzer parameteriza- tion to predict solution properties and salt solubility in the H-Na-K-Ca-Mg-NO,-H,O system at 298.15 K. Journal of Chemical & Engineering Data 63(3): 787-800. https:// doi.org/10.1021/acs.jced.7b00953 Lassin A, Christov C, André L, Azaroual M (2015) A thermodynamic model of aqueous elec- trolyte solution behavior and solid-liquid equilibrium in the Li-H-Na-K-Cl-OH-H,O system to very high concentrations (40 Molal) and from 0 to 250°C. American Journal of Science 315(3): 204-256. https://doi.org/10.2475/03.2015.02 Lassin A, Guignot S, Lach A, Christov C, André L, Madé B (2020) Modeling the solution prop- erties and mineral-solution equilibria in radionuclide-bearing aqueous nitrate systems. Ap- plication to binary and ternary systems containing U, ‘Th or lanthanides, at 25°C. Journal of Chemical & Engineering Data 65(7): 3613-3626. https://doi.org/10.1021/acs.jced.0c00180 Maliutin AS, Kovalenko NA, Uspenskaya IA (2020) Thermodynamic Properties and Phase Equilibria in the HO-HNO,—UO,(NO,), System. Moscow University Chemistry Bul- letin 75(2): 65-71. https://doi.org/10.3103/S002713142002009 1 Mikulin G (1968) Voprosy Fizicheskoi Khimii Electrolytov. Izd. Khimiya, St. Petersburg, 417 pp. Moller N, Christov C, Weare J (2006) Thermodynamic models of aluminumsilicate mineral solu- bility for application to enhanced geothermal systems. Proceedings 31% Workshop on Geo- thermal Reservoir Engineering. Stanford University, Stanford, California, January 30—Febru- ary 1, 8 pp. https://pangea.stanford.edu/ERE/pdf/ IGAstandard/SGW/2006/moller.pdf Moller N, Christov C, Weare J (2007) Thermodynamic model for predicting interactions of geothermal brines with hydrothermal aluminum silicate minerals. Proceedings 32" Work- shop on Geothermal Reservoir Engineering. Stanford University, Stanford, California, Jan- uary 22-24, 8 pp. https://geo.stanford.edu/ ERE/pdf/IGAstandard/SGW/2007/moller.pdf Ojkova T, Christov C, Mihov D (1999) Thermodynamic study of (NH,),SeO, (aq) and K,SeO, (aq) at the temperature 298.15 K. Monatshefte fir Chemie 130: 1061-1065. https://doi. org/10.1007/PL00010283 Park JH, Christov C, Ivanov A, Molina M (2009) On OH uptake by sea salt un- der humid Conditions. Geophysical Research Letters 36(2): LO2802. https://doi. org/10.1029/2008GL036160 Petrenko S, Pitzer K (1997) Thermodynamics of aqueous NaOH over the complete composi- tion range and to 523 K and 400 Mpa. The Journal of Physical Chemistry B 101(18): 3589-3595. https://doi.org/10.1021/jp963707 Pitzer KS (1973) Thermodynamics of Electrolytes. I. Theoretical Basis and General Equations. Journal of Physical Chemistry 77(2): 268-277. https://doi.org/10.1021/j100621a026 Pitzer KS (1991) Ion Interaction Approach: Theory and Data Correlation. Chapter 3 of Activ- ity Coefficients in Electrolyte Solutions. 2"4 edn. CRC Press, Boca Raton, Florida. Pitzer KS, Kim JJ (1974) Thermodynamics of Electrolytes. IV. Activity and Osmotic Coef- ficients for Mixed Electrolytes. Journal of the American Chemical Society 96(18): 5701- 5707. https://doi.org/10.1021/ja00825a004 406 Stanislav Donchev et al. / BioRisk 17: 389-406 (2022) Pitzer KS, Mayorga G (1973) Thermodynamics of electrolytes. II. Activity and osmotic coef- ficients for strong electrolytes with one or both ions univalent. Journal of Physical Chem- istry 77(19): 2300-2308. https://doi.org/10.1021/j100638a009 Pitzer KS, Mayorga G (1974) Thermodynamics of electrolytes. III. Activity and osmotic coef- ficients for 2—2 electrolytes. Journal of Solution Chemistry 3(7): 539-546. https://doi. org/10.1007/BF00648138 Rard JA (1987) Osmotic and activity coefficients of aqueous La(NO3)3 and densities and ap- parent molal volumes of aqueous Eu(NO3)3 at 25 °C. Journal of Chemical & Engineering Data 32: 92-98. https://doi.org/10.1021/je00047a026 Rard JA, Spedding FH (1981) Isopiestic determination of the activity coefficients of some aque- ous rare-earth electrolyte solut ions at 25 °C. 5. Dy(NO3)3Ho(NO3)3 and Lu(NO3)3. Jour- nal of Chemical & Engineering Data 26: 391-395. https://doi.org/10.1021/je00026a013 Rard JA, Shiers LE, Heiser DJ, Spedding FH (1977) Isopiestic determination of the activity coefh- cients of some aqueous rare earth electrolyte solutions at 25 °C. 3. The rare earth nitrates. Jour- nal of Chemical & Engineering Data 22(3): 337-347. https://doi.org/10.1021/je60074a015 Rard JA, Wijesinghe AM, Wolery TJ (2004) Review of the Thermodynamic Properties of Mg(NO,),(aq) and Their Representation with the Standard and Extended Ion-Interaction (Pitzer) Models at 298.15 K. Journal of Chemical & Engineering Data 49(5): 1127-1140. https://doi.org/10.1021/je0498681 Robinson R, Stokes R (1959) Electrolyte Solutions, 2"4 edn.; Butterworths, London. Sandia National Laboratories (2007) Qualification of thermodynamic data for geochemical modeling of mineral-water interactions in dilute systems (data0.ypf.R2) U.S. Department of Energy, ANL-WIS-GS-000003 REV 01. Trendafelov D, Prangova D, Nishev M, Christov C (1995a) Study of the conversion of BaSO, into BaCO, in the fourcomponent water-salt system BaSO,,+Na,CO, = BaCO,+Na,SO,,. Compt. rend. Acad. Bulg. Sci. 48: 39-41. Trendafelov D, Christov C, Balarew C, Karapetkova A (1995b) Study of the Conversion of CaSO, to CaCO, within the CaSO,+Na,CO, = CaCO,+Na,SO, fourcomponent water- salt system. Collection of Czechoslovak Chemical Contributions 60(12): 2107-2111. https://doi.org/10.1135/cccc19952107 Tsenov IT, Donchev S, Christov C (2022) Development of accurate chemical thermodynamic database for geochemical storage of nuclear waste. Part III: Models for predicting solution properties and solid-liquid equilibrium in cesium binary and mixed systems. In: Chankova S, Peneva V, Metcheva R, Beltcheva M, Vassilev K, Radeva G, Danova K (Eds) Current trends of ecology. BioRisk 17: 407-722. https://doi.org/10.3897/biorisk.17.77523 Wang P, Anderko A, Springer RD, Young RD (2006) Modeling phase equilibria and specia- tion in mixed-solvent electrolyte systems: II. Liquid-liquid equilibria and properties of associating electrolyte solutions. Journal of Molecular Liquids 125(1): 37-44. https://doi. org/10.1016/;.molliq.2005.11.030 Wijesinghe AM, Rard JA (2005) Conversion and Optimization of the Parameters from an Ex- tended Form of the Ion-Interaction Model for Ca(NO3)2(aq) and NaNO3(aq) to Those of the Standard Pitzer Model, and an Assessment of the Accuracy of the Parameter Tem- perature Representations. The Journal of Chemical Thermodynamics 37(11): 1196-1218. https://doi.org/10.1016/j.jct.2005.02.013