Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical review. B, Condensed matter

Exciton binding energy in the strong correlation limit of conjugated chains

Zhigang Shuai, Jean‐Luc Brédas, Swapan Kumar Pati, S. Ramasesha

Published Dec 15, 1998• 23 citations
Publisher

Abstract

By applying the numerically accurate symmetrized density-matrix renormalization-group method coupled with the extended Hubbard-Peierls model, we find that (i) the on-site Hubbard repulsion energy U dramatically reduces the binding energy of the lowest optically allowed ${1B}_{u}$ exciton; (ii) in the zero-dimerization limit, there exists a critical value of V at which the ${1B}_{u}$ exciton becomes bound; the critical value ${V}_{c}=2t$ is fully in agreement with the recent analytical results at the infinite-U limit by Gallagher and Mazumdar [Phys. Rev. B 56, 15 025 (1997)], furthermore, this critical value decreases appreciably for weaker on-site correlation strengths, when the dimerization amplitude $(\ensuremath{\delta})$ is nonzero. The present accurate numerical results contradict those obtained recently by Yu, Saxena, and Bishop [Phys. Rev. B 56, 3697 (1997)] both qualitatively and quantitatively. We also present first-order perturbation plus random-phase-approximation and single configuration-interaction analyses to rationalize the numerical calculations.

DOI
10.1103/physrevb.58.15329

Related papers

Binding energy of 1Bu singlet excitons in the one-dimensional extended Hubbard-Peierls model

Zhigang Shuai, Swapan Kumar Pati, W. P. Su +2 more

Using a symmetrized density-matrix renormalization-group formulation, we have investigated the electronic binding energy of the lowest optically allowed exciton (${1\mathrm{B}}_{\mathrm{u}}$) within an extended Hubbard-Peierls model (with parameters U, V, and \ensuremath{\delta}) for conjugated chains with N=80 sites. Three symmetries, the ${\mathrm{C}}_{2}$, spin parity, and electron-hole symmetries, have been applied to construct the projector operator. Analysis of our results on the ratios between exciton binding energies and exciton energies, sheds light on the current experimental and theoretical controversy on exciton binding for polyacetylene and polyparaphenylene vinylene. We show that in the absence of dimerization, the exciton binding energy is vanishingly small for U/t up to 5 and V/t up to 2.3; for a finite dimerization, it is only when V is large enough that the exciton gets to be significantly bound. This result questions the applicability to conjugated polymers of the strong correlation picture, in which the exciton binding energy is equal to V.

Physical review. B, Condensed matter•Jun 15, 1997•49 citations

Excitons and optical absorption in one-dimensional extended Hubbard models with short- and long-range interactions

Frank B. Gallagher, S. Mazumdar

We show that charge excitations within one-dimensional extended Hubbard models in the limit of large on-site Coulomb interaction can be mapped onto pairs of spin reversals from the ferromagnetic Ising-Heisenberg spin Hamiltonian. The exciton in the extended Hubbard models is equivalent to the bound two-magnon state of the spin Hamiltonian. From exact results for the spin Hamiltonian, for the intersite Coulomb interactions limited between nearest neighbors, the lowest optical exciton is separated in energy from the electron-hole continuum only for the intersite interaction larger than half the one-electron bandwidth. The commutation relationship between spin operators allows us to numerically investigate exciton binding energies and optical absorptions from the ground state in the limit of large on-site Coulomb interaction even for cases where analytic results are not available. For intersite Coulomb interaction limited between nearest neighbors, the binding energy of the lowest optical exciton changes very weakly for realistic Peierls bond dimerization, although there is a profound change in the optical absorption. While in the undimerized chain a single optical exciton can occur, two additional excitons that occur above the threshold of the electron-hole continuum become optically visible in the dimerized case. For long-range intersite Coulomb interactions, multiple excitons can occur even for the undimerized chain, and additional exciton states become optically visible upon dimerization. With the $1/r$ form for the intersite Coulomb interaction, our numerical calculations suggest that the formation of at least one optical exciton is unconditional, unlike the case of nearest-neighbor interaction. The binding energy of the lowest optical exciton again depends weakly on the dimerization for long-range interaction. Our theoretical results can explain the occurrence of two charge-transfer absorption bands in half-filled band segregated stack charge transfer solids below the spin-Peierls dimerization temperature. Our results also suggest that the lowest optical state of \ensuremath{\pi}-conjugated polymers, in which Coulomb interactions are long ranged, is necessarily an exciton.

Physical review. B, Condensed matter•Dec 15, 1997•68 citations

Density matrix renormalization group study of excitons in polydiacetylene chains

Gergely Barcza, Örs Legeza, Florian Gebhard +1 more

We study the elementary excitations of a model Hamiltonian for the $\ensuremath{\pi}$ electrons in polydiacetylene chains. In these materials, the bare band gap is only half the size of the observed single-particle gap and the binding energy of the exciton of 0.5 eV amounts to 20% of the single-particle gap. Therefore, exchange and correlations due to the long-range Coulomb interaction require a numerically exact treatment which we carry out using the density-matrix renormalization group method. Employing both the Hubbard-Ohno potential and the screened potential in one dimension, we reproduce the experimental results for the binding energy of the singlet exciton and its polarizability. Our results indicate that there are optically dark states below the singlet exciton, in agreement with experiment. In addition, we find a weakly bound second exciton with a binding energy of 0.1 eV. The energies in the triplet sector do not match the experimental data quantitatively, probably because we do not include polaronic relaxation effects.

arXiv:0908.4160•Physical Review B•Jan 5, 2010•14 citations

Perturbation theory for optical excitations in the one-dimensional extended Peierls–Hubbard model

Anja Grage, Florian Gebhard, Jörg Rissler

For the one-dimensional, extended Peierls–Hubbard model we calculate analytically the ground-state energy and the single-particle gap to second order in the Coulomb interaction for a given lattice dimerization. The comparison with numerically exact data from the density-matrix renormalization group shows that the ground-state energy is quantitatively reliable for Coulomb parameters as large as the bandwidth. The single-particle gap can almost triple from its bare Peierls value before substantial deviations appear. For the calculation of the dominant optical excitations, we follow two approaches. In Wannier theory, we perturb the Wannier exciton states to second order. In two-step perturbation theory, similar in spirit to the GW -BSE approach (where the Bethe–Salpeter equation is solved for quasi-particles which are calculated in the GW -approximation), we form excitons from dressed electron–hole excitations. We find the Wannier approach to be superior to the two-step perturbation theory. For singlet excitons, Wannier theory is applicable up to Coulomb parameters as large as half the bandwidth. For triplet excitons, second-order perturbation theory quickly fails completely.

arXiv:cond-mat/0505410•Journal of Statistical Mechanics Theory and Experiment•Aug 17, 2005•9 citations

Dimerization of polyacetylene using a distance-dependent Hubbard model and the density-matrix renormalization-group method

Marie-Bernadette Lepetit, G. M. Pastor

The dimerization of polyacetylene is studied using a distance-dependent Hubbard Hamiltonian for the valence $\ensuremath{\pi}$ electrons and approximating $\ensuremath{\sigma}$-electron contributions by a sum of pair interactions ${E}_{\ensuremath{\sigma}}{(r}_{\mathrm{ij}})$. The $\ensuremath{\pi}$-electron Coulomb repulsion $U$ is taken to be independent of the C-C distances ${r}_{\mathrm{ij}}$ and is varied systematically from the uncorrelated or H\"uckel limit all over to the strongly correlated or Heisenberg regime. For each value of $U$, ${E}_{\ensuremath{\sigma}}{(r}_{\mathrm{ij}})$, and the $\ensuremath{\pi}$-electron hopping integral ${t(r}_{\mathrm{ij}})$ are fitted to reproduce ab initio results on the ground-state and first-excited-state energy of the ethylene molecule. The one-dimensional many-body problem is solved numerically by means of the density-matrix renormalization-group (DMRG) method. The ground-state energy $E$ of the infinite open polymer is obtained by extrapolation considering chains of up to 300 atoms. The dimerization $\ensuremath{\delta}{=(r}_{i,i+1}\ensuremath{-}{r}_{i\ensuremath{-}1,i})/2$ and the average bond length ${R=(r}_{i,i+1}{+r}_{i\ensuremath{-}1,i})/2$ are determined by minimizing $E$. While $R$ is not significantly affected by the value of $U$, $\ensuremath{\delta}$ shows a remarkable nonmonotonic behavior. The differences between the delocalized regime, small $U/t$, and spin-Peierls behavior, large $U/t$, are discussed. The convergence properties of the DMRG algorithm are analyzed as a function of $U$, $\ensuremath{\delta}$, and the number states kept at each RG step. Criteria are examined in order to control the accuracy of the DMRG calculations in both weakly and strongly correlated situations.

Physical review. B, Condensed matter•Aug 15, 1997•21 citations

Electronic properties of the dimerized one-dimensional Hubbard model using lattice density-functional theory

Román López-Sandoval, G. M. Pastor

The dimerized one-dimensional Hubbard model is studied in the framework of lattice density-functional theory (LDFT). The single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ with respect to the lattice sites is considered as the basic variable. The corresponding interaction-energy functional $W[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ is defined by Levy's constrained search. Exact numerical results are obtained for $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23}),$ where ${\ensuremath{\gamma}}_{12}={\ensuremath{\gamma}}_{i,i+1}$ for odd i and ${\ensuremath{\gamma}}_{23}={\ensuremath{\gamma}}_{i,i+1}$ for even i are the nearest-neighbor density-matrix elements along the chain. The domain of representability of ${\ensuremath{\gamma}}_{\mathrm{ij}}$ and the functional dependence of $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23})$ are analyzed. A simple, explicit approximation to $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23})$ is proposed, which is derived from scaling properties of W, exact dimer results, and known limits. Using this approximation, LDFT is applied to determine ground-state properties and charge-excitation gaps of finite and infinite dimerized chains as a function of the Coulomb-repulsion strength $U/t$ and of the alternation $\ensuremath{\delta}t$ of the hopping integrals ${t}_{\mathrm{ij}}$ ${(t}_{\mathrm{ij}}=t\ifmmode\pm\else\textpm\fi{}\ensuremath{\delta}t).$ The accuracy of the method is demonstrated by comparison with available exact solutions and accurate numerical calculations. Goals and limitations of the present approach are discussed particularly concerning its ability to describe the crossover from weak to strong electron correlations.

arXiv:cond-mat/0210658•Physical review. B, Condensed matter•Jan 24, 2003•24 citations
By applying the numerically accurate symmetrized density-matrix renormalization-group method coupled with the extended Hubbard-Peierls model, we find that (i) the on-site Hubbard repulsion energy U dramatically reduces the binding energy of the lowest optically allowed 1Bu{1B}_{u}1Bu​ exciton; (ii) in the zero-dimerization limit, there exists a critical value of V at which the 1Bu{1B}_{u}1Bu​ exciton becomes bound; the critical value Vc=2t{V}_{c}=2tVc​=2t is fully in agreement with the recent analytical results at the infinite-U limit by Gallagher and Mazumdar [Phys. Rev. B 56, 15 025 (1997)], furthermore, this critical value decreases appreciably for weaker on-site correlation strengths, when the dimerization amplitude (\ensuremathδ)(\ensuremath{\delta})(\ensuremathδ) is nonzero. The present accurate numerical results contradict those obtained recently by Yu, Saxena, and Bishop [Phys. Rev. B 56, 3697 (1997)] both qualitatively and quantitatively. We also present first-order perturbation plus random-phase-approximation and single configuration-interaction analyses to rationalize the numerical calculations.
Using a symmetrized density-matrix renormalization-group formulation, we have investigated the electronic binding energy of the lowest optically allowed exciton (1Bu{1\mathrm{B}}_{\mathrm{u}}1Bu​) within an extended Hubbard-Peierls model (with parameters U, V, and \ensuremath{\delta}) for conjugated chains with N=80 sites. Three symmetries, the C2{\mathrm{C}}_{2}C2​, spin parity, and electron-hole symmetries, have been applied to construct the projector operator. Analysis of our results on the ratios between exciton binding energies and exciton energies, sheds light on the current experimental and theoretical controversy on exciton binding for polyacetylene and polyparaphenylene vinylene. We show that in the absence of dimerization, the exciton binding energy is vanishingly small for U/t up to 5 and V/t up to 2.3; for a finite dimerization, it is only when V is large enough that the exciton gets to be significantly bound. This result questions the applicability to conjugated polymers of the strong correlation picture, in which the exciton binding energy is equal to V.
We show that charge excitations within one-dimensional extended Hubbard models in the limit of large on-site Coulomb interaction can be mapped onto pairs of spin reversals from the ferromagnetic Ising-Heisenberg spin Hamiltonian. The exciton in the extended Hubbard models is equivalent to the bound two-magnon state of the spin Hamiltonian. From exact results for the spin Hamiltonian, for the intersite Coulomb interactions limited between nearest neighbors, the lowest optical exciton is separated in energy from the electron-hole continuum only for the intersite interaction larger than half the one-electron bandwidth. The commutation relationship between spin operators allows us to numerically investigate exciton binding energies and optical absorptions from the ground state in the limit of large on-site Coulomb interaction even for cases where analytic results are not available. For intersite Coulomb interaction limited between nearest neighbors, the binding energy of the lowest optical exciton changes very weakly for realistic Peierls bond dimerization, although there is a profound change in the optical absorption. While in the undimerized chain a single optical exciton can occur, two additional excitons that occur above the threshold of the electron-hole continuum become optically visible in the dimerized case. For long-range intersite Coulomb interactions, multiple excitons can occur even for the undimerized chain, and additional exciton states become optically visible upon dimerization. With the 1/r1/r1/r form for the intersite Coulomb interaction, our numerical calculations suggest that the formation of at least one optical exciton is unconditional, unlike the case of nearest-neighbor interaction. The binding energy of the lowest optical exciton again depends weakly on the dimerization for long-range interaction. Our theoretical results can explain the occurrence of two charge-transfer absorption bands in half-filled band segregated stack charge transfer solids below the spin-Peierls dimerization temperature. Our results also suggest that the lowest optical state of \ensuremath{\pi}-conjugated polymers, in which Coulomb interactions are long ranged, is necessarily an exciton.
We study the elementary excitations of a model Hamiltonian for the $\ensuremath{\pi}$ electrons in polydiacetylene chains. In these materials, the bare band gap is only half the size of the observed single-particle gap and the binding energy of the exciton of 0.5 eV amounts to 20% of the single-particle gap. Therefore, exchange and correlations due to the long-range Coulomb interaction require a numerically exact treatment which we carry out using the density-matrix renormalization group method. Employing both the Hubbard-Ohno potential and the screened potential in one dimension, we reproduce the experimental results for the binding energy of the singlet exciton and its polarizability. Our results indicate that there are optically dark states below the singlet exciton, in agreement with experiment. In addition, we find a weakly bound second exciton with a binding energy of 0.1 eV. The energies in the triplet sector do not match the experimental data quantitatively, probably because we do not include polaronic relaxation effects.
For the one-dimensional, extended Peierls–Hubbard model we calculate analytically the ground-state energy and the single-particle gap to second order in the Coulomb interaction for a given lattice dimerization. The comparison with numerically exact data from the density-matrix renormalization group shows that the ground-state energy is quantitatively reliable for Coulomb parameters as large as the bandwidth. The single-particle gap can almost triple from its bare Peierls value before substantial deviations appear. For the calculation of the dominant optical excitations, we follow two approaches. In Wannier theory, we perturb the Wannier exciton states to second order. In two-step perturbation theory, similar in spirit to the GW -BSE approach (where the Bethe–Salpeter equation is solved for quasi-particles which are calculated in the GW -approximation), we form excitons from dressed electron–hole excitations. We find the Wannier approach to be superior to the two-step perturbation theory. For singlet excitons, Wannier theory is applicable up to Coulomb parameters as large as half the bandwidth. For triplet excitons, second-order perturbation theory quickly fails completely.
The dimerization of polyacetylene is studied using a distance-dependent Hubbard Hamiltonian for the valence $\ensuremath{\pi}$ electrons and approximating $\ensuremath{\sigma}$-electron contributions by a sum of pair interactions ${E}_{\ensuremath{\sigma}}{(r}_{\mathrm{ij}})$. The $\ensuremath{\pi}$-electron Coulomb repulsion UUU is taken to be independent of the C-C distances rij{r}_{\mathrm{ij}}rij​ and is varied systematically from the uncorrelated or H\"uckel limit all over to the strongly correlated or Heisenberg regime. For each value of UUU, ${E}_{\ensuremath{\sigma}}{(r}_{\mathrm{ij}})$, and the $\ensuremath{\pi}$-electron hopping integral t(rij){t(r}_{\mathrm{ij}})t(rij​) are fitted to reproduce ab initio results on the ground-state and first-excited-state energy of the ethylene molecule. The one-dimensional many-body problem is solved numerically by means of the density-matrix renormalization-group (DMRG) method. The ground-state energy EEE of the infinite open polymer is obtained by extrapolation considering chains of up to 300 atoms. The dimerization $\ensuremath{\delta}{=(r}_{i,i+1}\ensuremath{-}{r}_{i\ensuremath{-}1,i})/2$ and the average bond length ${R=(r}_{i,i+1}{+r}_{i\ensuremath{-}1,i})/2$ are determined by minimizing EEE. While RRR is not significantly affected by the value of UUU, $\ensuremath{\delta}$ shows a remarkable nonmonotonic behavior. The differences between the delocalized regime, small U/tU/tU/t, and spin-Peierls behavior, large U/tU/tU/t, are discussed. The convergence properties of the DMRG algorithm are analyzed as a function of UUU, $\ensuremath{\delta}$, and the number states kept at each RG step. Criteria are examined in order to control the accuracy of the DMRG calculations in both weakly and strongly correlated situations.
The dimerized one-dimensional Hubbard model is studied in the framework of lattice density-functional theory (LDFT). The single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ with respect to the lattice sites is considered as the basic variable. The corresponding interaction-energy functional $W[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ is defined by Levy's constrained search. Exact numerical results are obtained for $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23}),$ where ${\ensuremath{\gamma}}_{12}={\ensuremath{\gamma}}_{i,i+1}$ for odd i and ${\ensuremath{\gamma}}_{23}={\ensuremath{\gamma}}_{i,i+1}$ for even i are the nearest-neighbor density-matrix elements along the chain. The domain of representability of ${\ensuremath{\gamma}}_{\mathrm{ij}}$ and the functional dependence of $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23})$ are analyzed. A simple, explicit approximation to $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23})$ is proposed, which is derived from scaling properties of W, exact dimer results, and known limits. Using this approximation, LDFT is applied to determine ground-state properties and charge-excitation gaps of finite and infinite dimerized chains as a function of the Coulomb-repulsion strength U/tU/tU/t and of the alternation $\ensuremath{\delta}t$ of the hopping integrals tij{t}_{\mathrm{ij}}tij​ ${(t}_{\mathrm{ij}}=t\ifmmode\pm\else\textpm\fi{}\ensuremath{\delta}t).$ The accuracy of the method is demonstrated by comparison with available exact solutions and accurate numerical calculations. Goals and limitations of the present approach are discussed particularly concerning its ability to describe the crossover from weak to strong electron correlations.