BZPEERReading, trust and collaboration for research
Menu

Discover

SearchFind papers, people, topics and opportunities from one place.GraphOpen maps around papers, topics and saved views.

Workspaces

HomeOpen the ranked research feed.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ireneusz W. Bulik

Ireneusz W. Bulik contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phcond-mat.str-elcond-mat.mtrl-sci

Trust snapshot

Quick read

Trust 21 - Emerging
11works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2015arXiv

Can single-reference coupled cluster theory describe static correlation?

While restricted single-reference coupled cluster theory truncated to singles and doubles (CCSD) provides very accurate results for weakly correlated systems, it usually fails in the presence of static or strong correlation. This failure is generally attributed to the qualitative breakdown of the reference, and can accordingly be corrected by using a multi-determinant reference, including higher-body cluster operators in the ansatz, or allowing symmetry breaking in the reference. None of these solutions are ideal; multi-reference coupled cluster is not black box, including higher-body cluster operators is computationally demanding, and allowing symmetry breaking leads to the loss of good quantum numbers. It has long been recognized that quasidegeneracies can instead be treated by modifying the coupled cluster ansatz. The recently introduced pair coupled cluster doubles (pCCD) approach is one such example which avoids catastrophic failures and accurately models strong correlations in a symmetry-adapted framework. Here we generalize pCCD to a singlet-paired coupled cluster model (CCD0) intermediate between coupled cluster doubles and pCCD, yielding a method that possesses the invariances of the former and much of the stability of the latter. Moreover, CCD0 retains the full structure of coupled cluster theory, including a fermionic wave function, antisymmetric cluster amplitudes, and well-defined response equations and density matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Combinations of coupled cluster, density functionals, and the random phase approximation for describing static and dynamic correlation, and van der Waals interactions

Contrary to standard coupled cluster doubles (CCD) and Brueckner doubles (BD), singlet-paired analogues of CCD and BD (denoted here as CCD0 and BD0) do not break down when static correlation is present, but neglect substantial amounts of dynamic correlation. In fact, CCD0 and BD0 do not account for any contributions from multielectron excitations involving only same-spin electrons at all. We exploit this feature to add---without introducing double counting, self-interaction, or increase in cost---the missing correlation to these methods via meta-GGA density functionals (TPSS and SCAN). Furthermore, we improve upon these CCD0+DFT blends by invoking range separation: the short- and long-range correlations absent in CCD0/BD0 are evaluated with DFT and the direct random phase approximation (dRPA), respectively. This corrects the description of long-range van der Waals forces. Comprehensive benchmarking shows that the combinations presented here are very accurate for weakly correlated systems, while also providing a reasonable description of strongly correlated problems without resorting to symmetry breaking.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pair extended coupled cluster doubles

The accurate and efficient description of strongly correlated systems remains an important challenge for computational methods. Doubly occupied configuration interaction (DOCI), in which all electrons are paired and no correlations which break these pairs are permitted, can in many cases provide an accurate account of strong correlations, albeit at combinatorial computational cost. Recently, there has been significant interest in a method we refer to as pair coupled cluster doubles (pCCD), a variant of coupled cluster doubles in which the electrons are paired. This is simply because pCCD provides energies nearly identical to those of DOCI, but at mean-field computational cost (disregarding the cost of the two-electron integral transformation). Here, we introduce the more complete pair extended coupled cluster doubles (pECCD) approach which, like pCCD, has mean-field cost and reproduces DOCI energetically. We show that unlike pCCD, pECCD also reproduces the DOCI wave function with high accuracy. Moreoever, pECCD yields sensible albeit inexact results even for attractive interactions where pCCD breaks down.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

Numerical results for ground state and excited state properties (energies, double occupancies, and Matsubara-axis self energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed node approximation, unrestricted coupled cluster theory, and multi-reference projected Hartree-Fock. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electron correlation in solids via density embedding theory

Density matrix embedding theory (Phys. Rev. Lett. 109, 186404 (2012)) and density embedding theory ((Phys. Rev. B 89, 035140 (2014)) have recently been introduced for model lattice Hamiltonians and molecular systems. In the present work, the formalism is extended to the ab initio description of infinite systems. An appropriate definition of the impurity Hamiltonian for such systems is presented and demonstrated in cases of 1, 2 and 3 dimensions, using coupled cluster theory as the impurity solver. Additionally, we discuss the challenges related to disentanglement of fragment and bath states. The current approach yields results comparable to coupled cluster calculations of infinite systems even when using a single unit cell as the fragment. The theory is formulated in the basis of Wannier functions but it does not require separate localization of unoccupied bands. The embedding scheme presented here is a promising way of employing highly accurate electronic structure methods for extended systems at a fraction of their original computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Seniority-based coupled cluster theory

Doubly occupied configuration interaction (DOCI) with optimized orbitals often accurately describes strong correlations while working in a Hilbert space much smaller than that needed for full configuration interaction. However, the scaling of such calculations remains combinatorial with system size. Pair coupled cluster doubles (pCCD) is very successful in reproducing DOCI energetically, but can do so with low polynomial scaling ($N^3$, disregarding the two-electron integral transformation from atomic to molecular orbitals). We show here several examples illustrating the success of pCCD in reproducing both the DOCI energy and wave function, and show how this success frequently comes about. What DOCI and pCCD lack are an effective treatment of dynamic correlations, which we here add by including higher-seniority cluster amplitudes which are excluded from pCCD. This frozen pair coupled cluster approach is comparable in cost to traditional closed-shell coupled cluster methods with results that are competitive for weakly correlated systems and often superior for the description of strongly correlated systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Density Matrix Embedding from Broken Symmetry Lattice Mean-Fields

Several variants of the recently proposed Density Matrix Embedding Theory (DMET) [G. Knizia and G. K-L. Chan, Phys. Rev. Lett. 109, 186404 (2012)] are formulated and tested. We show that spin symmetry breaking of the lattice mean-field allows precise control of the lattice and fragment filling while providing very good agreement between predicted properties and exact results. We present a rigorous proof that at convergence this method is guaranteed to preserve lattice and fragment filling. Differences arising from fitting the fragment one-particle density matrix alone versus fitting fragment plus bath are scrutinized. We argue that it is important to restrict the density matrix fitting to solely the fragment. Furthermore, in the proposed broken symmetry formalism, it is possible to substantially simplify the embedding procedure without sacrificing its accuracy by resorting to density instead of density matrix fitting. This simplified Density Embedding Theory (DET) greatly improves the convergence properties of the algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Particle-particle and quasiparticle random phase approximations: Connections to coupled cluster theory

We establish a formal connection between the particle-particle (pp) random phase approximation (RPA) and the ladder channel of the coupled cluster doubles (CCD) equations. The relationship between RPA and CCD is best understood within a Bogoliubov quasiparticle (qp) RPA formalism. This work is a follow-up to our previous formal proof on the connection between particle-hole (ph) RPA and ring-CCD. Whereas RPA is a quasibosonic approximation, CC theory is a correct bosonization in the sense that the wavefunction and Hilbert space are exactly fermionic. Coupled cluster theory achieves this goal by interacting the ph (ring) and pp (ladder) diagrams via a third channel that we here call "crossed-ring" whose presence allows for full fermionic antisymmetry. Additionally, coupled cluster incorporates what we call "mosaic" terms which can be absorbed into defining a new effective one-body Hamiltonian. The inclusion of these mosaic terms seems to be quite important. The pp-RPA an d qp-RPA equations are textbook material in nuclear structure physics but are largely unknown in quantum chemistry, where particle number fluctuations and Bogoliubov determinants are rarely used. We believe that the ideas and connections discussed in this paper may help design improved ways of incorporating RPA correlation into density functionals based on a CC perspective.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Proper and improper zero energy modes in Hartree-Fock theory and their relevance for symmetry breaking and restoration

We study the spectra of the molecular orbital Hessian (stability matrix) and random-phase approximation Hamiltonian of broken-symmetry Hartree-Fock solutions, focusing on zero eigenvalue modes. After all negative eigenvalues are removed from the Hessian by following their eigenvectors downhill, one is left with only positive and zero eigenvalues. Zero modes correspond to orbital rotations with no restoring force. These rotations determine states in the Goldstone manifold, which originates from a spontaneously broken continuous symmetry in the wave function. Zero modes can be classified as improper or proper according to their different mathematical and physical properties. Improper modes arise from symmetry breaking and their restoration always lowers the energy. Proper modes, on the other hand, correspond to degeneracies of the wave function, and their symmetry restoration does not necessarily lower the energy. We discuss how the RPA Hamiltonian distinguishes between proper and improper modes by doubling the number of zero eigenvalues associated with the latter. Proper modes in the Hessian always appear in pairs which do not double in RPA. We present several pedagogical cases exemplifying the above statements. The relevance of these results for projected Hartree-Fock methods is also addressed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Semilocal and Hybrid Meta-Generalized Gradient Approximations Based on the Understanding of the Kinetic-Energy-Density Dependence

We present a global hybrid meta-generalized gradient approximation (meta-GGA) with three empirical parameters, as well as its underlying semilocal meta-GGA and a meta-GGA with only one empirical parameter. All of them are based on the new meta-GGA resulting from the understanding of kinetic-energy-density dependence [J. Chem. Phys. 137, 051101 (2012)]. The obtained functionals show robust performances on the considered molecular systems for the properties of heats of formation, barrier heights, and noncovalent interactions. The pair-wise additive dispersion corrections to the functionals are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structural Phase Transitions of the Metal Oxide Perovskites SrTiO3, LaAlO3 and LaTiO3 Studied with a Screened Hybrid Functional

We have investigated the structural phase transitions of the transition metal oxide perovskites SrTiO$_{3}$, LaAlO$_{3}$ and LaTiO$_{3}$ using the screened hybrid density functional of Heyd, Scuseria and Ernzerhof (HSE06). We show that HSE06-computed lattice parameters, octahedral tilts and rotations, as well as electronic properties, are significantly improved over semilocal functionals. We predict the crystal field splitting ($Δ_{CF}$) resulting from the structural phase transition in SrTiO$_{3}$ and LaAlO$_{3}$ to be 3 meV and 10 meV, respectively, in excellent agreement with experimental results. HSE06 identifies correctly LaTiO$_{3}$ in the magnetic sates as a Mott insulator. Also, it predicts that the GdFeO$_{3}$-type distortion in non-magnetic LaTiO$_{3}$ will induce a large $Δ_{CF}$ of 410 meV. This large crystal-field splitting associated with the large magnetic moment found in the G-type antiferromagnetic state suggest that LaTiO$_{3}$ has an induced orbital order, which is confirmed by the visualisation of the highest occupied orbitals. These results strongly indicate that HSE06 is capable of efficiently and accurately modeling perovskite oxides, and promises to efficiently capture the physics at their heterointerfaces.

0 citations0 reviewsNo code linkedOpen access