10.8
CiteScore
 
5.3
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
Corrigendum
Current Issue
Editorial
Erratum
Full Length Article
Full lenth article
Original Article
Research article
Retraction notice
Review
Review Article
SPECIAL ISSUE: ENVIRONMENTAL CHEMISTRY
10.8
CiteScore
5.3
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
Corrigendum
Current Issue
Editorial
Erratum
Full Length Article
Full lenth article
Original Article
Research article
Retraction notice
Review
Review Article
SPECIAL ISSUE: ENVIRONMENTAL CHEMISTRY
View/Download PDF

Translate this page into:

Original article
01 2022
:16;
104421
doi:
10.1016/j.arabjc.2022.104421

Determination of the solubility of rivaroxaban (anticoagulant drug, for the treatment and prevention of blood clotting) in supercritical carbon dioxide: Experimental data and correlations

Department of Chemical Engineering, Faculty of Engineering, University of Kashan, 87317-53153 Kashan, Iran
Laboratory of Supercritical Fluids and Nanotechnology, University of Kashan, 87317-53153 Kashan, Iran
Modeling and Simulation Centre, Faculty of Engineering, University of Kashan, 87317-53153 Kashan, Iran

⁎Corresponding author. sodeifian@kashanu.ac.ir (Gholamhossein Sodeifian)

Disclaimer:
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.

Peer review under responsibility of King Saud University.

Abstract

Abstract

Supercritical processes are utilized in pharmaceutical, chemical, agricultural and food industries. To design the micro or nanoparticles formation processes via a supercritical fluid (SCF), it is necessary to have available the solubility data of a solid solute in a solvent such as supercritical carbon dioxide (SC- C O 2 ). In this article, the solubility of rivaroxaban (anticoagulant drug, for the treatment and prevention of blood clotting) in SC- C O 2 was investigated for the first time. Six semi-empirical models (Chrastil, Bartle, Jouyban, Garlapati and Madras, and Sodeifian) and two equation of states (EoSs) including Peng-Robinson (PR) and Soave-Redlich-Kwong (SRK) with twin-parametric van der Waals (vdW2) and covolume dependent (CVD) mixing rules were used to examine the correlation of drug solubility data. The best performance was obtained by semi-empirical Jouyban model (Average absolute relative deviation (AARD = 12.436 %)) and PR-vdW2 EoS (AARD = 5.104 %). Total and vaporization enthalpies were calculated 43.617 KJ.mol−1 and 63.288 KJ.mol−1, respectively.

Keywords

Supercritical carbon dioxide
Rivaroxaban
Solubility
Peng-Robinson
Sodeifian model
Soave-Redlich-Kwong

Nomenclature

a0-a7

Parameters for density-based models

a, b

EoS constants

AARD %

Percent of average absolute relative deviation

C s

Concentration of drug inside collection vial (µg/mL)

CVD

Covolume dependent mixing rule

H sol

Rivaroxaban solvation heat in kJ/mol

H sub

Rivaroxaban vaporization heat in kJ/mol

H tot

Total enthalpy of rivaroxaban in kJ/mol

K 12 & l 12

Interaction parameters of Eqs. (19) and (20)

Mij

Interaction parameters of Eqs. (27) and (28)

M C O 2 , Ms

Molar mass of CO2, and drug in g/mol

M S R

Mean square regression

M S E

Mean square residual

n C O 2

Mole of CO2

n solute

Mole of rivaroxaban drug

N

The number of experimental data

NIST

National Institute of Standards and Technology

P

Pressure (bar)

Pc

Critical pressure (bar)

Pref

Reference pressure (0.1 MPa)

P sub

Pure solid sublimation pressure (Pa)

PR

Peng − Robinson

Q

Number of independent variables in each equation

R

Universal gas constant in J/(mol K)

R 2

Correlation coefficient

S

Solubility of rivaroxaban in equilibrium state (g/L)

SRK

Soave-Redlich-Kwong

SSE

Sum of squares error

S S R

Regression sum of squares

S S T

Total sum of squares

SCF

Supercritical fluid

SC-CO2

Supercritical carbon dioxide

T

Temperature (K)

T b

Boiling temperature (K)

T r

Reduced temperature

T c

Critical temperature (K)

U

Expanded uncertainty

V s

Molar volume of rivaroxaban drug (cm3/mol)

V S ( L )

Volumes of the collection vial

V l ( L )

Volume of the sampling loop

vdW2

Twin-parametric van der Waals mixing rule

y2

Equilibrium molar fraction

y i c a l c

Calculated values of the solute solubility

y i exp

Experimental value of the solute solubility

Greek Symbols

ρ C O 2 , ρ 1

Density of SC-CO2 (kg/m3)

ω

Acentric factor

1

1 Introduction

Rivaroxaban is an oral, direct Factor Xa inhibitor, which is at an advanced stage of clinical development for prevention and treatment of thromboembolic disorders (Mueck et al., 2011). In fact, this anticoagulant drug can be used for the treatment and prevention of blood clotting. It was also evaluated for auxiliary anticipation after intense coronary disorder (Mega et al., 2012) and for thromboprophylaxis in intensely sick therapeutic patients (Cohen et al., 2011; Douxfils et al., 2012). The bioavailability of drugs after consumption depend on their solubility and solubility rate in an aqueous medium. For improvement these essential properties, the particle size reduction is an important tool (Sodeifian et al., 2020). Traditional methods for reducing particle size in the pharmaceutical industry include milling, evaporation, recrystallization, etc. (Sodeifian et al., 2017; Adachi and Lu, 1983). For this purpose, a supercritical fluid (SCF) process is useful in the pharmaceutical and biochemical industries compared to conventional methods. SCFs have properties between gas and liquid, such as viscosity, density and diffusivity, providing a wide application range for energy, solubility, polymers, extraction and etc. (Brunner, 2015; Sim Yeoh et al., 2013; Kiran, 2016; Sodeifian et al., 2020).

In general, using SCF technology, morphology, shape, and size of particles can be changed, which increases the efficiency and performance of the drug (Sodeifian et al., 2021). The common methods used in SCF technology can be mentioned RESS, SAS, GAS, RESOLVE ASES, SEDS, SFEE and PGSS (Chen et al., 2018; Montes et al., 2017; Jafari et al., 2015; Ciou et al., 2018; Lee et al., 2020; Kodama et al., 2018; Lee et al., 2018; Banchero, 2021; Aredo et al., 2021; Tokunaga et al., 2021; Razmimanesh et al., 2021). This is attributed to the supercritical conditions, as stated before, where properties of CO2 are between gas and liquid. This provides a suitable situation for CO2 solvent to achieve particle formation. Supercritical carbon dioxide (SC- C O 2 ), widely used as a solvent, is one of the most common SCFs and it has a low critical temperature (304 K) and pressure (73.8 bar). Other characteristics of SC- C O 2 include non-toxic, non-flammable, good chemical stability, inert, recyclable and low cost (Sodeifian et al., 2016; Sodeifian et al., 2018; Shamsipur et al., 2004; Sodeifian and Sajadian, 2017; Tabernero et al., 2012; Nasri, 2018; Coelho et al., 2016). To design and optimize drug processes such as micro/nanoparticles of drug particles, knowing the solubility of drug compounds in SCF is essential (Sang et al., 2017; Chen et al., 2017; Sodeifian et al., 2020).

Due to the experimental limitations like time consuming and high cost of experiments in various temperatures and pressures to measure the solubility of a drug in SC- C O 2 , the use of semi-empirical methods and mathematical modeling is valuable tools to investigate the solubility behavior (Hazaveie et al., 2020; Jouyban et al., 2005; Bartle et al., 1991). There are various methods for modeling and correlation of solubility data, which can be referred to EoSs with different mixing rules, semi-empirical models, solubility models and artificial neural networks (ANN) (Nasri, 2018; Huang et al., 2013; Sodeifian et al., 2020). For modeling solubility data using EoS, it is related to phase equilibria and thermodynamic properties. EoS-based models require a knowledge of physico-chemical properties of solid (drug). In the most cases, these properties are not available for pharmaceuticals and can be used group contribution methods to calculate them (Sodeifian et al., 2017; Farrokh-Niae et al., 2008; Yazdizadeh et al., 2011). Some of the common mixing rules are vdW1, vdW2, Kwak-Mansoori and Wong-Sandler. Furthermore, semi-empirical models are simpler than EoSs and only need temperature, pressure, and density of SC- C O 2 . These models do not require an additional information of solutes (Sim Yeoh et al., 2013; Garlapati and Madras, 2009; Jouyban et al., 2005; Hozhabr et al., 2014). In this work, semi-empirical models with different adjustable parameters from three to six such as Chrastil (Bian et al., 2016); Bian et al. (Bian et al., 2015); Khansary et al. (Asgarpour Khansary et al., 2015); Spark et al. (Sparks et al., 2008); Jouyban et al. (Jouyban et al., 2002), and Sodeifian et al. (Sodeifian et al., 2021) models were applied to compare the accuracy of models for correlating the experimental solubility data.

To the best of knowledge, there is no yet any report on the measurement and modeling of rivaroxaban drug in SC- C O 2 . Accordingly, for the first time, the solubility of rivaroxaban was measured in various temperatures (308–338 K) and pressures (120–270 bar). Then the solubility data was investigated and correlated by different EoSs (PR and SRK) (Sheikhi-Kouhsar et al., 2015; Jaubert and Privat, 2010) with twin-parametric van der Waals (vdW2) (Esmaeilzadeh et al., 2009) and covolume dependent (CVD) (Yazdizadeh et al., 2012) mixing rules; and with six semi-empirical equations based on density. The used models were evaluated by statistical analysis. Finally, total and vaporization enthalpies were estimated in SC- C O 2 for the first time.

2

2 Experimental

2.1

2.1 Materials

Methanol was purchased from Merck company (Germany). Pure CO2 (99.99 %) was provided by Fadak company (Kashan, Iran). Rivaroxaban with formulation C 19 H 18 C l N 3 O 5 S (molar mass = 435.90 g/mol) was prepared by Parsian company (Tehran, Iran). Some properties of the used materials are shown in Table 1.

Table 1 Some properties of the used materials.
Compound Structure MW (g/mol) CAS number Melting point Purity (mass%) and purification method by supplier
Rivaroxaban 435.90 366789–02-8 230 °C ≤99.99, HPLC
Carbon dioxide 44.01 124–38-9 99.99, GC
Methanol 32.04 67–56-1 99.98, GC

2.2

2.2 Experimental details

Experimental set up for measuring the solubility of rivaroxaban is shown in Fig. 1 (Sodeifian et al., 2019). Previous articles can be used to know more details (Sodeifian et al., 2021; Sodeifian et al., 2020).

Schematic of experimental process. E-1 CO2 cylinder; E-2, filter; E-3, refrigerator unit; E-4, air compressor; E-5, high pressure pump (Haskel model, MS-110, 0.33HP, air-driven liquid pump, USA); E-6, equilibrium cell; E-7, oven; E-8, six-port, two position valve (injection ring); E-9, control panel; E-10, syringe; E-11, collection vial; E-12, spectrophotometer UV–vis.
Fig. 1
Schematic of experimental process. E-1 CO2 cylinder; E-2, filter; E-3, refrigerator unit; E-4, air compressor; E-5, high pressure pump (Haskel model, MS-110, 0.33HP, air-driven liquid pump, USA); E-6, equilibrium cell; E-7, oven; E-8, six-port, two position valve (injection ring); E-9, control panel; E-10, syringe; E-11, collection vial; E-12, spectrophotometer UV–vis.

Components used include a CO2 storage cylinder, a refrigerator, a filter, a solubility column, a sample collector, and temperature and pressure measuring equipment. First, the refrigerator turned on and the temperature was set at 253 K. After the refrigerator temperature was stabilized, the CO2 in the tank was directed to the refrigerator. CO2 gas liquefied at 253 K. Liquid CO2 entered the reciprocating pump (Haskel model, MS-110, 0.33HP, air-driven liquid pump, USA) with the pressure in the tank (about 60 bar) to provide the desired pressure. Then fixed with a barometer up to ± 1 bar (U(p) = 1 bar). To maintain the experimental temperatures at the desired levels, a temperature-controlled oven with a temperature accuracy of u(T) = 0.1 K was utilized. A 70 mL equilibrium vessel (cell), containing glass beads of 2 mm in diameter and sintered filters (1 μ  m) on its both sides filled with 0.5 g of rivaroxaban drug for all experimental data, was placed in an oven to ensure that its temperature remains within ± 0.1 K. The cell was then saturated with high-pressure CO2 for 60 min to reach the equilibrium state that was determined by varying the contact times (Fig. 1S). After equilibrium, a 2-status 6-way port valve collected 600 μ L of saturated SC-CO2 and added it to a vial with 4 mL methanol. As a final step, the loop was washed with the solvent (1 mL) and the final volume of the solution was adjusted to 5 mL. Rivaroxaban absorbance in methanol was measured with a spectrophotometer, UNICO-4802 UV–Vis at λmax. The solubility tests were performed three times to determine averages. At the end of the experiments, 100 μg of the drug was dissolved in methanol. A number of standard solutions were then produced by diluting the initial solution to draw the calibration curve. The equilibrium molar fraction ( y 2 ), and solubility (S) at different temperatures and pressures in SC-CO2 were calculated as follows:

(1)
y 2 = n S o l u t e n s o l u t e + n C O 2

Two parameters n S o l u t e and n C O 2 are moles of solute (drug) and C O 2 , respectively that:

(2)
n S o l u t e = C S g L . V S L M S g m o l
(3)
n C O 2 = V l L . ρ g L M C O 2 g m o l

M S and M C O 2 are molecular weight of solute and C O 2 , respectively.

And

(4)
S g L = C S g l . V S ( L ) V l ( L )

C S is the solute concentration in the collection vial, V S ( L ) is volume of the collection vial ( 5 × 10 - 3 ), V l ( L ) is volume of the sampling loop ( 600 × 10 - 6 ). The expanded uncertainty of experimental solubility data is presented in Table 3.

3

3 Thermodynamic framework

To investigate the correlation of solubility of rivaroxaban in SC-CO2, the two EoSs (PR and SRK) with vdW2 and CVD mixing rules and six semi-empirical models based on density were used.

3.1

3.1 Semi-empirical models

In this paper, six semi-empirical were used to investigate the correlation of rivaroxaban solubility in SC-CO2. Two models of Chrastil and Belghait, relating a solute solubility to density and temperature of SCF, and other models including Sodeifian, Jouyban, Bartle and Jafari Nejad, relating solubility to temperature, pressure and density of SCF, were used. Table 2 shows the equations of these models.

Table 2 Representation of semi-empirical models equation used in this work.
Model Equation Number of parameters References
Chrastil L n S =  a 0 ln ρ + a 1 T + a 2 3 (Sridar et al., 2013)
Bartle L n ( y 2 P P r e f ) = a 0 + a 1 T + a 2 ( ρ - ρ r e f ) 3 (Reddy and Garlapati, 2019)
Jafari Nejad L n y 2 = a 0 + a 1 P 2 + a 2 T 2 + a 3 L n ρ 4 (Jafari Nejad et al., 2010)
Jouyban L n y 2 = a 0 + a 1 ρ + a 2 P 2 + a 3 PT + a 4 T P + a 5 L n ρ 6 (Si-Moussa et al., 2017)
Sodeifian L n y 2 = a 0 + a 1 P 2 T + a 2 L n ρ 1 T + a 3 ρ 1 L n ρ 1 + a 4 P L n T + a 5 L n ( ρ 1 ) T 6 (Sodeifian et al., 2019)
Belghait L n y 2 = a 0 + a 1 ρ C O 2 + a 2 ρ C O 2 2 + a 3 ρ C O 2 T + a 4 T + a 5 T 2 + a 6 L n ρ C O 2 + a 7 T 8 (Belghait et al., 2018)
Table 3 Solubility of rivaroxaban (component 2) in SC- C O 2 (component 1). The relative combined standard uncertainty and the experimental standard deviation expanded uncertainty are S ( y k ) = j = 1 n y j - y ¯ 2 n - 1 , U = k * ucombined, and.. u c o m b i n e d / y = i = 1 N P i u x i / x i 2
Temperature (K)a Pressure
(bar)a
Density of
SC- C O 2 (kg.m−3)
y 2 × 10 5 (Mole fraction) S (Equilibrium solubility) (g/L) Experimental standard deviation, S(ÿ) × (105) Expanded uncertainty of mole fraction (105 U)
308 120
150
180
210
240
270
769
817
849
875
896
914
0.393
0.542
0.697
1.079
1.200
1.603
0.0607
0.0889
0.1189
0.1895
0.2158
0.2942
0.011
0.014
0.020
0.041
0.051
0.061
0.027
0.037
0.051
0.095
0.115
0.141
318 120
150
180
210
240
270
661
744
791
824
851
872
0.304
0.440
0.554
0.973
1.316
1.731
0.0403
0.0657
0.0879
0.1610
0.2248
0.3031
0.012
0.021
0.012
0.015
0.010
0.073
0.026
0.046
0.034
0.052
0.063
0.165
328 120
150
180
210
240
270
509
656
725
769
802
829
0.203
0.345
0.455
0.770
1.490
1.886
0.0208
0.0455
0.0663
0.1189
0.2400
0.3138
0.010
0.011
0.013
0.019
0.021
0.032
0.022
0.025
0.031
0.051
0.078
0.104
338 120
150
180
210
240
270
388
557
652
710
751
783
0.104
0.268
0.365
0.679
1.635
2.062
0.0081
0.0300
0.0478
0.0969
0.2466
0.3242
0.001
0.002
0.011
0.028
0.021
0.025
0.005
0.013
0.028
0.064
0.084
0.104
Standard uncertainty u is herein set to u (T) = ±0.1 K and u (P) = ± 1 bar. The coverage factor, k = 2 corresponds to a confidence level of approximately 95 %.

3.2

3.2 Equation of state (EoS) models and mixing rules

The following formula is used to calculate the solubility of rivaroxaban in SC-CO2, assuming that SC-CO2 and rivaroxaban are compound 1 and 2, respectively:

(5)
y 2 = P 2 s u b T P φ 2 s a t φ 2 T , P , Y e x p v 2 s P - P 2 s u b T R T

P 2 s u b T is sublimation pressure of drug, φ 2 s a t is saturation fugacity coefficient of the solute, φ 2 is fugacity coefficient of the solute in supercritical phase and v 2 s is solid molar volume. If P s u b is very small, assuming that the saturation pressure of the solutes is approximately equal to one, the following equation can be used for calculating φ 2 with each EoS (Coimbra et al., 2006):

(6)
R T l n φ i = - R T l n Z + V P n i T , V , n j n i - R T V d V

Various methods can be used to estimate the physicochemical properties of drug compounds such as T c , P c , P s u b , and V s . In this study, Joback, Marrero and Pardillo, Grigoras, Fedors, Ambrose-Walton, and Lee-Kesler methods (Joback and Reid, 1987; Fedors, 1974; Poling et al., 2001) were used to calculate the critical properties. It should be noted that, using only a single EoS, the solubility of the compounds cannot be obtained. Mixing rules help us to solve this problem. In this research, PR and SRK EoSs with vdW2 and CVD mixing rules were used. PR EoS is written as follows:

(7)
P = R T V - b - a α V 2 + 2 b V - b 2

Where P, T, V and R are pressure, temperature, molar volume and gas constant, respectively. The parameters a and b calculated as follows:

(8)
a = 0.45724 R 2 T C 2 P C
(9)
b = 0.07780 R T C P C
(10)
α = ( 1 + k ( 1 - T r 1 2 ) ) 2
(11)
k = 0.37464 + 1.54226 ω - 0.26992 ω 2

Soave-Redlich-Kwong (SRK) EoS is written as follows:

(12)
P = R T V - b - a ( T ) V ( V + b )
(13)
a T = 0.42748 R 2 T C 2 P C × α ( T r , ω )
(14)
α T r , ω = [ 1 + m 1 - T r 0.5 ] 2
(15)
m = 0.480 + 1.574 ω - 0.176 ω 2
(16)
b = 0.08664 R T C P C

ω is acentric factor for each compound.

vdW2 and CVD mixing rules are explained as below:

A: vdW2 mixing rule:

(17)
a m = i j y i y j a ij
(18)
b m = i j y i y j b ij
(19)
a ij = a ii a jj ( 1 - K ij )
(20)
b ij = b i + b j 2 1 - l ij

Where lij and Kij are interaction parameters and given:

(21)
a i ^ = n a m n i T , P , n j i = 2 N j = 1 y i a ij
(22)
b i ^ = n b m n i TP n j i = 2 N j = 1 y i b ij

B: CVD mixing rule:

(23)
a = x i x j a ij ( b m b i ) M ij
(24)
b m = x i b i
(25)
a = ( a i a j
(26)
b ij = ( b i b j )

Where M ij is the interaction parameter and given:

(27)
i ^ = ( na m ) ( n i ) T , P , n j i = 2 j = 1 N x i a ij b m b ij M ij + b i b m - 1 i j x i x j a ij M ij b m b ij M ij - a m
(28)
b i ^ = ( n b m ) ( n i ) T , P , n j i = 2 j = 1 N x i b i

Where M ij is subjected to the following relations:

(29)
M ii = M jj = 1
(30)
M ij = M ji

4

4 Results and discussion

4.1

4.1 Experimental data

In this study, the solubility of rivaroxaban in SC-CO2 was investigated at various temperatures (308–338 K) and pressures (120–270 bar). Drug solubility in SC-CO2 is given in Table 3. All experiments were repeated triple to increase measurement accuracy. Using Spane-Wagner EoS (Span and Wagner, 1996), the density of C O 2 was estimated. Figs. 2 and 3 show the change in drug solubility with density and pressure, respectively. As shown in Fig. 2, the solubility of the drug increases with increasing density. Also, according to Fig. 3, the solubility of the drug increases with increasing pressure at a constant temperature, which is the result of increasing the density of C O 2 that increases the solubility of SC-CO2 at higher pressures (Khimeche et al., 2007; Foster et al., 1991; Sodeifian et al., 2021). Temperature has two different effects on solubility changes. There are graphs intersect at pressure of 225 bar. Before this pressure, the solubility decreases with increasing temperature and after this point, increases with increasing temperature. So at pressures less than 225 bar, with increasing temperature, solubility of drug decrease, and at pressures above 225 bar with increasing temperature the solubility increases. This type of solubility behavior is consistent with reports of other compounds by other researchers (Perrotin-Brunel et al., 2010; Sodeifian et al., 2018).

The solubility of rivaroxaban in SC-CO2 vs density of SC-CO2 at different temperatures.
Fig. 2
The solubility of rivaroxaban in SC-CO2 vs density of SC-CO2 at different temperatures.
The solubility of rivaroxaban in SC-CO2 vs pressure at different temperatures.
Fig. 3
The solubility of rivaroxaban in SC-CO2 vs pressure at different temperatures.

4.2

4.2 Semi-empirical models

As stated before, six semi-empirical models were used to correlate solubility data of rivaroxaban in SC-CO2. Average absolute relative deviation (AARD), adjusted correlation coefficient ( R adj ) and F-value by MATLAB function (program) were used for the capability evaluation and verification of these models (Table 4). The results of these models are given in Table 5.

Table 4 Statistical tools for the capability evaluation and verification of the used models.
Equation Parameters References
A A R D % = 100 N i = 1 N i y i c a l c - y i e x p y i e x p
  • N is the number of data base.

  • y i c a l c is calculated values of the solute solubility

  • y i e x p is experimental value of the solute solubility

(Esmaeili et al., 2019)
R a d j = R 2 - ( Q 1 - R 2 N - Q - 1 ) R 2 = 1 - S S E S S T
  • Q is the number of independent variables in each equation

  • R 2 is the correlation coefficient

  • S S E and S S T is the error sum and total sum of squares

(Sodeifian and Sajadian, 2019)
F - v a l u e = S S R Q S S E N - Q - 1 = M S R M S E
  • S S R is the regression sum of squares

  • M S R and M S E is the mean square regression and residual

(Ardestani et al., 2020)
Table 5 Result of semi-empirical models for rivaroxaban solubility in SC-CO2.
Model Chrastil Bartle Jafari Nejad Jouyban Sodeifian Belghait
parameter
a 0 9.508 15.912 −24.973 −13.986 −20.982 −168.862
a 1 −5246.281 −7612.323 0.002 −1.104 1.137 −0.265
a 2 −49.845 0.015 1.583E-6 1.182E-5 1.106 3.969e-5
a 3 1.812 3.232E-5 7.596e-4 6.705e-4
a 4 0.701 −0.017 0.068
a 5 2.962 −389.092 −8.536e-5
a 6 0.357
a 7 50314.666
AARD 25.201 24.828 13.788 12.463 14.297 19.188
Radj 0.959 0.960 0.921 0.969 0.923 0.936
F-value 90.285 91.603 44.410 90.54 45.266

As shown in Table 5, considering the values of AARD for each model (Chrasil (AARD = 25.201); Bartle (AARD = 24.828), Jafari Nejad (AARD = 13.788), Jouyban (AARD = 12.463), Sodeifian (AARD = 14.297), and Belghait (AARD = 19.188)), Jouyban model has the lowest value of AARD and Chrastil model has the highest value of AARD. Therefore, the Jouyban model shows a better correlation with solubility data. Fig. 4 demonstrates the solubility correlation obtained from six models.

Comparison of experimental (point) and calculated (line) values of rivaroxaban solubility.
Fig. 4
Comparison of experimental (point) and calculated (line) values of rivaroxaban solubility.

In terms of the number of adjustable parameters, the rivaroxaban solubility could be poorly correlated using models with three adjustable parameters, as compared to those with four, six and seven parameters, as shown in Table 5. In addition to AARD, Radj was calculated to compare models of different numbers of independent variables. A model with more adjustable parameters can be considered superior to other models with less parameters that, in addition to AARD decreasement, Radj also increases. This can be seen for Jouyban model.

Elaborating on the results presented in this research, it can be reasoned that the energy term (coefficient of the temperature term (a1 = ΔH/R)) in the Chrastil and Bartle et al. models eased the determination of vaporization heat (ΔHvap), total reaction heat (ΔHtotal), and solvation heat (ΔHsol) of the considered drug-CO2 systems based on the regressed energy parameters. Indeed, the solvation heat (ΔHsol) was defined as the difference between ΔHvap (Bartle et al. model) and ΔHtotal (Chrastil’s model). According to Table 6, wherein enthalpy of dissolution of rivaroxaban in SC-CO2 is presented, the corresponding values of ΔHtotal and Δ Hvap to the Chrastil and Bartle et al. models were 43.617 and 63.288 kJ.mol−1, respectively. Given the endothermic and exothermic natures of vaporization and solvation processes, respectively, the value of vaporization heat was observed to be larger than the total heat. Representing the difference between ΔHvap and ΔHtotal, the value of ΔHsol was evaluated as −19.671 kJ.mol−1.

Table 6 Approximated Δ H t o t a l , Δ H v a p and Δ H s o l for rivaroxaban.
Compound Δ H t o t a l ( K J m o l ) a Δ H v a p ( K J m o l ) b Δ H s o l ( K J m o l ) c
Rivaroxaban 43.617 63.288 −19.671

aObtained from the Chrastil model.

bObtained from the Bartle model.

c Δ H sol = Δ H total - Δ H vap .

4.3

4.3 Solubility correlation by EoS based model

Two mixing rules (vdW2 and CVD) were used to estimate the parameters of PR and SRK EoSs. The physicochemical properties of rivaroxaban were calculated with Ambrose, Joback, Marrero and Pardillo, Grigoras, Fedors and Lee-Kesler methods and shown in Table 7.

Table 7 The physicochemical properties of rivaroxaban.
Component T b ( K ) T c ( K ) P c ( b a r ) ω V s ( c m 3 m o l ) T ( K )
308 318 328 338
P s u b ( P a ) f
Rivaroxaban 586.298 a 751.094 b 22.837 c 1.111 d 296.5 e 0.0047 0.0169 0.0557 0.1685
C O 2 304.2 73.8 0.228

aEstimated by Marrero and Pardillo (Poling et al., 2001).

bEstimated by Joback (Poling et al., 2001).

cEstimated by Grigoras method (Poling et al., 2001).

dEstimated by Ambrose-Walton (Poling et al., 2001).

eEstimated by Fedors (Fedors, 1974).

fEstimated by lee-Kesler (Lee and Kesler, 1975).

Table 8 shows the results two EoSs (PR and SRK) with two mixing rules (vdW2 and CVD). Parameters k ij and l ij were calculated for vdW2 mixing rule and parameters M ij were calculated for CVD mixing rule. In addition to the above, the value of the 123 parameters for each of the EoS with two mixing rules are shown in Table 8. As shown in Table 8, the lowest AARD are for PR EoS + vdW2 (AARD%=5.104, T = 308 K), SRK EoS + vdW2 (AARD%=5.314, T = 308 K), PR EoS + CVD (AARD%=7.364, T = 308 K), and SRK EoS + CVD (AARD%=9.729, T = 308 K) were obtained, respectively. In all four cases, the value of AARD generally decreased by decreasing the temperature.

Table 8 Solubility correlation results of rivaroxaban for two different EoSs.
Model Temperature Mij kij lij Radj F-value AARD%
PR EoS + vdW2 308 0.441 0.343 0.986 90.18 5.104
318 0.504 0.443 0.989 1.141 8.012
328 0.571 0.557 0.974 47.713 15.53
338 0.650 0.701 0.966 35.528 14.732
SRK EoS + vdW2 308 0.466 0.369 0.986 88.662 5.314
318 0.522 0.462 0.989 110.17 8.091
328 0.527 0.410 0.755 4.317 15.117
338 0.701 0.787 0.873 9.033 20.634
PR EoS + CVD 308 1.030 0.926 16.048 7.364
318 1.092 0.8 5.447 15.826
328 1.183 0.37 1.396 22.286
338 1.339 0.842 −0.038 34.938
SRK EoS + CVD 308 1.040 0.825 6.308 9.729
318 1.100 0.687 3.231 17.577
328 1.207 0.575 0.378 31.97
338 1.982 1.649 −0.828 71.861

Parameters k ij and l ij for vdW2 are calculated as follows:

(31)
k ij = A T + B
(32)
l ij = C T + D

The values A, B, C, and D for k ij and l ij were calculated with data in Table 8 and linear regression fitting results (Figs. 5 and 6 for PR and SRK, respectively) are shown in Table 9.

The interaction parameters for the values of k ij and l ij in PR-vdW2.
Fig. 5
The interaction parameters for the values of k ij and l ij in PR-vdW2.
The interaction parameters for the values of k ij and l ij in SRK-vdW2.
Fig. 6
The interaction parameters for the values of k ij and l ij in SRK-vdW2.
Table 9 Adjustable parameters of rivaroxaban.
Model Interaction parameters
K i j l i j
A B C D
PR-vdW2 0.0069 −1.7001 0.0119 −3.3262
SRK-vdW2 0.0071 −1.7393 0.012 −3.3755

Fig. 7 shows plots of experimental data (point) and calculated values (line) for rivaroxaban drug with PR-vdW2 (a & b), SRK-vdW2 (c & d), PR-CVD (e & f), and SRK-CVD (g & h) at various temperatures relative to density and pressure, respectively. The correlation between experimental data and calculated values in different modes is well shown.

Comparison of experimental (point) and calculated (line) solubility of rivaroxaban at various temperatures. (a) & (b): PR-vdW2, (c) & (d): SRK-vdW2, (e) & (f): PR-CVD, and (g) & (h): SRK-CVD.
Fig. 7
Comparison of experimental (point) and calculated (line) solubility of rivaroxaban at various temperatures. (a) & (b): PR-vdW2, (c) & (d): SRK-vdW2, (e) & (f): PR-CVD, and (g) & (h): SRK-CVD.
Comparison of experimental (point) and calculated (line) solubility of rivaroxaban at various temperatures. (a) & (b): PR-vdW2, (c) & (d): SRK-vdW2, (e) & (f): PR-CVD, and (g) & (h): SRK-CVD.
Fig. 7
Comparison of experimental (point) and calculated (line) solubility of rivaroxaban at various temperatures. (a) & (b): PR-vdW2, (c) & (d): SRK-vdW2, (e) & (f): PR-CVD, and (g) & (h): SRK-CVD.

5

5 Conclusion

The solubility of rivaroxaban was determined at various temperatures and pressures in SC- C O 2 for the first time. The range of solubility of rivaroxaban was between 0.0104 and 0.2062 ( × 10 - 4 ) (in mole fraction). Maximum solubility of rivaroxaban at T = 338 K and P = 270 bar was 0.2062 × 10 - 4 (in mole fraction). Six semi-empirical models were used to examine the correlation of experimental data. Among the used models, the best performance was obtained for Jouyban model (AARD%=12.4635), and the lowest was for Chrastil model (AARD%=25.201). Then, two EoSs (PR and SRK) with vdW2 and CVD mixing rules were used to examine the correlation of solubility data. The best performance was obtained with PR + vdW2 (AARD%=5.104 at 308 K). In both cases of semi-empirical models and EoSs, the values of R adj , F-value, and AARD% were reported. Also, the amount of vaporization and total enthalpies were calculated by Bartle and Chrastil models. Using solubility data, important information can be obtained about the fabrication of micro/nano-scale particles using SCF to be used in the medical and pharmaceutical industries.

Acknowledgement

The authors would like to thank of the deputy of research, University of Kashan, under Grant Number (#Pajoohaneh 1401/19) for supporting this valuable project. Also, researchers would like to thank Parsian drug Company (Tehran), especially, Drs. Shojaei and Yosoufian, for providing the required drug API.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. , , . Supercritical fluid extraction with carbon dioxide and ethylene. Fluid Phase Equilib.. 1983;14:147-156.
    [CrossRef] [Google Scholar]
  2. , , , . Experimental measurement and thermodynamic modeling of capecitabine (an anticancer drug) solubility in supercritical carbon dioxide in a ternary system: effect of different cosolvents. J. Chem. Eng. Data. 2020;65:4762-4779.
    [CrossRef] [Google Scholar]
  3. , , , . Formation of edible oil-loaded beeswax microparticles using PGSS – Particles from Gas-Saturated Solutions. J. Supercritical Fluids. 2021;169:105106
    [CrossRef] [Google Scholar]
  4. , , , . Representing solute solubility in supercritical carbon dioxide: a novel empirical model. Chem. Eng. Res. Des.. 2015;93:355-365.
    [CrossRef] [Google Scholar]
  5. , . Supercritical carbon dioxide as a green alternative to achieve drug complexation with cyclodextrins. Pharmaceuticals. 2021;14:562.
    [CrossRef] [Google Scholar]
  6. , , , , . Solubilities of solids and liquids of low volatility in supercritical carbon dioxide. J. Phys. Chem. Ref. Data. 1991;20:713-756.
    [CrossRef] [Google Scholar]
  7. , , , , . Semi-empirical correlation of solid solute solubility in supercritical carbon dioxide: comparative study and proposition of a novel density-based model. C. R. Chim.. 2018;21:494-513.
    [CrossRef] [Google Scholar]
  8. , , , . A combined model for the solubility of different compounds in supercritical carbon dioxide. Chem. Eng. Res. Des.. 2015;104:416-428.
    [CrossRef] [Google Scholar]
  9. , , , . A five-parameter empirical model for correlating the solubility of solid compounds in supercritical carbon dioxide. Fluid Phase Equilib.. 2016;411:74-80.
    [CrossRef] [Google Scholar]
  10. , . Supercritical process technology related to energy and future directions – An introduction. J. Supercritical Fluids. 2015;96:11-20.
    [CrossRef] [Google Scholar]
  11. , , , , . Continuous nanonization of lonidamine by modified-rapid expansion of supercritical solution process. J. Supercritical Fluids. 2018;133:486-493.
    [CrossRef] [Google Scholar]
  12. , , , , . Experimental investigation for the solubility and micronization of pyridin-4-amine in supercritical carbon dioxide. J. CO2 Util.. 2017;18:173-180.
    [CrossRef] [Google Scholar]
  13. , , , . Measurement of solid solubility of warfarin in supercritical carbon dioxide and recrystallization study using supercritical antisolvent process. Adv. Powder Technol.. 2018;29:479-487.
    [CrossRef] [Google Scholar]
  14. , , , . Semi-empirical models and a cubic equation of state for correlation of solids solubility in scCO2: dyes and calix[4]arenes as illustrative examples. Fluid Phase Equilib.. 2016;426:37-46.
    [CrossRef] [Google Scholar]
  15. , , , . Extended-duration rivaroxaban thromboprophylaxis in acutely ill medical patients: MAGELLAN study protocol. J. Thromb. Thrombolysis. 2011;31:407-416.
    [CrossRef] [Google Scholar]
  16. , , , . Cubic equation-of-state correlation of the solubility of some anti-inflammatory drugs in supercritical carbon dioxide. Fluid Phase Equilib.. 2006;239:188-199.
    [CrossRef] [Google Scholar]
  17. , , , . Assessment of the impact of rivaroxaban on coagulation assays: Laboratory recommendations for the monitoring of rivaroxaban and review of the literature. Thromb. Res.. 2012;130:956-966.
    [CrossRef] [Google Scholar]
  18. , , , , . A data-driven model for predicting the effect of temperature on oil-water relative permeability. Fuel. 2019;236:264-277.
    [CrossRef] [Google Scholar]
  19. , , , . Calculation of the solid solubilities in supercritical carbon dioxide using a new Gex mixing rule. J. Supercritical Fluids. 2009;51:148-158.
    [CrossRef] [Google Scholar]
  20. , , , . A three-parameter cubic equation of state for prediction of thermodynamic properties of fluids. J. Chem. Thermodyn.. 2008;40:84-95.
    [CrossRef] [Google Scholar]
  21. , . A method for estimating both the solubility parameters and molar volumes of liquids. Polym. Eng. Sci.. 1974;14:147-154.
    [CrossRef] [Google Scholar]
  22. , , , . Significance of the crossover pressure in solid-supercritical fluid phase equilibria. Ind. Eng. Chem. Res.. 1991;30:1955-1964.
    [CrossRef] [Google Scholar]
  23. , , . Solubilities of solids in supercritical fluids using dimensionally consistent modified solvate complex models. Fluid Phase Equilib.. 2009;283:97-101.
    [CrossRef] [Google Scholar]
  24. , , , . Measurement and thermodynamic modeling of solubility of Tamsulosin drug (anti cancer and anti-prostatic tumor activity) in supercritical carbon dioxide. J. Supercritical Fluids. 2020;163:104875
    [CrossRef] [Google Scholar]
  25. , , , . Correlation of solute solubility in supercritical carbon dioxide using a new empirical equation. Chem. Eng. Res. Des.. 2014;92:2734-2739.
    [CrossRef] [Google Scholar]
  26. , , , . Correlation of solid solubilities of pharmaceutical compounds in supercritical carbon dioxide with solution model approach. J. Taiwan Inst. Chem. Eng.. 2013;44:349-358.
    [CrossRef] [Google Scholar]
  27. , , , , . Prediction of solute solubility in supercritical carbon dioxide: a novel semi-empirical model. Chem. Eng. Res. Des.. 2010;88:893-898.
    [CrossRef] [Google Scholar]
  28. , , , , . Gas-Antisolvent (GAS) crystallization of aspirin using supercritical carbon dioxide: experimental study and characterization. Ind. Eng. Chem. Res.. 2015;54:3685-3696.
    [CrossRef] [Google Scholar]
  29. , , . Relationship between the binary interaction parameters (kij) of the Peng-Robinson and those of the Soave–Redlich–Kwong equations of state: Application to the definition of the PR2SRK model. Fluid Phase Equilib.. 2010;295:26-37.
    [CrossRef] [Google Scholar]
  30. , , . Estimation of pure-component properties from group-contributions. Chem. Eng. Commun.. 1987;57:233-243.
    [CrossRef] [Google Scholar]
  31. , , , . Mathematical representation of solute solubility in supercritical carbon dioxide using empirical expressions. J. Supercritical Fluids. 2002;24:19-35.
    [CrossRef] [Google Scholar]
  32. Jouyban, A., Khoubnasabjafari, M., Acree, W.E., 2005. Mathematical representation of solute solubility in binary mixture of super- critical fluids by the Jouyban-Acree model. 3.
  33. , , , . Modeling the entrainer effects on solubility of solutes in supercritical carbon dioxide. Chem. Pharm. Bull.. 2005;53:290-295.
    [CrossRef] [Google Scholar]
  34. Khimeche, K., Alessi, P., Kikic, I., Dahmani, A., 2007. Solubility of diamines in supercritical carbon dioxide Experimental determination and correlation. 10.
  35. , . Supercritical fluids and polymers – The year in review – 2014. J. Supercritical Fluids. 2016;110:126-153.
    [CrossRef] [Google Scholar]
  36. , , , . Effect of the Z-isomer content on nanoparticle production of lycopene using solution-enhanced dispersion by supercritical fluids (SEDS) J. Supercritical Fluids. 2018;133:291-296.
    [CrossRef] [Google Scholar]
  37. , , . A generalized thermodynamic correlation based on three-parameter corresponding states. AIChE J.. 1975;21:510-527.
    [CrossRef] [Google Scholar]
  38. , , , . Microencapsulation of red palm oil as an oil-in-water emulsion with supercritical carbon dioxide solution-enhanced dispersion. J. Food Eng.. 2018;222:100-109.
    [CrossRef] [Google Scholar]
  39. , , , . Storage stability and degradation kinetics of bioactive compounds in red palm oil microcapsules produced with solution-enhanced dispersion by supercritical carbon dioxide: a comparison with the spray-drying method. Food Chem.. 2020;304:125427
    [CrossRef] [Google Scholar]
  40. , , , . Rivaroxaban in patients with a recent acute coronary syndrome. N. Engl. J. Med.. 2012;366:9-19.
    [CrossRef] [Google Scholar]
  41. , , , . Micronization of vanillin by rapid expansion of supercritical solutions process. J. CO2 Util.. 2017;21:169-176.
    [CrossRef] [Google Scholar]
  42. , , , . Rivaroxaban: population pharmacokinetic analyses in patients treated for acute deep-vein thrombosis and exposure simulations in patients with atrial fibrillation treated for stroke prevention. Clin. Pharmacokinet.. 2011;50:675-686.
    [CrossRef] [Google Scholar]
  43. , . Modified Wilson’s model for correlating solubilities in supercritical fluids of some polycyclic aromatic solutes. Polycyclic Aromat. Compd.. 2018;38:244-256.
    [CrossRef] [Google Scholar]
  44. , , , , . Solubility of 9-tetrahydrocannabinol in supercritical carbon dioxide. Experiments Modeling. 2010;5
    [Google Scholar]
  45. , , , . The properties of gases and liquids (5th ed.). New York: McGraw-Hill; .
  46. , , , . An investigation into Sunitinib malate nanoparticle production by US- RESOLV method: effect of type of polymer on dissolution rate and particle size distribution. J. Supercritical Fluids. 2021;170:105163
    [CrossRef] [Google Scholar]
  47. , , . Dimensionless empirical model to correlate pharmaceutical compound solubility in supercritical carbon dioxide. Chem. Eng. Technol.. 2019;42:2621-2630.
    [CrossRef] [Google Scholar]
  48. , , , , . Comparison and modelling of rutin solubility in supercritical carbon dioxide and subcritical 1,1,1,2-tetrafluoroethane. J. CO2 Util.. 2017;21:1-8.
    [CrossRef] [Google Scholar]
  49. Shamsipur M, Karami AR, Yamini Y, Sharghi H (2004) Solubilities of some 1-hydroxy-9,10-anthraquinone derivatives in supercritical carbon dioxide. 7.
  50. , , , . Modeling of ionic liquid+polar solvent mixture molar volumes using a generalized volume translation on the Peng-Robinson equation of state. Fluid Phase Equilib.. 2015;395:51-57.
    [CrossRef] [Google Scholar]
  51. , , , . Solubility measurement method and mathematical modeling in supercritical fluids. EJ. 2013;17:67-78.
    [CrossRef] [Google Scholar]
  52. , , , . Novel density-based model for the correlation of solid drugs solubility in supercritical carbon dioxide. C. R. Chim.. 2017;20:559-572.
    [CrossRef] [Google Scholar]
  53. , , , . Determination of solubility of Aprepitant (an antiemetic drug for chemotherapy) in supercritical carbon dioxide: empirical and thermodynamic models. J. Supercritical Fluids. 2017;128:102-111.
    [CrossRef] [Google Scholar]
  54. , , , , . Solubility measurement of an antihistamine drug (Loratadine) in supercritical carbon dioxide: assessment of qCPA and PCP-SAFT equations of state. Fluid Phase Equilib.. 2018;472:147-159.
    [CrossRef] [Google Scholar]
  55. , , , . Solubility measurement of a pigment (Phthalocyanine green) in supercritical carbon dioxide: Experimental correlations and thermodynamic modeling. Fluid Phase Equilib.. 2019;494:61-73.
    [CrossRef] [Google Scholar]
  56. , , , , . Experimental data and thermodynamic modeling of solubility of Sorafenib tosylate, as an anti-cancer drug, in supercritical carbon dioxide: Evaluation of Wong-Sandler mixing rule. J. Chem. Thermodyn.. 2020;142:105998
    [CrossRef] [Google Scholar]
  57. , , . Experimental measurement of solubilities of sertraline hydrochloride in supercriticalcarbon dioxide with/without menthol: data correlation. J. Supercritical Fluids. 2019;149:79-87.
    [CrossRef] [Google Scholar]
  58. , , , . Evaluation of the response surface and hybrid artificial neural network-genetic algorithm methodologies to determine extraction yield of Ferulago angulata through supercritical fluid. J. Taiwan Inst. Chem. Eng.. 2016;60:165-173.
    [CrossRef] [Google Scholar]
  59. , , , . Mathematical modelling for extraction of oil from Dracocephalum kotschyi seeds in supercritical carbon dioxide. Nat. Prod. Res.. 2018;32:795-803.
    [CrossRef] [Google Scholar]
  60. , , , . Solubility measurement of a chemotherapeutic agent (Imatinib mesylate) in supercritical carbon dioxide: assessment of new empirical model. J. Supercritical Fluids. 2019;146:89-99.
    [CrossRef] [Google Scholar]
  61. , , , . Prediction of solubility of sodium valproate in supercritical carbon dioxide: experimental study and thermodynamic modeling. J. Chem. Eng. Data. 2020;65:1747-1760.
    [CrossRef] [Google Scholar]
  62. , , , . Prediction of solubility of sunitinib malate (an anti-cancer drug) in supercritical carbon dioxide (SC–CO2): experimental correlations and thermodynamic modeling. J. Mol. Liq.. 2020;297:111740
    [CrossRef] [Google Scholar]
  63. , , , . Experimental measurement and thermodynamic modeling of Lansoprazole solubility in supercritical carbon dioxide: application of SAFT-VR EoS. Fluid Phase Equilib.. 2020;507:112422
    [CrossRef] [Google Scholar]
  64. , , , , . Measuring and modeling the solubility of an antihypertensive drug (losartan potassium, Cozaar) in supercritical carbon dioxide. J. Mol. Liq.. 2021;331:115745
    [CrossRef] [Google Scholar]
  65. , , , , . The solubility of Sulfabenzamide (an antibacterial drug) in supercritical carbon dioxide: evaluation of a new thermodynamic model. J. Mol. Liq.. 2021;335:116446
    [CrossRef] [Google Scholar]
  66. , , , , . Solubility of Amlodipine Besylate (Calcium Channel Blocker Drug) in supercritical carbon dioxide: measurement and correlations. J. Chem. Eng. Data. 2021;66:1119-1131.
    [CrossRef] [Google Scholar]
  67. , , . Investigation of essential oil extraction and antioxidant activity of Echinophora platyloba DC. using supercritical carbon dioxide. J. Supercritical Fluids. 2017;121:52-62.
    [CrossRef] [Google Scholar]
  68. , , , , . Experimental data and thermodynamic modeling of solubility of Azathioprine, as an immunosuppressive and anti-cancer drug, in supercritical carbon dioxide. J. Mol. Liq.. 2020;299:112179
    [CrossRef] [Google Scholar]
  69. , , . A new equation of state for carbon dioxide covering the fluid region from the triple-point temperature to 1100 K at pressures up to 800 MPa. J. Phys. Chem. Ref. Data. 1996;25:1509-1596.
    [CrossRef] [Google Scholar]
  70. , , , . Evaluation of density-based models for the solubility of solids in supercritical carbon dioxide and formulation of a new model. Chem. Eng. Sci.. 2008;63:4292-4301.
    [CrossRef] [Google Scholar]
  71. , , , . A new model for the solubility of dye compounds in supercritical carbon dioxide. Thermochim Acta. 2013;561:91-97.
    [CrossRef] [Google Scholar]
  72. , , , . Supercritical fluids for pharmaceutical particle engineering: methods, basic fundamentals and modelling. Chem. Eng. Process. Process Intensif.. 2012;60:9-25.
    [CrossRef] [Google Scholar]
  73. , , , . Microencapsulation of drug with enteric polymer Eudragit L100 for controlled release using the particles from gas saturated solutions (PGSS) process. J. Supercritical Fluids. 2021;167:105044
    [CrossRef] [Google Scholar]
  74. , , , . Thermodynamic modeling of solubilities of various solid compounds in supercritical carbon dioxide: effects of equations of state and mixing rules. J. Supercritical Fluids. 2011;55:861-875.
    [CrossRef] [Google Scholar]
  75. , , , . Applications of cubic equations of state for determination of the solubilities of industrial solid compounds in supercritical carbon dioxide: a comparative study. Chem. Eng. Sci.. 2012;71:283-299.
    [CrossRef] [Google Scholar]

Appendix A

Supplementary material

Supplementary data to this article can be found online at https://doi.org/10.1016/j.arabjc.2022.104421.

Appendix A

Supplementary material

The following are the Supplementary data to this article:

Supplementary data 1

Supplementary data 1

Show Sections