Source author record

Philip E. Hoggan

Philip E. Hoggan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

6 published item(s)

preprint2022arXiv

Quantum Monte Carlo method for metal-film catalysis: water addition to carbon monoxide adsorbed on Pt/Al(111), a route to hydrogen

Hydrogen production as a clean, sustainable replacement for fossil fuels is gathering pace. Doubling the capacity of Paris-CDG airport has been halted, even with the upcoming Olympic Games, until hydrogen powered planes can be used. It is thus timely to work on catalytic selective hydrogen production and optimise catalyst structure. Over 90 % of all chemical manufacture uses a solid catalyst. This work describes adsorption of carbon-monoxide (CO) on platinum thin films, supported by cheap Al(111). CO reacts with water to produce hydrogen (water-gas shift). Quantum Monte Carlo methods are the only ones accurate enough to investigate the early steps of this catalysed reaction at close-packed Pt/Al(111). Many chemical reactions involve bond-dissociation. This process is often the key to rate-limiting reaction steps at solid surfaces. Since bond-breaking is poorly described by Hartree-Fock and DFT methods, our embedded active site approach is used This work demonstrates a novel Quantum Monte Carlo (QMC) methodology. The water-gas shift reaction step studied is water addition to CO pre-adsorbed on a Pt-monolayer supported by Al(111). The water molecule is only partially dissociated. Its oxygen atom binds to CO giving adsorbed COOH and Pt-H. This concerted addition is rate-limiting. In subsequent steps, the adsorbed formate species (with acidic hydrogen) decomposes to carbon dioxide and, after proton migration to Pt-H, the clean product H$_2$ is obtained. The QMC activation barrier found is 64.8 $\pm$ 1.5 kJ/mol. Thus, QMC is shown to be encouraging for investigating similar catalytic systems.

preprint2022arXiv

Quantum Monte Carlo with ground-state input to investigate methane reactions on supported metal film catalysts

Nowadays, there is pressing demand for sustainable energy sources, or clean and 'green' fuel and hydrogen is a perfect candidate. It can be made by dissociating methane with the energy input compensated by metal-hydrogen bond formation. Industry uses Nickel supported on $γ$-alumina for this rate-limiting step. This work shows how to use Quantum Monte Carlo calculations to take (user-guided) information summarising chemistry of heterogeneous catalysis from a generic Jastrow factor, that allows electron correlation to be included in the density. This, in turn, gives accurate values of properties, in particular the activation barrier for stretching one methane C-H bond, before adsorbing CH3- and H- on vicinal nickel atoms in the close-packed (111) plane. The CH3- co-ordinates to a hollow site defined by three Ni atoms and, new to this work, another Ni-atom co-ordinates the H-atom from the dissociation.

preprint2014arXiv

Performance of numerical approximation on the calculation of two-center two-electron integrals over non-integer Slater-type orbitals using elliptical coordinates

The two-center two-electron Coulomb and hybrid integrals arising in relativistic and nonrelativistic ab-initio calculations of molecules are evaluated over the non-integer Slater-type orbitals via ellipsoidal coordinates. These integrals are expressed through new molecular auxiliary functions and calculated with numerical Global-adaptive method according to parameters of non-integer Slatertype orbitals. The convergence properties of new molecular auxiliary functions are investigated and the results obtained are compared with results found in the literature. The comparison for two-center twoelectron integrals is made with results obtained from one-center expansions by translation of wavefunction to same center with integer principal quantum number and results obtained from the Cuba numerical integration algorithm, respectively. The procedures discussed in this work are capable of yielding highly accurate two-center two-electron integrals for all ranges of orbital parameters.

preprint2011arXiv

Quantum Monte Carlo facing the Hartree-Fock symmetry dilemma: The case of hydrogen rings

When using Hartree-Fock (HF) trial wave functions in quantum Monte Carlo calculations, one faces, in case of HF instabilities, the HF symmetry dilemma in choosing between the symmetry-adapted solution of higher HF energy and symmetry-broken solutions of lower HF energies. In this work, we have examined the HF symmetry dilemma in hydrogen rings which present singlet instabilities for sufficiently large rings. We have found that the symmetry-adapted HF wave function gives a lower energy both in variational Monte Carlo and in fixed-node diffusion Monte Carlo. This indicates that the symmetry-adapted wave function has more accurate nodes than the symmetry-broken wave functions, and thus suggests that spatial symmetry is an important criterion for selecting good trial wave functions.

preprint2010arXiv

How specific exponential type orbitals recently became a viable basis set choice in NMR shielding tensor calculation

This paper advocates use of the atomic orbitals which have direct physical interpretation, i.e. Coulomb Sturmians and hydrogen-like orbitals. They are exponential type orbitals (ETOs). Their radial nodes are shown to be essential in obtaining accurate nuclear shielding tensors for NMR work. The present work builds on a 2003 French PhD and many numerical results were published by 2007. The improvements in this paper are noteworthy, the key being the actual basis function choice. Until 2008, their products on different atoms were difficult to manipulate for the evaluation of two-electron integrals. Coulomb resolutions provide an excellent approximation that reduces these integrals to a sum of one-electron overlap-like integral products that each involve orbitals on at most two centers. Such two-center integrals are separable in prolate spheroidal co-ordinates. They are thus readily evaluated. Only these integrals need to be re-evaluated to change basis functions. In this paper, a review of the translation procedures for Slater type orbitals (STO) and for Coulomb Sturmians follows that of the more recent application to ETOs of a particularly convenient Coulomb resolution.

preprint2010arXiv

Quantum Monte Carlo calculations of electronic excitation energies: the case of the singlet $n \to π^*$ (CO) transition in acrolein

We report state-of-the-art quantum Monte Carlo calculations of the singlet $n \to π^*$ (CO) vertical excitation energy in the acrolein molecule, extending the recent study of Bouabça {\it et al.} [J. Chem. Phys. {\bf 130}, 114107 (2009)]. We investigate the effect of using a Slater basis set instead of a Gaussian basis set, and of using state-average versus state-specific complete-active-space (CAS) wave functions, with or without reoptimization of the coefficients of the configuration state functions (CSFs) and of the orbitals in variational Monte Carlo (VMC). It is found that, with the Slater basis set used here, both state-average and state-specific CAS(6,5) wave functions give an accurate excitation energy in diffusion Monte Carlo (DMC), with or without reoptimization of the CSF and orbital coefficients in the presence of the Jastrow factor. In contrast, the CAS(2,2) wave functions require reoptimization of the CSF and orbital coefficients to give a good DMC excitation energy. Our best estimates of the vertical excitation energy are between 3.86 and 3.89 eV.