Quasiparticle self-consistent GW study of cuprates: electronic structure, model parameters, and the two-band theory for Tc.
Ontology highlight
ABSTRACT:
SUBMITTER: Jang SW
PROVIDER: S-EPMC4513282 | biostudies-other | 2015
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Low-order scaling GW implementations for molecules are usually restricted to approximations with diagonal self-energy. Here, we present an all-electron implementation of quasiparticle self-consistent GW for molecular systems. We use an efficient algorithm for the evaluation of the self-energy in imaginary time, from which a static non-local exchange-correlation potential is calculated via analytical continuation. By using a direct inversion of iterative subspace method, fast and stable convergence is achieved for almost all molecules in the GW100 database. Exceptions are systems which are associated with a breakdown of the single quasiparticle picture in the valence region. The implementation is proven to be starting point independent and good agreement of QP energies with other codes is observed. We demonstrate the computational efficiency of the new implementation by calculating the quasiparticle spectrum of a DNA oligomer with 1,220 electrons using a basis of 6,300 atomic orbitals in less than 4 days on a single compute node with 16 cores. We use then our implementation to study the dependence of quasiparticle energies of DNA oligomers consisting of adenine-thymine pairs on the oligomer size. The first ionization potential in vacuum decreases by nearly 1 electron volt and the electron affinity increases by 0.4 eV going from the smallest to the largest considered oligomer. This shows that the DNA environment stabilizes the hole/electron resulting from photoexcitation/photoattachment. Upon inclusion of the aqueous environment via a polarizable continuum model, the differences between the ionization potentials reduce to 130 meV, demonstrating that the solvent effectively compensates for the stabilizing effect of the DNA environment. The electron affinities of the different oligomers are almost identical in the aqueous environment.
| S-EPMC8446457 | biostudies-literature
Project description:We report an all-electron, atomic orbital (AO)-based, two-component (2C) implementation of the GW approximation (GWA) for closed-shell molecules. Our algorithm is based on the space-time formulation of the GWA and uses analytical continuation (AC) of the self-energy, and pair-atomic density fitting (PADF) to switch between AO and auxiliary basis. By calculating the dynamical contribution to the GW self-energy at a quasi-one-component level, our 2C-GW algorithm is only about a factor of 2-3 slower than in the scalar relativistic case. Additionally, we present a 2C implementation of the simplest vertex correction to the self-energy, the statically screened G3W2 correction. Comparison of first ionization potentials (IPs) of a set of 67 molecules with heavy elements (a subset of the SOC81 set) calculated with our implementation against results from the WEST code reveals mean absolute deviations (MAD) of around 70 meV for G0W0@PBE and G0W0@PBE0. We check the accuracy of our AC treatment by comparison to full-frequency GW calculations, which shows that in the absence of multisolution cases, the errors due to AC are only minor. This implies that the main sources of the observed deviations between both implementations are the different single-particle bases and the pseudopotential approximation in the WEST code. Finally, we assess the performance of some (partially self-consistent) variants of the GWA for the calculation of first IPs by comparison to vertical experimental reference values. G0W0@PBE0 (25% exact exchange) and G0W0@BHLYP (50% exact exchange) perform best with mean absolute deviations (MAD) of about 200 meV. Explicit treatment of spin-orbit effects at the 2C level is crucial for systematic agreement with experiment. On the other hand, eigenvalue-only self-consistent GW (evGW) and quasi-particle self-consistent GW (qsGW) significantly overestimate the IPs. Perturbative G3W2 corrections increase the IPs and therefore improve the agreement with experiment in cases where G0W0 alone underestimates the IPs. With a MAD of only 140 meV, 2C-G0W0@PBE0 + G3W2 is in best agreement with the experimental reference values.
| S-EPMC10501001 | biostudies-literature
Project description:Choosing optimal representation methods of atomic and electronic structures is essential when machine learning properties of materials. We address the problem of representing quantum states of electrons in a solid for the purpose of machine leaning state-specific electronic properties. Specifically, we construct a fingerprint based on energy decomposed operator matrix elements (ENDOME) and radially decomposed projected density of states (RAD-PDOS), which are both obtainable from a standard density functional theory (DFT) calculation. Using such fingerprints we train a gradient boosting model on a set of 46k G0W0 quasiparticle energies. The resulting model predicts the self-energy correction of states in materials not seen by the model with a mean absolute error of 0.14 eV. By including the material's calculated dielectric constant in the fingerprint the error can be further reduced by 30%, which we find is due to an enhanced ability to learn the correlation/screening part of the self-energy. Our work paves the way for accurate estimates of quasiparticle band structures at the cost of a standard DFT calculation.
| S-EPMC8813923 | biostudies-literature
Project description:An algorithm to perform stochastic generalized active space calculations, Stochastic-GAS, is presented, that uses the Slater determinant based FCIQMC algorithm as configuration interaction eigensolver. Stochastic-GAS allows the construction and stochastic optimization of preselected truncated configuration interaction wave functions, either to reduce the computational costs of large active space wave function optimizations, or to probe the role of specific electron correlation pathways. As for the conventional GAS procedure, the preselection of the truncated wave function is based on the selection of multiple active subspaces while imposing restrictions on the interspace excitations. Both local and cumulative minimum and maximum occupation number constraints are supported by Stochastic-GAS. The occupation number constraints are efficiently encoded in precomputed probability distributions, using the precomputed heat bath algorithm, which removes nearly all runtime overhead of GAS. This strategy effectively allows the FCIQMC dynamics to a priori exclude electronic configurations that are not allowed by GAS restrictions. Stochastic-GAS reduced density matrices are stochastically sampled, allowing orbital relaxations via Stochastic-GASSCF, and direct evaluation of properties that can be extracted from density matrices, such as the spin expectation value. Three test case applications have been chosen to demonstrate the flexibility of Stochastic-GAS: (a) the Stochastic-GASSCF [5·(6, 6)] optimization of a stack of five benzene molecules, that shows the applicability of Stochastic-GAS toward fragment-based chemical systems; (b) an uncontracted stochastic MRCISD calculation that correlates 96 electrons and 159 molecular orbitals, and uses a large (32, 34) active space reference wave function for an Fe(II)-porphyrin model system, showing how GAS can be applied to systematically recover dynamic electron correlation, and how in the specific case of the Fe(II)-porphyrin dynamic correlation further differentially stabilizes the 3Eg over the 5A1g spin state; (c) the study of an Fe4S4 cluster's spin-ladder energetics via highly truncated stochastic-GAS [4·(5, 5)] wave functions, where we show how GAS can be applied to understand the competing spin-exchange and charge-transfer correlating mechanisms in stabilizing different spin-states.
| S-EPMC8757470 | biostudies-literature
Project description:Individual regulatory proteins are typically charged with the simultaneous regulation of a battery of different genes. As a result, when one of these proteins is limiting, competitive effects have a significant impact on the transcriptional response of the regulated genes. Here we present a general framework for the analysis of any generic regulatory architecture that accounts for the competitive effects of the regulatory environment by isolating these effects into an effective concentration parameter. These predictions are formulated using the grand-canonical ensemble of statistical mechanics and the fold-change in gene expression is predicted as a function of the number of transcription factors, the strength of interactions between the transcription factors and their DNA binding sites, and the effective concentration of the transcription factor. The effective concentration is set by the transcription factor interactions with competing binding sites within the cell and is determined self-consistently. Using this approach, we analyze regulatory architectures in the grand-canonical ensemble ranging from simple repression and simple activation to scenarios that include repression mediated by DNA looping of distal regulatory sites. It is demonstrated that all the canonical expressions previously derived in the case of an isolated, non-competing gene, can be generalised by a simple substitution to their grand canonical counterpart, which allows for simple intuitive incorporation of the influence of multiple competing transcription factor binding sites. As an example of the strength of this approach, we build on these results to present an analytical description of transcriptional regulation of the lac operon.
| S-EPMC5501422 | biostudies-literature
Project description:A major challenge in understanding the cuprate superconductors is to clarify the nature of the fundamental electronic correlations that lead to the pseudogap phenomenon. Here we use ultrashort light pulses to prepare a non-thermal distribution of excitations and capture novel properties that are hidden at equilibrium. Using a broadband (0.5-2 eV) probe, we are able to track the dynamics of the dielectric function and unveil an anomalous decrease in the scattering rate of the charge carriers in a pseudogap-like region of the temperature (T) and hole-doping (p) phase diagram. In this region, delimited by a well-defined T*neq(p) line, the photoexcitation process triggers the evolution of antinodal excitations from gapped (localized) to delocalized quasiparticles characterized by a longer lifetime. The novel concept of photo-enhanced antinodal conductivity is naturally explained within the single-band Hubbard model, in which the short-range Coulomb repulsion leads to a k-space differentiation between nodal quasiparticles and antinodal excitations.
| S-EPMC4104437 | biostudies-literature
Project description:For the computational prediction of core electron binding energies in solids, two distinct kinds of modeling strategies have been pursued: the Δ-Self-Consistent-Field method based on density functional theory (DFT), and the GW method. In this study, we examine the formal relationship between these two approaches and establish a link between them. The link arises from the equivalence, in DFT, between the total energy difference result for the first ionization energy, and the eigenvalue of the highest occupied state, in the limit of infinite supercell size. This link allows us to introduce a new formalism, which highlights how in DFT─even if the total energy difference method is used to calculate core electron binding energies─the accuracy of the results still implicitly depends on the accuracy of the eigenvalue at the valence band maximum in insulators, or at the Fermi level in metals. We examine whether incorporating a quasiparticle correction for this eigenvalue from GW theory improves the accuracy of the calculated core electron binding energies, and find that the inclusion of vertex corrections is required for achieving quantitative agreement with experiment.
| S-EPMC10269324 | biostudies-literature
Project description:Organic-molecule fluorophores with emission wavelengths in the second near-infrared window (NIR-II, 1000-1700 nm) have attracted substantial attention in the life sciences and in biomedical applications because of their excellent resolution and sensitivity. However, adequate theoretical levels to provide efficient and accurate estimations of the optical and electronic properties of organic NIR-II fluorophores are lacking. The standard approach for these calculations has been time-dependent density functional theory (TDDFT). However, the size and large excitonic energies of these compounds pose challenges with respect to computational cost and time. In this study, we used the GW approximation combined with the Bethe-Salpeter equation (GW-BSE) implemented in many-body perturbation theory approaches based on density functional theory. This method was used to perform calculations of the excited states of two NIR molecular fluorophores (BTC980 and BTC1070), going beyond TDDFT. In this study, the optical absorption spectra and frontier molecular orbitals of these compounds were compared using TDDFT and GW-BSE calculations. The GW-BSE estimates showed excellent agreement with previously reported experimental results.
| S-EPMC8466807 | biostudies-literature
Project description:The G0W0, evGW0, evGW, and scGW0 approximations to many-body perturbation theory combined with the Bethe-Salpeter approach (BSE) are applied to calculate electronic and optical properties of the open-shell spinel ferrite ZnFe2 O4 . The effect of the various degrees of self-consistency is assessed by comparison to recent experimental results. Furthermore, the influence of the method for obtaining the ground-state wavefunction is studied, including the GGA functional PBE with and without an intermediate step using the COHSEX approximation, as well as PBE+U, where we try to minimize the influence of the Hubbard potential U. Best agreement for the optical band gap and the first maxima of the excitation spectrum is obtained with the evGW method based on a PBE+U wavefunction. This method is chosen and converged carefully to yield quantitative results for the optical spectra of four different magnetic structures and cation distributions of ZnFe2 O4 . With the results we provide a possible explanation for inconsistency in experimental results.
| S-EPMC7155046 | biostudies-literature
Project description:The notions of Bloch wave, crystal momentum, and energy bands are commonly regarded as unique features of crystalline materials with commutative translation symmetries. Motivated by the recent realization of hyperbolic lattices in circuit quantum electrodynamics, we exploit ideas from algebraic geometry to construct a hyperbolic generalization of Bloch theory, despite the absence of commutative translation symmetries. For a quantum particle propagating in a hyperbolic lattice potential, we construct a continuous family of eigenstates that acquire Bloch-like phase factors under a discrete but noncommutative group of hyperbolic translations, the Fuchsian group of the lattice. A hyperbolic analog of crystal momentum arises as the set of Aharonov-Bohm phases threading the cycles of a higher-genus Riemann surface associated with this group. This crystal momentum lives in a higher-dimensional Brillouin zone torus, the Jacobian of the Riemann surface, over which a discrete set of continuous energy bands can be computed.
| S-EPMC8442860 | biostudies-literature