Influence of Ethanol Parametrization on Diffusion Coefficients Using OPLS-AA Force Field.

Molecular dynamics simulations employing the all-atom optimized potential for liquid simulations (OPLS-AA) force field were performed for determining self-diffusion coefficients (D11) of ethanol and tracer diffusion coefficients (D12) of solutes in ethanol at several temperature and pressure conditions. For simulations employing the original OPLS-AA diameter of ethanol's oxygen atom (σOH), calculated and experimental diffusivities of protic solutes differed by more than 25%. To correct this behavior, the σOH was reoptimized using the experimental D12 of quercetin and of gallic acid in liquid ethanol as benchmarks. A substantial improvement of the calculated diffusivities was found by changing σOH from its original value (0.312 nm) to 0.306 nm, with average absolute relative deviations (AARD) of 3.71% and 4.59% for quercetin and gallic acid, respectively. The new σOH value was further tested by computing D12 of ibuprofen and butan-1-ol in liquid ethanol with AARDs of 1.55% and 4.81%, respectively. A significant improvement was also obtained for the D11 of ethanol with AARD = 3.51%. It was also demonstrated that in the case of diffusion coefficients of non-polar solutes in ethanol, the original σOH=0.312 nm should be used for better agreement with experiment. If equilibrium properties such as enthalpy of vaporization and density are estimated, the original diameter should be once again adopted.

Being Positive is Not Everything - Experimental and Computational Studies on the Selectivity of a Self-Assembled, Multiple Redox-State Receptor that Binds Anions with up to Picomolar Affinities.

The interaction of the self-assembled trinuclear ruthenium bowl 13+ , that displays three other accessible oxidation states, with oxo-anions is investigated. Using a combination of NMR and electrochemical experimental data, estimates of the binding affinities of 14+ , 15+ , and 16+ for both halide and oxo-anions were derived. This analysis revealed that, across the range of oxidation states of the host, both high anion binding affinities (>109  M-1 for specific guests bound to 16+ ) and high selectivities (a range of >107  M-1 ) were observed. As the crystal structure of binding of the hexafluorophosphate anion revealed that the host has two potential binding sites (named the α and ß pockets), the host-guest properties of both putative binding sites of the bowl, in all of its four oxidation states, were investigated through detailed quantum-based computational studies. These studies revealed that, due to the interplay of ion-ion interactions, charge-assisted hydrogen-bonding and anion-π interactions, binding to the α pocket is generally preferred, except for the case of the relatively large and lipophilic hexafluorophosphate anionic guest and the host in the highest oxidation states, where the ß pocket becomes relatively favourable. This analysis confirms that host-guest interactions involving structurally complex supramolecular architectures are driven by a combination of non-covalent interactions and, even in the case of charged binding pairs, simple ion-ion interactions alone cannot accurately define these recognition processes.

Solubilities of Amino Acids in Aqueous Solutions of Chloride or Nitrate Salts of Divalent (Mg2+ or Ca2+) Cations.

The solubilities of glycine, l-leucine, l-phenylalanine, and l-aspartic acid were measured in aqueous MgCl2, Mg(NO3)2, CaCl2,, and Ca(NO3)2 solutions with concentrations ranging from 0 to 2 mol/kg at 298.2 K. The isothermal analytical method was used combined with the refractive index measurements for composition analysis guaranteeing good accuracy. All salts induced a salting-in effect with a higher magnitude for those containing the Ca2+ cation. The nitrate anions also showed stronger binding with the amino acids, thus increasing their relative solubility more than the chloride anions. In particular, calcium nitrate induces an increase in the amino acid solubility from 2.4 (glycine) to 4.6 fold (l-aspartic acid) compared to the corresponding value in water. Amino acid solubility data in aqueous MgCl2 and CaCl2 solutions collected from the open literature were combined with that from this work, allowing us to analyze the relations between the amino acid structure and the salting-in magnitude.

Multifunctionality in an Ion-Exchanged Porous Metal-Organic Framework.

Porous robust materials are typically the primary selection of several industrial processes. Many of these compounds are, however, not robust enough to be used as multifunctional materials. This is typically the case of Metal-Organic Frameworks (MOFs) which rarely combine several different excellent functionalities into the same material. In this report we describe the simple acid-base postsynthetic modification of isotypical porous rare-earth-phosphonate MOFs into a truly multifunctional system, maintaining the original porosity features: [Ln(H3pptd)]·xSolvent [where Ln3+ = Y3+ (1) and (Y0.95Eu0.05)3+ (1_Eu)] are converted into [K3Ln(pptd)]·zSolvent [where Ln3+ = Y3+ (1K) and (Y0.95Eu0.05)3+ (1K_Eu)] by immersing the powder of 1 and 1_Eu into an ethanolic solution of KOH for 48 h. The K+-exchanged Eu3+-based material exhibits a considerable boost in CO2 adsorption, capable of being reused for several consecutive cycles. It can further separate C2H2 from CO2 from a complex ternary gas mixture composed of CH4, CO2, and C2H2. This high adsorption selectivity is, additionally, observed for other gaseous mixtures, such as C3H6 and C3H8, with all these results being supported by detailed theoretical calculations. The incorporation of K+ ions notably increases the electrical conductivity by 4 orders of magnitude in high relative humidity conditions. The conductivity is assumed to be predominantly protonic in nature, rendering this material as one of the best conducting MOFs reported to date.

Cinnamic Derivatives as Antitubercular Agents: Characterization by Quantitative Structure-Activity Relationship Studies.

Tuberculosis, caused by Mycobacterium tuberculosis (Mtb), remains one of the top ten causes of death worldwide and the main cause of mortality from a single infectious agent. The upsurge of multi- and extensively-drug resistant tuberculosis cases calls for an urgent need to develop new and more effective antitubercular drugs. As the cinnamoyl scaffold is a privileged and important pharmacophore in medicinal chemistry, some studies were conducted to find novel cinnamic acid derivatives (CAD) potentially active against tuberculosis. In this context, we have engaged in the setting up of a quantitative structure-activity relationships (QSAR) strategy to: (i) derive through multiple linear regression analysis a statistically significant model to describe the antitubercular activity of CAD towards wild-type Mtb; and (ii) identify the most relevant properties with an impact on the antitubercular behavior of those derivatives. The best-found model involved only geometrical and electronic CAD related properties and was successfully challenged through strict internal and external validation procedures. The physicochemical information encoded by the identified descriptors can be used to propose specific structural modifications to design better CAD antitubercular compounds.

Unravelling the Structure of Chemisorbed CO2 Species in Mesoporous Aminosilicas: A Critical Survey.

Chemisorbent materials, based on porous aminosilicas, are among the most promising adsorbents for direct air capture applications, one of the key technologies to mitigate carbon emissions. Herein, a critical survey of all reported chemisorbed CO2 species, which may form in aminosilica surfaces, is performed by revisiting and providing new experimental proofs of assignment of the distinct CO2 species reported thus far in the literature, highlighting controversial assignments regarding the existence of chemisorbed CO2 species still under debate. Models of carbamic acid, alkylammonium carbamate with different conformations and hydrogen bonding arrangements were ascertained using density functional theory (DFT) methods, mainly through the comparison of the experimental 13C and 15N NMR chemical shifts with those obtained computationally. CO2 models with variable number of amines and silanol groups were also evaluated to explain the effect of amine aggregation in CO2 speciation under confinement. In addition, other less commonly studied chemisorbed CO2 species (e.g., alkylammonium bicarbonate, ditethered carbamic acid and silylpropylcarbamate), largely due to the difficulty in obtaining spectroscopic identification for those, have also been investigated in great detail. The existence of either neutral or charged (alkylammonium siloxides) amine groups, prior to CO2 adsorption, is also addressed. This work extends the molecular-level understanding of chemisorbed CO2 species in amine-oxide hybrid surfaces showing the benefit of integrating spectroscopy and theoretical approaches.

Implicit solvent effects in the determination of Brønsted-Evans-Polanyi relationships for heterogeneously catalyzed reactions.

Heterogeneously catalyzed reactions take place at the catalyst surface where, depending on the conditions and process, the reacting molecules are either in the gas or liquid phase. In the latter case, computational heterogeneous catalysis studies usually neglect solvent effects. In this work, we systematically analyze how the electrostatic contribution to solvent effects influences the atomic structure of the reactants and products as well as the adsorption, activation, and reaction energy for the dissociation of water on several planar and stepped transition metal surfaces. The solvent effects were accounted for through an implicit model that describes the effect of electrostatics, cavitation, and dispersion on the interaction between the solute and solvent. The present study shows that the activation energy barriers are only slightly influenced by the inclusion of the electrostatic solvent effects accounted for in a continuum solvent approach, whereas the adsorption energies of the reactants or products are significantly affected. Encouragingly, the linear equations corresponding to the Brønsted-Evans-Polanyi relationships (BEPs), relating the activation energies for the dissociation reaction with a suitable descriptor, e.g. the adsorption energies of the products of the reaction on the difference surfaces, are similar in the presence or in the absence of the solvent. Despite the associated uncertainties, this suggests that BEP relationships derived without the implicit consideration of the solvent are still valid for predicting the activation energy barriers of catalytic reactions from a reaction descriptor.

Mechanisms of phase separation in temperature-responsive acidic aqueous biphasic systems.

The temperature responsive solubility of ionic liquids with 'bulky' polar regions, such as tributyltetradecyl phosphonium chloride ([P44414]Cl), in acidic aqueous solutions is elucidated through a combined experimental and computational approach. The temperature effect in the acidic aqueous biphasic system HCl/[P44414]Cl/H2O was characterised in the range 273 K to 373 K and was found to significantly deviate from the corresponding aqueous biphasic system with NaCl. A new transferable coarse grained MARTINI model for [P44414]Cl was developed, validated and applied to provide a molecular understanding of the experimental results. It is shown that the presence of large aliphatic moieties around the central phosphorus atoms of [P44414]Cl results in a decrease in the electrostatic repulsion between the cationic moieties, leading the [P44414]+ cation to present a behaviour conventionally associated with non-ionic surfactants. This difference in behaviour between HCl and NaCl was shown to result from the greater interaction of the hydronium cation with the micelle surface, thereby enhancing the [P44414]Cl aggregation.

Hydrogen bonding networks in gabapentin protic pharmaceutical salts: NMR and in silico studies.

Hydrogen bonds (HBs) play a key role in the supramolecular arrangement of crystalline solids and, although they have been extensively studied, the influence of their strength and geometry on crystal packing remains poorly understood. Here we describe the crystal structures of two novel protic gabapentin (GBP) pharmaceutical salts prepared with the coformers methanesulfonic acid (GBP:METHA) and ethanesulfonic acid (GBP:ETHA). This study encompasses experimental and computational electronic structure analyses of 1 H NMR chemical shifts (CSs), upon in silico HB cleavage. GBP:METHA and GBP:ETHA crystal packing comprise two main structural domains: an ionic layer (characterized by the presence of charge-assisted + NHGBP ⋯O-METHA/ETHA HB interactions) and a neutral layer generated in a different way for each salt, mainly due to the presence of bifurcated HB interactions. A comprehensive study of HB networks is presented for GBP:METHA, by isolating molecular fragments involved in distinct HB types (NH⋯O, OH⋯O, and CH⋯O) obtained from in silico disassembling of an optimized three-dimensional packing structure. Formation of HB leads to calculated 1 H NMR CS changes from 0.4 to ~5.8 ppm. This study further attempts to assess how 1 H NMR CS of protons engaged in certain HB are affected when other nearby HB, involving bifurcated or geminal/vicinal hydrogen atoms, are removed.

Flue gas adsorption on periodic mesoporous phenylene-silica: a DFT approach.

Periodic mesoporous organosilicas (PMOs) were suggested as potential adsorbents for CO2/CH4 separation because of their large affinities towards CO2 and low interaction with CH4. Herewith, we present a comprehensive computational study on the binding properties of flue gas species with the pore walls of periodic mesoporous phenylene-silica (Ph-PMO) for understanding the possible impact of other gaseous species in the CO2/CH4 separation. The calculations considered three exchange-correlation functionals (PBE, PBE-D2 and M06-2X) based on the density functional theory and the walls of the periodic mesoporous phenylene-silica were modelled within the cluster model approach. The components of the flue gas considered were the diatomic CO, H2, N2, O2 and NO molecules, the triatomic CO2, H2O, H2S and SO2 species, the tetratomic SO3 and NH3 gases and the pentatomic CH4 molecule. The calculated data demonstrate that the presence of H2O, SO2, NH3, H2S and SO3 is a significant threat to CO2 capture by Ph-PMO and suggest that the Ph-PMO material would present high selectivity for CO2 over CH4, CO, H2 or N2 adsorption. The adsorption behaviour of flue gas components in Ph-PMO can be directly related to the experimental proton affinities, basicities or even the polarizabilities of the gaseous molecules.

Mechanism of ionic-liquid-based acidic aqueous biphasic system formation.

Ionic-liquid-based acidic aqueous biphasic systems (IL-based AcABS) represent a promising alternative to the solvent extraction process for the recovery of critical metals, in which the substitution of the inorganic salt by an acid allows for a 'one-pot' approach to the leaching and separation of metals. However, a more fundamental understanding of AcABS formation remains wanting. In this work, the formation mechanisms of AcABS are elucidated through a comparison with traditional aqueous biphasic systems (ABS). A large screening of AcABS formation with a wide range of IL identifies the charge shielding of the cation as the primary structural driver for the applicability of an IL in AcABS. Through a systematic study of tributyltetradecylphosphonium chloride ([P44414]Cl) with various chloride salts and acids, we observed the first significant deviation to the cationic Hofmeister series reported for IL-based ABS. Furthermore, the weaker than expected salting-out ability of H3O+ compared to Na+ is attributed to the greater interaction of H3O+ with the [P44414]+ micelle surface. Finally, the remarkable thermomorphic properties of [P44414]Cl based systems are investigated with a significant increase in the biphasic region induced by the increase in the temperature from 298 K to 323 K. These finding allows for the extension of ABS to new acidic systems and highlights their versatility and tunability.

Structure of Chemisorbed CO(2) Species in Amine-Functionalized Mesoporous Silicas Studied by Solid-State NMR and Computer Modeling.

Two-dimensional (2D) solid-state nuclear magnetic resonance (SSNMR) experiments on samples loaded with 13C-labeled CO2, "under controlled partial pressures", have been performed in this work, revealing unprecedented structural details about the formation of CO2 adducts from its reaction with various amine-functionalized SBA-15 containing amines having distinct steric hindrances (e.g., primary, secondary) and similar loadings. Three chemisorbed CO2 species were identified by NMR from distinct carbonyl environments resonating at δC ≈ 153, 160, and 164 ppm. The newly reported chemisorbed CO2 species at δC ≈ 153 ppm was found to be extremely moisture dependent. A comprehensive 1H-based SSNMR study [1D 1H and 2D 1H-X heteronuclear correlation (HETCOR, X = 13C, 29Si) experiments] was performed on samples subjected to different treatments. It was found that all chemisorbed CO2 species are involved in hydrogen bonds (HBs) with either surface silanols or neighboring alkylamines. 1H chemical shifts up to 11.8 ppm revealed that certain chemisorbed CO2 species are engaged in very strong HBs. We effectively demonstrate that NMR may help in discriminating among free and hydrogen-bonded functional groups. 13C{14N} dipolar-recoupling NMR showed that the formation of carbonate or bicarbonate is excluded. Density functional theory calculations on models of alkylamines grafted into the silica surface assisted the 1H/13C assignments and validated various HB arrangements that may occur upon formation of carbamic acid. This work extends the understanding of the chemisorbed CO2 structures that are formed upon bonding of CO2 with surface amines and readily released from the surface by pressure swing.

Adsorption of CO on the rutile TiO2(110) surface: a dispersion-corrected density functional theory study.

The geometry, energy and stretching frequency of carbon monoxide on the rutile TiO2(110) surface for coverages between 0.125 and 1.5 ML are investigated by means of density functional theory calculations. Four different approaches were considered, namely, the PBE exchange-correlation functional and the PBE-D2, vdW-DF and vdW-DF2 methods incorporating van der Waals dispersion interactions of different theoretical complexity and empiricism. It is found that upon the increase of the surface coverage, the adsorption becomes less favorable due to lateral destabilizing interactions between adsorbed molecules. The preferred geometry for CO changes from an upright configuration at 0.125 ML to tilted configurations at 1.5 ML and the tilting of the C-O axis from the surface normal increases with the increase of the surface coverage. At 1 ML, all computational approaches predict alternate tilted configurations which contradict the interpretation of recent experimental infrared reflection-absorption spectroscopic findings suggesting upright CO geometries. Encouragingly, a very good agreement between calculated and experimental shifts of the C-O stretching frequency of adsorbed CO at different coverages with respect to gaseous CO species was reached.

Prediction of metallic nanotube reactivity for H2O activation.

The reactivity of metallic nanotubes toward the catalysis of water dissociation, a key step in the water gas shift reaction (WGSR), was analyzed through density functional theory (DFT) calculations. Water dissociation was studied on surfaces of nanotubes based on copper, gold and platinum, and also on platinum doped copper and gold nanotubes. Gold and copper nanotubes present activities that are similar to those of their corresponding extended surfaces but, in the case of the Pt(5,3) nanotube, a significant improvement in the activity is found when compared with the extended surfaces. In fact, the calculations predict the water dissociation to be spontaneous on Pt(5,3) with a low activation energy barrier. The platinum doping of gold and copper nanotubes leads to contrasting effects, i.e., with a slight increase of activity found on gold and a slight decrease of activity in the case of copper. The consideration of a Brönsted-Evans-Polanyi (BEP) relationship to estimate the activation energy barriers for the O-H bond break leads to a satisfactory agreement between estimated and explicitly calculated values which suggests the validity of the BEP relationship for qualitative predictions of the activities of metal nanotubes towards the water dissociation reaction.

Corrosion inhibition of copper in aqueous chloride solution by 1H-1,2,3-triazole and 1,2,4-triazole and their combinations: electrochemical, Raman and theoretical studies.

Triazoles are well-known organic corrosion inhibitors of copper. 1H-1,2,3-Triazole and 1,2,4-triazole, two very simple molecules with the only difference being the positions of the nitrogen atoms in the triazole ring, were studied in this work as corrosion inhibitors of copper in 50 mM NaCl solution using a set of electrochemical and analytical techniques. The results of electrochemical tests indicate that 1H-1,2,3-triazole exhibited superior inhibitor properties but could not suppress anodic copper dissolution at moderate anodic potentials (>+300 mV SCE), while 1,2,4-triazole, although it exhibited higher anodic currents, suppressed anodic copper dissolution at very anodic potentials. Density functional theory calculations were also performed to interpret the measured data and trends observed in the electrochemical studies. The computational studies considered either the inhibitors isolated in the gaseous phase or adsorbed onto Cu(111) surface models. From the calculations, the mechanisms of the inhibitive effects of both triazoles were established and plausible mechanisms of formation of the protective films on the Cu surface were proposed. The results of this study hold positive implications for research in the areas of catalysis, and copper content control in water purification systems.

Why are some cyano-based ionic liquids better glucose solvents than water?

Among different classes of ionic liquids (ILs), those with cyano-based anions have been of special interest due to their low viscosity and enhanced solvation ability for a large variety of compounds. Experimental results from this work reveal that the solubility of glucose in some of these ionic liquids may be higher than in water - a well-known solvent with enhanced capacity to dissolve mono- and disaccharides. This raises questions on the ability of cyano groups to establish strong hydrogen bonds with carbohydrates and on the optimal number of cyano groups at the IL anion that maximizes the solubility of glucose. In addition to experimental solubility data, these questions are addressed in this study using a combination of density functional theory (DFT) and molecular dynamics (MD) simulations. Through the calculation of the number of hydrogen bonds, coordination numbers, energies of interaction and radial and spatial distribution functions, it was possible to explain the experimental results and to show that the ability to favorably interact with glucose is driven by the polarity of each IL anion, with the optimal anion being dicyanamide.

Understanding the reactivity of metallic nanoparticles: beyond the extended surface model for catalysis.

Metallic nanoparticles (NPs) constitute a new class of chemical objects which are used in different fields as diverse as plasmonics, optics, catalysis, or biochemistry. The atomic structure of the NP and its size usually determine the chemical reactivity but this is often masked by the presence of capping agents, solvents, or supports. The knowledge of the structure and reactivity of isolated NPs is a requirement when aiming at designing NPs with a well-defined chemistry. Theoretical models together with efficient computational chemistry algorithms and parallel computer codes offer the opportunity to explore the chemistry of these interesting objects and to understand the effects of parameters such as size, shape and composition allowing one to derive some general trends.

Molecular simulation of the adsorption of methane in Engelhard titanosilicate frameworks.

Molecular simulations were carried out to elucidate the influence of structural heterogeneity and of the presence of extra-framework cations and water molecules on the adsorption of methane in Engelhard titanosilicates, ETS-10 and ETS-4. The simulations employed three different modeling approaches, (i) with fixed cations and water at their single crystal positions, (ii) with fixed cations and water at their optimized positions, and (iii) with mobile extra-framework cations and water molecules. Simulations employing the final two approaches provided a more realistic description of adsorption in these materials, and showed that at least some cations and water molecules are displaced from the crystallographic positions obtained from single crystal data. Upon methane adsorption in the case of ETS-10, the cations move to the large rings, while in the case of ETS-4, the water molecules and cations migrate to more available space in the larger 12-membered ring channels for better accommodation of the methane molecules. For ETS-4, we also considered adsorption in all possible pure polymorph structures and then combined these to provide an estimate of adsorption in a real ETS-4 sample. By comparing simulated adsorption isotherms to experimental data, we were able to show that both the mobility of extra-framework species and the structural heterogeneity should be taken into account for realistic predictions of adsorption in titanosilicate materials.

"Washing-out" ionic liquids from polyethylene glycol to form aqueous biphasic systems.

The molecular-level mechanisms behind the formation of aqueous biphasic systems (ABS) composed of ionic liquids (ILs) and polymers are hitherto not completely understood. For the first time, it is herein shown that polymer-IL-based ABS are a result of a "washing-out" phenomenon, and not of a salting-out effect of the IL over the polymer as assumed in the past few years. Novel evidence is herein provided by experimental results combined with molecular dynamics (MD) simulations and density functional theory (DFT) calculations.

Density functional theory study of the water dissociation on platinum surfaces: general trends.

We report a comparative periodic density functional theory study of the reaction of water dissociation on five platinum surfaces, e.g., Pt(111) Pt(100), Pt(110), Pt(211), and Pt(321). These surfaces were chosen to study the surface structural effects in the reaction of water dissociation. It was found that water molecules adsorb stronger on surfaces presenting low coordinated atoms in the surface. In the cases of the stepped Pt(110) and kinked Pt(321) surfaces, the activation energy barriers are smaller than the adsorption energies for the water molecule on the corresponding surfaces. Therefore, the calculations suggest that the dissociation reaction will take place preferentially at corner or edge sites on platinum particles with the (110) orientation. The inclusion of the results obtained in this work in previous derived BEP relationships confirms that the adsorption energy of the reaction products arises as the most appropriate descriptor for water dissociation on transition metal surfaces.