“Ion Solvation Spectra”: Free Energy Analysis of Solvation Structures of Multivalent Cations in Aprotic Solvents

Using advanced molecular dynamics free energy sampling techniquesboth classical and ab initiowe analyze the solvation structures of multivalent cations in aprotic solvents. In contrast to previous studies of mono- and bivalent ions in organic solvents, mainly performed using hybrid cluster-continuum quantum chemistry calculations that rely on the assumption of uniqueness of ion solvation free energies, here we find that monatomic bivalent cations may have multiple well-defined minima, as previously reported only for water, or plateaus of free energy with respect to the ion–solvent coordination. These observations are generalized in the concept of the “ion solvation spectrum“ to highlight the rich phenomenology related to ion solvation as opposed to the normally expected free energy profiles with a single coordination minimum. Specifically, we show that a single chemical species may exhibit a multiplicity of distinctly different electrochemical properties. Using one- and two-dimensional projections of the free energy landscape, we analyze the stability of ion solvation structures and reveal minimum free energy pathways for ion (de-)­solvation with low-dimensional approximations to associated kinetic barriers. Unexpectedly, we show that in some cases the process of opening the first ion solvation shell, by removing a solvent molecule, may actually drive the ion into a free energy basin with a higher coordination number. Our study highlights some deficiencies of conventional methodologies for studying ion solvation as a path to determine redox potentials and provides experimentally testable predictions

The Solvation Structure of Mg Ions in Dichloro Complex Solutions from First-Principles Molecular Dynamics and Simulated X‑ray Absorption Spectra

The knowledge of Mg solvation structure in the electrolyte is requisite to understand the transport behavior of Mg ions and their dissolution/deposition mechanism at electrolyte/electrode interfaces. In the first established rechargeable Mg-ion battery system [D. Aurbach et al. <i>Nature</i> <b>2000</b>, <i>407</i>, 724], the electrolyte is of the dichloro complex (DCC) solution family, Mg­(AlCl₂BuEt)₂/THF, resulting from the reaction of Bu₂Mg and EtAlCl₂ with a molar ratio of 1:2. There is disagreement in the literature regarding the exact solvation structure of Mg ions in such solutions, i.e., whether Mg²⁺ is tetra- or hexacoordinated by a combination of Cl^– and THF. In this work, theoretical insight into the solvation complexes present is provided based on first-principles molecular dynamics simulations (FPMD). Both Mg monomer and dimer structures are considered in both neutral and positively charged states. We found that, at room temperature, the Mg²⁺ ion tends to be tetracoordinated in the THF solution phase instead of hexacoordinated, which is the predominant solid-phase coordination. Simulating the X-ray absorption spectra (XAS) at the Mg K-edge by sampling our FPMD trajectories, our predicted solvation structure can be readily compared with experimental measurements. It is found that when changing from tetra- to hexacoordination, the onset of X-ray absorption should exhibit at least a 1 eV blue shift. We propose that this energy shift can be used to monitor changes in the Mg solvation sphere as it migrates through the electrolyte to electrolyte/electrode interfaces and to elucidate the mechanism of Mg dissolution/deposition

Exploration of the Detailed Conditions for Reductive Stability of Mg(TFSI)₂ in Diglyme: Implications for Multivalent Electrolytes

We reveal the general mechanisms of partial reduction of multivalent complex cations in conditions specific for the bulk solvent and in the vicinity of the electrified metal electrode surface and disclose the factors affecting the reductive stability of electrolytes for multivalent electrochemistry. Using a combination of <i>ab initio</i> techniques, we clarify the relation between the reductive stability of contact-ion pairs comprising a multivalent cation and a complex anion, their solvation structures, solvent dynamics, and the electrode overpotential. We found that for ion pairs with multiple configurations of the complex anion and the Mg cation whose available orbitals are partially delocalized over the molecular complex and have antibonding character, the primary factor of the reductive stability is the shape factor of the solvation sphere of the metal cation center and the degree of the convexity of a polyhedron formed by the metal cation and its coordinating atoms. We focused specifically on the details of Mg (II) bis­(trifluoro­methane­sulfonyl)­imide in diethylene glycol dimethyl ether (Mg­(TFSI)₂)/diglyme) and its singly charged ion pair, MgTFSI⁺. In particular, we found that both stable (MgTFSI)⁺ and (MgTFSI)⁰ ion pairs have the same TFSI configuration but drastically different solvation structures in the bulk solution. This implies that the MgTFSI/dyglyme reductive stability is ultimately determined by the relative time scale of the solvent dynamics and electron transfer at the Mg–anode interface. In the vicinity of the anode surface, steric factors and hindered solvent dynamics may increase the reductive stability of (MgTFSI)⁺ ion pairs at lower overpotential by reducing the metal cation coordination, in stark contrast to the reduction at high overpotential accompanied by TFSI decomposition. By examining other solute/solvent combinations, we conclude that the electrolytes with highly coordinated Mg cation centers are more prone to reductive instability due to the chemical decomposition of the anion or solvent molecules. The obtained findings disclose critical factors for stable electrolyte design and show the role of interfacial phenomena in reduction of multivalent ions

Tuning Semiconductor Band Edge Energies for Solar Photocatalysis via Surface Ligand Passivation

Semiconductor photocatalysts capable of broadband solar photon absorption may be nonetheless precluded from use in driving water splitting and other solar-to-fuel related reactions due to unfavorable band edge energy alignment. Using first-principles density functional theory and beyond, we calculate the electronic structure of passivated CdSe surfaces and explore the opportunity to tune band edge energies of this and related semiconductors via electrostatic dipoles associated with chemisorbed ligands. We predict substantial shifts in band edge energies originating from both the induced dipole at the ligand/CdSe interface and the intrinsic dipole of the ligand. Building on important induced dipole contributions, we further show that, by changing the size and orientation of the ligand’s intrinsic dipole moment via functionalization, we can control the direction and magnitude of the shifts of CdSe electronic levels. Our calculations suggest a general strategy for enabling new active semiconductor photocatalysts with both optimal opto-electronic, and photo- and electrochemical properties

Critical Factors in Computational Characterization of Hydrogen Storage in Metal–Organic Frameworks

Inconsistencies in high-pressure H2 adsorption data and a lack of comparative experiment–theory studies have made the evaluation of both new and existing metal–organic frameworks (MOFs) challenging in the context of hydrogen storage applications. In this work, we performed grand canonical Monte Carlo (GCMC) simulations in nearly 500 experimentally refined MOF structures to examine the variance in simulation results because of the equation of state, H2 potential, and the effect of density functional theory structural optimization. We find that hydrogen capacity at 77 K and 100 bar, as well as hydrogen 100-to-5 bar deliverable capacity, is correlated more strongly with the MOF pore volume than with the MOF surface area (the latter correlation is known as the Chahine’s rule). The tested methodologies provide consistent rankings of materials. In addition, four prototypical MOFs (MOF-74, CuBTC, ZIF-8, and MOF-5) with a range of surface areas, pore structures, and surface chemistries, representative of promising adsorbents for hydrogen storage, are evaluated in detail with both GCMC simulations and experimental measurements. Simulations with a three-site classical potential for H2 agree best with our experimental data except in the case of MOF-5, in which H2 adsorption is best replicated with a five-site potential. However, for the purpose of ranking materials, these two choices for H2 potential make little difference. More significantly, 100 bar loading estimates based on more accurate equations of state for the vapor–liquid equilibrium yield the best comparisons with the experiment

DataSheet1_Solar enhanced oxygen evolution reaction with transition metal telluride.PDF

The photo-enhanced electrocatalytic method of oxygen evolution reaction (OER) shows promise for enhancing the effectiveness of clear energy generation through water splitting by using renewable and sustainable source of energy. However, despite benefits of photoelectrocatalytic (PEC) water splitting, its uses are constrained by its low efficiency as a result of charge carrier recombination, a large overpotential, and sluggish reaction kinetics. Here, we illustrate that Nickel telluride (NiTe) synthesized by hydrothermal methods can function as an extremely effective photo-coupled electrochemical oxygen evolution reaction (POER) catalyst. In this study, NiTe was synthesized by hydrothermal method at 145°C within just an hour of reaction time. In dark conditions, the NiTe deposited on carbon cloth substrate shows a small oxygen evolution reaction overpotential (261 mV) at a current density of 10 mA cm–2, a reduced Tafel slope (65.4 mV dec−1), and negligible activity decay after 12 h of chronoamperometry. By virtue of its enhanced photo response, excellent light harvesting ability, and increased interfacial kinetics of charge separation, the NiTe electrode under simulated solar illumination displays exceptional photoelectrochemical performance exhibiting overpotential of 165 mV at current density of 10 mA cm-2, which is about 96 mV less than on dark conditions. In addition, Density Functional Theory investigations have been carried out on the NiTe surface, the results of which demonstrated a greater adsorption energy for intermediate -OH on the catalyst site. Since the -OH adsorption on the catalyst site correlates to catalyst activation, it indicates the facile electrocatalytic activity of NiTe owing to favorable catalyst activation. DFT calculations also revealed the facile charge density redistribution following intermediate -OH adsorption on the NiTe surface. This work demonstrates that arrays of NiTe elongated nanostructure are a promising option for both electrochemical and photoelectrocatalytic water oxidation and offers broad suggestions for developing effective PEC devices.</p

Understanding the Impact of Multi-Chain Ion Coordination in Poly(ether-Acetal) Electrolytes

Performant solid polymer electrolytes for battery applications usually have a low glass transition temperature and good ion solvation. Recently, to understand the success of PEO for solid-sate battery applications and explore alternatives, we have studied a series of polyacetals along with PEO, both from an experimental and a computational standpoint. We observed that even though the mechanism of transport may be more optimal in polyacetals, the lower glass transition temperature of the PEO-salt electrolyte system still makes it the best option, in this class of polymers, for battery applications. In this work, we explored the free-energy landscape of PEO and P(EO-MO) at various compositions and temperatures using metadynamics simulations to gain deeper insights into the various factors that affect the glass transition temperatures in these systems. In particular, we study the competition between intra- and inter-chain coordination of the cation in these systems that we had hypothesized in our previous work was responsible for the differences in the glass transition temperature. We observe that in PEO, the single-chain binding motif is thermodynamically more stable than the multi-chain binding motif, unlike P(EO-MO), where the opposite is true. We also show that multi-chain coordination, and the associated higher glass transition temperature, in P(EO-MO) is due to a larger strain energy for single-chain coordination that originates in the introduced OCO linkages (relative to PEO’s consistent OCCO linkages). Furthermore, the type of pathways to move from one transition state to another in the various systems do not change at higher concentrations though the relative probability of cation–anion coordinated states increases. Calculations at different temperatures to understand the entropic effect on the stability of these coordination environments reveal that as we increase the temperature, single-chain coordination becomes relatively more stable due to the entropic cost of multi-chain coordination, reducing the number of accessible states for the polymer. The various insights into the factors that affect glass transition temperature in these systems suggest design principles for polymer electrolyte systems with lower glass transition temperatures that need further research to compete with PEO at the same absolute battery working temperatures

Critical Factors in Computational Characterization of Hydrogen Storage in Metal–Organic Frameworks

Inconsistencies in high-pressure H2 adsorption data and a lack of comparative experiment–theory studies have made the evaluation of both new and existing metal–organic frameworks (MOFs) challenging in the context of hydrogen storage applications. In this work, we performed grand canonical Monte Carlo (GCMC) simulations in nearly 500 experimentally refined MOF structures to examine the variance in simulation results because of the equation of state, H2 potential, and the effect of density functional theory structural optimization. We find that hydrogen capacity at 77 K and 100 bar, as well as hydrogen 100-to-5 bar deliverable capacity, is correlated more strongly with the MOF pore volume than with the MOF surface area (the latter correlation is known as the Chahine’s rule). The tested methodologies provide consistent rankings of materials. In addition, four prototypical MOFs (MOF-74, CuBTC, ZIF-8, and MOF-5) with a range of surface areas, pore structures, and surface chemistries, representative of promising adsorbents for hydrogen storage, are evaluated in detail with both GCMC simulations and experimental measurements. Simulations with a three-site classical potential for H2 agree best with our experimental data except in the case of MOF-5, in which H2 adsorption is best replicated with a five-site potential. However, for the purpose of ranking materials, these two choices for H2 potential make little difference. More significantly, 100 bar loading estimates based on more accurate equations of state for the vapor–liquid equilibrium yield the best comparisons with the experiment

Critical Factors in Computational Characterization of Hydrogen Storage in Metal–Organic Frameworks

Inconsistencies in high-pressure H2 adsorption data and a lack of comparative experiment–theory studies have made the evaluation of both new and existing metal–organic frameworks (MOFs) challenging in the context of hydrogen storage applications. In this work, we performed grand canonical Monte Carlo (GCMC) simulations in nearly 500 experimentally refined MOF structures to examine the variance in simulation results because of the equation of state, H2 potential, and the effect of density functional theory structural optimization. We find that hydrogen capacity at 77 K and 100 bar, as well as hydrogen 100-to-5 bar deliverable capacity, is correlated more strongly with the MOF pore volume than with the MOF surface area (the latter correlation is known as the Chahine’s rule). The tested methodologies provide consistent rankings of materials. In addition, four prototypical MOFs (MOF-74, CuBTC, ZIF-8, and MOF-5) with a range of surface areas, pore structures, and surface chemistries, representative of promising adsorbents for hydrogen storage, are evaluated in detail with both GCMC simulations and experimental measurements. Simulations with a three-site classical potential for H2 agree best with our experimental data except in the case of MOF-5, in which H2 adsorption is best replicated with a five-site potential. However, for the purpose of ranking materials, these two choices for H2 potential make little difference. More significantly, 100 bar loading estimates based on more accurate equations of state for the vapor–liquid equilibrium yield the best comparisons with the experiment

Mechanistic Advantages of Organotin Molecular EUV Photoresists

Extreme ultraviolet (EUV)-induced radiation exposure chemistry in organotin–oxo systems, represented by the archetypal [(R–Sn)12O14(OH)6]­(A)2 cage, has been investigated with density functional theory. Upholding existing experimental evidence of Sn–C cleavage-dominant chemistry, computations have revealed that either electron attachment or ionization can single-handedly trigger tin–carbon bond cleavage, partially explaining the current EUV sensitivity advantage of metal oxide systems. We have revealed that tin atoms at different parts of the molecule react differently to ionization and electron attachment and have identified such selectivity as a result of local coordination chemistry instead of the macro geometry of the molecule. An ionization–deprotonation pathway has also been identified to explain the observed evolution of an anion conjugate acid upon exposure and anion mass dependence in resist sensitivity