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 D. O'Regan

David D. O'Regan 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-scicond-mat.str-elphysics.chem-phphysics.comp-phBiological Physicsquant-phBiomoleculescond-mat.mes-hallcond-mat.othercond-mat.stat-mechphysics.app-phphysics.atm-clus

Catalog footprint

What is connected

15works
12topics
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

15 published item(s)

preprint2020arXiv

Equipartition of Energy Defines the Size-Thickness Relationship in Liquid-Exfoliated Nanosheets

Liquid phase exfoliation is a commonly used method to produce 2D nanosheets from a range of layered crystals. However, such nanosheets display broad size and thickness distributions and correlations between area and thickness, issues that limit nanosheet application potential. To understand the factors controlling the exfoliation process, we have liquid-exfoliated 11 different layered materials, size-selecting each into fractions before using AFM to measure the nanosheet length, width, and thickness distributions for each fraction. The resultant data show a clear power-law scaling of nanosheet area with thickness for each material. We have developed a simple nonequilibrium thermodynamics-based model predicting that the power-law prefactor is proportional to both the ratios of in-plane-tearing/out-of-plane-peeling energies and in-plane/out-of-plane moduli. By comparing the experimental data with the modulus ratio calculated from first-principles, we find close agreement between experiment and theory. This supports our hypothesis that energy equipartition holds between nanosheet tearing and peeling during sonication-assisted exfoliation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

First-principles Hubbard U and Hund's J corrected approximate density-functional theory predicts an accurate fundamental gap in rutile and anatase TiO2

Titanium dioxide (TiO$_2$) presents a long-standing challenge for approximate Kohn-Sham density-functional theory (KS-DFT), as well as to its Hubbard-corrected extension, DFT+$U$. We find that a previously proposed extension of first-principles DFT+$U$ to incorporate a Hund's $J$ correction, termed DFT+$U$+$J$, in combination with parameters calculated using a recently proposed linear-response theory, predicts fundamental band-gaps accurate to well within the experimental uncertainty in rutile and anatase TiO$_2$. Our approach builds upon established findings that Hubbard correction to both titanium $3d$ and oxygen $2p$ subspaces in TiO$_2$, symbolically giving DFT+$U^{d,p}$, is necessary to achieve acceptable band-gaps using DFT+$U$. This requirement remains when the first-principles Hund's $J$ is included. We also find that the calculated gap depends on the correlated subspace definition even when using subspace-specific first-principles $U$ and $J$ parameters. Using the simplest reasonable correlated subspace definition and underlying functional, the local density approximation, we show that high accuracy results from using a relatively uncomplicated form of the DFT+$U$+$J$ functional. For closed-shell systems such as TiO$_2$, we describe how various DFT+$U$+$J$ functionals reduce to DFT+$U$ with suitably modified parameters, so that reliable band gaps can be calculated for rutile and anatase with no modifications to a conventional DFT+$U$ code.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Reproducibility in $G_0W_0$ Calculations for Solids

Ab initio many-body perturbation theory within the $GW$ approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn-Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange correlation functional ``starting point'', are used to construct $G$ and $W$, and then perturbatively corrected by the resultant $GW$ self-energy. In practice, there are several ways to construct the $GW$ self-energy, and these can lead to variations in predicted quasiparticle energies. For example, for ZnO and TiO$_2$, reported $GW$ fundamental gaps can vary by more than 1 eV. In this work, we address the convergence and key approximations in contemporary $G_0W_0$ calculations, including frequency-integration schemes and the treatment of the Coulomb divergence in the exact-exchange term. We study several systems,and compare three different $GW$ codes: BerkeleyGW, Abinit and Yambo. We demonstrate, for the first time, that the same quasiparticle energies for systems in the condensed phase can be obtained with different codes, and we provide a comprehensive assessment of implementations of the $GW$ approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

XDFT: an efficient first-principles method for neutral excitations in molecules

State-of-the-art methods for calculating neutral excitation energies are typically demanding and limited to single electron-hole pairs and their composite plasmons. Here we introduce excitonic density-functional theory (XDFT) a computationally light, generally applicable, first-principles technique for calculating neutral excitations based on generalized constrained DFT. In order to simulate an M-particle excited state of an N-electron system, XDFT automatically optimizes a constraining potential to confine N-M electrons within the ground-state Kohn-Sham valence subspace. We demonstrate the efficacy of XDFT by calculating the lowest single-particle singlet and triplet excitation energies of the well-known Thiel molecular test set, with results which are in excellent agreement with time-dependent DFT. Furthermore, going beyond the capability of adiabatic time-dependent DFT, we show that XDFT can successfully capture double excitations. Overall our method makes optical gaps, excition bindings and oscillator strengths readily accessible at a computational cost comparable to that of standard DFT. As such, XDFT appears as an ideal candidate to work within high-throughput discovery frameworks and within linear-scaling methods for large systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inapplicability of exact constraints and a minimal two-parameter generalization to the DFT+$U$ based correction of self-interaction error

In approximate density functional theory (DFT), the self-interaction error is an electron delocalization anomaly associated with underestimated insulating gaps. It exhibits a predominantly quadratic energy-density curve that is amenable to correction using efficient, constraint-resembling methods such as DFT + Hubbard $U$ (DFT+$U$). Constrained DFT (cDFT) enforces conditions on DFT exactly, by means of self-consistently optimized Lagrange multipliers, and while its use to automate error corrections is a compelling possibility, we show that it is limited by a fundamental incompatibility with constraints beyond linear order. We circumvent this problem by utilizing separate linear and quadratic correction terms, which may be interpreted either as distinct constraints, each with its own Hubbard $U$ type Lagrange multiplier, or as the components of a generalized DFT+$U$ functional. The latter approach prevails in our tests on a model one-electron system, $H_2^+$, in that it readily recovers the exact total-energy while symmetry-preserving pure constraints fail to do so. The generalized DFT+$U$ functional moreover enables the simultaneous correction of the total-energy and ionization potential or the correction of either together with the enforcement of Koopmans condition. For the latter case, we outline a practical, approximate scheme by which the required pair of Hubbard parameters, denoted as U1 and U2, may be calculated from first-principles.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimization of constrained density functional theory

Constrained density functional theory (cDFT) is a versatile electronic structure method that enables ground-state calculations to be performed subject to physical constraints. It thereby broadens their applicability and utility. Automated Lagrange multiplier optimisation is necessary for multiple constraints to be applied efficiently in cDFT, for it to be used in tandem with geometry optimization, or with molecular dynamics. In order to facilitate this, we comprehensively develop the connection between cDFT energy derivatives and response functions, providing a rigorous assessment of the uniqueness and character of cDFT stationary points while accounting for electronic interactions and screening. In particular, we provide a new, non-perturbative proof that stable stationary points of linear density constraints occur only at energy maxima with respect to their Lagrange multipliers. We show that multiple solutions, hysteresis, and energy discontinuities may occur in cDFT. Expressions are derived, in terms of convenient by-products of cDFT optimization, for quantities such as the dielectric function and a condition number quantifying ill-definition in multi-constraint cDFT.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Supercell convergence of charge-transfer energies in pentacene molecular crystals from constrained DFT

Singlet fission (SF) is a multi-exciton generation process that could be harnessed to improve the efficiency of photovoltaic devices. Experimentally, systems derived from the pentacene molecule have been shown to exhibit ultrafast SF with high yields. Charge-transfer (CT) configurations are likely to play an important role as intermediates in the SF process in these systems. In molecular crystals, electrostatic screening effects and band formation can be significant in lowering the energy of CT states, enhancing their potential to effectively participate in SF. In order to simulate these, it desirable to adopt a computational approach which is acceptably accurate, relatively inexpensive, which and scales well to larger systems, thus enabling the study of screening effects. We propose a novel, electrostatically-corrected constrained Density Functional Theory (cDFT) approach as a low-cost solution to the calculation of CT energies in molecular crystals such as pentacene. Here we consider an implementation in the context of the ONETEP linear-scaling DFT code, but our electrostatic correction method is in principle applicable in combination with any constrained DFT implementation, also outside the linear-scaling framework. Our newly developed method allows us to estimate CT energies in the infinite crystal limit, and with these to validate the accuracy of the cluster approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Renormalization of myoglobin-ligand binding energetics by quantum many-body effects

We carry out a first-principles atomistic study of the electronic mechanisms of ligand binding and discrimination in the myoglobin protein. Electronic correlation effects are taken into account using one of the most advanced methods currently available, namely a linear-scaling density functional theory (DFT) approach wherein the treatment of localized iron 3d electrons is further refined using dynamical mean-field theory (DMFT). This combination of methods explicitly accounts for dynamical and multi-reference quantum physics, such as valence and spin fluctuations, of the 3d electrons, whilst treating a significant proportion of the protein (more than 1000 atoms) with density functional theory. The computed electronic structure of the myoglobin complexes and the nature of the Fe-O2 bonding are validated against experimental spectroscopic observables. We elucidate and solve a long standing problem related to the quantum-mechanical description of the respiration process, namely that DFT calculations predict a strong imbalance between O2 and CO binding, favoring the latter to an unphysically large extent. We show that the explicit inclusion of many body-effects induced by the Hund's coupling mechanism results in the correct prediction of similar binding energies for oxy- and carbonmonoxymyoglobin.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Importance of many body effects in the kernel of hemoglobin for ligand binding

We propose a mechanism for binding of diatomic ligands to heme based on a dynamical orbital selection process. This scenario may be described as bonding determined by local valence fluctuations. We support this model using linear-scaling first-principles calculations, in combination with dynamical mean-field theory, applied to heme, the kernel of the hemoglobin metalloprotein central to human respiration. We find that variations in Hund's exchange coupling induce a reduction of the iron 3d density, with a concomitant increase of valence fluctuations. We discuss the comparison between our computed optical absorption spectra and experimental data, our picture accounting for the observation of optical transitions in the infrared regime, and how the Hund's coupling reduces, by a factor of five, the strong imbalance in the binding energies of heme with CO and O_2 ligands.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ligand Discrimination in Myoglobin from Linear-Scaling DFT+U

Myoglobin modulates the binding of diatomic molecules to its heme group via hydrogen-bonding and steric interactions with neighboring residues, and is an important benchmark for computational studies of biomolecules. We have performed calculations on the heme binding site and a significant proportion of the protein environment (more than 1000 atoms) using linear-scaling density functional theory and the DFT+U method to correct for self-interaction errors associated with localized 3d states. We confirm both the hydrogen-bonding nature of the discrimination effect (3.6 kcal/mol) and assumptions that the relative strain energy stored in the protein is low (less than 1 kcal/mol). Our calculations significantly widen the scope for tackling problems in drug design and enzymology, especially in cases where electron localization, allostery or long-ranged polarization influence ligand binding and reaction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Generalized Wannier functions: a comparison of molecular electric dipole polarizabilities

Localized Wannier functions provide an efficient and intuitive means by which to compute dielectric properties from first principles. They are most commonly constructed in a post-processing step, following total-energy minimization. Nonorthogonal generalized Wannier functions (NGWFs) [Skylaris et al., Phys. Rev. B 66, 035119 11 (2002); Skylaris et al., J. Chem. Phys. 122, 084119 (2005)] may also be optimized in situ, in the process of solving for the ground-state density. We explore the relationship between NGWFs and orthonormal, maximally localized Wannier functions (MLWFs) [Marzari and Vanderbilt, Phys. Rev. B 56, 12847 (1997); Souza, Marzari, and Vanderbilt, ibid. 65, 035109 (2001)], demonstrating that NGWFs may be used to compute electric dipole polarizabilities efficiently, with no necessity for post-processing optimization, and with an accuracy comparable to MLWFs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Linear-scaling DFT+U with full local orbital optimization

We present an approach to the DFT+U method (Density Functional Theory + Hubbard model) within which the computational effort for calculation of ground state energies and forces scales linearly with system size. We employ a formulation of the Hubbard model using nonorthogonal projector functions to define the localized subspaces, and apply it to a local-orbital DFT method including in situ orbital optimization. The resulting approach thus combines linear-scaling and systematic variational convergence. We demonstrate the scaling of the method by applying it to nickel oxide nano-clusters with sizes exceeding 7,000 atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Vanadium dioxide : A Peierls-Mott insulator stable against disorder

Vanadium dioxide undergoes a first order metal-insulator transition at 340 K. In this work, we develop and carry out state of the art linear scaling DFT calculations refined with non-local dynamical mean-field theory. We identify a complex mechanism, a Peierls-assisted orbital selection Mott instability, which is responsible for the insulating M$_1$ phase, and furthermore survives a moderate degree of disorder.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Subspace representations in ab initio methods for strongly correlated systems

We present a generalized definition of subspace occupancy matrices in ab initio methods for strongly correlated materials, such as DFT+U and DFT+DMFT, which is appropriate to the case of nonorthogonal projector functions. By enforcing the tensorial consistency of all matrix operations, we are led to a subspace projection operator for which the occupancy matrix is tensorial and accumulates only contributions which are local to the correlated subspace at hand. For DFT+U in particular, the resulting contributions to the potential and ionic forces are automatically Hermitian, without resort to symmetrization, and localized to their corresponding correlated subspace. The tensorial invariance of the occupancies, energies and ionic forces is preserved. We illustrate the effect of this formalism in a DFT+U study using self-consistently determined projectors.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Projector self-consistent DFT+U using non-orthogonal generalized Wannier functions

We present a formulation of the density-functional theory + Hubbard model (DFT+U) method that is self-consistent over the choice of Hubbard projectors used to define the correlated subspaces. In order to overcome the arbitrariness in this choice, we propose the use of non-orthogonal generalized Wannier functions (NGWFs) as projectors for the DFT+U correction. We iteratively refine these NGWF projectors and, hence, the DFT+U functional, such that the correlated subspaces are fully self-consistent with the DFT+U ground-state. We discuss the convergence characteristics of this algorithm and compare ground-state properties thus computed with those calculated using hydrogenic projectors. Our approach is implemented within, but not restricted to, a linear-scaling DFT framework, opening the path to DFT+U calculations on systems of unprecedented size.

0 citations0 reviewsNo code linkedOpen access