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

Lubos Mitas

Lubos Mitas 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.str-elphysics.atom-phquant-phcond-mat.othercond-mat.quant-gashep-latnucl-th

Catalog footprint

What is connected

22works
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

22 published item(s)

preprint2022arXiv

A new generation of effective core potentials from correlated and spin-orbit calculations: selected heavy elements

We introduce new correlation consistent effective core potentials (ccECPs) for the elements I, Te, Bi, Ag, Au, Pd, Ir, Mo, and W with $4d$, $5d$, $6s$ and $6p$ valence spaces. These ccECPs are given as a sum of spin-orbit averaged relativistic effective potential (AREP) and effective spin-orbit (SO) terms. The construction involves several steps with increasing refinements from more simple to fully correlated methods. The optimizations are carried out with objective functions that include weighted many-body atomic spectra, norm-conservation criteria, and spin-orbit splittings. Transferability tests involve molecular binding curves of corresponding hydride and oxide dimers. The constructed ccECPs are systematically better and in a few cases on par with previous effective core potential (ECP) tables on all tested criteria and provide a significant increase in accuracy for valence-only calculations with these elements. Our study confirms the importance of the AREP part in determining the overall quality of the ECP even in the presence of sizable spin-orbit effects. The subsequent quantum Monte Carlo (QMC) calculations point out the importance of accurate trial wave functions which in some cases (mid series transition elements) require treatment well beyond single-reference.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Assessing the accuracy of compound formation energies with quantum Monte Carlo

Accurately predicting the formation energy of a compound, which describes its thermodynamic stability, is a key challenge in materials physics. Here, we employ many-body quantum Monte Carlo (QMC) with single-reference trial functions to compute the formation energy of two electronically disparate compounds, the intermetallic VPt$_2$ and the semiconductor CuI, for which standard density functional theory (DFT) predictions using both the Perdew-Burke Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) density functional approximations deviate markedly from available experimental values. For VPt$_2$, we find an agreement between QMC, SCAN, and PBE0 estimates, which therefore remain in disagreement with the much less exothermic experimental value. For CuI, the QMC result agrees with neither SCAN nor PBE pointing towards DFT exchange-correlation biases, likely related to the localized Cu $3d$ electrons. Compared to the behavior of some density functional approximations within DFT, spin-averaged QMC exhibits a smaller but still appreciable deviation when compared to experiment. The QMC result is slightly improved by incorporating spin-orbit corrections for CuI and solid I$_2$, so that experiment and theory are brought into imperfect but reasonable agreement within about 120~meV/atom.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Electronic structure of $\boldsymbolα$-RuCl$_3$ by fixed-node and fixed-phase diffusion Monte Carlo methods

Layered material $α$-RuCl$_3$ has caught wide attention due to its possible realization of Kitaev's spin liquid and its electronic structure that involves the interplay of electron-electron correlations and spin-orbit effects. Several DFT$+U$ studies have suggested that both electron-electron correlations and spin-orbit effects are crucial for accurately describing the band gap. This work studies the importance of these two effects using fixed-node and fixed-phase diffusion Monte Carlo calculations both in spin-averaged and explicit spin-orbit formalisms. In the latter, the Slater-Jastrow trial function is constructed from two-component spin-orbitals using our recent quantum Monte Carlo (QMC) developments and thoroughly tested effective core potentials. Our results show that the gap in the ideal crystal is already accurately described by the spin-averaged case, with the dominant role being played by the magnetic ground state with significant exchange and electron correlation effects. We find qualitative agreement between hybrid DFT, DFT+$U$, and QMC. In addition, QMC results agree very well with available experiments, and we identify the values of exact Fock exchange mixing that provide comparable gaps. Explicit spin-orbit QMC calculations reveal that the effect of spin-orbit coupling on the gap is minor, of the order of 0.2 eV, which corresponds to the strength of the spin-orbit of the Ru atom.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Weighted nodal domain averages of eigenstates for quantum Monte Carlo and beyond

We study the nodal properties of many-body eigenstates of stationary Schrödinger equation that affect the accuracy of real-space quantum Monte Carlo calculations. In particular, we introduce weighted nodal domain averages that provide a new probe of nodal surfaces beyond the usual expectations. Particular choices for the weight function reveal, for example, that the difference between two arbitrary fermionic eigenvalues is given by the nodal hypersurface integrals normalized by overlaps with the bosonic ground state of the given Hamiltonian. Noninteracting and fully interacting Be atom with corresponding almost exact and approximate wave functions are used to illustrate several aspects of these concepts. Variational formulations that employ different weights are proposed for prospective improvement of nodes in variational and fixed-node diffusion Monte Carlo calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accurate atomic correlation and total energies for correlation consistent effective core potentials

Very recently, we introduced a set of correlation consistent effective core potentials (ccECPs) constructed within full many-body approaches. By employing significantly more accurate correlated approaches we were able to reach a new level of accuracy for the resulting effective core Hamiltonians. We also strived for simplicity of use and easy transferability into a variety of electronic structure methods in quantum chemistry and condensed matter physics. Here, as a reference for future use, we present exact or nearly-exact total energy calculations for these ccECPs. The calculations cover H-Kr elements and are based on the state-of-the-art configuration interaction (CI), coupled-cluster (CC), and quantum Monte Carlo (QMC) calculations with systematically eliminated/improved errors. In particular, we carry out full CI/CCSD(T)/CCSDT(Q) calculations with cc-pVnZ with up to n=6 basis sets and we estimate the complete basis set limits. Using combinations of these approaches, we achieved an accuracy of $\approx$ 1-10 mHa for K-Zn atoms and $\approx$ 0.1-0.3 mHa for all other elements $-$ within about 1% or better of the ccECP total correlation energies. We also estimate the corresponding kinetic energies within the feasible limit of full CI calculations. In order to provide data for QMC calculations, we include fixed-node diffusion Monte Carlo energies for each element that give quantitative insights into the fixed-node biases for single-reference trial wave functions. The results offer a clear benchmark for future high accuracy calculations in a broad variety of correlated wave function methods such as CI and CC as well is in stochastic approaches such as real space sampling QMC.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Many-body electronic structure of LaScO$_3$ by real space quantum Monte Carlo

We present real space quantum Monte Carlo (QMC) calculations of the scandate LaScO$_3$ that proved to be challenging for traditional electronic structure approaches due to strong correlation effects resulting in inaccurate band gaps from DFT and $GW$ methods when compared with existing experimental data. Besides calculating an accurate QMC band gap corrected for supercell size biases and in agreement with numerous experiments, we also predict the cohesive energy of the crystal using the standard fixed-node QMC without any empirical or non-variational parameters. We show that promotion (optical) gap and fundamental gap agree with each other illustrating a clear absence of significant excitonic effects in the ideal crystal. We obtained these results in perfect consistency in two independent tracks that employ different basis sets (plane wave vs. localized gaussians), different codes for generating orbitals (\textsc{Quantum Espresso} vs. \textsc{Crystal}), different QMC codes (\textsc{Qmcpack} vs. \textsc{Qwalk}) and different high-accuracy pseudopotentials (ccECPs vs. Troullier-Martins) presenting the maturity and consistency of QMC methodology and tools for studies of strongly correlated problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion Quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing accuracy. Advances in real space methods include techniques for accurate computation of band gaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods including GW and density functional based techniques. To provide an improved foundation for these calculations we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum Monte Carlo for Noncovalent Interactions: An Efficient Protocol Attaining Benchmark Accuracy

Reliable theoretical predictions of noncovalent interaction energies, which are important e.g. in drug-design and hydrogen-storage applications, belong to longstanding challenges of contemporary quantum chemistry. In this respect, the fixed-node diffusion Monte Carlo (FN-DMC) is a promising alternative to the commonly used "gold standard" coupled-cluster CCSD(T)/CBS method for its benchmark accuracy and favourable scaling, in contrast to other correlated wave function approaches. This work is focused on the analysis of protocols and possible tradeoffs for FN-DMC estimations of noncovalent interaction energies and proposes an efficient yet accurate computational protocol using simplified explicit correlation terms with a favorable O(N^3) scaling. It achieves an excellent agreement (mean unsigned error ~0.2 kcal/mol) with respect to the CCSD(T)/CBS data on a number of complexes, including benzene/hydrogen,T-shape benzene dimer, stacked adenine-thymine and a set of small noncovalent complexes A24. The high accuracy and reduced computational costs predestinate the reported protocol for practical interaction energy calculations of large noncovalent complexes, where the CCSD(T)/CBS is prohibitively expensive.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum Monte Carlo with Variable Spins

We investigate the inclusion of variable spins in electronic structure quantum Monte Carlo, with a focus on diffusion Monte Carlo with Hamiltonians that include spin-orbit interactions. Following our previous introduction of fixed-phase spin-orbit diffusion Monte Carlo (FPSODMC), we thoroughly discuss the details of the method and elaborate upon its technicalities. We present a proof for an upper-bound property for complex nonlocal operators, which allows for the implementation of T-moves to ensure the variational property. We discuss the time step biases associated with our particular choice of spin representation. Applications of the method are also presented for atomic and molecular systems. We calculate the binding energies and geometry of the PbH and Sn$_2$ molecules, as well as the electron affinities of the 6$p$ row elements in close agreement with experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spin-Orbit Interactions in Electronic Structure Quantum Monte Carlo

We develop generalization of the fixed-phase diffusion Monte Carlo method for Hamiltonians which explicitly depend on particle spins such as for spin-orbit interactions. The method is formulated in zero variance manner and is similar to treatment of nonlocal operators in commonly used static- spin calculations. Tests on atomic and molecular systems show that it is very accurate, on par with the fixed-node method. This opens electronic structure quantum Monte Carlo methods to a vast research area of quantum phenomena in which spin-related interactions play an important role.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fixed-Node Diffusion Monte Carlo of Lithium Systems

We study lithium systems over a range of number of atoms, e.g., atomic anion, dimer, metallic cluster, and body-centered cubic crystal by the diffusion Monte Carlo method. The calculations include both core and valence electrons in order to avoid any possible impact by pseudo potentials. The focus of the study is the fixed-node errors, and for that purpose we test several orbital sets in order to provide the most accurate nodal hyper surfaces. We compare our results to other high accuracy calculations wherever available and to experimental results so as to quantify the the fixed-node errors. The results for these Li systems show that fixed-node quantum Monte Carlo achieves remarkably high accuracy total energies and recovers 97-99 % of the correlation energy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fixed-node errors in quantum Monte Carlo: interplay of electron density and node nonlinearities

We elucidate the origin of large differences (two-fold or more) in the fixed-node errors between the first- vs second-row systems for single-configuration trial wave functions in quantum Monte Carlo calculations. This significant difference in the fixed-node biases is studied across a set of atoms, molecules, and also Si, C solid crystals. The analysis is done over valence isoelectronic systems that share similar correlation energies, bond patterns, geometries, ground states, and symmetries. We show that the key features which affect the fixed-node errors are the differences in electron density and the degree of node nonlinearity. The findings reveal how the accuracy of the quantum Monte Carlo varies across a variety of systems, provide new perspectives on the origins of the fixed-node biases in electronic structure calculations of molecular and condensed systems, and carry implications for pseudopotential constructions for heavy elements

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Study of dipole moments of LiSr and KRb molecules by quantum Monte Carlo methods

Heteronuclear dimers are of significant interest to experiments seeking to exploit ultracold polar molecules in a number of novel ways including precision measurement, quantum computing, and quantum simulation. We calculate highly accurate Born-Oppenheimer total energies and electric dipole moments as a function of internuclear separation for two such dimers, LiSr and KRb. We apply fully-correlated, high-accuracy quantum Monte Carlo methods for evaluating these molecular properties in a many-body framework. We use small-core effective potentials combined with multi-reference Slater-Jastrow trial wave functions to provide accurate nodes for the fixed-node diffusion Monte Carlo method. For reference and comparison, we calculate the same properties with Hartree-Fock and with restricted Configuration Interaction methods, and carefully assess the impact of the recovered many-body correlations on the calculated quantities. For LiSr we find a highly nonlinear dipole moment curve, which may make this molecule's dipole moment tunable through vibrational state control.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Atomic Fermi gas at the unitary limit by quantum Monte Carlo methods: Effects of the interaction range

We calculate the ground-state properties of unpolarized two-component Fermi gas by the diffusion quantum Monte Carlo (DMC) methods. Using an extrapolation to the zero effective range of the attractive two-particle interaction, we find $E/E_{\rm free}$ to be 0.212(2), 0.407(2), 0.409(3) and 0.398(3) for 4, 14, 38 and 66 atoms, respectively. Our results indicate that the dependence of the total energy on the effective range is sizable and the extrapolation is therefore quite important. In order to test the quality of nodal surfaces and to estimate the impact of the fixed-node approximation we perform released-node DMC calculations for 4 and 14 atoms. Analysis of the released-node and the fixed-node results suggests that the main sources of the fixed-node errors are long-range correlations which are difficult to sample in the released-node approaches due to the fast growth of the bosonic noise. Besides energies, we evaluate the two-body density matrix and the condensate fraction. We find that the condensate fraction for the 66 atom system converges to 0.56(1) after the extrapolation to the zero interaction range.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Impact of the Electron Density on the Fixed-Node Errors in Quantum Monte Carlo

We analyze the effect of increasing charge density on the Fixed Node Errors in Diffusion Monte Carlo by comparing FN-DMC calculations of the total ground state energy on a 4 electron system done with a Hartree-Fock based trial wave function to calculations by the same method on the same system using a Configuration Interaction based trial wave function. We do this for several different values of nuclear charge, Z. The Fixed Node Error of a Hartree-Fock trial wave function for a 4 electron system increases linearly with increasing nuclear charge.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Precision benchmark calculations for four particles at unitarity

The unitarity limit describes interacting particles where the range of the interaction is zero and the scattering length is infinite. We present precision benchmark calculations for two-component fermions at unitarity using three different ab initio methods: Hamiltonian lattice formalism using iterated eigenvector methods, Euclidean lattice formalism with auxiliary-field projection Monte Carlo, and continuum diffusion Monte Carlo with fixed and released nodes. We have calculated the ground state energy of the unpolarized four-particle system in a periodic cube as a dimensionless fraction of the ground state energy for the non-interacting system. We obtain values 0.211(2) and 0.210(2) using two different Hamiltonian lattice representations, 0.206(9) using Euclidean lattice, and an upper bound of 0.212(2) from fixed-node diffusion Monte Carlo. Released-node calculations starting from the fixed-node result yield a decrease of less than 0.002 over a propagation of 0.4/E_F in Euclidean time, where E_F is the Fermi energy. We find good agreement among all three ab initio methods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Variational Monte Carlo for spin-orbit interacting systems

Recently, a diffusion Monte Carlo algorithm was applied to the study of spin dependent interactions in condensed matter. Following some of the ideas presented therein, and applied to a Hamiltonian containing a Rashba-like interaction, a general variational Monte Carlo approach is here introduced that treats in an efficient and very accurate way the spin degrees of freedom in atoms when spin orbit effects are included in the Hamiltonian describing the electronic structure. We illustrate the algorithm on the evaluation of the spin-orbit splittings of isolated carbon and lead atoms. In the case of the carbon atom, we investigate the differences between the inclusion of spin-orbit in its realistic and effective spherically symmetrized forms. The method exhibits a very good accuracy in describing the small energy splittings, opening the way for a systematic quantum Monte Carlo studies of spin-orbit effects in atomic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Applications of quantum Monte Carlo methods in condensed systems

The quantum Monte Carlo methods represent a powerful and broadly applicable computational tool for finding very accurate solutions of the stationary Schroedinger equation for atoms, molecules, solids and a variety of model systems. The algorithms are intrinsically parallel and are able to take full advantage of the present-day high-performance computing systems. This review article concentrates on the fixed-node/fixed-phase diffusion Monte Carlo method with emphasis on its applications to electronic structure of solids and other extended many-particle systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Electronic structure quantum Monte Carlo

Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Wave functions for quantum Monte Carlo calculations in solids: Orbitals from density functional theory with hybrid exchange-correlation functionals

We investigate how the fixed-node diffusion Monte Carlo energy of solids depends on single-particle orbitals used in Slater--Jastrow wave functions. We demonstrate that the dependence can be significant, in particular in the case of 3d transition-metal compounds, which we adopt as examples. We illustrate how exchange-correlation functionals with variable exact-exchange component can be exploited to reduce the fixed-node errors. On the basis of these results we argue that the fixed-node quantum Monte Carlo provides a variational approach for optimization of effective hamiltonians with parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Structural (B1 to B8) Phase Transition in MnO under Pressure: Comparison of All-electron and Pseudopotential Approaches

We employ the density functional theory to study a structural transition of MnO from B1 (rocksalt) to B8 (NiAs) structures that was observed experimentally at pressures around 100 GPa. We utilize all-electron description as well as norm-conserving pseudopotentials and demonstrate that these two approaches can significantly differ in quantitative predictions. We explicitly show that even small-core pseudopotentials exhibit transferability inaccuracies for quantities sensitive to the energy differences between high- and low-spin polarizations of valence electrons.

0 citations0 reviewsNo code linkedOpen access