We extend the symmetrized density-matrix renormalization-group method to compute the dynamic nonlinear optical coefficients for long chains. By computing correction vectors in the appropriate symmetry subspace, we obtain the dynamic polarizabilities, ${\ensuremath{\alpha}}_{\mathrm{ij}}(\ensuremath{\omega})$, and third-order polarizabilities ${\ensuremath{\gamma}}_{\mathrm{ijkl}}(\ensuremath{\omega},\ensuremath{\omega},\ensuremath{\omega})$ of the Hubbard and ``U-V'' chains in an all transpolyacetylene geometry, with and without dimerization. We rationalize the behavior of $\overline{\ensuremath{\alpha}}$ and $\overline{\ensuremath{\gamma}}$ on the basis of the low-lying excitation gaps in the system.
We extend the density-matrix renormalization-group (DMRG) method to exploit parity,
C2 (rotation by \ensuremath{\pi}), and electron-hole symmetries of model Hamiltonians. We demonstrate the power of this method by obtaining the lowest-energy states in all eight symmetry subspaces of Hubbard chains with up to 50 sites. The ground-state energy, optical gap, and spin gap of regular U=4t and U=6t Hubbard chains agree very well with exact results. This development extends the scope of the DMRG method and allows future applications to study of optical properties of low-dimensional conjugated polymeric systems. \textcopyright{} 1996 The American Physical Society.
We extend the density matrix renormalization group method to exploit Parity,
C2 (rotation by
π) and electron-hole symmtries of model Hamiltonians. We demonstrate the power of this method by obtaining the lowest energy states in all the eight symmetry subspaces of Hubbard chains with upto 50 sites. The ground-state energy, optical gap and spin gap of regular
U=4t and
U=6t Hubbard chains agree very well with exact results. This development extends the scope of the DMRG method and allows future applications to study of optical properties of low-dimensional conjugated polymeric systems.
The density-matrix renormalization group (DMRG) method is used for a comparative study of low-lying excitations in trans-polyacetylene (t-PA) and transversely substituted t-PA (TS-t-PA). We have employed the Pariser-Parr-Pople model Hamiltonian which incorporates long-range electronic correlations to model these systems. We find some fundamental differences in the excited states of the t-PA and TS-t-PA. We find that the lowest two-photon allowed excited state in TS-t-PA is not made up of two triplet excitons and the gap to this state is nonzero even for undimerized chains in the thermodynamic limit. Contrary to earlier results for the Hubbard model, we find that the lowest two-photon state is always below the first optically allowed state in all the systems studied here making TS-t-PA systems only weakly fluorescent materials. Nonresonant tumbling averaged linear and third harmonic generation optic coefficients of TS-t-PA systems are also much smaller than that of t-PA.
A density-matrix renormalization-group (DMRG) method for calculating dynamical properties and excited states in low-dimensional lattice quantum many-body systems is presented. The method is based on an exact variational principle for dynamical correlation functions and the excited states contributing to them. This dynamical DMRG is an alternate formulation of the correction vector DMRG but is both simpler and more accurate. The finite-size scaling of spectral functions is discussed and a method for analyzing the scaling of dense spectra is described. The key idea of the method is a size-dependent broadening of the spectrum. The dynamical DMRG and the finite-size scaling analysis are demonstrated on the optical conductivity of the one-dimensional Peierls-Hubbard model. Comparisons with analytical results show that the spectral functions of infinite systems can be reproduced almost exactly with these techniques. The optical conductivity of the Mott-Peierls insulator is investigated and it is shown that its spectrum is qualitatively different from the simple spectra observed in Peierls (band) insulators and one-dimensional Mott-Hubbard insulators.
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
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
U, ${E}_{\ensuremath{\sigma}}{(r}_{\mathrm{ij}})$, and the $\ensuremath{\pi}$-electron hopping integral
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
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.