How to Predict the pK a of Any Compound in Any Solvent

Molecular Electrochemical Catalysis of CO-to-Formaldehyde Conversion with a Cobalt Complex

Formox, a highly energy-intensive process, currently serves as the primary source of formaldehyde (HCHO), for which there is a crucial and steadily growing chemical demand. The alternative electrochemical production of HCHO from C1 carbon sources such as CO2 and CO is still in its early stages, with even the few identified cases lacking mechanistic rationalization. In this study, we demonstrate that cobalt phthalocyanine (CoPc) immobilized on multiwalled carbon nanotubes (MW-CNTs) constitutes an excellent electrocatalytic system for producing HCHO with productivity through the direct reduction of CO, the two-electron reduction product of CO2. By carefully adjusting both the pH and the applied potential, we identified conditions that enable the production of HCHO with a partial current density of 0.64 mA cm-2 (17.5% Faradaic efficiency, FE) and a total FE of 61.2% for the liquid products (formaldehyde and methanol). A reduction mechanism is proposed.

Exploration of Vitamin B6-Based Redox-Active Pyridinium Salts Towards the Application in Aqueous Organic Flow Batteries

Pyridoxal hydrochloride, a vitamin B6 vitamer, was synthetically converted to a series of diverse redox-active benzoyl pyridinium salts. Cyclic voltammetry studies demonstrated redox reversibility under basic conditions, and two of the most promising salts were subjected to laboratory-scale flow battery tests involving galvanostatic cycling at 10 mM in 0.1 M NaOH. In these tests, the battery was charged completely, corresponding to the transfer of two electrons to the electrolyte, but no discharge was observed. Both CV analysis and electrochemical simulations confirmed that the redox wave observed in the experimental voltammograms corresponds to a two-electron process. To explain the irreversibility in the battery tests, we conducted bulk electrolysis with the benzoyl pyridinium salts, affording the corresponding benzylic secondary alcohols. Computational studies suggest that the reduction proceeds in three consecutive steps: first electron transfer (ET), then proton-coupled electron transfer (PCET) and finally proton transfer (PT) to give the secondary alcohol. 1H NMR deuterium exchange studies indicated that the last PT step is not reversible in 0.1 M NaOH, rendering the entire redox process irreversible. The apparent reversibility observed in CV at the basic media likely arises from the slow rate of the PT step at the timescale of the measurement.

Increasing the polarity of β-lapachone does not affect its binding capacity with bovine plasma protein

Despite ortho-quinones showing several biological and pharmacological activities, there is still a lack of biophysical characterization of their interaction with albumin - the main carrier of different endogenous and exogenous compounds in the bloodstream. Thus, the interactive profile between bovine serum albumin (BSA) with β-lapachone (1) and its corresponding synthetic 3-sulfonic acid (2, under physiological pH in the sulphonate form) was performed. There is one main binding site of albumin for both β-lapachones (n ≈ 1) and a static fluorescence quenching mechanism was proposed. The Stern-Volmer constant (KSV) values are 104 M-1, indicating a moderate binding affinity. The enthalpy (-3.41 ± 0.45 and - 8.47 ± 0.37 kJ mol-1, for BSA:1 and BSA:2, respectively) and the corresponding entropy (0.0707 ± 0.0015 and 0.0542 ± 0.0012 kJ mol-1 K-1) values indicate an enthalpically and entropically binding driven. Hydrophobic interactions and hydrogen bonding are the main binding forces. The differences in the polarity of 1 and 2 did not change significantly the affinity to albumin. In addition, the 1,2-naphthoquinones showed a similar binding trend compared with 1,4-naphthoquinones.

Accurate Computation of Aqueous pKas of Biologically Relevant Organic Acids: Overcoming the Challenges Posed by Multiple Conformers, Tautomeric Equilibria, and Disparate Functional Groups with the Fully Black-Box pK-Yay Method

The use of static electronic structure calculations to compute solution-phase pKas offers a great advantage in that a macroscopic bulk property could be computed via microscopic computations involving very few molecules. There are various sources of errors in the quantum chemical calculations though. Overcoming these errors to accurately compute pKas of a plethora of acids is an active area of research in physical chemistry pursued by both computational as well as experimental chemists. We recently developed the pK-Yay method in our attempt to accurately compute aqueous pKas of strong and weak acids. The method is fully black-box, computationally inexpensive, and is very easy for even a nonexpert to use. However, the method was thus far tested on very few molecules (only 16 in all). Herein, in order to assess the future applicability of pK-Yay, we study the effect of multiple conformers, the presence of tautomers under equilibrium, and the impact of a wide variety of functional groups (derivatives of acetic acid with substituents at various positions, dicarboxylic acids, aromatic carboxylic acids, amines and amides, phenols and thiols, and fluorine bearing organic acids). Starting with more than 1000 conformers and tautomers, this study establishes that overall errors of ∼ 1.0 pKa units are routinely obtained for a majority of the molecules. Larger errors are noted in cases where multiple charges, intramolecular hydrogen bonding, and several ionizable functional groups are simultaneously present. An important conclusion to emerge from this work is that, the computed pKas are insensitive (difference <0.5) to whether we consider multiple conformers/tautomers or only choose the most stable conformer/tautomer. Further, pK-Yay captures the stereoelectronic effects arising due to differing axial vs equatorial pattern, and is useful to predict the dominant acid-base equilibrium in a system featuring several equilibria. Overall, pK-Yay may be employed in several chemical applications featuring organic molecules and biomonomers.

Measurement of the pKa Values of Organic Molecules in Aqueous-Organic Solvent Mixtures by 1H NMR without External Calibrants

Aqueous-organic solvent mixtures are commonly used for reactions or analyses, where the components of a system are insoluble in pure water. The acid dissociation constant is an important property to measure in these media as it determines the charge state, solubility, and reactivity of a molecule. While NMR spectroscopy is an established tool for the measurement of pKa in water, its use in aqueous-organic solvents is greatly hindered by the requirement for external calibrants on which a working pH scale is set. Such calibrants include buffer solutions, "anchor" molecules with known pKa values, and pH electrodes that have undergone lengthy calibration procedures in the solvent mixture of interest. However, such calibrations are often inconvenient to perform, while literature pKa data covering the required range may not be available at the solvent composition or the temperature of interest. Here, we present a method to determine pKa in aqueous-organic solvents directly by NMR. We first determine pKa of an organic acid such as 2,6-dihydroxybenzoic acid (2,6-DHB) by measuring its 1H chemical shift as a function of concentration along a concentration gradient using chemical shift imaging (CSI). Using 2,6-DHB as a reference, we then determine pKa of less acidic molecules in single CSI experiments via the variation of their 1H chemical shifts along pH gradients. As proof of concept, we determine the pKa values of organic acids and bases up to pKa 10 in 50% (v/v) 1-propanol/water, 50% (v/v) dimethyl sulfoxide/water, and 30% (v/v) acetonitrile/water and obtain good agreement with the literature values.

Density functional theory and machine learning for electrochemical square-scheme prediction: an application to quinone-type molecules relevant to redox flow batteries

Proton-electron transfer (PET) reactions are rather common in chemistry and crucial in energy storage applications. How electrons and protons are involved or which mechanism dominates is strongly molecule and pH dependent. Quantum chemical methods can be used to assess redox potential (Ered.) and acidity constant (pKa) values but the computations are rather time consuming. In this work, supervised machine learning (ML) models are used to predict PET reactions and analyze molecular space. The data for ML have been created by density functional theory (DFT) calculations. Random forest regression models are trained and tested on a dataset that we created. The dataset contains more than 8200 quinone-type organic molecules that each underwent two proton and two electron transfer reactions. Both structural and chemical descriptors are used. The HOMO of the reactant and LUMO of the product participating in the oxidation reaction appeared to be strongly associated with Ered.. Trained models using a SMILES-based structural descriptor can efficiently predict the pKa and Ered. with a mean absolute error of less than 1 and 66 mV, respectively. Good prediction accuracy of R2 > 0.76 and >0.90 was also obtained on the external test set for Ered. and pKa, respectively. This hybrid DFT-ML study can be applied to speed up the screening of quinone-type molecules for energy storage and other applications.

Calculating the Aqueous pKa of Phenols: Predictions for Antioxidants and Cannabinoids

We aim to develop a theoretical methodology for the accurate aqueous pK_a prediction of structurally complex phenolic antioxidants and cannabinoids. In this study, five functionals (M06-2X, B3LYP, BHandHLYP, PBE0, and TPSS) and two solvent models (SMD and PCM) were combined with the 6-311++G(d,p) basis set to predict pK_a values for twenty structurally simple phenols. None of the direct calculations produced good results. However, the correlations between the calculated Gibbs energy difference of each acid and its conjugate base, ΔGaq(BA)°=ΔGaqA-°-ΔGaq(HA)°, and the experimental aqueous pK_a values had superior predictive accuracy, which was also tested relative to an independent set of ten molecules of which six were structurally complex phenols. New correlations were built with twenty-seven phenols (including the phenols with experimental pK_a values from the test set), which were used to make predictions. The best correlation equations used the PCM method and produced mean absolute errors of 0.26-0.27 pK_a units and R² values of 0.957-0.960. The average range of predictions for the potential antioxidants (cannabinoids) was 0.15 (0.25) pK_a units, which indicates good agreement between our methodologies. The new correlation equations could be used to make pK_a predictions for other phenols in water and potentially in other solvents where they might be more soluble.

Temperature-Controlled Syngas Production via Electrochemical CO2 Reduction on a CoTPP/MWCNT Composite in a Flow Cell

The mixture of CO and H₂, known as syngas, is a building block for many substantial chemicals and fuels. Electrochemical reduction of CO₂ and H₂O to syngas would be a promising alternative approach for its synthesis due to negative carbon emission footprint when using renewable energy to power the reaction. Herein, we present temperature-controlled syngas production by electrochemical CO₂ and H₂O reduction on a cobalt tetraphenylporphyrin/multiwalled carbon nanotube (CoTPP/MWCNT) composite in a flow cell in the temperature range of 20-50 °C. The experimental results show that for all the applied potentials the ratio of H₂/CO increases with increasing temperature. Interestingly, at -0.6 V _RHE and 40 °C, the H₂/CO ratio reaches a value of 1.2 which is essential for the synthesis of oxo-alcohols. In addition, at -1.0 V _RHE and 20 °C, the composite shows very high selectivity toward CO formation, reaching a Faradaic efficiency of ca. 98%. This high selectivity of CO formation is investigated by density functional theory modeling which underlines that the potential-induced oxidation states of the CoTPP catalyst play a vital role in the high selectivity of CO production. Furthermore, the stability of the formed intermediate species is evaluated in terms of the pK_a value for further reactions. These experimental and theoretical findings would provide an alternative way for syngas production and help us to understand the mechanism of molecular catalysts in dynamic conditions.

Stability of Dibromo-Dipyrromethene Complexes Coordinated with B, Zn, and Cd in Solutions of Various Acidities

The spectral luminescent properties of dipyrromethenates halogenated with bromine on both ends of the long axis and coordinated using boron fluoride, zinc, or cadmium in neutral ethanol and acidified with hydrochloric acid solutions were studied. The constants of the acid-base equilibrium of the complexes in the proton-donor solvents in the ground and excited states was determined. The mechanisms of complex protonation were discussed, depending on the structure of the compounds. The electronic structures of the neutral and protonated compounds were modeled and analyzed based on the quantum-chemical method. The structures and spectral-luminescence properties were calculated using the SMD model of ethanol solvent using the TD-DFT theory with the B3LYP functional and the composite def2-SVP/def2-TZVP/def2-TZVPP_ECP basis sets, depending on the atomic number of the elements.

An Accurate Approach for Computational pKa Determination of Phenolic Compounds

Computational chemistry is a valuable tool, as it allows for in silico prediction of key parameters of novel compounds, such as pKa. In the framework of computational pKa determination, the literature offers several approaches based on different level of theories, functionals and continuum solvation models. However, correction factors are often used to provide reliable models that adequately predict pKa. In this work, an accurate protocol based on a direct approach is proposed for computing phenols pKa. Importantly, this methodology does not require the use of correction factors or mathematical fitting, making it highly practical, easy to use and fast. Above all, DFT calculations performed in the presence two explicit water molecules using CAM-B3LYP functional with 6-311G+dp basis set and a solvation model based on density (SMD) led to accurate pKa values. In particular, calculations performed on a series of 13 differently substituted phenols provided reliable results, with a mean absolute error of 0.3. Furthermore, the model achieves accurate results with -CN and -NO₂ substituents, which are usually excluded from computational pKa studies, enabling easy and reliable pKa determination in a wide range of phenols.

Universal Trends between Acid Dissociation Constants in Protic and Aprotic Solvents

pKa values in non-aqueous solvents are of critical importance in many areas of chemistry. Our knowledge is, despite their relevance, still limited to the most fundamental properties and few pKa values in the most common solvents. Taking advantage of a recently introduced computationally efficient procedure we computed the pKa values of 182 compounds in 21 solvents. This data set is used to establish for the first time universal trends between all solvents. Our computations indicate, that the total charge of the molecule and the charge of the acidic group combined with the Kamlet-Taft solvatochromic parameters are sufficient to predict pKa values with at least semi- quantitative accuracy. We find, that neutral acids such as alcohols are strongly affected by the solvent properties. This is contrasted by cationic acids like ammonium ions whose pKa is often almost completely independent from the choice of solvent.

Modeling pKa of the Brønsted Bases as an Approach to the Gibbs Energy of the Proton in Acetonitrile

A simple but efficient computational approach to calculate pK_a in acetonitrile for a set of phosphorus, nitrogen, and carbon bases was established. A linear function that describes relations between the calculated ΔG’_a.sol(BH⁺) and pK_a values was determined for each group of bases. The best model was obtained through the variations in the basis set, in the level of theory (density functionals or MP2), and in the continuum solvation model (IPCM, CPCM, or SMD). The combination of the IPCM/B3LYP/6-311+G(d,p) solvation approach with MP2/6-311+G(2df,p)//B3LYP/6-31G(d) gas-phase energies provided very good results for all three groups of bases with R² values close to or above 0.99. Interestingly, the slopes and the intercepts of the obtained linear functions showed significant deviations from the theoretical values. We made a linear plot utilizing all the conducted calculations and all the structural variations and employed methods to prove the systematic nature of the intercept/slope dependence. The interpolation of the intercept to the ideal slope value enabled us to determine the Gibbs energy of the proton in acetonitrile, which amounted to −258.8 kcal mol⁻¹. The obtained value was in excellent agreement with previously published results.