A Peierls Transition in Long Polymethine Molecular Wires: Evolution of Molecular Geometry and Single-Molecule Conductance

Molecules capable of mediating charge transport over several nanometers with minimal decay in conductance have fundamental and technological implications. Polymethine cyanine dyes are fascinating molecular wires because up to a critical length, they have no bond-length alternation (BLA) and their electronic structure resembles a one-dimensional free-electron gas. Beyond this threshold, they undergo a symmetry-breaking Peierls transition, which increases the HOMO-LUMO gap. We have investigated cationic cyanines with central polymethine chains of 5-13 carbon atoms (Cy3⁺-Cy11⁺). The absorption spectra and crystal structures show that symmetry breaking is sensitive to the polarity of the medium and the size of the counterion. X-ray crystallography reveals that Cy9·PF₆ and Cy11·B(C₆F₅)₄ are Peierls distorted, with high BLA at one end of the π-system, away from the partially delocalized positive charge. This pattern of BLA distribution resembles that of solitons in polyacetylene. The single-molecule conductance is essentially independent of molecular length for the polymethine salts of Cy3⁺-Cy11⁺ with the large B(C₆F₅)₄^- counterion, but with the PF₆^- counterion, the conductance decreases for the longer molecules, Cy7⁺-Cy11⁺, because this smaller anion polarizes the π-system, inducing a symmetry-breaking transition. At higher bias (0.9 V), the conductance of the shorter chains, Cy3⁺-Cy7⁺, increases with length (negative attenuation factor, β = -1.6 nm^-1), but the conductance still drops in Cy9⁺ and Cy11⁺ with the small polarizing PF₆^- counteranion.

Electronic charge density analysis of Li-doped polyacetylene: molecular vs periodic descriptions and nature of Li-to-chain bonding

A detailed analysis of the electronic structure and charge distribution around the trigonal site of Li-doped polyacetylene is reported using finite chain and periodic descriptions of the polymer. Atoms-in-molecules (AIM) analysis is done to characterize the nature of the bond between Li and the polymer backbone through the location of the bond critical points and computation of the total charge on the atomic basins around the doping site. We find that the Li atom donates practically one electron to the π-system, in accordance with the classical Su-Schriffer-Heeger model. However, despite that the Li atom is equidistant from the three closest C atoms in the geometric soliton, a single Li-C bond critical point is found. The AIM quantitative analysis of the electronic density reveals that the Li(+) ion is immersed into the polymer π-cloud in a way that resembles a metallic bonding interaction.

A theoretical investigation of electric properties of L-arginine phosphate monohydrate including environment polarization effects

The dipole moment (μ), linear polarizability (α), and first hyperpolarizability (β(tot)) of the asymmetric unit of L-arginine phosphate (LAP) monohydrate crystal are investigated using the supermolecule approach in combination with an iterative electrostatic polarization scheme. Environment polarization effects are attained by assuring the convergence of the dipole moment of LAP embedded in the polarization field of the surrounding molecules whose atomic sites are treated as point charges. The results obtained show that in the presence of the embedding charges, the value of μ is increased by 9% but the static values of α and β(tot) are decreased, respectively, by 3% and 13%, as compared with the isolated situation. The MP2/6-311+G(d) model predicts for the in-crystal dipole moment the converged value of 33 D, in good concordance with the available experimental result of 32 D. Our estimates for the converged results of α and β(tot) are, respectively, 22.51×10(-24) and 5.01×10(-30) esu. Dispersion effects are found to have a small impact on the nonlinear optical responses of LAP in the visible region. In addition, MP2/6-311G results obtained for β(tot) by using isolated and embedded LAP dimers show that crystal packing effects have a significant contribution of the electrostatic interactions. Our results suggest that the role of the crystal environment is to minimize the effects of the intermolecular interactions in the electric properties. That is, μ and β(tot) gain a more additive character in the presence of the field of the embedding charges. This is specially marked for β(tot).

On the accurate calculation of polarizabilities and second hyperpolarizabilities of polyacetylene oligomer chains using the CAM-B3LYP density functional

The polarizability and second hyperpolarizability of polyacetylene oligomer chains of increasing size up to C(24)H(26) were investigated by means of the Coulomb-attenuating method (CAM-B3LYP) using response theory. It was found that this long-range corrected density functional removes to large parts the overestimation observed for standard methods and in many cases provides results close to those of coupled cluster calculations. A direct comparison to experimentally observed dynamic hyperpolarizabilities is made to estimate the accuracy of the method. A basis set study revealed a noticeable contribution of diffuse orbitals to the hyperpolarizability also for larger oligomers. Furthermore, CAM-B3LYP is also confirmed to provide molecular geometries close to experimentally observed structures, especially for longer chain lengths.

Dynamical simulations of charged soliton transport in conjugated polymers with the inclusion of electron-electron interactions

We present numerical studies of the transport dynamics of a charged soliton in conjugated polymers under the influence of an external time-dependent electric field. All relevant electron-phonon and electron-electron interactions are nearly fully taken into account by simulating the monomer displacements with classical molecular dynamics and evolving the wave function for the pi electrons by virtue of the adaptive time-dependent density matrix renormalization group simultaneously and nonadiabatically. It is found that after a smooth turn on of the external electric field the charged soliton is accelerated at first up to a stationary constant velocity as one entity consisting of both the charge and the lattice deformation. An Ohmic region (6 mV/A</=E(0)</=12 mV/A) where the stationary velocity increases linearly with the electric field strength is observed. The relationship between electron-electron interactions and charged soliton transport is also investigated in detail. We find that the dependence of the stationary velocity of a charged soliton on the on-site Coulomb interactions U and the nearest-neighbor interactions V is due to the extent of delocalization of the charged soliton defect.

Theoretical Investigation of Static Characterization on Nonlinear Elementary Excitations in trans-Polyacetylene

Semiempirical quantum chemical studies on neutral and positively charged H(CH)nH homologues have been performed for systems with n up to 101, where different kinds of nonlinear excitation are found with increasing chain length. The Pariser-Parr-Pople (PPP) model has been employed and solved with the density matrix renormalization group (DMRG) method. The geometrical and electronic distortions induced by defects are obtained and compared with previous theoretical work, indicating that an adequate account of the electron correlation is essential for describing such systems. The structural distortion of a charged soliton (half-width calculated as L = 13) is shown to be more extended than that of a neutral soliton (L = 6); the geometrical distortion is even more extended in a polaron. In bipolarons, our calculations show that the coupling of the soliton-antisoliton pair might be longer ranged than expected. The phase transition from a bipolaron to a separated soliton-antisoliton pair occurs when n is close to 100.

Periodic Hartree-Fock and density functional theory calculations for Li-doped polyacetylene chains

We have performed periodic restricted Hartree-Fock/6-31G** and B3LYP6-31G** density functional theory calculations on Li-doped trans-polyacetylene at various dopant concentrations, using C(2m)H(2m)Li2 unit cells (m = 7-14). Except for maintaining P1 rod symmetry the geometry was completely optimized for both uniform and nonuniform doping structures. In addition to geometry we obtain atomic charges, along with soliton formation and dopant binding energies, as well as band structures and densities of states. A thorough analysis of the band structure and density of states, as a function of dopant concentration, is presented. We also characterize the complex nature of the binding interaction between Li and the polyacetylene chain.

Theoretical study of the lowest π→π* excitation energies for neutral and doped polyenes

In earlier theoretical studies, it has been widely noticed that the electron correlation effect played an important role in determining the excitation energies of low-lying pi-->pi(*) excited states for neutral polyenes and their radical cations and dications. In this paper, neutral and doped polyene oligomers of medium to large sizes are investigated with the Pariser-Parr-Pople model, and the pi-electron correlation effect is fully taken into consideration by virtue of the density-matrix renormalization group method. The excitation properties in the polymer limit are also obtained by exponential extrapolation from the finite oligomers. The reasonable agreement of our results with the available experimental observations and advanced ab initio calculations is witnessed. It is also observed that while charge doping can significantly lower the exciting energy, the odd-charged oligomers show lower excitation energies than the even-charged ones.