Researcher profile

David R. Bowler

David R. Bowler contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
4topics
4close 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)

preprint2022arXiv

Large-scale DFT methods for calculations of materials with complex structures

Large-scale density functional theory (DFT) calculations provide a powerful tool to investigate the atomic and electronic structure of materials with complex structures. This article reviews a large-scale DFT calculation method, the multi-site support function (MSSF) method, in the CONQUEST code. MSSFs are linear combinations of the basis functions which belong to a group of atoms in a local region. The method can reduce the computational time while preserving accuracy. The accuracy of MSSFs has been assessed for bulk Si, Al, Fe and NiO and hydrated DNA, which demonstrate the applicability of the MSSFs for varied materials. The applications of MSSFs on large systems with several thousand atoms, which have complex interfaces and non-periodic structures, indicate that the MSSF method is promising for precise investigations of materials with complex structures.

preprint2020arXiv

Blue moon ensemble simulation of aquation free energy profiles applied to mono and bifunctional platinum anticancer drugs

Aquation free energy profiles of neutral cisplatin and cationic monofunctional derivatives, including triaminochloroplatinum(II) and cis-diammine(pyridine)chloroplatinum(II), were computed using state of the art thermodynamic integration, for which temperature and solvent were accounted for explicitly using density functional theory based canonical molecular dynamics (DFT-MD). For all the systems the "inverse-hydration" where the metal center acts as an acceptor of hydrogen bond has been observed. This has motivated to consider the inversely bonded solvent molecule in the definition of the reaction coordinate required to initiate the constrained DFT-MD trajectories. We found that there exists little difference in free enthalpies of activations, such that these platinum-based anticancer drugs are likely to behave the same way in aqueous media. Detailed analysis of the microsolvation structure of the square-planar complexes, along with the key steps of the aquation mechanism are discussed.

preprint2020arXiv

Insight into the charge density wave gap from contrast inversion in topographic STM images

Charge density waves (CDWs) are understood in great details in one dimension, but they remain largely enigmatic in two dimensional systems. In particular, numerous aspects of the associated energy gap and the formation mechanism are not fully understood. Two long standing riddles are the amplitude and position of the CDW gap with respect to the Fermi level ($E_F$) and the frequent absence of CDW contrast inversion (CI) between opposite bias scanning tunneling microscopy (STM) images. Here, we find compelling evidence that these two issues are intimately related. Combining density functional theory and STM to analyse the CDW pattern and modulation amplitude in 1$T$-TiSe$_2$, we find that CI takes place at an unexpected negative sample bias because the CDW gap opens away from $E_F$, deep inside the valence band. This bias becomes increasingly negative as the CDW gap shifts to higher binding energy with electron doping. This study shows the importance of CI in STM images to identify periodic modulations with a CDW and to gain valuable insight into the CDW gap, whose measurement is notoriously controversial.

preprint2020arXiv

Large scale and linear scaling DFT with the CONQUEST code

We survey the underlying theory behind the large-scale and linear scaling DFT code, Conquest, which shows excellent parallel scaling and can be applied to thousands of atoms with exact solutions, and millions of atoms with linear scaling. We give details of the representation of the density matrix and the approach to finding the electronic ground state, and discuss the implementation of molecular dynamics with linear scaling. We give an overview of the performance of the code, focussing in particular on the parallel scaling, and provide examples of recent developments and applications.

preprint2020arXiv

Polar morphologies from first principles: PbTiO$_3$ films on SrTiO$_3$ substrates and the $p(2 \times Λ)$ surface reconstruction

Low dimensional structures comprised of ferroelectric (FE) PbTiO$_3$ (PTO) and quantum paraelectric SrTiO$_3$ (STO) are hosts to complex polarization textures such as polar waves, flux-closure domains and polar skyrmion phases. Density functional theory (DFT) simulations can provide insight into this order, but, are limited by the computational effort needed to simulate the thousands of required atoms. To relieve this issue, we use the novel multi-site support function (MSSF) method within DFT to reduce the solution time for the electronic groundstate whilst preserving high accuracy. Using MSSFs, we simulate thin PTO films on STO substrates with system sizes $>2000$ atoms. In the ultrathin limit, the polar wave texture with cylindrical chiral bubbles emerges as an intermediate phase between full flux closure domains and in-plane polarization. This is driven by an internal bias field born of the compositionally broken inversion symmetry in the [001] direction. Since the exact nature of this bias field depends sensitively on the film boundary conditions, this informs a new principle of design for manipulating chiral order on the nanoscale through the careful choice of substrate, surface termination or use of overlayers. Antiferrodistortive (AFD) order locally interacts with these polar textures giving rise to strong FE/AFD coupling at the PbO terminated surface driving a $p(2 \times Λ)$ surface reconstruction. This offers another pathway for the local control of ferroelectricity.

preprint2020arXiv

The pseudoatomic orbital basis: electronic accuracy and soft-mode distortions in ABO$_3$ perovskites

The perovskite oxides are known to be susceptible to structural distortions over a long wavelength when compared to their parent cubic structures. From an ab initio simulation perspective, this requires accurate calculations including many thousands of atoms; a task well beyond the remit of traditional plane wave-based density functional theory (DFT). We suggest that this void can be filled using the methodology implemented in the large-scale DFT code, CONQUEST, using a local pseudoatomic orbital (PAO) basis. Whilst this basis has been tested before for some structural and energetic properties, none have treated the most fundamental quantity to the theory, the charge density $n(\mathbf{r})$ itself. An accurate description of $n(\mathbf{r})$ is vital to the perovskite oxides due to the crucial role played by short-range restoring forces (characterised by bond covalency) and long range coulomb forces as suggested by the soft-mode theory of Cochran and Anderson. We find that modestly sized basis sets of PAOs can reproduce the plane-wave charge density to a total integrated error of better than 0.5% and provide Bader partitioned ionic charges, volumes and average charge densities to similar degree of accuracy. Further, the multi-mode antiferroelectric distortion of PbZrO$_3$ and its associated energetics are reproduced by better than 99% when compared to plane-waves. This work suggests that electronic structure calculations using efficient and compact basis sets of pseudoatomic orbitals can achieve the same accuracy as high cutoff energy plane-wave calculations. When paired with the CONQUEST code, calculations with high electronic and structural accuracy can now be performed on many thousands of atoms, even on systems as delicate as the perovskite oxides.