Analysis of the Ordering Effects in Anthraquinone Thin Films and Its Potential Application for Sodium Ion Batteries

Probing the range of applicability of structure- and energy-adjusted QM/MM link bonds III: QM/MM MD simulations of solid-state systems at the example of layered carbon structures

The previously introduced workflow to achieve an energetically and structurally optimized description of frontier bonds in quantum mechanical/molecular mechanics (QM/MM)-type applications was extended into the regime of computational material sciences at the example of a layered carbon model systems. Optimized QM/MM link bond parameters at HSEsol/6-311G(d,p) and self-consistent density functional tight binding (SCC-DFTB) were derived for graphitic systems, enabling detailed investigation of specific structure motifs occurring in graphene-derived structures v i a quantum-chemical calculations. Exemplary molecular dynamics (MD) simulations in the isochoric-isothermic (NVT) ensemble were carried out to study the intercalation of lithium and the properties of the Stone-Thrower-Wales defect. The diffusivity of lithium as well as hydrogen and proton adsorption on a defective graphene surface served as additional example. The results of the QM/MM MD simulations provide detailed insight into the applicability of the employed link-bond strategy when studying intercalation and adsorption properties of graphitic materials.

Characterization of the Coordination and Solvation Dynamics of Solvated Systems─Implications for the Analysis of Molecular Interactions in Solutions and Pure H2O

The characterization of solvation shells of atoms, ions, and molecules in solution is essential to relate solvation properties to chemical phenomena such as complex formation and reactivity. Different definitions of the first-shell coordination sphere from simulation data can lead to potentially conflicting data on the structural properties and associated ligand exchange dynamics. The definition of a solvation shell is typically based on a given threshold distance determined from the respective solute-solvent pair distribution function g(r) (i.e., GC). Alternatively, a nearest neighbor (NN) assignment based on geometric properties of the coordination complex without the need for a predetermined cutoff criterion, such as the relative angular distance (RAD) or the modified Voronoi (MV) tessellation, can be applied. In this study, the effect of different NN algorithms on the coordination number and ligand exchange dynamics evaluated for a series of monatomic ions in aqueous solution, carbon dioxide in aqueous and dichloromethane solutions, and pure liquid water has been investigated. In the case of the monatomic ions, the RAD approach is superior in achieving a well separated definition of the first solvation layer. In contrast, the MV algorithm provides a better separation of the NNs from a molecular point of view, leading to better results in the case of solvated CO2. When analyzing the coordination environment in pure water, the cutoff-based GC framework was found to be the most reliable approach. By comparison of the number of ligand exchange reactions and the associated mean ligand residence times (MRTs) with the properties of the coordination number autocorrelation functions, it is shown that although the average coordination numbers are sensitive to the different definitions of the first solvation shell, highly consistent estimates for the associated MRT of the solvated system are obtained in the majority of cases.

Perylenetetracarboxylic Diimide Composite Electrodes as Organic Cathode Materials for Rechargeable Sodium-Ion Batteries: A Joint Experimental and Theoretical Study

The organic semiconductor 3,4,9,10-perylenetetracarboxylic diimide (PTCDI), a widely used industrial pigment, has been identified as a diffusion-less Na-ion storage material, allowing for exceptionally fast charging/discharging rates. The elimination of diffusion effects in electrochemical measurements enables the assessment of interaction energies from simple cyclic voltammetry experiments through the theoretical work of Laviron and Tokuda. In this work, the two N-substituted perylenes, N,N'-dimethyl-3,4,9,10-perylenetetracarboxylic diimide (Me2PTCDI) and N,N'-diphenyl-3,4,9,10-perylenetetracarboxylic diimide (Ph2PTCDI), as well as the parent molecule 3,4,9,10-perylenetetracarboxylic diimide (H2PTCDI) are investigated as thin-film composite electrodes on carbon fibers for sodium-ion batteries. The composite electrodes are analyzed with Raman spectroscopy. Interaction parameters are extracted from cyclic voltammetry measurements. The stability and rate capability of the three PTCDI derivatives are examined through galvanostatic measurements in sodium-ion half-cell batteries and the influence of the interactions on those parameters is evaluated. In addition, self-consistent charge density function tight binding calculations of the different PTCDI systems interacting with graphite have been carried out. The results show that the binding motif displays notable deviations from an ideal ABA stacking, especially for the neutral state. In addition, data obtained for the electron-transfer integrals show that the difference in performance between different PTCDI thin-film batteries cannot be solely explained by the electron-transfer properties and other factors such as H-bonding have to be considered.

Towards hybrid quantum mechanical/molecular mechanical simulations of Li and Na intercalation in graphite - force field development and DFTB parametrisation

In this work a previously established QM/MM simulation protocol for the treatment of solid-state interfaces was extended towards the treatment of layered bulk materials enabling for instance investigation of metal intercalation in graphitic carbon materials. In order to study the intercalation of Li in graphite, new density functional tight binding (DFTB) parameters for Li have been created. Molecular dynamics (MD) simulations at constant temperatures (273.15, 298.15 and 323.15 K) have been carried out to assess the performance of the presented DFTB MD simulation approach. The intercalation of variable lithium and sodium content was investigated via z-distribution functions and analysis of the diffusivity in the direction parallel to the graphene plane. Both the calculated diffusion coefficients and the activation energy in case of lithium are in good agreement with experimental data. The comparison of the QM/MM MD simulation results provide detailed insights into the structural and dynamical properties of intercalated metal ions.

How to use a rotating ring-disc electrode (RRDE) subtraction method to investigate the electrocatalytic oxygen reduction reaction?

When studying electrochemical oxygen reduction reactions in homogeneous media, special attention must be given to the significant background activity present with conventional electrode materials. The intrinsic electrocatalytic activity of different materials can be investigated using complementary methods, such as the rotating ring-disc electrode (RRDE) technique and chronoamperometric electrolysis with product quantification. This report presents a detailed investigation of the electrocatalytic ability of hydroxy anthraquinone derivatives and riboflavin towards hydrogen peroxide (H2O2) production via a novel RRDE subtraction method together with chronoamperometric electrolysis. Qualitative trends linking the two methods were obtained, such as a higher excess current correlating with both higher productivity and selectivity. As such, a valuable tool is provided to increase the understanding of the electrocatalytic ability of homogeneous solutions toward improving the oxygen reduction reaction.

Nanostructural Organization in a Biredox Ionic Liquid

Ionic liquids generally display peculiar structural features that impact their physical properties, such as the formation of polar and apolar domains. Recently, ionic liquids functionalized with anthraquinone and TEMPO redox groups were shown to increase the energy storage performance of supercapacitors, but their structure has not yet been characterized. In this work, we use polarizable molecular dynamics to study the nanostructuration of such biredox ionic liquids. We show that TEMPO nitroxyl functions tend to aggregate, while the anthraquinone groups favor stacked arrangements. The latter eventually percolate through the whole liquid, which sheds some light on the mechanisms at play within biredox ionic liquid-based supercapacitors.

Substantial Na-Ion Storage at High Current Rates: Redox-Pseudocapacitance through Sodium Oxide Formation

Batteries and supercapacitors, both governed by electrochemical processes, operate by different electrochemical mechanisms which determine their characteristic energy and power densities. Battery materials store large amounts of energy by ion intercalation. Electrical double-layer capacitors store charge through surface-controlled ion adsorption which leads to high power and rapid charging, but much smaller amounts of energy stored. Pseudocapacitive materials offer the promise to combine these properties by storing charge through surface-controlled, battery-like redox reactions but at high rates approaching those of electrochemical double-layer capacitors. This work compares the pseudo-capacitive charge storage characteristics of self-organized titanium dioxide (TiO_2-x) nanotubes (NTs) to flat TiO_2-x surface films to further elucidate the proposed charge storage mechanism within the formed surface films. By comparing TiO_2-x NTs to flat TiO_2-x surface films, having distinctively different oxide mass and surface area ratios, it is shown that NaO₂ and Na₂O₂ formation, which constitutes the active surface film material, is governed by the metal oxide bulk. Our results corroborate that oxygen diffusion from the lattice oxide is key to NaO₂ and Na₂O₂ formation.

Direct Electrochemical CO2 Capture Using Substituted Anthraquinones in Homogeneous Solutions: A Joint Experimental and Theoretical Study

Electrochemical capture of carbon dioxide (CO2) using organic quinones is a promising and intensively studied alternative to the industrially established scrubbing processes. While recent studies focused only on the influence of substituents having a simple mesomeric or nucleophilicity effect, we have systematically selected six anthraquinone (AQ) derivatives (X-AQ) with amino and hydroxy substituents in order to thoroughly study the influence thereof on the properties of electrochemical CO2 capture. Experimental data from cyclic voltammetry (CV) and UV-Vis spectroelectrochemistry of solutions in acetonitrile were analyzed and compared with innovative density functional tight binding computational results. Our experimental and theoretical results provide a coherent explanation of the influence of CO2 on the CV data in terms of weak and strong binding nomenclature of the dianions. In addition to this terminology, we have identified the dihydroxy substituted AQ as a new class of molecules forming rather unstable [X-AQ-(CO2) n ]2- adducts. In contrast to the commonly used dianion consideration, the results presented herein reveal opposite trends in stability for the X-AQ-CO2 •- radical species for the first time. To the best of our knowledge, this study presents theoretically calculated UV-Vis spectra for the various CO2-AQ reduction products for the first time, enabling a detailed decomposition of the spectroelectrochemical data. Thus, this work provides an extension of the existing classification with proof of the existence of X-AQ-CO2 species, which will be the basis of future studies focusing on improved materials for electrochemical CO2 capture.

Anthraquinone and its derivatives as sustainable materials for electrochemical applications - a joint experimental and theoretical investigation of the redox potential in solution

Anthraquinone (AQ) has long been identified as a highly promising lead structure for various applications in organic electronics. Considering the enormous number of possible substitution patterns of the AQ lead structure, with only a minority being commercially available, a systematic experimental screening of the associated electrochemical potentials represents a highly challenging and time consuming task, which can be greatly enhanced via suitable virtual pre-screening techniques. In this work the calculated electrochemical reduction potentials of pristine AQ and 12 hydroxy- or/and amino-substituted AQ derivatives in N,N-dimethylformamide have been correlated against newly measured experimental data. In addition to the calculations performed using density functional theory (DFT), the performance of different semi-empirical density functional tight binding (DFTB) approaches has been critically assessed. It was shown that the SCC DFTB/3ob parametrization in conjunction with the COSMO solvation model provides a highly adequate description of the electrochemical potentials also in the case of the two-fold reduced species. While the quality in the correlation against the experimental data proved to be slightly inferior compared to the employed DFT approach, the highly advantageous cost-accuracy ratio of the SCC DFTB/3ob/COSMO framework has important implications in the formulation of hierarchical screening strategies for materials associated with organic electronics. Based on the observed performance, the low-cost method provides sufficiently accurate results to execute efficient pre-screening protocols, which may then be followed by a DFT-based refinement of the best candidate structures to facilitate a systematic search for new, high-performance organic electronic materials.