Modeling Solvation in Supercritical CO2

A Molecular Dynamics Study of the Solvation Properties of Sugars in Supercritical Carbon Dioxide

Among the various strategies used to enhance the solvation power of supercritical carbon dioxide (scCO2), the use of CO2-philic compounds has been extensively studied over the recent two decades. Given the biocompatibility of this medium, extraction technologies based on scCO2 are particularly attractive, and a molecular-level understanding of intermolecular interactions is crucial for optimizing processing conditions. Functionalized sugars and cyclic oligosaccharides, such as cyclodextrins, can be rendered soluble in scCO2, opening new avenues for vectorization strategies and supramolecular chemistry in this medium. To support the exploration of CO2-philic compounds relevant to these research goals, we conducted a molecular dynamics investigation into the solvation properties of cyclodextrins functionalized with CO2-philic groups. We thoroughly analyzed the key solute-solvent interactions and their influence on the cavity shape. Additionally, we provided insights into the solvation behavior of peracetylated α and β-glucose across different regions of the carbon dioxide phase diagram. We were able to confirm the importance of the well-known (acetyl)C-O⋯C(CO2) interaction, as the most important signature of CO2-philicity of carbonyl compounds. Depending on the substituent, this interaction can be assisted by a cooperative (methyl)2HCH⋯O(CO2) intermolecular bond. In cyclodextrins, conformational flexibility, with a possible change in the conformation of some pyranose units, was observed in the macromolecular structure. On the other hand, these structural modifications were not present for α- and β-glucose.

Deciphering the Role of Vinylene Carbonate in Shaping CO2 Cluster Growth and Stability via Wavelet-Enhanced Microwave Spectroscopy

This study presents the first integration of wavelet denoising with resonator-enhanced microwave spectroscopy to explore the formation and characteristics of gas-phase heterogeneous vinylene carbonate-(CO2)1-5 clusters. Through this innovative approach, faint spectral lines were detected with efficiency, enhancing our understanding of complex interaction patterns that govern the growth and stability of these molecular clusters. It was found that CO2 aggregates preferentially around the vinylene carbonate molecule through C···O tetrel bonds at the carbonyl and ether groups, altering the typical self-aggregation topology of free CO2 clusters. This leads to more symmetrical growth around the vinylene carbonate monomer. Our findings provide insights into the solvation mechanisms of compounds in supercritical CO2 thereby advancing the understanding of subnanoscale CO2 aggregation patterns.

Singlet Oxygen Photophysics: From Liquid Solvents to Mammalian Cells

Molecular oxygen, O2, has long provided a cornerstone for studies in chemistry, physics, and biology. Although the triplet ground state, O2(X3Σg-), has garnered much attention, the lowest excited electronic state, O2(a1Δg), commonly called singlet oxygen, has attracted appreciable interest, principally because of its unique chemical reactivity in systems ranging from the Earth's atmosphere to biological cells. Because O2(a1Δg) can be produced and deactivated in processes that involve light, the photophysics of O2(a1Δg) are equally important. Moreover, pathways for O2(a1Δg) deactivation that regenerate O2(X3Σg-), which address fundamental principles unto themselves, kinetically compete with the chemical reactions of O2(a1Δg) and, thus, have practical significance. Due to technological advances (e.g., lasers, optical detectors, microscopes), data acquired in the past ∼20 years have increased our understanding of O2(a1Δg) photophysics appreciably and facilitated both spatial and temporal control over the behavior of O2(a1Δg). One goal of this Review is to summarize recent developments that have broad ramifications, focusing on systems in which oxygen forms a contact complex with an organic molecule M (e.g., a liquid solvent). An important concept is the role played by the M+•O2-• charge-transfer state in both the formation and deactivation of O2(a1Δg).

The Structure of Carbon Dioxide at the Air-Water Interface and its Chemical Implications

The efficient reduction of CO2 into valuable products is a challenging task in an international context marked by the climate change crisis and the need to move away from fossil fuels. Recently, the use of water microdroplets has emerged as an interesting reaction media where many redox processes which do not occur in conventional solutions take place spontaneously. Indeed, several experimental studies in microdroplets have already been devoted to study the reduction of CO2 with promising results. The increased reactivity in microdroplets is thought to be linked to unique electrostatic solvation effects at the air-water interface. In the present work, we report a theoretical investigation on this issue for CO2 using first-principles molecular dynamics simulations. We show that CO2 is stabilized at the interface, where it can accumulate, and that compared to bulk water solution, its electron capture ability is larger. Our results suggest that reduction of CO2 might be easier in interface-rich systems such as water microdroplets, which is in line with early experimental data and indicate directions for future laboratory studies. The effect of other relevant factors which could play a role in CO2 reduction potential is discussed.

Substantial breakdown of the hydrogen-bonding network, local density inhomogeneities and fluid-liquid structural transitions in supercritical octanol-1: A molecular dynamics investigation

Molecular dynamics simulations have been employed to explore the hydrogen-bonding structure and dynamics in supercritical octanol-1 at a near-critical temperature and up to high densities and pressures. A substantial breakdown of the hydrogen-bonding network when going from ambient-liquid to supercritical conditions is revealed. The fraction of the non-hydrogen bonded molecules significantly increases in supercritical octanol-1, and a substantial decrease in the intermittent hydrogen-bond lifetime is observed. This behavior is also reflected on the maximum local density augmentation, which is comparable to the values obtained for non-polar and non-hydrogen bonded fluids. The existence of a structural transition from an inhomogeneous fluid phase to a soft-liquid one at densities higher than 2.0 ρc is also revealed. At higher densities, a significant change in the reorientational relaxation process is observed, reflected on the significant increase in the ratio of the Legendre reorientational times τ1R/τ2R. The latter becomes much higher than the value predicted by the Debye model of diffusive reorientation and the corresponding ratio for ambient liquid octanol-1. The non-polar tail of octanol-1 under supercritical conditions reorients more slowly in comparison with the polar tail. Interestingly, the opposite behavior is observed for the ambient liquid, further verifying the strong effect of the breakdown of the hydrogen bonding network on the properties of supercritical octanol-1. In accordance with the above-mentioned findings, the static dielectric constant of supercritical octanol-1 is very low even at high densities and pressures, comparable to the values obtained for non-polar and non-hydrogen bonded fluids.

Conformational equilibria in acrolein-CO2: the crucial contribution of n → π* interactions unveiled by rotational spectroscopy

Using gas phase Fourier-transform microwave spectroscopy complemented by theoretical analysis, this study delivers a comprehensive depiction of the physical origin of the 'n → π* interaction' between CO2 and acrolein, one of the most reactive aldehydes. Three distinct isomers of the acrolein-CO2 complex, linked through a C⋯O tetrel bond (or n → π* interaction) and a C-H⋯O hydrogen bond, have been unambiguously identified in the pulsed jet. Relative intensity measurements allowed estimation on the population ratio of the three isomers to be T1/T2/C1 ≈ 25/5/1. Advanced theoretical analyses were employed to elucidate the intricacies of the noncovalent interactions within the examined complex. This study not only sheds light on the molecular underpinnings of n → π* interactions but also paves the way for future exploration in carbon dioxide capture and utilization, leveraging the fundamental principles uncovered in the study of acrolein-carbon dioxide interactions.

Design of Molecularly Imprinted Polymers Using Supercritical Carbon Dioxide Technology

The design and development of affinity polymeric materials through the use of green technology, such as supercritical carbon dioxide (scCO2), is a rapidly evolving field of research with vast applications across diverse areas, including analytical chemistry, pharmaceuticals, biomedicine, energy, food, and environmental remediation. These affinity polymeric materials are specifically engineered to interact with target molecules, demonstrating high affinity and selectivity. The unique properties of scCO2, which present both liquid- and gas-like properties and an accessible critical point, offer an environmentally-friendly and highly efficient technology for the synthesis and processing of polymers. The design and the synthesis of affinity polymeric materials in scCO2 involve several strategies. Commonly, the incorporation of functional groups or ligands into the polymer matrix allows for selective interactions with target compounds. The choice of monomer type, ligands, and synthesis conditions are key parameters of material performance in terms of both affinity and selectivity. In addition, molecular imprinting allied with co-polymerization and surface modification are commonly used in these strategies, enhancing the materials' performance and versatility. This review aims to provide an overview of the key strategies and recent advancements in the design of affinity polymeric materials using scCO2.

Decellularization of cartilage microparticles: Effects of temperature, supercritical carbon dioxide and ultrasound on biochemical, mechanical, and biological properties

One of the approaches to restoring the structure of damaged cartilage tissue is an intra-articular injection of tissue-engineered medical products (TEMPs) consisting of biocompatible matrices loaded with cells. The most interesting are the absorbable matrices from decellularized tissues, provided that the cellular material is completely removed from them with the maximum possible preservation of the structure and composition of the natural extracellular matrix. The present study investigated the mechanical, biochemical, and biological properties of decellularized porcine cartilage microparticles (DCMps) obtained by techniques, differing only in physical treatments, such as freeze-thaw cycling (Protocol 1), supercritical carbon dioxide fluid (Protocol 2) and ultrasound (Protocol 3). Full tissue decellularization was achieved, as confirmed by the histological analysis and DNA quantification, though all the resultant DCMps had reduced glycosaminoglycans (GAGs) and collagen. The elastic modulus of all DCMp samples was also significantly reduced. Most notably, DCMps prepared with Protocol 3 significantly outperformed other samples in viability and the chondroinduction of the human adipose-derived stem cells (hADSCs), with a higher GAG production per DNA content. A positive ECM staining for type II collagen was also detected only in cartilage-like structures based on ultrasound-treated DCMps. The biocompatibility of a xenogenic DCMps obtained with Protocol 3 has been confirmed for a 6-month implantation in the thigh muscle tissue of mature rats (n = 18). Overall, the results showed that the porcine cartilage microparticles decellularized by a combination of detergents, ultrasound and DNase could be a promising source of scaffolds for TEMPs for cartilage reconstruction.

Rotational spectrum of anisole-CO2: Cooperative C···O tetrel bond and CH···O hydrogen bond

Rotational spectrum of the 1:1 anisole-CO₂ complex has been investigated using a pulsed jet Fourier transform microwave spectrometer supplemented with quantum chemical calculations. In the pulsed jet, only one isomer has been observed which is characterized by a dominant C···O tetrel bond and two CH···O_CO2 weak hydrogen bonds. Different theoretical methods predict different orders of relative energies of plausible conformations. The experimental observation is most consistent with the theoretical estimation at the B3LYP-D3(BJ)/6-311++G(d,p) level of theory. Johnson's non-covalent interaction, quantum theory of atoms in molecules and natural bond orbital analyses have been applied to better understand the nature of non-covalent interactions at play in the anisole-CO₂ complex.

Biocompatible and functional properties of a microdispersed tissue-specific 3D matrix from decellularized porcine cartilage

In contrast to decellularization of soft tissues for use as tissue-specific matrices in the creation of tissue-engineered constructs, decellularization of cartilage tissue requires several processing techniques, which can negatively affect the biocompatibility and functional properties of the native extracellular matrix (ECM).

Objective: to study the biocompatible and functional properties of microdispersed tissue-specific 3D matrix from a porcine cartilage that is decellularized by sequential use of chemical, physical and enzymatic techniques.

Materials and methods. For decellularization, microdispersed cartilage particles (MCPs), obtained by cryomilling, were incubated in detergent solutions (sodium dodecyl sulfate and Triton X-100), then treated with supercritical carbon dioxide (scCO2) with 10% ethanol and DNase I. The Ames test (Salmonella typhimurium reverse mutation assay) was used to determine the genotoxicity of decellularized microdispersed cartilage particles (dMCPs). Local and general toxic effects, as well as resorption of dMCPs were studied in vivo on sexually mature outbred rats. Decellularized MCP specimens (10 mg) were implanted into the thigh muscle tissue. Viability of human adipose-derived mesenchymal stem/stromal cells (hAdMSCs), when cultured on dMCPs, was analyzed by in vivo microscopy, stained with fluorescent Calcein AM dye. Cell metabolic activity was assessed using PrestoBlue™ Cell Viability Reagent.

Results. It has been proven that porcine dMCPs implanted in rat muscle after treatment with scCO2 do not exhibit local and general toxic effects, and do not show genotoxicity and negative effects on the reproductive system of animals. After 6 months of in vivo experiment, most (87%) of the implanted decellularized cartilage was resorbed. It was shown that the resulting matrices are able to support adhesion and proliferation of hAdMSCs. Conclusion. Porcine dMCP specimens are suitable for biocompatible medical products in terms of local and general toxic effects, genotoxicity and reproductive toxicity, and can be used as a matrix for creating cell- and tissue-engineered cartilage constructs.

Competitive tetrel bond and hydrogen bond in benzaldehyde-CO2: characterization via rotational spectroscopy

The rotational spectrum of the 1 : 1 benzaldehyde-CO₂ complex has been investigated using pulsed-jet Fourier transform microwave spectroscopy complemented with quantum chemical calculations. Two isomers, both characterized by one C⋯O tetrel bond (n → π* interaction) and one C-H⋯O hydrogen bond (n → σ* interaction), have been observed in the pulsed jet. Competition between the tetrel bond and the hydrogen bond has been disclosed by natural bond orbital analysis: isomer I is characterized by one dominating OC_CO2⋯O tetrel bond (12.6 kJ mol^-1) and a secondary (C-H)_formyl⋯O hydrogen bond (2.2 kJ mol^-1); by contrast, in isomer II the (C-H)_phenyl⋯O hydrogen bond (7.6 kJ mol^-1) becomes the dominant bond, while the OC_CO2⋯O tetrel bond (5.8 kJ mol^-1) becomes much weaker with respect to that of isomer I. Using intensity measurements the relative population ratio of the two isomers was estimated to be N_I/N_II ≈ 2/1.

Ultra-fast single-crystal polymerization of large-sized covalent organic frameworks

In principle, polymerization tends to produce amorphous or poorly crystalline materials. Efficiently producing high-quality single crystals by polymerization in solvent remains as an unsolved issue in chemistry, especially for covalent organic frameworks (COFs) with highly complex structures. To produce μm-sized single crystals, the growth time is prolonged to >15 days, far away from the requirements in practical applications. Here, we find supercritical CO2 (sc-CO2) accelerates single-crystal polymerization by 10,000,000 folds, and produces two-dimensional (2D) COF single crystals with size up to 0.2 mm within 2~5 min. Although it is the fastest single-crystal polymerization, the growth in sc-CO2 leads to not only the largest crystal size of 2D COFs, but also higher quality with improved photoconductivity performance. This work overcomes traditional concept on low efficiency of single-crystal polymerization, and holds great promise for future applications owing to its efficiency, industrial compatibility, environmental friendliness and universality for different crystalline structures and linkage bonds.

Characterizing hydrogen and tetrel bonds in clusters of CO2 with carboxylic acids

The interaction between carbon dioxide and planar carboxylic acids has been investigated through the analysis of the microwave spectrum of the acrylic acid·CO2 complex and quantum chemical modeling of the R-COOH·(CO2)1,16 clusters, where R = H, CH2CH. As regards the 1 : 1 compounds, two species, involving the s-cis and s-trans conformers of acrylic acid were observed. For both of them, a similar bidentate interaction arises between the carbonyl group of CO2 and the carboxylic group of the organic acid, leading to the formation of a planar six-membered ring. The binding energy is estimated to be De ≃ 21 kJ mol-1, 1/3 being the energy contributions of the tetrel to hydrogen bonds, respectively. In the 1 : 16 clusters, the ring arrangement is broken, allowing for the interaction of the acid with several CO2 molecules. The CO2 molecules completely surround formic acid, whereas, in the case of acrylic acid, they tend to avoid the allyl chain.

Influence of structure and solubility of chain transfer agents on the RAFT control of dispersion polymerisation in scCO2

Reversible addition-fragmentation chain transfer (RAFT) dispersion polymerisation of methyl methacrylate (MMA) is performed in supercritical carbon dioxide (scCO2) with 2-(dodecylthiocarbonothioylthio)-2-methylpropionic acid (DDMAT) present as chain transfer agent (CTA) and surprisingly shows good control over PMMA molecular weight. Kinetic studies of the polymerisation in scCO2 also confirm these data. By contrast, only poor control of MMA polymerisation is obtained in toluene solution, as would be expected for this CTA which is better suited for acrylates. In this regard, we select a range of CTAs and use them to determine the parameters that must be considered for good control in dispersion polymerisation in scCO2. A thorough investigation of the nucleation stage during the dispersion polymerisation reveals an unexpected "in situ two-stage" mechanism that strongly determines how the CTA works. Finally, using a novel computational solvation model, we identify a correlation between polymerisation control and degree of solubility of the CTAs. All of this ultimately gives rise to a simple, elegant and counterintuitive guideline to select the best CTA for RAFT dispersion polymerisation in scCO2.

pH-Responsive Nanostructures Based on Surface Active Fatty Acid-Protic Ionic Liquids for Imiquimod Delivery in Skin Cancer Topical Therapy

For topical treatment of skin cancer, the design of pH-responsive nanocarriers able to selectively release the drug in the tumor acidic microenvironment represents a reliable option for targeted delivery. In this context, a series of newly synthesized surface-active fatty acid-protic ionic liquids (FA-PILs), based on tetramethylguanidinium cation and different natural hydrophobic fatty acid carboxylates, have been investigated with the aim of developing a pH-sensitive nanostructured drug delivery system for cutaneous administration in the skin cancer therapy. The capability of FA-PILs to arrange in micelles when combined with each other and with the non-ionic surfactant d-α-Tocopherol polyethylene glycol succinate (vitamin E TPGS) as well as their ability to solubilize imiquimod, an immuno-stimulant drug used for the treatment of skin cancerous lesions, have been demonstrated. The FA-PILs-TPGS mixed micelles showed pH-sensitivity, suggesting that the acidic environment of cancer cells can trigger nanostructures’ swelling and collapse with consequent rapid release of imiquimod and drug cytotoxic potential enhancement. The in vitro permeation/penetration study showed that the micellar formulation produced effective imiquimod concentrations into the skin exposed to acid environment, representing a potential efficacious and selective drug delivery system able to trigger the drug release in the tumor tissues, at lower and less irritating drug concentrations.

Conformational equilibria of pharmaceuticals in supercritical CO2, IR spectroscopy and quantum chemical calculations

In this work we demonstrate a self-consistent effective technique of analyzing the conformational equilibria of active pharmaceutical ingredient (API) molecules dissolved in supercritical carbon dioxide in a wide range of thermodynamic parameters of state. This approach can be useful for pharmaceutics when the crystalline forms of pharmaceuticals with a high purity degree and desirable polymorphism are produced using CO₂-based supercritical fluids technologies. Within this approach we use a combination of quantum chemical calculations and in situ IR spectroscopy. Quantum chemical calculations allow us to perform the initial conformational search and to determine the energy characteristics of the most stable conformers of API and the energy barriers of transitions between them. IR spectroscopy gives the information on the equilibrium of the most stable conformers of pharmaceuticals dissolved in scCO₂ in the thermodynamic parameter range of interest. Finally we validate our approach by applying it to the study of carbamazepine dissolved in scCO₂ being in permanent contact with an excess of crystalline carbamazepine as an example. The conformational search for carbamazepine molecules in scCO₂ was also performed using molecular dynamics simulation for comparison with the results obtained by the technique presented in this paper.

Insight into the Intermolecular Interaction and Free Radical Polymerizability of Methacrylates in Supercritical Carbon Dioxide

High pressure in situ Fourier transfer infrared/near infrared technology (HP FTIR/NIR) along with theoretical calculation of density functional theory (DFT) method was employed. The solvation behaviors and the free radical homopolymerization of methyl methacrylate (MMA), methacrylate acid (MAA), trifluoromethyl methacrylate (MTFMA) and trifluoromethyl methacrylate acid (TFMAA) in scCO2 were systematically investigated. Interestingly, the previously proposed mechanism of intermolecular-interaction dynamically-induced solvation effect (IDISE) of monomer in scCO2 is expected to be well verified/corroborated in view that the predicted solubility order of the monomers in scCO2 via DFT calculation is ideally consistent with that observed via HP FTIR/NIR. It is shown that MMA and MAA can be easily polymerized, while the free radical polymerizability of MTFMA is considerably poor and TFMAA cannot be polymerized via the free radical initiators. The α trifluoromethyl group (-CF3) may effectively enhance the intermolecular hydrogen bonding and restrain the diffusion of the monomer in scCO2. More importantly, the strong electron-withdrawing inductive effect of -CF3 to C=C may distinctly decrease the atomic charge of the carbon atom in the methylene (=CH2). These two factors are believed to be predominantly responsible for the significant decline of the free radical polymerizability of MTFMA and the other alkyl 2-trifluoromethacrylates in scCO2.

The use of process simulation in supercritical fluids applications

Modelling and simulation from micro- to macro-scale are needed to attain a broader commercialization of supercritical technologies.

Structural inhomogeneity as a factor promoting the homogenous catalysis of CO2 hydrogenation by (PMe3)4RuH2

During homogenous catalysis by organometallic complexes, the dissociation of a ligand to produce an unsaturated site on the metal center is often invoked as the first step of activation, especially when photo-excitation is involved. In this theoretical study, we demonstrated that under mild conditions, a thermodynamically unstable yet dynamically favorable active intermediate could be produced by the inhomogeneity of the solvent distribution around the catalyst rather than by ligand dissociation. This occurred at the end of the first catalytic cycle when the product was eliminated. The empty site was immediately filled by one of the additive molecules aggregated around the reaction center even when the intermediate complex was unstable, producing a transient and more active catalyst. This process accounted for the accelerated reaction rate observed in the landmark CO₂ hydrogenation catalyzed by (PMe₃)₄RuH₂ in supercritical CO₂ when H₂O, MeOH, or HNMe₂ was added. This also suggests a new way to exploit the structural inhomogeneity around an organometallic complex for the design of superior catalysts.

Electronic Interactions in Iminophosphorane Superbase Complexes with Carbon Dioxide

Iminophosphoranes or phosphazenes are an important class of compounds with increasing use in synthetic organic chemistry as neutral organic superbases exhibiting low nucleophilicity. Their electronic structure and therefore their properties strongly depend on substitution, but there have been very few theoretical studies devoted to this topic, and more specifically to the formation of electron donor-acceptor complexes of iminophosphoranes with electrophiles. In this work, we have investigated the interaction with carbon dioxide at different ab initio levels. Carbon dioxide usually behaves as a Lewis acid and the reaction with iminiphosphoranes has been described as a nonconventional aza-Wittig process leading to isocyanates. The reaction can be conducted in supercritical CO₂ conditions (carbon dioxide acts as both solvent and reactant), which is a promising strategy in the context of green chemistry. Our calculations have been carried out at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ level for model systems and at the M06-2X/6-611+G(d,p) level for a larger species used in experiments. The electronic interactions and the interaction energies are analyzed and discussed in detail using the natural bond orbital method. Proton affinities and gas-phase basicities are provided as well.