Pia
Mueller
*a,
Aikaterini
Vriza
a,
Adam D.
Clayton
a,
Oliver S.
May
bc,
Norman
Govan
b,
Stuart
Notman
b,
Steven V.
Ley
c,
Thomas W.
Chamberlain
a and
Richard A.
Bourne
*a
aInstitute of Process Research and Development, School of Chemistry & School of Chemical and Process Engineering, University of Leeds, LS2 9JT, UK. E-mail: p.mueller@leeds.ac.uk; r.a.bourne@leeds.ac.uk
bDefence Science and Technology Lab, Porton Down, Salisbury, SP4 0JQ, UK
cYusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK
First published on 10th January 2023
Automated platforms allow for fast and efficient optimisation of single and multi-objective chemical systems. Herein, we report the application of automated optimisation platforms for the chemical screening of sulfide oxidations in flow. The identification of different optima for each substrate highlights the requirement of rapid, individual experimental optimisation for successful compound screening.
Only a few instances are published that combine chemical or solvent screening with algorithm supported optimisation.12,13 More recently, Kershaw et al. extended multi-objective optimisations to mixed variables screening, including ligands and solvents, allowing continuous and discrete variables to be optimised simultaneously, for one substrate.14 Simon et al. screened several nitro aromatic substrates in heterogeneous hydrogenation to evaluate the tolerance and stability of the catalysts by varying conditions in set steps, finding enough data to give an informed benchmarking. While these studies exceed a simple, one-point screening, they are not focussed on finding individual optima for different substrates.12
For this study the oxidation of organic sulfides to sulfoxides is chosen due to its appeal in pharmaceutical chemistry (Fig. 1).15 Furthermore, oxidation of organic sulfides containing β-chloro or β-hydroxyl substituents is widely used in sulfur-mustard decontamination studies.16
Whilst standard experimental studies have previously been able to find optima for one objective by changing one variable at a time,17 self-optimized studies are able to correlate variables in a coherent functionality to multiple objectives.18,19 A wide range of heterogeneous and homogenous catalytic studies20,21 and electrochemical reaction routes22,23 have been performed on the selective oxidation of organic sulfides achieving fast, selective oxidation in minutes or even seconds at close to room temperatures.
In this study we selected a homogenous phosphotungstic acid catalyst, which has been shown to exhibit catalytic activity for selective oxidations, owing to the presence of Brønsted acid sites and proton channelling.24 Whilst not being the catalyst with highest reactivity, this comparatively low-cost catalyst demonstrates good stability, enables a green solvent system to be used (water with butanol as organic sulfide carrier),17,25 and can potentially be recycled via phase separation.26 In this study a more complex experimental set-up is chosen to highlight the versatility of the optimisation of automated continuous flow platforms. Usually, these optimisations are performed in tubular reactors, assuming plug flow conditions, using one liquid phase or 2 liquid phases as slugs.27,28 Herein, the biphasic reaction with varying phase ratios requires active mixing to increase the extent of contact between the phases. For this purpose, mini-continuously stirred tank reactor (CSTR) cells were found effective.29 Using these cells for automated optimisations in this way constitutes a new field which has the potential to expand the breadth of flow chemistry conditions currently accessible.30,31
The experimental study was performed on a previously described flow platform using an on-line GCMS fitted with an automated liquid sampling system for determining conversion and selectivity (Fig. 1b). Flow rates, reaction temperature and sampling were monitored by the MATLAB interface, which used the calibrated mass spectrum peak areas for analysis.19 (more information can be found in the ESI,† S1.2–S1.4).
Data collection was automatically performed during a time span of ∼24 hours. The stability of the system was tested using diphenylsulfide (DPS). Following 10 replications a standard deviation of 1.7% of the yield was observed (ESI,† S2.3). Each compound optimisation was initialised using a Latin hypercube (LHC, size = 9, for 2n + 1, n = 4 parameter) for maximal experimental distribution.32 Subsequent experiments were then designed using a Bayesian Optimisation with an Adaptive Expected Improvement algorithm (BOAEI).
The oxidation of methyl phenyl sulfide (TA) is a common model reaction in literature due to its commercial availability and low toxicity. As the first model substrate, the catalytic system was able to selectively oxidise TA to the corresponding sulfoxide (2a).
The optimisation converged on long residence times (9–10 min) and a temperature above 75 °C (Fig. 2). An excess of hydrogen peroxide of 1.5 to 2.5 equivalents to TA was seen as advantageous. This is on average slightly below what can be found in literature (typically used equiv. of 2.5 or an excess), indicating that the presence of a catalyst and water as a solvent increases the reactivity.33 Compared to the reaction rates the decomposition of hydrogen peroxide was negligible in our system. Similar 3D scatter plots as in Fig. 2 for the other 3 substrates can be found in the ESI,† S2.1. The optimisations for all 4 compounds were successful in finding a global optimum defined by a maximal yield which could not be improved by 5 subsequent runs (see ESI,† S2.1). However, in the given boundaries the system was unable to achieve full conversion for DPS (2b) and CEPS (2d). For DPS, we propose this is due to the 2nd phenyl ring sterically hindering the reaction and decreasing water solubility (see ESI,† S3), leading to lower reactivity in this bi-phasic system. For CEPS the low yield was in contrast caused by the formation of side products. The mass spectrometry identified the elimination product, allyl phenylsulfide, and its oxidised sulfoxides and sulfone. This suggests a parallel oxidation of the eliminated sulfide. Additional signals were seen for some experimental runs but due to the low analytical response could not be quantified or qualified (ESI,† S1.4).
With the more complex formation of side products, the BOAEI algorithm correlated with the need to perform an increased number of experimental iterations to find global optima than in the previous optimisations of HEPS and TA.
The accuracy of the trained surrogate models by the BOAEI algorithm for each sulfide was tested by predicting the response (yield) for each experiment and comparing the result to the experiments in a Pareto plot (S2.4†). A close agreement between BAOEI models and experiments were found. Due to this synergy, a simulation of the bounds was performed to show the surface response of each parameter on the yield. Fig. 3 shows the response for the optimisation parameters paired in comparison for each sulfide. The 4-dimensionality was broken down by only showing the response at the fixed optima of the other 2 parameters (for further explanations see ESI,† S2.5).
The surface response for HEPS and at lower yields CEPS, show a localised optimum, while for TA the optimum is clearly at the upper bounds for the temperature and residence time. Whereas with DPS there is a band in the optimisation which balances between the desired oxidation and over-oxidation. In accordance with rate laws, temperature and residence time impact the yield significantly (up to 60% deviation in yield in the given bounds). In comparison, the influence of equivalents of hydrogen peroxide and catalyst on the sulfide is lower (maximal 30%). Interestingly, the zone in which the hydrogen peroxide equivalent leads to optimal yields is quite narrow but moves from 2 for HEPS and TA, over 2.6 for DPS and to above 2.9 for CEPS.
This is likely correlated with the lower selectivity of the reactions. The catalyst ratio shows a slight increase in optimal yields between 1 and 1.5%. The ratio of the catalyst stream to the sulfide did not show a clear trend as catalytic activity was counteracted by the dilution of the system through the 3rd stream.
Focusing on the optima the 4 varied process conditions lead to a multi-dimensional problem often visualised using a parallel coordinate plot. Fig. 4 shows such a graph for the four screened sulfides. The conditions necessary for each optimum (dark lines) and their standard deviation (light lines) are shown. Small deviations relate to high impact of variation as seen for residence time and temperature for TA, HEPS and CEPS and are in alignment with the algorithm simulated surface responses (Fig. 4). The difference in reactivity of DPS and CEPS can be easily seen by comparing the associated total conversion with the yield of the formed sulfoxide. While DPS has a low conversion, CEPS is able to convert in the given bounds to up to 85%, but due to the less selective nature leading only to low optimal yields as well. The phosphotungstic acid catalyst is known for selective oxidation by providing additional acid sites. For the non-chlorinated sulfides, selectivities of 94% were seen at the optimum, for CEPS a maximum selectivity of 80% was found at the highest yields. A likely explanation is that the formation of hydrochloric acid leads to a further decrease in the pH making the hydrogen peroxide more likely to decompose or over oxidise to the corresponding sulfone. An indication of this is the significant formation of ethyl-phenyl-sulfone seen in the GCMS trace. Comparing our novel full chemical space screening method towards common best point screenings, the increase in chemical reaction understanding is clear. When only looking at the optima, the traditional approach of finding and optimising for a model substrate, in our case TA, and testing the best conditions for DPS and HEPS resulted in final yields of 32% and 70% respectively. Being 5–25% below the optimised screening defined global optima, making it a compelling case as a viable alternative to traditional screenings.
The real advantage of the proposed system is the use of a fully automated system equipped with a smart algorithm finding optima in less than 20 experiments for each substrate, which were continuously performed without human interference. This autonomous approach reduced experimental time and allowed for data collection overnight.
Advances in automated optimisation allow routine tasks to be handled by the equipment and gain a deeper process understanding without prior knowledge. This approach was utilised in the work to successfully screen the catalysed oxidation of four sulfides ranging from relatively simple model compound thioanisole to mustard simulant chloroethylphenylsulfide. The experimental focus lays on the complexity of the process, using organic and aqueous streams mixed in mini CSTR cascades varying in total 4 parameters in the optimisation. The study found different optima for the oxidation of each substrate in less than 25 experiments and a total of less than 24 hours of experimental time for each of them, showcasing the power of modern Bayesian algorithm in finding global optima. While temperature and residence time proved to have a clear and significant influence on all compounds, optimal equivalents of hydrogen peroxide were found for each species and did increase in correlation with the amount of formed by-products per sulfide. In the absence of β-chloro substituents, the catalyst achieved selective oxidisation to sulfoxide. However, no clear trend emerged since varying the catalyst concentration resulted in a dilution effect which impacted the total reaction yield. Overall, the automated screening delivered reliable global optima with up to 25% more yield than using the one-point screening based on the model compound optimisation. An increased process understanding could be gained by simulating the surface response and parameter sensitivity of optimal points. Herein, an approach is suggested that is superior to the traditional screening approach, providing a faster and more reliable way to detect the optimal reaction conditions for four structurally diverse sulfides.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2re00552b |
This journal is © The Royal Society of Chemistry 2023 |