Search

Understanding hydrogen-metal interactions is important in various fields of surface science, including the aqueous corrosion of metals. The interaction between atomic H and a Mg surface is a key process for the formation of sub-surface Mg hydride, which may play an important role in Mg aqueous corrosion. In the present work, we performed first-principles Density Functional Theory (DFT) calculations to study the mechanisms for hydrogen adsorption and crystalline Mg hydride formation under aqueous conditions. The Electron Localisation Function (ELF) is found to be a promising indicator for predicting stable H adsorption in the Mg surface. It is found that H adsorption and hydride layer formation is dominated by high ELF adsorption sites. Our calculations suggest that the on-surface adsorption of atomic H, OH radicals and atomic O could enhance the electron localisation at specific sites in the sub-surface region, thus forming effective H traps locally. This is predicted to result in the formation of a thermodynamically stable sub-surface hydride layer, which is a potential precursor of the crucial hydride corrosion product of magnesium.

We investigate the phase behavior of model ternary triacylglycerol blends, comprising triolein (C57H104O6, OOO), tripalmitin (C51H98O6, PPP) and tristearin (C57H110O6, SSS), building upon extensive characterisation of single and binary mixtures, in order to rigorously map the thermal transitions of model natural ‘fats’. A combination of calorimetry, X-ray diffraction, and FTIR spectroscopy is employed to determine crystallisation and melting temperatures and identify the corresponding phases in the complex ternary system. We recover the eutectic behaviour of SSS-PPP blends and the invariability of OOO neat transitions, and resolve the complex β′ + β ternary surface, reflecting the roles of unsaturation and polymorphism of its constituents. Our results provide a representation of the OOO:PPP:SSS:temperature phase behaviour into a triangular prism, consistent with binary pair-wise data, which can inform a range of food science, cosmetic, pharmaceutical and cleaning applications that depend strongly on the physical-chemistry of such multicomponent ‘triglycerides’.

In this work, we develop a photocatalyst wherein nitrogen and phosphorus co-doped carbon quantum dots are scaffolded onto TiO2 nanoparticles (NPCQD/TiO2), denoted as NPCT hereafter. The developed NPCT photocatalyst exhibits an enhanced visible light photocatalytic hydrogen production of 533 μmol h−1 g−1 compared to nitrogen doped CQD/TiO2 (478 μmol h−1 g−1), phosphorus doped CQD/TiO2 (451 μmol h−1 g−1) and pure CQD/TiO2 (427 μmol h−1 g−1) photocatalysts. The enhanced photocatalytic activity of the NPCT photocatalyst is attributed to the excellent synergy between NPCQDs and TiO2 nanoparticles, which results in the creation of virtual energy levels, a decrease in work function and suppressed recombination rates, thereby increasing the lifetime of photogenerated electrons. A detailed mechanism is proposed for the enhancement in visible light hydrogen production by the NPCT photocatalyst from the experimental results, Mott–Schottky plots and ultraviolet photoelectron spectroscopy results. Further, first-principles density functional theory (DFT) simulations are carried out which predict the decrease in the work function and band gap, and the increase in the density of states of NPCT as the factors responsible for the observed enhancement in visible light photocatalytic hydrogen production.

Room Temperature Ionic Liquids (RTILs) attract significant interest in nanotribology. However, their microscopic lubrication mechanism is still under debate. Here, using non-equilibrium molecular dynamics simulations, we investigate the lubrication performance of ultra-thin (0.5 GPa, while bi-layers sustain pressures in the 1–2 GPa range. The compression of these films results in monolayers that can sustain loads of several GPa without significant loss in their lubrication performance. Surprisingly, in such ultra-thin films the imidazolium rings show higher orientational in-plane disorder, with and the rings adopting a tilted orientation with respect to the gold surface. The friction force and friction coefficient of the monolayers depends strongly on the structure of the gold plates, with the friction coefficient being four times higher for monolayers confined between Au(100) surfaces than for more compact Au(111) surfaces. We show that the general behaviour described here is independent of whether the metallic surfaces are modelled as polarizable or non-polarizable surfaces and speculate on the nature of this unexpected conclusion.

Activation of CO2 is the first step towards its reduction to more useful chemicals. Here we systematically investigate the CO2 activation mechanism on Cu3X (X is a first-row transition metal atom) using density functional theory computations. The CO2 adsorption energies and the activation mechanisms depend strongly on the selected dopant. The dopant electronegativity, the HOMO-LUMO gap and the overlap of the frontier molecular orbitals control the CO2 dissociation efficiency. Our calculations reveal that early transition metal-doped (Sc, Ti, V) clusters exhibit a high CO2 adsorption energy, a low activation barrier for its dissociation, and a facile regeneration of the clusters. Thus, early transition metal-doped copper clusters, particularly Cu3Sc, may be efficient catalysts for the carbon capture and utilization process. ispartof: PHYSICAL CHEMISTRY CHEMICAL PHYSICS vol:23 issue:38 pages:21738-21747 ispartof: location:England status: published

The local activation overpotential describes the electrostatic potential shift away from equilibrium at an electrode/electrolyte interface. This electrostatic potential is not entirely satisfactory for describing the reaction kinetics of a mixed ionic–electronic conducting (MIEC) solid-oxide cell (SOC) electrode where charge transfer occurs at the electrode–gas interface. Using the theory of the electrostatic potential at the MIEC–gas interface as an electrochemical driving force, charge transfer at the ceria–gas interface has been modelled based on the intrinsic dipole potential of the adsorbate. This model gives a physically meaningful reason for the enhancement in electrochemical activity of a MIEC electrode as the steam and hydrogen pressure is increased in both fuel cell and electrolysis modes. This model was validated against operando XPS data from previous literature to accurately predict the outer work function shift of thin film Sm0.2Ce0.8O1.9 in a H2/H2O atmosphere as a function of overpotential.

WO3 photoanodes offer rare stability in acidic media, but are limited by their selectivity for oxygen evolution over parasitic side reactions, when employed in photoelectrochemical (PEC) water splitting. Herein, this is remedied via the modification of nanostructured WO3 photoanodes with surface decorated PdO as an oxygen evolution co-catalyst (OEC). The photoanodes and co-catalyst particles are grown using an up-scalable aerosol assisted chemical vapour deposition (AA-CVD) route, and their physical properties characterised by X-ray diffraction (XRD), Raman spectroscopy, X-ray photoelectron spectroscopy (XPS), scanning electron microscopy (SEM), high-resolution transmission electron microscopy (HR-TEM) and UV-vis absorption spectroscopy. Subsequent PEC and transient photocurrent (TPC) measurements showed that the use of a PdO co-catalyst dramatically increases the faradaic efficiency (FE) of water oxidation from 52% to 92%, whilst simultaneously enhancing the photocurrent generation and charge extraction rate. The Pd oxidation state was found to be critical in achieving these notable improvements to the photoanode performance, which is primarily attributed to the higher selectivity towards oxygen evolution when PdO is used as an OEC and the formation of a favourable junction between WO3 and PdO, that drives band bending and charge separation.

Ab initio electronic excited state calculations are necessary for the quantitative study of photochemical reactions, but their accurate computation on classical computers is plagued by prohibitive resource scaling. The Variational Quantum Deflation (VQD) is an extension of the quantum-classical Variational Quantum Eigensolver (VQE) algorithm for calculating electronic excited state energies, and has the potential to address some of these scaling challenges using quantum computers. However, quantum computers available in the near term can only support a limited number of quantum circuit operations, so reducing the quantum computational cost in VQD methods is critical to their realisation. In this work, we investigate the use of adaptive quantum circuit growth (ADAPT-VQE) in excited state VQD calculations, a strategy that has been successful previously in reducing the resources required for ground state energy VQE calculations. We also invoke spin restrictions to separate the recovery of eigenstates with different spin symmetry to reduce the number of calculations and accumulation of errors in computing excited states. We created a quantum eigensolver emulation package - Quantum Eigensolver Building on Achievements of Both quantum computing and quantum chemistry (QEBAB) - for testing the proposed adaptive procedure against two existing VQD methods that use fixed-length quantum circuits: UCCGSD-VQD and k-UpCCGSD-VQD. For a lithium hydride test case we found that the spin-restricted adaptive growth variant of VQD uses the most compact circuits out of the tested methods by far, consistently recovers adequate electron correlation energy for different nuclear geometries and eigenstates while isolating the singlet and triplet manifold. This work is a further step towards developing techniques which improve the efficiency of hybrid quantum algorithms for excited state quantum chemistry, opening up the possibility of exploiting real quantum computers for electronic excited state calculations sooner than previously anticipated.

The expansion of lithium-ion batteries from consumer electronics to larger-scale transport and energy storage applications has made understanding the many mechanisms responsible for battery degradation increasingly important. The literature in this complex topic has grown considerably; this perspective aims to distil current knowledge into a succinct form, as a reference and a guide to understanding battery degradation. Unlike other reviews, this work emphasises the coupling between the different mechanisms and the different physical and chemical approaches used to trigger, identify and monitor various mechanisms, as well as the various computational models that attempt to simulate these interactions. Degradation is separated into three levels: the actual mechanisms themselves, the observable consequences at cell level called modes and the operational effects such as capacity or power fade. Five principal and thirteen secondary mechanisms were found that are generally considered to be the cause of degradation during normal operation, which all give rise to five observable modes. A flowchart illustrates the different feedback loops that couple the various forms of degradation, whilst a table is presented to highlight the experimental conditions that are most likely to trigger specific degradation mechanisms. Together, they provide a powerful guide to designing experiments or models for investigating battery degradation.

Biomolecular conformational transitions are usually modeled as barrier crossings in a free energy landscape. The transition paths connect two local free energy minima and transition path times (TPT) are the actual durations of the crossing events. The simplest model employed to analyze TPT and to fit empirical data is that of a stochastic particle crossing a parabolic barrier. Motivated by some disagreement between the value of the barrier height obtained from the TPT distributions as compared to the value obtained from kinetic and thermodynamic analyses, we investigate here TPT for barriers which deviate from the symmetric parabolic shape. We introduce a continuous set of potentials, that starting from a parabolic shape, can be made increasingly asymmetric by tuning a single parameter. The TPT distributions obtained in the asymmetric case are very well-fitted by distributions generated by parabolic barriers. The fits, however, provide values for the barrier heights and diffusion coefficients which deviate from the original input values. We show how these findings can be understood from the analysis of the eigenvalues spectrum of the Fokker-Planck equation and highlight connections with experimental results., Comment: We dedicate this paper to the memory of our dear friend and colleague Carlo Vanderzande

We use a variety of computational methods to characterize and compare the hydrogen atom transfer (HAT) and epoxidation reaction pathways for oxidation of cyclohexene by an iron(iv)-oxo complex. Previous B3LYP calculations have led to predictions that both alcohol (from the HAT route) and epoxide should be formed in similar amounts, which was not in agreement with experiment where only the HAT product was observed. We show here that ab initio calculations which can take both static and dynamic correlation into account are needed to explain the experimentally observed dominance of the HAT process. Since these systems do not have very strong multireference character we have also tested different flavours of local coupled cluster methods. We suggest that further improvements are necessary before they can provide highly accurate results for these systems. ispartof: PHYSICAL CHEMISTRY CHEMICAL PHYSICS vol:22 issue:41 pages:23908-23919 ispartof: location:England status: published

Planar and quasi-planar boron clusters with a disk-like shape are investigated in search of common bonding characteristics. Methods used involve molecular orbital calculations based on Density Functional Theory (DFT), and valence bond partitioning using Adaptive Natural Density Partitioning (AdNDP) analysis. For high-symmetry cases the proposed bonding schemes are confirmed using the group-theoretical induction method. The focus is on the electron occupation of delocalized in-plane 3-center and 4-center bonds. For disks consisting of concentric rings this inner electron count is found to be equal to a multiple of the vertex count of the inner polygon. For two concentric rings the multiplying factor is four, for three concentric rings it is eight. The appropriate bonding schemes are presented which explain these results. Some giant clusters with two hexagonal holes are also discussed. ispartof: PHYSICAL CHEMISTRY CHEMICAL PHYSICS vol:21 issue:2 pages:729-735 ispartof: location:England status: published

The radiative cooling of highly excited carbon cluster cations of sizes N = 8, 10, 13-16 has been studied in an electrostatic storage ring. The cooling rate constants vary with cluster size from a maximum at N = 8 of 2.6 × 104 s-1 and a minimum at N = 13 of 4.4 × 103 s-1. The high rates indicate that photon emission takes place from electronically excited ions, providing a strong stabilizing cooling of the molecules. ispartof: PHYSICAL CHEMISTRY CHEMICAL PHYSICS vol:21 issue:3 pages:1587-1596 ispartof: location:England status: published

We report the high-pressure properties of two heterometallic perovskite-type metal-organic frameworks (MOFs) templated by dimethylammonium (NH2(CH3)2, DMA+) with the general formula [DMA]MI0.5CrIII0.5(HCOO)3, where MI = Na+ (DMANaCr) and K+ (DMAKCr). The high-pressure Raman scattering studies show crystal instabilities in the 4.0-4.4 GPa and 2.0-2.5 GPa ranges for DMANaCr and DMAKCr, respectively. The mechanism is similar in the two compounds and involves strong deformation of the metal-formate framework, especially pronounced for the subnetwork of CrO6 octahedra, accompanied by substantial compressibility of the DMA+ cations. Comparison with previous high-pressure Raman studies of sodium-chromium heterometallic MOFs show that the stability depends on the templated cation and increases as follows: ammonium < imidazolium < DMA+. Density functional theory (DFT) calculations are performed to get a better understanding of the structural properties leading to the existence of phase transitions. We calculate the energy of the hydrogen bonds (HBs) between the DMA+ cation and the metal formate cage, revealing a stronger interaction in the DMAKCr compound due to a HB arrangement that primarily involves the energetically preferred bonding to KO6 octahedra. This material however also has a smaller structural tolerance factor (TF) and a higher vibrational entropy than DMANaCr. This indicates a more flexible crystal structure, explaining the lower phase transition pressure, as well as the previously observed phase transition at 190 K, which is absent in the DMANaCr compound. The DFT high-pressure simulations show the largest contraction to be along the trigonal axis, leading to a minimal distortion of the HBs formed between the DMA+ cations and the metal-formate sublattice.

The recent discovery of extremely red-shifted chlorophyll f pigments in both photosystem I (PSI) and photosystem II has led to the conclusion that chlorophyll f plays a role not only in the energy transfer, but also in the charge separation processes [Nürnberg et al., Science, 2018, 360, 1210–1213]. We have employed ultrafast transient infrared absorption spectroscopy to study the contribution of far-red light absorbing chlorophyll f to energy transfer and charge separation processes in far-red light-grown PSI (FRL-PSI) from the cyanobacterium Chroococcidiopsis thermalis PCC 7203. We compare the kinetics and spectra of FRL-grown PSI excited at 670 nm and 740 nm wavelengths to those of white light-grown PSI (WL-PSI) obtained at 675 nm excitation. We report a fast decay of excited state features of chlorophyll a and complete energy transfer from chlorophyll a to chlorophyll f in FRL-PSI upon 670 nm excitation, as indicated by a frequency shift in a carbonyl absorption band occurring within a 1 ps timescale. While the WL-PSI measurements support the assignment of initial charge separation to A−1+˙A0−˙ [Di Donato et al., Biochemistry, 2011, 50, 480–490] from the kinetics of a distinct cation feature at 1710 cm−1, in the case of FRL-PSI, small features at 1715 cm−1 from the chlorophyll cation are present from sub-ps delays instead, supporting the replacement of the A−1 pigment with chlorophyll f. Comparisons of nanosecond spectra show that charge separation proceeds with 740 nm excitation, which selectively excites chlorophyll f, and modifications in specific carbonyl absorption bands assigned to P700+˙ minus P700 and A1−˙ minus A1 indicate dielectric differences of FRL-PSI compared to WL-PSI in one or both of the two electron transfer branches of FRL-PSI.

The lowest-lying isomer of the B50 boron cluster is confirmed to have a quasi-planar shape with two hexagonal holes. By applying a topological (leap-frog) dual operation followed by boron capping, we demonstrated that such a quasi-planar structure actually comes from the smallest elongated B102-, and its high thermodynamic stability is due to its inherent disk aromaticity arising from its 32 valent π electrons that fully occupy a disk configuration of [(1σ)2(1π)4(1δ)4(2σ)2(1φ)4(2π)4(1γ)4(2δ)4(1η)4]. The aromatic character of the quasi-planar B50 is further supported by a strong diatropic magnetic current flow. The sudden appearance of a quasi-planar B50 again points out that the growth pattern of pure boron clusters is still far from being completely understood. ispartof: PHYSICAL CHEMISTRY CHEMICAL PHYSICS vol:21 issue:13 pages:7039-7044 ispartof: location:England status: published

Physical chemistry, chemical physics 21(7), 3721 - 3733 (2019). doi:10.1039/C8CP05063E, Published by RSC Publ., Cambridge

The skin surface, our first barrier against the external environment, is covered by the sebum oil, a lipid film composed of sebaceous and epidermal lipids, which is important in the regulation of the hydration level of our skin. Here, we investigate the pathways leading to the transfer of epidermal lipids from the skin lipid bilayer to the sebum. We show that the sebum triglycerides, a major component of sebum, interact strongly with the epidermal lipids and extract them from the bilayer. Using microsecond time scale molecular dynamics simulations, we identify and quantify the free energy associated with the skin lipid extraction process.

We present results of combined experimental photoelectron spectroscopy and theoretical modeling studies of solvated dicarboxylate species (-O2C(CH2)2CO2-) in complex with Na+ and K+ metal cations. These ternary clusters serve as simple models for the investigation of aqueous ion/solute specific effects that play an important role in biological systems. The experimental characterization of these systems was performed in the presence of up to six solvating waters. In both Na+ and K+ cases, we observe the presence of one major broad band that gradually shifts to higher electron binding energy (EBE) with an increasing number of waters. In the Na+ case further detailed analysis of experimental spectra was performed using ab initio calculations. In particular, we have identified the structures of the lowest energy clusters whose EBE values match well the major band in the experimental spectra. Our results show that evolution of an aqueous solvation shell emphasizes the coordination of the negatively charged carboxylate groups accompanied by simultaneous interaction with metal cations. Calculations also indicate that in the solvation range investigated experimentally (up to 6 waters), Na+ retains direct contact with the dicarboxylate species, i.e. a contact ion-pair (CIP) complex. Preliminary modeling studies show evidence of an alternative solvent separated ion-pair complex once the solvation range approaches 8 waters, however its energy still remains above that of (∼7-8 kcal/ mol-1) the CIP complex. At a higher number of waters (n = 3 for Na+ and n = 5 for K+), the experimental spectra also show the development of a weak low energy band. Its origin cannot be precisely identified. Our calculations in the Na+ case point out the existence of a quaternary complex consisting of Na+, H2O, OH- and a singly protonated dicarboxylate anion (HO2C(CH)2CO2-). Such a complex appears to be stabilized in the solvation range corresponding to the appearance of the low EBE band and does match its peak, even though the energy of such a complex is fairly high compared to the ternary structure. ispartof: PHYSICAL CHEMISTRY CHEMICAL PHYSICS vol:20 issue:46 pages:29051-29060 ispartof: location:England status: published

Zinc oxide (ZnO) nanorods (NRs) have been demonstrated as a promising platform for enhanced fluorescence-based sensing. It is, however, desirable to achieve a tuneable fluorescence enhancement with these platforms so that the fluorescence output can be adjusted based on the real need. Here we show that the fluorescence enhancement can be tuned by changing the diameter of the ZnO nanorods, simply controlled by potassium chloride (KCl) concentration during synthesis, using arrays of previously developed aligned NRs (a.k.a. aligned NR forests) and nanoflowers (NFs). Combining the experimental results obtained from ZnO nanostructures with controlled morphology and computer-aided verification, we show that the fluorescence enhancement factor increases when ZnO NRs become thicker. The fluorescence enhancement factor of NF arrays is shown to have a much stronger dependency on the rod diameter than that of aligned NR arrays. We prove that the morphology of nanostructures, which can be controlled, can be an important factor for fluorescence enhancement. Our (i) effort towards understanding the structure–property relationships of ZnO nanostructured arrays and (ii) demonstration on tuneable fluorescence enhancement by nanostructure engineering can provide some guidance towards the rational design of future fluorescence amplification platforms potentially for bio-sensing.

Deep eutectic solvents (DESs) constitute a new class of ionic solvents that has been developing at a fast pace in recent years. Since these solvents are commonly suggested as green alternatives to organic solvents, it is important to understand their physical properties. In particular, polarity plays an important role in solvation phenomena. In this work, the polarity of different families of DESs was studied through solvatochromic responses of UV-vis absorption probes. Kamlet-Taft α, β, π* and ETN parameters were evaluated using different solvatochromic probes, as 2,6-dichloro-4-(2,4,6-triphenyl-N-pyridino)-phenolate (Reichardt's betaine dye 33), 4-nitroaniline, and N,N-diethyl-4-nitroaniline for several families of DESs based on cholinium chloride, dl-menthol and a quaternary ammonium salt ([N4444]Cl). In addition, a study to understand the difference in polarity properties between DESs and the corresponding ILs, namely ILs based on cholinium cation and carboxylic acids as anions ([Ch][Lev], [Ch][Gly] and [Ch][Mal]), was carried out. The chemical structure of the hydrogen bond acceptor (HBA) in a DES clearly controls the dipolarity/polarizability afforded by the DES. Moreover, Kamlet-Taft parameters do not vary much within the family, but they differ among families based on different HBA, either for DESs containing salts ([Ch]Cl or [N4444]Cl) or neutral compounds (dl-menthol). A substitution of the HBD was also found to play an important role in solvatochromic probe behaviour for all the studied systems.

The atomic contributions to valence electronic structure for 37 ionic liquids (ILs) are identified using a combination of variable photon energy XPS, resonant Auger electron spectroscopy (RAES) and a subtraction method. The ILs studied include a diverse range of cationic and anionic structural moieties. We introduce a new parameter for ILs, the energy difference between the energies of the cationic and anionic highest occupied fragment orbitals (HOFOs), which we use to identify the highest occupied molecular orbital (HOMO). The anion gave rise to the HOMO for 25 of the 37 ILs studied here. For 10 of the ILs, the energies of the cationic and anionic HOFOs were the same (within experimental error); therefore, it could not be determined whether the HOMO was from the cation or the anion. For two of the ILs, the HOMO was from the cation and not from the anion; consequently it is energetically more favourable to remove an electron from the cation than the anion for these two ILs. In addition, we used a combination of area normalisation and subtraction of XP spectra to produce what are effectively XP spectra for individual ions; this was achieved for seven cations and 14 anions.

Enhanced conductivity in YSZ films has been of substantial interest over the last decade. In this paper we examine the effects of substrate lattice mismatch and film thickness on the strain in YSZ films and the resultant effect on the conductivity. 8 mol% YSZ films have been grown on MgO, Al2O3, LAO and NGO substrates, thereby controlling the lattice mismatch at the film/substrate interface. The thickness of the films was varied to probe the interfacial contribution to the transport properties, as measured by impedance spectroscopy and tracer diffusion. No enhancement in the transport properties of any of the films was found over single crystal values, and instead the effects of lattice strain were found to be minimal. The interfaces of all films were more resistive due to a heterogeneous distribution of grain boundaries, and no evidence for enhanced transport down dislocations was found.

Efficient enhancement of photoluminescence in rare-earth activated upconversion materials is of great significance for their practical applications in various fields. In this work{,} three-dimensional mesoporous gold films were fabricated by a low-cost and facile dealloying approach to improve the upconversion photoluminescence efficiency. The mesoporous Au films exhibit good chemical stability{,} large-area uniformity and abundant distribution of porous nanospaces. Varying the time of the dealloying process leads to modification of pore size distribution{,} surface roughness and residual Ag content{,} resulting in effective tuning the wavelength of the broadband localized surface plasmon resonance (LSPR). Enhancement factors were identified to be a function of the dealloying time. With the optimized upconversion photoluminescence enhancement a 41-fold increase was achieved with the mesoporous gold substrate which had been dealloyed for 8 days. These results pave the way to overcome the limitation of poor upconversion efficiency for widespread practical applications in life sciences and energy applications.

The principal difference between 1-benzyl-3-methyl-imidazolium triflimide [BzC1im][NTf2] and an equimolar mixture of benzene and dimethylimidazolium triflimide [C1C1im][NTf2] is that in the former the benzene moieties are tied to the imidazolium ring, while in the latter they move independently. We use femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy (OHD-RIKES) and molecular simulations to explore some properties of these two systems. The Kerr spectra show small differences in the spectral densities; the simulations also show very similar environments for both the imidazolium rings and the phenyl or benzene parts of the molecules. The low frequency vibrational densities of states are also similar in the model systems. In order to perform the simulations we developed a model for the [BzC1im](+) cation and found that the barriers to rotation of the two parts of the molecule are low.

Hybrid perovskite solar cells were first reported in 2009 by Miyasaka and co-workers (J. Am. Chem. Soc., 2009,131, 6050) and sparked an exciting change in photovoltaics research, a field that has been growing rapidly ever since.

In the past decade, ionic liquids (ILs) have been the focus of intensive research regarding their use as potential and alternative solvents in many chemical applications. Targeting their effectiveness, recent investigations have attempted to establish polarity scales capable of ranking ILs according to their chemical behaviours. However, some major drawbacks have been found since polarity scales only report relative ranks because they depend on the set of probe dyes used, and they are sensitive to measurement conditions, such as purity levels of the ILs and procedures employed. Due to all these difficulties it is of crucial importance to find alternative and/or predictive methods and to evaluate them as a priori approaches capable of providing the chemical properties of ILs. Furthermore, the large number of ILs available makes their experimental characterization, usually achieved by a trial and error methodology, burdensome. In this context, we firstly evaluated COSMO-RS, COnductor-like Screening MOdel for Real Solvents, as an alternative tool to estimate the hydrogen-bond basicity of ILs. After demonstrating a straight-line correlation between the experimental hydrogen-bond basicity values and the COSMO-RS hydrogen-bonding energies in equimolar cation-anion pairs, an extended scale for the hydrogen-bond accepting ability of IL anions is proposed here. This new ranking of the ILs' chemical properties opens the possibility to pre-screen appropriate ILs (even those not yet synthesized) for a given task or application.