Exploiting Supramolecular Interactions to Control Isomer Distributions in Reduced-Symmetry [Pd2L4]4+ Cages.

High-symmetry metallosupramolecular architectures (MSAs) have been exploited for a range of applications including molecular recognition, catalysis, and drug delivery. Recently, there have been increasing efforts to enhance those applications by generating reduced-symmetry MSAs. Here we report our attempts to use supramolecular (dispersion and hydrogen-bonding) forces and solvophobic effects to generate isomerically pure [Pd2(L)4]4+ cage architectures from a family of new reduced-symmetry ditopic tripyridyl ligands. The reduced-symmetry tripyridyl ligands featured either solvophilic polyether chains, solvophobic alkyl chains, or amino substituents. We show using NMR spectroscopy, high-performance liquid chromatography, X-ray diffraction data, and density functional theory calculations that the combination of dispersion forces and solvophobic effects does not provide any control of the [Pd2(L)4]4+ isomer distribution with mixtures of all four cage isomers (HHHH, HHHT, cis-HHTT, or trans-HTHT, where H = head and T = tail) obtained in each case. More control was obtained by exploiting hydrogen-bonding interactions between amino units. While the cage assembly with a 3-amino-substituted tripyridyl ligand leads to a mixture of all four possible isomers, the related 2-amino-substituted tripyridyl ligand generated a cis-HHTT cage architecture. Formation of the cis-HHTT [Pd2(L)4]4+ cage was confirmed using NMR studies and X-ray crystallography.

Comparison of Taboo Search Methods for Atomic Cluster Global Optimization with a Basin-Hopping Algorithm.

The basin-hopping algorithm (BHA) allows for the efficient exploration of atomic cluster potential energy surfaces by random perturbations in configuration space, followed by energy minimizations. Here, the taboo search method is incorporated to prevent the search from revisiting recently visited regions of the search space. Two taboo search modes are implemented, one mode resets the search to random coordinates upon encountering the taboo region, while the other simply rejects any proposed move into the taboo region. These two modes are tested and compared on a variety of potential energy surfaces─several clusters where atomic interactions are described by the Lennard-Jones potential, and Au55 where a semi-empirical tight binding potential is used to describe atomic interactions. Some differences in performance between the two taboo search modes were noted for LJ38 and Au55, with the mode that rejects all hops into the taboo region performing better, offering a means to improve the efficiency of the BHA for multifunnel systems. However, both taboo search modes failed to significantly improve performance on multifunnel systems where more than two funnels were present in the system.

Quantitative Assessment of Ligand Substituent Effects on σ- and π-Contributions to Fe-N Bonds in Spin Crossover FeII Complexes.

The effect of para-substituent X on the electronic structure of sixteen tridentate 4-X-(2,6-di(pyrazol-1-yl))-pyridine (bppX ) ligands and the corresponding solution spin crossover [FeII (bppX )2 ]2+ complexes is analysed further, to supply quantitative insights into the effect of X on the σ-donor and π-acceptor character of the Fe-NA (pyridine) bonds. EDA-NOCV on the sixteen LS complexes revealed that neither ΔEorb,σ+π (R2 =0.48) nor ΔEorb,π (R2 =0.31) correlated with the experimental solution T1/2 values (which are expected to reflect the ligand field imposed on the iron centre), but that ΔEorb,σ correlates well (R2 =0.82) and implies that as X changes from EDG→EWG (Electron Donating to Withdrawing Group), the ligand becomes a better σ-donor. This counter-intuitive result was further probed by Mulliken analysis of the NA atomic orbitals: NA (px ) involved in the Fe-N σ-bond vs. the perpendicular NA (pz ) employed in the ligand aromatic π-system. As X changes EDG→EWG, the electron population on NA (pz ) decreases, making it a better π-acceptor, whilst that in NA (px ) increases, making it a better σ-bond donor; both increase ligand field, and T1/2 as observed. In 2016, Halcrow, Deeth and co-workers proposed an intuitively reasonable explanation of the effect of the para-X substituents on the T1/2 values in this family of complexes, consistent with the calculated MO energy levels, that M→L π-backdonation dominates in these M-L bonds. Here the quantitative EDA-NOCV analysis of the M-L bond contributions provides a more complete, coherent and detailed picture of the relative impact of M-L σ-versus π-bonding in determining the observed T1/2 , refining the earlier interpretation and revealing the importance of the σ-bonding. Furthermore, our results are in perfect agreement with the ΔE(HS-LS) vs. σp + (X) correlation reported in their work.

Halimane Diterpenes in the Alpine Daisy Celmisia viscosa: Absolute Configuration, 2,6-Dideoxyhexopyran-3-ulosides, Conformational Flexibility, and Intraspecific Variation.

The methyl-migrated bicyclic skeleton of the halimane diterpenes has been found in a wide range of organisms, including flowering plants, liverworts, marine animals, and bacteria. The discovery of halima-1(10),14-dien-13-ol (3) from the Aotearoa New Zealand endemic alpine daisy Celmisia viscosa is now reported. The full configuration was assigned for the first time by X-ray crystallography, enantiomeric to that of a liverwort isolate. The absolute configuration at C-5 of the halimane is opposite to that at C-5 of the labdane epimanool (1) found in some C. viscosa specimens. Two new 2,6-dideoxyhexopyran-3-uloside halimane derivatives (4 and 5) were also found, and the absolute configuration of 5 was determined by 1H NMR analysis of the Mosher esters. Line broadening in the 13C NMR spectra of these halim-1(10)-enes was due to conformational exchange in the decalin ring A, as shown by molecular modeling and DFT calculations. 1H NMR and GC analyses of leaf extracts of individual plants from across the full geographic range of C. viscosa revealed intraspecific variation of diterpenes: 37 samples had halimadienol as the main diterpene in large amounts and 2 specimens had predominantly epimanool, again in large amounts. Three other viscid (sticky leaved) Celmisia species also contained diterpenes, but none was detectable in four nonviscid Celmisia species.

Development of a Structural Comparison Method to Promote Exploration of the Potential Energy Surface in the Global Optimization of Nanoclusters.

A structural comparison method (SCM) was created to quantify the structural diversity of nanoclusters and was implemented into a global optimization algorithm to evaluate structural diversity between generated clusters on the fly and promote exploration of the potential energy surface. The SCM evaluated topological differences between clusters using the common neighbor analysis and provided a numerical measure of similarity between the two clusters. The SCM was implemented into a genetic algorithm by integrating it into a new structure + energy fitness operator such that structural diversity of clusters in the population and their energies were used to assign fitness values to clusters. The efficiency of the genetic algorithm with this new fitness operator was benchmarked against several Lennard-Jones clusters (LJ38, LJ75, and LJ98) known to be difficult cases for global optimization algorithms. For LJ38 and LJ75, this new structure + energy fitness operator performed equally well or better than the energy fitness operator. However, the efficiency of locating the global minimum of LJ98 decreased using this new structure + energy fitness operator. Further analysis of the genetic algorithm with this fitness operator showed that the algorithm did indeed promote exploration of the PES of LJ98 as desired but hindered refinement of clusters, preventing it from locating the global minimum even if the energy funnel of the global minimum had been located.

A theoretical investigation of 38-atom CuPd clusters: the effect of potential parameterisation on structure and segregation.

Understanding the structure of bimetallic clusters is increasingly important due to their emerging practical applications. Herein we investigate the structure of 38-atom CuPd clusters using a genetic algorithm with cluster energies described by the semi-empirical Gupta potential. Selected clusters are then refined with density functional theory. Three different parameterisations of the Gupta potential are used and their performance assessed to understand what features of bulk and surfaces are necessary to capture for accurate description of small clusters. Three general regions of motif stability exist; for the Pd majority clusters (Pd38 to Cu4Pd34) the truncated octahedron is most stable, while for clusters of intermediate compositions (Cu5Pd33 to Cu25Pd13) a "pancake" icosahedron is most stable, and for the Cu majority clusters (Cu26Pd12 to Cu38) again the truncated octahedron is most stable. CuPd clusters tend to segregate to a Cu-core, Pd-shell structure if possible, and at higher Cu compositions, the Pd segregates to the faces of the cluster. Using multiple parameterisations of the Gupta potential ensures the full variety of possible structures is found, and improves the search for the most stable CuPd clusters.

Quantitative and Chemically Intuitive Evaluation of the Nature of M-L Bonds in Paramagnetic Compounds: Application of EDA-NOCV Theory to Spin Crossover Complexes.

To improve understanding of M-L bonds in 3d transition metal complexes, analysis by energy decomposition analysis and natural orbital for chemical valence model (EDA-NOCV) is desirable as it provides a full, quantitative and chemically intuitive ab initio description of the M-L interactions. In this study, a generally applicable fragmentation and computational protocol was established and validated by using octahedral spin crossover (SCO) complexes, as the transition temperature (T1/2 ) is sensitive to subtle changes in M-L bonding. Specifically, EDA-NOCV analysis of Fe-N bonds in five [FeII (Lazine )2 (NCBH3 )2 ], in both low-spin (LS) and paramagnetic high-spin (HS) states led to: 1) development of a general, widely applicable, corrected M+L6 fragmentation, tested against a family of five LS [FeII (Lazine )3 ](BF4 )2 complexes; this confirmed that three Lazine are stronger ligands (ΔEorb,σ+π =-370 kcal mol-1 ) than 2 Lazine +2 NCBH3 (=-335 kcal mol-1 ), as observed. 2) Analysis of Fe-L bonding on LS→HS, reveals more ionic (ΔEelstat ) and less covalent (ΔEorb ) character (ΔEelstat :ΔEorb 55:45 LS→64:36 HS), mostly due to a big drop in σ (ΔEorb,σ ↓50 %; -310→-145 kcal mol-1 ), and a drop in π contributions (ΔEorb,π ↓90 %; -30→-3 kcal mol-1 ). 3) Strong correlation of observed T1/2 and ΔEorb,σ+π , for both LS and HS families (R2 =0.99 LS, R2 =0.95 HS), but no correlation of T1/2 and ΔΔEorb,σ+π (LS-HS) (R2 =0.11). Overall, this study has established and validated an EDA-NOCV protocol for M-L bonding analysis of any diamagnetic or paramagnetic, homoleptic or heteroleptic, octahedral transition metal complex. This new and widely applicable EDA-NOCV protocol holds great promise as a predictive tool.

Does the Reaction of Cyclopropyl Acid Chlorides and Imines To Form 1,3-Oxazin-4-enone Heterocycles Proceed via a Ketene or an N-Acyl-iminium Mechanism?

The mechanism of the reaction between cyclopropyl acid chlorides and imines to form 1,3-oxazin-4-enones was probed through physical and computational experiments. The data gathered strongly support the reaction proceeding through an N-acyl iminium intermediate mechanism rather than a ketene intermediate mechanism.

Predictable Substituent Control of CoIII/II Redox Potential and Spin Crossover in Bis(dipyridylpyrrolide)cobalt Complexes.

A family of five easily prepared tridentate monoanionic 2,5-dipyridyl-3-(R1)-4-(R2)-pyrrolide anions (dppR1,R2)-, varying in the nature of the R1 and R2 substituents [R1, R2 = CN, Ph; CO2Et, CO2Et; CO2Me, 4-Py; CO2Me, Me; Me, Me], has been used to generate the analogous family of neutral [CoII(dppR1,R2)2] complexes, two of which are structurally characterized at both 100 and 298 K. Both the oxidation and spin states of these complexes can be switched in response to appropriate external stimuli. All complexes, except [CoII(dppMe,Me)2], exhibit gradual spin crossover (SCO) in the solid state, and SCO activity is observed for three complexes in CDCl3 solution. The cobalt(II) centers in the low spin (LS) complexes are Jahn-Teller tetragonally compressed along the pyrrolide-Co-pyrrolide axis. The complexes in their high spin (HS) states are more distorted than in the LS states, as is also usually the case for SCO active iron(II) complexes. The reversible CoIII/II redox potentials are predictably tuned by choice of substituents R1 and R2, from -0.95 (Me,Me) to -0.45 (CN,Ph) V vs Fc+/Fc, with a linear correlation observed between E1/2(CoIII/II) and the Swain-Lupton parameters of the pyrrolide substituents.

Elucidation of temperature-programmed desorption of high-coverage hydrogen on Pt(211), Pt(221), Pt(533) and Pt(553) based on density functional theory calculations.

In this work we compute high-coverage hydrogen adsorption energies and geometries on the stepped platinum surfaces Pt(211) and Pt(533) which contain a (100)-step type and the Pt(221) and Pt(553) surface with a (111) step edge. We discuss these results in relation to ultra-high-vacuum temperature programmed desorption (TPD) data to elucidate the origin of the desorption features. Our results indicated that on surfaces with a (100)-step type, two distinct ranges of adsorption energy for the step and terrace are observed, which mirrors the TPD spectra for which we find a clear separation of the desorption peaks. For the (111) step type, the TPD spectra show much less separation of the step and terrace features, which we assign to the low individual adsorption energies for H atoms on this step edge. From our results we obtain a much clearer understanding of the surface-hydrogen bonding at high coverages and the origin of the different TPD features present for the two step types studied.

Mechanistic Evaluation of Bioorthogonal Decaging with trans-Cyclooctene: The Effect of Fluorine Substituents on Aryl Azide Reactivity and Decaging from the 1,2,3-Triazoline.

Bioorthogonal prodrug activation/decaging strategies need to be selective, rapid and release the drug from the masking group upon activation. The rates of the 1,3-dipolar cycloaddition between a trans-cyclooctene (TCO) and a series of fluorine-substituted azido-PABC self-immolative spacers caging two model drugs, and subsequent release from the 1,2,3-triazoline are reported. As the number of fluorine substituents on the PABC linker increases from one to four, the rate of cycloaddition increases by almost one order of magnitude. Using a combination of fluorescence, 1H/19F NMR, and computational experiments, we have been able to determine how substituents on the PABC ring can influence the degradation rates and also the product distribution of the 1,2,3-triazoline. We have also been able to determine how these substituents influence the rate of imine hydrolysis and 1,6-self-immolation decaging rates of the generated anilines. The NMR and computational studies demonstrate that fluorine substituents on the aromatic ring lower the transition state energy required for converting the triazoline to the imine or aziridine intermediates via extrusion of diatomic nitrogen, and that in the case of a tetrafluoro substituted aromatic ring, it is the imine hydrolysis and 1,6-self-immolation that is rate-limiting. This knowledge further enhances the understanding of factors which influence the stability of triazolines, and enables potential applications of fluorinated aromatics, in particular, perfluorinated aromatics, in synthetic chemistry and sustained-release drug delivery systems.

A Simple Method of Predicting Spin State in Solution.

A simple method, using density functional theory (DFT), of predicting spin-state in advance of synthesis is reported. Specifically, an excellent correlation is observed between the switching temperatures (T1/2) in CDCl3 solution of five spin-crossover (SCO)-active [FeII(Lazine)2(NCBH3)2] complexes and the DFT-calculated (and observed) 15N NMR chemical shift (δNA) of the five different azine-substituted 1,2,4-triazole ligands employed, Lazine = 4-(4-methylphenyl)-3-phenyl-5-(azine)-1,2,4-triazole, where azine = pyridine, pyridazine, 4-pyrimidine, pyrazine, and 2-pyrimidine. To test the generality of this finding, DFT was also employed to readily predict the δNA values for a family of 16 literature ligands, known as bppX,Y [X,Y-substituted 2,6-(pyrazol-1-yl)pyridines], which have produced 16 SCO-active [FeII(bppX,Y)2](Z)2 complexes (Z = BF4 or in one case PF6) in (CD3)2CO solution: again an excellent correlation was found between the computed δNA and the observed T1/2. These correlations represent a key advance in the field, as they allow a simple DFT calculation on a modified ligand to be used to reliably predict, before synthesis of the ligand or complex, the T1/2 that would result from that modification. Achieving such easily predictable tuning of T1/2, and hence of spin-state, is a significant step forward in the field of SCO and also has big implications in many other fields in which spin-state is key, including catalysis, metallo-enzyme modeling studies, and host-guest chemistry.

Self-Assembly of Cyclohelicate [M3 L3 ] Triangles Over [M4 L4 ] Squares, Despite Near-Linear Bis-terdentate L and Octahedral M.

Self-assembly of 1:1:2 MII (BF4 )2 (M=Zn or Fe), pyrazine-2,5-dicarbaldehyde (1) and 2-(2-aminoethyl)pyridine gave trimetallic triangle architectures rather than the anticipated tetrametallic [2×2] squares. Options for the nontrivial synthesis of 1 are considered, and synthetic details provided for both preferred routes. Rare cyclohelicate triangle architectures are observed for the pair of structurally characterized yellow-brown [Zn3 L3 ](BF4 )6 and dark green [Fe3 L3 ](BF4 )6 complexes of the neutral bis-terdentate Schiff base L. In order to form these pyrazine-edged triangles, the octahedral metal ions-with all 6 N-donors provided by the terdentate binding pockets of two L-are located 0.4-0.5 Šout of the plane of the bridging pyrazines, towards the center of the triangle. Density functional theory calculations confirm that simple particle counting entropic arguments, which predict triangles over squares, are correct here, with the triangles shown to be energetically favored over the corresponding squares. However, importantly, DFT analysis of these and related triangle versus square systems also show that vibrational contributions to entropy dominate and may significantly influence the preferred architecture, such that simple particle counting cannot in general be reliably employed to predict the observed architecture.

Heterometallic [MnPtn(L)2n](x+) Macrocycles from Dichloromethane-Derived Bis-2-pyridyl-1,2,3-triazole Ligands.

There is continued interest in the assembly of heterometallic metallosupramolecular assemblies due to the potential for new structural types and/or interesting chemical and physical properties. Two novel methylene-linked bis-2-pyridyl-1,2,3-triazole ligands (1a, 1b) were synthesized through CuAAC "click" chemistry between dichloromethane (DCM), sodium azide, and the appropriate ethynylpyridine, giving the shortest possible "regular" bis-pyridyl-1,2,3-triazole. The first example of bis-coordination of two 2-pyridyl-1,2,3-triazoles (Rpytri) around one Pt(II) center is reported, giving exclusive formation of head-to-tail [Pt(Rpytri)2](2+) complexes with vacant binding sites suitable for complexation with other metals and the formation of heterometallic assemblies. The concentration-dependent formation of an equilibrium mixture of a heterometallic [Pd2Pt2(L)4](8+) [4 + 4] square and [Pd3Pt3(L)6](12+) [6 + 6] hexagon was observed, and at lower concentrations ([reactants] = 1.5 mM) hexamer formation was negligible. The [Pt(L)2](2+) building block could also be utilized in the synthesis of a concentration-independent [Cu2Pt2(L)4](6+) metallomacrocyle. These compounds were characterized with (1)H, (13)C, and (1)H DOSY NMR spectroscopies, elemental analysis, mass spectrometry, and in some cases X-ray crystallography.

Enabling electrochemical reduction of nitrogen to ammonia at ambient conditions through rational catalyst design.

Commercial design of a sustainable route for on-site production of ammonia represents a potential economic and environmental breakthrough. In an analogous process to the naturally occurring enzymatic mechanism, synthesis of ammonia could be achieved in an electrochemical cell, in which electricity would be used to reduce atmospheric nitrogen and water into ammonia at ambient conditions. To date, such a process has not been realized due to slow kinetics and low faradaic efficiencies. Although progress has been made in this regard, at present there exists no device that can produce ammonia efficiently from air and water at room temperature and ambient pressure. In this work, a scheme is presented in which electronic structure calculations are used to screen for catalysts that are stable, active and selective towards N2 electro-reduction to ammonia, while at the same time suppressing the competing H2 evolution reaction. The scheme is applied to transition metal nitride catalysts. The most promising candidates are the (100) facets of the rocksalt structures of VN and ZrN, which show promise of producing ammonia in high yield at low onset potentials.

Predictable electronic tuning of FeII and RuII complexes via choice of azine: correlation of ligand pKa with Epa(MIII/II) of complex.

Five new mononuclear ruthenium(II) tris-ligated complexes have been synthesised, varying through the choice of azine in the family of 3-azinyl-4-(4-methylphenyl)-5-phenyl-4H-1,2,4-triazole ligands (Lazine): [Ru(Lpyridine)](PF6)2 (1), [Ru(Lpyridazine)](PF6)2 (2), [Ru(L4-pyrimidine)](PF6)2 (3), [Ru(Lpyrazine)](PF6)2 (4), [Ru(L2-pyrimidine)](PF6)2 (5). Three of them, 1·2MeCN·Et2O, 3·2MeCN·Et2O and 4·2MeCN, have been structurally characterised, confirming the presence of the meridional isomer, as was previously reported for the FeII analogues. Cyclic voltammetry studies, in dry CH3CN vs. Ag/0.01 M AgNO3, show that all five RuII complexes undergo a reversible RuIII/RuII process, with the midpoint potential (Em) increasing from 0.87 to 1.18 V as the azine is changed: pyridine < pyridazine < 2-pyrimidine < 4-pyrimidine < pyrazine. A strong inverse linear correlation (R2 = 0.98) is found between the RuIII/RuII redox potential and the calculated HOMO orbital energies, which is consistent with the expectation that it is easier to oxidise (lower Em) a metal ion with a higher HOMO orbital energy. The same trend was reported earlier for the family of analogous FeII complexes, albeit at lower values of Em in all cases. In addition, the ionisation potentials of the RuII complexes, as well as those of the other group 8 analogues (FeII and OsII), showed a linear relationship with Epa. As the MIII/II redox potentials of a family of complexes has been previously reported to correlate with ligand pKa values, a computational protocol to calculate, in silico, the pKa of the Lazine family of ligands was developed. A strong linear relationship was found between the readily calculated pKa of the Lazine ligand and the Epa of the MII complex, for all three families of complexes (R2 = 0.98).

Calculations of the effect of catalyst size and structure on the electrocatalytic reduction of CO2 on Cu nanoclusters.

The structure and catalytic properties of Cu nanoclusters of sizes between 55 and 147 atoms were examined to understand if small Cu clusters could provide enhancement over traditional catalysts for the electrocatalysis of CO2 to CO and carbon-based fuels, such as CH4 and CH3OH, compared to bulk Cu surfaces and large Cu nanoparticles. Clusters studied included Cu55, Cu78, Cu101, Cu124, and Cu147, the structures of which were determined using global optimisation. The majority of Cu clusters examined were icosahedral, including the perfect closed-shell, partial-shell, elongated and distorted icosahedral clusters. Free energy diagrams for the reduction of CO2 showed the potential required for the formation of CO is notably smaller for all cluster sizes considered, relative to Cu(111). Less variation is observed for the limiting potential for the formation of CH4 and CH3OH. However, it was found that clusters that are either a distorted motif or contain vacancy defects yielded the best activity and provide an interesting synthesis target for future experiments.

Dynamic sampling of liquid metal structures for theoretical studies on catalysis.

Liquid metals have recently emerged as promising catalysts that can outcompete their solid counterparts for many reactions. Although theoretical modelling is extensively used to improve solid-state catalysts, there is currently no way to capture the interactions of adsorbates with a dynamic liquid metal. We propose a new approach based on ab initio molecular dynamics sampling of an adsorbate on a liquid catalyst. Using this approach, we describe time-resolved structures for formate adsorbed on liquid Ga-In, and for all intermediates in the methanol oxidation pathway on Ga-Pt. This yields a range of accessible adsorption energies that take into account the at-temperature motion of the liquid metal. We find that a previously proposed pathway for methanol oxidation on Ga-Pt results in unstable intermediates on a dynamic liquid surface, and propose that H desorption must occur during the path. The results showcase a more accurate way to treat liquid metal catalysts in this emerging field.

Protecting group free glycosylation: one-pot stereocontrolled access to 1,2-trans glycosides and (1→6)-linked disaccharides of 2-acetamido sugars.

Unprotected 2-acetamido sugars may be directly converted into their oxazolines using 2-chloro-1,3-dimethylimidazolinium chloride (DMC), and a suitable base, in aqueous solution. Freeze drying and acid catalysed reaction with an alcohol as solvent produces the corresponding 1,2-trans-glycosides in good yield. Alternatively, dissolution in an aprotic solvent system and acidic activation in the presence of an excess of an unprotected glycoside as a glycosyl acceptor, results in the stereoselective formation of the corresponding 1,2-trans linked disaccharides without any protecting group manipulations. Reactions using aryl glycosides as acceptors are completely regioselective, producing only the (1→6)-linked disaccharides.

Size-dependent trends in the hydrogen evolution activity and electronic structure of MoS2 nanotubes.

The thermodynamics of hydrogen evolution on MoS2 nanotubes is studied for the first time using periodic density functional theory calculations to obtain hydrogen adsorption free energies (ΔG Hads ) on pristine nanotubes and those with S-vacancy defects. Armchair and zigzag MoS2 nanotubes of different diameters, ranging from 12 to 22 Å, are examined. The H adsorption energy is observed to become more favourable (lower ΔG Hads ) as nanotube diameter decreases, with ΔG Hads values ranging from 1.82 to 1.39 eV on the pristine nanotubes, and from 0.03 to -0.30 eV at the nanotube S-vacancy defect sites. An ideal thermoneutral ΔG Hads value of nearly 0 eV is observed at the S-vacancy site on nanotubes around 20 to 22 Å in diameter. For the pristine nanotubes, density of states calculations reveal that electron transfer from S to Mo occurs during H adsorption, and the energy gap between these two states yields a highly reliable linear correlation with ΔG Hads , where a smaller gap leads to a more favourable hydrogen adsorption. For the S-vacancy defect site the H adsorption resembles that on a pure metallic surface, meaning that a traditional d-band centre model can be applied to explain the trends in ΔG Hads . A linear relation between the position of the Mo d-states and ΔG Hads is found, with d-states closer to the Fermi level leading to strong hydrogen adsorption. Overall this work highlights the relevance of MoS2 nanotubes as promising hydrogen evolution catalysts and explains trends in their activity using the energies of the electronic states involved in binding hydrogen.