Ab Initio Study of the Stabilization of Multiply Charged Anions in Water

Unraveling the mystery of solvation-dependent fluorescence of fluorescein dianion using computational study

Fluorescein, one of the brightest fluorescent dye molecules, is a widely used fluorophore for various applications from biomedicine to industry. The dianionic form of fluorescein is responsible for its high fluorescence quantum yield. Interestingly, the molecule was found to be nonfluorescent in the gas phase. This characteristic is attributed to the photodetachment process, which out-competes the fluorescence emission in the gas phase. In this work, we show that the calculated vertical and adiabatic detachment energies of fluorescein dianion in the gas and solvent phases account for the drastic differences observed in their fluorescence characteristics. The functional dependence of these detachment energies on the dianion's microsolvation was systematically investigated. The performance of different solvent models was also assessed. The higher thermodynamic stability of fluorescein dianion over the monoanion doublet in the solvent phase plays a crucial role in quenching photodetachment and activating the radiative channel with a high fluorescence quantum yield.

Enriching Metal-Oxygen Species and Phosphate Modulating of Active Sites for Robust Electrocatalytical CO2 Reduction

Direct electrochemical reduction of CO2 (CO2 RR) into value-added chemicals is a promising solution to reduce carbon emissions. The activity of CO2 RR is influenced deeply by the reaction microenvironment and electronic properties of the catalysts. Herein, the surface PO4 3- anions are tuned to modulate the local microenvironment and the electronic properties of the indium-based catalyst with abundant metal-oxygen species enabling efficient electrochemical conversion of CO2 to HCOO- . Indium nanoparticles coupled with PO4 3- anions (PO4 3- -In NPs) achieve a high selectivity of HCOO- up to 91.4% at a low potential of -0.98 V versus reversible hydrogen electrode (versus RHE) and a high HCOO- partial current density of 279.3 mA cm-2 at -1.1 V versus RHE in the electrochemical flow cell. In situ and ex situ characterizations confirm the PO4 3- anions keep stable on the surface of indium during CO2 RR, accelerating the generation of OCHO* intermediate. From density functional theory calculations, PO4 3- anions enrich the metal-oxygen species on the substrate to optimize the electronic structure of the catalysts and induce a local microenvironment with massive K+ ions on the interface, thus reducing the activation energy barrier of CO2 RR.

Time-resolved terahertz-Raman spectroscopy reveals that cations and anions distinctly modify intermolecular interactions of water

The solvation of ions changes the physical, chemical and thermodynamic properties of water, and the microscopic origin of this behaviour is believed to be ion-induced perturbation of water's hydrogen-bonding network. Here we provide microscopic insights into this process by monitoring the dissipation of energy in salt solutions using time-resolved terahertz-Raman spectroscopy. We resonantly drive the low-frequency rotational dynamics of water molecules using intense terahertz pulses and probe the Raman response of their intermolecular translational motions. We find that the intermolecular rotational-to-translational energy transfer is enhanced by highly charged cations and is drastically reduced by highly charged anions, scaling with the ion surface charge density and ion concentration. Our molecular dynamics simulations reveal that the water-water hydrogen-bond strength between the first and second solvation shells of cations increases, while it decreases around anions. The opposite effects of cations and anions on the intermolecular interactions of water resemble the effects of ions on the stabilization and denaturation of proteins.

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 CO₂ species, which may form in aminosilica surfaces, is performed by revisiting and providing new experimental proofs of assignment of the distinct CO₂ species reported thus far in the literature, highlighting controversial assignments regarding the existence of chemisorbed CO₂ 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 ¹³C and ¹⁵N NMR chemical shifts with those obtained computationally. CO₂ models with variable number of amines and silanol groups were also evaluated to explain the effect of amine aggregation in CO₂ speciation under confinement. In addition, other less commonly studied chemisorbed CO₂ 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 CO₂ adsorption, is also addressed. This work extends the molecular-level understanding of chemisorbed CO₂ species in amine-oxide hybrid surfaces showing the benefit of integrating spectroscopy and theoretical approaches.

Ab Initio Studies of Calcium Carbonate Hydration

Ab initio simulations of large hydrated calcium carbonate clusters are challenging due to the existence of multiple local energy minima. Extensive conformational searches around hydrated calcium carbonate clusters (CaCO3·nH2O for n = 1-18) were performed to find low-energy hydration structures using an efficient combination of Monte Carlo searches, density-functional tight binding (DFTB+) method, and density-functional theory (DFT) at the B3LYP level, or Møller-Plesset perturbation theory at the MP2 level. This multilevel optimization yields several low-energy structures for hydrated calcium carbonate. Structural and energetics analysis of the hydration of these clusters revealed a first hydration shell composed of 12 water molecules. Bond-length and charge densities were also determined for different cluster sizes. The solvation of calcium carbonate in bulk water was investigated by placing the explicitly solvated CaCO3·nH2O clusters in a polarizable continuum model (PCM). The findings of this study provide new insights into the energetics and structure of hydrated calcium carbonate and contribute to the understanding of mechanisms where calcium carbonate formation or dissolution is of relevance.

ACIDITY AND BASICITY OF ALKOXIDES: ATTAINMENT OF HYDRATION CLUSTERS USING THE AGOA METHODOLOGY

The optimized geometries of the methoxy, ethoxy, iso-propoxy and tert-butoxy anions were calculated using the B3LYP/6-311++G(d, p) theoretical level. Thereby, these structures were used to explore their hydration clusters by means of the AGOA methodology. Analyzing the molecular electrostatic potential of the solute and using the TIP4P model to locate the water molecules over the alkoxide anions, it was observed that the energies of the hydration clusters decrease in order from methoxy to the tert-butoxy. These results corroborate the experimental acidity of the correspondent conjugate alcohols, indicating that the AGOA methodology can be used as an efficient computational methodology to appraise the structures of molecules in aqueous media.

On the stabilization of the carbonate dianion by sulfur dioxide

The stabilization in the gas phase of the carbonate dianion [Formula: see text] by SO2 molecules is being investigated. The geometries of various isomers of [Formula: see text] (SO2)n and [Formula: see text] (SO2)n, for n = 1–4, have been optimized by the B3PW91/6−311+G(3df) method. Single-point CCSD and CCSD(T) energies at the DFT-optimized geometries were obtained for n = 1–3, using the 6−311+G(d) basis set. For n = 1 and 2, the monoanionic clusters are adiabatically more stable than the dianionic ones. However, starting at n = 3, they become less stable. The CCSD adiabatic electron detachment energy of the dianionic cluster switches from −0.39 eV for n = 2 to +0.20 eV for n = 3. The vertical electron detachment energy turns positive at n = 2, with a CCSD value of 1.35 eV. Several of the less stable dianionic, and most of the monoionic clusters, are characterized by the transfer of an oxygen atom from CO3 to SO2, forming [Formula: see text] or [Formula: see text] units, owing to [Formula: see text] + CO2 being more stable than [Formula: see text] + SO2. For the stabilization of the sulfate dianion by stepwise hydration, studied both experimentally and theoretically by other groups, a minimum of three water molecules was required.

Molecular dynamics study of the solvation of calcium carbonate in water

We performed molecular dynamics simulations of diluted solutions of calcium carbonate in water. To this end, we combined and tested previous polarizable models. The carbonate anion forms long-living hydrogen bonds with water and shows an amphiphilic character, in which the water molecules are expelled in a region close to its C(3) symmetry axis. The calcium cation forms a strongly bound ion pair with the carbonate. The first hydration shell around the CaCO(3) pair is found to be very similar to the location of the water molecules surrounding CaCO(3) in ikaite, the hydrated mineral.

Hydrogen bonding mediated by key orbital interactions determines hydration enthalpy differences of phosphate water clusters

Electronic structure calculations have been carried out to provide a molecular interpretation for dihydrogen phosphate stability in water relative to that of metaphosphate. Specifically, hydration enthalpies of biologically important metaphosphate and dihydrogen phosphate with one to three waters have been computed with second-order Møller-Plesset perturbation and density functional theory (B3LYP) with up to the aug-cc-pvtz basis set and compared to experiment. The inclusion of basis set superposition error corrections and supplemental diffuse functions are necessary to predict hydration enthalpies within experimental uncertainty. Natural bond orbital analysis is used to rationalize underlying hydrogen bond configurations and key orbital interactions responsible for the experimentally reported difference in hydration enthalpies between metaphosphate and dihydrogen phosphate. In general, dihydrogen phosphate forms stronger hydrogen bonds compared to metaphosphate due to a greater charge transfer or enhanced orbital overlap between the phosphoryl oxygen lone pairs, n(O), and the antibonding O-H bond of water. Intramolecular distal lone pair repulsion with the donor n(O) orbital of dihydrogen phosphate distorts symmetric conformations, which improves n(O) and sigma*(O-H) overlap and ultimately the hydrogen bond strength. Unlike metaphosphate, water complexed to dihydrogen phosphate can serve as both a hydrogen bond donor and a hydrogen bond acceptor, which results in cooperative charge transfer and a reduction of the energy gap between n(O) and sigma*(O-H), leading to stronger hydrogen bonds. This study offers insight into how orbital interactions mediate hydrogen bond strengths with potential implications on the understanding of the kinetics and mechanism in enzymatic phosphoryl transfer reactions.

Probing the intrinsic features and environmental stabilization of multiply charged anions

Multiply charged anions (MCAs) represent exotic, highly energetic species in the gas-phase due to their propensity to undergo unimolecular decay via electron loss or ionic fragmentation. There is considerable fundamental interest in these systems since they display novel potential energy surfaces that are characterized by Coulomb barriers. Over recent years, considerable progress has been made in understanding the factors that affect the stability, decay pathways and reactivity of gas-phase MCAs, mainly as a result of the application of electrospray ionization as a generic technique for transferring solution-phase MCAs into the gas-phase for detailed characterization. We review contemporary work in this field, focusing on the factors that control the intrinsic stability of MCAs, both as isolated gas-phase ions, and on their complexation with solvent molecules and counter-ions. While studies of MCAs are primarily of fundamental interest, several classes of important biological ions are commonly observed as MCAs in the gas-phase (e.g. oligonucleotides, sugars). Recent results for biologically relevant ions are emphasised, since a fundamental understanding of the properties of gas-phase MCAs will be highly valuable for developing further analytical methods to study these important systems.

Gas-Phase Stability of Tetrahedral Multiply Charged Anions: A Conceptual and Computational DFT Study

Multiply charged anions (MCA's) are unstable relative to electron autoejection; however, the repulsive Coulomb barrier (RCB) provides electronic stability. In view of their interest in biological systems, the behavior of isolated AsO(4)(3-), PO(4)(3-), SO(4)(2-), and SeO(4)(2-) in the gas phase and in solution has been studied. To calculate the RCB values, the electrostatic and point charge model-two methods currently used in the literature-are applied, together with a recently introduced Conceptual Density Functional Theory (DFT) based approach. The relative stability of the above-mentioned MCA's is compared. The trends of the RCB are analyzed by including analogous compounds from the second and third row and by passing from dianionic to trianionic systems. Considering the effect of solvent, using the SCI-PCM solvent model, the evolution of the RCB when passing to higher dielectric constants is evaluated. The RCB is related to the properties of the system as polarizability/softness. Both a numerical and a conceptual correlation between the RCB and the global softness is found.

Dianionic tetraborates do exist as stable entities

To date, B(6)H(6)(2-) and some of its derivatives are the smallest members of the closo-borates that have been synthesized and analyzed in condensed phases. In contrast, no stable dianionic tetraborate has yet been observed, either in solution or solids or in the gas phase. In this work, the gas-phase stability of dianionic tetraborates B(4)X(4)(2-) (X = H, CN, NC, or BO) is investigated with ab initio methods. For this objective, the geometries of the dianions are optimized, the electronic stability is tested, and various fragmentation channels are studied. In agreement with previous examinations, tetrahedral isomers of all examined tetraborates have been found to represent geometrically stable isomers and to exhibit a triplet electronic ground state. However, these isomers are electronically unstable, i.e., their additional electrons are not bound. Furthermore, new D(2)(d)()-symmetric isomers of B(4)X(4)(2-) (X = H, CN, NC, or BO) have been identified that have a closed-shell singlet ground state and are lower in energy than their tetrahedral counterparts. Moreover, B(4)(CN)(4)(2-) and B(4)(BO)(4)(2-) represent stable gas-phase dianions and are predicted to be observable in suitable experiments. The electronic properties and geometries of these dianions are discussed in detail and explained in terms of the electrostatic repulsion of the excess electrons and the aromaticity of the dianions.

Cyclic carbon cluster dianions and their aromaticity

Cyclic carbon cluster dianions (CC(2))(2-)(n)(n = 3-6) are investigated by ab initio methods with regard to their geometric properties, electronic stability, and aromaticity. The unique wheel-like structures of these dianions consist of a n-membered carbon ring, where a C(2) unit is attached to each carbon atom. All investigated dianions represent stable gas-phase dianions. While the smallest member of this family (CC(2))(2-)(3) is clearly aromatic, the aromatic character decreases rapidly with increasing ring size. The geometries and the aromaticity of the cyclic clusters (CC(2))(2-)(n)(n = 3-6) can be nicely explained using resonance structure arguments.