We present in this contribution the first application of local correlation in the context of multicomponent methods. Multicomponent approaches allow for the targeted simulation of electrons together with other Fermions (most commonly protons) as quantum particles. These methods have become increasingly popular over the last years, particularly for the description of nuclear quantum effects (in strong hydrogen bonds, proton tunneling, and many more). However, most implementations are still based on canonical formulations of wave function theory, which we know for decades to be computationally inefficient for capturing dynamical correlation effects. Local correlation approaches, particularly with the use of pair natural orbitals (PNOs), enable asymptotically linear scaling of computational costs with very little impact on the overall accuracy. In this context, the efficient use of density fitting approximations in the integral calculation proves essential. We start by discussing our implementation of density-fitted NEO-MP2 and NEO-PNO-LMP2, upgrading the electronic correlation treatment up to PNO local coupled cluster level of theory. Several challenging examples are provided to benchmark the method in terms of accuracy as well as computational cost scaling. Following appropriate protocols, anharmonic corrections to localized X-H stretches can be applied routinely with little computational overhead.
Methods are described which can be used for an analytical description of the strong electron correlations in the planes of the Cu–O based high-Tc superconducting materials. One is a projection technique which is applied in order to calculate the photoelectron and inverse photoelectron spectra for the cases of half filling and hole doping of the Cu–O planes. Good agreement with the results of exact diagonalization of (CuO2)4 clusters is obtained, proving the usefulness of the method. Another one is the concept of string states which can be applied to study the t–J model. Also in this case a number of numerical findings can be reproduced and improved by this analytic method.
In a previous publication, we have discussed the usage of tensor decomposition in the canonical polyadic (CP) tensor format for electronic structure methods. There, we focused on two-electron integrals and second order Møller-Plesset perturbation theory (MP2). In this work, we discuss the CP format for Coupled Cluster (CC) theory and present a pilot implementation for the Coupled Cluster Doubles method. We discuss the iterative solution of the CC amplitude equations using tensors in CP representation and present a tensor contraction scheme that minimizes the effort necessary for the rank reductions during the iterations. Furthermore, several details concerning the reduction of complexity of the algorithm, convergence of the CC iterations, truncation errors, and the choice of threshold for chemical accuracy are discussed.
The influence of Zn and Ni impurities on the normal-state pseudogap of underdoped high-Tc cuprate superconductors is studied using exact diagonalization of effective $t\ensuremath{-}J\text{-like}$ Hamiltonians describing low energy electronic excitations of the CuO2 plane with some of the copper ions replaced with Zn/Ni. The Ni case Hamiltonian has been obtained by a sequence of approximations from a more complete model involving Cu 3d, Ni 3d, and O 2p orbitals. Our main findings are: (i) The width ${\mathrm{\ensuremath{\Omega}}}_{\mathrm{PG}}$ of the pseudogap occurring in the many body density of states, and manifesting itself also in the c-axis infrared conductivity, decreases with increasing Zn concentration as a consequence of a suppression of short range spin correlations. (ii) In the case of one hole and one Ni impurity, the hole is---for realistic values of the model parameters---weakly bound to the Ni site. This causes a slight increase of ${\mathrm{\ensuremath{\Omega}}}_{\mathrm{PG}}$ with respect to the pure case. (iii) Based on this result and further results for 1--2 holes and 1--2 Ni impurities, we suggest that in the real Ni substituted CuO2 plane ${\mathrm{\ensuremath{\Omega}}}_{\mathrm{PG}}$ is larger than in the pure case due to the binding of the doped holes to the Ni sites and effective underdoping. Our findings clarify the trends observed in the c-axis infrared conductivity data of Zn and Ni substituted (Sm,Nd)${\mathrm{Ba}}_{2}{\mathrm{Cu}}_{3}{\mathrm{O}}_{7\ensuremath{-}\ensuremath{\delta}}$ crystals.
Many properties of the two-dimensional Hubbard model have been explored for the model parameters appropriate for strongly correlated electronic systems, especially cuprate superconductors.Most of the calculations are done using a well-established dynamical cluster quantum Monte Carlo method.Using this method, we investigate the effect of long-range hoppings on superconductivity with and without the presence of phonons on small clusters.The superconducting transition temperature, T c , is found to generally decrease with a negative next nearest neighbor hopping, t ′ .In the presence of the Holstein phonons, a finite t ′ enhances T c in the under-doped region for the hole-doped system, consistent with band structure calculations and experiment.The validity of the spin-susceptibility-mediated pairing in this model is studied and found to yield symmetries other than d-wave when a finite t ′ is considered.A new numerical algorithm for solving the embedded cluster problems is introduced, and used to calculate the thermodynamic properties of the model on larger clusters, especially those associated with the quantum critical behavior at finite doping.Our results suggest that the quantum critical point (QCP) which separates the Fermi liquid and pseudogap regions, is the second order terminus of the line of first order phase separation transition in the limit when a positive t ′ goes to zero.For small t ′ , the superconducting dome is centered at the QCP, suggesting that charge fluctuations might have a role in the pairing mechanism in this model.I am very thankful to University of Cincinnati for providing me with the opportunity to pursue my education and the support I received during my program.Part of this thesis is a result of my work as a researcher at Louisiana State University.I am also thankful to them for their support.I would like to thank my adviser, Mark Jarrell, who taught me how to be an effective researcher, provided me with the opportunities to interact with other professionals and fully supported me financially for most of the duration of the
We have studied electron correlations in the doped two-dimensional (2D) Hubbard model by using the coupled-cluster method (CCM) to investigate whether or not the method can be applied to correct the independent particle approximations actually used in ab initio band calculations. The double excitation version of the CCM, implemented using the approximate coupled pair (ACP) method, accounts for most of the correlation energies of the 2D Hubbard model in the weak $(U/t\ensuremath{\simeq}1)$ and the intermediate U/t regions $(U/t\ensuremath{\simeq}4).$ The error in this case is always less than 1 The ACP approximation gets less accurate for large U/t $(U/t\ensuremath{\simeq}8)$ and near half-filling. Further incorporation of electron correlation effects is necessary in this region. The accuracy depends neither on the system size, nor the gap between the lowest unoccupied level and the highest occupied level in the finite size effect. Hence, the CCM may be favorably applied to ab initio band calculations on metals as well as semiconductors and insulators.