Anhydrous Proton Conduction Through a Chemically Robust Electrolyte Enabling a High-Temperature Non-Precious Metal Catalyzed Fuel Cell.

Fuel cells offer great promise for portable electricity generation, but their use is currently limited by their low durability, excessive operating temperatures, and expensive precious metal electrodes. It is therefore essential to develop fuel cell systems that can perform effectively using more robust electrolyte materials, at reasonable temperatures, with lower-cost electrodes. Recently, proton exchange membrane fuel cells have attracted attention due to their generally favorable chemical stability and quick start-up times. However, in most membrane materials, water is required for proton conduction, severely limiting operational temperatures. Here, for the first time it is demonstrated that when acidified, PAF-1 can conduct protons at high temperatures, via a unique framework diffusion mechanism. It shows that this acidified PAF-1 material can be pressed into pellets with high proton conduction properties even at high temperatures and pellet thickness, highlighting the processibility, and ease of use of this material. Furthermore, a fuel cell is shown with high power density output is possible using a non-precious metal copper electrode. Acid-doped PAF-1 therefore represents a significant step forward in the potential for a broad-purpose fuel cell due to it being cheap, robust, efficient, and easily processible.

Artificial synthesis of covalent triazine frameworks for local structure and property determination.

Here we show an 'artificial synthesis' method for covalent triazine framework (CTF) materials, enabling localised structural features to be incorporated that result directly from the acid-catalysed synthetic protocol that would otherwise not be captured. This advancement will enable prediction and design of new CTF materials with targeted properties.

Investigation of structure and dynamics in a photochromic molecular crystal by NMR crystallography.

A photochromic anil, N-(3,5-di-t-butylsalicylidene)-4-amino-pyridine, has been studied by single-crystal X-ray diffraction, multinuclear magic-angle spinning NMR, and first-principles density functional theory (DFT) calculations. Interpretation of the solid-state NMR data on the basis of calculated chemical shifts confirms the structure is primarily composed of molecules in the ground-state enol tautomer, whereas thermally activated cis-keto and photoisomerised trans-keto states exist as low-level defects with populations that are too low to detect experimentally. Variable temperature 13 C NMR data reveal evidence for solid-state dynamics, which is found to be associated with fast rotational motion of t-butyl groups and 180° flips of the pyridine ring, contrasting the time-averaged structure obtained by X-ray diffraction. Comparison of calculated chemical shifts for the full crystal structure and an isolated molecule also reveals evidence for an intermolecular hydrogen bond involving the pyridine ring and an adjacent imine carbon, which facilitates the flipping motion. The DFT calculations also reveal that the molecular conformation in the crystal structure is very close to the energetic minimum for an isolated molecule, indicating that the ring dynamics arise as a result of considerable steric freedom of the pyridine ring and which also allows the molecule to adopt a favourable conformation for photochromism.

Rapidly accessible "click" rotaxanes utilizing a single amide hydrogen bond templating motif.

The synthesis of hydrogen bond templated rotaxanes using the CuAAC click reaction has been achieved in yields of up to 47%, employing near stoichiometric equivalents of macrocycle and readily prepared azide and alkyne half-axle components. Interlocked structure formation has been confirmed by NMR spectroscopy and mass spectrometry. Density functional theory calculations support 1H NMR spectroscopic analysis that the macrocycle resides over the amide of the axle component, rather than the newly formed triazole, as a result of more favourable hydrogen bond interactions.

The synthesis of a pyridine-N-oxide isophthalamide rotaxane utilizing supplementary amide hydrogen bond interactions.

The synthesis of a pyridine-N-oxide containing rotaxane, not requiring an additional ionic template, has been achieved in 32% yield. Successful rotaxane formation is dependent upon the structure of the isophthalamide macrocycle used, an observation which has been rationalised by a combination of NMR spectroscopy, X-ray crystallography and computational modelling.

Syntheses, structures, and comparison of the photophysical properties of cyclometalated iridium complexes containing the isomeric 1- and 2-(2'-pyridyl)pyrene ligands.

Two cyclometalated iridium complexes of the form IrL2(acac) have been synthesized, where L is either of the isomeric ligands 1- or 2-(2'-pyridyl)pyrene (1-pypyrH or 2-pypyrH). These complexes have been investigated in terms of their photophysical behavior and, although both complexes exhibit similar pure radiative lifetimes, they have substantially different observed phosphorescence lifetimes and quantum yields. Moreover, the observed phosphorescence lifetimes and quantum yields of both complexes, as well as the absorption spectra of Ir(1-pypyr)2(acac), exhibit a strong solvent dependence, while there is essentially no solvatochromism in the emission spectra of either complex. Single-crystal X-ray diffraction studies of both ligands and both iridium complexes reveal structural differences between the two isomers. The crystal structures of the ligands, supported by density functional theory (DFT) modeling, show that a twist is present between the pyridyl and pyrenyl rings in 1-pypyrH, but is absent in 2-pypyrH, which leads to the requirement for more unusual cyclometalation conditions for 1-pypyrH. Furthermore, it is suggested that the strained structure of Ir(1-pypyr)2(acac) provides access to a facile nonradiative excited state deactivation pathway, which leads to the higher value of knr for this isomer. DFT, TD-DFT, and ΔSCF calculations have been conducted to investigate further the photophysical properties of the complexes, allowing a detailed comparison of the two isomers. We find that Tamm-Dancoff Approximation TD-DFT with the CAM-B3LYP functional provides the best agreement between experimentally and theoretically determined transition energies, performing better than the more common combination of TD-DFT with B3LYP, the reasons for which are outlined. We also highlight some difficulties with performing optimization calculations on oxidized complexes to assess electrochemical data.

Overcoming low orbital overlap and triplet instability problems in TDDFT.

Low orbital overlap and triplet instability problems in time-dependent density functional theory (TDDFT) are investigated for a new benchmark set, encompassing challenging singlet and triplet excitation energies of local, charge-transfer, and Rydberg character. The low orbital overlap problem is largely overcome for both singlet and triplet states by the use of a Coulomb-attenuated functional. For all the categories of functional considered, however, errors associated with triplet instability problems plague high overlap excitations, as exemplified by the excited states of acenes and polyacetylene oligomers. Application of the Tamm-Dancoff approximation reduces these errors for both singlet and triplet states, while leaving low-overlap excitations unaffected. The study illustrates the synergy between overlap and stability and highlights the success of a combined, Coulomb-attenuated Tamm-Dancoff approach.

On the evaluation of the non-interacting kinetic energy in density functional theory.

The utility of both an orbital-free and a single-orbital expression for computing the non-interacting kinetic energy in density functional theory is investigated for simple atomic systems. The accuracy of both expressions is governed by the extent to which the Kohn-Sham equation is solved for the given exchange-correlation functional and so special attention is paid to the influence of finite Gaussian basis sets. The orbital-free expression is a statement of the virial theorem and its accuracy is quantified. The accuracy of the single-orbital expression is sensitive to the choice of Kohn-Sham orbital. The use of particularly compact orbitals is problematic because the failure to solve the Kohn-Sham equation exactly in regions where the orbital has decayed to near-zero leads to unphysical behaviour in regions that contribute to the kinetic energy, rendering it inaccurate. This problem is particularly severe for core orbitals, which would otherwise appear attractive due to their formally nodeless nature. The most accurate results from the single-orbital expression are obtained using the relatively diffuse, highest occupied orbitals, although special care is required at orbital nodes.

Shielding constants and chemical shifts in DFT: influence of optimized effective potential and Coulomb-attenuation.

The influence of the optimized effective potential (OEP) and Coulomb-attenuation on shielding constants and chemical shifts is investigated for three disparate categories of molecule: main group, hydrogen bonded, and transition metal systems. Expanding the OEP in the orbital basis leads to physically sensible exchange-correlation potentials; OEP generalized gradient approximation results provide some indication of the accuracy of the expansion. OEP uncoupled magnetic parameters from representative hybrid and Coulomb-attenuated functionals can be a dramatic improvement over conventional results; both categories yield similar accuracy. Additional flexibility is introduced by expanding the OEP in an extensive even-tempered basis set, but this leads to the well-known problem of unphysical, oscillatory potentials. Smooth potentials are recovered through the use of a smoothing norm, but deficiencies in the procedure are highlighted for transition metal complexes. The study reiterates the importance of the OEP procedure in magnetic response calculations using orbital-dependent functionals, together with the need for careful attention to ensure physically sensible potentials. It also illustrates the utility of Coulomb-attenuated functionals for computing short-range molecular properties.

Mechanistic insights into the cytochrome P450-mediated oxidation and rearrangement of littorine in tropane alkaloid biosynthesis.

During the biosynthesis of certain tropane alkaloids, littorine (1) is rearranged to hyoscyamine (3). Recent evidence indicates that this isomerisation is a two-step process in which the first step is an oxidation/rearrangement to give hyoscyamine aldehyde (2). This step is catalysed by CYP80F1, a cytochrome P450 enzyme, which was recently identified from the plant Hyoscyamus niger; CYP80F1 also catalyses the hydroxylation of littorine at the 3'-position. The mechanisms of the reactions catalysed by CYP80F1 were probed with synthetic deutero and arylfluoro analogues of 1. Measurement of the primary kinetic isotope effects indicates that C3' hydrogen abstraction is the rate-limiting step for the oxidation/rearrangement of natural littorine, and for the 3'-hydroxylation reaction of the unnatural S enantiomer of littorine. The character of the intermediates in the oxidation/rearrangement and hydroxylation reaction was probed with the use of arylfluorinated analogues of (R)-littorine (natural stereoisomer) and (S)-littorine (unnatural stereoisomer) as substrates for CYP80F1. The relative conversions of ortho-, meta- and para-fluorolittorine analogues were used to obtain information on the likely intermediacy of either a benzylic radical or benzylic carbocation intermediate. The data suggest that hydroxylation takes place via a benzylic carbocation intermediate, whereas the product profile arising from rearrangement is more consistent with a benzylic radical intermediate.

Time-dependent density functional theory calculations of near-edge X-ray absorption fine structure with short-range corrected functionals.

We report calculations of core excitation energies and near-edge X-ray absorption fine structure (NEXAFS) spectra computed with time-dependent density functional theory (TDDFT). TDDFT with generalized gradient approximation and standard hybrid exchange-correlation functionals is known to underestimate core excitation energies. This failure is shown to be associated with the self-interaction error at short interelectronic distances. Short-range corrected hybrid functionals are shown to reduce the error in the computed core excitation energies for first and second row nuclei in a range of molecules to a level approaching that observed in more traditional excited states calculations in the ultraviolet region. NEXAFS spectra computed with the new functionals agree well with experiment and the pre-edge features in the NEXAFS spectra of plastocyanin are correctly predicted.

Adiabatic connection forms in density functional theory: H2 and the He isoelectronic series.

Full configuration interaction (FCI) data are used to quantify the accuracy of approximate adiabatic connection (AC) forms in describing two challenging problems in density functional theory--the singlet ground state potential energy curve of H(2) in a restricted formalism and the energies of the helium isoelectronic series, H(-) to Ne(8+). For H(2), an exponential-based form yields a potential energy curve that is virtually indistinguishable from the FCI curve, eliminating the unphysical barrier to dissociation observed previously with a [1,1]-Pade-based form and with the random phase approximation. For the helium isoelectronic series, the Pade-based form gives the best overall description, followed by the exponential form, with errors that are orders of magnitude smaller than those from a standard hybrid functional. Particular attention is paid to the limiting behavior of the AC forms with increasing bond distance in H(2) and increasing atomic number in the isoelectronic series; several forms describe both limits correctly. The study illustrates the very high quality results that can be obtained using exchange-correlation functionals based on simple AC forms, when near-exact data are used to determine the parameters in the forms.

Fractional Electron Loss in Approximate DFT and Hartree-Fock Theory.

Plots of electronic energy vs electron number, determined using approximate density functional theory (DFT) and Hartree-Fock theory, are typically piecewise convex and piecewise concave, respectively. The curves also commonly exhibit a minimum and maximum, respectively, in the neutral → anion segment, which lead to positive DFT anion HOMO energies and positive Hartree-Fock neutral LUMO energies. These minima/maxima are a consequence of using basis sets that are local to the system, preventing fractional electron loss. Ground-state curves are presented that illustrate the idealized behavior that would occur if the basis set were to be modified to enable fractional electron loss without changing the description in the vicinity of the system. The key feature is that the energy cannot increase when the electron number increases, so the slope cannot be anywhere positive, meaning frontier orbital energies cannot be positive. For the convex (DFT) case, the idealized curve is flat beyond a critical electron number such that any additional fraction of an electron added to the system is unbound. The anion HOMO energy is zero. For the concave (Hartree-Fock) case, the idealized curve is flat up to some critical electron number, beyond which it curves down to the anion energy. A minimum fraction of an electron is required before any binding occurs, but beyond that, the full fraction abruptly binds. The neutral LUMO energy is zero. Approximate DFT and Hartree-Fock results are presented for the F → F(-) segment, and results approaching the idealized behavior are recovered for highly diffuse basis sets. It is noted that if a DFT calculation using a highly diffuse basis set yields a negative LUMO energy then a fraction of an electron must bind and the electron affinity must be positive, irrespective of whether an electron binds experimentally. This is illustrated by calculations on Ne → Ne(-).

Assessment of Tuning Methods for Enforcing Approximate Energy Linearity in Range-Separated Hybrid Functionals.

A range of tuning methods, for enforcing approximate energy linearity through a system-by-system optimization of a range-separated hybrid functional, are assessed. For a series of atoms, the accuracy of the frontier orbital energies, ionization potentials, electron affinities, and orbital energy gaps is quantified, and particular attention is paid to the extent to which approximate energy linearity is actually achieved. The tuning methods can yield significantly improved orbital energies and orbital energy gaps, compared to those from conventional functionals. For systems with integer M electrons, optimal results are obtained using a tuning norm based on the highest occupied orbital energy of the M and M + 1 electron systems, with deviations of just 0.1-0.2 eV in these quantities, compared to exact values. However, detailed examination for the carbon atom illustrates a subtle cancellation between errors arising from nonlinearity and errors in the computed ionization potentials and electron affinities used in the tuning.

A new water···Na+ coordination motif in an unexpected diatrizoic acid disodium salt crystal form.

The disodium salt of diatrizoic acid crystallises as a tetragonal channel hydrate structure. One of the incorporated water molecules is coordinated to three individual sodium cations in a unique geometry.

Influence of Triplet Instabilities in TDDFT.

Singlet and triplet vertical excitation energies from time-dependent density functional theory (TDDFT) can be affected in different ways by the inclusion of exact exchange in hybrid or Coulomb-attenuated/range-separated exchange-correlation functionals; in particular, triplet excitation energies can become significantly too low. To investigate these issues, the explicit dependence of excitation energies on exact exchange is quantified for four representative molecules, paying attention to the effect of constant, short-range, and long-range contributions. A stability analysis is used to verify that the problematic TDDFT triplet excitations can be understood in terms of the ground state triplet instability problem, and it is proposed that a Hartree-Fock stability analysis should be used to identify triplet excitations for which the presence of exact exchange in the TDDFT functional is undesirable. The use of the Tamm-Dancoff approximation (TDA) significantly improves the problematic triplet excitation energies, recovering the correct state ordering in benzoquinone; it also affects the corresponding singlet states, recovering the correct state ordering in naphthalene. The impressive performance of the TDA is maintained for a wide range of molecules across representative functionals.

TDDFT diagnostic testing and functional assessment for triazene chromophores.

A simple diagnostic test based on orbital overlap [M. J. G. Peach et al., J. Chem. Phys., 2008, 128, 044118] may be used to help judge the reliability of excitation energies in time-dependent density functional theory (TDDFT) when using generalized gradient approximation (GGA) and hybrid functionals. Orbital plots are used to illustrate the test for a model tripeptide and for 4-(N,N-dimethylamino)benzonitrile, which are representative of systems containing low- and high-overlap charge-transfer excitations. The scheme is then applied to a series of triazene chromophores in solvent, highlighting the relationship between overlap and oscillator strength and its implications for theoretical absorption spectra. No low-overlap excitations are observed with a hybrid functional; a single one is identified using a GGA. To assess the diagnostic test and to judge functional performance, gas phase triazene TDDFT excitations are compared with correlated ab initio values. The diagnostic test correctly identifies two low-overlap problematic GGA excitations. However, it does not identify another problematic excitation where the electron is excited to a spatially extended orbital, which necessarily has reasonable overlap with the occupied orbital; an improved diagnostic quantity is required for such cases. The best agreement between TDDFT and correlated ab initio excitations is obtained using a Coulomb-attenuated functional; the errors are significantly smaller than from the GGA and hybrid functionals. The study provides further support for the high quality excitations from Coulomb-attenuated functionals, negating the need for diagnostic tests.

Excitation energies in density functional theory: an evaluation and a diagnostic test.

Electronic excitation energies are determined using the CAM-B3LYP Coulomb-attenuated functional [T. Yanai et al. Chem. Phys. Lett. 393, 51 (2004)], together with a standard generalized gradient approximation (GGA) and hybrid functional. The degree of spatial overlap between the occupied and virtual orbitals involved in an excitation is measured using a quantity Lambda, and the extent to which excitation energy errors correlate with Lambda is quantified. For a set of 59 excitations of local, Rydberg, and intramolecular charge-transfer character in 18 theoretically challenging main-group molecules, CAM-B3LYP provides by far the best overall performance; no correlation is observed between excitation energy errors and Lambda, reflecting the good quality, balanced description of all three categories of excitation. By contrast, a clear correlation is observed for the GGA and, to a lesser extent, the hybrid functional, allowing a simple diagnostic test to be proposed for judging the reliability of a general excitation from these functionals--when Lambda falls below a prescribed threshold, excitations are likely to be in very significant error. The study highlights the ambiguous nature of the term "charge transfer," providing insight into the observation that while many charge-transfer excitations are poorly described by GGA and hybrid functionals, others are accurately reproduced.