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

Tim Gould

Tim Gould 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.chem-phphysics.comp-phcond-mat.otherphysics.atom-phcond-mat.mes-hallquant-phcond-mat.stat-mechcond-mat.str-elphysics.class-ph

Catalog footprint

What is connected

26works
10topics
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

26 published item(s)

preprint2022arXiv

Effect of exchange-correlation and core-electron approximations on the calculated superconducting transition temperature of palladium hydride

Realistic prediction of the superconducting transition temperature (Tc) for PdH is a long-standing challenge, because it depends on robust calculations of the electron and phonon band structures to obtain the electron-phonon scattering matrix element. To date, first-principles calculations employing density functional theory (DFT) have been based on selected exchange-correlation and core-electron approximations. Incorporating anharmonicity produced a more realistic value of Tc that nevertheless still disagreed strongly with experiment unless adjustable parameters were introduced. Here we consider how the value of Tc predicted using DFT in the harmonic approximation depends on the DFT scheme employed. The rationale for this work is that unless the calculation of Tc within the harmonic approximation is robust, albeit incorrect, there is not a solid foundation for incorporating anharmonicity meaningfully. Six combinations of exchange-correlation approximation (LDA, PBE, PBEsol) and core-electron approximation (PAW, USPP) were tested. Following a carefully systematic methodology, the calculated Tc was found to vary by a factor exceeding two across the tested DFT schemes. This suggests strongly that "standard" DFT, even including anharmonicity, is not reliable for PdH, implying that a higher-rung method will be needed to calculate a realistic lattice constant and phonon band structure, and so predict Tc accurately.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Combining density functional theories to correctly describe the energy, lattice structure and electronic density of functional oxide perovskites

Functional oxide perovskites are the pillar of cutting-edge technological applications. Density functional theory (DFT) simulations are the theoretical methods of choice to understand and design perovskite materials. However, tests on the reliability of DFT to describe fundamental properties of oxide perovskites are scarce and mostly ill-defined due to a lack of rigorous theoretical benchmarks for solids. Here, we present a quantum Monte Carlo benchmark study of DFT on the archetypal perovskite BaTiO$_{3}$ (BTO). It shows that no DFT approximation can simultaneously reproduce the energy, structure, and electronic density of BTO. Traditional protocols to select DFT approximations are empirical and fail to detect this shortcoming. An approach combining two different non-empirical DFT schemes, "SCAN" and "HSE06", is able to holistically describe BTO with accuracy. Combined DFT approaches should thus be considered as a promising alternative to standard methods for simulating oxide perovskites.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Density driven correlations in ensemble density functional theory: insights from simple excitations in atoms

Ensemble density functional theory extends the usual Kohn-Sham machinery to quantum state ensembles involving ground- and excited states. Recent work by the authors [Phys. Rev. Lett. 119, 243001 (2017); 123, 016401 (2019)] has shown that both the Hartree-exchange and correlation energies can attain unusual features in ensembles. Density-driven(DD) correlations -- which account for the fact that pure-state densities in Kohn-Sham ensembles do not necessarily reproduce those of interacting pure states -- are one such feature. Here we study atoms (specifically $S$--$P$ and $S$--$S$ transitions) and show that the magnitude and behaviour of DD correlations can vary greatly with the variation of the orbital angular momentum of the involved states. Such estimations are obtained through an approximation for DD correlations built from relevant exact conditions Kohn-Sham inversion, and plausible assumptions for weakly correlated systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Strong correlation and charge localisation in Kohn-Sham theories with fractional orbital occupations

We study static correlation and delocalisation errors and show that even methods with good energies can yield significant delocalization errors that affect the density, leading to large errors in predicting {\em e.g.} dipole moments. We illustrate this point by comparing existing state-of-art approaches with an accurate exchange correlation functional based on a generalised valence-bond ansatz, in which orbitals and fractional occupations are treated as variational parameters via an optimized effective potential (OEP). We show that the OEP exhibits step and peak features which, similar to the exact Kohn-Sham (KS) potential of DFT, are crucial to prevent charge delocalization. We further show that the step is missing in common approximations within reduced density matrix functional theory resulting in delocalization errors comparable to those found in DFT approximations. Finally, we explain the delocalization error as coming from an artificial mixing of the ground state with a charge-transfer excited state which is avoided if occupation numbers exhibit discontinuities.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermodynamics from first principles: correlations and nonextensivity

The standard formulation of thermostatistics, being based on the Boltzmann-Gibbs distribution and logarithmic Shannon entropy, describes idealized uncorrelated systems with extensive energies and short-range interactions. In this letter, we use the fundamental principles of ergodicity (via Liouville's theorem), the self-similarity of correlations, and the existence of the thermodynamic limit to derive generalized forms of the equilibrium distribution for long-range-interacting systems. Significantly, our formalism provides a justification for the well-studied nonextensive thermostatistics characterized by the Tsallis distribution, which it includes as a special case. We also give the complementary maximum entropy derivation of the same distributions by constrained maximization of the Boltzmann-Gibbs-Shannon entropy. The consistency between the ergodic and maximum entropy approaches clarifies the use of the latter in the study of correlations and nonextensive thermodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A simple self-interaction correction to RPA-like correlation energies

The random phase approximation (RPA) is exact for the exchange energy of a many-electron ground state, but RPA makes the correlation energy too negative by about 0.5 eV/electron. That large short-range error, which tends to cancel out of iso-electronic energy differences, is largely corrected by an exchange-correlation kernel, or (as in RPA+) by an additive local or semilocal correction. RPA+ is by construction exact for the homogeneous electron gas, and it is also accurate for the jellium surface. RPA+ often gives realistic total energies for atoms or solids in which spin-polarization corrections are absent or small. RPA and RPA+ also yield realistic singlet binding energy curves for H2 and N2, and thus RPA+ yields correct total energies even for spin-unpolarized atoms with fractional spins and strong correlation, as in stretched H2 or N2. However, RPA and RPA+ can be very wrong for spin-polarized one-electron systems (especially for stretched H2+), and also for the spin-polarization energies of atoms. The spin-polarization energy is often a small part of the total energy of an atom, but important for ionization energies, electron affinities, and the atomization energies of molecules. Here we propose a computationally efficient generalized RPA+ (gRPA+) that changes RPA+ only for spin-polarized systems by making gRPA+ exact for all one-electron densities, in the same simple semilocal way that the correlation energy densities of many meta-generalized gradent approximations are made self-correlation free. By construction, gRPA+ does not degrade the exact RPA+ description of jellium. gRPA+ is found to greatly improve upon RPA and RPA+ for the ionization energies and electron affinities of light atoms. Many versions of RPA with an approximate exchange-correlation kernel fail to be exact for all one-electron densities, and they can also be self-interaction corrected in this way.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Faraday-cage screening reveals intrinsic aspects of the van der Waals attraction

General properties of the recently observed screening of the van der Waals (vdW) attraction between a silica substrate and silica tip by insertion of graphene are predicted using basic theory and first-principles calculations. Results are then focused on possible practical applications, as well as an understanding of the nature of vdW attraction, considering recent discoveries showing it competing against covalent and ionic bonding. The traditional view of the vdW attraction as arising from pairwise-additive London dispersion forces is considered using Grimme's "D3" method, comparing results to those from Tkatchenko's more general many-body dispersion (MBD) approach, all interpreted in terms of Dobson's general dispersion framework. Encompassing the experimental results, MBD screening of the vdW force between two silica bilayers is shown to scale up to medium separations as 1.25 de/d, where d is the bilayer separation and de its equilibrium value, depicting antiscreening approaching and inside de. Means of unifying this correlation effect with those included in modern density functionals are urgently required.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

$C_6$ coefficients and dipole polarizabilities for all atoms and many ions in rows 1-6 of the periodic table

Using time-dependent density functional theory (tdDFT) with exchange kernels we calculate and test imaginary frequency-dependent dipole polarizabilities for all atoms and many ns in rows 1-6 of the periodic table. These are then integrated over frequency to produce $C_6$ coefficients. Results are presented under different models: straight tdDFT calculations using two different kernels, "benchmark" tdDFT calculations corrected by more accurate quantum chemical and experimental data, and "benchmark" tdDFT with frozen orbital anions. Parametrisations are presented for 411+ atoms and ions, allowing results to be easily used by other researchers. A curious relationship, $C_{6,XY}\propto [α_X(0)α_Y(0)]^{0.73}$ is found between $C_6$ coefficients and static polarizabilities $α(0)$. The relationship $C_{6,XY}=2C_{6,X}C_{6,Y}/[α_X/α_YC_{6,Y}+α_Y/α_XC_{6,X}]$ is tested and found to work well ($<5$\% errors) in about 80\% of cases, but can break down badly ($>30$\% errors) in a small fraction of cases.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Graphene as a p-type metal for ultimate miniaturization

We report macroscopic sheets of highly conductive bilayer graphene with exceptionally high hole concentrations of ~ $10^{15}$ $cm^{-2}$ and unprecedented sheet resistances of 20-25 Ω per square over macroscopic scales, and obtained in-situ over a thin cushion of molecular oxygen on a silicon substrate. The electric and electronic properties of this specific configuration remain stable upon thermal anneals and months of exposure to air. We further report a complementary ab-initio study, predicting an enhancement of graphene adhesion energy of up to a factor 20, also supported by experimental fracture tests. Our results show that the remarkable properties of graphene can be realized in a reliable fashion using a high-throughput process. In addition to providing exceptional material properties, the growth process we employed is scalable to large areas so that the outstanding conduction properties of graphene can be harnessed in devices fabricated via conventional semiconductor manufacturing processes. We anticipate that the approach will provide the necessary scalability and reliability for future developments in the graphene nanoscience and technology fields, especially in areas where further miniaturization is hampered by size effects and electrical reliability of classical conductors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How polarizabilities and $C_6$ coefficients actually vary with atomic volume

In this work we investigate how atomic $C_6$ coefficients and static dipole polarizabilities $α$ scale with effective volume. We show, using confined atoms covering rows 1-5 of the periodic table, that $C_6/C_6^R\approx (V/V^R)^{p_Z}$ and $α/α^R\approx (V/V^R)^{p'_Z}$ (for volume $V=\int dr \frac{4π}{3}r^3 n(r)$) where $C_6^R$, $α^R$ and $V^R$ are the reference values and effective volume of the free atom. The scaling exponents $p_Z$ and $p'_Z$ vary substantially as a function of element number $Z=N$, in contrast to the standard "rule of thumb" that $p_Z=2$ and $p'_Z=1$. Remarkably, We find that the polarizability and $C_6$ exponents $p'$ and $p$ are related by $p'\approx p-0.615$ rather than the expected $p'\approx p/2$. Results are largely independent of the form of the confining potential (harmonic, cubic and quartic potentials are considered) and kernel approximation, justifying this analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Layer Response Theory: Energetics of layered materials from semi-analytic high-level theory

We present a readily computable semi-analytic Layer Response Theory (LRT) for analysis of cohesive energetics involving two-dimensional layers such as BN or graphene. The theory approximates the Random Phase Approximation (RPA) correlation energy. Its RPA character ensures that the energy has the correct van der Waals asymptotics for well-separated layers, in contrast to simple pairwise atom-atom theories which fail qualitatively for layers with zero electronic energy gap. At the same time our theory is much less computationally intensive than the full RPA energy. It also gives accurate correlation energies near to the binding minimum, in contrast to Lifshitz-type theory. We apply our LRT theory successfully to graphite and to BN, and to a graphene-BN heterostructure.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Locality of correlation in density functional theory

The Hohenberg-Kohn density functional was long ago shown to reduce to the Thomas-Fermi approximation in the non-relativistic semiclassical (or large-$Z$) limit for all matter, i.e, the kinetic energy becomes local. Exchange also becomes local in this limit. Numerical data on the correlation energy of atoms supports the conjecture that this is also true for correlation, but much less relevant to atoms. We illustrate how expansions around large particle number are equivalent to local density approximations and their strong relevance to density functional approximations. Analyzing highly accurate atomic correlation energies, we show that the correlation energy tends to $-A_c Z ln Z + B_c Z$ as $Z$ tends to infinity, where $Z$ is the atomic number, $A_c$ is known, and we estimate $B_c$ to be about 37 millihartrees. The local density approximation yields $A_c$ exactly, but a very incorrect value for $B_c$, showing that the local approximation is less relevant for correlation alone. This limit is a benchmark for the non-empirical construction of density functional approximations. We conjecture that, beyond atoms, the leading correction to the local density approximation in the large-$Z$ limit generally takes this form, but with $B_c$ a functional of the TF density for the system. The implications for construction of approximate density functionals are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A comment on "Interlayer interactions in graphites" [Chen et al., Sci. Rep. 3, 3046 (2013)]

Determining the material properties of layered systems like graphite and bigraphene from \emph{ab initio} calculations is very difficult. This is mostly due to the complex van der Waals forces which help bind the layers. Recently, Chen~\emph{et al.} [Chen et al., Sci. Rep. 3, 3046 (2013)] reported a novel approach for extracting geometry dependent energetic properties of general, layered graphitic systems from periodic graphite calculations on AA, AB and ABC graphite. Unfortunately, their analysis suffered from a number of technical and theoretical flaws which make their results unreliable for predicting energetic properties. We propose that their conclusions in this regard should be reassessed, or reanalysed using more appropriate van der Waals theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Atomic correlation energies and the generalized gradient approximation

Careful extrapolation of atomic correlation energies suggests that $E_c$ tends to $-AZ\log{Z} + BZ$ as $Z$ tends to infinity, where $Z$ is the atomic number, $A$ is known, and $B$ is about 38 milliHartrees. The coefficients roughly agree with those of the high-density limit of the real-space construction of the generalized gradient approximation. An asymptotic coefficient, missed by previous derivations, is included in a revised approximation. The exchange is also corrected, reducing atomic errors considerably.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

How many-body effects modify the van der Waals interaction between graphene sheets

Undoped graphene (Gr) sheets at low temperatures are known, via Random Phase Approximation (RPA) calculations, to exhibit unusual van der Waals (vdW) forces. Here we show that graphene is the first known system where effects beyond the RPA make qualitative changes to the vdW force. For large separations, $D \gtrsim 10$nm where only the $π_z$ vdW forces remain, we find the Gr-Gr vdW interaction is substantially reduced from the RPA prediction. Its $D$ dependence is very sensitive to the form of the long-wavelength many-body enhancement of the velocity of the massless Dirac fermions, and may provide independent confirmation of the latter via direct force measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kohn-Sham potentials in exact density-functional theory at non-integer electron numbers

Within exact electron density-functional theory, we investigate Kohn-Sham (KS) potentials, orbital energies, and non-interacting kinetic energies of the fractional ions of Li, C and F. We use quantum Monte Carlo densities as input, which are then fitted, interpolated at non-integer electron numbers $N$, and inverted to produce accurate KS potentials $v_s^N(r)$. We study the dependence of the KS potential on $N$, and in particular we numerically reproduce the theoretically predicted spatially constant discontinuity of $v_s^N(r)$ as $N$ passes through an integer. We further show that, for all the cases considered, the inner orbital energies and the non-interacting kinetic energy are nearly piecewise linear functions of $N$. This leads us to propose a simple approximation of the KS potential $v_s^N(r)$ at any fractional electron number $N$ which uses only quantities of the systems with the adjacent integer electron numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The flexible nature of exchange, correlation and Hartree physics: resolving "delocalization" errors in a 'correlation free' density functional

By exploiting freedoms in the definitions of 'correlation', 'exchange' and 'Hartree' physics in ensemble systems we better generalise the notion of 'exact exchange' (EXX) to systems with fractional occupations functions of the frontier orbitals, arising in the dissociation limit of some molecules. We introduce the Linear EXX ("LEXX") theory whose pair distribution and energy are explicitly \emph{piecewise linear} in the occupations $f^σ_{i}$. {\hi}We provide explicit expressions for these functions for frontier $s$ and $p$ shells. Used in an optimised effective potential (OEP) approach it yields energies bounded by the piecewise linear 'ensemble EXX' (EEXX) energy and standard fractional optimised EXX energy: $E^{EEXX}\leq E^{LEXX} \leq E^{EXX}$. Analysis of the LEXX explains the success of standard OEP methods for diatoms at large spacing, and why they can fail when both spins are allowed to be non-integer so that "ghost" Hartree interactions appear between \emph{opposite} spin electrons in the usual formula. The energy $E^{LEXX}$ contains a cancellation term for the spin ghost case. It is evaluated for H, Li and Na fractional ions with clear derivative discontinuities for all cases. The $p$-shell form reproduces accurate correlation-free energies of B-F and Al-Cl. We further test LEXX plus correlation energy calculations on fractional ions of C and F and again shows both derivative discontinuities and good agreement with exact results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Assessment of range-separated time-dependent density-functional theory for calculating C6 dispersion coefficients

We assess a variant of linear-response range-separated time-dependent density-functional theory (TDDFT), combining a long-range Hartree-Fock (HF) exchange kernel with a short-range adiabatic exchange-correlation kernel in the local-density approximation (LDA) for calculating isotropic C6 dispersion coefficients of homodimers of a number of closed-shell atoms and small molecules. This range-separated TDDFT tends to give underestimated C6 coefficients of small molecules with a mean absolute percentage error of about 5%, a slight improvement over standard TDDFT in the adiabatic LDA which tends to overestimate them with a mean absolute percentage error of 8%, but close to time-dependent Hartree-Fock which has a mean absolute percentage error of about 6%. These results thus show that introduction of long-range HF exchange in TDDFT has a small but beneficial impact on the values of C6 coefficients. It also confirms that the present variant of range-separated TDDFT is a reasonably accurate method even using only a LDA-type density functional and without adding an explicit treatment of long-range correlation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Binding and interlayer force in the near-contact region of two graphite slabs: experiment and theory

Via a novel experiment, Liu \emph{et al.} [Phys. Rev. B, {\bf 85}, 205418 (2012)] estimated the graphite binding energy, specifically the cleavage energy, an important physical property of bulk graphite. We re-examine the data analysis and note that within the standard Lennard-Jones model employed, there are difficulties in achieving internal consistency in the reproduction of the graphite elastic properties. By employing similar models which guarantee consistency with the elastic constant, we find a wide range of model dependent binding energy values from the same experimental data. We attribute some of the difficulty in the determination of the binding energy to: i) limited theoretical understanding of the van der Waals dispersion of graphite cleavage, ii) the mis-match between the strong bending stiffness of the graphite-SiO$_2$ cantilever and the weak asymptotic inter-layer forces that are integrated over to produce the binding energy. We find, however, that the data does support determination of a maximum inter-layer force that is relatively model independent. We conclude that the peak force per unit area is $1.1 \pm 0.15$GPa for cleavage, and occurs at an inter-layer spacing of $0.377\pm 0.013$nm.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dispersion corrections in graphenic systems: a simple and effective model of binding

We combine high-level theoretical and \emph{ab initio} understanding of graphite to develop a simple, parametrised force-field model of interlayer binding in graphite, including the difficult non-pairwise-additive coupled-fluctuation dispersion interactions. The model is given as a simple additive correction to standard density functional theory (DFT) calculations, of form $ΔU(D)=f(D)[U^{vdW}(D)-U^{DFT}(D)]$ where $D$ is the interlayer distance. The functions are parametrised by matching contact properties, and long-range dispersion to known values, and the model is found to accurately match high-level \emph{ab initio} results for graphite across a wide range of $D$ values. We employ the correction on the difficult bigraphene binding and graphite exfoliation problems, as well as lithium intercalated graphite LiC$_6$. We predict the binding energy of bigraphene to be 0.27 J/m^2, and the exfoliation energy of graphite to be 0.31 J/m^2, respectively slightly less and slightly more than the bulk layer binding energy 0.295 J/m^2/layer. Material properties of LiC$_6$ are found to be essentially unchanged compared to the local density approximation. This is appropriate in view of the relative unimportance of dispersion interactions for LiC$_6$ layer binding.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Beyond the RPA on the cheap: improved correlation energies with the efficient "Radial Exchange Hole" kernel

The "ACFD-RPA" correlation energy functional has been widely applied to a variety of systems to successfully predict energy differences, and less successfully predict absolute correlation energies. Here we present a parameter-free exchange-correlation kernel that systematically improves absolute correlation energies, while maintaining most of the good numerical properties that make the ACFD-RPA numerically tractable. The "RXH" kernel is constructed to approximate the true exchange kernel via a carefully weighted, easily computable radial averaging. Correlation energy errors of atoms with two to eighteen electrons show a thirteenfold improvement over the RPA and a threefold improvement over the related "PGG" kernel, for a mean absolute error of 13mHa or 5%. The average error is small compared to all but the most difficult to evaluate kernels. van der Waals $C_6$ coefficients are less well predicted, but still show improvements on the RPA, especially for highly polarisable Li and Na.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Correlation energies beyond the random-phase approximation: ISTLS applied to spherical atoms and ions

The inhomogeneous Singwi, Tosi, Land and Sjolander (ISTLS) correlation energy functional of Dobson, Wang and Gould [PRB {\bf 66} 081108(R) (2008)] has proved to be excellent at predicting correlation energies in semi-homogeneous systems, showing promise as a robust `next step' fifth-rung functional by using dynamic correlation to go beyond the limitations of the direct random-phase approximation (dRPA), but with similar numerical scaling with system size. In this work we test the functional on fourteen spherically symmetric, neutral and charged atomic systems and find it gives excellent results (within 2mHa/$e^-$ except Be) for the absolute correlation energies of the neutral atoms tested, and good results for the ions (within 4mHa/$e^-$). In all cases it performs better than the dRPA. When combined with the previous successes, these new results point to the ISTLS functional being a prime contender for high-accuracy, benchmark DFT correlation energy calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantum Continuum Mechanics Made Simple

In this paper we further explore and develop the quantum continuum mechanics (CM) of [Tao \emph{et al}, PRL{\bf 103},086401] with the aim of making it simpler to use in practice. Our simplifications relate to the non-interacting part of the CM equations, and primarily refer to practical implementations in which the groundstate stress tensor is approximated by its Kohn-Sham version. We use the simplified approach to directly prove the exactness of CM for one-electron systems via an orthonormal formulation. This proof sheds light on certain physical considerations contained in the CM theory and their implication on CM-based approximations. The one-electron proof then motivates an approximation to the CM (exact under certain conditions) expanded on the wavefunctions of the Kohn-Sham (KS) equations. Particular attention is paid to the relationships between transitions from occupied to unoccupied KS orbitals and their approximations under the CM. We also demonstrate the simplified CM semi-analytically on an example system.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Efficient, long-range correlation from occupied wavefunctions only

We use continuum mechanics [Tao \emph{et al}, PRL{\bf 103},086401] to approximate the dynamic density response of interacting many-electron systems. Thence we develop a numerically efficient exchange-correlation energy functional based on the Random Phase Approximation (dRPA). The resulting binding energy curve $E(D)$ for thin parallel metal slabs at separation $D$ better agrees with full dRPA calculations than does the Local Density Approximation. We also reproduce the correct non-retarded van der Waals (vdW) power law $E(D)\aeq -C_{5/2}D^{-5/2}$ as $D\to\infty$, unlike most vdW functionals.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

van der Waals dispersion power laws for cleavage, exfoliation and stretching in multi-scale, layered systems

Layered and nanotubular systems that are metallic or graphitic are known to exhibit unusual dispersive van der Waals (vdW) power laws under some circumstances. In this letter we investigate the vdW power laws of bulk and finite layered systems and their interactions with other layered systems and atoms in the electromagnetically non-retarded case. The investigation reveals substantial difference between `cleavage' and `exfoliation' of graphite and metals where cleavage obeys a $C_2 D^{-2}$ vdW power law while exfoliation obeys a $C_3 \log(D/D_0) D^{-3}$ law for graphitics and a $C_{5/2} D^{-5/2}$ law for layered metals. This leads to questions of relevance in the interpretation of experimental results for these systems which have previously assumed more trival differences. Furthermore we gather further insight into the effect of scale on the vdW power laws of systems that simultaneously exhibit macroscopic and nanoscopic dimensions. We show that, for metallic and graphitic layered systems, the known "unusual" power laws can be reduced to standard or near standard power laws when the effective scale of one or more dimension is changed. This allows better identification of the systems for which the commonly employed `sum of $C_6 D^{-6}$' type vdW methods might be valid such as layered bulk to layered bulk and layered bulk to atom.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

A theoretical and semiemprical correction to the long-range dispersion power law of stretched graphite

In recent years intercalated and pillared graphitic systems have come under increasing scrutiny because of their potential for modern energy technologies. While traditional \emph{ab initio} methods such as the LDA give accurate geometries for graphite they are poorer at predicting physicial properties such as cohesive energies and elastic constants perpendicular to the layers because of the strong dependence on long-range dispersion forces. `Stretching' the layers via pillars or intercalation further highlights these weaknesses. We use the ideas developed by [J. F. Dobson et al, Phys. Rev. Lett. {\bf 96}, 073201 (2006)] as a starting point to show that the asymptotic $C_3 D^{-3}$ dependence of the cohesive energy on layer spacing $D$ in bigraphene is universal to all graphitic systems with evenly spaced layers. At spacings appropriate to intercalates, this differs from and begins to dominate the $C_4 D^{-4}$ power law for dispersion that has been widely used previously. The corrected power law (and a calculated $C_3$ coefficient) is then unsuccesfully employed in the semiempirical approach of [M. Hasegawa and K. Nishidate, Phys. Rev. B {\bf 70}, 205431 (2004)] (HN). A modified, physicially motivated semiempirical method including some $C_4 D^{-4}$ effects allows the HN method to be used successfully and gives an absolute increase of about $2-3%$ to the predicted cohesive energy, while still maintaining the correct $C_3 D^{-3}$ asymptotics.

0 citations0 reviewsNo code linkedOpen access