Controlling the product crystal size distribution by ... - ACS Publications

1. INTRODUCTION. 26. In the manufacturing of Active Pharmaceutical Ingredients (API), crystallization is. 27 of vital importance and is often the crit...
10 downloads 8 Views 3MB Size
Subscriber access provided by UNIV OF DURHAM

Article

Controlling the product crystal size distribution by strategic application of ultrasonication Kiran A. Ramisetty, and Åke C. Rasmuson Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.7b01619 • Publication Date (Web): 20 Feb 2018 Downloaded from http://pubs.acs.org on February 20, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Crystal Growth & Design is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

1

Controlling the product crystal size distribution by strategic application of ultrasonication

2

Kiran A. Ramisetty, Åke C. Rasmuson*

3

Department of Chemical and Environmental Science, Synthesis and Solid State Pharmaceutical

4

Centre (SSPC), Bernal Institute, University of Limerick, Limerick, Ireland.

5

*Email: [email protected]

6

ABSTRACT:

7

In this work, different strategies of ultrasonication (continuous, single-pulse and multiple-pulse)

8

are compared for control of the product crystal size distribution of three model API compounds:

9

piracetam, paracetamol, and ibuprofen. Experiments have been performed in 0.5L and 3L scale

10

continuously recorded by FTIR and FBRM. Irrespective of the sonication operating mode,

11

sonication in general produced smaller sized crystals with a more narrow size distribution than a

12

normal cooling crystallization process. A multiple pulse sonication mode, in particular, was

13

capable of delivering more narrow size distributions. Sonication power per unit mass of solution

14

does not appear to be a relevant scaling-up parameter.

15 16 17 18 19 20 21 22 23 24 25

1 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

26

1. INTRODUCTION

27

In the manufacturing of Active Pharmaceutical Ingredients (API), crystallization is

28

of vital importance and is often the critical process technology to ensure the product purity,

29

and to deliver the required crystal size, shape, surface area and bulk density; of importance

30

for downstream processing and end-use. Traditionally, often the solid product from the

31

API manufacturing is exposed to a milling step, which adds another unit operation to the

32

manufacturing is energy consuming and may lead to undesirable or even unacceptable

33

alternations of essential product characteristics like changes in crystal structure or

34

dehydration. A desirable alternative would be directly in the crystallization step to produce

35

the right size distribution that meets the requirements of the drug formulation. Thus,

36

processes that can deliver a narrow size distribution around a particular mean size below

37

100 µm is of significant interest.

38

Careful monitoring of the solution concentration using an in situ ATR-FTIR probe

39

is becoming popular and gaining a considerable attention as it allows to precisely control

40

crystallization kinetics and thus the final product quality.1,2,3 During the API

41

manufacturing process, size and shape engineered crystals are usually obtained either in

42

batch or continuous mode via concentration feedback control strategy, anti-solvent

43

addition, optimal cooling with or without the addition of seed crystals, etc.4 Under this

44

context, feedback control strategy of direct nucleation control (DNC) seems to be a

45

promising technique to control final crystal size.5,6 The DNC approach relies on a protocol

46

to directly control the apparent onset of nucleation or in situ fines removal by switching

47

between heating and cooling cycles or solvent or anti-solvent addition to the nucleated

48

supersaturation solution until a desired number of FBRM counts is reached. Additionally,

49

DNC provides a model-free approach to design and does not require concentration

50

measurement, or information on nucleation or growth kinetics to design an operating

51

curve. Despite the advantages of these controlling strategies, they usually produce larger

52

crystals with narrow crystal size distribution (CSD) depending on the method applied and

53

thus require the additional step of milling (to reduce the crystal size) to meet size

54

requirements for pharmaceutical processing in various dosage forms. It is worth to mention

55

here, that dry milling is unsuitable for potent compounds for reasons like human exposure, 2 ACS Paragon Plus Environment

Page 2 of 37

Page 3 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

56

undesirable polymorphic changes, and product loss. Although the wet milling can mitigate

57

some of these issues, the high shear force and breakage due to attrition usually produce

58

bimodal CSD.

59

Technologies such as spray drying, supercritical anti-solvent, sono-crystallization,

60

sono-milling can surpass the limitations discussed above and provide an elegant way to

61

produce

62

sonocrystallization pose some advantages over the other techniques. Recent advances in

63

the design of ultrasonic probes that can operate either in-situ or transfer ultrasound

64

externally (or indirect sonication) via a metal column or the steel wall of the crystallizer

65

which is attached to an ultrasonicator probe, makes this technique versatile and easy to

66

scale-up at least at the pilot plant scale level. Current technologies and methods allow the

67

design of a sonicated crystallizer that handles large batch volume (required by industries)

68

that are comparable to any conventional crystallizers.7 Though indirect sonication is

69

expensive, it enjoys the advantage of no metal scavenging, a much-needed design property

70

for API crystallization as it avoids the incorporation of impurities in the final product. In

71

sonocrystallization, nucleation occurs within the supersaturated solution, induced by

72

transient cavitation when an ultrasound of sufficient amplitude is applied.8 The

73

mechanisms involved in sonocrystallization have been discussed elsewhere.9,10

smaller

sized

crystals

with

targeted

CSD

directly.

Amongst

these,

74

Generally, nucleation can be enhanced by increasing the supersaturation, addition

75

of appropriate anti-solvents or by adjusting the cooling rate. Alternatively, nucleation can

76

be created/enhanced by introducing an external field such as ultrasound. Conceptually,

77

ultrasonication increases the nucleation rate as it effectively creates cavitation, which

78

enhances the creation of a large number of the nucleus via heterogeneous primary

79

nucleation. Earlier studies confirm introducing ultrasound to a conventional crystallizer

80

will drastically modify the nucleation and the growth kinetics.11 Sonocrystallization is

81

considered a proven technique to produce smaller crystals; however, the strategies required

82

to produce uniform-sized crystals are not clear. As discussed above, product size and CSD

83

can be controlled simultaneously by identifying the right control strategy for

84

sonocrystallization.

3 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 37

85

The purpose of this work is to propose and test different ultrasonication strategies

86

to precisely control the product CSD concerning mean size and distribution width. The

87

strategies directly rely on ‘when and how’ the ultrasound is applied while operating the

88

batch cooling crystallizer. To execute this approach, we systematically modified the mode

89

of ultrasonication and determined the product size distribution, and use various PAT

90

technologies to trace the behaviour of the process.

91

piracetam, paracetamol, and ibuprofen were used in this work, and the objectives were to

92

produce a narrow size distribution in the lower size range (20 – 50 µm), potentially

93

allowing for downstream milling to be circumvented.

94

2. MATERIALS AND METHODS

95

2.1. Materials. Piracetam (Oxo Industries) recrystallized before use; Isopropanol (99.5%),

96

Acetaminophen (98.0-102.0%), Ibuprofen, and Acetonitrile (99.8%), purchased from

97

Sigma Aldrich and used as obtained

98

2.2 Equipment

99

Figure.1 shows the schematic representation of the experimental setup and image analysis

100

procedures used to determine the CSD. The experiments were performed in 2 different

101

glass-jacketed reactors from Lenz Laborglas (GmbH & Co.KG) with working volumes of

102

0.5L (CR1) and 3L (CR2). To prevent vortexing and mixing dead zones, and to ensure

103

complete mixing, two beavertail (Cowie) baffles with a baffle clearance of 0.5 cm placed

104

inside the crystallizers. A 45 degree pitched blade impeller with 4 blades placed

105

diametrically at the center of CR1 and CR2, respectively, with a clearance of 16.5mm for

106

CR1 and 26mm for CR2 from the bottom of the reactor. Precision controlled overhead

107

stirrer (IKA Eurostar 100) used to adjust the impeller speed. All the experiments performed

108

at 500 rpm in CR1 and 464 rpm in CR2. The agitation speed was chosen carefully to

109

ensure and maintain the same agitation power per unit mass (ɛ = 0.211W/kg). The

110

temperature of the solution mixture measured with the help of a PT 100 thermocouple,

111

controlled by LAUDA E300 circulator and Wintherm Software. An ultrasonic processor

112

(Qsonica Q500) with a maximum power of 500W and with operating a fixed frequency of

4 ACS Paragon Plus Environment

Three model API compounds,

Page 5 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

113

20 kHz used (VWR International Ltd) to study the effect of ultrasound on the process of

114

crystallisation.

115 116

Specifications Volume Dimensions Impeller Diameter Agitation speed Power per unit mass (ε) Number of Baffles Impeller Clearance from bottom of the reactor 117

Crystallizer-1 (CR1) 0.5(L) Diameter = 100 mm, Height = 110 mm 40mm 500rpm ɛ = 0.211W/kg 2 + (3 in-situ Probes)

Crystallizer-2 (CR2) 3(L) Diameter = 150 mm, Height = 230 mm 60mm 464rpm ɛ = 0.211W/kg 2 + (3 in-situ Probes)

16.5mm

26mm

Table 1, Crystallizer specifications Most of the experiments were performed with

118

2.3. Experimental Methodology.

119

Piracetam; the procedures of these experiments are described in detail here. Some

120

experiments were performed with Acetaminophen and Ibuprofen and the conditions are

121

described in conjunction with presenting the corresponding results in chapter 3.5. A

122

solution saturated at 40 oC of Form III Piracetam in isopropanol with a concentration of

123

32g/L was prepared and heated to 45 oC, and maintained for 30 min at 500-rpm agitation to

124

ensure complete dissolution.12 The solution was then cooled down to 40 oC, maintained for

125

30min and then cooled down to 0 oC at specified cooling rates. Three different ultrasonic

126

power levels: 25, 43 and 60 W were used, and the power delivered to the liquid per unit

127

volume (W/kg) was determined by the calorimetric method.

128

5 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

129 130

Figure 1. Schematic representation of experimental setup and CSD analysis

131

This work mainly focuses on the influence of sonication on the final crystal size

132

distribution. The ultrasound is applied in a cooling crystallizer in different modes: (i) CU –

133

continuous ultrasonication, (ii) SPU - single-pulse ultrasonication and (iii) MPU – multi-

134

pulse ultrasonication. In CU, ultrasound was applied throughout the crystallization time. In

135

SPU the ultronication is turned on after 2500 s of cooling at 0.5 ℃/min when the solution

136

is close to the metastable limit (MSZW as determined in the absence of sonication in CC

137

experiments). The influence of the time of ultrasonication on the crystal size distribution

138

was studied for three different US times: 2, 7 and 12 min at a constant sonication power of

139

60 W/kg. The influence of the US power on the CSD was studied at a constant US time of

140

2 minutes and with three different US power, 60, 103 and 144 W/kg. In MPU, ultrasound

141

was applied to the crystallizer at certain and specific time intervals (2min on and 5min off

142

per cycle), beginning from ∆Tmax to the end of the crystallization experiment (see section

143

3.3 for more details).

144

Introducing ultrasound into industrial crystallizers can create undesirable noise in

145

addition to other technical difficulties such as crystal contamination by metal scavenging,

146

attrition of the crystals affecting product quality and thermal fluctuations in the suspension

147

etc..13,14,15. Ultrasonication transfer energy to the solution leading to a rise in temperature

148

that tends to upset the intended cooling profile. This effect was observed in the 6 ACS Paragon Plus Environment

Page 6 of 37

Page 7 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

149

introductory CU sonocrystallization experiments. The heat generated by sonication in CU

150

mode created a temperature upset at least by +5oC from the actual set temperature profile at

151

a cooling rate of 1oC/min. In the sonocrystallization experiments performed with 0.1

152

o

153

of 0.5 oC/min was chosen as it provided a balance between sonication time and solution

154

temperature control (the temperature difference between the linear profile and the solution

155

temperature is < 0.3 oC). Additionally, in order to bypass the temperature upset created by

156

sonication, an equilibration time was provided at the beginning of each sonication

157

experiment (not shown), which allows the temperature to reach the set point or the

158

saturation temperature. At some stage and only for demonstration purpose, experimental

159

results obtained for different cooling rates and US power are briefly discussed. Only for

160

demonstration purpose, experimental results obtained for different cooling rates and US

161

power are briefly discussed.

162

2.4. Analytical Technique. The concentration of Piracetam was monitored online with the

163

Attenuated Total Reflectance Fourier Transform Infrared in-situ probe (ATR-FTIR, model

164

iC 10, DiComp probe with AgX halide fibre). The spectra with a resolution of 4 cm-1 were

165

collected in the spectral region of 600 - 1890 cm-1 and analysed using iC IR software v.4.3.

166

FTIR spectra showed a characteristic peak for Piracetam at 1684 cm-1 and for isopropanol

167

950 cm-1 (see Figure 2).16 The dissolved solute concentration was obtained using a relation

168

that correlates the ratio of peak area of piracetam to the peak area of the isopropanol. This

169

relation was obtained by dissolving 20 g of Piracetam in 500 ml of isopropanol at 0oC. This

170

solution was then heated up to 40oC at a very low heating rate of 0.1oC/min. It is assumed

171

that, this low heating rate provides the sufficient time for the piracetam to completely reach

172

the saturation at the studied temperature range. The ratio of the peak area of piracetam to

173

the peak area of isopropanol, PAp/PAipa recorded during this experiment was then plotted

174

against the solubility to obtain a relation that correlates the dissolved concentration as a

175

function of PAp/PAipa (see the polynomial expression in Figure.2 (B)). This relation was

176

later used in the cooling crystallization experiments in order to determine the concentration

177

from the PAp/PAipa (See in Figure.2 (C)). Additionally, IR absorbance (A.U) is sensitive to

178

temperature and in general, they linearly increase with a decrease in temperature. Thus,

179

during the cooling crystallization experiments, after converting the PAp/PAipa to

C/min the sonication time became prohibitedly long. Thus, an intermideate cooling rate

7 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

180

concentration, a correction was applied to the concentration in order to subtract the

181

temperature effect (see in Figure.2 (D)). This was done by subtracting the concentration by

182

a value m*T. Where T is the temperature and m is the slope that measures the change in

183

concentration with respect to temperature.

184 185

Figure 2. A) ATR- FTIR spectra (600-1900cm-1) for isopropanol and a saturated solution of

186

Piracetam in isopropanol. B) The correlation between solubility v’s peak area ratio. C)

187

Calculated solubility from peak area ratio. D) Temperature correction on concentration

188

profile in typical experiment.

189

Focused Beam Reflectance Measurement (FBRM, model S400A) manufactured by

190

Metler-Toledo was used to monitor the chord length distribution (CLD) and change in

191

count rate with the measurement time of 10 s using the software iCFBRM. The nucleation

192

or cloud point was determined with the help of FBRM and ATR-FTIR. The FBRM and

193

ATR-FTIR probes were placed at a certain distance (60mm) from the ultrasonic probe. In

194

Figure 3, the initial rise in FBRM counts when the ultrasonication starts is due to the

195

cavitation bubbles generated. However, this contribution to the number counts is negligible

196

compared to the following crystal nucleation contribution. For an accurate measurement of

197

the CSD, optical microscopy was (Olympus IX53) used. Crystals suspended in antisolvent 8 ACS Paragon Plus Environment

Page 8 of 37

Page 9 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

198

or saturated solution in a glass cuvette (Starna Scientific) were dispersed in an ultrasonic

199

bath for 3 s and then placed under a microscope. The glass cuvette image was divided into

200

regular frames based on the magnification of the microscopic lens to avoid repeated

201

measurement of the same crystals. Image analysis software called ImageJ was then used to

202

analyze the images to determine the crystal size distribution (CSD). The Ferret diameter

203

was used as the characteristic linear dimension of the crystal. The ImageJ macro program

204

was used to analyze a number of images by which more than 20,000 crystals from each

205

experiment were analyzed. MS-Excel and Minitab software were used to analyze the CSD

206

and is in this paper presented as number based percentage distribution.

Point of nucleation Cavity bubbles shown in FBRM counts

207 208

Figure 3. Observed cavity bubble size smaller than crystal size in FBRM count rate.

209

To calibrate the FBRM count rate in separate experiments, the solution was seeded

210

by a well-defined number of crystals of known size.17 Different size range of seed crystals

211

were carefully isolated by sieving using 125, 250, 300, 355 and 500 µm mesh sieve plates

212

using dry sieving followed by wet sieving. The volumetric shape factor (0.93) was

213

determined by counting and weighing the seed crystals, which was the same for all seed 9 ACS Paragon Plus Environment

Crystal Growth & Design

214

crystals regardless of their size. Using this shape factor (f) and the crystal density (ρc) the

215

number of crystals (N) added per unit volume of solution (V) is known from the equation. =

(

 )   

216

For example, for 250-300µm sieved crystals of mass 1.5 g and 3 g, respectively, the

217

corresponding number concentration becomes 7.58 × 107 and 1.52 × 108 (#/m3). In Figure

218

4, the relation between the actual number concentration and the FBRM count rate at the

219

CR1 conditions is shown. A power law model is fitted to the data and the exponent

220

becomes very close to unity. This correlation is used to convert the count rate during the

221

experiments to the actual number of concentrations. This correlation is only applicable to

222

CR1 conditions. 1.E+10 Actual Number of Crystals/m3

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 37

y = 567274x0.984 R² = 0.9858

1.E+9

1.E+8

1.E+7 10

100

223

Counts/s

1000

10000

224

Figure 4. Calibration of FBRM counts with an actual number of crystals per unit volume

225

for CR1.

226

The crystal structure of the product has been examined by PXRD. In all experiments only

227

the stable Form III with nice hexagonal block shape crystals is found.

228

3. RESULTS AND DISCUSSION.

10 ACS Paragon Plus Environment

Page 11 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

229

Figure 5 shows the product crystal size distribution and the temperature as a function of the

230

time profile during the cooling crystallization of piracetam at three different modes of

231

ultrasonication. For comparison purposes, we also show the results obtained from

232

conventional crystallization (CC-no ultrasound) experiments. All experiments started at 40

233

o

234

few experiments where the cooling rate is 0.1 or 1 oC/min. For the cases of SPU and MPU,

235

the arrows in Figure 5 show the time at which ultrasound is introduced to the system and

236

how the temperature of the solution increases due to the ultrasonication (on-cycle), and

237

rapidly decreases during the off-cycle. The rapid cooling is the result of the temperature

238

control system trying to bring the solution temperature back to the linear cooling profile.

239

The crystal size distribution shown in Figure 5 points to a major conclusion that switching

240

sonication mode can produce crystals with different size and size distribution. Below the

241

manuscript is structured to detail how the crystal size distribution is influenced by the

242

mode of the sonication, and how the outcome compares to a normal cooling process, and to

243

explain why certain experimental conditions produce more narrow crystals size

244

distributions, for instance SPU, and MPU can produce smaller crystals with a a more

245

narrow size distribution compared to the CU mode.

246 247

Figure 5. The temperature profiles at various conditions of ultrasound modes, 1. CC, 2.

C and the solution was cooled to 0 oC with a constant cooling rate of 0.5 oC/min, besides a

11 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

248

CU, 3. SPU, 4. MPU and their corresponding final CSD. (Arrow marks show the point

249

where the ultrasound applied for SPU and MPU)

250

3.1 Continuous Ultrasonication (CU):

251

Figure 6 shows the concentration vs temperature curve obtained from CU experiments

252

performed with different ultrasonication power at a constant cooling rate of 0.5 oC/min.

253

The results obtained from the conventional cooling crystallization experiments performed

254

at different cooling rates are shown for comparison. Quite clearly, the depletion in

255

concentration begins at an earlier temperature and the concentration rapidly decreases

256

towards the solubility line in the sonicated experiments. Theoretically this difference in

257

saturation temperature and the temperature at which nucleation occurs is defined as the

258

metastable zone width (MSZW). The MSZW obtained from Figure 6 is plotted against

259

product crystal mean size for different operating conditions in Figure 7. Figure 7 shows

260

that, ultrasonication significantly decreases the MSZW. With conventional cooling

261

crystallization performed at a cooling rate of 0.1, 0.5 and 1oC/min, the MSZW was 14.5,

262

19.4 and 23.4 oC, respectively (Table 4). This is as expected, much higher compared to the

263

MSZW of 8.1, 7.2 and 5.1 oC observed in the sonicated experiments performed with a

264

cooling rate of 0.5 oC/min and ultrasound power of 60, 103 and 144 kg/m3, respectively.

265

The decrease in MZW, reflects that nucleation in a supersaturated solution is stimulated by

266

the sonication, and increasingly so with increasing sonication power.18,19 The MSZW can

267

be altered by changing the ultrasonication power; MSZW of the sonicated crystallization

268

performed at 0.5 oC/min is significantly lower than that of a conventional cooling

269

crystallization even when performed at a higher (1 oC/min) or lower cooling rate (0.1

270

o

271

conditions (no ultrasound) with a very low cooling rate still produced a higher MSZW,

272

compared to the CU experiments. Increase in the creation of collapsing cavitational events

273

with increase in power of ultrasound would be the possible reason for an earlier nucleation

274

at increasing power of ultrasonication. Even with the low US frequency of 20kHz it is

275

possible to reduce the energy barrier for the nucleation with a high US power.

C/min). As evidenced from Figure 6, the cooling crystallization performed at silent

12 ACS Paragon Plus Environment

Page 12 of 37

Page 13 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

276 277

Figure 6. Concentration vs temperature in the case of CC and CU experiments.

CC

CU

278 279

Figure 7: MSZW and its corresponding mean crystal size at CU and CC

280

experiments: From table 1, exp no 1-6 in CC-mode and exp no 7-15 in CU-mode are

281

included in this graph for comprehensive relation between mean size and MSZW.

13 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

282 283

Figure 8. CSD for CC and CU experiments at various conditions and microscopic image of

284

crystals (A: 0.5oC/min for CC, B: 0.5oC/min, for CU 60W/kg, (scale bar is 500µm))

285

Figure 8 shows the influence of the cooling rate in a CC process and of the sonication

286

power in a CU process on the final crystal size distribution. The crystal size and the

287

distribution width (D90-D10/D50) (CSDW) obtained using the microscopic image analysis

288

techniques are given in Table 4. Increasing the cooling rate and sonication power decrease

289

the mean size and the CSDW, but with sonication, the product crystals are much smaller

290

and the distribution is much more narrow. The reported hypothesis is that solvent

291

evaporation at the hotspots created by cavity implosion, produces a massive number of

292

nucleation sites which results in a reducing MSZW and product crystal size20,21. It has

293

been reported that the cavitation bubbles act similarly to seed crystals or gas bubbles22.

294

Since the nucleation time distribution basically controls the product crystal size

295

distribution, an increase in nucleation rate due to ultrasonication leads to smaller product

296

crystals when using ultrasound.23 CU experiments performed at a constant cooling rate

297

show a moderate decrease in crystal size and a fairly significant decrease in CSDW with

298

increasing US power from 60W/kg to 144W/kg: for CSDW from 1.72 to 1.45 at 0.1 14 ACS Paragon Plus Environment

Page 14 of 37

Page 15 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

299

o

C/min, from 1.65 to 1.32 at 0.5 oC/min, and from 1.55 to 1.35 at 1 oC/min. CU

300

experiments performed at constant US power showed a decrease in mean crystal size and

301

CSDW at an increasing cooling rate, but the influence of cooling rate appears to decrease

302

as the sonication power increases. In CC experiments, the mean crystal size decreases from

303

133 to 70 µm as the cooling rate increases from 0.1 to 1oC/min, and the CSDW becomes

304

higher (>2) for the highest cooling rate (1oC/min). Continuous ultrasonication at 144 W/kg

305

and a cooling rate of 1 oC/min) produced the smallest sized crystals (~25 um) and a CSDW

306

of 1.35.

307

The reduced mean size is a consequence of ultrasonication generating a higher

308

number of crystals to share the total product mass as is illustrated in Figure 9. In the CC

309

process, a higher cooling rate leads to an increased number of crystals, and ultrasonication

310

substantially increases the number. The graphs reveal that the number generation includes

311

two different regions: a very fast primary nucleation and a following slower secondary

312

nucleation, in accordance with previous reports.24,25 From straight-line approximations of

313

the slopes, we have estimated the corresponding nucleation rates and they are collected in

314

Table 2. Obviously, the primary nucleation rates are much higher in the CU process as

315

compared to the CC process even though the supersaturation is much lower (Figure 6).

316

Also, the secondary nucleation rates are clearly higher for the CU process. For instance, the

317

CC experiments performed with a cooling rate of 0.5 oC/min exhibited a low nucleation

318

rate, lower by several orders of magnitude when compared to CU experiments. Irrespective

319

of the cooling rate, the number of crystals involved or the nucleation rate in CC

320

experiments is always remarkably lower than the experiments performed with CU.

321

Ultrasound continuously generates a vast number of cavition bubbles that burst and collide

322

within the solution which promotes the formation of crystal nuclei.These nucleation results

323

are more detailed in the Supplementary Information, in relation to Figure S1, clearly

324

showing the effect of supersaturation on the nucleation rate and how it is influenced by the

325

sonication. The nucleation rate increases with increasing sonication power. These

326

observations commensurate with the results shown in Figure 8; the higher the nucleation

327

rate the lower is the crystal size. In spite of a much lower D50 value, the CSDW observed

328

in CU crystallization is lower, i.e., the distribution width relative to the mean size is 15 ACS Paragon Plus Environment

Crystal Growth & Design

329

narrower. The Table 2 also reveals that the nucleation rates – primary and secondary –

330

increase with increasing sonication power. As shown in Table 4 the metastable zone width

331

becomes narrower, i.e. ultrasonication does induce primary nucleation at lower

332

supersaturation than the normal MSZW limit. 1.E+10 9.E+09 Number of Crystals(#)/m3

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 37

1 °C/min, CC

8.E+09 0.5 °C/min, CC

7.E+09

0.1 °C/min, CC

6.E+09 5.E+09

60W/kg, 0.5 °C/min, CU

4.E+09

103W/kg, 0.5 °C/min, CU

3.E+09

144W/kg, 0.5 °C/min, CU

2.E+09 1.E+09 0.E+00 0

5000

10000

15000

20000

25000

Time, s

333 334

Figure 9. Comparison of nucleation rates in CC and CU conditions at varying cooling rates

335

and power of ultrasound. The slope of the trend lines showing empirical nucleation rates. Bprim(#/m3/s) 0.1 oc/min 0.5 oc/min 1 oc/min 60W/kg, 0.5oc/min 103W/kg, 0.5oc/min 144W/kg, 0.5oc/min

336

CC CC CC CU CU CU

5.67E+05 1.01E+06 2.36E+06 6.72E+06 8.11E+06 8.66E+06

Bsec (#/m3/s) 1.10E+04 1.70E+05 4.53E+05 3.88E+05 6.37E+05 9.14E+05

Table 2: Nucleation rates in CU and CC experiments calculated from the slopes of

337

the lines presented in Figure 9.

338

3.2 Single Pulse Ultrasonication (SPU):

339

SPU experiments were performed at different power of ultrasound and different

340

sonation time at a constant cooling rate of 0.5 oC/min. The effect of time of ultrasonication 16 ACS Paragon Plus Environment

Page 17 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

341

on the crystal size distribution was studied for three different US times: 2, 7 and 12 min at

342

a constant sonication power of 60 W/Kg. The influence of the US power on the CSD was

343

studied using three different US power, 60, 103 and 144 W/Kg, at a constant US time of 2

344

minutes. Figure 10 shows the crystal size distribution as a function of US time and power

345

in SPU experiments. The higher the US power, the lower the mean size and the CSDW

346

(Table 4). As mentioned earlier, when ultrasound is applied at the maximum

347

supersaturation (MSZW), there is a rapid supersaturation consumption as a consequence of

348

crystal growth of the nuclei formed. Throughtout the rest of the process, the concentration

349

in the solution remains close to the solubility because of the big crystal surface area

350

available. Increase in the power of ultrasound reduced the final mean crystal size and

351

narrowed the crystal size distribution. The width of the distribution is 1.44, 1.24 and 1.19

352

at ultrasonic powers of 60, 103 and 144W/Kg. The influence of the sonication time is not

353

entirely clear, but overall does not appear to be strong. US time of 7 min produced the

354

smallest crystals (the mean sizes produced are 53, 36, and 44 µm for 2, 7 and 12 min of

355

sonication), while the distribution width increases with sonication time ( 1.44, 1.50, and

356

2.12, respectively), as shown Figure 10.

357

17 ACS Paragon Plus Environment

Crystal Growth & Design

358

Figure 10. The comparison of CSD for all SPU experiments at the varying condition

359 360

of US time and power.

9.E+9

Number of crystals(#)/m3

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 37

8.E+9

B

a)

7.E+9

C

D

B C

6.E+9 5.E+9

B

D C

4.E+9 3.E+9

2min SPU 144W/kg 2min SPU 103W/kg 2min SPU 60W/kg

2.E+9 1.E+9 0.E+0 2000

A

2500

3000

361

3500 Time, s

18 ACS Paragon Plus Environment

4000

4500

D

Page 19 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

362 363

Figure 11a) Effect of ultrasonic power and b) time on the nucleation rate in an SPU

364

experiment at a supersaturation ratio of 2.8. The lines indicating the time of sonication

365

from point A (∆Tmax).

366

Bprim(#/m3/s) 1.93E+07 1.43E+07 1.05E+07 1.04E+07 0.98E+07

144W/kg, 2min 103W/kg, 2min 60W/kg, 2min 60W/kg, 7min 60W/kg, 12min 367 368

Table 3: Empirical nucleation rates calculated from the slopes of the lines presented in Figure 11.

369

Figure 11 shows the influence of US power and time on the crystal number

370

generation with time of crystallization. The crystal number here refers to the number of

371

particles (of size > 1 µm) in the suspension as obtained from the FBRM counts.

372

Irrespective of the operating conditions, a rapid and linearly increasing crystal number is

373

followed by a smaller or more pronounced decrease. The initial rapid rise in the crystal 19 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

374

number characterizes the primary nucleation (point A to B in Figure 11), whereas the latter

375

is attributed to a partial dissolution. Secondary processes like secondary nucleation and

376

crystal aggregation (all these features are shown in Figure.S2) will influence the rest of the

377

curve. From the early straight-line part of the number, increase with time the nucleation

378

rates have been determined and are given in Table 3. Obviously, the nucleation rate clearly

379

increases with increasing sonication power, but the influence of the sonication time in SPU

380

is very low. We also note that the nucleation rates in general are higher than those of the

381

CU process given in Table 2. The likely reason for this is that a higher supersaturation has

382

been built up in the SPU experiments before the nucleation commence.

383

At lower US power of 60 W/kg and 103 W/kg the process ends with an increase in

384

the number of crystals (point C to Point D in Figure 11a), while at 144 W/kg the number

385

decreases during the second half of the process. The lower the US power, the lower the

386

primary nucleation rate is and thus the lower is the number of crystals generated in the

387

system during the early period. A lower number generated early may lead to a higher

388

supersaturation later on depending on the cooling rate, thus leading to a stronger secondary

389

nucleation during the second half of the process. When the US power is higher, an exact

390

opposite trend can be realized. At high US power, the primary nucleation rate is high and

391

the number of crystals generated becomes high. With increasing US power the number of

392

cavitational collapses in the solution increases.26 The solution contains a higher number of

393

crystals and surface area leading to lower supersaturation during the second half of the

394

experiment, and thus lower secondary nucleation rate. This means crystal size is inversely

395

proportional to the US power as shown in Figure .10 and also in Table 1 The decrease in

396

the number of particles can be correlated with the frequently encountered phenomenon in

397

this particular system: crystal aggregation that occurs in the absence of sonication27. A

398

similar behaviour of decrease in the number of crystals is observed in experiments

399

performed by longer US time (Figure. 11b). This explains why both the mean size of the

400

final crystals and the CSDW obtained with low US power are higher than the

401

corresponding values for 144 W/kg.

402

As shown in Figure 11a), the nucleation continues beyond the sonication time (2

403

min) for additional 5 minutes and the primary nucleation basically end at the same time 20 ACS Paragon Plus Environment

Page 20 of 37

Page 21 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

404

regardless of the sonication power. As shown in figure 11b) for the experiments with

405

different sonication time, not only the nucleation rate is about the same but also the time

406

period of nucleation. At 2 minutes sonication, the nucleation continues for another 5 min,

407

at 7 minutes sonication it ends at about the same time as the sonication ends and at 12

408

minutes sonication the nucleation ends 5 minutes before the sonication ends. This suggests

409

a mechanism where the effect of the sonication is primarily to initiate the nucleation

410

process, and the ones this initiation has occurred nucleation continues probably up to the

411

point where the supersaturation has decreased sufficiently for this type of nucleation to

412

terminate.28,29 This resembles the behaviour of an avalanche and is actually a description

413

introduced by Mersmann perhaps 30 years ago30. A reasonably short sonication pulse

414

initiates the nucleation, which then proceed and spread out through the solution without the

415

need for further sonication.

416

3.3 Multiple Pulse Ultrasonication (MPU):

417

In Figure 5, it shows one example of the final crystal size distribution and how the

418

temperature varies with the ultrasonication in an MPU experiment. The first sonication

419

pulse is applied at the normal metastable zone limit of 19 oC at the concentration used (first

420

arrow in Figure 5), and the sonication instantly leads to an increased solution temperature

421

by approximately 5 oC (during the US on-cycle period). During the following US off-cycle

422

period, the solution temperature is forced back towards the predetermined linear cooling

423

profile (see Figure 5). It takes approximately 5 minutes for the solution temperature to

424

reach the set linear cooling profile, and based on this observation, the MPU cycle times

425

were fixed to 7 minutes: 2 minutes for the on - cycle period and 5 minutes for the off-cycle

426

period. Thus, at the end of the US off-cycle period, the next US pulse was introduced, and

427

this was repeated in total five times in a row while cooling the solution from 19 oC to 0 oC

428

at a cooling rate of 0.5 oC/min.

21 ACS Paragon Plus Environment

Crystal Growth & Design

2.5

2 (D90 -D10)/D50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 37

1.5

1 CC SPU

0.5

CU MPU

0 0

429 430

20

40

60

80

100

120

140

Crystal Size, µ m

Figure 12. Plot of CDSW vs crystal size for all experiments in Table 4

431

Figure 12 summarizes the characteristics of the product from MPU process experiments in

432

comparison to those of other processing modes and includes experiments in a larger scale.

433

Obviously, the MPU mode is capable of producing small crystals width a narrow size

434

distribution. The crystals obtained from SPU and CU could produce equally small crystals,

435

but the CSDW tends to be higher. This illustrates that by changing the mode of

436

ultrasonication, it is possible to alter the product size distribution to meet specific needs of

437

small crystals with a narrow distribution

438

In Figure 13, the evolution of the process during an MPU experiment is illustrated.

439

Because of the first ultrasonication pulse at the metastable limit, nucleation occurs, and the

440

concentration instantly decreases because of crystal growth. The supersaturation dropped

441

from 31 g/L of solution (S = 2.8; point A) to 20 g/L (S = 1.18; point B) in Figure 13a. In

442

addition, the temperature increases from 20℃ to 25.6 oC, partly due to the release of

443

enthalpy of crystallization (as evidenced from the increase in the number of particles – see

444

Figure 13 b) but primarily because of the ultrasonication as is evidenced from the

445

temperature profile shown in Figure 5. A larger generation of crystals is confirmed by the

446

FBRM counts shown in Figure 13b. During the first US off-cycle period (point B to C) the

447

solution temperature is forced to approach the previously set linear cooling profile (as 22 ACS Paragon Plus Environment

Page 23 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

448

discussed earlier in section 3.1 and 3.2), and there is a slight decrease in the number of

449

particles. The rise in the number of particles and increase in temperature during the US on-

450

cycle period and the drop in the number of particles and decrease in temperature during the

451

US off-cycle period, is repeated for every US pulse, but is most pronounced during the first

452

cycle. During the US on-cycle period, the rise in temperature favours the dissolution of the

453

smaller crystals and any generated secondary nuclei. Deora et.al reported that in

454

dissolution studies, ultrasound leads to faster dissolution of crystals than by thermal

455

treatment because of the heat generation near collapsing cavitational events.26 However, If

456

the supersaturation is sufficiently high, sonication obviously trigger nucleation to occur, in

457

spite of the simultaneous increase in temperature. In MPU mode, nucleation primarily

458

occurs during the first pulse where the supersaturation is very high. This high

459

supersaturation is however rapidly consumed by crystal growth. During the second and

460

third pulse the supersaturation is low and primary nucleation tends to be much weaker or

461

non-existent, while some dissolution of fine crystals will occur due to the temperature

462

increase. This is evidenced from Figure.13b, exposing the increase and decrease in FBRM

463

count during the US on-cycle and off-cycle period, respectively. Despite the fluctuation in

464

FBRM counts, overall (see the linear trendline in Figure 13 b) there is an increase in the

465

number of larger sized particles with respect to time. This effect is more pronounced or at

466

least clearly visible with medium sized crystals (5-10 and 10-20µm) than with the smaller

467

ones (1-3 and 3-5 µm). On the other hand, larger crystals (20-30 and 30-40µm) are the

468

least affected by sonication.31 Instead, their number tends to increase slightly due to the

469

earlier discussed Ostwald ripening effect. This is also supported by the decrease in the

470

CLD (Chord Length Distribution) width (C90-C10) with each US on-cycle period (see the

471

annotation symbol

472

the annotation symbol in figure 13b) due to dissolution of fines and growing of larger

473

crystals at the expense of smaller ones, respectively.

474

According to figure 13a), the solution temperature does not increase enough for it to

475

become undersaturated. We have investigated whether a substantial temperature increase

476

can be recorded locally close to the US probe where the sonication field strength should be

477

the highest. However, we only find about 0.5℃ higher temperature here as compared to the

in figure 13b) followed by an increase during the off-cycle period (see

23 ACS Paragon Plus Environment

Crystal Growth & Design

478

average temperature, and the main reason for this is probably that there is a significant

479

flow of liquid passing by the probe. This confirms the observations stated by Gielena

480

et.al.32 Accordingly, this reduction in the number of very small crystals is either due to

481

normal Oswald ripening perhaps reinforced by a particularly low stability of the fresh

482

crystals nucleated in the sonication field, or due to dissolution because of the very extreme

483

conditions close to the collapsing sonication cavities.

484

40 35

Concentration, g/L

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 37

A

30 25

B

20 15

C

10

F

E

Solubility(g/L)

D

2min, MPU, 144W/Kg

5 0 0 485

10

20 30 Temperature, oC

486

24 ACS Paragon Plus Environment

40

50

Page 25 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

487 488

489 490

Figure .13a) Desupersaturation of MPU experiment b) Change in FBRM count rate with

491

the time of the experiment c) corresponding CLD at specified moments ( )

492

The above-made arguments were further confirmed by measuring the mean crystal size and

493

the CLD (Chord Length Distribution) of the crystals obtained at the end of each US on25 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 37

494

cycle periods ( ). The CSDW decreases with each sonication pulse, whereas the mean size

495

of the crystals tends to increase. This behaviour is typical of a system where the growth of

496

the larger crystal occurs at the expense of smaller crystals, and as such show similarities

497

with the ADNC or DNC (Direct Nucleation Control) techniques.33

498

3.4 Effect of volume of reactor

499

For the scaling-up investigation, Single Pulse Ultrasonication (SPU) and Multiple Pulse

500

Ultrasonication (MPU) experiments were carried out also in 3L scale (CR2) using the same

501

cooling rate (0.5 ℃/min). The same sonication probe (20 kHz) was used as in the 0.5 L

502

experiments and thus, US power/mass is much lower (144 W/kg for CR1 and 24W/kg for

503

CR2). For comparison purpose, also conventional cooling crystallizations

504

performed in the larger volume. Figure 14 shows a comparison of the concentration versus

505

temperature in CR1 and CR2 operated in SPU mode and in the CC mode. Obviously, the

506

principal influence of the SPU in the larger scale is quite similar as in the small scale. The

507

ultrasonication leads to an earlier onset of nucleation and an increase in temperature of the

508

solution. However, the increase in temperature in CR2 is less than in CR1 since the

509

sonication power per unit mass is 6 times lower.

510

As shown in Table 4, in the absence of sonication (CC mode), the product size distribution

511

is essentially independent of crystallization volume – the mean size and the CSDW value

512

are essentially unchanged. Using SPU the mean size from CR2 at 24 W/kg operates lower

513

than that from CR1 at higher specific US power (60W/kg), while the CSDW is essentially

514

unchanged. In addition, a trend of decreasing mean size with increasing US time is found

515

in CR2, but there is no tendency of the CSDW to increase. As stated above for the CR1

516

experiments, primary nucleation continue beyond the SPU sonication being turned off. The

517

fact that the size distribution from CR2 is relatively similar to that from CR1 at only 40 %

518

of the specific power input suggests that a significant part of the primary nucleation is not

519

only taking place after that the sonication has been turned off but is also spreading in space

520

from an initial burst directly generated by the sonication (because otherwise, the product

521

size from CR2 should be larger than that from CR1). This further supports an analogy with

522

the behaviour of an avalanche. 26 ACS Paragon Plus Environment

(CC) were

Page 27 of 37

40 35 Concentration, g/L

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

30 25 20 15

Solubility in g/L 0.5 °C/min,2min SPU, CR1 0.5 °C/min,2min SPU, CR2 0.5 °C/min, CC, CR1 0.5 °C/min, CC, CR2

10 5 0 0

10

523

20 30 Temperature, oC

40

50

524

Figure 14. Comparison of desupersaturation in CR1 and CR2 crystallizers for ultrasonic

525

and non-ultrasonic one pulse mode.

526

Unfortunately, the number of actual crystals present in CR2 cannot be calculated using the

527

empirical correlation shown in Figure 4, since this is only valid for the flow conditions in

528

the small crystallizer.

529

Multipulse sonication (MPU) was also applied in the same way as in CR1, i.e. 2 min of

530

sonication followed by 5 min without ultrasound. As was the case in the 0.5 L crystallizer

531

(please see section 3.1 and 3.2), the MPU mode produced a narrower CSD distribution than

532

the SPU mode, and in the CR2 crystallizer the low mean size is maintained but not further

533

reduced from that obtained by SPU, Figure 15. The rise in temperature observed during the

534

2 min of ultrasonication is again less than that in CR1. The solution concentration stayed

535

close to the solubility curve already after two successive pulses of ultrasound. Ultrasound

536

applied after the first two pulses introduced Ostwald ripening effect as observed earlier in

537

CR1 (Figure 13b).

27 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

538 539

Figure 15. Comparison of crystal size distribution of SPU and MPU mode of ultrasound in

540

CR2.

541

3.5 Effect of sonication on Paracetamol/water and Ibuprofen/Isopropanol systems.

542

MPU crystallization experiments were performed with two other model compounds;

543

paracetamol and ibuprofen with a solution concentration corresponding to the saturation

544

temperature of 50 and 20 oC, respectively. For paracetamol, 17 g of paracetamol was

545

dissolved in 500mL (CR1) of water at 50 oC (solubility of paracetamol: 34 g/L at 50 oC).34

546

The solution was heated up to 55 oC (and maintained at this temperature for 10 min) to

547

dissolve the crystals completely and then cooled and maintained at the saturation

548

temperature of 50 oC. The solution was then allowed to cool at a constant cooling rate of

549

0.5 oC /min until it reached the final temperature of 5 oC. The first pulse of US was applied

550

precisely at the metastable zone limit at 41.35 oC (equal to ∆Tmax = 8.65 oC determined

551

experimentally in this work in the absence of the US) for 2 minutes and then followed by a

552

5-minute off-time. The experiment then involves nine additional equal US pulses.

553

The experiments with Ibuprofen was performed by dissolving 63 g of Ibuprofen in 500ml

554

(CR1) of acetonitrile at 20 oC.35 The solution was heated up to 25 oC and maintained at

555

this temperature for 10 min. The solution was then allowed to cool at a cooling rate of 0.5

556

o

C/min until it reached 0 oC. Five sonication pulses were introduced, the first pulse of US 28 ACS Paragon Plus Environment

Page 28 of 37

Page 29 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

557

was applied precisely at the metastable limit, i.e. 15 oC (equal to ∆Tmax = 5 oC ) and

558

throughout the experiment the US on-time period was 2 min followed by a US off-time

559

period of 5 min. For comparison, experiments were also performed in the absence of US

560

(CC experiments) for both compounds.

561

Figure 16 shows the microscopic images and the crystal size distribution of the final

562

crystals obtained from MPU experiments, and the results are compared with the CC mode

563

experiments in Table 4. Clearly, the MPU mode produced much smaller crystals and more

564

narrow distributions than the CC mode. For paracetamol the mean size decreases from 158

565

µm (CC) to 45 µm (MPU), and the CSDW from 1.5 (CC) to 1.12 MPU). For ibuprofen the

566

mean size decreased from 83 µm (CC) to 21 µm (MPU), and the CSDW from 1.45 (CC) to

567

1.14 (MPU). As opposed to previous work (Bhangu et al)23, only the stable polymorph of

568

paracetamol was obtained in the present study.

569 570 571 572

29 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

a)

b)

573 574

Figure 16.a) Crystal size distribution comparison of Paracetamol and Ibuprofen b)

575

microscopic image of final crystals of paracetamol and Ibuprofen from MPU mode.

576

From the results of the present work, it appears as if the principal influence of sonication

577

on the product crystal size distribution is the same for different systems.A single-pulse

578

(SPU) and a multiple pulse sonication mode (MPU) obviously can produce small crystals

579

with a narrow size distribution and this finding seems to hold true for all the three different

580

systems studied in this work. The detailed sonication scheme needs to be further

581

researched. It also appears to be quite clear that scaling-up to obtain an equal product

582

crystal size do not require equal US power/unit volume.

583

Conclusions

584

In this work it is shown that ultrasonication can be used during batch cooling

585

crystallization to control the product size distribution into low mean sizes and narrow

586

distributions. There appears to be no benefit from the size distribution point of view in

587

continuously exposing the solution to sonication, but short sonication pulses strategically

588

used will be more efficient, reducing energy consumption and increasing sonication probe

589

life time. In single pulse sonication mode, sonication power is an important factor in 30 ACS Paragon Plus Environment

Page 30 of 37

Page 31 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

590

influencing the product properties but sonication time much less so. Multiple pulse

591

sonication appears to more consistently produce small crystals with a narrow particle size

592

distribution, but more work is needed to detail the optimum sonication strategy. The

593

nucleation generated by the sonication appears to follow an avalanche mechanism. The

594

sonication initiates the nucleation, which then progress by itself in time and space. As a

595

consequence scaling-up of a sonicated process will not require the same sonication power

596

per unit volume to reach the same product size distribution.

597

Supporting Information: The supporting information available for free of charge on the

598

ACS publication site at DOI:

599

Empirical Nucleation rate at measured corresponding supersaturation ratio curves for CC

600

and CUexperiments, FBRM count data and corresponding CLD at pinned locations in the

601

process of crystallization at SPU conditions.

602

Funding: This work has been supported by the Science Foundation Ireland (Grants SFI

603

SSPC2 12/RC/2275)

604

REFERENCES

605 606

(1)

Duffy, D.; Barrett, M.; Glennon, B. Novel, Calibration-Free Strategies for Supersaturation

607

Control in Antisolvent Crystallization Processes. Cryst. Growth Des. 2013, 13, 3321–

608

3332.

609

(2)

Aamir, E.; Nagy, Z. K.; Rielly, C. D. Optimal Seed Recipe Design for Crystal Size

610

Distribution Control for Batch Cooling Crystallisation Processes. Chem. Eng. Sci. 2010,

611

65, 3602–3614.

612

(3)

Yu, Z. Q.; Chow, P. S.; Tan, R. B. H. Application of Attenuated Total Reflectance-Fourier

613

Transform Infrared (ATR-FTIR) Technique in the Monitoring and Control of Anti-

614

Solvent Crystallization. Ind. Eng. Chem. Res. 2006, 45, 438–444.

615 616

(4)

Bakar, M. R. A.; Nagy, Z. K.; Rielly, C. D. Seeded Batch Cooling Crystallization with Temperature Cycling for the Control of Size Uniformity and Polymorphic Purity of 31 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Sulfathiazole Crystals. Org. Process Res. Dev. 2009, 13, 1343–1356.

617 618

(5)

Abu Bakar, M. R.; Nagy, Z. K.; Saleemi, A. N.; Rielly, C. D.; Bakar, M. R. A.; Nagy, Z.

619

K.; Saleemi, A. N.; Rielly, C. D.; Abu Bakar, M. R.; Nagy, Z. K.; et al. The Impact of

620

Direct Nucleation Control on Crystal Size Distribution in Pharmaceutical Crystallization

621

Processes. Cryst. Growth Des. 2009, 9, 1378–1384.

622

(6)

Rev. Chem. Biomol. Eng. 2012, 3, 55–75.

623 624

Nagy, Z. K.; Braatz, R. D. Advances and New Directions in Crystallization Control. Annu.

(7)

Ruecroft, G.; Hipkiss, D.; Ly, T.; Maxted, N.; Cains, P. W.; Technology, C.; Ox, O.;

625

Kingdom, U.; Ruecroft, G.; Hipkiss, D.; et al. Sonocrystallization : The Use of Ultrasound

626

for Improved Industrial Crystallization Abstract : Org. Process Res. Dev. 2005, 9, 923–

627

932.

628

(8)

Hightower, J. R. Ultrasonics, the Low- and High-Intensity Applications, Dale Ensminger,

629

Marcel Dekker, Inc., New York, N. Y.(1973). 570 pages.$24.50. AIChE J. 1973, 19, 879–

630

880.

631

(9)

Sonocrystallization. Ultrason. Sonochem. 2011, 18, 345–355.

632 633

Nalajala, V. S.; Moholkar, V. S. Investigations in the Physical Mechanism of

(10)

Kaur Bhangu, S.; Ashokkumar, M.; Lee, J. Ultrasound Assisted Crystallization of

634

Paracetamol: Crystal Size Distribution and Polymorph Control. Cryst. Growth Des. 2016,

635

16, 1934–1941.

636

(11)

Habit. Ultrason. Sonochem. 2007, 14, 143–152.

637 638

(12)

(13)

Nii, S.; Takayanagi, S. Growth and Size Control in Anti-Solvent Crystallization of Glycine with High Frequency Ultrasound. Ultrason. Sonochem. 2014, 21, 1182–1186.

641 642

Maher, A.; Croker, D.; Rasmuson, Å. C.; Hodnett, B. K. Solubility of Form III Piracetam in a Range of Solvents. J. Chem. Eng. Data 2010, 55, 5314–5318.

639 640

Bund, R. K.; Pandit, a B. Sonocrystallization: Effect on Lactose Recovery and Crystal

(14)

Hatkar, U. N.; Gogate, P. R. Chemical Engineering and Processing : Process 32 ACS Paragon Plus Environment

Page 32 of 37

Page 33 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

643

Intensification Process Intensification of Anti-Solvent Crystallization of Salicylic Acid

644

Using Ultrasonic Irradiations. Chem. Eng. Process. Process Intensif. 2012, 57–58, 16–24.

645

(15)

Sonocrystallization. J. Cryst. Growth 2010, 312, 961–966.

646 647

Boels, L.; Wagterveld, R. M.; Mayer, M. J.; Witkamp, G. J. Seeded Calcite

(16)

Barrett, M.; Hao, H.; Maher, A.; Hodnett, K.; Glennon, B.; Croker, D. In Situ Monitoring

648

of Supersaturation and Polymorphic Form of Piracetam during Batch Cooling

649

Crystallization. Org. Process Res. Dev. 2011, 15, 681–687.

650

(17)

Kim, J.-W.; Kim, J.; Lee, K.-D.; Koo, K.-K. Evaluation of Nucleation Rate by in-Situ

651

Focused Beam Reflectance Measurement in an Unseeded Batch Cooling Crystallization.

652

Cryst. Res. Technol. 2013, 48, 1097–1105.

653

(18)

Bouyarmane, H.; Gouza, A.; Masse, S.; Saoiabi, S.; Saoiabi, A.; Coradin, T.; Laghzizil, A.

654

Nanoscale Conversion of Chlorapatite into Hydroxyapatite Using Ultrasound Irradiation.

655

Colloids Surfaces A Physicochem. Eng. Asp. 2016, 495, 187–192.

656

(19)

with Gassing of Adipic Acid. Chem. Eng. Sci. 2010, 65, 1016–1027.

657 658

Wohlgemuth, K.; Ruether, F.; Schembecker, G. Sonocrystallization and Crystallization

(20)

Bund, R. K.; Pandit, A. B. Rapid Lactose Recovery from Paneer Whey Using

659

Sonocrystallization: A Process Optimization. Chem. Eng. Process. Process Intensif. 2007,

660

46, 846–850.

661

(21)

Miyasaka, E.; Takai, M.; Hidaka, H.; Kakimoto, Y.; Hirasawa, I. Effect of Ultrasonic

662

Irradiation on Nucleation Phenomena in a Na2HPO4.12H2O Melt Being Used as a Heat

663

Storage Material. Ultrason. Sonochem. 2006, 13, 308–312.

664

(22)

Jordens, J.; Gielen, B.; Braeken, L.; Van Gerven, T. Determination of the Effect of the

665

Ultrasonic Frequency on the Cooling Crystallization of Paracetamol. Chem. Eng. Process.

666

Process Intensif. 2014, 84, 38–44.

667

(23)

Kaur Bhangu, S.; Ashokkumar, M.; Lee, J. Ultrasound Assisted Crystallization of

668

Paracetamol: Crystal Size Distribution and Polymorph Control. Cryst. Growth Des. 2016,

669

16, 1934–1941. 33 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

670

(24)

Guo, Z.; Jones, a. G.; Hao, H.; Patel, B.; Li, N. Effect of Ultrasound on the

671

Heterogeneous Nucleation of BaSO[sub 4] during Reactive Crystallization. J. Appl. Phys.

672

2007, 101, 54907.

673

(25)

Zisu, B.; Sciberras, M.; Jayasena, V.; Weeks, M.; Palmer, M.; Dincer, T. D.

674

Sonocrystallisation of Lactose in Concentrated Whey. Ultrason. Sonochem. 2014, 21,

675

2117–2121.

676

(26)

Deora, N. S.; Misra, N. N.; Deswal, A.; Mishra, H. N.; Cullen, P. J.; Tiwari, B. K.

677

Ultrasound for Improved Crystallisation in Food Processing. Food Eng. Rev. 2013, 5, 36–

678

44.

679

(27)

Kim, J.-M.; Chang, S.-M.; Kim, K.-S.; Chung, M.-K.; Kim, W.-S. Acoustic Influence on

680

Aggregation and Agglomeration of Crystals in Reaction Crystallization of Cerium

681

Carbonate. Colloids Surfaces A Physicochem. Eng. Asp. 2011, 375, 193–199.

682

(28)

Solutions. Chem. Eng. Res. Des. 2015, 97, 68–76.

683 684

Kim, H. J.; Yeo, S. Do. Liquid Antisolvent Crystallization of Griseofulvin from Organic

(29)

Kitayama, H.; Yoshimura, Y.; So, M.; Sakurai, K.; Yagi, H.; Goto, Y. A Common

685

Mechanism Underlying Amyloid Fibrillation and Protein Crystallization Revealed by the

686

Effects of Ultrasonication. Biochim. Biophys. Acta - Proteins Proteomics 2013, 1834,

687

2640–2646.

688

(30)

Mersmann, A. Crystallization Technology Handbook; 2001.

689

(31)

Lakerveld, R.; Verzijden, N. G.; Kramer, H.; Jansens, P.; Grievink, J. Application of

690

Ultrasound for Start-up of Evaporative Batch Crystallization of Ammonium Sulfate in a

691

75-L Crystallizer. AIChE J. 2011, 57, 3367–3377.

692

(32)

Gielen, B.; Kusters, P.; Jordens, J.; Thomassen, L. C. J.; Van Gerven, T.; Braeken, L.

693

Energy Efficient Crystallization of Paracetamol Using Pulsed Ultrasound. Chem. Eng.

694

Process. Process Intensif. 2017, 114, 55–66.

695 696

(33)

Saleemi, A. N.; Steele, G.; Pedge, N. I.; Freeman, A.; Nagy, Z. K. Enhancing Crystalline Properties of a Cardiovascular Active Pharmaceutical Ingredient Using a Process 34 ACS Paragon Plus Environment

Page 34 of 37

Page 35 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

697

Analytical Technology Based Crystallization Feedback Control Strategy. Int. J. Pharm.

698

2012, 430, 56–64.

699

(34)

Eng. Data 1999, 44, 1391–1395.

700 701 702

Granberg, R. A.; Rasmuson, Å. C. Solubility of Paracetamol in Pure Solvents. J. Chem.

(35)

Khalifeh, I. THERMODYNAMIC EVALUATION OF IBUPROFEN SOLUBILITY IN AQUEOUS AND NON-AQUEOUS COSOLVENT SYSTEMS, Purdue University, 2000.

703

35 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

704 705

Table 4: Experimental results at varying conditions on crystal size, MZSW and CSDW.

706 36 ACS Paragon Plus Environment

Page 36 of 37

Page 37 of 37

For Table of Contents Use Only

707 708

Manuscript Title:

709

Controlling the product crystal size distribution by strategic application of ultrasonication

710

Author List:

711

Kiran A. Ramisetty, Åke C. Rasmuson

712

TOC graphic, and synopsis. 2.5

2

(D90-D10)/D50

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

1.5

1 CC SPU

0.5

CU MPU

0 0

713 714

20

40

60

80

100

120

140

Crystal Size, µ m

Changing the ultrasonication mode can tailor the crystal size and its distribution width.

715

37 ACS Paragon Plus Environment