Chiral Symmetry Breaking and Deracemization of Sodium Chlorate in

Dec 12, 2017 - The chiral symmetry breaking and deracemization of sodium chlorate was investigated when varying the agitation speed and cooling rate i...
1 downloads 3 Views 547KB Size
Subscriber access provided by University of Florida | Smathers Libraries

Article

Chiral Symmetry Breaking and Deracemization of Sodium Chlorate in Turbulent Flow JaekyuKyu Ahn, Do Hyun Kim, Gérard Coquerel, and Woo-Sik Kim Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.7b01247 • Publication Date (Web): 12 Dec 2017 Downloaded from http://pubs.acs.org on December 21, 2017

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 free 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 accessible to all readers and 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 32 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

Chiral Symmetry Breaking and Deracemization of Sodium Chlorate in Turbulent Flow

Jaekyu Ahn1, Do Hyun Kim2, Gerard Coquerel3 and Woo-Sik Kim1*

1

Functional Crystallization Center, Department of Chemical Engineering, Kyung Hee University, Kyungki-do, Korea

2

3

Department of Chemical Engineering, KAIST, Daejeon, Korea SMS laboratory EA3233, University of Rouen Normandy, F-76821 Mont Saint Aignan CEDEX, France

* Corresponding author: Woo-Sik Kim ([email protected])

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

Abstract The chiral symmetry breaking and deracemization of sodium chlorate was investigated when varying the agitation speed and cooling rate in cooling crystallization. At a low agitation speed, almost zero crystal enantiomeric excess occurred at the induction period. When increasing the agitation speed, the crystal enantiomeric excess at the induction period (called the ‘initial CEE’) increased due to the promotion of chiral symmetry breaking. The chiral symmetry breaking at the induction period (called the ‘initial chiral symmetry breaking’) also varied with the cooling rate. At a low cooling rate of 0.0738oC/min, the initial CEE reached up to about 90% and was rapidly reduced when increasing the cooling rate. The experiments also showed enhanced deracemization of the chiral crystals when increasing the initial CEE. Thus, complete deracemization was achieved when the initial CEE was over 60%. The influence of the agitation speed and cooling rate on the initial CEE originated from secondary nucleation depending on the supersaturation at the induction period (called the ‘induction supersaturation’). Using a nucleation rate equation, the initial CEE was found to correlate well with the induction supersaturation. Also, varying the final setting temperature and agitator confirmed that secondary nucleation was significantly involved in the chiral symmetry breaking at the induction period.

Key words: chiral symmetry breaking; deracemization; induction supersaturation; Eve crystals; secondary nucleation

ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32 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. INTRODUCTION The isolation of chiral enantiomers is currently a special area of interest for research on fine chemicals and pharmaceuticals, as the biological and pharmacological actions of two enantiomers can differ distinctly despite identical non-vectorial physico-chemical properties, such as: density, solubility in non-chiral solvents, pKa, melting points, hardness, etc. Thus, many studies have investigated chiral symmetry breaking and deracemization to obtain homochiral crystals. Sodium chlorate is an interesting chiral substance, being optically inactive in a solution and generating optically active Lform and D-form crystals via crystallization. Since Kipping and Pope first demonstrated chiral symmetry breaking over 100 years ago by adding D-glucose to the crystallization of sodium chlorate,1 many studies have used sodium chlorate as a model substance for investigating the mechanism and effective ways for chiral symmetry breaking and deracemization. The milestone study on chiral symmetry breaking was reported by Kondepudi et al.2. While symmetric chiral crystals were produced with unstirred evaporation crystallization, complete chiral symmetry breaking was only obtained when stirring the crystallization, supposedly due to secondary nucleation by the stirring. They also used a stochastic kinetic model to explain the chiral symmetry breaking mechanism via secondary nucleation.3 Many subsequent studies have also examined chiral symmetry breaking. For example, McBride and Carter suggested the Eve effect (authors called it the ‘Adam effect’) on chiral symmetry breaking, where the initial ‘Eve’ crystal is formed and grows to a certain size, and then hundreds of offspring crystals are cloned via secondary nucleation in stirred crystallization.4 Plus, using the Metcalfe and Ottino model, Cartwright et al. developed an advection-mediated nonlinear autocatalytic process to describe chiral symmetry breaking by secondary nucleation in a fluid flow.5,6 Viedma also explained chiral symmetry breaking in terms of supersaturation in stirred crystallization.7 With a low supersaturation, the primary nucleation slowly generates a few homo-chiral crystals, which then become mother crystals breeding many daughter crystals via secondary nucleation, resulting in chiral symmetry breaking. However, with a high superstation, the primary nucleation generates stochastically a high number of two crystal forms with almost equal kinetics, resulting in a nearly symmetric amount of chiral crystals. Similar results were also reported by Szurgot and Szurgot.8

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

Meanwhile, Botsaris studied the influence of the seed and supersaturation on chiral symmetry breaking in crystallization.9,10 When using seeds, chiral symmetry breaking was always obtained via secondary nucleation, regardless of stirring the crystallization. However, the chiral symmetry breaking varied according to the seed size and supersaturation. The secondary nucleation was enhanced when increasing the seed size, thereby increasing the chiral symmetry breaking.9 Thus, the secondary nucleation of sodium chlorate was suggested to be proportional to L4 (L=seed size).11 Meanwhile, when increasing the supersaturation, spontaneous primary nucleation became more predominant in the seed crystallization, thereby diminishing the chiral symmetry breaking of the product crystals. In contrast, Qian and Botsaris observed complete chiral symmetry breaking with unseeded crystallization, even at a high supersaturation inducing spontaneous primary nucleation.9 Yet, while this result is the opposite of other studies, their study did not measure the chirality of the crystals generated at the induction period of nucleation in the unseeded crystallization. Moreover, the crystallization time was long at about 18 days and the chiral symmetry breaking was only analyzed at end of the crystallization. Thus, it is possible that the unseeded crystallization achieved complete chiral symmetry breaking via deracemization during the long crystallization, even though symmetric chiral crystals were initially formed at a high supersaturation. Most previous studies have found that seed crystals duplicate crystals of the same form via secondary nucleation. Notwithstanding, the formation of crystals dissimilar to the seed via secondary nucleation has been observed in stirred tank and oscillatory crystallizers due to the fluid shear.12 Niinomi et al. explained chiral symmetry breaking in unseeded crystallization at a high supersaturation using the phase transformation of achiral metastable crystals.13,14 At a high supersaturation, the primary nucleation quickly generates achiral metastable crystals, which are optically inactive. These metastable crystals are then slowly transformed into optically active chiral crystals due to the solubility difference between metastable crystals and chiral crystals. During this phase transformation, chiral symmetry breaking was induced via secondary nucleation. Other influences on chiral symmetry breaking in crystallization have also been investigated, such as ultrasonication, magnetic fields, and additives.15-17 While most of the above-mentioned studies of chiral symmetry breaking were carried out in a non-

ACS Paragon Plus Environment

Page 4 of 32

Page 5 of 32 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

equilibrium where nucleation plays an important role, a quasi-equilibrium where classical nucleation was not involved has been adopted for the enantiomeric enrichment of chiral crystals by deracemization. In the case of deracemization, the milestone study was reported by Viedma.18 In this study, a symmetric mixture of chiral crystals was completely shifted to a pure enantiomer using glass beads. The explanation provided for this enantiomeric enrichment by glass beads was attrition that produces tiny crystal fragments. These tiny fragments are then dissolved in the solution due to fluctuation of the equilibrium and consumed for the growth of large crystals, like Ostwald ripening.1922

They could also be directly incorporated on the surface of crystals of the same handedness. This

process was known as: Viedma ripening.20 Meanwhile, Uwaha and Katsuno suggested that chiral clusters were formed by grinding, which caused a significant initial chiral imbalance, resulting in enantiomeric enrichment.23

To promote deracemization, temperature cycling and ultra-sonication

have both been applied.24-27 According to Coquerel, deracemization was significantly accelerated by temperature cycling, as the dissolution of minor crystals and growth of major crystals were both promoted.24,25 Plus, ultra-sonication was more effective for deracemization than glass beads, as the mechanical impact of the cavitation by ultrasonic waves was stronger than the attrition by glass beads.26,27 Most previous studies have demonstrated that secondary nucleation was induced by the fluid flow in stirred crystallization, causing symmetry breaking. However, chiral symmetry breaking was invariably confirmed based on the chirality of crystals after a long crystallization of several days and weeks, and little attention has been paid to the nucleation in the early stage of crystallization, which was practically critical for determining chiral symmetry breaking. Moreover, the deracemization of chiral crystals during stirred crystallization has seldom been studied. Accordingly, the present study systematically investigated the influence of a turbulent flow on chiral symmetry breaking via secondary nucleation in the early stage of crystallization and the working mechanism of a turbulent flow on chiral symmetry breaking. In addition, the influence of a turbulent flow on the deracemization of chiral crystals during crystallization was also studied.

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 6 of 32

2. EXPERIMENT The sodium chlorate (NaClO3, >98.0%) was purchased from Duksan Chemical Co. (Korea) and used in the experiment without further purification. The feed solution was prepared by the dissolution of 96g of sodium chlorate in 100 g of water, which was saturated at 20 ºC, and then filtered using microfilter paper. A Rushton mixing tank was used as the crystallizer (called a mixing tank (MT) crystallizer). The MT crystallizer made of glass was equipped with an agitator and four-baffle for a turbulent eddy flow and thermal jacket on the outer wall for temperature control of the crystallization, as shown in Fig. 1. The working volume of the crystallizer was about 300 ml. Initially, the feed solution of sodium chlorate in the crystallizer was maintained at 26 ºC for 20 min. It was then cooled at a constant cooling rate to a final setting temperature of 13 ºC. Thereafter, the crystallizer was maintained at the final setting temperature for up to 100 hrs for further crystallization. During the crystallization, the agitation speed was varied from 300rpm to 1000rpm to adjust the hydrodynamic intensity of the turbulent flow, and the cooling rate was changed from 0.074 ºC/min to 0.59 ºC/min. For reproducibility, the crystallization was repeated at least three times under the same conditions. The induction point of nucleation was visually measured. Suspension samples were taken intermittently, quickly filtered using a vacuum pump, and dried in a desiccator. The chirality of the crystals and crystal sizes of both enantiomers (L- and D-forms of sodium chlorate) were analyzed using a polarized optical microscope (Olympus, BX53M, Japan). Over 500 crystals were counted from each sample to analyze the chirality and crystal sizes of both forms.

3. RESULTS and DISCUSSION The influence of the turbulent flow, generated by the propeller agitation in the mixing tank (MT) crystallizer, on the chiral symmetry breaking and deracemization of the sodium chlorate crystals was investigated when varying the agitation speed, as shown in Fig. 2. The chirality of the crystals was indexed based on the crystal enantiomeric excess (CEE), defined as

CEE =

.      .     .      .    

× 100(%)

ACS Paragon Plus Environment

(1)

Page 7 of 32 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

In eq. (1), the major crystals indicated the enantiomeric crystals in excess, while the minor crystals indicated the other enantiomeric crystals. With an agitation speed of 300rpm, the initial nucleation of the two enantiomeric crystals was almost symmetric. Thus, the CEE at the induction period (here, called the ‘initial CEE’) was close to zero. However, this crystal chirality slowly shifted to the major form by deracemization during the crystallization. Yet, no further deracemization occurred after reaching a CEE of 25%. When increasing the agitation speed, the initial CEE was also enhanced. Moreover, the deracemization of the chiral crystals was also promoted. As a result, complete deracemization was achieved with an agitation speed of 1000rpm. It should be mentioned that, in the present study, the induction point of nucleation was detected by a slight color change of the clear solution. However, at this point, the crystals were too tiny to clearly discriminate their chirality using a polarized optical microscope. Thus, the first sample was taken one minute after the induction point to analyze the crystal enantiomeric excess. This result was counted as the CEE at the induction period (initial CEE). The influence of the agitation on the chiral symmetry breaking at the induction period (called ‘initial chiral symmetry breaking’) was summarized in Fig. 3(a). The initial CEE rapidly increased from almost zero to 70% when increasing the agitation speed from 300rpm to 1000rpm. This initial chiral symmetry breaking was due to secondary nucleation by agitation, as suggested by Kondepudi and Viedma.2,3,7 Thus, if chiral L- and D-form crystals were spontaneously generated by a high primary nucleation rate with equal stochastics, two symmetrical chiral crystals population would be produced in the induction period. This symmetric formation of chiral crystals was usually observed with a high supersaturation above S=1.5.7 However, if the primary nucleation rate were slow with a low supersaturation, only a few enantiomeric crystals were initially generated. These enantiomeric crystals then became ‘Eve’ crystals duplicating offspring crystals with the same form via secondary nucleation, before crystals with the other form were nucleated and grown. Thus, the form of the Eve crystal population in the induction period dominated the crystal population over the other form, resulting in initial chiral symmetry breaking. Moreover, the secondary nucleation from the Eve crystals suppressed the crystal formation of the other form by reducing the supersaturation, also

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

contributing to the initial chiral symmetry breaking. According to Viedma, the primary nucleation rate was 2~3 homochiral nuclei/min under a low supersaturation. Thus, the interval between the crystal nucleation of one form and the other form was long enough for the induction of secondary nucleation from the Eve crystals, resulting in initial chiral symmetry breaking.7 In the present study, the promotion of secondary nucleation by agitation enhanced the initial chiral symmetry breaking, thereby increasing the initial CEE. This influence of the agitation speed on the initial symmetry breaking was also clearly reflected in the supersaturation at the induction period (called ‘induction supersaturation’), as shown in Fig. 3(b). Since the induction of nucleation was promoted by the turbulent flow, the primary nucleation was induced at a lower supersaturation when increasing the agitation speed and the primary nucleation rate was lower. Thus, the interval between the nucleation of the two forms increased, providing a higher chance for the secondary nucleation from the Eve crystals at the induction period. Therefore, the initial chiral symmetry breaking was enhanced when increasing the agitation speed. It should be mentioned that the induction supersaturation in the present study varied at around S=1.04 with the agitation speed, which was consistent with the supersaturation for chiral symmetry breaking by secondary nucleation reported by Viedma.7 In addition, it was interesting to note from Figs. 2 and 3(a) that the deracemization of the chiral crystals during the crystallization was related to the initial chiral symmetry breaking. Thus, when the initial chiral symmetry breaking was low with a CEE of 3%, the deracemization was also limited below a final CEE of 28%. When the initial chiral symmetry breaking was enhanced by the agitation, the limit of the deracemization in the crystallization was also increased. Thus, total deracemization of the chiral crystals was achieved when the initial CEE was 72% at 1000rpm. This experimental result suggested that the initial chiral symmetry breaking (initial CEE) was an important factor determining the deracemization of chiral crystals in crystallization and a certain initial level of chiral symmetry breaking is required for the complete deracemization of chiral crystals. The crystal size changes of the major and minor forms during deracemization were monitored, as shown in Fig. 4. At 300rpm, the major and minor crystals grew together during the deracemization. Thus, the crystal size difference between the major and minor forms was small. When

ACS Paragon Plus Environment

Page 8 of 32

Page 9 of 32 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

increasing the agitation speed, the major crystal size increased faster than the minor crystal size due to enantiomeric enrichment of the major form, resulting in a larger difference in the crystal sizes between the two forms at the end of the crystallization. On further increase of the agitation speed up to 1000rpm, the minor crystal size was rapidly reduced after reaching 65µm in the early stage, and eventually disappeared due to complete deracemization, whereas the major crystal size continuously increased. Interestingly, the crystal size profiles for the two forms revealed that for all the agitation conditions the initial major crystals were always larger than the initial minor crystals due to the Eve crystals. That is, the form of enantiomeric crystals generated by secondary nucleation from Eve crystals began to grow and consume the supersaturation earlier than the other form of enantiomeric crystals. Thus, the form of enantiomeric crystals nucleated first in the induction period were larger in size and higher in population than the other form of enantiomeric crystals. Furthermore, this high population and large size of the major crystals might have been a driving force for enantiomeric enrichment in the deracemization.

The deracemization of chiral crystals were also reflected in the crystal size distributions of two forms (Supporting Information Figure S1).

In this figure, it was

clearly observed that the crystal size distributions of two forms were shifted to large scale by the deracemization. However, the crystal population of major form increased whereas the crystal population was significantly reduced (Figure S1(a) and (b)). So, CEE at each class of crystal size was enhanced by the deracemization at 700 rpm of agitation speed (Figure S1(c)). When increasing the agitation speed to 1000 rpm, the crystal population of major form at the induction period was much higher than that of minor form at the induction period (Figure S1(d)), as estimated as high initial CEE in Figure 2. Due to complete deracemization, the no crystal of minor form remained at the end of crystallization, resulting in 100% of CEE at each class of crystal size (Figure S1(f)). This change of crystal size distribution was well consistent with change of CEE in the deracemization For a further examination of the influence of agitation on the deracemization of chiral

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

crystals, the agitation speed was changed during the crystallization, as shown in Fig 5. The crystallization was initiated with an agitation speed of 700rpm and cooling rate of 0.197oC/min. The agitation speed was then suddenly increased to 1200rpm after the CEE reached its limit under these crystallization conditions. As shown in Fig. 5(a), the initial chiral crystals were generated at about 35% of CEE and they were deracemized to about 67% of CEE at the end of crystallization (called the ‘final CEE’), which was similar to the deracemization at 700rpm in Fig. 2, even though the agitation was suddenly increased from 700rpm to 1200rpm during the deracemization. Furthermore, when comparing with the complete deracemization at 1000rpm in Fig 2, the final CEE in this deracemization was much lower. It was also interesting that the crystal sizes of both forms increased continuously, even though the deracemization did not proceed further (Fig. 5(b)). Here, the crystalsize-increase of the two forms after suddenly increasing the agitation speed from 700rpm to 1200rpm was slightly faster than that without a sudden increase of the agitation speed, as shown in Fig. 4(c), due to enhancement of the deracemization by a more rigorous turbulent flow. However, when the agitation speed was suddenly changed from 700rpm to 1200rpm at the induction period of nucleation, the deracemization was markedly different, as shown in Fig. 5(c). While a similar initial chiral symmetry breaking was obtained at 700rpm, the shift of chiral crystals to homo-chiral crystals was continuously amplified by the deracemization at 1200rpm and complete enantiomeric enrichment was achieved. Also, the crystal-size profiles of the two forms corresponded well to those in the complete deracemization, as shown in Fig. 5(d). From these results, it would seem that the agitation did not modify the deracemization if the crystals of the two forms grew above a certain size. Thus, the crystal size was identified as an influential factor determining deracemization. The turbulent flow was also involved in deracemization; a more rigorous turbulent flow attained a faster and higher deracemization. In Fig 6, the cooling rate was changed to investigate its influence on the initial chiral symmetry breaking and deracemization, as the induction supersaturation varied with the cooling rate. At a high cooling rate of 0.59 ºC/min, symmetrical chiral crystals appeared at the induction period and were partially deracemized up to about 40% of CEE after 100 hrs. When decreasing the cooling rate,

ACS Paragon Plus Environment

Page 10 of 32

Page 11 of 32 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

the initial chiral symmetry breaking was clearly amplified and reached an initial CEE to 90% at a cooling rate of 0.0738 ºC/min. Correspondingly, the deracemization was enhanced when decreasing the cooling rate. Thus, complete deracemization was achieved within 10 hours at a cooling rate below 0.15oC/min. These experimental results could be consistently explained in terms of the secondary nucleation depending on the induction supersaturation. When the induction supersaturation was reduced by decreasing the cooling rate, this produced a slow nucleation rate. As a result, the interval between nucleation of the two forms was enough for secondary nucleation from the Eve crystals, resulting in an increase of the initial chiral symmetry breaking. This description of the initial chiral symmetry breaking was also supported when the induction temperature and induction supersaturation were varied according to the cooling rate, as shown in Fig. 7. The induction supersaturation with a low cooling rate was low enough for secondary nucleation predominantly dictating chiral symmetry breaking at the induction period. When increasing the induction supersaturation, the crystal generation at the induction period was predominantly dictated by the primary nucleation, producing symmetrical chiral crystals, as suggested by Viedma (Fig. 7(a)).7 Therefore, the initial CEE was close to symmetrical chiral crystals with a high cooling rate of 0.59oC/min. The crystal-size profiles of the two forms were also monitored during the crystallization, as shown in Fig. 8. In case of complete deracemization with a low cooling rate -e.g. below 0.015 ºC/min- the minor crystals grew initially and then were completely dissolved out after reaching a maximum size, whereas the major crystals grew continuously. However, with a high cooling rate above 0.20 /min, both the major and minor crystals grew together. The effect of the induction supersaturation on chiral symmetry breaking was further examined when varying the final setting temperature, as shown in Fig. 9. In this case, different induction supersaturations were designed by different final setting temperatures without changing the crystallization conditions. Thus, at fixed crystallization conditions, including the agitation speed (700rpm), cooling rate (0.295 ºC/min), and initial temperature (26 ºC), three different final setting temperatures of 13oC, 16.3oC and 18.1oC were applied. When the final temperature was set at 13ºC, the induction of nucleation occurred at 14.7ºC. Therefore, 14.7 ºC was counted as the induction

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

supersaturation. However, when the final temperatures were set at 16.3ºC and 18.1ºC, the induction of nucleation occurred after the solution was cooled to the final setting temperature, respectively. Thus, the induction supersaturations were counted as 16.3ºC and 18.1ºC, respectively, for these two final setting temperatures. Here, it is interesting to note that the induction times for the three different final setting temperatures were similar at about 40 min, as shown in Fig. 9(a). At a high induction supersaturation of 14.7oC (final setting temperature of 13oC), the initial CEE was low as 10% due to the high primary nucleation at the induction period, as shown in Fig. 9(b). Meanwhile, when the induction supersaturation was reduced to 16.3oC and 18.1oC, the secondary nucleation was dominantly involved in the chiral symmetry breaking at the induction period increased, resulting in high initial CEEs. In this experiment, it was also confirmed that the initial chiral symmetry breaking was closely related to determining the deracemization. As such, when the initial CEE was low as 10% at an induction supersaturation of 14.7oC, the deracemization was limited to about 40 % of final CEE. The limit of deracemization increased when increasing the initial CEE. Thus, 100% deracemization was attained when the initial CEE was 67 % at an induction supersaturation of 18.1oC. The relationship of the induction supersaturation to the initial chiral symmetry breaking (initial CEE), and the initial chiral symmetry breaking (initial CEE) to deracemization (final CEE) are described in Fig. 10. Based on classical nucleation theory, Volmer’s equation was used for the primary nucleation rate expressed in terms of the supersaturation (S), and the interval between the crystal nucleation of the two forms was assumed as inversely proportional to the nucleation rate. Thus, the initial chiral symmetry breaking was correlated with the induction supersaturation on a plot of the initial CEE to 1/(lnS)2, as shown in Fig. 10(a). The plot revealed that as the induction supersaturation was reduced, the interval between the nucleation of the two forms (1/(lnS)2) increased, predicting a high initial CEE. Also, Fig. 10(b) suggested that the initial chiral symmetry breaking (initial CEE) dictated directly the deracemization of the chiral crystals (final CEE); a higher initial chiral symmetry breaking resulted in higher deracemization. An initial chiral symmetry breaking of at least 55% was required for complete deracemization.

ACS Paragon Plus Environment

Page 12 of 32

Page 13 of 32 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

Secondary nucleation was induced by collisions between mother crystals and between the mother crystals and the agitator in a turbulent flow. The collision efficiency depended on the pattern of the fluid motion and hardness of the agitator. When the fluid motion was vigorous and effective, the collision frequency of crystals and fluid shear would increase to induce high secondary nucleation. Also, harder agitator materials would be better for secondary nucleation by collision. Thus, three kinds of agitator; a steel propeller, steel impeller, and plastic propeller, were applied to the crystallization. Fig. 11 showed the influence of the agitator type and agitator material on the initial chiral symmetry breaking and deracemization. The steel agitators (steel propeller and steel impeller) were more effective for secondary nucleation by collision than the plastic propeller due to hardness of agitators. Thus, the initial chiral symmetry breaking by the steel agitators was much higher than that by the plastic propeller. Between the steel agitators (steel propeller and steel impeller), the steel propeller produced a higher initial chiral symmetry breaking than the steel impeller, as the strong radial and axial fluid motion provided by the propeller agitator was more effective for collision than the radial flow mainly generated by the impeller agitator. These experimental results also suggested that secondary nucleation was directly involved in the initial chiral symmetry breaking, and turbulent agitation was an influential parameter for chiral symmetry breaking and deracemization.

4. CONCLUSION In the cooling crystallization of sodium chlorate, the induction supersaturation, depending itself on the agitation and cooling rate, was identified as a key factor determining chiral symmetry breaking and deracemization. At a low induction supersaturation, secondary nucleation was significantly involved in the initial chiral symmetry breaking due to slow primary nucleation. However, when increasing the induction supersaturation, the contribution of secondary nucleation in the induction period was significantly reduced, resulting in the symmetric formation of chiral crystals due to high spontaneous primary nucleation. Agitation of the supersaturated solution promoted the induction of nucleation and a slow cooling rate induced nucleation at a high temperature (corresponding to a low supersaturation). Therefore, the initial chiral symmetry breaking was

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

enhanced when increasing the agitation speed and decreasing the cooling rate. Deracemization was predominantly determined by the initial chiral symmetry breaking. When it began with low initial CEE, the deracemization of chiral crystals was limited to a low final CEE. When increasing the initial CEE, the deracemization was also enhanced. Thus, 100 % deracemization was attained when the initial CEE was above 55 %. The initial chiral symmetry breaking was also examined by varying the final setting temperature and agitator. The final setting temperature variations confirmed that the induction supersaturation was a critical factor for the initial chiral symmetry breaking. When using three different agitators, the initial chiral symmetry breaking was markedly influenced by the turbulent flow. The correlation of the initial chiral symmetry breaking with the induction supersaturation suggested that secondary nucleation played a key role in chiral symmetry breaking.

Acknowledgement This work was supported by the Engineering Research Center of the Excellence Program of the Korean Ministry of Science, ICT & Future Planning (MSIP)/National Research Foundation of Korea (NRF) (Grant NRF-2014R1A5A1009799)

ACS Paragon Plus Environment

Page 14 of 32

Page 15 of 32 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

REFERENCES (1) Kipping, W. S.; Pope, W. J. LXIII-Enantiomorphism, J. Chem. SOC. Trans. 1898, 73 606-617. (2) Kondepudi, D. K.; Kaufman, R. J.; Singh N. Chiral symmetry breaking in sodium chlorate crystallization. Science. 1990, 250 (4983), 975-976. (3) Kondepudi, D. K.; Bullock, K. L.; Digits, J. A.; Hall, J. K.; Miller, J. M. Kinetics of chiral symmetry breaking in crystallization. J. Am. Chem. Soc. 1993, 115 (22), 10211-10216. (4) McBride, J. M.; Carter, R. L. Spontaneous Resolution by Stirred Crystallization. Angew. Chem., Int. Ed. Engl. 1991, 30 (3), 293-295. (5) Cartwright, J. H. E.; Garcia-Ruiz, J. M.; Piro, O.; Sainz-Dias, C. I.; Tuval, I. Chiral Symmetry Breaking during Crystallization: An Advection-Mediated Nonlinear Autocatalytic Process. Phys. Rev. Lett. 2004, 93 (3), 035502(4). (6) Metcalfe, G.; Ottino, J. M. Autocatalitic Processes in Mixing Flows, Phys. Rev. Lett. 1994, 72, 2875-2878. (7) Viedma, C. Experimental evidence of chiral symmetry breaking in crystallization from primary nucleation. J. Cryst. Growth. 2004, 261 (1), 118-121. (8) Szurgot, J.; Szurgot, M. Enantiomorphism in sodium chlorate crystals. Cryst. Res. Technol. 1995, 30 (1), 71-79. (9) Qian, R.-Y.; Botsaris, G. D. Nuclei breeding from a chiral crystal seed of NaClO3. Chem. Eng. Sci. 1998, 53 (9), 1745-1756. (10) Denk, E. G.; Botsaris, G. D. Fundamental Studies in Secondary Nucleation, J. Crystal Growth, 1972, 13/14, 493-499. (11) Mayer, C.; Lacmann, R. Secondary Nucleation of Sodium Chlorate Studied with the Aid of Asymmetric Crystallization. Chem. Eng. Technol. 1997, 20 (9), 633-640. (12) Callahan, C. J.; Ni, X.-W. An investigation into the effect of mixing on the secondary nucleation of sodium chlorate in a stirred tank and an oscillatory baffled crystallizer. CrystEngComm. 2014, 16, 690-697. (13) Niinomi, H.; Yamazaki, T.; Harada, S.; Ujihara, T.; Miura, H.; Kimura, Y.; Kuribayashi, T.; Uwaha, M.; Tsukamoto, K. Achiral Metastable Crystals of Sodium Chlorate Forming Prior to Chiral Crystals in Solution Growth. Cryst. Growth Des. 2013, 13 (12), 5188–5192. (14) Niinomi,H.; Miura, H.; Kimura, Y.; Uwaha, M.; Katsuno, H.; Harad, S.; Ujihara, T.; Tsukamoto,

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

K. Emergence and Amplification of Chirality via Achiral–Chiral Polymorphic Transformation in Sodium Chlorate Solution Growth. Cryst. Growth Des. 2014, 14 (7), 3596–3602. (15) Song, Y.; Chen, W.; Chen, X. Ultrasonic Field Induced Chiral Symmetry Breaking of NaClO3 Crystallization. Cryst. Growth Des. 2008, 8 (5), 1448-1450. (16) Zorzano, M.-P.; Susana, O.-E.; Marta, R.-B.; Cesar, M.-S.; Sabino, V.-V. Enantioselective Crystallization of Sodium Chlorate in the Presence of Racemic Hydrophobic Amino Acids and Static Magnetic Fields. Challenges 2014, 5 (1) , 175-192. (17) Alexander, A. J. Crystallization of Sodium Chlorate with d-Glucose Co-Solute Is Not Enantioselective. Cryst. Growth Des.2008, 8 (8), 2630–2632. (18) Viedma, C. Chiral symmetry breaking during crystallization: complete chiral purity induced by nonlinear autocatalysis and recycling. Phys. Rev. Lett. 2005, 94 (6), 065504(4). (19) Hein, J. E.; Cao, B. H.; Viedma, C.; Kellogg, R. M.; Blackmond, D. G. Pasteur’s Tweezers Revisited: On the Mechanism of Attrition-Enhanced Deracemization and Resolution of Chiral Conglomerate Solids. J. Am. Chem. Soc. 2012, 134 (30), 12629–12636. (20) Noorduin, W. L.; Enckevort, W. J.; Meekes, H.; Kaptein, B.; Kellogg, R. M.; Tully, J. C.; McBride, J. M.; Vlieg, E. The driving mechanism behind attrition-enhanced deracemization. Angew. Chem. Int. Ed. 2010, 49 (45), 8435-8438. (21) Viedma, C. Chiral symmetry breaking and complete chiral purity by thermodynamic-kinetic feedback near equilibrium: implications for the origin of biochirality. Astrobiology 2007, 7 (2), 312-319. (22) Cartwright, J. H. E.; Piro, O.; Tuval, I. Ostwald ripening, chiral crystallization, and the commonancestor effect. Phys. Rev. Lett. 2007, 98 (16), 165501(4). (23) Uwaha, M.; Katsuno, H. Mechanism of chirality conversion by grinding crystals:Ostwald ripening vs crystallization of chiral clusters. J. Phys. Soc. Jpn. 2009, 78 (2), 023601(4). (24) Suwannasang, K.; Flood, A. E.; Rougeot, C.; Coquerel, G. Using programmed heating–cooling cycles with racemization in solution for complete symmetry breaking of a conglomerate forming system. Cryst. Growth Des. 2013, 13 (8), 3498-3504. (25) Suwannasang, K.; Flood, A. E.; Rougeot, C.; Coquerel, G. Use of programmed damped temperature cycles for the deracemization of a racemic suspension of a conglomerate forming system. Org. Process Res. Dev. 2017, 21 (4), 623-630. (26) Rougeot, C.; Guillen, F.; Plaquevent, J.-C.; Coquerel, G. Ultrasound-enhanced deracemization: toward the existence of agonist effects in the interpretation of spontaneous symmetry breaking.

ACS Paragon Plus Environment

Page 16 of 32

Page 17 of 32 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

Cryst. Growth Des. 2015, 15 (5), 2151-2155. (27) Xiouras, C.; Aeken, J. V.; Panis, J.; Horst, J. H.;

Gerven, T. V.;

Stefanidis, G. D. Attrition-

enhanced deracemization of NaClO3: comparison between ultrasonic and abrasive grinding. Cryst. Growth Des. 2015, 15 (11), 5476-5484.

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

List of Figures Figure 1. Experiment apparatus for cooling crystallization of sodium chlorate in mixing tank (MT) crystallizer Figure 2. Deracemization of sodium chlorate when varying agitation speed in cooling crystallization. The cooling rate and final setting temperature were fixed at 0.197oC/min and 13oC, respectively. A stainless steel propeller was used for the agitation. Figure 3. (a) Chiral symmetry breaking in induction period and (b) supersaturation in induction period when varying agitation speed in cooling crystallization. The agitation speed and final setting temperature were fixed at 700rpm and 13oC, respectively. A stainless steel propeller was used for the agitation. Figure 4. Crystal-size profiles of two forms during deracemization with various agitation speeds. (a) 300rpm agitation speed, (b) 500rpm agitation speed, (c) 700rpm agitation speed, and (d) 1000rpm agitation speed. The cooling rate and final setting temperature were fixed at 0.197oC/min and 13oC, respectively. A stainless steel propeller was used for the agitation. Figure 5. Deracemization of sodium chlorate with sudden change of agitation during crystallization. (a) CEE profile with sudden change of agitation speed when CEE reached limit of deracemization, (b) crystal-size profiles of two forms with sudden change of agitation speed when CEE reached limit of deracemization, (c) CEE profile with sudden change of agitation speed after induction period, (b) crystal-size profiles of two forms with sudden change of agitation speed after induction period. The cooling rate and final setting temperature were fixed at 0.197oC/min and 13oC, respectively. A stainless steel propeller was used for the agitation Figure 6. Deracemization of sodium chlorate when varying cooling rate in cooling crystallization. The agitation speed and final setting temperature were fixed at 700rpm and 13oC, respectively. A stainless steel propeller was used for the agitation. Figure 7. (a) Supersaturation in induction period and (b) chiral symmetry breaking in induction period when varying cooling rate in cooling crystallization. The agitation speed and final setting temperature were fixed at 700rpm and 13oC, respectively. A stainless steel propeller was used for the agitation. Figure 8. Crystal-size profiles of two forms during deracemization with various cooling rates. (a) 0.074oC/min cooling rate, (b) 0.148oC/min cooling rate, (c) 0.197oC/min cooling rate, (d) 0.295oC/min cooling rate, and (e) 0.59oC/min cooling rate. The agitation speed and final setting temperature were fixed at 700rpm and 13oC, respectively. A stainless steel propeller

ACS Paragon Plus Environment

Page 18 of 32

Page 19 of 32 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

was used for the agitation. Figure 9. Deracemization of sodium chlorate with various final setting temperatures in cooling crystallization. (a) Temperature profiles and (b) CEE profiles. The agitation speed and cooling rate were fixed at 700rpm and 0.295oC/min, respectively. A stainless steel propeller was used for the agitation Figure 10. (a) Correlation of initial CEE to induction supersaturation, and (b) correlation of final CEE to initial CEE. Figure 11. Influence of agitator type and material on deracemization. (a) CEE profile with different agitator types and materials, and (b) initial CEE with different agitator types and materials. The agitation speed, cooling rate, and final setting temperature were fixed at 700rpm, 0.295oC/min and 13oC, respectively.

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

Figure 1. Experiment apparatus for cooling crystallization of sodium chlorate in mixing tank (MT) crystallizer

ACS Paragon Plus Environment

Page 20 of 32

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

Crystal Enantiomeric Excess [%]

Page 21 of 32

100

80

60

40

20 1000 RPM 700 RPM 500 RPM 300 RPM

0 0

20

40

60

80

100

Time [hr]

Figure 2. Deracemization of sodium chlorate when varying agitation speed in cooling crystallization. The cooling rate and final setting temperature were fixed at 0.197oC/min and 13oC, respectively. A stainless steel propeller was used for the agitation.

ACS Paragon Plus Environment

Crystal Growth & Design

CEE in Induction Period [%]

100

80

60

40

20

0 200

400

600

800

1000

Agitation Speed [RPM]

20

1.10

1.08 16 1.06

12

1.04

1.02 8 1.00 200

400

600

800

Supersaturation in Induction Period [-]

(a)

Temperature in Induction Period [oC]

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 32

1000

Agitation Speed [RPM]

(b)

Figure 3. (a) Chiral symmetry breaking in induction period and (b) supersaturation in induction period when varying agitation speed in cooling crystallization. The agitation speed and final setting temperature were fixed at 700rpm and 13oC, respectively. A stainless steel propeller was used for the agitation.

ACS Paragon Plus Environment

Page 23 of 32

100

Crystal Size [µm]

Crystal Size [µm]

100

80

60

Major form Minor form

40

0

20

40

60

80

80

60

Major form Minor form

40

0

100

20

40

60

Time [hr]

Time [hr]

(a)

(b)

80

100

100 100

Crystal Size [µm]

80

Crystal Size [µm]

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

80

60

Major form Minor form

40

0

20

40

60

80

100

60

40

20

Major form Minor form

0 0

4

8

12

16

20

Time [hr]

Time [hr]

(c)

(d)

Figure 4. Crystal-size profiles of two forms during deracemization with various agitation speeds. (a) 300rpm agitation speed, (b) 500rpm agitation speed, (c) 700rpm agitation speed, and (d) 1000rpm agitation speed. The cooling rate and final setting temperature were fixed at 0.197oC/min and 13oC, respectively. A stainless steel propeller was used for the agitation.

ACS Paragon Plus Environment

100

Page 24 of 32

1200 rpm

700 rpm

80

60

40

20

0 0

20

40

60

80

Time [hr]

(a)

100

Crystal Size [µm]

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 Enantiomeric Excess [%]

Crystal Growth & Design

700 rpm

1200 rpm

80

60

Major form Minor form

40

0

20

40

60

80

Time [hr]

(b)

Figure 5. Deracemization of sodium chlorate with sudden change of agitation during crystallization. (a) CEE profile with sudden change of agitation speed when CEE reached limit of deracemization, (b) crystal-size profiles of two forms with sudden change of agitation speed when CEE reached limit of deracemization, (c) CEE profile with sudden change of agitation speed after induction period, (b) crystal-size profiles of two forms with sudden change of agitation speed after induction period. The cooling rate and final setting temperature were fixed at 0.197oC/min and 13oC, respectively. A stainless steel propeller was used for the agitation.

ACS Paragon Plus Environment

Page 25 of 32

Crystal Enantiomeric Excess [%]

100

80

60

40

20

700 rpm

1200 rpm

0 0

2

4

6

8

10

12

Time [hr]

(c)

80

Major form Minor form 60

Crystal Size [µm]

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

40

20

700 rpm

1200 rpm

0

0

2

4

6

8

Time [hr]

(d)

Figure 5. Continued

ACS Paragon Plus Environment

10

12

Crystal Growth & Design

100

Crystal Enantiomeric Excess [%]

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 32

80

60

40

0.074 oC/min 0.148 oC/min 0.197 oC/min 0.295 oC/min 0.59 oC/min

20

0 0

20

40

60

80

100

Time [hr]

Figure 6. Deracemization of sodium chlorate when varying cooling rate in cooling crystallization. The agitation speed and final setting temperature were fixed at 700rpm and 13oC, respectively. A stainless steel propeller was used for the agitation.

ACS Paragon Plus Environment

1.10

1.08

16 1.06

12 1.04

8 1.02

4

1.00

0.0

0.1

0.2

0.3

0.4

0.5

Supersaturation in Induction Period [-]

20

0.6

Cooling Rate [oC/min]

(a)

100

CEE in Induction Period [%]

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

Temperature in Induction Period [oC]

Page 27 of 32

80

60

40

20

0 0.0

0.1

0.2

0.3

0.4

0.5

0.6

Cooling Rate [oC/min]

(b)

Figure 7. (a) Supersaturation in induction period and (b) chiral symmetry breaking in induction period when varying cooling rate in cooling crystallization. The agitation speed and final setting temperature were fixed at 700rpm and 13oC, respectively. A stainless steel propeller was used for the agitation.

ACS Paragon Plus Environment

Crystal Growth & Design

80

80

60

60

Crystal Size [µm]

100

Crystal Size [mm]

100

40 20

40

20

Major form Minor form

0 0

1

2

Major form Minor form

0

3

4

5

6

7

8

0

2

4

Time [hr]

6

8

10

12

Time [hr]

(a)

(b) 120

100

Crystal Size [µm]

100

Crystal Size [µm]

80

60

Major form Minor form

40

0

20

40

60

80

80

60

Major form Minor form

40

100

0

20

Time [hr]

40

60

80

100

Time [hr]

(c)

(d) 120

100

Crystal Size [µm]

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 28 of 32

80

60

Major form Minor form 40 0

20

40

60

80

100

Time [hr]

(e) Figure 8. Crystal-size profiles of two forms during deracemization with various cooling rates. (a) 0.074oC/min cooling rate, (b) 0.148oC/min cooling rate, (c) 0.197oC/min cooling rate, (d) 0.295oC/min cooling rate, and (e) 0.59oC/min cooling rate. The agitation speed and final setting temperature were fixed at 700rpm and 13oC, respectively. A stainless steel propeller was used for the agitation.

ACS Paragon Plus Environment

Page 29 of 32

28

1: 18.1 oC 2: 16.3 oC 3: 13 oC

Temperature [oC]

24

20

1 2 16

3 12 0

10

20

30

40

50

60

70

Time [min]

(a)

Crystal Enantiomeric Excess [%]

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

100

80

60

40

20

18.1 oC 16.3 oC 13 oC

0 0

10

20

30

40

50

60

70

Time [min]

(b)

Figure 9. Deracemization of sodium chlorate with various final setting temperatures in cooling crystallization. (a) Temperature profiles and (b) CEE profiles. The agitation speed and cooling rate were fixed at 700rpm and 0.295oC/min, respectively. A stainless steel propeller was used for the agitation.

ACS Paragon Plus Environment

Crystal Growth & Design

6 5

Ln(initial CEE) [%]

4 3 2 1

Final setting temp. Agitation speed Cooling rate

0 0

2000

4000

6000

8000

10000

12000

2

1/(lnS)

(a)

100 80

Final CEE [%]

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 30 of 32

60 40 20 Cooling rate Agitation speed Finial setting temp.

0 0

20

40

60

80

100

Initial CEE [%]

(b)

Figure 10. (a) Correlation of initial CEE to induction supersaturation, and (b) correlation of final CEE to initial CEE.

ACS Paragon Plus Environment

Page 31 of 32

Steel propeller Steel impeller Plastic propeller

Crystal Enantiomeric Excess [%]

100

80

60

40

20

0 0

20

40

60

80

100

120

Time [hr]

(a)

100

CEE in Induction Period [%]

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

80

60

40

20

0 Steel Propeller

Steel Impeller

Plastic Propeller

Agitator Type

(b)

Figure 11. Influence of agitator type and material on deracemization. (a) CEE profile with different agitator types and materials, and (b) initial CEE with different agitator types and materials. The agitation speed, cooling rate, and final setting temperature were fixed at 700rpm, 0.295oC/min and 13oC, respectively.

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 32 of 32

Graphical Abstract

For Table of Contents Use Only

Title: Chiral Symmetry Breaking and Deracemization of Sodium Chlorate in Turbulent Flow Authors: Jaekyu Ahn, Do Hyun Kim, Gerard Coquerel and Woo-Sik Kim* Primary nucleation

Secondary nucleation

Chiral non-symmetry

Low S

Primary nucleation

Chiral symmetry

High S

§ The chiral symmetry breaking in the induction period was caused by the secondary nucleation in the crystallization. § The secondary nucleation in the induction period was predominantly controlled by the turbulent flow supersaturation. § The supersaturation in the induction was markedly dependent on the agitation speed and cooling rate in the crystallization.

ACS Paragon Plus Environment