CHINESE JOURNAL OF ENERGETIC MATERIALS
+Advanced Search
en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。
参考文献 1
CAOXin‑mao, LIFu‑ping . HMX and its application[M]. Beijing:Weapon Industry Press,1993:104-106.
参考文献 2
CAODuan‑lin, WANGJian‑long . A study on the recovery of small quantity of RDX in the producing of HMX with acetic anhydride[J]. Journal of North China Institute of Technology, 1996, 17(3): 279-282.
参考文献 3
YELing, WUYu‑zhong, LIZhi‑hua . Study on transformation of α‑HMX into β‑HMX with spent acetic acid[J]. Chinese Journal of Explosives & Propellants, 2000, 23(1): 38-39, 42.
参考文献 4
LIQiao‑ling, YEYu‑peng . A novel process for the refining of HMX[J]. Acta Armamentarii, 2002, 23(4): 555-557.
参考文献 5
ZHAORui‑xian, LIGuang‑ming, GAOTian‑ping, et al . Study on RDX granulate stage process[J]. Chinese Journal of Explosives & Propellants, 2003, 26(4): 67-70.
参考文献 6
ZHANGMin‑hua, DOUMao‑bin, WANGMeng‑yan,et al . Study on the solubility parameter of supercritical carbon dioxide system by molecular dynamics simulation[J]. Journal of Molecular Liquids, 2017, 248: 322-329.
参考文献 7
LUOYan‑long, WANGRun‑guo, WANGWei, et al . Molecular dynamics simulation insight into two‑component solubility parameters of graphene and thermodynamic compatibility of graphene and styrene butadiene rubber[J]. Journal of Physical Chemistry C,2017, 121(18): 10163-10173.
参考文献 8
GoharshadiE K, AkhlamadiG, MahdizaehS J . Investigation of graphene oxide nanosheets dispersion in water based on solubility parameters: a molecular dynamics simulation study[J]. RSC Advance,2015, 5(129): 106421-106430.
参考文献 9
LeachA R . Molecular modelling: principles and applications[M]. Beijing: World Book Inc, 2003.
参考文献 10
ChoiC S, BennemaP . The crystal structure of cyclotrimethylenetrinitramine[J]. Acta Crystallogr B, 1972, 28: 2857-2862.
参考文献 11
BartonA F M . Solubility parameters[J]. Chemical Reviews, 1975, 75(6): 731-753.
参考文献 12
SUNYe‑bin, HUIJun‑ming, CAOXin‑mao . Military Use Blended Explosives[M]. Beijing: Weapon Industry Press, 1995: 242-251.
参考文献 13
LIJing, JINShao‑hua, LANGuan‑chao, et al . Morphology control of 3‑nitro‑1,2,4‑triazole‑5‑one (NTO) by molecular dynamics simulation[J]. Crystengcomm, 20(40): 6252-6260.
参考文献 14
LANGuan‑chao, JINShao‑hua, LIJing, et al .The study of external growth environments on the crystal morphology of epsilon‑HNIW by molecular dynamics simulation[J]. Journal of Materials Science,2018, 53(18): 12921-12936.
参考文献 15
LIJing, JINShao‑hua, LANGuan‑chao, et al . Molecular dynamics simulations on miscibility, glass transition temperature and mechanical properties of PMMA/DBP binary system[J]. Journal of Molecular Graphics & Modelling,2018, 84: 182-188.
参考文献 16
LANGuan‑chao, JINShao‑hua, LIJing, et al . Miscibility, glass transition temperature and mechanical properties of NC/DBP binary systems by molecular dynamics[J]. Propellants, Explosives, Pyrotechnics,2018, 43(6): 559-567.
参考文献 17
AlmarriF, HaqN, AlanaziF K, et al . Solubilityand thermodynamic function of vitamin D3 in different mono solvent[J].Journal of molecular liquids,2017, 229: 477-481.
参考文献 18
ShakeelF, ImranM, Abida, HaqN, et al . Solubility and thermodynamic/solvation behavior of 6‑phenyl‑4,5‑dihydropyridazin‑3(2H)‑one in different(Transcutol‑water) mixtures[J].Journal of molecular liquids,2017, 230: 511-517.
参考文献 19
CHENLi‑zhen, ZHANGTian‑bei, LIMan, et al . Solubility of 2,2',4,4',6,6'‑hexanitrostilbene in binary solvent of N, N‑dimethylformamide and acetonitrile[J]. Journal of Chemical Thermodynamics, 2016, 99: 99-104.
参考文献 20
LANGuan‑chao, WANGJian‑long, CHENLi‑zhen, et al . Measurement and correlation of the solubility of 3,4‑bis(3‑nitrofurazan‑4‑yl)furoxan (DNTF) in different solvents[J].Journal of Chemical Thermodynamics, 2015, 89: 264-269.
参考文献 21
ApelblatA, ManzurolaE .Solubilities of manganese, cadmium, mercury and lead acetates in water from T=278.15 K to T=340.15 K[J]. Journal of Chemical Thermodynamics,2001, 33(2): 147-153.
参考文献 22
ApelblatA, ManzurolaE . Solubilities of L‑aspartic, DL‑aspartic, DL‑glutamic, p‑hydroxybenzoic, o‑anisic, p‑anisic, and itaconic acids in water from T=278 K to T=345 K[J]. Journal of Chemical Thermodynamics,1997, 29(12): 1527-1533.
参考文献 23
ApelblatA, ManzurolaE . Solubilities of o‑acetylsalicylic, 4‑aminosalicylic, 3,5‑dinitrosalicylic, and p‑toluic acid, and magnesium‑DL‑aspartate in water from T=(278 to 348) K[J]. Journal of Chemical Thermodynamics,1999, 31(1): 85-91.
参考文献 24
ZHOULi, ZHANGPei‑pei, YANGGuang‑de, et al .Solubility of chrysin in ethanol and water mixtures[J]. Journal of Chemical and Engineering Data,2014, 59(7): 2215-2220.
参考文献 25
CHENZhao‑guo, YANGWen‑ge, HUYong‑hong, et al . Measurement and correlation for the solubility of dimethyl 1,4‑cyclohexanedione‑2,5‑dicarboxylate in different solvents at temperatures from (278.15 to 323.15) K[J]. Journal of Chemical and Engineering Data,2011, 56(5):2726-2729.
参考文献 26
AcreeJrW E .Comments concerning ‘model for solubility estimation in mixed solvent system’[J].International Journal of Pharmaceutics,1996, 127(1): 27-30.
参考文献 27
Jouyban‑GharamalekiA, HanaeeJ . A general model from theoretical cosolvency models[J].International Journal of Pharmaceutics,1997, 152(2): 247-250.
参考文献 28
LIUBao‑shu, SUNHua, WANGJing‑kang, et al . Solubility of disodium 5′‑guanylate heptahydrate in aqueous methanol mixtures[J].Food Chemistry,2011,128(1): 218-221.
参考文献 29
JouybanA . Review of the cosolvency models for predicting solubility of drugs in water‑cosolvent mixtures[J].Journal of Pharmacy and Pharmaceutical Sciences,2008, 11(1): 32-57.
参考文献 30
WANGShui,QINLi‑ying, ZHOUZhi‑mao, et al . Solubility and solution thermodynamics of betaine in different pure solvents and binary mixtures[J]. Journal of Chemical and Engineering Data,2012,57(8): 2128-2135.
参考文献 31
WANGGuan, WANGYong‑li, You‑guangMA, et al . Determination and correlation of cefuroxime acid solubility in (acetonitrile‑water) mixtures[J].Journal of Chemical Thermodynamics, 2014, 77: 144-150.
参考文献 32
ZHANGHui, YINQiu‑xiang, LIUZeng‑kun, et al . Measurement and correlation of solubility of dodecanedioic acid in different pure solvents from T=(288.15 to 323.15) K[J]. Journal of Chemical Thermodynamics,2014, 68: 270-274.
参考文献 33
DelgadoD R, HolguinA R, AlmanzaO, et al . Solubility and preferential solvation of meloxicam in ethanol‑water mixtures[J]. Fluid Phase Equilibria,2011, 305(1): 80-85.
参考文献 34
HolguinA R, DelgadoD R, MartinezF, et al . Solution thermodynamics and preferential solvation of meloxicam in propylene glycol‑water mixtures[J]. Journal of Solution Chemistry, 2011, 40(12):1987-1999.
目录 contents

    Abstract

    To research the dissolution characteristics of Hexogen (RDX) in ethyl acetate–water binary system, the solubility of RDX at temperatures ranging from 298.15 to 338.15 K was measured The experimental data were fitted using Apelblat equation, CNIBS/R‑K model and Jouyban‑Acree equation. The standard enthalpy of dissolution, entropy of dissolution and Gibbs free energy were calculated. To study the dissolution mechanism, the solubility parameters of RDX in ethyl acetate‑water binary mixed solvent were calculated by molecular dynamics simulation method. Results show that the solubility of RDX increases with the increase of temperature and water content. The fitting values of empirical equation are basically consistent with the experimental ones. The experimental solubility data, model parameters, thermodynamic properties and solubility parameters provide basic data and models for the recovery process of RDX and HMX.

    摘要

    为了研究黑索今(RDX)在乙酸乙酯‑水二元体系中的溶解特性,测定了RDX在298.15~338.15 K温度下的溶解度。分别采用Apelblat方程、CNIBS/R‑K模型及Jouyban‑Acree方程对实验数据进行拟合。算出了标准溶解焓、标准溶解熵及吉布斯自由能。为了研究溶解机理,采用分子动力学模拟方法计算了RDX在二元混合溶剂中的溶度参数。结果表明,RDX的溶解度随着温度和水含量的升高而增大。经验方程拟合值与实验值基本保持一致。实验所得的溶解度数据,模型参数、热力学性质、溶度参数都为RDX和HMX的回收过程提供了基础数据和模型。

  • 1 Introduction

    Octogen(HMX), with eight‑membered ring of nitramine structure, is a kind of preferable single explosive judged from high density, high detonation velocity and good thermal stability. It has wide applications in rocket propellant, ammunition and weapons. HMX can be prepared in many methods[1], among which acetic anhydride method with simple process, easy operation and production safety is popularized. In the production process of HMX by acetic anhydride method, a mass of “waste explosive” is produced, and the so‑called “waste explosive” is a mixture containing RDX, HMX and multiple by products, which precipitates from waste acid and waste solvents[2]. The interactions between RDX and HMX of “waste explosive” are mainly consist of van der Waals(vdW) force and electrostatic, which can be broken easily that makes it possible to separate the RDX and HMX. Since 1960s, the research of separation of RDX and HMX has entered a period of rapid development[2,3]. The representative methods are solvent method, crystallization with waste acid or nitric acid, mechanical separating method and complexation with cyclohexanone or dimethylformamide.

    Nowadays, crystallization with ethyl acetate (EtOAc) that has mature recovery technology and obvious difference in solubility of RDX and HMX is selected to separate the RDX and HMX. Then HMX of high quality is obtained by rotating crystal method[4] and RDX of multi‑granularity can be obtained through evaporative crystallization with acetone and water[5].

    The solubility is an extremely significant thermodynamic parameter for crystallization and calculation of other thermodynamic parameters. Because EtOAc is used in the separation of RDX and HMX, the measurement of the solubility of EtOAc is very important. In addition, the moisture content of raw material is pretreated to be below 25 percent, to leave out the drying process, during which we must consider the influence of moisture on solubility, so it is necessary to measure the solubility of RDX in ethyl acetate‑water mixed solvents with different mole ratios. The solubility results will support the design and operation of the HMX recycling process.

    In addition, molecular dynamics(MD) with provision of atomic‑scale information has become a powerful simulation technique. Some researchers have simulated the solubility parameter[6,7,8] by MD method to explore the effect of components on solubility.

    In the present work, we report the solubility of RDX in ethyl acetate‑water mixture measured with a gravimetric method and the solubility parameter calculated using MD simulation.This is quite useful in the study of recycling process of RDX and HMX.

  • 2 Materials and methodology

  • 2.1 Materials

    RDX (99.7%) was supplied by Research Institute of Gansu Yin Guang Chemical Industry Group Co.Ltd. Ethyl acetate (99.5%) without further purification was purchased from Sinopharm Chemical Reagent Co.Ltd. Deionized water was provided by laboratory‑made.

  • 2.2 Solubility measurements

    40 g RDX and 200 g EtOAc were added into a 500 mL flask with three necks under 300 r·min-1 with an agitator made from polytetrafluoroethylene at temperatures ranging from 298.15 K to 338.15 K. A thermostatic water bath(Shanghai Bilon Precision Instrument Co. Ltd., China) was used to control the temperature with an uncertainty of 0.05 K. A condenser was used to prevent evaporation of solvent during the experiments. After 120 min, the agitation was switched off and the solution was stood for 120 min. Then 5 mL of solution was sampled into a bottle weighed before hand, and the mass of solution could be calculated after the total mass was measured. Similarly, the mass of solute could be calculated after the solvent evaporated completely. Each mass was determined using an analytical balance (Mettler Toledo XS 105DU, Switzerland) with an uncertainty of 0.0001 g. The same solubility experiments were repeated three times to obtain the mean values. The operations above were repeated with the mole ratio of water ranging from 0.0000 to 0.1280.

  • 2.3 Computation

    MD simulation in present work was performed using Forcite module and Amorphous Cell in Materials Studio 6.0 package[9]. The single crystal structure of RDX was obtained from the Cambridge Structural Database(CCDC: 705291)[10]. All the geometry was optimized by Conjugate gradient algorithm in COMPASS force field. The comparison of the experimental lattice parameters and the optimized values of RDX are listed in Table 1. It can be illustrated that the relative error between the optimized lattice parameters and the experimental values is within 5%, demonstrating that COMPASS force field is suitable for the simulation of RDX.

    Table 1 Comparison of the experimental lattice parameters and the optimized values of RDX

    lattice

    parameter

    experimental

    value[10]

    optimized value

    relative

    error/%

    a13.18213.450-2.033
    b11.57411.2832.514
    c10.70910.2254.520
    β/(°)90.00090.0000.000

    The solubility parameter (δ) of a material is defined as square root of its cohesive energy density, and the molar cohesive energy is the energy associated with all the molecular interactions in a mole of the material[11]. According to the solubility parameter similarity rule, if the difference of solute and solvent is small, the solute will prefer to dissolve into the solvent[12].Therefore, the solubility parameter of RDX and ethyl acetate‑water is calculated to evaluate the miscibility. The supercell of RDX and AC models of pure EtOAc and EtOAc/water system are displayed in Fig.1.

    Fig.1
                            The super cell of RDX and amorphous cells of the EtOAc and EtOAc/water

    Fig.1 The super cell of RDX and amorphous cells of the EtOAc and EtOAc/water

    With COMPASS force field, a MD simulation of 200 psis was conducted under NPT (constant number of particles, pressure, and temperature) ensemble to equilibrate the AC structure of pure EtOAc and EtOAc/water binary system, and then another MD simulation of 200 psis was performed under NVT ensemble to calculate the solubility parameter. The periodic boundary condition was applied to all the simulations. Andersen thermostat[13] and Berendsen barostat[14] were selected to control the temperature and pressure, respectively. The Verlet velocity time integration method[15] with a time step of 1 fs was adopted to integrate the Newtonian equation of motion. The electrostatic interactions were calculated by the Ewald method[16] with an accuracy of 0.0001 kcal·mol-1, and the van der Waals interactions were calculated by the atom‑based method with a cutoff distance of 15.5 Å.

  • 3 Results and discussion

  • 3.1 Determination of solubility

    The mole fraction of water in the mixed solvent was defined by Eq.(1) and the saturated mole fraction solubility of RDX in the binary solvent mixtures was calculated by Eq.(2)[17,18].

    x0=m3/M3m2/M2+m3/M3
    (1)
    x=m1/M1m1/M1+m2/M2+m3/M3
    (2)

    where m 1,m 2 and m 3 represent the mass of RDX, ethyl acetate and water, respectively, M 1,M 2 and M 3 are the molar mass of RDX, ethyl acetate and water, respectively, x 0 and x are the mole fraction of water and RDX, respectively.

  • 3.2 Correlation of the solubility values

    Since solid‑liquid equilibrium is usually not available, the correlation and prediction schemes are frequently utilized. The Apelblat equation, CNIBS/R‑K and Jouyban‑Acree models are used to correlate the experimental values. The root‑mean‑square deviation(RMSD) of each solvent is used to evaluate the fitting results of the correlation equation. RMSD is defined as Eq.(3)[19,20]:

    RMSD=1Ni=1Nxcal-xexp212
    (3)
  • 3.2.1 Modified Apelblat model

    The relationship between the mole fraction solubility and temperatureis generally modeled by Apelblat equation[21,22,23,24] deduced from the Clausius‑ Clapeyron equation, a semi‑empirical equation, which can describe the solid‑liquid equilibrium precisely. The equation can be expressed as Eq.(4):

    lnx=A+BT+ClnT
    (4)

    where A, B and C are the model parameters, and T denotes the absolute temperature. The constants A and B represent the variation in the solution activity coefficient and provide an indication of the influence of non‑ideal solution on the solubility of solute; the parameter C reflects the effect of temperature on the enthalpy of fusion[25].

    The correlated solubility curves are presented in Fig.2. The model parameters, correlation coefficient (R 2) together with RMSD of Apelblat equation are displayed in Table 2. The Apelblat equation obtained from experimental results describes the relationship between the solubility of RDX and the temperature precisely. The R 2 is close to 1 with atiny RMSD, resulting in that each fitting curve passes through all the experimental points.

    Fig.2
                            The mole fraction solubility of RDX in ethyl acetate‑ water binary mixed solvents correlated by Apelblat equation

    Fig.2 The mole fraction solubility of RDX in ethyl acetate‑ water binary mixed solvents correlated by Apelblat equation

    Table 2 Model parameters A, B, C, R2 and RMSD of Apelblat equation

    x 0 A B C R 2 109×RMSD
    0.0000-158.35745253.060023.79710.99688.8363
    0.0239-84.27041738.241512.86220.99735.6023
    0.0467-85.41101798.676713.03230.99765.5914
    0.0684-87.66601894.021413.37820.99922.2630
    0.0892-91.15981986.310513.93990.99961.2730
    0.1090-88.39481828.099613.55381.00000.1813
    0.1280-81.33461591.404212.46690.99932.8186
  • 3.2.2 CNIBS/R‑K model

    The CNIBS/R‑K model is one of the theoretical models used for correlating the relationship between the solute solubility and the concentration of the binary solvents[26,27] as Eq.(5):

    lnx=B0+B1x0+B2x02+B3x03+B4x04
    (5)

    where B 0, B 1, B 2, B 3 and B 4 are the model parameters.

    The correlated solubility curves are presented in Fig.3, respectively. The model parameters,R 2, together with RMSD of CNIBS/R‑K model are listed in Table 3. The fitting curves pass through every experimental point, which means that the CNIBS/R‑K model obtained for experimental results describes the relationship between the solubility of RDX and the concentration of the binary solvents precisely.

    Fig.3
                            The mole fraction solubility of RDX in ethyl acetate‑water binary mixed solvents correlated by CNIBS/R‑K equation

    Fig.3 The mole fraction solubility of RDX in ethyl acetate‑water binary mixed solvents correlated by CNIBS/R‑K equation

    Table 3 Model parameters B0-B4, R2 and RMSD of CNIBS/R‑K model

    T / K B 0 B 1 B 2 B 3 B 4 R 2 1011×RMSD
    298.15-5.1807-1.5317110.3606-1437.90656116.91440.995318.8180
    303.15-5.05750.116745.3821-635.05692925.86350.996015.7598
    308.15-4.92400.322945.8469-665.12022944.65741.00000.0040
    313.15-4.8230-0.486457.4424-649.11242589.37140.99934.4262
    318.15-4.69471.1474-24.3906400.4854-1478.25630.99925.2310
    323.15-4.5906-0.303942.9317-380.56941297.39670.99919.9144
    328.15-4.50580.77919.241162.4654-499.22310.997940.7612
    333.15-4.36090.57475.2307132.9940-819.16920.99991.7660
    338.15-4.23600.8955-9.4206316.0742-1533.36320.99992.1997
  • 3.2.3 Jouyban‑Acree Model

    Jouyban‑Acree model is one of the versatile models to describe the solubility on both solvent compositions and temperature for binary mixed solvents[28,29,30]. The expression of this model can be described by Eq.(6):

    lnx=x0lnx0+x1lnx1+x0x1i=1NJix0-x1iT
    (6)

    where Ji is a model constant, T is the absolute temperature, and x 0 and x 1 represent the initial mole fractions of compositions of the binary solvent. When N=2, Eq.(6) can be simplified as Eq.(7):

    lnx=A0+A11T+A2lnT+A3x0+A4x01T+A5x021T+A6x031T+A7x041T+A8x0lnT
    (7)

    where A 0-A 8 are empirical model parameters which can be obtained by least‑squares analysis. The model parameter values of A 0-A 8, R 2 and RMSD are displayed in Table 4. The three‑dimensional diagram of x, x 0 and T is shown in Fig.4. The points represent the experimentalvalues, and the surface represents the results fitted by Jouyban‑Acree model. All the experimental values are almost on the surface, indicating that the experimental results are fitted well by the Jouyban‑Acree model. In addition, the values of R 2 are very close to 1 and those of RMSD are tiny. Therefore, the Jouyban‑Acree model is a suitable equation to correlate the experimental solubility values of RDX in ethyl acetate‑water system. Also, the solubility of RDX in ethyl acetate‑water system at a random temperature or concentration can be calculated by the simplified Jouyban‑Acree model obtained from this study.

    Fig.4
                            The mole fraction solubility of RDX in ethyl acetate‑water binary mixed solvents correlated by Jouyban‑Acree model

    Fig.4 The mole fraction solubility of RDX in ethyl acetate‑water binary mixed solvents correlated by Jouyban‑Acree model

    Table 4 Model parameters A0-A8, R2 and RMSD of Jouyban‑Acree model

    parametersvaluesparametersvalues
    A 0 -110.1876 A 6 -10000.8545
    A 1 479.4046 A 7 6578.7368
    A 2 16.6790 A 8 3.4612
    A 3 -19.4866 R 2 0.9989
    A 4 149.8984104×RMSD1.0516
    A 5 4113.8900
  • 3.3 Thermodynamic properties of the solution

    Some thermodynamic properties such as the standard enthalpy of dissolution (Δdissln H, kJ·mol-1), standard entropy of dissolution (Δdissln S, J·K-1·mol-1) and the standard Gibbs free energy(Δdissln G, kJ·mol-1) can be calculated when the solubility of RDX in ethyl acetate‑water mixed solvents at different temperatures is confirmed. According to the Van′t Hoff analysis, the apparent enthalpy change of solution can be related to the temperature and the solubility as Eq.(8)[31]:

    ΔdisslnHR=-lnx1/TP
    (8)

    Over a limited temperature interval (298.15 K to 338.15 K), the heat capacity change of solution may be assumed to be constant. Hence, the values ofH and S would be valid for the mean temperature,T mean(318.15 K). Thus, combined with the Apelblat model, the Δdissln H, Δdissln S, Δdissln G can be calculated by Eq.(9)-Eq.(11), respectively[32].

    ΔdisslnH=-RB-CTmean
    (9)
    ΔdisslnS=RlnxlnT+lnx=RA+C1+lnTmean
    (10)
    ΔdisslnG=ΔdisslnH-TmeanΔdisslnS
    (11)

    The results of the standard Gibbs energy, enthalpy, of dissolution and entropy of dissolution are shown in Table 5, together with %ξ H and %ξ S. %ξ H and %ξ S represent the relative contribution to the standard Gibbs energy made by enthalpy and entropy in the dissolution process[33,34], as Eq.(12)-Eq.(13).

    Table 5 Standard enthalpy of dissolution (ΔdisslnH), entropy of dissolution (ΔdisslnS), and Gibbs energy (ΔdisslnG) at the mean temperature (318.15 K) together with %ξH and %ξS

    x 0

    Δdissln H

    /kJ·mol-1

    Δdissln S

    /J·K-1·mol-1

    Δdissln G

    /kJ·mol-1

    %ξ H %ξ S
    0.000019.271821.376612.470873.9226.08
    0.023919.570122.535312.400573.1926.81
    0.046719.517522.616112.322273.0626.94
    0.068419.639823.315912.221872.5827.42
    0.089220.358325.849612.134271.2328.77
    0.109020.652327.129512.021170.5329.47
    0.128019.745224.717711.881371.5228.48
    %ξH=ΔdisslnHΔdisslnH+TΔdisslnS×100
    (12)
    %ξS=TΔdisslnSΔdisslnH+TΔdisslnS×100
    (13)

    Table 5 shows that the values of Δdissln H are positive in the binary solvent mixtures, indicating that the dissolution of RDX is an endothermic process. What's more, %ξ H is larger than %ξ S for each solvent, indicating that the main contributor to the standard molar Gibbs energy of dissolution is the enthalpy rather than the entropy.

  • 3.4 Computational parameters

    The solubility parameter (δ) of RDX and EtOAc/water is calculated by MD method to research why RDX solubility increases with the increasing water ratio of EtOAc /water binary system. The calculated solubility parameter of RDX and EtOAc /water is displayed in Table 6.

    Table 6 The solubility parameters of RDX and solvent

    RDXsolvent
    0.00000.02390.04670.06840.08920.10900.1280
    solubility parameter/(J·cm-3)1/2 30.20818.31718.51118.58818.68118.78419.00919.330

    The solubility parameters of RDX and ethyl acetate are 30.208 (J·cm-3)1/2 and 18.317 (J·cm-3)1/2 respectively. The simulation results agree with the reported results of 31.8 (J·cm-3)1/2 and 18.6(J·cm-3)1/2[12], demonstrating that this method of computation is valid and reliable. From the calculation results we can deduce that with the increase of water, the solubility parameter of solvent is increasing and the difference of solubility parameter between RDX and solvent is decreasing. This phenomenon demonstrates that RDX solubility should increase with the mole fraction water in ethyl acetate‑water binary system theoretically, which has a good agreement with the experimental results.

  • 4 Conclusions

    (1) The solubility values of RDX increase with an increase of water concentration and temperature as a nonlinear function.

    (2) Three equations (Apelblat, CNIBS/R‑K, and Jouyban‑Acree) are used for the correlation of the experimental data, and all the models agree well with the experimental values.

    (3) Some important thermodynamic properties such as the standard enthalpy of dissolution, standard entropy of dissolution and standard Gibbs free energy of dissolution have been calculated. It can be deduced that the dissolution process of RDX in ethyl acetate‑water binary system is endothermic because of the positive value of Δdissln H.

    (4) The solubility parameter is calculated by MD simulations to verify the reliability of the experimental results. The difference of the solubility parameter between RDX and mixed solvents decreases with an increase of the mole fraction of water. So, RDX prefer to dissolve into ethyl acetate‑water binary system with increasing the mole fraction of water, which reflects a good consistency with the experimental results.

    (责编:高 毅)

  • References

    • 1

      CAO Xin‑mao, LI Fu‑ping . HMX and its application[M]. Beijing:Weapon Industry Press,1993:104-106.

    • 2

      CAO Duan‑lin, WANG Jian‑long . A study on the recovery of small quantity of RDX in the producing of HMX with acetic anhydride[J]. Journal of North China Institute of Technology, 1996, 17(3): 279-282.

    • 3

      YE Ling, WU Yu‑zhong, LI Zhi‑hua . Study on transformation of α‑HMX into β‑HMX with spent acetic acid[J]. Chinese Journal of Explosives & Propellants, 2000, 23(1): 38-39, 42.

    • 4

      LI Qiao‑ling, YE Yu‑peng . A novel process for the refining of HMX[J]. Acta Armamentarii, 2002, 23(4): 555-557.

    • 5

      ZHAO Rui‑xian, LI Guang‑ming, GAO Tian‑ping, et al . Study on RDX granulate stage process[J]. Chinese Journal of Explosives & Propellants, 2003, 26(4): 67-70.

    • 6

      ZHANG Min‑hua, DOU Mao‑bin, WANG Meng‑yan,et al . Study on the solubility parameter of supercritical carbon dioxide system by molecular dynamics simulation[J]. Journal of Molecular Liquids, 2017, 248: 322-329.

    • 7

      LUO Yan‑long, WANG Run‑guo, WANG Wei, et al . Molecular dynamics simulation insight into two‑component solubility parameters of graphene and thermodynamic compatibility of graphene and styrene butadiene rubber[J]. Journal of Physical Chemistry C,2017, 121(18): 10163-10173.

    • 8

      Goharshadi E K, Akhlamadi G, Mahdizaeh S J . Investigation of graphene oxide nanosheets dispersion in water based on solubility parameters: a molecular dynamics simulation study[J]. RSC Advance,2015, 5(129): 106421-106430.

    • 9

      Leach A R . Molecular modelling: principles and applications[M]. Beijing: World Book Inc, 2003.

    • 10

      Choi C S, Bennema P . The crystal structure of cyclotrimethylenetrinitramine[J]. Acta Crystallogr B, 1972, 28: 2857-2862.

    • 11

      Barton A F M . Solubility parameters[J]. Chemical Reviews, 1975, 75(6): 731-753.

    • 12

      SUN Ye‑bin, HUI Jun‑ming, CAOXin‑mao . Military Use Blended Explosives[M]. Beijing: Weapon Industry Press, 1995: 242-251.

    • 13

      LI Jing, JIN Shao‑hua, LAN Guan‑chao, et al . Morphology control of 3‑nitro‑1,2,4‑triazole‑5‑one (NTO) by molecular dynamics simulation[J]. Crystengcomm, 20(40): 6252-6260.

    • 14

      LAN Guan‑chao, JIN Shao‑hua, LI Jing, et al .The study of external growth environments on the crystal morphology of epsilon‑HNIW by molecular dynamics simulation[J]. Journal of Materials Science,2018, 53(18): 12921-12936.

    • 15

      LI Jing, JIN Shao‑hua, LAN Guan‑chao, et al . Molecular dynamics simulations on miscibility, glass transition temperature and mechanical properties of PMMA/DBP binary system[J]. Journal of Molecular Graphics & Modelling,2018, 84: 182-188.

    • 16

      LAN Guan‑chao, JIN Shao‑hua, LI Jing, et al . Miscibility, glass transition temperature and mechanical properties of NC/DBP binary systems by molecular dynamics[J]. Propellants, Explosives, Pyrotechnics,2018, 43(6): 559-567.

    • 17

      Almarri F, Haq N, Alanazi F K, et al . Solubilityand thermodynamic function of vitamin D3 in different mono solvent[J].Journal of molecular liquids,2017, 229: 477-481.

    • 18

      Shakeel F, Imran M, Abida, Haq N, et al . Solubility and thermodynamic/solvation behavior of 6‑phenyl‑4,5‑dihydropyridazin‑3(2H)‑one in different(Transcutol‑water) mixtures[J].Journal of molecular liquids,2017, 230: 511-517.

    • 19

      CHEN Li‑zhen, ZHANG Tian‑bei, LI Man, et al . Solubility of 2,2',4,4',6,6'‑hexanitrostilbene in binary solvent of N, N‑dimethylformamide and acetonitrile[J]. Journal of Chemical Thermodynamics, 2016, 99: 99-104.

    • 20

      LAN Guan‑chao, WANG Jian‑long, CHEN Li‑zhen, et al . Measurement and correlation of the solubility of 3,4‑bis(3‑nitrofurazan‑4‑yl)furoxan (DNTF) in different solvents[J].Journal of Chemical Thermodynamics, 2015, 89: 264-269.

    • 21

      Apelblat A, Manzurola E .Solubilities of manganese, cadmium, mercury and lead acetates in water from T=278.15 K to T=340.15 K[J]. Journal of Chemical Thermodynamics,2001, 33(2): 147-153.

    • 22

      Apelblat A, Manzurola E . Solubilities of L‑aspartic, DL‑aspartic, DL‑glutamic, p‑hydroxybenzoic, o‑anisic, p‑anisic, and itaconic acids in water from T=278 K to T=345 K[J]. Journal of Chemical Thermodynamics,1997, 29(12): 1527-1533.

    • 23

      Apelblat A, Manzurola E . Solubilities of o‑acetylsalicylic, 4‑aminosalicylic, 3,5‑dinitrosalicylic, and p‑toluic acid, and magnesium‑DL‑aspartate in water from T=(278 to 348) K[J]. Journal of Chemical Thermodynamics,1999, 31(1): 85-91.

    • 24

      ZHOU Li, ZHANG Pei‑pei, YANG Guang‑de, et al .Solubility of chrysin in ethanol and water mixtures[J]. Journal of Chemical and Engineering Data,2014, 59(7): 2215-2220.

    • 25

      CHEN Zhao‑guo, YANG Wen‑ge, HU Yong‑hong, et al . Measurement and correlation for the solubility of dimethyl 1,4‑cyclohexanedione‑2,5‑dicarboxylate in different solvents at temperatures from (278.15 to 323.15) K[J]. Journal of Chemical and Engineering Data,2011, 56(5):2726-2729.

    • 26

      AcreeJr W E .Comments concerning ‘model for solubility estimation in mixed solvent system’[J].International Journal of Pharmaceutics,1996, 127(1): 27-30.

    • 27

      Jouyban‑Gharamaleki A, Hanaee J . A general model from theoretical cosolvency models[J].International Journal of Pharmaceutics,1997, 152(2): 247-250.

    • 28

      LIU Bao‑shu, SUN Hua, WANG Jing‑kang, et al . Solubility of disodium 5′‑guanylate heptahydrate in aqueous methanol mixtures[J].Food Chemistry,2011,128(1): 218-221.

    • 29

      Jouyban A . Review of the cosolvency models for predicting solubility of drugs in water‑cosolvent mixtures[J].Journal of Pharmacy and Pharmaceutical Sciences,2008, 11(1): 32-57.

    • 30

      WANG Shui,QIN Li‑ying, ZHOU Zhi‑mao, et al . Solubility and solution thermodynamics of betaine in different pure solvents and binary mixtures[J]. Journal of Chemical and Engineering Data,2012,57(8): 2128-2135.

    • 31

      WANG Guan, WANG Yong‑li, MA You‑guang, et al . Determination and correlation of cefuroxime acid solubility in (acetonitrile‑water) mixtures[J].Journal of Chemical Thermodynamics, 2014, 77: 144-150.

    • 32

      ZHANG Hui, YIN Qiu‑xiang, LIU Zeng‑kun, et al . Measurement and correlation of solubility of dodecanedioic acid in different pure solvents from T=(288.15 to 323.15) K[J]. Journal of Chemical Thermodynamics,2014, 68: 270-274.

    • 33

      Delgado D R, Holguin A R, Almanza O, et al . Solubility and preferential solvation of meloxicam in ethanol‑water mixtures[J]. Fluid Phase Equilibria,2011, 305(1): 80-85.

    • 34

      Holguin A R, Delgado D R, Martinez F, et al . Solution thermodynamics and preferential solvation of meloxicam in propylene glycol‑water mixtures[J]. Journal of Solution Chemistry, 2011, 40(12):1987-1999.

LIJing

机 构: 北京理工大学材料学院, 北京 100081

Affiliation: School of Materials Science and Engineering, Beijing Institute of Technology, Beijing 100081, China

邮 箱:360018908@qq.com

Profile: LI Jing(1989-), female, doctoral candidate, majoring in the study of energetic materials. e‑mail:360018908@qq.com

JINShao‑hua

机 构: 北京理工大学材料学院, 北京 100081

Affiliation: School of Materials Science and Engineering, Beijing Institute of Technology, Beijing 100081, China

XUZi‑shuai

机 构: 甘肃银光化学工业集团有限公司, 甘肃 白银 730900

Affiliation: Research Institute of Gansu Yin Guang Chemical Industry Group Co.Ltd, Baiyin 730900, China

WUNa‑na

机 构: 甘肃银光化学工业集团有限公司, 甘肃 白银 730900

Affiliation: Research Institute of Gansu Yin Guang Chemical Industry Group Co.Ltd, Baiyin 730900, China

LANGuan‑chao

机 构: 北京理工大学材料学院, 北京 100081

Affiliation: School of Materials Science and Engineering, Beijing Institute of Technology, Beijing 100081, China

CHENShu‑sen

机 构: 北京理工大学材料学院, 北京 100081

Affiliation: School of Materials Science and Engineering, Beijing Institute of Technology, Beijing 100081, China

WANGDong‑xu

机 构: 北京理工大学机电学院, 北京 100081

Affiliation: School of Mechatronical Engineering, Beijing Institute of Technology, Beijing 100081, China

角 色:通讯作者

Role:Corresponding author

邮 箱:dadadolindsay@126.com

Profile: WANG Dong‑xu(1987-), male, doctor, majoring in the study of energetic materials. e‑mail:dadadolindsay@126.com

lattice

parameter

experimental

value[10]

optimized value

relative

error/%

a13.18213.450-2.033
b11.57411.2832.514
c10.70910.2254.520
β/(°)90.00090.0000.000
html/hncl/CJEM2018194/alternativeImage/eebf16ab-83bb-4ddd-b8ac-af2f89be9009-F001.jpg
html/hncl/CJEM2018194/alternativeImage/eebf16ab-83bb-4ddd-b8ac-af2f89be9009-F002.jpg
x 0 A B C R 2 109×RMSD
0.0000-158.35745253.060023.79710.99688.8363
0.0239-84.27041738.241512.86220.99735.6023
0.0467-85.41101798.676713.03230.99765.5914
0.0684-87.66601894.021413.37820.99922.2630
0.0892-91.15981986.310513.93990.99961.2730
0.1090-88.39481828.099613.55381.00000.1813
0.1280-81.33461591.404212.46690.99932.8186
html/hncl/CJEM2018194/alternativeImage/eebf16ab-83bb-4ddd-b8ac-af2f89be9009-F003.jpg
T / K B 0 B 1 B 2 B 3 B 4 R 2 1011×RMSD
298.15-5.1807-1.5317110.3606-1437.90656116.91440.995318.8180
303.15-5.05750.116745.3821-635.05692925.86350.996015.7598
308.15-4.92400.322945.8469-665.12022944.65741.00000.0040
313.15-4.8230-0.486457.4424-649.11242589.37140.99934.4262
318.15-4.69471.1474-24.3906400.4854-1478.25630.99925.2310
323.15-4.5906-0.303942.9317-380.56941297.39670.99919.9144
328.15-4.50580.77919.241162.4654-499.22310.997940.7612
333.15-4.36090.57475.2307132.9940-819.16920.99991.7660
338.15-4.23600.8955-9.4206316.0742-1533.36320.99992.1997
html/hncl/CJEM2018194/alternativeImage/eebf16ab-83bb-4ddd-b8ac-af2f89be9009-F004.jpg
parametersvaluesparametersvalues
A 0 -110.1876 A 6 -10000.8545
A 1 479.4046 A 7 6578.7368
A 2 16.6790 A 8 3.4612
A 3 -19.4866 R 2 0.9989
A 4 149.8984104×RMSD1.0516
A 5 4113.8900
x 0

Δdissln H

/kJ·mol-1

Δdissln S

/J·K-1·mol-1

Δdissln G

/kJ·mol-1

%ξ H %ξ S
0.000019.271821.376612.470873.9226.08
0.023919.570122.535312.400573.1926.81
0.046719.517522.616112.322273.0626.94
0.068419.639823.315912.221872.5827.42
0.089220.358325.849612.134271.2328.77
0.109020.652327.129512.021170.5329.47
0.128019.745224.717711.881371.5228.48
RDXsolvent
0.00000.02390.04670.06840.08920.10900.1280
solubility parameter/(J·cm-3)1/2 30.20818.31718.51118.58818.68118.78419.00919.330

Table 1 Comparison of the experimental lattice parameters and the optimized values of RDX

Fig.1 The super cell of RDX and amorphous cells of the EtOAc and EtOAc/water

Fig.2 The mole fraction solubility of RDX in ethyl acetate‑ water binary mixed solvents correlated by Apelblat equation

Table 2 Model parameters A, B, C, R2 and RMSD of Apelblat equation

Fig.3 The mole fraction solubility of RDX in ethyl acetate‑water binary mixed solvents correlated by CNIBS/R‑K equation

Table 3 Model parameters B0-B4, R2 and RMSD of CNIBS/R‑K model

Fig.4 The mole fraction solubility of RDX in ethyl acetate‑water binary mixed solvents correlated by Jouyban‑Acree model

Table 4 Model parameters A0-A8, R2 and RMSD of Jouyban‑Acree model

Table 5 Standard enthalpy of dissolution (ΔdisslnH), entropy of dissolution (ΔdisslnS), and Gibbs energy (ΔdisslnG) at the mean temperature (318.15 K) together with %ξH and %ξS

Table 6 The solubility parameters of RDX and solvent

image /

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

  • References

    • 1

      CAO Xin‑mao, LI Fu‑ping . HMX and its application[M]. Beijing:Weapon Industry Press,1993:104-106.

    • 2

      CAO Duan‑lin, WANG Jian‑long . A study on the recovery of small quantity of RDX in the producing of HMX with acetic anhydride[J]. Journal of North China Institute of Technology, 1996, 17(3): 279-282.

    • 3

      YE Ling, WU Yu‑zhong, LI Zhi‑hua . Study on transformation of α‑HMX into β‑HMX with spent acetic acid[J]. Chinese Journal of Explosives & Propellants, 2000, 23(1): 38-39, 42.

    • 4

      LI Qiao‑ling, YE Yu‑peng . A novel process for the refining of HMX[J]. Acta Armamentarii, 2002, 23(4): 555-557.

    • 5

      ZHAO Rui‑xian, LI Guang‑ming, GAO Tian‑ping, et al . Study on RDX granulate stage process[J]. Chinese Journal of Explosives & Propellants, 2003, 26(4): 67-70.

    • 6

      ZHANG Min‑hua, DOU Mao‑bin, WANG Meng‑yan,et al . Study on the solubility parameter of supercritical carbon dioxide system by molecular dynamics simulation[J]. Journal of Molecular Liquids, 2017, 248: 322-329.

    • 7

      LUO Yan‑long, WANG Run‑guo, WANG Wei, et al . Molecular dynamics simulation insight into two‑component solubility parameters of graphene and thermodynamic compatibility of graphene and styrene butadiene rubber[J]. Journal of Physical Chemistry C,2017, 121(18): 10163-10173.

    • 8

      Goharshadi E K, Akhlamadi G, Mahdizaeh S J . Investigation of graphene oxide nanosheets dispersion in water based on solubility parameters: a molecular dynamics simulation study[J]. RSC Advance,2015, 5(129): 106421-106430.

    • 9

      Leach A R . Molecular modelling: principles and applications[M]. Beijing: World Book Inc, 2003.

    • 10

      Choi C S, Bennema P . The crystal structure of cyclotrimethylenetrinitramine[J]. Acta Crystallogr B, 1972, 28: 2857-2862.

    • 11

      Barton A F M . Solubility parameters[J]. Chemical Reviews, 1975, 75(6): 731-753.

    • 12

      SUN Ye‑bin, HUI Jun‑ming, CAOXin‑mao . Military Use Blended Explosives[M]. Beijing: Weapon Industry Press, 1995: 242-251.

    • 13

      LI Jing, JIN Shao‑hua, LAN Guan‑chao, et al . Morphology control of 3‑nitro‑1,2,4‑triazole‑5‑one (NTO) by molecular dynamics simulation[J]. Crystengcomm, 20(40): 6252-6260.

    • 14

      LAN Guan‑chao, JIN Shao‑hua, LI Jing, et al .The study of external growth environments on the crystal morphology of epsilon‑HNIW by molecular dynamics simulation[J]. Journal of Materials Science,2018, 53(18): 12921-12936.

    • 15

      LI Jing, JIN Shao‑hua, LAN Guan‑chao, et al . Molecular dynamics simulations on miscibility, glass transition temperature and mechanical properties of PMMA/DBP binary system[J]. Journal of Molecular Graphics & Modelling,2018, 84: 182-188.

    • 16

      LAN Guan‑chao, JIN Shao‑hua, LI Jing, et al . Miscibility, glass transition temperature and mechanical properties of NC/DBP binary systems by molecular dynamics[J]. Propellants, Explosives, Pyrotechnics,2018, 43(6): 559-567.

    • 17

      Almarri F, Haq N, Alanazi F K, et al . Solubilityand thermodynamic function of vitamin D3 in different mono solvent[J].Journal of molecular liquids,2017, 229: 477-481.

    • 18

      Shakeel F, Imran M, Abida, Haq N, et al . Solubility and thermodynamic/solvation behavior of 6‑phenyl‑4,5‑dihydropyridazin‑3(2H)‑one in different(Transcutol‑water) mixtures[J].Journal of molecular liquids,2017, 230: 511-517.

    • 19

      CHEN Li‑zhen, ZHANG Tian‑bei, LI Man, et al . Solubility of 2,2',4,4',6,6'‑hexanitrostilbene in binary solvent of N, N‑dimethylformamide and acetonitrile[J]. Journal of Chemical Thermodynamics, 2016, 99: 99-104.

    • 20

      LAN Guan‑chao, WANG Jian‑long, CHEN Li‑zhen, et al . Measurement and correlation of the solubility of 3,4‑bis(3‑nitrofurazan‑4‑yl)furoxan (DNTF) in different solvents[J].Journal of Chemical Thermodynamics, 2015, 89: 264-269.

    • 21

      Apelblat A, Manzurola E .Solubilities of manganese, cadmium, mercury and lead acetates in water from T=278.15 K to T=340.15 K[J]. Journal of Chemical Thermodynamics,2001, 33(2): 147-153.

    • 22

      Apelblat A, Manzurola E . Solubilities of L‑aspartic, DL‑aspartic, DL‑glutamic, p‑hydroxybenzoic, o‑anisic, p‑anisic, and itaconic acids in water from T=278 K to T=345 K[J]. Journal of Chemical Thermodynamics,1997, 29(12): 1527-1533.

    • 23

      Apelblat A, Manzurola E . Solubilities of o‑acetylsalicylic, 4‑aminosalicylic, 3,5‑dinitrosalicylic, and p‑toluic acid, and magnesium‑DL‑aspartate in water from T=(278 to 348) K[J]. Journal of Chemical Thermodynamics,1999, 31(1): 85-91.

    • 24

      ZHOU Li, ZHANG Pei‑pei, YANG Guang‑de, et al .Solubility of chrysin in ethanol and water mixtures[J]. Journal of Chemical and Engineering Data,2014, 59(7): 2215-2220.

    • 25

      CHEN Zhao‑guo, YANG Wen‑ge, HU Yong‑hong, et al . Measurement and correlation for the solubility of dimethyl 1,4‑cyclohexanedione‑2,5‑dicarboxylate in different solvents at temperatures from (278.15 to 323.15) K[J]. Journal of Chemical and Engineering Data,2011, 56(5):2726-2729.

    • 26

      AcreeJr W E .Comments concerning ‘model for solubility estimation in mixed solvent system’[J].International Journal of Pharmaceutics,1996, 127(1): 27-30.

    • 27

      Jouyban‑Gharamaleki A, Hanaee J . A general model from theoretical cosolvency models[J].International Journal of Pharmaceutics,1997, 152(2): 247-250.

    • 28

      LIU Bao‑shu, SUN Hua, WANG Jing‑kang, et al . Solubility of disodium 5′‑guanylate heptahydrate in aqueous methanol mixtures[J].Food Chemistry,2011,128(1): 218-221.

    • 29

      Jouyban A . Review of the cosolvency models for predicting solubility of drugs in water‑cosolvent mixtures[J].Journal of Pharmacy and Pharmaceutical Sciences,2008, 11(1): 32-57.

    • 30

      WANG Shui,QIN Li‑ying, ZHOU Zhi‑mao, et al . Solubility and solution thermodynamics of betaine in different pure solvents and binary mixtures[J]. Journal of Chemical and Engineering Data,2012,57(8): 2128-2135.

    • 31

      WANG Guan, WANG Yong‑li, MA You‑guang, et al . Determination and correlation of cefuroxime acid solubility in (acetonitrile‑water) mixtures[J].Journal of Chemical Thermodynamics, 2014, 77: 144-150.

    • 32

      ZHANG Hui, YIN Qiu‑xiang, LIU Zeng‑kun, et al . Measurement and correlation of solubility of dodecanedioic acid in different pure solvents from T=(288.15 to 323.15) K[J]. Journal of Chemical Thermodynamics,2014, 68: 270-274.

    • 33

      Delgado D R, Holguin A R, Almanza O, et al . Solubility and preferential solvation of meloxicam in ethanol‑water mixtures[J]. Fluid Phase Equilibria,2011, 305(1): 80-85.

    • 34

      Holguin A R, Delgado D R, Martinez F, et al . Solution thermodynamics and preferential solvation of meloxicam in propylene glycol‑water mixtures[J]. Journal of Solution Chemistry, 2011, 40(12):1987-1999.