Assessment of two hybrid van der Waals density functionals for covalent and non-covalent binding of molecules

Two hybrid van der Waals density functionals (vdW-DFs) are constructed using 25%, Fock exchange with i) the consistent-exchange vdW-DF-cx functional and ii) with the vdW-DF2 functional. The ability to describe covalent and non-covalent binding properties of molecules are assessed. For properties related to covalent binding, atomization energies (G2-1 set), molecular reaction energies (G2RC set), as well as ionization energies (G21IP set) are benchmarked against experimental reference values. We find that hybrid-vdW-DF-cx yields results that are rather similar to those of the standard non-empirical hybrid PBE0 [JCP 110, 6158 (1996)]. Hybrid vdW-DF2 follows somewhat different trends, showing on average significantly larger deviations from the reference energies, with a MAD of 14.5 kcal/mol for the G2-1 set. Non-covalent binding properties of molecules are assessed using the S22 benchmark set of non-covalently bonded dimers and the X40 set of dimers of small halogenated molecules, using wavefunction-based quantum chemistry results for references. For the S22 set, hybrid-vdW-DF-cx performs better than standard vdW-DF-cx for the mostly hydrogen-bonded systems. Hybrid-vdW-DF2 offers a slight improvement over standard vdW-DF2. Similar trends are found for the X40 set, with hybrid-vdW-DF-cx performing particularly well for binding involving the strongly polar hydrogen halides, but poorly for systems with tiny binding energies. Our study of the X40 set reveals both the potential of mixing Fock exchange with vdW-DF, but also highlights shortcomings of the hybrids constructed here. The solid performance of hybrid-vdW-DF-cx for covalent-bonded systems, as well as the strengths and issues uncovered for non-covalently bonded systems, makes this study a good starting point for developing even more precise hybrid vdW-DFs

Nonequilibrium thermodynamics of interacting tunneling transport: variational grand potential, density-functional formulation, and nature of steady-state forces

The standard formulation of tunneling transport rests on an open-boundary modeling. There, conserving approximations to nonequilibrium Green function or quantum-statistical mechanics provide consistent but computational costly approaches; alternatively, use of density-dependent ballistic-transport calculations [e.g., Phys. Rev. B 52, 5335 (1995)], here denoted `DBT', provide computationally efficient (approximate) atomistic characterizations of the electron behavior but has until now lacked a formal justification. This paper presents an exact, variational nonequilibrium thermodynamic theory for fully interacting tunneling and provides a rigorous foundation for frozen-nuclei DBT calculations as a lowest order approximation to an exact nonequilibrium thermodynamics density functional evaluation. The theory starts from the complete electron nonequilibrium quantum statistical mechanics and I identify the operator for the nonequilibrium Gibbs free energy. I demonstrate a minimal property of a functional for the nonequilibrium thermodynamic grand potential which thus uniquely identifies the solution as the exact nonequilibrium density matrix. I also show that a uniqueness-of-density proof from a closely related study [Phys. Rev. B 78, 165109 (2008)] makes it possible to provide a single-particle formulation based on universal electron-density functionals. I illustrate a formal evaluation of the thermodynamics grand potential value which is closely related to the variation in scattering phase shifts and hence to Friedel density oscillations. This paper also discusses the difference between the here-presented exact thermodynamics forces and the often-used electrostatic forces. Finally the paper documents an inherent adiabatic nature of the thermodynamics forces and observes that these are suited for a nonequilibrium implementation of the Born-Oppenheimer approximation.Comment: 37 pages, 3 Figure

van der Waals density functionals built upon the electron-gas tradition: Facing the challenge of competing interactions

The theoretical description of sparse matter attracts much interest, in particular for those ground-state properties that can be described by density functional theory (DFT). One proposed approach, the van der Waals density functional (vdW-DF) method, rests on strong physical foundations and offers simple yet accurate and robust functionals. A very recent functional within this method called vdW-DF-cx [K. Berland and P. Hyldgaard, Phys. Rev. B 89, 035412] stands out in its attempt to use an exchange energy derived from the same plasmon-based theory from which the nonlocal correlation energy was derived. Encouraged by its good performance for solids, layered materials, and aromatic molecules, we apply it to several systems that are characterized by competing interactions. These include the ferroelectric response in PbTiO$_3$, the adsorption of small molecules within metal-organic frameworks (MOFs), the graphite/diamond phase transition, and the adsorption of an aromatic-molecule on the Ag(111) surface. Our results indicate that vdW-DF-cx is overall well suited to tackle these challenging systems. In addition to being a competitive density functional for sparse matter, the vdW-DF-cx construction presents a more robust general purpose functional that could be applied to a range of materials problems with a variety of competing interactions

Understanding adhesion at as-deposited interfaces from ab initio thermodynamics of deposition growth: thin-film alumina on titanium carbide

We investigate the chemical composition and adhesion of chemical vapour deposited thin-film alumina on TiC using and extending a recently proposed nonequilibrium method of ab initio thermodynamics of deposition growth (AIT-DG) [Rohrer J and Hyldgaard P 2010 Phys. Rev. B 82 045415]. A previous study of this system [Rohrer J, Ruberto C and Hyldgaard P 2010 J. Phys.: Condens. Matter 22 015004] found that use of equilibrium thermodynamics leads to predictions of a non-binding TiC/alumina interface, despite the industrial use as a wear-resistant coating. This discrepancy between equilibrium theory and experiment is resolved by the AIT-DG method which predicts interfaces with strong adhesion. The AIT-DG method combines density functional theory calculations, rate-equation modelling of the pressure evolution of the deposition environment and thermochemical data. The AIT-DG method was previously used to predict prevalent terminations of growing or as-deposited surfaces of binary materials. Here we extent the method to predict surface and interface compositions of growing or as-deposited thin films on a substrate and find that inclusion of the nonequilibrium deposition environment has important implications for the nature of buried interfaces.Comment: 8 pages, 6 figures, submitted to J. Phys.: Condens. Matte

Minimum Thermal Conductivity of Superlattices

The phonon thermal conductivity of a multilayer is calculated for transport perpendicular to the layers. There is a cross over between particle transport for thick layers to wave transport for thin layers. The calculations shows that the conductivity has a minimum value for a layer thickness somewhat smaller then the mean free path of the phonons.Comment: new results added, to appear in PR

Implication of the overlap representation for modelling generalized parton distributions

Based on a field theoretically inspired model of light-cone wave functions, we derive valence-like generalized parton distributions and their double distributions from the wave function overlap in the parton number conserved s-channel. The parton number changing contributions in the t-channel are restored from duality. In our construction constraints of positivity and polynomiality are simultaneously satisfied and it also implies a model dependent relation between generalized parton distributions and transverse momentum dependent parton distribution functions. The model predicts that the t-behavior of resulting hadronic amplitudes depends on the Bjorken variable x_Bj. We also propose an improved ansatz for double distributions that embeds this property.Comment: 15 pages, 8 eps figure

The effect of inelastic processes on tunneling

We study an electron that interacts with phonons or other linear or nonlinear excitations as it resonantly tunnels. The method we use is based on mapping a many-body problem in a large variational space exactly onto a one-body problem. The method is conceptually simpler than previous Green's function approaches, and allows the essentially exact numerical solution of much more general problems. We solve tunneling problems with transverse channels, multiple sites coupled to phonons, and multiple phonon degrees of freedom and excitations.Comment: 12 pages, REVTex, 4 figures in compressed tar .ps forma

Zero Frequency Current Noise for the Double Tunnel Junction Coulomb Blockade

We compute the zero frequency current noise numerically and in several limits analytically for the coulomb blockade problem consisting of two tunnel junctions connected in series. At low temperatures over a wide range of voltages, capacitances, and resistances it is shown that the noise measures the variance in the number of electrons in the region between the two tunnel junctions. The average current, on the other hand, only measures the mean number of electrons. Thus, the noise provides additional information about transport in these devices which is not available from measuring the current alone.Comment: 33 pages, 10 figure

Phylogeography reveals a potential cryptic invasion in the Southern Hemisphere of Ceratophyllum demersum, New Zealand's worst invasive macrophyte

Ceratophyllum demersum (common hornwort) is presently considered the worst invasive submerged aquatic macrophyte in New Zealand. We explored the global phylogeographic pattern of the species, based on chloroplast and nuclear DNA, in order to identify the origin of the invasive populations in New Zealand and to clarify if there were multiple introductions. The phylogeographic study identified geographically differentiated gene pools in North America, tropical Asia, Australia, and South Africa, likely native to these regions, and a recent dispersal event of a Eurasian-related haplotype to North America, New Zealand, Australia, and South Africa. At least two different invasive genotypes of this Eurasian-related haplotype have been found in New Zealand. One genotype is closely related to genotypes in Australia and South Africa, while we could not trace the closest relatives of the other genotype within our C. demersum sample set. Contrasting spectra of genetic distances in New Zealand and in a region within the native range (Denmark), suggest that the invasive population was founded by vegetative reproduction, seen as low genetic distances among genotypes. We also discovered the introduction of the same Eurasian-related haplotype in Australia and South Africa and that a cryptic invasion may be occurring in these continents

Accurate Nonempirical Range-Separated Hybrid van der Waals Density Functional for Complex Molecular Problems, Solids, and Surfaces

We introduce a new, general-purpose, range-separated hybrid van der Waals density functional termed vdW-DF2-ahbr within the nonempirical vdW-DF method [Hyldgaard, et al. J. Phys. Condens. Matter 32, 393001 (2020)]. It combines a correlation from vdW-DF2 with a screened Fock exchange that is fixed by a new model of exchange effects in the density-explicit vdW-DF-b86r or rev-vdW-DF2 functional [Hamada, Phys. Rev. B 89, 121103(R) (2014)]. The new vdW-DF2-ahbr prevents spurious exchange binding and has a small-density-gradient form set from many-body perturbation analysis. It is accurate for bulk as well as layered materials, and it systematically and significantly improves the performance of the present vdW-DFs for molecular problems. Importantly, vdW-DF2-ahbr also outperforms present-standard (dispersion-corrected) range-separated hybrids on a broad collection of noncovalent-interaction benchmark sets, while at the same time successfully mitigating the density-driven errors that often affect the description of molecular transition states and isomerization calculations. vdW-DF2-ahbr furthermore improves on state-of-the-art density-functional-theory approaches by succeeding at challenging problems. For example, it (1) correctly predicts both the substrate structure and the site preference for CO adsorption on Pt(111), (2) it outperforms existing nonempirical vdW-DFs for the description of CO2 adsorption in both a functionalized and in a simple metal-organic framework, and (3) it is highly accurate for the set of base-pair interactions in a model of DNA assembly