BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David R. Bowler

David R. Bowler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.comp-phphysics.chem-phcond-mat.mes-hallBiological Physicscond-mat.softmath-phmath.MP

Catalog footprint

What is connected

21works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

21 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Generalized canonical purification for density matrix minimization

A Lagrangian formulation for the constrained search for the $N$-representable one-particle density matrix based on the McWeeny idempotency error minimization is proposed, which converges systematically to the ground state. A closed form of the canonical purification is derived for which no a posteriori adjustement on the trace of the density matrix is needed. The relationship with comparable methods are discussed, showing their possible generalization through the hole-particle duality. The appealing simplicity of this self-consistent recursion relation along with its low computational complexity could prove useful as an alternative to diagonalization in solving dense and sparse matrix eigenvalue problems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Raoult's law revisited: accurately predicting equilibrium relative humidity points for humidity control experiments

The equilibrium relative humidity values for a number of the most commonly used precipitants in biological macromolecule crystallisation have been measured using a new humidity control device. A simple argument in statistical mechanics demonstrates that the saturated vapour pressure of a solvent is proportional to its mole fraction in an ideal solution (Raoult's Law). The same argument can be extended to the case where solvent and solute molecules are of different size.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A linear scaling DFT study of the growth of a new \{105\} facet layer on a Ge hut cluster

Ge hut clusters, which appear during heteroepitaxy of Ge on Si(001), are prototypical examples of islands formed through strain. Experimentally, complete facets are observed to form rapidly, though the mechanism is unknown. We model the growth of new faces on Ge hut clusters, using linear scaling DFT. We build realistic small huts on substrates, and show that the growth of \{105\} facets proceeds from top-to-bottom, even for these small huts. The growth of the facet is driven by the reconstruction on the \{105\} facet, and nucleates at the boundaries of the facet. We thus resolve any ambiguities in kinetic models of hut cluster growth.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Linear Scaling Density Matrix Real Time TDDFT: Propagator Unitarity \& Matrix Truncation

Real time, density matrix based, time dependent density functional theory proceeds through the propagation of the density matrix, as opposed to the Kohn-Sham orbitals. It is possible to reduce the computational workload by imposing spatial cut-off radii on sparse matrices, and the propagation of the density matrix in this manner provides direct access to the optical response of very large systems, which would be otherwise impractical to obtain using the standard formulations of TDDFT. Following a brief summary of our implementation, along with several benchmark tests illustrating the validity of the method, we present an exploration of the factors affecting the accuracy of the approach. In particular we investigate the effect of basis set size and matrix truncation, the key approximation used in achieving linear scaling, on the propagator unitarity and optical spectra. Finally we illustrate that, with an appropriate density matrix truncation range applied, the computational load scales linearly with the system size and discuss the limitations of the approach.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimized multi-site local orbitals in the large-scale DFT program CONQUEST

We introduce numerical optimization of multi-site support functions in the linear-scaling DFT code CONQUEST. Multi-site support functions, which are linear combinations of pseudo-atomic orbitals on a target atom and those neighbours within a cutoff, have been recently proposed to reduce the number of support functions to the minimal basis while keeping the accuracy of a large basis [J. Chem. Theory Comput., 2014, 10, 4813]. The coefficients were determined by using the local filter diagonalization (LFD) method [Phys. Rev. B, 2009, 80, 205104]. We analyse the effect of numerical optimization of the coefficients produced by the LFD method. Tests on crystalline silicon, a benzene molecule and hydrated DNA systems show that the optimization improves the accuracy of the multi-site support functions with small cutoffs. It is also confirmed that the optimization guarantees the variational energy minimizations with multi-site support functions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

DFT Studies of Adsorption of Benzoic Acid on the Rutile (110) Surface: Modes and Patterns

Adsorption of benzoic acid on the (110) surface of rutile, both unreconstructed and (1x2)- reconstructed ones, at saturation coverage of one benzoate per two adjacent five-coordinated Ti5c sites has been studied using DFT simulations, as implemented through the Vienna Ab initio Simulations Package (VASP). In order to study the effects of hydrogen bonding and Van der Waals forces in influencing the relative stabilities of different adsorbate overlayer superstructures, these studies were performed through Local Density Approximation (LDA), Generalized Gradient Approximation(GGA) and DFT-D2. Through the calculations, it was found out that although the optimized structures corresponded with the proposed models for the experimental results, the relative energetic stabilities of different overlayer structures have shown some differences with the experimental results. In the GGA calculations, the overlayer structures involving the benzene rings aligned along the (001)-direction were shown to be more stable than those aligned perpendicular to it, regardless of whether the benzoates were arranged in (1x1) or (1x2) symmetries on the (1x1)-unreconstructed surface, or on the (1x2)-reconstructed surface. These are in contrast to the experimental studies using STM, whose observations revealed that configurations with the benzene rings aligned along the (110)-direction when adsorbed on the (1x2)-reconstructed surface.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

DSSC Anchoring Groups: A Surface Dependent Decision

Electrodes in dye sensitised solar cells (DSSCs) are typically nanocrystalline anatase TiO2 with a majority (101) surface exposed. Generally the sensitising dye employs a carboxylic anchoring moiety through which it adheres to the TiO2 surface. Recent interest in exploiting the properties of differing TiO2 electrode morphologies, such as rutile nanorods exposing the (110) surface and anatase electrodes with high percentages of the (001) surface exposed, begs the question of whether this anchoring strategy is best, irrespective of the majority surface exposed. Here we address this question by presenting density functional theory calculations contrasting the binding properties of two promising anchoring groups, phosphonic acid and boronic acid, to that of carboxylic acid. Anchor-electrode interactions are studied for the pro- totypical anatase (101) surface, along with the anatase (001) and rutile (110) surfaces. Finally the effect of using these alternative anchoring groups to bind a typical coumarin dye (NKX- 2311) to these TiO2 substrates is examined. Significant differences in the binding properties are found depending on both the anchor and surface, illustrating that the choice of anchor is necessarily dependent upon the surface exposed in the electrode. In particular the boronic acid is found to show the potential to be an excellent anchor choice for electrodes exposing the anatase (001) surface.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficient Calculations with Multisite Local Orbitals in the large-scale DFT Code CONQUEST

Multisite local orbitals, which are formed from linear combinations of pseudo-atomic orbitals from a target atom and its neighbor atoms, have been introduced in the large-scale density functional theory calculation code CONQUEST. Multisite local orbitals correspond to local molecular orbitals so that the number of required local orbitals can be minimal. The multisite support functions are determined by using the localized filter diagonalization (LFD) method [Phys. Rev. B 2009, 80, 205104]. Two new methods, the double cutoff method and the smoothing method, are introduced to the LFD method to improve efficiency and stability. The Hamiltonian and overlap matrices with multisite local orbitals are constructed by efficient sparse-matrix multiplications in CONQUEST. The investigation of the calculated energetic and geometrical properties and band structures of bulk Si, Al and DNA systems demonstrate the accuracy and the computational efficiency of the present method. The representability of both occupied and unoccupied band structures with the present method has been also confirmed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hydrogen adsorption and diffusion around Si(001)/Si(110) corners in nanostructures

While the diffusion of hydrogen on silicon surfaces has been relatively well characterised both experimentally and theoretically, the diffusion around corners between surfaces, as will be found on nanowires and nanostructures, has not been studied. Motivated by nanostructure fabrication by Patterned Atomic Layer Epitaxy (PALE), we present a density functional theory (DFT) study of the diffusion of hydrogen around the edge formed by the orthogonal (001) and (110) surfaces in silicon. We find that the barrier from (001) to (110) is approximately 0.3 eV lower than from (110) to (001), and that it is comparable to diffusion between rows on the clean surface, with no significant effect on the hydrogen patterns at the growth temperatures used.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Intrinsic Oxygen Vacancy and Extrinsic Aluminium Dopant Interplay: A Route to the Restoration of Defective TiO$_2$

Density functional theory (DFT) and DFT corrected for on-site Coulomb interactions (DFT+U) calculations are presented on Aluminium doping in bulk TiO$_2$ and the anatase (101) surface. Particular attention is paid to the mobility of oxygen vacancies throughout the doped TiO$_2$ lattice, as a means by which charge compensation of trivalent dopants can occur. The effect that Al doping of TiO$_2$ electrodes has in dye sensitised solar cells is explained as a result of this mobility and charge compensation. Substitutional defects in which one Al3+ replaces one Ti4+ are found to introduce valence band holes, while intrinsic oxygen vacancies are found to introduce states in the band-gap. Coupling two of these substitutional defects with an oxygen vacancy results in exothermic defect formation which maintain charge neutrality. Nudged elastic band calculations have been performed to investigate the formation of these clustered defects in the (101) surface by oxygen vacancy diffusion, with the resulting potential energy surface suggesting energetic gains with small diffusion barriers. Efficiency in- creases observed in dye sensitised solar cells as a result of aluminium doping of TiO$_2$ electrodes are investigated by adsorbing the tetrahydroquinoline C2-1 chromophore on the defective surfaces. Adsorption on the clustered extrinsic Al3+ and intrinsic oxygen vacancy defects are found to behave as if adsorbed on a clean surface, with vacancy states not present, while adsorption on the oxygen vacancy results in a down shift of the dye localised states within the band-gap and defect states being present below the conduction band edge. Aluminium doping therefore acts as a benign dopant for 'cleaning' TiO$_2$ through oxygen vacancy diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stable and Efficient Linear Scaling First-Principles Molecular Dynamics for 10,000+ atoms

The recent progress of linear-scaling or O(N) methods in the density functional theory (DFT) is remarkable. We expect that first-principles molecular dynamics (FPMD) simulations based on DFT can now treat more realistic and complex systems using the O(N) technique. However, very few examples of O(N) FPMD simulations exist so far and the information for the accuracy or reliability of the simulations is very limited. In this paper, we show that efficient and robust O(N) FPMD simulations are now possible by the combination of the extended Lagrangian Born-Oppenheimer molecular dynamics method, which was recently proposed by Niklasson et al (Phys. Rev. Lett. 100, 123004 (2008)), and the density matrix method as an O(N) technique. Using our linear-scaling DFT code Conquest, we investigate the reliable calculation conditions for the accurate O(N) FPMD and demonstrate that we are now able to do actual and reliable self-consistent FPMD simulation of a very large system containing 32,768 atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Investigating individual arsenic dopant atoms in silicon using low-temperature scanning tunnelling microscopy

We study sub-surface arsenic dopants in a hydrogen terminated Si(001) sample at 77 K, using scanning tunnelling microscopy and spectroscopy. We observe a number of different dopant related features that fall into two classes, which we call As1 and As2. When imaged in occupied states the As1 features appear as anisotropic protrusions superimposed on the silicon surface topography, and have maximum intensities lying along particular crystallographic orientations. In empty-state images the features all exhibit long-range circular protrusions. The images are consistent with buried dopants that are in the electrically neutral (D0) charge state when imaged in filled states, but become positively charged (D+) through electrostatic ionisation when imaged under empty state conditions, similar to previous observations of acceptors in GaAs. Density functional theory (DFT) calculations predict that As dopants in the third layer of the sample induce two states lying just below the conduction band edge, which hybridize with the surface structure creating features with the surface symmetry consistent with our STM images. The As2 features have the surprising characteristic of appearing as a protrusion in filled state images and an isotropic depression in empty state images, suggesting they are negatively charged at all biases. We discuss the possible origins of this feature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Understanding the role of ions and water molecules in the NaCl dissolution process

The dissolution of NaCl in water is one of the most common everyday processes, yet it remains poorly understood at the molecular level. Here we report the results of an extensive density functional theory study in which the initial stages of NaCl dissolution have been examined at low water coverages. Our specific approach is to study how the energetic cost of moving an ion or a pair of ions to a less coordinated site at the surface of various NaCl crystals varies with the number of water molecules adsorbed on the surface. This "microsolvation" approach allows us to study the dependence of the defect energies on the number of water molecules in the cluster and thus to establish when and where dissolution becomes favorable. Moreover, this approach allows us to understand the roles of the individual ions and water molecules in the dissolution process. Consistent with previous work we identify a clear preference for dissolution of Cl ions over Na ions. However, the detailed information obtained here leads to the conclusion that the process is governed by the higher affinity of the water molecules to Na ions than to Cl ions. The Cl ions are released first as this exposes more Na ions at the surface creating favorable adsorption sites for water. We discuss how this mechanism is likely to be effective for other alkali halides.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Linear-scaling implementation of exact exchange using localized numerical orbitals and contraction reduction integrals

We present enhancements to the computational efficiency of exact exchange calculations using the density matrix and local support functions. We introduce a numerical method which avoids the explicit calculation the four-center two-electron repulsion integrals and reduces the prefactor scaling by a factor N, where N is the number of atoms within the range of the exact exchange Hamiltonian. This approach is based on a contraction-reduction scheme, and takes advantage of the discretization space which enables the direct summation over the support functions in a localized space. Using the sparsity property of the density matrix, the scaling of the prefactor can be further reduced to reach asymptotically O(N).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Van der Waals density functionals applied to solids

The van der Waals density functional (vdW-DF) of Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)] is a promising approach for including dispersion in approximate density functional theory exchange-correlation functionals. Indeed, an improved description of systems held by dispersion forces has been demonstrated in the literature. However, despite many applications, standard general tests on a broad range of materials are lacking. Here we calculate the lattice constants, bulk moduli, and atomization energies for a range of solids using the original vdW-DF and several of its offspring. We find that the original vdW-DF overestimates lattice constants in a similar manner to how it overestimates binding distances for gas phase dimers. However, some of the modified vdW functionals lead to average errors which are similar to those of PBE or better. Likewise, atomization energies that are slightly better than from PBE are obtained from the modified vdW-DFs. Although the tests reported here are for "hard" solids, not normally materials for which dispersion forces are thought to be important, we find a systematic improvement in cohesive properties for the alkali metals and alkali halides when non-local correlations are accounted for.

0 citations0 reviewsNo code linkedOpen access