Risk Considerations on Developing a Continuous ... - ACS Publications

Jun 15, 2017 - Here, we focus on how to develop a continuous crystallization system, from the perspective of process engineering, and the related risk...
0 downloads 5 Views 7MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Article pubs.acs.org/OPRD

Risk Considerations on Developing a Continuous Crystallization System for Carbamazepine Xiaochuan Yang, David Acevedo, Adil Mohammad, Naresh Pavurala, Huiquan Wu, Alex L. Brayton, Ryan A. Shaw, Mark J. Goldman, Fan He, Shuaili Li, Robert J. Fisher, Thomas F. O’Connor, and Celia N. Cruz* Office of Pharmaceutical Quality, CDER, FDA, 10903 New Hampshire Avenue, Silver Spring, Maryland 20993-0002, United States S Supporting Information *

ABSTRACT: Continuous manufacturing (CM) is an emerging technology in the pharmaceutical manufacturing sector, and the understanding of the impact on product quality is currently evolving. As the final purification and isolation step, crystallization has a significant impact on the final physicochemical properties of drug substance and is considered a critical process step in achieving the continuous manufacturing of drug substances. Although many publications previously focused on various innovative techniques to continuously make crystals with desired properties, engineering difficulties such as system design, automation, and integration with process analytical technology (PAT) tools have not been thoroughly discussed. Here, we focus on how to develop a continuous crystallization system, from the perspective of process engineering, and the related risk considerations on product quality. Specifically, we designed and built an automated two-stage mixed suspension mixed product removal (MSMPR) crystallization platform for a model compound (carbamazepine, CBZ) that exhibits multiple polymorphs. The crystallization process includes the integration of PAT tools (online Raman microscopy and focused beam reflectance microscopy, FBRM) for real-time monitoring. A series of case studies were done to evaluate the performance of the continuous system and PAT tools. Specifically, the drawing schemes, slurry transport, and variations on process variables are considered as the three key risk areas for continuous crystallization process development. Our proof-of-concept continuous crystallization system uses feedback/feedforward controls to achieve constant levels in crystallizers, a centralized automation program coded in LabVIEW, and PAT monitoring for polymorphs and particle size distribution (Raman and FBRM). To the best of our knowledge, this is also the first study on continuous crystallization of carbamazepine for form III and its polymorphic transition (between form II and form III). lizers,15−19 which are suitable for faster kinetics and better temperature control. However, tubular crystallizers are more prone to have crystal fouling, clogging, settling issues, and poor mixing, which are all hurdles for industrial implementation.20 In addition to these two common setups, there are several novel crystallization approaches that potentially can be used for continuous crystallization, such as microfluidics,18 impinging jet,21 film crystallization,22 electro-spraying,23 nanospraying,24 and crystallization within porous media.25−27 Research efforts mentioned above mainly focus on novel technologies or methodologies with little discussion from the perspective of process robustness and product quality. In this paper, we would like to emphasize the importance of several scientific and engineering considerations from the perspective of quality risk management, specifically what are the potential quality risks when developing a two-stage MSMPR continuous crystallization system, and how to minimize these risks through process design, monitoring, and control. To demonstrate these concepts, a laboratory scale continuous crystallization platform is presented as a proof-of-concept for the development and implementation of automation, monitoring, and control strategies.

1. INTRODUCTION The pharmaceutical industry nowadays is undergoing a shift from batch to continuous manufacturing.1−3 While an end-toend integration process from synthesis to final dosage forms still remains challenging, continuous manufacturing of drug product part was demonstrated successfully with the commercial production of Orkambi tablets by Vertex Pharmaceuticals4,5 and Prezista tablets by Janssen Pharmaceutica.6 As for an integrated line on continuous manufacturing of drug substance part, no commercial manufacturing line has been reported yet, and there are still several roadblocks to be solved. One major roadblock is a continuous crystallization step, though several advances have been made in flow chemistry for pharmaceutical drug substances.7−9 Crystallization is a vital step in drug substance manufacturing process, which not only purifies APIs, but also solidifies APIs from liquid phase for downstream processing. Research efforts on continuous crystallization have increased significantly over the past few years. Currently, there are two popular approaches to achieve continuous crystallization. One approach is the cascade of MSMPR crystallizers, where one or more well-stirred crystallizers are connected in sequence for multiple-stage operations.10−14 This design is suitable for systems that require longer residence times, but its strong agitation and nonuniform temperature profile may cause problems for particle size and polymorphism control. The other approach is tubular crystalThis article not subject to U.S. Copyright. Published 2017 by the American Chemical Society

Received: March 30, 2017 Published: June 15, 2017 1021

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

Figure 1. Process and instrumentation diagram for the continuous crystallizer.

Table 1. Continuous Crystallizer Equipment Description number

name(s)

CR 1

Crystallizer 1

CR 2

Crystallizer 2

P1 P2 P3 P4

Pump 1, feed pump Pump 2, intermediate pump Pump 3, product pump Pump 4, anti-solvent pump, dosing pump feed tank product tank anti-solvent tank

V1 V2 V3

description A thermostatted, stirred tank (500 mL) produced by Mettler Toledo as part of the EasyMax 402 Synthesis Workstation. EasyMax 402 allows for temperature and impeller speed control. A thermostatted, stirred tank (500 mL) produced by Mettler Toledo as part of the EasyMax 402 Synthesis Workstation. EasyMax 402 allows for temperature and impeller speed control. A peristaltic pump produced by Ismatec as part of the Reglo Analog line (ISM829). A peristaltic pump produced by Ismatec as part of the Reglo Analog line (ISM829). A peristaltic pump produced by Ismatec as part of the Reglo Analog line (ISM829). A 50 mL (alternatively 10 mL) syringe pump produced by Mettler Toledo as an accessory to EasyMax 402. A bottle (which may be used with or without a sealed top) that feeds Crystallizer 1 via the feed pump. A bottle (which may be used with or without a sealed top) that draws from Crystallizer 2 via the product pump. A bottle (which may be used with or without a sealed top) that feeds Crystallizer 1 or Crystallizer 2 via the antisolvent pump.

2. EXPERIMENTAL SECTION 2.1. Materials. Carbamazepine (CBZ) was chosen as the model compound throughout the various studies performed in this work. CBZ possesses four polymorphs (I, II, III, and IV).28 Polymorph III is the most stable form under room temperature. The commercial CBZ purchased from Sigma-Aldrich was tested by X-ray powder diffraction (XRPD), and the results showed polymorph III. Ethanol 200 proof (Decon Laboratories) was used as solvent throughout this work. The solubility of CBZ in ethanol at 295.80 K is 2.2 g/mL ethanol. Deionized water (Millipore System) was used as antisolvent for the polymorphic transformation experiments. Air was used for the level indicator and went through an ethanol buffer bottle to get saturated with ethanol first before entering the system. 2.2. Experimental Setup. The two-stage MSMPR system and associated instrumentation is shown in Figure 1; the equipment contained is described in Table 1. Feed solution from the feed tank (V 1) is fed through the feed pump (P 1) into crystallizer 1 (Cr 1). Similarly, slurry may be removed from crystallizer 1 by the intermediate pump (P 2) and sent to crystallizer 2 (Cr 2) before being removed by the product pump (P 3) to a product tank (V 2). In addition, antisolvent from the antisolvent tank (V 3) may be pumped to either crystallizer 1 or crystallizer 2 via the dosing pump (P 4), not shown in the figure. Crystallizer temperature, crystallizer stirring speed, and antisolvent flow are controlled by Mettler Toledo’s EasyMax 402 System. The temperature in the crystallizer was measured using a Pt100 thermocouple. The system uses a Peltier device for cooling and electric system for heating the water that passes through the

jacket system. Pumps P 1 through P 3 are controlled manually through dials on the pump housing or automatically through the control system described as below. 2.3. Automation Control System. The feedback and feedforward control system for the experimental setup is briefly discussed in this section. For detailed engineering documentation, see the Supporting Information (SI). The control strategy for flow and level control system can also be seen in Figure 1: (1) Computer control of peristaltic pumps P 1, P 2, and P 3 via FC01, FC02, and FC03; (2) Level indication in Cr 1 and Cr 2; (3) Level control of Cr 1 and Cr 2 via manipulation of P 1 and P 3. In this system total process throughput is established by the flow controller on P 2 (FC02), level in Cr 1 is manipulated via FC01, and level in Cr 2 is manipulated via FC03. This configuration ensures level disturbances do not propagate between crystallizers. This software front-end for the flow and level control system has been implemented via LabVIEW, a development environment by National Instruments. 2.3.1. Pump Flow Control. The analog pumps are controlled using LabVIEW via a data acquisition unit to allow communication between the pumps and LabVIEW. To facilitate this communication, custom circuitry has been designed, tested, and constructed. This circuitry allows LabVIEW to send the following signals to the pumps: • Remote condition (manual or computer control); • Run status (on or off); • Pump direction (forward or reverse); 1022

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

collection frequency, the filtration system has two parallel vacuum filter setups; while one of these setups is actively collecting sample, the other setup is being harvested and reestablished. At the end of the collection period, the setups are switched via two valves controlled by the control system implemented in LabVIEW. This system is illustrated in Figure 3. This type of setup is inspired by industrial practices in catalytic bed reactors or pressure swing adsorption.29 2.5. Off-Line Characterization. The particle size distribution (PSD) and the polymorphic form of samples obtained throughout the various case studies were determined using offline characterization techniques. The PSD was measured using a Sympatec HELOS laser diffractometer with a RODOS disperser. An optimal concentration of 0.1% was setup as the trigger condition and a feed rate of 50% and 1 bar. The HELOS sensor can detect particles in the range of 0.1−875 μm. X-ray powder diffraction (XRPD) was used to identify the polymorph obtained throughout the continuous process. The CBZ polymorphs form II and III (form III is more stable than form II) can be easily identified through their XRPD spectra. Each polymorph has characteristic peaks at distinct values of 2θ.27 The XRD was performed by a Bruker D8 Advance (Bruker AXS, Madison, Wisconsin) diffractometer equipped with the LYNXEYE scintillating detector and Cu Kα (λ = 1.5405 Å) using a voltage of 40 kV and an anode current of 40 mA. Before measurement, the instrument functionality was checked using corundum as an external standard. About 500 mg of sample was placed in the sample holder, and three replicate diffractograms of each sample were collected over 2θ range of 4−40° with an increment of 0.00251 at 0.2 s per step (7181 total steps). Sample holder was rotated during run to minimize preferred orientation and to get average diffractogram of the sample. The XRPD operation, data collection, and data analysis were achieved through Diffract.Suite (V2.2). 2.6. Process Analytical Technology (PAT) Tools. Ensuring that a process stays within a state of control is crucial to achieving consistent quality. Interrogating the quality of the materials in process via process monitoring and PAT tools is a key element of establishing a state of control. Raman spectrophotometer and focused beam reflectance measurement (Mettler Toledo, FBRM) probes were chosen as PAT tools to monitor the process and the product quality. Figure 4 illustrates

• Pump speed. A LabVIEW program has been constructed, and a userinterface for pump flow control has been designed. The pump calibration has been implemented in the LabVIEW pump control. 2.3.2. Level Indication and Control. Figure 2 illustrates the bubbler level indication scheme for each of the crystallizer tanks.

Figure 2. Schematic of bubbler level indication.

A bubbler operates by sending a stream of air to the base of a liquid containing tank. The pressure difference between the gas in the bubbler and the headspace of the tank gives an indication of the liquid level, via the following equation, where ΔP is the differential pressure, ρ is the liquid density, g is the force due to gravity, and h is the height. ΔP = ρgh ⇒ h =

ΔP ρg

In order for the differential pressure signal to be successfully integrated into LabVIEW, it is filtered and amplified via a strain gauge transmitter. The filtered and amplified signal from transmitter is sent to LabVIEW via the data acquisition unit. Using this signal, LabVIEW adjusts the pump flow rates via a PI controller to maintain the level at its desired set point. 2.4. Semicontinuous Filtration Setup. The semicontinuous filtration device filters in a traditional batch fashion but does so intermittently. The designed filtration system targets a collection frequency of one residence time. To achieve this

Figure 3. Periodic batch filtration system schematic. 1023

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

Figure 4. A photo of the continuous crystallization system. The electronics of the automation control system are mostly mounted on the back rack. In the Easymax system, the Raman and FBRM probes were located in the Cr 2 (the left crystallizer). Contrary to Figure 1, the flow direction is from right to left.

Figure 5. Ishikawa diagram that shows possible causes for process failure in a continuous crystallization process.

the final setup of the whole continuous crystallization system: a two-stage continuous crystallization unit with engineered flow and level control, implemented downstream semicontinuous filtration, and integrated PAT tools. 2.6.1. Focused Beam Reflectance Measurement (FBRM). FBRM is based on a large sample of 2D projections of the particles with random rotations.30,31 A count of these events as a function of the chord length yields a chord length distribution (CLD), which can be correlated to particle size distribution in real time. However, significant factors can impact this correlation such as crystal shape and orientation.33,34 Because the particles are scanned in one dimension as they flow past the probe window, the CLD count is biased toward lower size ranges as the shape tends away from ideal spherical. For needle like crystals or rod-like crystals the dominant chord length will be closer to its minor axis length, and as a result, FBRM is not very sensitive to axial growth of particles.32 Thus, without detailed knowledge of particle shape, CLD cannot be mapped onto a crystal size distribution with sufficient confidence.31 The system must also be properly calibrated for accuracy. The CLD can, however, track

how one-dimensional size and count change during a process.33,34 The data collected can be trended over time. Focusing on the qualitative analysis, it can capture events such as nucleation, crystal growth, agglomeration, and breakage though tracking of CLD evolution.12,35 Therefore, the utility of FBRM as a quantitative PAT tool for this system needs to be demonstrated with increased efforts in calibration and correlation development for the particulate system. The FBRM was installed and integrated with iControl, allowing for simultaneous collection and display of both the parameters of EasyMax operation (i.e., temperature andstir speed) with FBRM data. With FBRM as a PAT, trends related to particle size distribution can be tracked real-time to better understand, and eventually optimize, and control the crystallization system. 2.6.2. Raman Spectroscopy. Raman spectroscopy is used to distinguish polymorphs in the system. The Raman analyzer used is the Raman Rxn2 system by the Kaiser Optical Systems, Inc.: a four channel analyzer that operates sequentially allowing both fast analysis per channel and programmable channel interrogation. It supports both off-line analysis and on-line analysis 1024

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

Figure 6. Stir speed experimental results as presented by FBRM raw data. The purple line is the stir speed. Other lines for counts represent different ranges of chord sizes.

method. Measurements were collected and saved every 5 s for the primary distribution mode. First, a solution of carbamazepine (CBZ) was kept at 20 °C in CR 2 and stirred at a constant stir speed of 300 rpm The FBRM was tested at different locations. No significant difference was observed at different locations, as long as the probe tip was sufficiently submerged. The impact of stir speed on the FBRM was evaluated by varying the stir speed between 100 and 500 rpm in CR 2 shown in Figure 6. Under these saturated conditions, there is no crystal growth; i.e., the crystal size distribution should be almost constant over time (Ostwald ripening and crystal breakage due to agitation may exist). From the experiment, it is clear that the FBRM is quite responsive to changes in stir speed. As explained above, the FBRM measures the numbers of particles per rotation of the laser beam. Therefore, the number of particles observed by the FBRM varies as the mixing behavior in the vessel changes. This can be observed clearly when significant step changes are applied; a significant decrease in counts for all bin ranges is observed when the stir speed decreases to 100 rpm. Low stir speeds increase the probability of settling of crystals in the system, which results in a significant decrease in particle counts. The use of FBRM as a tool to evaluate the impact of stir speed on homogeneity of the system has been presented in literature.36 Nonetheless, this is a qualitative observation for a fixed volume, solid concentration, and crystal size distribution. The interpretation of the changes observed, in terms of particle size and impact on product quality, will be evaluated with additional development. 3.1.2. Raman. Form II and form III CBZ were made in-house and tested with off-line Raman spectroscopy (Figure 7). The characteristic peaks of the two polymorphs of CBZ correspond closely to the peak locations that have been characterized in literature.37,38 This is especially useful for monitoring the transition from form III to form II when antisolvent (water) is added to the system. The peak area ratio of the peak at 1050 cm−1 to the peak at 1027 cm−1 is used to distinguish form II and form III,38,39 as shown in Figure 7. In addition, form II also has two small characteristic peaks at 257 and 391 cm−1. The peak ratio will be used throughout this work since the two small

mode with a 785 nm laser beam. The initialization of Raman requires about 1 h warming up time. Proper coverage of the system such as aluminum foils is required because the Raman probe is sensitive to visible light.

3. RESULTS AND DISCUSSION Unlike many other unit operations in pharmaceutical manufacturing process, continuous crystallization (excluding some novel methods mentioned in the introduction) requires handling a mixture of both solid and liquid phases. The two-phase flow creates many challenges to process robustness, such as clogging in tubing and settling of large particles in low-velocity flow zones. These process failure modes, if not properly addressed, may affect product quality. Here, as shown in Figure 5, we will discuss these potential process risks and our related findings in the following sections. To achieve a robust continuous crystallization process, we have identified six areas of focus during process development: (1) Drawing scheme; (2) slurry transport; (3) automation system for level control; (4) PAT tools; (5) materials; (6) process variables. Based on our experience, Areas 1, 2, and 6 should be considered as the key risk areas for continuous crystallization process design and are discussed in detail in the following sections. Without appropriate considerations on these three areas, a fully functional integrated system is hardly possible. Meanwhile, Areas 3, 4, and 5 are crucial additions to the system to ensure that the system is operated in a monitored state-ofcontrol. Considerations for Area 3 have been described in the above Experimental Setup section and in the Supporting Information. In addition, regarding Area 4, the details of the PAT method development and validation will be discussed in a separate manuscript. For Area 5, more studies will be conducted on the effect of seeded crystals (different polymorphs and size distributions) on the process dynamics of continuous crystallization later. 3.1. PAT Tools. 3.1.1. Focused Beam Reflectance Measurement (FBRM). A series of batch experiments were conducted with the FBRM to determine the sensitivity and precision of the 1025

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

Figure 7. Off-line Raman spectrum of CBZ forms II and form III powders. The dashed squares highlight the characteristic peaks.

characteristic peaks (i.e., 257 and 391 cm−1) are at the end of the spectrum where the measurement noise increases. 3.1.3. Case Study: Polymorph Transition in a Single-Stage MSMPR Setup. To evaluate the ability of the system to monitor process status and detect disturbances in real-time with the PAT tools (specifically polymorphic change from form III to form II in this case study), an experimental run using antisolvent in a singlestage MSMPR setup is conducted. Raman and FBRM were implemented into the crystallizer to monitor real-time particle counts and polymorph transition from form III to form II. The system started at 20 °C with form III partially dissolved in ethanol (a total concentration of 4.0 g/100 mL). The system was stirred under 300 rpm and kept at 20 °C. After the system reached a steady state according to FBRM data (after ∼40 min), both the inlet (only water, no CBZ ethanol solution) and outlet pumps were initiated simultaneously. The flow rate of water (antisolvent) in and the slurry out were set to the same target to achieve a residence time of 30 min (the crystallizer runs at a set volume of 400 mL and thus the flow rate is 13.3 mL/min). The process was monitored by FBRM as shown in Figure 8.

Figure 9. Dynamic profile obtained for characteristic peaks for CBZ form III and II during a continuous MSMPR operation. Blue (1050 cm−1, form III); green (1027 cm−1, form II); red (257 cm−1, form II); light pink (391 cm−1, form II). The Y-axis is peak height. Continuous operation and antisolvent addition were initiated after 40 min (red dashed line).

ratio near 1040 cm−1, between 1050 and 1027 cm−1 was clearly seen. When form II is formed in the system, it is expected that the peak at 1027 cm−1 would increase and 1050 cm−1 decrease, which is clearly observed in Figure 9; in addition, the characteristic peaks of form II at 257 and 391 cm−1 increased. However, it should be noted that, for polymorphic systems like CBZ which lacks major characteristic peak differences, selection of PAT tool with adequate sensitivity is crucial in CM process development. If needed, signal processing or multivariate analysis should be considered to aid detection and monitoring. The off-line analysis of product 1, product 2, and Cr 2 was conducted to further characterize the properties and completeness of the polymorph transition. From the XRPD patterns in Figure 10, we can see that both product 1 and product 2 and Cr 2 are mainly form II crystals as shown from the characteristic peaks at 2θ values of 8.7 and 13.3.28 There are some small peaks in the XRPD pattern of product 1 (one residence time) indicating the existence of some form III crystals. Figure 11 shows that the PSD of crystals (product 2) becomes broader after the antisolvent addition. It is possible that the antisolvent addition created supersaturation very fast to generate a large amount of nuclei or small crystals and therefore the PSD tends to be broader, or because the shape of form II is needle-like (CBZ crystals has been well-characterized in literature and it has been observed from microscopy that form II results in needle-like crystals.40,41 It should be noted that both FBRM and Raman methods will be developed and validated in later studies to better monitor the process dynamics with full details. 3.2. Drawing Scheme for Intermediate Crystal Collection. In continuous crystallization, crystals should be collected during the process, without disturbing the state of control (steady state in some cases). Therefore, the design for collecting the crystals in continuous crystallization should be carefully considered. For example, some questions should be asked: (1) where in the crystallizers should the crystals be collected? (2) Should the drawing be continuous or periodic? (3) How should the liquid level in the crystallizers be controlled? Regarding the collecting methodologies, there are several common setups.12,14,36 One example is continuous skimming of

Figure 8. FBRM trends of counts in various ranges during adding antisolvent to the system.

The experiment included two residence times to ensure full transformation of the polymorphs in the vessel. Samples from the product crystals after one residence time, product crystals after two residence times and crystals in the crystallizer itself after two residence times were collected, filtered, and denoted as product 1, product 2, and Cr 2. The Raman spectra were also monitored during the process, three of which are shown in Supplement. After one residence time (∼01:20:00 in Figure 9), some obvious changes can be observed for the characteristic peaks. After two residence times (∼1:50:00 in Figure 9), a change of peak area 1026

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

Figure 10. XRPD results of product 1, product 2, and Cr 2 comparing with pure form II and form III samples.

suspension, due to the settling of larger crystals. A single-stage continuous crystallization system was used to test different drawing schemes from the top (see Supporting Information). For each test, crystallization was run for approximately one residence time (∼30 min), so the collection tank has the same amount of suspension as the crystallizer tank. The collection efficiency was defined by percentage as (weight of crystals in collection tank)/(weight of crystals in crystallizer tank); results are shown in Table 2. By visual inspection, the suspension in the Table 2. Collection Efficiency by Skimming at Different Stirring Speeds RPM

Figure 11. Comparison of the particle size distribution of sampling drawn before (red) and after adding antisolvent (blue, product 2).

300 500

the top layer of the suspension, as shown in Figure 12 (left). This setup requires that the pump handle both slurry and air and that level control is achieved by the position of the tubing end. However, in terms of the number and particle size of crystals, material from the top layer may not be representative of the bulk

500 500

collection efficiency by continuous skimming 22% 19% by periodic drawing 49% by bottom drawing 73%

crystallizer tank appeared uniform, indicating that settling of large crystals may not be an issue during the experiment. To further investigate settling, the stirring speed was increased from 300 to 500 in the second experiment, but no improvement on collection efficiency was observed. The uniformity of the suspension was further studied using FBRM as qualitative tool and will be addressed in later sections (section 3.1.3). The result suggests that settling of large crystals was not the cause here for low collection efficiency. Periodic drawing is another widely used setup as shown in Figure 12 (right).12,42,43 In a periodic drawing scheme (or alternatively intermittent, cyclic, or pulsatile), the liquid level is allowed to rise above that of the draw tube in the skimming

Figure 12. Schematic for skimming (left) and periodic collection (right) experiments. 1027

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

Figure 13. (a) Schematic of bottom draw by metal guide and (b) particle size distribution in collection tank (red square) and crystallizer (blue circle).

scheme, with normal inlet flow rate with pneumatic system disengaged. After a set amount of time or volume, the pneumatic system is engaged by reintroduction of seal to the tank. In a relative short amount of time, the liquid is transferred via the pressure gradient and reduced back to the level of the draw tube, at which point the seal on the tank is disengaged. The additional 15% of the tank volume or ∼13% of the total overflow volume is similar to those cited in literature for lab-scale continuous crystallization system.44 At 500 rpm, the periodic drawing (49%) is able to collect more than twice the amount of crystal solid as that by skimming. This could be due to the draw tube being more submerged in the tank at the onset of the drawing by the pneumatic system (location), or the longer duration it has in collecting the suspension (dynamics). However, it was observed that the mean particle size in the collection tank was slightly smaller than that in the crystallizer tank (see the Supporting Information). In addition, during nonflow periods in periodic drawing, the liquid residue evaporates to crystallize solid CBZ at the tubing end. After accumulating over a sufficiently long period, the solid may clog the flow or fall out as product, affecting the process robustness or product quality in the tank. The third configuration evaluated was the “bottom draw” scheme (Figure 13a), where a metal guide tube is inserted along the side to the bottom of the tank. The experiment performed with continuous drawing led to a collection efficiency of 73% and almost identical particle size distributions (PSD) in the two tanks as shown in Figure 13b. This represented a significant improvement from skimming and periodic draws by capturing the majority of solid crystals and removing inefficiencies of drawing over an inconsistent flow. In addition, the efficiency was improved by the reduction of clogging of the outlet tube (addressed later in this article). The “bottom draw” was selected for the final material collection design and integrated with the flow and level control. Further improvement could be performed in the design of the continuous system by using gravity in favor of the continuous flow. Bottom draw by gravity is a drawing scheme which can be continuous or periodic. Unlike the previous three methods, this setup normally does not require pumping or pressure difference to transfer slurry. Instead, the slurry naturally flows into the second crystallizer by gravity without energy consumption. This setup was not studied since it requires customized glassware, which is not congruent with the existing experimental apparatus. An interesting question to consider is whether a representative (i.e., entire PSD) collection is strictly necessary to mitigate product quality risk. If the material balance is achieved and the PSD of crystals reaches a dynamic equilibrium within the

crystallizers, the crystallization may be in a state of control and run continuously. An example of this may be a drawing scheme that consistently pulls larger crystals but allows nucleation and small crystal growth to occur simultaneously. In this case, the process yield and efficiency will be maximized. However, achieving this may require complex design to only extract large crystals and comprehensive understanding toward crystallization kinetics (or even process modeling/simulation tools). Also, any potential disruptions in flow that would allow a different population of crystal particle size to exit the vessel would need to be avoided. Nonetheless, the main attribute of an efficient collect scheme is the consistent performance throughout the continuous operation with respect to particle properties of interest. 3.3. Slurry Transport. After drawing the crystals from the vessels, the next question here is how to handle the slurry transport to another stage or filtration? Slurry transport may encounter several issues. Clogging and settling of solids are major concerns since they may fully stop flow in the tubing, force the process variables out of design space, and require the whole system to shut down. Furthermore, the heat loss to the environment can cause either additional crystal growth during transfer and therefore clogging/settling in the tubing, or crystals to dissolve and thus lower yields for the process. In this section, we will discuss possible solutions to these risks. Maintaining consistent flow of the suspension and preventing clogging is a major challenge to the robust operation of continuous crystallization. The crystal sizes obtained in these experiments are mostly 10−300 μm, while the tubing used is 2.49 mm ID. However, the connection fittings reduce the inner diameter to approximate 1 mm. Fouling occurs generally at these connection points, causing a buildup of solid. Once particulate flow is obstructed, a solid plug forms, further restricting flow and creating of air bubbles downstream. To mitigate this risk, possible solutions are to remove unnecessary fittings, enlarge tubing inner diameters, increase flow velocity, selecting material of construction of tubing inner walls, etc. However, larger tubing inner diameters will decrease the flow velocity and cause settling (another failure mode) of particles if the total flow rate is kept constant. The system was improved via process design, by removing unnecessary connections and reducing clogging risk without increasing settling risk. After the system upgrades, in a 2h-run continuous crystallization experiment, only one clogging event was observed which was easily addressed with manual tools within 2 min. Overall, clogging can be a serious problem for process safety and product quality. Prevention via process design and control of process parameters and step-by-step instructions 1028

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

⎡ ⎤ ⎛ 1310 ⎞ v horizontal = 1.5⎢2 × 9.8⎜ − 1⎟ × 2.49 × 10−3⎥ ⎝ 789 ⎠ ⎣ ⎦

on handling clogging events should be a key consideration in the control strategy.45 Another potentially process risk factor is settling of crystals. Settling can be seen clearly through the transparent tubing via crystal solids flow at the bottom of the tubing (Figure 14),

⎛ 100 × 10−6 ⎞1/6 ×⎜ ⎟ = 0.0283 m/s ⎝ 2.49 × 10−3 ⎠

In addition, the minimum transport velocity (i.e., for heterogeneous suspensions) in vertical tubing is calculated by the following equation:47,48 vvertical =

Re p =

g (s − 1)d p2

for Re p < 1.4

18ν

d pvvertical ν

where Rep is the particle Reynolds number, and ν is the fluid kinematic viscosity (= viscosity/density). A value of 0.001095 Pa· s was used for the viscosity of ethanol in the calculation. So the minimum vertical transport velocity is

(

9.8 vvertical =

1310 789

)

− 1 (100 × 10−6)2

(

18

0.001095 789

)

= 0.00259 m/s

Re p = 0.187 < 1.4 The actual velocity should be designed to be larger than both minimum transport velocities νhorizontal and νvertical. In our case, the actual velocity in the tubing is vactual =

Figure 14. Solid stream visible at the bottom of the tubing during transfer.

13.3 mL/min 13.3 × 10−6 /60 = = 0.0452 m/s 2 2 πr 2.49 × 10−3 π 2

(

)

Finally, process risk due to slurry transport concerns heat loss to the environment. While the slurry temperature is wellcontrolled in the crystallizers, there may be significant cooling during transport between thermostat vessels. This can exacerbate the issues of clogging and settling, if a stream cools sufficiently to promote growth on the tubing inner wall. The best solution is to keep the transport distances as short as possible. In addition, insulation in form of a fiber glass wrap was implemented for the transfer line between the feed tank and CR 1, as this was the transition with the greatest temperature gradient (e.g., feed solution T = 50 °C). Due to loss of visual access to the connections, other methods of detecting clogging are also suggested to be implemented (e.g., periodic checks and verification at the end of the run). 3.4. Process Variables. Potential risks from process variables, such as temperature and flow rate, should be considered. Robust control of process variables is critical to avoid process induced disturbances, which may result in significant variations in the product quality. For example, an increase on the temperature can affect the solute concentration and crystal size distribution by dissolving fine crystals in the system.49 Therefore, it is important to choose an appropriate temperature control strategy. Commercially available external temperature baths with programmed simple loops to control the process temperature can be connected to the crystallizer with jackets, which allows for efficient controllability. In our work, the temperature and the agitation speed were controlled using the programmed control algorithm built within the iControl platform. This software allows for implementation

creating an observable gradient. Several potential risks to quality due to crystal settling can be (1) risk to material uniformity, (2) risk of degradation of material adhered to the walls, (3) perturbations in the transfer lines can initiate a bolus of slurry concentrated in crystal solids, and a potential for clogs downstream. The nature of the flow is a function of factors such as velocity, solid loading, particle size, and physical properties of the system.46 These physical phenomena should be sufficiently understood in order to avoid settling problems. Fluid dynamics modeling and heuristic equations can be used for the design and optimization of slurry transport.47,48 The minimum transport velocity (i.e., for heterogeneous suspensions) in horizontal tubing can be calculated by the modified Durand equation:47,48 ⎛ d p ⎞1/6 v horizontal = F[2g (s − 1)D]⎜ ⎟ ⎝D⎠

where F is an empirical constant that varies between 0.4 and 1.5 (the maximum value of 1.5 was used in the calculation), and D is the tubing diameter, s is the ratio of particle and fluid densities (1310 kg/m3 for CBZ and 789 kg/m3 for ethanol), dp is the particle diameter (100 μm was used in the calculation), and g is the gravity constant 9.8 m/s2. Therefore, 1029

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

of cooling and stir speed ramps in the EasyMax platform (refer to section 2.2). The minimum temperature suggested by the manufacturer for the configuration developed in this system is −10 °C. The temperature controller in the EasyMax platform maintained the temperature within ±0.2 °C at the operating conditions studied. Other cooling mediums such as cryostat or silicone oil should be considered if the desired operating temperature is out of the suggested range. Variations in temperature may affect the nucleation and crystal growth kinetics and thus alter the critical quality attributes of the product crystals during process. For our study, the effect of introduced disturbances in temperature on the product quality will be evaluated in later studies. Nonetheless, responsive control of process variables, such as temperature, will enable faster response in the case of disturbances that may occur during the continuous crystallization process. The agitation speed is also a potential risk factor since it impacts the mixing conditions of the process; variations in the agitation speed can affect the crystal growth and withdrawal of crystals, among others. The isokinetic withdrawal of the suspension from the vessel is critical for achieving consistent crystal properties. A homogeneous suspension should be achieved to avoid any settling of crystals in the system. The settling of crystals in the crystallizer may result in variations on the crystal growth since fine crystals have a higher probability to growth. The deposition of crystals at the bottom could also cause aggregation or agglomeration of particles. Hence, it is important to determine an agitation speed region in which internal particle classification effects in the operating volume are minimized. A series of batch cooling experiments were performed in which a saturated solution was cooled down from 40 to 5 °C at 0.5 °C per minute and 300 rpm. For a minimum of 2 h holding period or until equilibrium was reached, the position and agitation speed was varied following the methodology described by Hou et al.36 The agitation speed was varied among 150, 300, and 450 rpm. The volume was fixed to 400 mL since it is the operating volume used in the continuous crystallization process. Figure 15 shows the average CLDs obtained at both the top and the bottom of the crystallizer. The particle counts change significantly throughout the entire chord length between the top and the bottom of the crystallizer when the agitation speed was set to 150 rpm The results indicate that significant settling of particles occurred at this condition. However, little differences were observed between the top and the bottom CLDs for the scenarios in which the agitation speed was set to 300 and 450 rpm, as observed in Figure 15b−c. The agitation speeds between 300 and 450 rpm demonstrated sufficient mixing to avoid differences on the CLDs. Although the approximated power input per unit volume for 300 to 450 rpm is 0.05−0.2 kW per m3, the Reynolds number for stirred vessels was approximately 3 × 103 for 300 rpm and 104 for 450 rpm, indicating that the flow patterns were both in the turbulent region.45 Therefore, agitation speed is a potential risk factor that should be well-studied and controlled, since it can affect product quality (e.g., particle size) and process efficiency. Another potential risk factor is the working volume (or the suspension level). The suspension level can change due to variations on the flow rate due to clogging, pump malfunction, or flow of air, among others. Hence, an efficient control strategy that could react to changes on the suspension level needs to be implemented to avoid variations that will affect the product properties. The development requires the design and implementation of a level indicator, efficient tuning of the pumps, and a

Figure 15. Comparison between the CLD of suspensions withdrawn from the top or bottom of the crystallizer at different agitation speeds: (a) 150, (b) 300, and (c) 450 rpm.

software for automation. Hydrostatic, radar, and ultrasonic level indicators, among others, can be implemented for continuous monitoring. The indicator should be capable of providing accurate monitoring for suspensions of high solid concentration. For this work, a bubbler was used as level indicator for the controller (refer to section 2.3.2). Tuning of the pumps and level control is necessary for the development of an efficient control framework. For our work, a level controller (feedback and feedforward controls) was developed and implemented as described in section 2 and Supporting Information. The performance of controller was improved by tuning the PID control loops via open-loop studies. The characteristic gain, time constant, and dead time can be evaluated from the dynamic response as a result of the step increase. A set of tuning rules, such as Ziegler-Nichols, can be used to estimate appropriate tuning parameters for the PID controllers. As Figure 16 illustrates, the controller is satisfactory in keeping the level around 105 mL (±5 mL) with setting the P2 at the set point (7.5 mL/min) and varying the P1 in real time through feedback control. Further fine-tuning of controller parameters may be able to further reduce the level oscillations. Overall, process variables such as agitation speed, suspension level, and vessel temperature can have significant impact in product quality, often with a multivariate relationship. Therefore, 1030

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

Figure 16. Closed loop response of level controller around its set point. The pump flow rate is shown in the top figure, and the level measurement and set point (105 mL) are shown in the bottom figure.

implementing a control strategy to monitor and reduce variability in process parameters will be the key to achieving consistent quality of the crystals and process efficiency. 3.5. Case Study: Robustness of Continuous Platform. To evaluate the robustness of the system, long-run experiments of continuous crystallization were performed on the two-stage MSMPR system. The results of two long-run experiments are discussed here. Specifically, we picked a scenario without clogging events (Experiment 1) and another one with a clogging event (Experiment 2) to demonstrate the monitoring ability of PAT tools. It should be noted that we continue to optimize our system as discussed above. Currently we have seen relatively low occurrence of clogging in the system. Since only one FBRM and one Raman mobile station were available, the FBRM was placed into the first crystallizer for particle size distribution, and the Raman probe was put into the second crystallizer for polymorphism monitoring. The first crystallizer is mainly for primary/secondary nucleation, and the second crystallizer is mainly for crystal growth, so monitoring the particle size distribution in Cr 1 can better capture process dynamics and disturbances. Due to the proximity to final product collection, monitoring the polymorphs in Cr 2 was selected as a good preliminary surrogate for monitoring product polymorphic forms. Polymorphism change during transport in the tubing from Cr 2 to filtration (∼10 s) is low likelihood. The system was precharged with 4.0 g/100 mL of CBZ ethanol solution in both crystallizers. The P2 flow rate from Cr 1 (400 mL) to Cr 2 (400 mL) is 13.3 mL/min for a residence time of 30 min each stage (in total, 2−3 h was planned for an initial experiment, and normally the system needs 2−4 residence times to reach a relatively steady state). The P1 and P3 flow rate were controlled by the automation system as described above to maintain the liquid levels of Cr 1 and Cr 2. Both crystallizers were stirred at 300 rpm. The Cr 1 temperature was set at 15 °C, and Cr 2 temperature was set at 5 °C. To initiate nucleation, form III crystals (in-house made) were seeded. FBRM and Raman data on two continuous run experiments with the same conditions were collected). Due to the limited volume of the feed solution, the longest run time achieved so far is 2 h and 17 min for Experiment 2 (Figure 17). However, we anticipate it can be run continuously for a much longer period. The FBRM data can be used for qualitative analysis of the dynamic behavior for the continuous process. In Experiment 1, the FBRM data (Figure 17a) seems to indicate that the system

Figure 17. FBRM data showing dynamic profile and relatively steady (RS) operation for case studies where (a) no clogging and (b) one clogging event occurred.

achieved a relatively steady operation at about ∼20 min after seeding. However, in Experiment 2, there was a clogging in tubing between Cr 1 and Cr 2, for about 10 min (around 45 min) as shown in Figure 17b). The system was able to detect the clogging event, which caused significant crystal number growth in Cr 1 as the FBRM trends show. After the clogging was detected and resolved, the numbers of crystals started to decrease gradually. It took approximately two residence times (∼80 min) for the system to return to a relatively stable operation after the 10 min clog event. However, the square weighted mean crystal size obtained for Experiments 1 and 2 is close to 200 μm which shows that the clogging event has a small impact on the crystal size after it reaches a relatively stable state, as characterized by process conditions (temperature, stirring rate, residence time, etc.). The Raman spectra (Figure 18) in both experiments suggest form III consistently, and the offline XRPD analysis also confirmed all products were form III. In addition, the large variations in Raman spectra of Experiment 2 also suggest the clogging event and its development as discussed above. In both experiments, during a fixed input volume (200 mL) of feed solution, yield was calculated as the ratio between the weight of the crystals collected and the amount of CBZ in the 200 mL feed solution supposed to crystallize at 5 °C (CBZ solubility at 5 °C = 16.1 ± 0.2 g/kg solvent): yield = crystals collected during input of 200 mL feed solution amount of CBZ in 200 mL feed solution supposed to crystallize at 5 °C

Experiment 1’s yield is 44.1%, and Experiment 2’s yield is 45.0%. Given the short residence time (totally 60 min), the results seem reasonable compared to multiple-stage continuous crystallization process without recycle.35 For the purpose of a 1031

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

times to study system dynamics and reproducibility of process startup. Issues with clogging and particle settling have been mostly addressed and are detectable by the PAT. Next steps include process variation studies, definition of system deliverables for product quality (form, particle size, morphology), method development, and validation for PAT, in order to evaluate more advanced control strategies, for example, an automatic system for statistically detecting clogging issues, robustness of level control, and modeling/simulation applications to predict system behaviors.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.oprd.7b00130. Detailed description on the tuning and development of the control system. The calibration procedure of the bubbler is explained in detail. Also, figures showing more details of the experimental setup are provided (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Figure 18. Real-time Raman spectra obtained for studies where (a) no clogging and (b) one clogging event occurred.

Xiaochuan Yang: 0000-0001-7003-444X Notes

The authors declare no competing financial interest.

benchtop two-stage continuous crystallizer to be used for demonstrating control strategy concepts, ∼45% yield is considered adequate. However, continuous crystallization in industry may have additional stages and/or recycle streams to increase yield. As the number of stages increases, the residence time will be increased, and thus the total yield will be increased. As a result, the solid density in the slurry will be higher at the last few stages, and the material transport would need to address more challenges to avoid potential clogs. This case study illustrated that, by mitigating the risks through the above considerations, we were able to run the continuous crystallization for a period of time. Our LabVIEW-based feedforward and feedback control system was able to maintain a relatively steady operation for the crystallizers. In addition, the PAT tools (Raman and FBRM) successfully detected when there was a clogging event in the system.



ACKNOWLEDGMENTS We appreciate the collaboration and funding from FDA-MIT David Koch Practice School Internship Program, and the help on instrumentation from Mettler Toledo. This publication only reflects the views of the authors and should not be construed to represent FDA’s views or policies.



REFERENCES

(1) Mascia, S.; Heider, P. L.; Zhang, H.; Lakerveld, R.; Benyahia, B.; Barton, P. I.; Braatz, R. D.; Cooney, C. L.; Evans, J. M. B.; Jamison, T. F.; Jensen, K. V.; Myerson, A. S.; Trout, B. L. Angew. Chem., Int. Ed. 2013, 52, 12359−12363. (2) Lee, S. L.; O’Connor, T. F.; Yang, X.; Cruz, C. N.; Chatterjee, S.; Madurawe, R. D.; Moore, C. M. V.; Yu, L. X.; Woodcock, J. J. Pharm. Innov. 2015, 10, 191−199. (3) Byrn, S.; Futran, M.; Thomas, H.; Jayjock, E.; Maron, N.; Meyer, R. F.; Myerson, A. S.; Thien, M. P.; Trout, B. L. J. Pharm. Sci. 2015, 104 (3), 792−802. (4) Vertex Case Study. Available from: http://www.gea.com/global/ en/stories/vertex-case-study.jsp (cited April 4, 2016). (5) Smith, D. E. Quality by Design (QbD) for the Continuous Manufacturing of Solid Oral Dosage Forms. Available from: http://pqri. org/wp-content/uploads/2015/11/Embiata-Smith.pdf (April 4, 2016). (6) FDA Approves Tablet Production on Janssen Continuous Manufacturing Line. Available from: http://www.pharmtech.com/fdaapproves-tablet-production-janssen-continuous-manufacturing-line (April 13, 2016). (7) Pimparkar, K.; Yen, B.; Goodell, J. R.; Martin, V. I.; Lee, W. H.; Porco, J. A.; Beeler, A. B.; Jensen, K. F. J. Flow Chem. 2011, 1, 53−55. (8) Goodell, J. R.; McMullen, J. P.; Zaborenko, N.; Maloney, J. R.; Ho, C.-X.; Jensen, K. F.; Porco, J. A., Jr.; Beeler, A. B. J. Org. Chem. 2009, 74, 6169−6180. (9) Heider, P. L.; Born, S. C.; Basak, S.; Benyahia, B.; Lakerveld, R.; Zhang, H.; Hogan, R.; Buchbinder, L.; Wolfe, A.; Mascia, S.; Evans, J. M. B.; Jamison, T. F.; Jensen, K. F. Org. Process Res. Dev. 2014, 18, 402−409.

4. CONCLUSION Continuous manufacturing, as an emerging technology in the pharmaceutical industry, presents many challenges for process and product development. It is essential to understand and consider the risks and how they can be mitigated via process dynamics understanding and process analytics tools capable of detecting changes at a relevant time scale. In this work, a proofof-concept two-stage MSMPR continuous crystallization system was developed with integration of two PAT tools (Raman and FBRM) to monitor the process in real time. This research scale system can be used to evaluate concepts of control strategy development and process risks. Currently, the system is successfully demonstrated to perform two-stage cooling crystallization of CBZ. Preliminary studies indicate that a level control system with bottom-draw suspension removal configuration and alternating filtration setup can provide a basis of operation. However, the current setup requires further optimization, particularly around defining characteristic run 1032

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033

Organic Process Research & Development

Article

(10) Wong, S. Y.; Tatusko, A. P.; Trout, B. L.; Myerson, A. S. Cryst. Growth Des. 2012, 12 (11), 5701−5707. (11) Quon, J. L.; Zhang, H.; Alvarez, A.; Evans, J.; Myerson, A. S.; Trout, B. L. Cryst. Growth Des. 2012, 12 (6), 3036−3044. (12) Peña, R.; Nagy, Z. K. Cryst. Growth Des. 2015, 15 (9), 4225−4236. (13) Zhang, H.; Lakerveld, R.; Heider, P. L.; Tao, M.; Su, M.; Testa, C. J.; D’Antonio, A. N.; Barton, P. L.; Braatz, R. D.; Trout, B. L.; Myerson, A. S.; Jensen, K. F.; Evans, J. M. B. Cryst. Growth Des. 2014, 14 (5), 2148−2157. (14) Yang, Y.; Song, L.; Nagy, Z. K. Cryst. Growth Des. 2015, 15 (12), 5839−5848. (15) Wong, S. Y.; Cui, Y.; Myerson, A. S. Cryst. Growth Des. 2013, 13 (6), 2514−2521. (16) Alvarez, A. J.; Myerson, A. S. Cryst. Growth Des. 2010, 10 (5), 2219−2228. (17) Lawton, S.; Steele, G.; Shering, P. Org. Process Res. Dev. 2009, 13 (6), 1357−1363. (18) Zhang, S.; Ferte, N.; Candoni, N.; Veesler, S. Org. Process Res. Dev. 2015, 19 (12), 1837−1841. (19) Yang, H.; Yu, X.; Raval, V.; Makkawi, Y.; Florence, A. Cryst. Growth Des. 2016, 16, 875−886. (20) Majumder, A.; Nagy, Z. K. Cryst. Growth Des. 2015, 15 (3), 1129− 1140. (21) Liu, W. J.; Ma, C. Y.; Wang, X. Z. Procedia Eng. 2015, 102, 499− 507. (22) Yazdanpanah, N.; Ferguson, S. T.; Myerson, A. S.; Trout, B. L. Cryst. Growth Des. 2016, 16 (1), 285−296. (23) Wang, M.; Rutledge, G. C.; Myerson, A. S.; Trout, B. L. J. Pharm. Sci. 2012, 101 (3), 1178−1188. (24) Yang, X.; Myerson, A. S. CrystEngComm 2015, 17 (4), 723−728. (25) Yang, X.; Ong, T. C.; Michaelis, V. K.; Heng, S.; Griffin, R. G.; Myerson, A. S. CrystEngComm 2015, 17 (31), 6044−6052. (26) Acevedo, D.; Ling, J.; Chadwick, K.; Nagy, Z. K. Cryst. Growth Des. 2016, 16 (8), 4263−4271. (27) Yang, X.; Ong, T. C.; Michaelis, V. K.; Heng, S.; Huang, J.; Griffin, R. G.; Myerson, A. S. CrystEngComm 2014, 16 (39), 9345−9352. (28) Grzesiak, A. L.; Lang, M.; Kim, K.; Matzger, A. J. Pharm. Sci. 2003, 92 (11), 2260−2271. (29) Gomes, V. G.; Yee, K. W. K. Sep. Purif. Technol. 2002, 28, 161− 171. (30) Preikschat, F. K.; Preikschat, E. United States Patent 4871251, 1989. (31) Kail, N.; Marquardt, W.; Briesen, H. Chem. Eng. Sci. 2009, 64, 984−1000. (32) Leyssens, T.; Baudry, C.; Escudero-Hernandez, M. L. Org. Process Res. Dev. 2011, 15, 413−426. (33) Whelan, J.; Murphy, E.; Pearson, A.; Jeffers, P.; Kieran, P.; McDonnell, S.; Glennon, B. Bioprocess Biosyst. Eng. 2012, 35 (6), 963− 975. (34) Barthe, S. C.; Grover, M. A.; Rousseau, R. W. Cryst. Growth Des. 2008, 8 (9), 3316−3322. (35) Yu, Z. Q.; Chow, P. S.; Tan, R. B. H. Org. Process Res. Dev. 2008, 12, 646−654. (36) Hou, G.; Power, G.; Barrett, M.; Glennon, B.; Morris, G.; Zhao, Y. Cryst. Growth Des. 2014, 14, 1782−1793. (37) Hickey, M. B.; Peterson, M. L.; Scoppettuolo, L. A.; Morrisette, S. L.; Vetter, A.; Guzman, H.; Remenar, J. F.; Zhang, Z.; Tawa, M. D.; Haley, S. Eur. J. Pharm. Biopharm. 2007, 67 (1), 112−119. (38) Tian, F.; Zeitler, J. A.; Strachan, C. J.; Saville, D. J.; Gordon, K. C.; Rades, T. J. Pharm. Biomed. Anal. 2006, 40 (2), 271−280. (39) O’Brien, L. E.; Timmins, P.; Williams, A. C.; York, P. J. Pharm. Biomed. Anal. 2004, 36, 335−340. (40) Cruz Cabeza, A. J.; Day, G. M.; Motherwell, W. D. S.; Jones, W. Chem. Commun. 2007, 16, 1600−16002. (41) Rustichelli, C.; Gamberini, G.; Ferioli, V.; Gamberini, M. C.; Ficarra, R.; Tommasini, S. J. Pharm. Biomed. Anal. 2000, 23, 41−54. (42) Polster, C. S.; Cole, K. P.; Burcham, C. L.; Campbell, B. M.; Frederick, A. K.; Hansen, M. M.; Harding, M.; Heller, M. R.; Miller, M.

T.; Phillips, J. L.; Pollock, P. M.; Zaborenko, N. Org. Process Res. Dev. 2014, 18, 1295−1309. (43) Cole, K. P.; Campbell, B. M.; Forst, M. B.; McClary Groh, J.; Hess, M.; Johnson, M. D.; Miller, R. D.; MItchell, D.; Polster, C. S.; Reizman, B. J.; Rosemeyer, M. Org. Process Res. Dev. 2016, 20, 820−830. (44) Ferguson, S.; Ortner, F.; Quon, J.; Peeva, L.; Livingston, A.; Trout, B. L.; Myerson, A. S. Cryst. Growth Des. 2014, 14 (2), 617−627. (45) Myerson, A. Handbook of Industrial Crystallization; ButterworthHeinemann, 2002. (46) Weber, M. Liquid-solid flow. Thermopedia, A-to-Z Guide to Thermodynamics, Heat & Mass Transfer, and Fluids Engineering, 2010. Available from: http://thermopedia.com/. (47) Heywood, N. I. Chem. Eng. Prog. 1999, 95 (9), 21−41. (48) Shook, C. A. Slurry Pipeline Flow. Processing of Solid-Liquid Suspensions; Butterworth-Heinemann: Boston, 1993. (49) Nagy, Z. K.; Aamir, E.; Rielly, C. D. Cryst. Growth Des. 2011, 11, 2205−2219.

1033

DOI: 10.1021/acs.oprd.7b00130 Org. Process Res. Dev. 2017, 21, 1021−1033