Vibrational analyses oftrans‐polyacetylene based onabinitiosecond‐order Mo/ller–Plesset perturbation calculations oftrans‐oligoenes

Molecular Structures and Vibrational Spectra of trans- and cis-Polyacetylene and Their Oligoenes Revisited Using Density Functional Theory Calculations

Polyacetylene, the most representative synthetic conducting polymer, has attracted much attention because it exhibits high conductivity upon doping. In this paper, molecular structures, electronic excitation energies, and Raman and infrared spectra were calculated using density functional theory for trans- and cis-oligoenes with various chain lengths up to the number of C═C bonds (n) of 100 and trans- and cis-polyacetylenes under one-dimensional periodic boundary condition. The harmonic vibrational frequencies obtained at the B3LYP/6-311G(d,p) level were scaled by the scaling factors determined with respect to the anharmonic vibrational frequencies using the B2PLYP method, in which the coefficients of the functional were optimized for trans-oligoenes. The calculated infrared and Raman frequencies reproduce reasonably well the observed frequencies for trans- and cis-polyacetylene. Based on the chain-length dependence of the calculated Raman spectra of trans-oligoenes, we proposed the possibility of longer conjugated trans-segments observed in the resonance Raman spectra of trans-polyacetylene excited at longer wavelengths of 647.1 and 1064 nm. We also elucidated the origin of the excitation-wavelength dependence of the resonance Raman spectra of trans-polyacetylene and the structure of isomerization intermediates from cis-form to trans-form. In addition, the previous assignments of Raman and infrared spectra of trans- and cis-polyacetylene were reexamined in the present study based on the chain-length dependence of the spectra.

Exploiting the Momentum Distribution in Atomically Confined Plasmonic Fields by Inelastic Scatterings

The utilization of atomically confined plasmonic fields has revolutionized the imaging technique. According to the fundamental position-momentum uncertainty principle, such a narrow spatial distribution certainly leads to a broad momentum distribution in the fields, which has however been overlooked. Here we propose a novel exploitation for the momentum distribution by adaptively satisfying the conservation law of momentum in inelastic Raman scatterings in periodic systems, providing a unique optical means of directly measuring the whole phonon dispersions. The proposed technique is particularly useful for measuring phonon dispersions of low-dimensional hydrogen-rich materials, which are completely inaccessible via other techniques. The numerical results for a single all-trans polyacetylene chain demonstrate that all phonon dispersion branches can be conclusively measured from their Raman images for the first time. Our findings highlight a unique advantage of the emerging momentum-based nanophotonics and open the door for exploiting highly confined plasmonic fields in another dimension.

DFT studies on the structural and vibrational properties of polyenes

Detailed density functional theory (DFT) calculations on the structure and harmonic frequencies of model all-trans and all-cis polyenes were undertaken. For the first time, we report on the convergence of selected B3LYP/6-311++G** and BLYP/6-311++G** calculated structural parameters resulting from a systematic increase in polyene size (chains containing 2 to 14 C = C units). The limiting values of the structural parameters for very long chains were estimated using simple three-parameter empirical formulae. BLYP/6-311++G** calculated ν(C = C) and ν(C-C) frequencies for all-trans and all-cis polyenes containing up to 14 carbon-carbon double bonds were used to estimate these values for very long chains. Correction of raw, unscaled vibrational data was performed by comparing theoretical and experimental wavenumbers for polyenes chains containing 3 to 12 conjugated C = C units with both ends substituted by tert-butyl groups. The corrected ν(C = C) and ν(C-C) wavenumbers for all-trans molecules were used to estimate the presence of 9 - 12 C = C units in all-trans polyene pigment in red coral.

Ground State Geometries of Polyacetylene Chains from Many-Particle Quantum Mechanics

Due to the crucial role played by electron correlation, the accurate determination of ground state geometries of π-conjugated molecules is still a challenge for many quantum chemistry methods. Because of the high parallelism of the algorithms and their explicit treatment of electron correlation effects, Quantum Monte Carlo calculations can offer an accurate and reliable description of the electronic states and of the geometries of such systems, competing with traditional quantum chemistry approaches. Here, we report the structural properties of polyacetylene chains H-(C₂H₂)(N)-H up to N = 12 acetylene units, by means of Variational Monte Carlo (VMC) calculations based on the multi-determinant Jastrow Antisymmetrized Geminal Power (JAGP) wave function. This compact ansatz can provide for such systems an accurate description of the dynamical electronic correlation as recently detailed for the 1,3-butadiene molecule [J. Chem. Theory Comput. 2015 11 (2), 508-517]. The calculated Bond Length Alternation (BLA), namely the difference between the single and double carbon bonds, extrapolates, for N → ∞, to a value of 0.0910(7) Å, compatible with the experimental data. An accurate analysis was able to distinguish between the influence of the multi-determinantal AGP expansion and of the Jastrow factor on the geometrical properties of the fragments. Our size-extensive and self-interaction-free results provide new and accurate ab initio references for the structures of the ground state of polyenes.

π-Conjugation in trans-1,3-butadiene: static and dynamical electronic correlations described through quantum Monte Carlo

We investigate the effects of the static and dynamical electronic correlations on the level of conjugation of the trans-1,3-butadiene molecule through Quantum Monte Carlo methods applied to an Antisymmetrized Geminal Power (AGP) wave function, with a Jastrow factor similar to the Gutzwiller ansatz. The degree of conjugation is measured through the convergence of the structural properties of 1,3-butadiene and in particular of the Bond Length Alternation (BLA), that is the difference between the lengths of the single and double carbon bonds. After verifying the different roles of the Fermionic AGP part of our wave function and of the Jastrow factor in recovering electronic correlation, we study the effects of a constrained Active Space AGP (AGPAS), similar to that used in the Complete Active Space (CAS) representation. Through this AGPAS, we are able to identify the effect of the limited active space on the degree of conjugation, showing that in the limit of infinite active space the structural properties converge exactly to those of the atomic AGP, giving a BLA for 1,3-butadiene around 0.1244(5) Å.

First-principles calculations on anharmonic vibrational frequencies of polyethylene and polyacetylene in the Gamma approximation

The frequencies of the infrared- and/or Raman-active (k=0) vibrations of polyethylene and polyacetylene are computed by taking account of the anharmonicity in the potential energy surfaces (PESs) and the resulting phonon-phonon couplings explicitly. The electronic part of the calculations is based on Gaussian-basis-set crystalline orbital theory at the Hartree-Fock and second-order Møller-Plesset (MP2) perturbation levels, providing one-, two-, and/or three-dimensional slices of the PES (namely, using the so-called n-mode coupling approximation with n=3), which are in turn expanded in the fourth-order Taylor series with respect to the normal coordinates. The vibrational part uses the vibrational self-consistent field, vibrational MP2, and vibrational truncated configuration-interaction (VCI) methods within the Gamma approximation, which amounts to including only k=0 phonons. It is shown that accounting for both electron correlation and anharmonicity is essential in achieving good agreement (the mean and maximum absolute deviations less than 50 and 90 cm(-1), respectively, for polyethylene and polyacetylene) between computed and observed frequencies. The corresponding values for the calculations including only one of such effects are in excess of 120 and 300 cm(-1), respectively. The VCI calculations also reproduce semiquantitatively the frequency separation and intensity ratio of the Fermi doublet involving the nu(2)(0) fundamental and nu(8)(pi) first overtone in polyethylene.

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.

Application of a Novel Linear/Exponential Hybrid Force Field Scaling Scheme to the Longitudinal Raman Active Mode of Polyyne

The properties of an infinite carbon chain (polyyne), an allotropic form of elemental carbon, are of importance in materials science as well as astronomy. The Raman active longitudinal optical (LO) frequencies are calculated with first-principles methods for oligoynes and polyyne and compared with experiments. Since traditional force constant scaling schemes fail in this case, we introduced a linear/exponential scaling scheme based on the exponential behavior of the carbon-carbon bond stretching force constant couplings in quasi-one-dimensional conjugated chains. The LO Raman active frequency is predicted at 1870-1877 cm-1. Our results provides further evidence for the assignment of the characteristic Raman peaks near 1850 cm-1 of the recently discovered long linear carbon chains encapsulated inside multiwalled or double-walled carbon nanotubes.

Molecular vibrations of [n]oligoacenes (n=2−5 and 10) and phonon dispersion relations of polyacene

As model compounds for nanosize carbon clusters, the phonon dispersion curves of polyacene are constructed based on density functional theory calculations for [n]oligoacenes (n=2-5, 10, and 15). Complete vibrational assignments are given for the observed Fourier-transform infrared and Raman spectra of [n]oligoacenes (n=2-5). Raman intensity distributions by the 1064-nm excitation are well reproduced by the polarizability-approximation calculations for naphthalene and anthracene, whereas several bands of naphthacene and pentacene at 1700-1100 cm(-1) are calculated to be enhanced by the resonance Raman effect. It is found from vibronic calculations that the coupled a(g) modes between the Kekulé deformation and joint CC stretching give rise to the Raman enhancements of the Franck-Condon type, and that the b(3g) mode corresponding to the graphite G mode is enhanced by vibronic coupling between the (1)L(a)((1)B(1u)) and (1)B(b)((1)B(2u)) states. The phonon dispersion curves of polyacene provide a uniform foundation for understanding molecular vibrations of the oligoacenes in terms of the phase difference. The mode correlated with the defect-sensitive D mode of the bulk carbon networks is also found for the present one-dimensional system.

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.

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.

An Efficient Fragment-Based Approach for Predicting the Ground-State Energies and Structures of Large Molecules

An efficient fragment-based approach for predicting the ground-state energies and structures of large molecules at the Hartree-Fock (HF) and post-HF levels is described. The physical foundation of this approach is attributed to the "quantum locality" of the electron correlation energy and the HF total energy, which is revealed by a new energy decomposition analysis of the HF total energy proposed in this work. This approach is based on the molecular fractionation with conjugated caps (MFCC) scheme (Zhang, D. W.; Zhang, J. Z. H. J. Chem. Phys. 2003, 119, 3599), by which a macromolecule is partitioned into various capped fragments and conjugated caps formed by two adjacent caps. We find that the MFCC scheme, if corrected by the interaction between non-neighboring fragments, can be used to predict the total energy of large molecules only from energy calculations on a series of small subsystems. The approach, named as energy-corrected MFCC (EC-MFCC), computationally achieves linear scaling with the molecular size. Our test calculations on a broad range of medium- and large molecules demonstrate that this approach is able to reproduce the conventional HF and second-order Moller-Plesset perturbation theory (MP2) energies within a few millihartree in most cases. With the EC-MFCC optimization algorithm described in this work, we have obtained the optimized structures of long oligomers of trans-polyacetylene and BN nanotubes with up to about 400 atoms, which are beyond the reach of traditional computational methods. In addition, the EC-MFCC approach is also applied to estimate the heats of formation for a series of organic compounds. This approach provides an appealing approach alternative to the traditional additivity rules based on either bond or group contributions for the estimation of thermochemical properties.

Raman spectroscopy of polyconjugated molecules and materials: confinement effect in one and two dimensions

The effect of the confinement of pi electrons in one- and two-dimensional domains is illustrated with several examples ranging from linear polyene chains to planar molecules with honeycomb structure. Theoretical computations and specific Raman experiments on molecular materials demonstrate that a molecular approach provides a unified key to the interpretation of the Raman response both of linear polyconjugated polymers (polyacetylene) and of nanostructured graphitic materials.

Origin of high second- and third-order nonlinear optical response in ammonio/borato diphenylpolyene zwitterions: the remarkable role of polarized aromatic groups

We present a quantum-chemical analysis of the molecular structure and second- and third-order polarizabilities in a series of promising nonlinear optical (NLO) chromophores, the zwitterionic ammonio/borato diphenylpolyenes, R3N+Ph(C=C)nPhB-R3, whose synthesis has been reported recently. The molecular geometries are obtained via MP2/6-31G optimization, while the NLO properties are calculated with the INDO Hamiltonian using the sum-over-states and finite-field real-space methods. The real-space approach allows the direct evaluation of the NLO-active segments of the molecules, while the sum-over-states results illustrate the virtual excitations and charge-transfer pathways that are essential in the NLO response. Both methods highlight the remarkable and unexpected result that it is the strongly polarized phenylene groups that play the key role in generating a high NLO response.

Experimental CC stretching phonon dispersion curves and electron phonon coupling in polyene derivatives

We analyze the infrared and Raman spectra (both experimentally and with the aid of quantum chemical calculations) of a series of polyenals which provide us with the fortunate case of a set of polyene chains with one of the end groups consisting of a C=O group which not only does take part in the conjugation but also pulls electrons from the chain making the whole system highly polar, thus affecting the vibrational transition moments. In the following we show, for the first time, that it is possible to derive experimental phonon dispersion curves and these prove to be different for each chain length. We support our experimental findings with Density Functional Theory quantum chemical calculations which reproduce with sufficient accuracy the IR and Raman spectral pattern and at the same time help in disentangling the assignment of the fine structure observed in the experimental spectra.