Researcher profile

Nicholas M. Harrison

Nicholas M. Harrison contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - Baseline
4works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

4 published item(s)

preprint2022arXiv

A transferable prediction model of molecular adsorption on metals based on adsorbate and substrate properties

Surface adsorption is one of the fundamental processes in numerous fields, including catalysis, environment, energy and medicine. The development of an adsorption model which provides an effective prediction of binding energy in minutes has been a long term goal in surface and interface science. The solution has been elusive as identifying the intrinsic determinants of the adsorption energy for various compositions, structures and environments is non-trivial. We introduce a new and flexible model for predicting adsorption energies to metal substrates. The model is based on easily computed, intrinsic properties of the substrate and adsorbate. It is parameterised using machine learning based on first-principles calculations of probe molecules (e.g., H$_2$O, CO$_2$, O$_2$, N$_2$) adsorbed to a range of pure metal substrates. The model predicts the computed dissociative adsorption energy to metal surfaces with a correlation coefficient of 0.93 and a mean absolute error of 0.77 eV for the large database of molecular adsorption energies provided by \textit{Catalysis-Hub.org} which have a range of 15 eV. As the model is based on pre-computed quantities it provides near-instantaneous estimates of adsorption energies and it is sufficiently accurate to eliminate around 90\% of candidates in screening study of new adsorbates. The model, therefore, significantly enhances current efforts to identify new molecular coatings in many applied research fields.

preprint2013arXiv

Comment on "First-principles study of the influence of (110)-oriented strain on the ferroelectric properties of rutile TiO2" [arXiv:1106.2820]

In a recent article, Grünebohm et al. [Phys. Rev. B 84 132105 (2011), arXiv:1106.2820] report that they fail to reproduce the A2u ferroelectric instability of TiO2 in the rutile structure calculated with density functional theory within the PBE-GGA approximation by Montanari et al. [Chem. Phys. Lett 364, 528 (2002)]. We demonstrate that this disagreement arises from an erroneous treatment of Ti 3s and 3p semi-core electrons as core in their calculations. Fortuitously the effect of the frozen semi-core pseudopotential cancels the phonon instability of the PBE exchange-correlation, and the combination yields phonon frequencies similar to the LDA harmonic values. Grünebohm et al. also attempted and failed to reproduce the soft acoustic phonon mode instability under (110) strain reported by Mitev et al. [Phys. Rev. B 81 134303 (2010)]. For this mode the combination of PBE-GGA and frozen semi-core yields a small imaginary frequency of 9.8i. The failure of Grünebohm et al. to find this mode probably arose from numerical limitations of the geometry optimization approach in the presence of a shallow double well potential; the optimization method is not suitable for locating such instabilities.

preprint2013arXiv

Linear-scaling time-dependent density-functional theory in the linear response formalism

We present an implementation of time-dependent density-functional theory (TDDFT) in the linear response formalism enabling the calculation of low energy optical absorption spectra for large molecules and nanostructures. The method avoids any explicit reference to canonical representations of either occupied or virtual Kohn-Sham states and thus achieves linear-scaling computational effort with system size. In contrast to conventional localised orbital formulations, where a single set of localised functions is used to span the occupied and unoccupied state manifold, we make use of two sets of in situ optimised localised orbitals, one for the occupied and one for the unoccupied space. This double representation approach avoids known problems of spanning the space of unoccupied Kohn-Sham states with a minimal set of localised orbitals optimised for the occupied space, while the in situ optimisation procedure allows for efficient calculations with a minimal number of functions. The method is applied to a number of medium sized organic molecules and a good agreement with traditional TDDFT methods is observed. Furthermore, linear scaling of computational cost with system size is demonstrated on a system of carbon nanotubes.

preprint2012arXiv

Ab initio complex band structure of conjugated polymers: Effects of hydrid DFT and GW schemes

The non-resonant tunneling regime for charge transfer across nanojunctions is critically dependent on the so-called β parameter, governing the exponential decay of the current as the length of the junction increases. For periodic materials, this parameter can be theoretically evaluated by computing the complex band structure (CBS) -- or evanescent states -- of the material forming the tunneling junction. In this work we present the calculation of the CBS for organic polymers using a variety of computational schemes, including standard local, semilocal, and hybrid-exchange density functionals, and many-body perturbation theory within the GW approximation. We compare the description of localization and β parameters among the adopted methods and with experimental data. We show that local and semilocal density functionals systematically underestimate the β parameter, while hybrid-exchange schemes partially correct for this discrepancy, resulting in a much better agreement with GW calculations and experiments. Self-consistency effects and self-energy representation issues of the GW corrections are discussed together with the use of Wannier functions to interpolate the electronic band-structure.