Source author record

Carlos A. Jiménez-Hoyos

Carlos A. Jiménez-Hoyos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

16 published item(s)

preprint2022arXiv

Coupled cluster and perturbation theories based on a cluster mean-field reference applied to strongly correlated spin systems

We introduce perturbation and coupled-cluster theories based on a cluster mean-field reference for describing the ground state of strongly-correlated spin systems. In cluster mean-field, the ground state wavefunction is written as a simple tensor product of optimized cluster states. The cluster language and the mean-field nature of the ansatz allows for a straightforward improvement which uses perturbation theory and coupled-cluster to account for inter-cluster correlations. We present benchmark calculations on the 1D chain and 2D square $J_1-J_2$ Heisenberg model, using cluster mean-field, perturbation theory and coupled-cluster. We also present an extrapolation scheme that allows us to compute thermodynamic limit energies accurately. Our results indicate that, with sufficiently large clusters, the correlated methods (cPT2, cPT4 and cCCSD) can provide a relatively accurate description of the Heisenberg model in the regimes considered, which suggests that the methods presented can be used for other strongly-correlated systems. Some ways to improve upon the methods presented in this work are discussed.

preprint2021arXiv

Ground States of Heisenberg Spin Clusters from Projected Hartree-Fock Theory

We apply Projected Hartree-Fock theory (PHF) for approximating ground states of Heisenberg spin clusters. Spin-rotational, point-group and complex-conjugation symmetry are variationally restored from a broken-symmetry mean-field reference, where the latter corresponds to a product of local spin states. A fermionic formulation of the Heisenberg model furnishes a conceptual connection to PHF applications in quantum chemistry and detailed equations for a self-consistent field optimization of the reference state are provided. Different PHF variants are benchmarked for ground-state energies and spin-pair correlation functions of antiferromagnetic spin rings and three different polyhedra, with various values of the local spin-quantum number s. The low computational cost and the compact wave-function representation make PHF a promising complement to existing approaches for ground states of molecular spin clusters, particularly for large s and moderately large N. The present work may also motivate future explorations of more accurate post-PHF methods for Heisenberg spin clusters.

preprint2020arXiv

Approaching the Full Configuration Interaction Ground State from an Arbitrary Wavefunction with Gradient Descent and Quasi-Newton Algorithms

We consider gradient descent and quasi-Newton algorithms to optimize the full configuration interaction (FCI) ground state wavefunction starting from an arbitrary reference state $|0 \rangle$. We show that the energies obtained along the optimization path can be evaluated in terms of expectation values of $|0 \rangle$, thus avoiding explicit storage of intermediate wavefunctions. This allows us to find the energies after the first few steps of the FCI algorithm for systems much larger than what standard, deterministic FCI codes can handle at present. We show an application of the algorithm with reference wavefunctions constructed as linear combinations of non-orthogonal determinants.

preprint2020arXiv

Minimal matrix product states and generalizations of mean-field and geminal wavefunctions

Simple wavefunctions of low computational cost but which can achieve qualitative accuracy across the whole potential energy surface (PES) are of relevance to many areas of electronic structure as well as to applications to dynamics. Here, we explore a class of simple wavefunctions, the minimal matrix product state (MMPS), that generalizes many simple wavefunctions in common use, such as projected mean-field wavefunctions, geminal wavefunctions, and generalized valence bond states. By examining the performance of MMPSs for PESs of some prototypical systems, we find that they yield good qualitative behavior across the whole PES, often significantly improving on the aforementioned ansätze.

preprint2020arXiv

On a dual representation of the Goldstone manifold

An intrinsic wavefunction with a broken continuous symmetry can be rotated with no energy penalty leading to an infinite set of degenerate states known as a Goldstone manifold. In this work, we show that a dual representation of such manifold exists that is sampled by an infinite set of non-degenerate states. A proof that both representations are equivalent is provided. From the work of Peierls and Yoccoz (Proc. Phys. Soc. A {\bf 70}, 381 (1957)), it is known that collective states with good symmetries can be obtained from the Goldstone manifold using a generator coordinate trial wavefunction. We show that an analogous generator coordinate can be used in the dual representation; we provide numerical evidence using an intrinsic wavefunction with particle number symmetry-breaking for the electronic structure of the Be atom and one with $\hat{S}^z$ symmetry-breaking for a H$_5$ ring. We discuss how the dual representation can be used to evaluate expectation values of symmetry-projected states when the norm $|\langle Φ| \hat{P}^q | Φ\rangle|$ becomes very small.

preprint2016arXiv

A practical guide to density matrix embedding theory in quantum chemistry

Density matrix embedding theory (DMET) provides a theoretical framework to treat finite fragments in the presence of a surrounding molecular or bulk environment, even when there is significant correlation or entanglement between the two. In this work, we give a practically oriented and explicit description of the numerical and theoretical formulation of DMET. We also describe in detail how to perform self-consistent DMET optimizations. We explore different embedding strategies with and without a self-consistency condition in hydrogen rings, beryllium rings, and a sample S$_{\text{N}}$2 reaction. The source code for the calculations in this work can be obtained from \url{https://github.com/sebwouters/qc-dmet}.

preprint2015arXiv

A cluster-based mean-field and perturbative description of strongly correlated fermion systems. Application to the 1D and 2D Hubbard model

We introduce a mean-field and perturbative approach, based on clusters, to describe the ground state of fermionic strongly-correlated systems. In cluster mean-field, the ground state wavefunction is written as a simple tensor product over optimized cluster states. The optimization of the single-particle basis where the cluster mean-field is expressed is crucial in order to obtain high-quality results. The mean-field nature of the ansatz allows us to formulate a perturbative approach to account for inter-cluster correlations; other traditional many-body strategies can be easily devised in terms of the cluster states. We present benchmark calculations on the half-filled 1D and (square) 2D Hubbard model, as well as the lightly-doped regime in 2D, using cluster mean-field and second-order perturbation theory. Our results indicate that, with sufficiently large clusters or to second-order in perturbation theory, a cluster-based approach can provide an accurate description of the Hubbard model in the considered regimes. Several avenues to improve upon the results presented in this work are discussed.

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.

preprint2014arXiv

Multi-reference symmetry-projected variational approximation for the ground state of the doped one-dimensional Hubbard model

A multi-reference configuration mixing scheme is used to describe the ground state, characterized by well defined spin and space group symmetry quantum numbers as well as doping fractions $N_{e}/N_{sites}$, of one dimensional Hubbard lattices with nearest-neighbor hopping and periodic boundary conditions. Within this scheme, each ground state is expanded in a given number of nonorthogonal and variationally determined symmetry-projected configurations. The results obtained for the ground state and correlation energies of half-filled and doped lattices with 30, 34 and 50 sites, compare well with the exact Lieb-Wu solutions as well as with the ones obtained with other state-of-the-art approximations. The structure of the intrinsic symmetry-broken determinants resulting from the variational procedure is interpreted in terms of solitons whose translational and breathing motions can be regarded as basic units of quantum fluctuations. It is also shown that in the case of doped 1D lattices, a part of such fluctuations can also be interpreted in terms of polarons. In addition to momentum distributions, both spin-spin and density-density correlation functions are studied as functions of doping. The spectral functions and density of states, computed with an ansatz whose quality can be well-controlled by the number of symmetry-projected configurations used to approximate the $N_{e} \pm 1$ electron systems, display features beyond a simple quasiparticle distribution, as well as spin-charge separation trends.

preprint2014arXiv

Polyradical character and spin frustration in fullerene molecules: An ab initio non-collinear Hartree--Fock study

Most {\em ab initio} calculations on fullerene molecules have been carried out based on the paradigm of the Hückel model. This is consistent with the restricted nature of the independent-particle model underlying such calculations, even in single-reference-based correlated approaches. On the other hand, previous works on some of these molecules using model Hamiltonians have clearly indicated the importance of short-range inter-atomic spin-spin correlations. In this work, we consider {\em ab initio} non-collinear Hartree--Fock (HF) solutions for representative fullerene systems: the bowl, cage, ring, and pentagon isomers of C$_{20}$, and the larger C$_{30}$, C$_{36}$, C$_{60}$, C$_{70}$, and C$_{84}$ fullerene cages. In all cases but the ring we find that the HF minimum corresponds to a truly non-collinear solution with a torsional spin density wave. Optimized geometries at the generalized HF (GHF) level lead to fully symmetric structures, even in those cases where Jahn-Teller distortions have been previously considered. The nature of the GHF solutions is consistent with the $π$-electron space becoming polyradical in nature: each $p$-orbital remains effectively singly occupied. The spin frustration, induced by the pentagon rings in an otherwise anti-ferromagnetic background, is minimized at the HF level by aligning the spins in non-collinear arrangements. The long-range magnetic ordering observed is reminiscent of the character of broken symmetry HF solutions in polyacene systems.

preprint2014arXiv

Symmetry-projected Wave Functions in Quantum Monte Carlo Calculations

We consider symmetry-projected Hartree--Fock trial wave functions in constrained-path Monte Carlo (CPMC) calculations. Previous CPMC calculations have mostly employed Hartree--Fock (HF) trial wave functions, restricted or unrestricted. The symmetry-projected HF approach results in a hierarchy of wave functions with increasing quality: the more symmetries that are broken and restored in a self-consistent manner, the higher the quality of the trial wave function. This hierarchy is approximately maintained in CPMC calculations: the accuracy in the energy increases and the statistical variance decreases when further symmetries are broken and restored. Significant improvement is achieved in CPMC with the best symmetry-projected trial wave functions over those from simple HF. We analyze and quantify the behavior using the two-dimensional repulsive Hubbard model as an example. In the sign-problem-free region, where CPMC can be made exact but a constraint is deliberately imposed here, spin-projected wave functions remove the constraint bias. Away from half-filling, spatial symmetry restoration in addition to that of the spin lead to highly accurate results from CPMC. Since the computational cost of symmetry-projected HF trial wave functions in CPMC can be made to scale algebraically with system size, this provides a potentially general approach for accurate calculations in many-fermion systems.

preprint2014arXiv

Variational description of the ground state of the repulsive two-dimensional Hubbard model in terms of nonorthogonal symmetry-projected Slater determinants

The few determinant (FED) methodology, introduced in our previous works for 1D lattices, is here adapted for the repulsive two-dimensional Hubbard model at half-filling and with finite doping fractions. Within this configuration mixing scheme, a given ground state with well defined spin and space group quantum numbers, is expanded in terms of a nonorthogonal symmetry-projected basis determined through chains of variation-after projection calculations. The results obtained for the ground state and correlation energies of half-filled and doped 4$\times$4, 6$\times$6, 8$\times$8, and 10$\times$10 lattices, as well as momentum distributions and spin-spin correlation functions in small lattices, compare well with those obtained using other state-of-the-art approximations. The structure of the intrinsic determinants resulting from the variational strategy is interpreted in terms of defects that encode information on the basic units of quantum fluctuations in the considered 2D systems. The varying nature of the underlying quantum fluctuations, reflected in a transition to a stripe regime for increasing onsite repulsions, is discussed using the intrinsic determinants belonging to a 16$\times$4 lattice with 56 electrons. Such a transition is further illustrated by computing spin-spin and charge-charge correlation functions with the corresponding multireference FED wave functions. In good agreement with previous studies, the analysis of the pairing correlation functions reveals a weak enhancement of the extended $s$-wave and d$_{x^2-y^2}$ pairing modes. Given the quality of results here reported together with those previously obtained for 1D lattices and the parallelization properties of the FED scheme, we believe that symmetry projection techniques are very well suited for building ground state wave functions of correlated electronic systems, regardless of their dimensionality.

preprint2013arXiv

Excited electronic states from a variational approach based on symmetry-projected Hartree--Fock configurations

Recent work from our research group has demonstrated that symmetry-projected Hartree--Fock (HF) methods provide a compact representation of molecular ground state wavefunctions based on a superposition of non-orthogonal Slater determinants. The symmetry-projected ansatz can account for static correlations in a computationally efficient way. Here we present a variational extension of this methodology applicable to excited states of the same symmetry as the ground state. Benchmark calculations on the C$_2$ dimer with a modest basis set, which allows comparison with full configuration interaction results, indicate that this extension provides a high quality description of the low-lying spectrum for the entire dissociation profile. We apply the same methodology to obtain the full low-lying vertical excitation spectrum of formaldehyde, in good agreement with available theoretical and experimental data, as well as to a challenging model $C_{2v}$ insertion pathway for BeH$_2$. The variational excited state methodology developed in this work has two remarkable traits: it is fully black-box and will be applicable to fairly large systems thanks to its mean-field computational cost.

preprint2013arXiv

Multi-component symmetry-projected approach for molecular ground state correlations

The symmetry-projected Hartree--Fock ansatz for the electronic structure problem can efficiently account for static correlation in molecules, yet it is often unable to describe dynamic correlation in a balanced manner. Here, we consider a multi-component, systematically-improvable approach, that accounts for all ground state correlations. Our approach is based on linear combinations of symmetry-projected configurations built out of a set of non-orthogonal, variationally optimized determinants. The resulting wavefunction preserves the symmetries of the original Hamiltonian even though it is written as a superposition of deformed (broken-symmetry) determinants. We show how short expansions of this kind can provide a very accurate description of the electronic structure of simple chemical systems such as the nitrogen and the water molecules, along the entire dissociation profile. In addition, we apply this multi-component symmetry-projected approach to provide an accurate interconversion profile among the peroxo and bis($μ$-oxo) forms of [Cu$_2$O$_2$]$^{2+}$, comparable to other state-of-the-art quantum chemical methods.

preprint2013arXiv

Multi-reference symmetry-projected variational approaches for ground and excited states of the one-dimensional Hubbard model

We present a multi-reference configuration mixing scheme for describing ground and excited states, with well defined spin and space group symmetry quantum numbers, of the one-dimensional Hubbard model with nearest-neighbor hopping and periodic boundary conditions. Within this scheme, each state is expanded in terms of non-orthogonal and variationally determined symmetry-projected configurations. The results for lattices up to 30 and 50 sites compare well with the exact Lieb-Wu solutions as well as with results from other state-of-the-art approximations. In addition to spin-spin correlation functions in real space and magnetic structure factors, we present results for spectral functions and density of states computed with an ansatz whose quality can be well-controlled by the number of symmetry-projected configurations used to approximate the systems with $N_{e}$ and $N_{e} \pm 1$ electrons. The intrinsic symmetry-broken determinants resulting from the variational calculations have rich structures in terms of defects that can be regarded as basic units of quantum fluctuations. Given the quality of the results here reported, as well as the parallelization properties of the considered scheme, we believe that symmetry-projection techniques, which have found ample applications in nuclear structure physics, deserve further attention in the study of low-dimensional correlated many-electron systems.

preprint2012arXiv

Symmetry-projected variational approach for ground and excited states of the two-dimensional Hubbard model

We present a symmetry-projected configuration mixing scheme to describe ground and excited states, with well defined quantum numbers, of the two-dimensional Hubbard model with nearestneighbor hopping and periodic boundary conditions. Results for the half-filled 2{\times}4, 4{\times}4, and 6{\times}6 lattices, as well as doped 4 {\times} 4 systems, compare well with available results, both exact and from other state-of-the-art approximations. We report spectral functions and density of states obtained from a well-controlled ansatz for the (Ne {\pm} 1)-electron system. Symmetry projected methods have been widely used for the many-body nuclear physics problem but have received little attention in the solid state community. Given their relatively low (mean-field) computational cost and the high quality of results here reported, we believe that they deserve further scrutiny.