DFT+U and quantum Monte Carlo study of electronic and optical properties of AgNiO2 and AgNi1-xCoxO2 delafossite.

As the only semimetallic d10-based delafossite, AgNiO2 has received a great deal of attention due to both its unique semimetallicity and its antiferromagnetism in the NiO2 layer that is coupled with a lattice distortion. In contrast, other delafossites such as AgCoO2 are insulating. Here we study how the electronic structure of AgNi1-xCoxO2 alloys vary with Ni/Co concentration, in order to investigate the electronic properties and phase stability of the intermetallics. While the electronic and magnetic structure of delafossites have been studied using density functional theory (DFT), earlier studies have not included corrections for strong on-site Coulomb interactions. In order to treat these interactions accurately, in this study we use Quantum Monte Carlo (QMC) simulations to obtain accurate estimates for the electronic and magnetic properties of AgNiO2. By comparison to DFT results we show that these electron correlations are critical to account for. We show that Co doping on the magnetic Ni sites results in a metal-insulator transition near x ∼0.33, and reentrant behavior near x ∼ 0.66.

Hydrogen separation with a graphenylene monolayer: Diffusion Monte Carlo study.

We performed fixed-node diffusion Monte Carlo (DMC) calculations to investigate structural and energetic properties of graphenylene (GPNL), a two-dimensional network of sp2-bonded carbon atoms with large near-circular pores, and its H2 separation performance for gas mixtures. We have found that the energetic stability of a GPNL monolayer is comparable to that of γ-graphyne, as evidenced by its large cohesive energy of 6.755(3) eV/atom. Diffusion barriers of several gas molecules, including hydrogen, through a GPNL membrane were determined from the analysis of their adsorption energies depending on the adsorption distance, which led to our estimation for hydrogen selectivity with respect to other target molecules. DMC hydrogen selectivity of a GPNL monolayer was found to be exceptionally high at 300 K, as high as 1010-1011 against CO and N2 gases. This, along with high hydrogen permeance due to its generic pore structure, leads us to conclude that GPNL is a promising membrane to be used as a high-performance hydrogen separator from gas mixtures. We find that when compared to our DMC results, DFT calculations tend to overestimate H2 selectivity, which is mostly due to their inaccurate description of short-range repulsive interactions.

The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations.

In this work, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew-Burke-Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd-Scuseria-Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.

Metastable Metallic Phase of a Bilayer Blue Phosphorene Induced by Interlayer Bonding and Intralayer Charge Redistributions.

We have carried out diffusion Monte Carlo calculations for an A1B-1-stacked bilayer blue phosphorene to find that it undergoes a semiconductor-metal transition as the interlayer distance decreases. While the most stable bilayer structure is a semiconducting one with two monolayers coupled through a weak van der Waals interaction, the metallic bilayer at a shorter interlayer distance is found to be only metastable. This is in contrast to a recent theoretical prediction based on a random phase approximation that the metallic phase would be the most stable bilayer configuration of blue phosphorene. Our analysis of charge density distributions reveals that the metastable metallic phase is induced by interlayer chemical bonding and intralayer charge redistributions. This study enriches our understanding of interlayer binding of a blue phosphorene and contributes to the establishment of correct energetic order between its different phases, which will be essential in devising an experimental pathway for a metallic phosphorene.

Importance of Van der Waals Interactions in Hydrogen Adsorption on a Silicon-carbide Nanotube Revisited with vdW-DFT and Quantum Monte Carlo.

Density functional theory (DFT) is a valuable tool for calculating adsorption energies toward designing materials for hydrogen storage. However, dispersion forces being absent from the local/semi-local theory, it remains unclear as to how the consideration of van der Waals (vdW) interactions affects such calculations. For the first time, we applied diffusion Monte Carlo (DMC) to evaluate the adsorption characteristics of a hydrogen molecule on a (5,5) armchair silicon-carbide nanotube (H2-SiCNT). Within the DFT framework, we benchmarked various exchange-correlation functionals, including those recently developed for treating dispersion or vdW interactions. We found that the vdW-corrected DFT methods agree well with DMC, whereas the local (semilocal) functional significantly over (under)-binds. Furthermore, we fully optimized the H2-SiCNT geometry within the DFT framework and investigated the correlation between the structure and charge density. The vdW contribution to the adsorption was found to be non-negligible at ∼1 kcal/mol per hydrogen molecule, which amounts to 9-29% of the ideal adsorption energy required for hydrogen storage applications.

Adsorption of a single Pt atom on graphene: spin crossing between physisorbed triplet and chemisorbed singlet states.

Diffusion Monte Carlo (DMC) calculations have been performed to study the adsorption of a single Pt atom on pristine graphene. We obtain the adsorption energy curves of a single Pt atom adsorbed at three different adsorption sites (bridge, on-top, hollow) as functions of the vertical distance from a graphene surface for both spin singlet and triplet states. The bridge-site adsorption in a singlet spin state is found to be energetically most stable, which is consistent with previous theoretical predictions. As the Pt atom moves away from a graphene surface, spin triplet states are favored over spin singlet states for all three adsorption sites, reflecting that the ground state of an isolated Pt atom is in a spin triplet state. Furthermore, our DMC calculations reveal local-minimum features in the triplet region which is understood to be due to van der Waals interaction between the Pt atom and graphene. This provides a comprehensive understanding for a spin crossing from a physisorbed triplet state to a chemisorbed singlet state in the adsorption process of a single Pt atom on graphene.

Quantum Monte Carlo benchmarking of large noncovalent complexes in the L7 benchmark set.

We have used diffusion Monte Carlo (DMC) to perform calculations on the L7 benchmark set. DMC is a stochastic numerical integration scheme in real-space and part of a larger set of quantum Monte Carlo methods. The L7 set was designed to test the ability of electronic structure methods to include dispersive interactions. While the agreement between DMC and quantum-chemical state-of-the-art methods is excellent for some of the structures, there are significant differences in others. In contrast to wavefunction-based quantum chemical methods, DMC is a first-principle many-body method with the many-body wavefunction evolving in real space. It includes explicitly all electron-electron interactions and is relatively insensitive to the size of the basis set.

Competition between Hückel's Rule and Jahn-Teller Distortion in Small Carbon Rings: A Quantum Monte Carlo Study.

Structural properties and energetics of carbon rings are studied with the diffusion Monte Carlo (DMC) method. Our DMC-based geometry optimization reveals that both polyynic C4n and cumulenic C4n + 2 rings exhibit bond length alternations for n ≥ 3, which is understood to be due to Jahn-Teller distortions. The bond length alternation even in a cumulenic (4n + 2) carbon ring was experimentally observed in a recently synthesized C18 molecule. From a comparison of the DMC cohesive energies of C4n with those of C4n + 2, we present a comprehensive picture of the competition between Hückel's rule and Jahn-Teller distortion in small carbon rings; the former is more dominant than the latter for n < 5 where C4n + 2 rings are more stable than C4n, while C4n rings are as stable as C4n + 2 for n < 5 where dimerization effects due to Jahn-Teller distortion are more important.

Giant Anisotropy of Gilbert Damping in Epitaxial CoFe Films.

Tailoring Gilbert damping of metallic ferromagnetic thin films is one of the central interests in spintronics applications. Here we report a giant Gilbert damping anisotropy in epitaxial Co_{50}Fe_{50} thin films with a maximum-minimum damping ratio of 400%, determined by broadband spin-torque ferromagnetic resonance as well as inductive ferromagnetic resonance. We conclude that the origin of this damping anisotropy is the variation of the spin orbit coupling for different magnetization orientations in the cubic lattice, which is further corroborated from the magnitude of the anisotropic magnetoresistance in Co_{50}Fe_{50}.

QMCPACK: an open source ab initio quantum Monte Carlo package for the electronic structure of atoms, molecules and solids.

QMCPACK is an open source quantum Monte Carlo package for ab initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wavefunctions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary-field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit and graphical processing unit systems. We detail the program's capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://qmcpack.org.

Nature of Interlayer Binding and Stacking of sp-sp2 Hybridized Carbon Layers: A Quantum Monte Carlo Study.

α-Graphyne is a two-dimensional sheet of sp-sp2 hybridized carbon atoms in a honeycomb lattice. While the geometrical structure is similar to that of graphene, the hybridized triple bonds give rise to electronic structure that is different from that of graphene. Similar to graphene, α-graphyne can be stacked in bilayers with two stable configurations, but the different stackings have very different electronic structures: one is predicted to have gapless parabolic bands, and the other, a tunable band gap which is attractive for applications. In order to realize applications, it is crucial to understand which stacking is more stable. This is difficult to model, as the stability is a result of weak interlayer van der Waals interactions which are not well captured by density functional theory (DFT). We have used quantum Monte Carlo simulations that accurately include van der Waals interactions to calculate the interlayer binding energy of bilayer graphyne and to determine its most stable stacking mode. Our results show that interlayer bindings of sp- and sp2-bonded carbon networks are significantly underestimated in a Kohn-Sham DFT approach, even with an exchange-correlation potential corrected to include, in some approximation, van der Waals interactions. Finally, our quantum Monte Carlo calculations reveal that the interlayer binding energy difference between the two stacking modes is only 0.9(4) meV/atom. From this we conclude that the two stable stacking modes of bilayer α-graphyne are almost degenerate with each other, and both will occur with about the same probability at room temperature unless there is a synthesis path that prefers one stacking over the other.

Cohesion energetics of carbon allotropes: quantum Monte Carlo study.

We have performed quantum Monte Carlo calculations to study the cohesion energetics of carbon allotropes, including sp(3)-bonded diamond, sp(2)-bonded graphene, sp-sp(2) hybridized graphynes, and sp-bonded carbyne. The computed cohesive energies of diamond and graphene are found to be in excellent agreement with the corresponding values determined experimentally for diamond and graphite, respectively, when the zero-point energies, along with the interlayer binding in the case of graphite, are included. We have also found that the cohesive energy of graphyne decreases systematically as the ratio of sp-bonded carbon atoms increases. The cohesive energy of γ-graphyne, the most energetically stable graphyne, turns out to be 6.766(6) eV/atom, which is smaller than that of graphene by 0.698(12) eV/atom. Experimental difficulty in synthesizing graphynes could be explained by their significantly smaller cohesive energies. Finally, we conclude that the cohesive energy of a newly proposed graphyne can be accurately estimated with the carbon-carbon bond energies determined from the cohesive energies of graphene and three different graphynes considered here.

Effects of 3He impurities on the superfluid response of the 4He monolayer on a C20 molecule.

The path-integral Monte Carlo calculations have been performed to investigate the effects of (3)He impurities on structural and superfluid properties of the (4)He monolayer on a single C(20) molecule. According to our previous study, the helium monolayer exhibits different quantum states for different numbers of (4)He adatoms and is completed to form a commensurate solid where nanoscale supersolidity can be realized through the activation of mobile vacancy states. We first observe that different structures for different numbers of helium atoms are mostly preserved with the replacement of a few (4)He atoms with the same number of (3)He atoms, whether the helium layer is a fluid or a solid. However, the substitution of (3)He impurities is found to have different effects on the superfluid response of the helium layer, depending on its quantum state. For a partially-filled fluid layer the superfluid fraction decreases monotonically with the increasing (3)He concentration, which can be understood in terms of the suppression of exchange couplings among (4)He atoms due to the presence of (3)He impurities. On the other hand, the substitution of a few (3)He impurity atoms may increase the superfluid fraction of a near-complete monolayer that is in a crystalline solid state. The enhancement of superfluidity in a solid layer is interpreted to be due to interstitial and vacancy defects promoted by larger quantum fluctuations of lighter (3)He atoms. This provides strong evidence that the (4)He monolayer on C(20) shows the vacancy-based supersolidity near its completion.

Commensurate-incommensurate transition of 4He adsorbed on a single C60 molecule.

Path-integral Monte Carlo calculations have been performed to study (4)He adsorption on a single C(60) molecule. Helium corrugations on the fullerene molecular surface are incorporated with the (4)He-C(60) interaction described by the sum of all (4)He-C interatomic pair potentials. Radial density distributions show a layer-by-layer growth of (4)He with the first adlayer being located at a distance of ~6.3 Å from the center of the C(60) molecule. The monolayer shows different quantum states as the number of (4)He adatoms N varies. For N = 32, we find a commensurate solid, with each of the 32 adsorption sites on the molecular surface being occupied by a single (4)He atom. Various domain-wall structures are observed as more (4)He atoms are added and the first layer crystallizes into an incommensurate solid when it is completely filled. This commensurate-incommensurate transition of the helium monolayer is found to be accompanied by re-entrant superfluid response at a low temperature of 0.31 K with the superfluidity being totally quenched at N = 32, 44, and 48. Finally, the different quantum states observed in the helium monolayer around C(60) are compared with phase diagrams proposed for the corresponding layer on a graphite surface.