Exact and Ubiquitous Condition for Solid-State Deracemization in Vitro and in Nature.

Solid-state deracemization is the amplification of an enantiomeric excess in suspensions of conglomerate-forming chiral compounds. Although numerous chemical and biochemical compounds deracemize, its governing mechanism has remained elusive. We introduce a novel formulation of the classical population-based model of deracemization through temperature cycles to prove that suspensions deracemize whenever a simple and ubiquitous condition is met: crystal dissolution must be faster than crystal growth. Such asymmetry is a known principle of crystallization, hence explaining the generality of deracemization. Through both experiments and a theoretical analysis, we demonstrate that this condition applies even for very small temperature cycles and for random temperature fluctuations. These findings establish solid-state deracemization as an attractive route to the manufacture of enantiopure products and as a plausible pathway toward the emergence of homochirality in nature.

Monitoring Aqueous Sucrose Solutions Using Droplet Microfluidics: Ice Nucleation, Growth, Glass Transition, and Melting.

Freezing and freeze-drying processes are commonly used to extend the shelf life of drug products and to ensure their safety and efficacy upon use. When designing a freezing process, it is beneficial to characterize multiple physicochemical properties of the formulation, such as nucleation rate, crystal growth rate, temperature and concentration of the maximally freeze-concentrated solution, and melting point. Differential scanning calorimetry has predominantly been used in this context but does have practical limitations and is unable to quantify the kinetics of crystal growth and nucleation. In this work, we introduce a microfluidic technique capable of quantifying the properties of interest and use it to investigate aqueous sucrose solutions of varying concentration. Three freeze-thaw cycles were performed on droplets with 75-µm diameters at cooling and warming rates of 1 °C/min. During each cycle, the visual appearance of the droplets was optically monitored as they experienced nucleation, crystal growth, formation of the maximally freeze-concentrated solution, and melting. Nucleation and crystal growth manifested as increases in droplet brightness during the cooling phase. Heating was associated with a further increase as the temperature associated with the maximally freeze-concentrated solution was approached. Heating beyond the melting point corresponded to a decrease in brightness. Comparison with the literature confirmed the accuracy of the new technique while offering new visual data on the maximally freeze-concentrated solution. Thus, the microfluidic technique presented here may serve as a complement to differential scanning calorimetry in the context of freezing and freeze-drying. In the future, it could be applied to a plethora of mixtures that undergo such processing, whether in pharmaceutics, food production, or beyond.

Assessment of Potential and Techno-Economic Performance of Solid Sorbent Direct Air Capture with CO2 Storage in Europe.

Direct air capture with CO2 storage (DACCS) is among the carbon dioxide removal (CDR) options, with the largest gap between current deployment and needed upscaling. Here, we present a geospatial analysis of the techno-economic performance of large-scale DACCS deployment in Europe using two performance indicators: CDR costs and potential. Different low-temperature heat DACCS configurations are considered, i.e., coupled to the national power grid, using waste heat and powered by curtailed electricity. Our findings reveal that the CDR potential and costs of DACCS systems are mainly driven by (i) the availability of energy sources, (ii) the location-specific climate conditions, (iii) the price and GHG intensity of electricity, and (iv) the CO2 transport distance to the nearest CO2 storage location. The results further highlight the following key findings: (i) the limited availability of waste heat, with only Sweden potentially compensating nearly 10% of national emissions through CDR, and (ii) the need for considering transport and storage of CO2 in a comprehensive techno-economic assessment of DACCS. Finally, our geospatial analysis reveals substantial differences between regions due to location-specific conditions, i.e., useful information elements and consistent insights that will contribute to assessment and feasibility studies toward effective DACCS implementation.

Life Cycle Assessment of Direct Air Carbon Capture and Storage with Low-Carbon Energy Sources.

Direct air carbon capture and storage (DACCS) is an emerging carbon dioxide removal technology, which has the potential to remove large amounts of CO2 from the atmosphere. We present a comprehensive life cycle assessment of different DACCS systems with low-carbon electricity and heat sources required for the CO2 capture process, both stand-alone and grid-connected system configurations. The results demonstrate negative greenhouse gas (GHG) emissions for all eight selected locations and five system layouts, with the highest GHG removal potential in countries with low-carbon electricity supply and waste heat usage (up to 97%). Autonomous system layouts prove to be a promising alternative, with a GHG removal efficiency of 79-91%, at locations with high solar irradiation to avoid the consumption of fossil fuel-based grid electricity and heat. The analysis of environmental burdens other than GHG emissions shows some trade-offs associated with CO2 removal, especially land transformation for system layouts with photovoltaics (PV) electricity supply. The sensitivity analysis reveals the importance of selecting appropriate locations for grid-coupled system layouts since the deployment of DACCS at geographic locations with CO2-intensive grid electricity mixes leads to net GHG emissions instead of GHG removal today.

Description of Adsorption in Liquid Chromatography under Nonideal Conditions.

A thermodynamically consistent description of binary adsorption in reversed-phase chromatography is presented, accounting for thermodynamic nonidealities in the liquid and adsorbed phases. The investigated system involves the adsorbent Zorbax 300SB-C18, as well as phenetole and 4- tert-butylphenol as solutes and methanol and water as inert components forming the eluent. The description is based on adsorption isotherms, which are a function of the liquid-phase activities, to account for nonidealities in the liquid phase. Liquid-phase activities are calculated with a UNIQUAC model established in this work, based on experimental phase equilibrium data. The binary interaction in the adsorbed phase is described by the adsorbed solution theory, assuming an ideal (ideal adsorbed solution theory) or real (real adsorbed solution theory) adsorbed phase. Implementation of the established adsorption model in a chromatographic code achieves a quantitative description of experimental elution profiles, with feed compositions exploiting the entire miscible region, and involving a broad range of different eluent compositions (methanol/water). The quantitative agreement of the model and experimental data serves as a confirmation of the underlying physical (thermodynamic) concepts and of their applicability to a broad range of operating conditions.

Molecular-dynamics simulations of urea nucleation from aqueous solution.

Despite its ubiquitous character and relevance in many branches of science and engineering, nucleation from solution remains elusive. In this framework, molecular simulations represent a powerful tool to provide insight into nucleation at the molecular scale. In this work, we combine theory and molecular simulations to describe urea nucleation from aqueous solution. Taking advantage of well-tempered metadynamics, we compute the free-energy change associated to the phase transition. We find that such a free-energy profile is characterized by significant finite-size effects that can, however, be accounted for. The description of the nucleation process emerging from our analysis differs from classical nucleation theory. Nucleation of crystal-like clusters is in fact preceded by large concentration fluctuations, indicating a predominant two-step process, whereby embryonic crystal nuclei emerge from dense, disordered urea clusters. Furthermore, in the early stages of nucleation, two different polymorphs are seen to compete.

Ripening of Semiconductor Nanoplatelets.

Ostwald ripening describes how the size distribution of colloidal particles evolves with time due to thermodynamic driving forces. Typically, small particles shrink and provide material to larger particles, which leads to size defocusing. Semiconductor nanoplatelets, thin quasi-two-dimensional (2D) particles with thicknesses of only a few atomic layers but larger lateral dimensions, offer a unique system to investigate this phenomenon. Experiments show that the distribution of nanoplatelet thicknesses does not defocus during ripening, but instead jumps sequentially from m to (m + 1) monolayers, allowing precise thickness control. We investigate how this counterintuitive process occurs in CdSe nanoplatelets. We develop a microscopic model that treats the kinetics and thermodynamics of attachment and detachment of monomers as a function of their concentration. We then simulate the growth process from nucleation through ripening. For a given thickness, we observe Ostwald ripening in the lateral direction, but none perpendicular. Thicker populations arise instead from nuclei that capture material from thinner nanoplatelets as they dissolve laterally. Optical experiments that attempt to track the thickness and lateral extent of nanoplatelets during ripening appear consistent with these conclusions. Understanding such effects can lead to better synthetic control, enabling further exploration of quasi-2D nanomaterials.

A low-energy chilled ammonia process exploiting controlled solid formation for post-combustion CO2 capture.

A new ammonia-based process for CO2 capture from flue gas has been developed, which utilizes the formation of solid ammonium bicarbonate to increase the CO2 concentration in the regeneration section of the process. Precipitation, separation, and dissolution of the solid phase are realized in a dedicated process section, while the packed absorption and desorption columns remain free of solids. Additionally, the CO2 wash section applies solid formation to enable a reduction of the wash water consumption. A rigorous performance assessment employing the SPECCA index (Specific Primary Energy Consumption for CO2 Avoided) has been implemented to allow for a comparison of the overall energy penalty between the new process and a standard ammonia-based capture process without solid formation. A thorough understanding of the relevant solid-solid-liquid-vapor phase equilibria and an accurate modeling of them have enabled the synthesis of the process, and have inspired the development of the optimization algorithm used to screen a wide range of operating conditions in equilibrium-based process simulations. Under the assumptions on which the analysis is based, the new process with controlled solid formation achieved a SPECCA of 2.43 MJ kgCO2-1, corresponding to a reduction of 17% compared to the process without solid formation (with a SPECCA of 2.93 MJ kgCO2-1). Ways forward to confirm this significant improvement, and to increase the accuracy of the optimization are also discussed.

On the potential of phase-change adsorbents for CO2 capture by temperature swing adsorption.

We investigate the potential of a class of recently discovered metal-organic-framework materials for their use in temperature swing adsorption (TSA) processes for CO2 capture; the particularity of the considered materials is their reversible and temperature dependent step-shaped CO2 adsorption isotherm. Specifically, we present a comprehensive modeling study, where the performance of five different materials with step-shaped isotherms [McDonald et al., Nature, 2015, 519, 303] in a four step TSA cycle is assessed. The specific energy requirement of the TSA process operated with these materials is lower than for a commercial 13X zeolite, and a smaller temperature swing is required to reach similar levels of CO2 purity and recovery. The effect of a step in the adsorption isotherm is illustrated and discussed, and design criteria that lead to an optimal and robust operation of the considered TSA cycle are identified. The presented criteria could guide material scientists in designing novel materials whose step position is tailored to specific CO2 separation tasks.

Overcoming time scale and finite size limitations to compute nucleation rates from small scale well tempered metadynamics simulations.

Condensation of a liquid droplet from a supersaturated vapour phase is initiated by a prototypical nucleation event. As such it is challenging to compute its rate from atomistic molecular dynamics simulations. In fact at realistic supersaturation conditions condensation occurs on time scales that far exceed what can be reached with conventional molecular dynamics methods. Another known problem in this context is the distortion of the free energy profile associated to nucleation due to the small, finite size of typical simulation boxes. In this work the problem of time scale is addressed with a recently developed enhanced sampling method while contextually correcting for finite size effects. We demonstrate our approach by studying the condensation of argon, and showing that characteristic nucleation times of the order of magnitude of hours can be reliably calculated. Nucleation rates spanning a range of 10 orders of magnitude are computed at moderate supersaturation levels, thus bridging the gap between what standard molecular dynamics simulations can do and real physical systems.

Urea homogeneous nucleation mechanism is solvent dependent.

The composition of the mother phase plays a primary role in crystallization processes, affecting both crystal nucleation and growth. In this work, the influence of solvents on urea nucleation has been investigated by means of enhanced sampling molecular dynamics simulations. We find that, depending on the solvent, the nucleation process can either follow a single-step or a two-step mechanism. While in methanol and ethanol a single-step nucleation process is favored, in acetonitrile a two-step process emerges as the most likely nucleation pathway. We also find that solvents have a minor impact on polymorphic transitions in the early stages of urea nucleation. The impact of finite size effects on the free energy surfaces is systematically considered and discussed in relation to the simulation setup.

Modelling the stochastic behaviour of primary nucleation.

We study the stochastic nature of primary nucleation and how it manifests itself in a crystallisation process at different scales and under different operating conditions. Such characteristics of nucleation are evident in many experiments where detection times of crystals are not identical, despite identical experimental conditions, but instead are distributed around an average value. While abundant experimental evidence has been reported in the literature, a clear theoretical understanding and an appropriate modelling of this feature is still missing. In this contribution, we present two models describing a batch cooling crystallisation, where the interplay between stochastic nucleation and deterministic crystal growth is described differently in each. The nucleation and growth rates of the two models are estimated by a comprehensive set of measurements of paracetamol crystallisation from aqueous solution in a 1 mL vessel [Kadam et al., Chemical Engineering Science, 2012, 72, 10-19]. Both models are applied to the cooling crystallisation process above under different operating conditions, i.e. different volumes, initial concentrations, cooling rates. The advantages and disadvantages of the two approaches are illustrated and discussed, with particular reference to their use across scales of nucleation rate measured in very small crystallisers.

Flue gas CO2 mineralization using thermally activated serpentine: from single- to double-step carbonation.

Carbon dioxide capture and utilization by mineralization seeks to combine greenhouse gas emission control with the production of value-added materials in the form of solid carbonates. This experimental work demonstrates that the world's most abundant mineralization precursor, the magnesium (Mg) silicate serpentine, in its thermally activated, partially dehydroxylated form can be carbonated without the use of chemical additives at process temperatures (T) below 90 °C and CO2 partial pressures (pCO2) below 1 bar. A first series of single-step batch experiments was performed varying the temperature and slurry density to systematically assess the precipitation regime of the relevant Mg-carbonates and the fate of silicon (Si) species in solution. The results suggested that the reaction progress was hindered by a passivating layer of re-precipitated silica or quartz, as well as by equilibrium limitations. Concurrent grinding proved effective in tackling the former problem. A double-step strategy proved successful in addressing the latter problem by controlling the pH of the solution. This is achieved by continuously removing the Mg from the dissolution reactor and letting it precipitate at a higher T and a lower pCO2 in a separate reactor, thus yielding a combined T-pCO2-swing-the working principle of a new flue gas mineralization route is presented herein. Simulations and experiments of the different individual steps of the process are reported, in order to make an assessment of its feasibility.

Kinetics and Thermodynamics of Lactose Mutarotation through Chromatography.

The mutarotation kinetics and thermodynamics of the reaction α-lactose ⇌ ß-lactose have been measured in dilute solutions using liquid chromatography without any derivatization step, using a C18 column and pure water as the mobile phase. The effect of temperature (0.5-45 °C) of the starting powder composition (α-lactose-rich or ß-lactose-rich powder) and of the solvent composition (water with up to 35% weight fraction of seven organic solvents) has been experimentally investigated. Increasing the temperature leads to faster kinetics, following an Arrhenius model, and to slightly decreasing concentration-based equilibrium ratio. Conversely, increasing the weight fraction of organic solvent at 25 °C resulted in slower kinetics and smaller concentration-based equilibrium ratio. The starting powder composition is shown not to influence the kinetics or thermodynamics of the process. The corresponding parameter estimation problem is thoroughly discussed, taking into account the small difference in response factors of the lactose diastereomers.

HPLC method development for the quantification of a mixture of reacting species: The case of lactose.

Preparative and analytical chromatography are impaired by analytes that undergo a chemical reaction during the chromatographic separation, leading to peak distortion and systematic errors during the subsequent quantification phase. The pitfalls are highlighted through a combination of analytical results and numerical simulations. Two different quantification strategies for partially overlapping and reacting peaks are compared. A novel method development strategy based on the valley-to-peak ratio instead of the more common resolution is proposed. The method has been used to experimentally investigate the chromatographic behavior of a mutarotating sugar, lactose. The separation of the unprotected lactose isomers, α and ß, has been optimized using a C18 column and pure water as the mobile phase. Phase dewetting phenomena during method development have also been studied and discussed.

On Comparing Packed Beds and Monoliths for CO2 Capture from Air Through Experiments, Theory, and Modeling.

This study compares the performance of amine-functionalized γ-alumina sorbents in the form of 3 mm γ-alumina pellets and of a γ-alumina wash-coated monolith for CO2 capture for direct air capture (DAC). Breakthrough experiments were conducted on the two contactors to analyze the adsorption kinetics and performance for different gas feeds. A constant pattern analysis revealed dominant mass transfer resistances in the gas film and in the pores, with axial dispersion also observed, particularly at higher concentrations. A 1D, physical model was used to fit the experiments and thus to estimate mass transfer and axial dispersion coefficients, which appear to be consistent with the hypotheses derived from constant pattern analysis. A dual kinetic model to describe mass transfer was found to better describe the tail behavior in the monolith, whereas a pseudo-first-order model was sufficient to describe breakthroughs on packed beds. A substantial two-order magnitude decrease in mass transfer coefficients was noted when reducing the feed concentration from 5.6% to 400 ppm CO2, thus underscoring the significant mass transfer limitations observed in DAC. Comparison between the contactors revealed notably higher mass transfer coefficients in the monolith compared to the packed beds, which are attributed to shorter diffusion lengths and lower equilibrium capacity. While the faster mass transfer coefficients observed in the monolith experiments led to reduced specific energy consumption and increased adsorption productivity compared to the packed bed at 400 ppm, no significant improvement was observed for the same process at the higher concentration of 5.6% CO2 in the feed. This finding highlights the need to tailor the contactor design to the specific gas separation requirements. This research contributes to the understanding and quantification of mass transfer kinetics at DAC concentrations in both packed bed and monolith contactors. It demonstrates the crucial role of the contactor in DAC systems and the importance of optimizing the adsorption step to enhance productivity and DAC performance.

On Solute Recovery and Productivity in Chiral Resolution through Solid-State Deracemization by Temperature Cycling.

Temperature cycling represents an effective means for the deracemization of chiral compounds that crystallize as conglomerates and racemize in solution. In such a process, a suspension enriched in the desired enantiomer is converted into an enantiopure one through periodic cycles of crystal dissolution and crystal growth. We show that performing temperature cycling at higher temperatures leads to faster deracemization and, consequently, higher productivity. However, this comes at the cost of lower recovery, as the solution contains potentially relevant amounts of solute due to the higher solubility at an elevated temperature. In this work, we introduce and compare two process variants that mitigate this issue. The first involves temperature cycling, followed by linear cooling, whereas the second is based on merging the temperature cycles and cooling crystallization. Experiments carried out with the chiral compound N-(2-methylbenzylidene)-phenylglycine amide show that the former variant is faster than the latter, and it is easier to design and implement. In this process, the choice of an appropriate cooling rate is essential to avoid nucleation of the undesired enantiomer.

Conceptual Validation of Stochastic and Deterministic Methods To Estimate Crystal Nucleation Rates.

This work presents a generalized framework to assess the accuracy of methods to estimate primary and secondary nucleation rates from experimental data. The crystallization process of a well-studied model compound was simulated by means of a novel stochastic modeling methodology. Nucleation rates were estimated from the simulated data through multiple methods and were compared with the true values. For primary nucleation, no method considered in this work was able to estimate the rates accurately under general conditions. Two deterministic methods that are widely used in the literature were shown to overpredict rates in the presence of secondary nucleation. This behavior is shared by all methods that extract rates from deterministic process attributes, as they are insensitive to primary nucleation if secondary nucleation is sufficiently fast. Two stochastic methods were found to be accurate independent of whether secondary nucleation is present, but they underestimated rates in the case where a large number of primary nuclei are formed. We hence proposed a criterion to probe the accuracy of stochastic methods for arbitrary data sets, thus providing the theoretical foundations required for their rational use. Finally, we showed how both primary and secondary nucleation rates can be inferred from the same set of detection time data by combining deterministic and stochastic considerations.

Thermodynamics Explains How Solution Composition Affects the Kinetics of Stochastic Ice Nucleation.

The freezing of aqueous solutions is of great relevance to multiple fields, yet the kinetics of ice nucleation, its first step, remains poorly understood. The literature focuses on the freezing of microdroplets, and it is unclear if those findings can be generalized and extended to larger volumes such as those used in the freezing of biopharmaceuticals. To this end, we study ice nucleation from aqueous solutions of ten different compositions in vials at the milliliter scale. The statistical analysis of the approximately 6,000 measured nucleation events reveals that the stochastic ice nucleation kinetics is independent of the nature and concentration of the solute. We demonstrate this by estimating the values of the kinetic parameters in the nucleation rate expression for the selected solution compositions, and we find that a single set of parameters can describe quantitatively the nucleation behavior in all solutions. This holds regardless of whether the nucleation rate is expressed as a function of the chemical potential difference, of the water activity difference, or of the supercooling. While the chemical potential difference is the thermodynamically correct driving force for nucleation and hence is more accurate from a theoretical point of view, the other two expressions allow for an easier implementation in mechanistic freezing models in pharmaceutical manufacturing.

Process Performance and Operational Challenges in Continuous Crystallization: A Study of the Polymorphs of L-Glutamic Acid.

The crystallization of the two polymorphs of l-glutamic acid (LGA) is carried out in a continuous crystallization process, and its performance according to different criteria is evaluated. The study aims at identifying suitable operating conditions for producing either αLGA or ßLGA with a high polymorphic purity. To this end, we investigate the process both from a theoretical perspective and through experiments using either a single stirred-tank crystallizer or a cascade of two stirred-tank crystallizers in series. In terms of theory, we extend the MSMPR-based steady-state stability analysis of Farmer et al. (Farmer, T. C. et al. AIChE J.2016, 62, 3505-3514) by accounting for the possibility of a nonrepresentative withdrawal of the solid phase from the crystallizer. Additionally, the process is simulated using population balance equations, thereby investigating the effect of operating conditions on polymorphic purity, yield, and productivity. Guided by the model-based conclusions, we identified suitable operating conditions and experimentally tested them. The experimental campaign has demonstrated that ßLGA could be successfully and continuously produced in both process configurations according to the theory with performance as expected, whereas that was not possible for αLGA. The difference between the two stems from different operational challenges, whose consequence is that steady-state operation is attained in the case of ßLGA but not in that of αLGA. In the former case, the needle-like ßLGA crystals, which exhibit no agglomeration, tend to be only slightly oversampled; in the latter case, the prismatic αLGA crystals undergo major agglomeration and hence are very difficult to suspend and effectively withdraw from the crystallizer.