The sensitive aspects of modelling polymer-ceramic composite solid-state electrolytes using molecular dynamics simulations.

Solid-state composite electrolytes have arisen as one of the most promising materials classes for next-generation Li-ion battery technology. These composites mix ceramic and solid-polymer ion conductors with the aim of combining the advantages of each material. The ion-transport mechanisms within such materials, however, remain elusive. This knowledge gap can to a large part be attributed to difficulties in studying processes at the ceramic-polymer interface, which are expected to play a major role in the overall ion transport through the electrolyte. Computational efforts have the potential of providing significant insight into these processes. One of the main challenges to overcome is then to understand how a sufficiently robust model can be constructed in order to provide reliable results. To this end, a series of molecular dynamics simulations are here carried out with a variation of certain structural (surface termination and polymer length) and pair potential (van der Waals parameters and partial charges) models of the Li7La3Zr2O12 (LLZO) poly(ethylene oxide) (PEO) system, in order to test how sensitive the outcome is to each variation. The study shows that the static and dynamic properties of Li-ion are significantly affected by van der Waals parameters as well as the surface terminations, while the thickness of the interfacial region - where the structure-dynamic properties are different as compared to the bulk-like regime - is the same irrespective of the simulation setup.

Theoretical investigation of solvent and oxidation/deprotonation effects on the electronic structure of a mononuclear Ru-aqua-polypyridine complex in aqueous solution.

Mononuclear polypyridine ruthenium (Ru) complexes can catalyze various reactions, including water splitting, and can also serve as photosensitizers in solar cells. Despite recent progress in their synthesis, accurately modeling their physicochemical properties, particularly in solution, remains challenging. Herein, we conduct a theoretical investigation of the structural and electronic properties of a mononuclear Ru-aqua polypyridine complex in aqueous solution, considering five of its possible oxidation/protonation states species: [RuII(H2O)(py)(bpy)2]2+, [RuII(OH)(py)(bpy)2]+, [RuIII(H2O)(py)(bpy)2]3+, [RuIII(OH)(py)(bpy)2]2+ and [RuIV(O)(py)(bpy)2]2+, where py = pyridine and bpy = 2,2'-bipyridine. At first, we investigate the impact of proton-coupled and non-coupled electron transfer reactions on the geometry and electronic structure of the complexes in vacuum and in solution, using an implicit solvent model. Then, using a sequential multiscale approach that combines quantum mechanics and molecular mechanics (S-QM/MM), we examine the explicit solvent effects on the electronic excitations of the complexes, and compare them with the experimental results. The complexes were synthesized, and their absorption spectra measured in aqueous solution. To accurately describe the QM interactions between the metal center and the aqueous ligand in the MM simulations, we developed new force field parameters for the Ru atom. We analyze the solvent structure around the complexes and account for its explicit influence on the polarization and electronic excitations of the complexes. Notably, accounting for the explicit solvent polarization effects of the first solvation shells is essential to correctly describe the energy of the electronic transitions, and the explicit treatment of the hydrogen bonds at the QM level in the excitation calculations improves the accuracy of the description of the metal-to-ligand charge-transfer bands. Transition density matrix analysis is used to characterize all electronic transitions in the visible and ultraviolet ranges according to their charge-transfer (CT) character. This study elucidates the electronic structure of those ruthenium polypyridyl complexes in aqueous solution and underscores the importance of precisely describing solvent effects, which can be achieved employing the S-QM/MM method.

Unveiling the impact of exchange-correlation functionals on the description of key electronic properties of non-fullerene acceptors in organic photovoltaics.

Non-fullerene electron acceptors have emerged as promising alternatives to traditional electron-acceptors in the active layers of organic photovoltaics. This is due to their tunable energy levels, optical response in the visible light spectrum, high electron mobility, and photochemical stability. In this study, the electronic properties of two representative non-fullerene acceptors, ITIC and Y5, have been calculated within the framework of density functional theory using a range of hybrid and non-hybrid density functionals. Screened range-separated hybrid (SRSH) approaches were also tested. The results are analyzed in light of the previously reported experimental outcomes. Specifically, we have calculated the oxidation and reduction potentials, fundamental and optical gaps, the highest occupied molecular orbital and lowest unoccupied molecular orbital energies, and exciton binding energies. Additionally, we have investigated the effects of the medium dielectric constant on these properties employing a universal implicit solvent model. It was found that hybrid functionals generally perform poorly in predicting oxidation potentials, while non-hybrid functionals tend to overestimate reduction potentials. The inclusion of a large Hartree-Fock contribution to the global or long range was identified as the source of inaccuracy for many hybrid functionals in predicting both redox potentials and the fundamental and optical gaps. Corroborating with the available literature, ∼50% of all tested functionals predicted very small exciton binding energies, within the range of ±0.1 eV, that become even smaller by increasing the dielectric constant of the material. Finally, the OHSE2PBE and tHCTHhyb functionals and the optimal tuning SRSH approach emerged as the best-performing methods, with good accuracy in the description of the electronic properties of interest.

On the Conformation of Dimeric Acceptors and Their Polymer Solar Cells with Efficiency over 18 .

The determination of molecular conformations of oligomeric acceptors (OAs) and their impact on molecular packing are crucial for understanding the photovoltaic performance of their resulting polymer solar cells (PSCs) but have not been well studied yet. Herein, we synthesized two dimeric acceptor materials, DIBP3F-Se and DIBP3F-S, which bridged two segments of Y6-derivatives by selenophene and thiophene, respectively. Theoretical simulation and experimental 1D and 2D NMR spectroscopic studies prove that both dimers exhibit O-shaped conformations other than S- or U-shaped counter-ones. Notably, this O-shaped conformation is likely governed by a distinctive "conformational lock" mechanism, arising from the intensified intramolecular π-π interactions among their two terminal groups within the dimers. PSCs based on DIBP3F-Se deliver a maximum efficiency of 18.09 %, outperforming DIBP3F-S-based cells (16.11 %) and ranking among the highest efficiencies for OA-based PSCs. This work demonstrates a facile method to obtain OA conformations and highlights the potential of dimeric acceptors for high-performance PSCs.

Unraveling the acid-base characterization and solvent effects on the structural and electronic properties of a bis-bidentate bridging ligand.

Understanding the interactions and the solvent effects on the distribution of several species in equilibrium and how it can influence the 1H-NMR properties, spectroscopy (UV-vis absorption), and the acid-base equilibria can be especially challenging. This is the case of a bis-bidentate bridging ligand bis(2-pyridyl)-benzo-bis(imidazole), where the two pyridyl and four imidazolyl nitrogen atoms can be protonated in different ways, depending on the solvent, generating many isomeric/tautomeric species. Herein, we report a combined theoretical-experimental approach based on a sequential quantum mechanics/molecular mechanics procedure that was successfully applied to describe in detail the acid-base characterization and its effects on the electronic properties of such a molecule in solution. The calculated free-energies allowed the identification of the main species present in solution as a function of the solvent polarity, and its effects on the magnetic shielding of protons (1H-NMR chemical shifts), the UV-vis absorption spectra, and the acid-base equilibrium constants (pKas) in aqueous solution. Three acid-base equilibrium constants were experimentally/theoretically determined (pKa1 = 1.3/1.2, pKa2 = 2.1/2.2 and pKa5 = 10.1/11.3) involving mono-deprotonated and mono-protonated cis and trans species. Interestingly, other processes with pKa3 = 3.7 and pKa4 = 6.0 were also experimentally determined and assigned to the protonation and deprotonation of dimeric species. The dimerization of the most stable neutral species was investigated by Monte Carlo simulations and its electronic effects were considered for the elucidation of the UV-vis absorption bands, revealing transitions mainly with the charge-transfer characteristic and involving both the monomeric species and the dimeric species. The good matching of the theoretical and experimental results provides an atomistic insight into the solvent effects on the electronic properties of this bis-bidentate bridging ligand.

Understanding the lithiation limits of high-capacity organic battery anodes by atomic charge derivative analysis.

The superlithiation of organic anodes is a promising approach for developing the next generation of sustainable Li-ion batteries with high capacity. However, the lack of fundamental understanding hinders its faster development. Here, a systematic study of the lithiation processes in a set of dicarboxylate-based materials is carried out within the density functional theory formalism. It is demonstrated that a combined analysis of the Li insertion reaction thermodynamics and the conjugated-moiety charge derivative enables establishing the experimentally observed maximum storage, thus allowing an assessment of the structure-function relationships also.

Small Organic Molecule Based on Benzothiadiazole for Electrocatalytic Hydrogen Production.

A small organic molecule 2,1,3-benzothiadiazole-4, 7-dicarbonitrile (BTDN) is assessed for electrocatalytic hydrogen evolution on glassy carbon electrode and shows a hydrogen production Faradaic efficiency of 82% in the presence of salicylic acid. The key catalytic intermediates of reduced species BTDN-• and protonated intermediates are characterized or hypothesized by using various spectroscopic methods and density functional theory (DFT)-based calculations. With the experimental and theoretical results, a catalytic mechanism of BTDN for electrocatalytic H2 evolution is proposed.

Understanding carbon dioxide capture on metal-organic frameworks from first-principles theory: The case of MIL-53(X), with X = Fe3+, Al3+, and Cu2.

Metal-organic frameworks (MOFs) constitute a class of three-dimensional porous materials that have shown applicability for carbon dioxide capture at low pressures, which is particularly advantageous in dealing with the well-known environmental problem related to the carbon dioxide emissions into the atmosphere. In this work, the effect of changing the metallic center in the inorganic counterpart of MIL-53 (X), where X = Fe3+, Al3+, and Cu2+, has been assessed over the ability of the porous material to adsorb carbon dioxide by means of first-principles theory. In general, the non-spin polarized computational method has led to adsorption energies in fair agreement with the experimental outcomes, where the carbon dioxide stabilizes at the pore center through long-range interactions via oxygen atoms with the axial hydroxyl groups in the inorganic counterpart. However, spin-polarization effects in connection with the Hubbard corrections, on Fe 3d and Cu 3d states, were needed to properly describe the metal orbital occupancy in the open-shell systems (Fe- and Cu-based MOFs). This methodology gave rise to a coherent high-spin configuration, with five unpaired electrons, for Fe atoms leading to a better agreement with the experimental results. Within the GGA+U level of theory, the binding energy for the Cu-based MOF is found to be Eb = -35.85 kJ/mol, which is within the desirable values for gas capture applications. Moreover, it has been verified that the adsorption energetics is dominated by the gas-framework and internal weak interactions.

Spectroscopic Fingerprints of Intermolecular H-Bonding Interactions in Carbon Nitride Model Compounds.

The effect of intermolecular H-bonding interactions on the local electronic structure of N-containing functional groups (amino group and pyridine-like N) that are characteristic of polymeric carbon nitride materials p-CN(H), a new class of metal-free organophotocatalysts, was investigated. Specifically, the melamine molecule, a building block of p-CN(H), was characterized by X-ray photoelectron (XPS) and near edge X-ray absorption fine structure (NEXAFS) spectroscopy. The molecule was studied as a noninteracting system in the gas phase and in the solid state within a H-bonded network. With the support of DFT simulations of the spectra, it was found that the H-bonds mainly affect the N 1s level of the amino group, leaving the N 1s level of the pyridine-like N mostly unperturbed. This is responsible for a reduction of the chemical shift between the two XPS N 1s levels relative to free melamine. Consequently, N K-edge NEXAFS resonances involving the amino N 1s level also shift to lower photon energies. Moreover, the solid-state absorption spectra showed significant modification/quenching of resonances related to transitions from the amino N 1s level to σ* orbitals involving the NH2 termini.

Assessing the electrochemical properties of polypyridine and polythiophene for prospective applications in sustainable organic batteries.

Conducting polymers are being considered promising candidates for sustainable organic batteries mainly due to their fast electron transport properties and high recyclability. In this work, the key properties of polythiophene and polypyridine have been assessed through a combined theoretical and experimental study focusing on such applications. A theoretical protocol has been developed to calculate redox potentials in solution within the framework of the density functional theory and using continuous solvation models. Here, the evolution of the electrochemical properties of solvated oligomers as a function of the length of the chain is analyzed and then the polymer properties are estimated via linear regressions using ordinary least square. The predicted values were verified against our electrochemical experiments. This protocol can now be employed to screen a large database of compounds in order to identify organic electrodes with superior properties.

Electrolyte decomposition on Li-metal surfaces from first-principles theory.

An important feature in Li batteries is the formation of a solid electrolyte interphase (SEI) on the surface of the anode. This film can have a profound effect on the stability and the performance of the device. In this work, we have employed density functional theory combined with implicit solvation models to study the inner layer of SEI formation from the reduction of common organic carbonate electrolyte solvents (ethylene carbonate, propylene carbonate, dimethyl carbonate, and diethyl carbonate) on a Li metal anode surface. Their stability and electronic structure on the Li surface have been investigated. It is found that the CO producing route is energetically more favorable for ethylene and propylene carbonate decomposition. For the two linear solvents, dimethyl and diethyl carbonates, no significant differences are observed between the two considered reduction pathways. Bader charge analyses indicate that 2 e- reductions take place in the decomposition of all studied solvents. The density of states calculations demonstrate correlations between the degrees of hybridization between the oxygen of adsorbed solvents and the upper Li atoms on the surface with the trend of the solvent adsorption energies.

Crafting ferromagnetism in Mn-doped MgO surfaces with p-type defects.

We have employed first-principles calculations based on density functional theory (DFT) to investigate the underlying physics of unusual magnetism in Mn-doped MgO surface. We have studied two distinct scenarios. In the first one, two Mn atoms are substitutionally added to the surface, occupying the Mg sites. Both are stabilized in the Mn[Formula: see text] valence state carrying a local moment of 4.3 [Formula: see text] having a high-spin configuration. The magnetic interaction between the local moments display a very short-ranged characteristic, decaying very quickly with distance, and having antiferromagnetic ordering lower in energy. The energetics analysis also indicates that the Mn ions prefer to stay close to each other with an oxygen atom bridging the local interaction. In the second scenario, we started exploring the effect of native defects on the magnetism by crafting both Mg and O vacancies, which are p- and n-type defects, respectively. It is found that the electrons and holes affect the magnetic interaction between Mn ions in a totally different manner. The n-type defect leads to very similar magnetism, with the AFM configuration being energetically preferred. However, in the presence of Mg vacancy, the situation is quite different. The Mn atoms are further oxidized, giving rise to mixed Mn(d) ionic states. As a consequence, the Mn atoms couple ferromagnetically, when placed in the close configuration, and the obtained electronic structure is coherent with the double-exchange type of magnetic interaction. To guarantee the robustness of our results, we have benchmarked our calculations with three distinct theory levels, namely DFT-GGA, DFT-GGA+U and DFT-hybrid functionals. On the surface, the Mg vacancy displays lower formation energy occurring at higher concentrations. Therefore, our model systems can be the basis to explain a number of controversial results regarding transition metal doped oxides.

Inducing molecular orientation in solution-processed thin films of fluorene-bithiophene-based copolymer: thermal annealing vs. solvent additive.

A deep understanding of the factors influencing the morphology of thin films based on conjugated polymers is essential to boost their performance in optoelectronic devices. Herein, we investigated the electronic structure and morphology of thin films of the copolymer poly(9,9-dioctyl-fluorenyl-co-bithiophene) (F8T2) in its pristine form as well as samples processed with the solvent additive 1,8-diiodooctane (DIO) or post-processed through thermal annealing treatment. Measurements were carried out using angle-resolved S K-edge NEXAFS (near-edge X-ray absorption fine structure) in total electron yield (TEY) and fluorescence yield (FY) detection modes. Two main transitions were observed at the S 1s NEXAFS spectra: S 1s → π* and S 1s → σ* (S-C). The observed dichroism pointed to a face-on orientation of the conjugated backbone, which was significantly increased for F8T2 films processed with DIO. Resonant Auger decay spectra were obtained and analyzed using the core-hole clock (CHC) method. An enhancement in the charge transfer process was observed for thermally annealed films, especially for samples processed with DIO, corresponding to an increase in film ordering. Furthermore, the investigated films were characterized using X-ray photoelectron spectroscopy, attesting to the presence of the thiophene unit in the samples and demonstrating that some of its sulfur atoms were positively polarized in the F8T2 films. All these experimental findings were compared with molecular dynamics (MD) simulations of film evaporation with and without DIO. The use of MD, together with mathematical modeling, was able to explain the major effects found in the experiments, including the polarization of sulfur atoms. The simultaneous use of powerful spectroscopic techniques and theoretical methods shed light on key aspects linking film morphology with fabrication procedures.

C60-mediated hydrogen desorption in Li-N-H systems.

Hydrogen desorption from a LiH + NH(3) mixture is very difficult due to the formation of the stable LiNH(4) compound. Using cluster models and first-principles theory, we demonstrate that the C(60) molecule can in fact significantly improve the thermodynamics of ammonia-mediated hydrogen desorption from LiH due to the stabilization of the intermediate state, LiNH(4). The hydrogen desorption following the path of LiNH(4)-C(60) → LiNH(3)-C(60) + 1/2H(2) is exothermic. Molecular dynamic simulations show that this reaction can take place even at room temperature (300 K). In contrast, the stable LiNH(4) compound cannot desorb hydrogen at room temperature in the absence of C(60). The introduction of C(60) also helps to restrain the NH(3) gas which is poisonous in proton exchange membrane fuel cell applications.

Exploring Metastable Phases During Lithiation of Organic Battery Electrode Materials.

In this work, the Li-ion insertion mechanism in organic electrode materials is investigated through the lens of atomic-scale models based on first-principles theory. Starting with a structural analysis, the interplay of density functional theory with evolutionary and potential-mapping algorithms is used to resolve the crystal structure of the different (de)lithiated phases. These methods elucidate different lithiation reaction pathways and help to explore the formation of metastable phases and predict one- or multi-electron reactions, which are still poorly understood for organic intercalation electrodes. The cathode material dilithium 2,5-oxyterephthalate (operating at 2.6 V vs. Li/Li+) is investigated in depth as a case study, owing to its rich redox chemistry. When compared with recent experimental results, it is demonstrated that metastable phases with peculiar ring-ring molecular interactions are more likely to be controlling the redox reactions thermodynamics and consequently the battery voltage.

Theoretical study of C60 as catalyst for dehydrogenation in LiBH4.

Complex light metal hydrides possess many properties which make them attractive as a storage medium for hydrogen, but typically catalysts are required to lower the hydrogen desorption temperature and to facilitate hydrogen uptake in the form of a reversible reaction. The overwhelming focus in the search for catalysing agents has been on compounds containing titanium, but the precise mechanism of their actions remains somewhat obscure. A recent experiment has now shown that fullerenes (C(60)) can also act as catalysts for both hydrogen uptake and release in lithium borohydride (LiBH(4)). In an effort to understand the involved mechanism, we have employed density functional theory to carry out a detailed study of the interaction between this complex metal hydride and the carbon nanomaterial. Considering a stepwise reduction of the hydrogen content in LiBH(4), we find that the presence of C(60) can lead to a substantial reduction of the involved H-removal energies. This effect is explained as a consequence of the interaction between the BH(x)( - ) complex and the C(60) entity.

Role of catalysts in dehydrogenation of MgH2 nanoclusters.

A fundamental understanding of the role of catalysts in dehydrogenation of MgH(2) nanoclusters is provided by carrying out first-principles calculations based on density functional theory. It is shown that the transition metal atoms Ti, V, Fe, and Ni not only lower desorption energies significantly but also continue to attract at least four hydrogen atoms even when the total hydrogen content of the cluster decreases. In particular, Fe is found to migrate from the surface sites to the interior sites during the dehydrogenation process, releasing more hydrogen as it diffuses. This diffusion mechanism may account for the fact that a small amount of catalysts is sufficient to improve the kinetics of MgH(2), which is essential for the use of this material for hydrogen storage in fuel-cell applications.

Significance of self-trapping on hydrogen diffusion.

The diffusion rate of hydrogen in Nb was calculated using ab initio molecular dynamics simulations. At low temperatures the hydrogen is strongly trapped in a local strain field which is caused by the elastic response of the lattice. At elevated temperatures, the residence time (τ) of hydrogen in an interstitial site is not sufficient for fully developing the local strain field. This unbinding of the interstitial hydrogen and the strain field increases the hopping rate (1/τ) at elevated temperatures (>400 K). These results call for a revision of the conceptual framework of diffusion of hydrogen in transition metals at elevated temperatures.

Li-decorated metal-organic framework 5: a route to achieving a suitable hydrogen storage medium.

A significant improvement in molecular hydrogen uptake properties is revealed by our ab initio calculations for Li-decorated metal-organic framework 5. We have found that two Li atoms are strongly adsorbed on the surfaces of the six-carbon rings, one on each side, carrying a charge of +0.9e per Li atom. Each Li can cluster three H(2) molecules around itself with a binding energy of 12 kJ (mol H(2))(-1). Furthermore, we show from ab initio molecular dynamics simulations with a hydrogen loading of 18 H(2) per formula unit that a hydrogen uptake of 2.9 wt % at 200 K and 2.0 wt % at 300 K is achievable. To our knowledge, this is the highest hydrogen storage capacity reported for metal-organic framework 5 under such thermodynamic conditions.

Carbon nanomaterials as catalysts for hydrogen uptake and release in NaAlH4.

A synergistic approach involving experiment and first-principles theory not only shows that carbon nanostructures can be used as catalysts for hydrogen uptake and release in complex metal hydrides such as sodium alanate, NaAlH(4), but also provides an unambiguous understanding of how the catalysts work. Here we show that the stability of NaAlH(4) originates with the charge transfer from Na to the AlH(4) moiety, resulting in an ionic bond between Na(+) and AlH(4)(-) and a covalent bond between Al and H. Interaction of NaAlH(4) with an electronegative substrate such as carbon fullerene or nanotube affects the ability of Na to donate its charge to AlH(4), consequently weakening the Al-H bond and causing hydrogen to desorb at lower temperatures as well as facilitating the absorption of H(2) to reverse the dehydrogenation reaction. In addition, based on our experimental observations and theoretical calculations it appears the curvature of the carbon nanostructure plays a role in the catalytic process. Ab initio molecular dynamics simulation further reveals the time evolution of the charge transfer process.