Researcher profile

M. J. Paterson

M. J. Paterson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
2topics
3close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

6 published item(s)

preprint2014arXiv

Applying Monte Carlo configuration interaction to transition metal dimers: exploring the balance between static and dynamic correlation

We calculate potential curves for transition metal dimers using Monte Carlo configuration interaction (MCCI). These results, and their associated spectroscopic values, are compared with experimental and computational studies. The multireference nature of the MCCI wavefunction is quantified and we estimate the important orbitals. We initially consider the ground state of the chromium dimer. Next we calculate potential curves for Sc$_{2}$ where we contrast the lowest triplet and quintet states. We look at the molybdenum dimer where we compare non-relativistic results with the partial inclusion of relativistic effects via effective core potentials, and report results for scandium nickel.

preprint2014arXiv

Characterizing a configuration interaction excited state using natural transition geminals

We introduce natural transition geminals as a means to qualitatively understand a transition where double excitations are important. The first two $A_{1}$ singlet states of the CH cation are used as an initial example. We calculate these states with configuration interaction singles (CIS) and state-averaged Monte Carlo configuration interaction (SA-MCCI). For each method we compare the important natural transition geminals with the dominant natural transition orbitals. We then compare SA-MCCI and full configuration interaction (FCI) with regards to the natural transition geminals using the beryllium atom. We compare using the natural transition geminals with analyzing the important configurations in the CI expansion to give the dominant transition for the beryllium atom and the carbon dimer. Finally we calculate the natural transition geminals for two electronic excitations of formamide.

preprint2014arXiv

State-averaged Monte Carlo configuration interaction applied to electronically excited states

We introduce state-averaging into the method of Monte Carlo configuration interaction (SA-MCCI) to allow the stable and efficient calculation of excited states. We show that excited potential curves for H$_{3}$, including a crossing with the ground state, can be accurately reproduced using a small fraction of the FCI space. A recently introduced error measure for potential curves [J. P. Coe and M. J. Paterson, J. Chem. Phys., 137, 204108 (2012)] is shown to also be a fair approach when considering potential curves for multiple states. We demonstrate that potential curves for LiF using SA-MCCI agree well with the FCI results and the avoided crossing occurs correctly. The seam of conical intersections for CH$_{2}$ found by Yarkony [J. Chem. Phys., 104, 2932 (1996)] is used as a test for SA-MCCI and we compare potential curves from SA-MCCI with FCI results for this system for the first three triplet states. We then demonstrate the improvement from using SA-MCCI on the dipole of the $2$ $^{1}A_{1}$ state of carbon monoxide. We then look at vertical excitations for small organic molecules up to the size of butadiene where the SA-MCCI energies and oscillator strengths are compared with CASPT2 values [M. Schreiber, M. R. Silva-Junior, S. P. A. Sauer, and W. Thiel, J. Chem. Phys., 128, 134110 (2008)]. We finally see if the SA-MCCI results for these excitation energies can be improved by using MCCIPT2 with approximate natural orbitals when the PT2 space is not onerously large.

preprint2013arXiv

Development of Monte Carlo configuration interaction: Natural orbitals and second-order perturbation theory

Approximate natural orbitals are investigated as a way to improve a Monte Carlo configuration interaction (MCCI) calculation. We introduce a way to approximate the natural orbitals in MCCI and test these and approximate natural orbitals from MP2 and QCISD in MCCI calculations of single-point energies. The efficiency and accuracy of approximate natural orbitals in MCCI potential curve calculations for the double hydrogen dissociation of water, the dissociation of carbon monoxide and the dissociation of the nitrogen molecule are then considered in comparison with standard MCCI when using full configuration interaction as a benchmark. We also use the method to produce a potential curve for water in an aug-cc-pVTZ basis. A new way to quantify the accuracy of a potential curve is put forward that takes into account all of the points and that the curve can be shifted by a constant. We adapt a second-order perturbation scheme to work with MCCI (MCCIPT2) and improve the efficiency of the removal of duplicate states in the method. MCCIPT2 is tested in the calculation of a potential curve for the dissociation of nitrogen using both Slater determinants and configuration state functions.

preprint2013arXiv

Monte Carlo configuration interaction applied to multipole moments, ionisation energies and electron affinities

The method of Monte Carlo configuration interaction (MCCI) [1,2] is applied to the calculation of multipole moments. We look at the ground and excited state dipole moments in carbon monoxide. We then consider the dipole of NO, the quadrupole of the nitrogen molecule and of BH. An octupole of methane is also calculated. We consider experimental geometries and also stretched bonds. We show that these non-variational quantities may be found to relatively good accuracy when compared with FCI results, yet using only a small fraction of the full configuration interaction space. MCCI results in the aug-cc-pVDZ basis are seen to generally have reasonably good agreement with experiment. We also investigate the performance of MCCI when applied to ionisation energies and electron affinities of atoms in an aug-cc-pVQZ basis. We compare the MCCI results with full configuration-interaction quantum Monte Carlo [3,4] and `exact' non-relativistic results [3,4]. We show that MCCI could be a useful alternative for the calculation of atomic ionisation energies however electron affinities appear much more challenging for MCCI. Due to the small magnitude of the electron affinities their percentage errors can be high, but with regards to absolute errors MCCI performs similarly for ionisation energies and electron affinities.

preprint2012arXiv

Calculations of Potential Energy Surfaces Using Monte Carlo Configuration Interaction

We apply the method of Monte Carlo configuration interaction (MCCI) to calculate ground-state potential energy curves for a range of small molecules and compare the results with full configuration interaction. We show that the MCCI potential energy curve can be calculated to relatively good accuracy, as quantified using the non-parallelity error, using only a very small fraction of the FCI space. In most cases the potential curve is of better accuracy than its constituent single-point energies. We finally test the MCCI program on systems with basis sets beyond full configuration interaction: a lattice of fifty hydrogen atoms and ethylene. The results for ethylene agree fairly well with other computational work while for the lattice of fifty hydrogens we find that the fraction of the full configuration interaction space we were able to consider appears to be too small as, although some qualitative features are reproduced, the potential curve is less accurate.