Structural and Electronic Properties of Polyacetylene and Polyyne from Hybrid and Coulomb-Attenuated Density Functionals

Bond length alternation of π-conjugated polymers predicted by the Fermi-Löwdin orbital self-interaction correction method

π-conjugated polymers have been used in a wide range of practical applications, partly due to their unique properties that originate in the delocalization of electrons through the polymer backbone. The level of delocalization can be characterized by the induced bond length alternation (BLA), with shorter BLA connected with strong delocalization and vice versa. The accurate description of this structural parameter can be considered a benchmark for testing the capability of different electronic structure methods for self-interaction error (SIE) removal and electron correlation inclusion. Density functional theory (DFT), in its local or semi-local flavors, suffers from SIE and, thus, underestimates the BLA compared to self-interaction-free methods. In this work, we utilize the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method for one-electron self-interaction removal to characterize the BLA of five oligomers with increasing length extrapolated to the polymeric limit. We compare the self-interaction-free BLA to several DFT approximations, Møller-Plesset second-order perturbation theory (MP2), and the BLA obtained with the domain based local pair natural orbital CCSD(T) [DLPNO-CCSD(T)] approximation. Our findings show that FLOSIC corrects for the small BLA given by (semi-)local DFT approximations, but it tends to overcorrect with respect to CAM-B3LYP, MP2, and DLPNO-CCSD(T).

Comparative rotatory power of bent and twisted polyynes

Linear polyynes of the formula C18 H2 (symmetry D∞h ) were bent in silico by progressively introducing CCC angles less than 180°. The bent structures (symmetry C2v ) were then twisted by introducing torsion angles across the CCCC segments by as much as 60°. The gyration tensors of these 19 structures (linear, bent, and twisted) were computed by linear response methods. Bending is massively generative of optical activity in oriented structures, even achiral structures, whereas twisting in conjunction with bending, serves to linearize the molecules and diminish maximally observable optical activity. This computational exercise is intended to unbind the infelicitous linkage of optical activity and chirality, which is only meaningful in isotropic media. Although bent structures are not optically active in solution-the spatial average of the optical activity is necessarily zero-solution measurements that deliver the spatial averages are a special class of measurements, albeit the overwhelmingly most common chiroptical measurements, that prejudice our common understanding of how π-conjugated structures generate gyration. Bending is far more effective than twisting at generating optical activity along some directions for oriented structures. The respective contributions from the transition electric dipole-magnetic dipole polarizability and the transition electric dipole-electric quadrupole polarizability are compared.

The vertical excitation energies and a lifetime of the two lowest singlet excited states of the conjugated polyenes from C2 to C22: Ab initio, DFT, and semiclassical MNDO-MD simulations

Electronic excited states in the series of polyene molecules were explored. Optimal ground-state geometry was used for the evaluation of vertical excitation energies. Results of a chosen set of functionals were compared to post-HF methods (EOM-CCSD, NEVPT2, CASPT2, and MRCI). In addition, the semiempirical OM2/MNDO method using MRCISD computational level was confronted with the above-mentioned techniques. Despite the fact that the first excited state has a significant double-excitation character some functionals were able to qualitatively determine the correct state order (where the lowest excited state has a A g - character). The most successful functionals in transition energies predictions were PBE, TPSS and BLYP in Tamm-Dancoff approach (TDA), which had the smallest root-mean-square deviation (RMSD) scoring towards the experimental values. Regarding RMSD scoring, the OM2/MNDO method performed fairly well, too. Besides absorption spectra, lifetimes of the first two excited states were estimated based on a stochastic approach exploring a swarm of OM2/MNDO hopping dynamics using the Tully fewest switch algorithm for each molecule. The longest lifetime of the first excited state (S₁ ) was found for decapentaene (about 5 ps). Further elongation of the conjugated chain caused a mild decrease of this value to ca 1.5 ps for docosaundecaene.

Theoretical Study on Non-Linear Optics Properties of Polycyclic Aromatic Hydrocarbons and the Effect of Their Intercalation with Carbon Nanotubes

Results of a theoretical study devoted to comparing NLO (non-linear optics) responses of derivatives of tetracene, isochrysene, and pyrene are reported. The static hyperpolarizability β, the dipole moment μ, the HOMO and LUMO orbitals, and their energy gap were calculated using the CAM-B3LYP density functional combined with the cc-pVDZ basis set. The para-disubstituted NO2-tetracene-N(CH3)2 has the highest NLO response, which is related to a large intramolecular charge transfer. Adding vinyl groups to the para-disubstituted NO2-tetracene-N(CH3)2 results in an increase in the NLO responses. We further investigated the effect of the intercalation of various push-pull molecules inside an armchair single-walled carbon nanotube. The intercalation leads to increased NLO responses, something that depends critically on the position of the guest molecule and/or on functionalization of the nanotube by donor and attractor groups.

Low-Lying Excited States of Natural Carotenoids Viewed by Ab Initio Methods

Low-lying excited states of carotenoids (the optically dark 2A_g^- and bright 1B_u⁺) are deeply involved in energy transfer processes in photosynthetic antennas, such as light harvesting and non-photochemical quenching. Though any ab initio modeling of these phenomena has to rely on precise energies of the carotenoid electronic states, the accurate evaluation of these states remains a challenging problem due to their different natures. The paper aims to study the accuracy of the excitation energies of the low-lying excited states of certain open- and closed-chain carotenoids obtained by a state-of-the-art multireference approach for electronic structure calculation. Here, density matrix renormalization group SCF (DMRGSCF) and a perturbative approach based on driven similarity renormalization group second-order multireference perturbation theory (DSRG-MRPT2) were used to treat the static and dynamic correlation, respectively. Nuclear geometries of the electronic states were optimized with DFT-based approaches. It is demonstrated that spin-flip TD-DFT can replace multiconfigurational methods for the geometry optimization of the 2A_g^- state but not for the calculation of the excitation energy. Adiabatic excitation energies to the 1B_u⁺ state were shown to be within a margin of 1000 cm^-1 with an appropriate flow parameter value. Adiabatic excitation energies to the 2A_g^- state for the open-chain carotenoids lie within a range of experimental values (taking into account the broad range of experimental estimates); for the closed-chain ones, the error does not exceed 2000 cm^-1. Ab initio stationary (1A_g^- → 1B_u⁺) and transient (2A_g^- → 1B_u⁺) absorption spectra were modeled for violaxanthin and lycopene, and these spectra showed good agreement with the experimental ones both in terms of the vibronic structure and the transition energies.

Examining fundamental and excitation gaps at the thermodynamic limit: A combined (QTP) DFT and coupled cluster study on trans-polyacetylene and polyacene

Interest in ab initio property prediction of π-conjugated polymers for technological applications places significant demand on “cost-effective” and conceptual computational methods, particularly effective, one-particle theories. This is particularly relevant in the case of Kohn–Sham Density Functional Theory (KS-DFT) and its new competitors that arise from correlated orbital theory, the latter defining the QTP family of DFT functionals. This study presents large, ab initio equation of motion-coupled cluster calculations using the massively parallel ACESIII to target the fundamental bandgap of two prototypical organic polymers, trans-polyacetylene (tPA) and polyacene (Ac), and provides an assessment of the new quantum theory project (QTP) functionals for this problem. Further results focusing on the 1 A g (1 A g), 1 B u (1 B2 u), and 3 B u (3 B2 u) excited states of tPA (Ac) are also presented. By performing calculations on oligomers of increasing size, extrapolations to the thermodynamic limit for the fundamental and all excitation gaps, as well as estimations of the exciton binding energy, are made. Thermodynamic-limit results for a combination of “optimal” and model geometries are presented. Calculated results for excitations that are adequately described using a single-particle model illustrate the benefits of requiring a KS-DFT functional to satisfy the Bartlett ionization potential theorem.

DFT investigation on the application of pure and doped X12N12 (X = B and Al) fullerene-like nano-cages toward the adsorption of temozolomide

The sensitivity of pure and doped X₁₂N₁₂ (X = B and Al) fullerene-like nano-cages (FLNs) toward the anti-cancer drug temozolomide (TMZ) is probed herein at DFT/M06-2X-D3/6-311G(d,p) theoretical level in both gas phase and water. A noticeable affinity of the FLNs toward TMZ was observed along with the negative gas-phase adsorption energies -1.37 and -2.09 eV for the most stable configurations of pure B₁₂N₁₂ and Al₁₂N₁₂ pristines, respectively. Considerable charge transfer from TMZ to the FLNs was also revealed via NBO analysis and the Hirshfeld atomic charges, making the dipole moment vector of the molecular complexes to be oriented from the nano-cages to the TMZ moiety. Furthermore, a percentage decrease in the HOMO-LUMO energy gap (ΔE _g) of 38.09 and 17.72% was obtained for the B₁₂N₁₂ and Al₁₂N₁₂ nano-cages, respectively. The percentage change in ΔE _g was found to be reduced upon doping and solvation of the FLNs. Finally, a recovery time in vacuum ultraviolet light of 1.06 s is found for the complex with pure B₁₂N₁₂, which in addition to the above-mentioned parameters make this boron nitride cage the best sensor for TMZ, among the FLNs considered in the present work.

Confinement on the optical response in h-BNCs: Towards highly efficient SERS-active 2D substrates

Several experimental and theoretical studies have shown that 2D hybrid structures formed by boron, nitrogen and carbon atoms (h-BNCs) possess a highly tunable linear and non-linear optical responses. Recent advances towards the controlled synthesis of these unique structures have motivated an important number of experimental and theoretical work. In this work, the confinement on the optical response induced by boron-nitride (BN) strings in h-BNC 2D structures is investigated using time-dependent density functional theory (TDDFT) and electron density response properties. The number of surrounding BN strings (N_BN) necessary to "isolate" the optical modes of a carbon nanoisland (nanographene) from the remaining substrate has been characterized in two different nanoisland models: benzene and pyrene. It was found that for N_BN ≥ 3 the excitation wavelengths of the optically active modes remain constant and the changes in the transition densities, the ground to excited state density differences and their associated electron deformation orbitals are negligible and strongly confined within the carbon nanoisland. Using a water molecule as model system, Raman enhancement factors of 10 [6] for the water vibrational modes are obtained when these electromagnetic "hot spots" are activated by an external electromagnetic field. The high tunability of the optical absorption bands of nanographenes through changes in size and morphology makes h-BNCs be perfect materials to construct platforms for surface enhancement Raman spectroscopy (SERS) for a wide range of laser sources.

Dominant Role of Quantum Anharmonicity in the Stability and Optical Properties of Infinite Linear Acetylenic Carbon Chains

Carbyne, an infinite-length straight chain of carbon atoms, is supposed to undergo a second order phase transition from the metallic bond-symmetric cumulene (═C═C═)_∞ toward the distorted insulating polyyne chain (-C≡C-)_∞ displaying bond-length alternation. However, recent synthesis of ultra long carbon chains (∼6000 atoms, [Nat. Mater., 2016, 15, 634]) did not show any phase transition and detected only the polyyne phase, in agreement with previous experiments on capped finite carbon chains. Here, by performing first-principles calculations, we show that quantum-anharmonicity reduces the energy gain of the polyyne phase with respect to the cumulene one by 71%. The magnitude of the bond-length alternation increases by increasing temperature, in stark contrast with a second order phase transition, confining the cumulene-to-polyyne transition to extremely high and unphysical temperatures. Finally, we predict that a high temperature insulator-to-metal transition occurs in the polyyne phase confined in insulating nanotubes with sufficiently large dielectric constant due to a giant quantum-anharmonic bandgap renormalization.

Excited States of Xanthophylls Revisited: Toward the Simulation of Biologically Relevant Systems

Xanthophylls are a class of oxygen-containing carotenoids, which play a fundamental role in light-harvesting pigment-protein complexes and in many photoresponsive proteins. The complexity of the manifold of the electronic states and the large sensitivity to the environment still prevent a clear and coherent interpretation of their photophysics and photochemistry. In this Letter, we compare cutting-edge ab initio methods (CC3 and DMRG/NEVPT2) with time-dependent DFT and semiempirical CI (SECI) on model keto-carotenoids and show that SECI represents the right compromise between accuracy and computational cost to be applied to real xanthophylls in their biological environment. As an example, we investigate canthaxanthin in the orange carotenoid protein and show that the conical intersections between excited states and excited-ground states are mostly determined by the effective bond length alternation coordinate, which is significantly tuned by the protein through geometrical constraints and electrostatic effects.

TD-DFT Simulation and Experimental Studies of a Mirrorless Lasing of Poly[(9,9-dioctylfluorenyl-2,7-diyl)-co-(1,4-diphenylene-vinylene-2-methoxy-5-{2-ethylhexyloxy}-benzene)]

In this work, we investigate the TD-DFT simulation, optical, and mirrorless laser properties of conjugated polymer (CP) Poly[(9,9-dioctylfluorenyl-2,7-diyl)-co-(1,4-diphenylene-vinylene-2-methoxy-5-{2-ethylhexyloxy}-benzene)], also known as (PFO-co-PPV-MEHB) or ADS125GE. TD-DFT calculations were performed for three monomer units with truncated tails using time-dependent density functional theory (TD-DFT) calculations. The calculations showed a highest occupied and lowest unoccupied molecular orbital (HOMO-LUMO) structure and a very high oscillator strength of 6.434 for the singlet-singlet transition at 374.43 nm. Experimentally, the absorption and fluorescence spectra were examined at various concentrations in verity of solvents, such as benzene, toluene, and hexane. The experimental results obtained in hexane were comparable with theoretical UV-VIS spectra calculated under vacuum. Amplified spontaneous emission (ASE) spectra peaked at approximately 509 nm for CO PFO-co-PPV-MEHB in solution and were obtained at suitable concentrations and pump energies. Additionally, the photochemical stability of this CP and coumarin (C510) were compared. Time-resolved spectroscopy (TRS) studies with a sub-nanosecond resolution were performed for the CO under various pump energies. These results showed the excited state dynamics and single-pass optical gain of CO PFO-co-PPV-MEHB.

Effect of varying the TD-lc-DFTB range-separation parameter on charge and energy transfer in a model pentacene/buckminsterfullerene heterojunction

Atomistic modeling of energy and charge transfer at the heterojunction of organic solar cells is an active field with many remaining outstanding questions owing, in part, to the difficulties in performing reliable photodynamics calculations on very large systems. One approach to being able to overcome these difficulties is to design and apply an appropriate simplified method. Density-functional tight binding (DFTB) has become a popular form of approximate density-functional theory based on a minimal valence basis set and neglect of all but two center integrals. We report the results of our tests of a recent long-range correction (lc) [A. Humeniuk and R. Mitrić, J. Chem. Phys. 143, 134120 (2015)] for time-dependent (TD) lc-DFTB by carrying out TD-lc-DFTB fewest switches surface hopping calculations of energy and charge transfer times using the relatively new DFTBABY [A. Humeniuk and R. Mitrić, Comput. Phys. Commun. 221, 174 (2017)] program. An advantage of this method is the ability to run enough trajectories to get meaningful ensemble averages. Our interest in the present work is less in determining exact energy and charge transfer rates than in understanding how the results of these calculations vary with the value of the range-separation parameter (R_lc = 1/μ) for a model organic solar cell heterojunction consisting of a gas-phase van der Waals complex P/F made up of a single pentacene (P) molecule together with a single buckminsterfullerene (F) molecule. The default value of R_lc = 3.03 a₀ is found to be much too small as neither energy nor charge transfer is observed until R_lc ≈ 10 a₀. Tests at a single geometry show that the best agreement with high-quality ab initio spectra is obtained in the limit of no lc (i.e., very large R_lc). A plot of energy and charge transfer rates as a function of R_lc is provided, which suggests that a value of R_lc ≈ 15 a₀ yields the typical literature (condensed-phase) charge transfer time of about 100 fs. However, energy and charge transfer times become as high as ∼300 fs for R_lc ≈ 25 a₀. A closer examination of the charge transfer process P^*/F → P⁺/F^- shows that the initial electron transfer is accompanied by a partial delocalization of the P hole onto F, which then relocalizes back onto P, consistent with a polaron-like picture in which the nuclei relax to stabilize the resultant redistribution of charges.

Electrochemical and Solvent-Mediated Visible-to-Near-Infrared Spectroscopic Switching of Benzoselenadiazole Fluorophores

A series of electronic push-pull, pull-pull, and push fluorophores has been prepared from a benzoselenadiazole core so that their spectroscopic, electrochemical, spectro-electrochemical, and spectro-electrofluorescence properties could be examined. The emission wavelengths and fluorescence quantum yields (Φ_fl ) of the N,N-dimethyl fluorophores were contingent on the solvent polarity and they ranged from 615 to 850 nm in aprotic solvents. The positive solvatochromism and the quenched Φ_fl in polar solvents were consistent with an intramolecular charge-transfer state (ICT). Meanwhile, a locally excited state (LE) was assigned in nonpolar solvents from the blue-shifted emission and high Φ_fl . The N,N-dimethylamine fluorophores examined could be both electrochemically oxidized and reduced, whereas the symmetric dinitro pull-pull derivative could be only reversibly reduced. Courtesy of their electrochemical reversibility, the fluorophores could reversibly change color from yellow to blue with an applied potential in addition to switching off their emission. The absorption of the electrochemically generated intermediates of the N,N-dimethyl derivatives spanned 500 nm over the visible and the NIR regions. The colors could be switched for upwards of two hours with applied potential, illustrating their potential use as electroactive materials in electrochromic devices.

Optical gap and fundamental gap of oligoynes and carbyne

The optoelectronic properties of various carbon allotropes and nanomaterials have been well established, while the purely sp-hybridized carbyne remains synthetically inaccessible. Its properties have therefore frequently been extrapolated from those of defined oligomers. Most analyses have, however, focused on the main optical transitions in UV-Vis spectroscopy, neglecting the frequently observed weaker optical bands at significantly lower energies. Here, we report a systematic photophysical analysis as well as computations on two homologous series of oligoynes that allow us to elucidate the nature of these weaker transitions and the intrinsic photophysical properties of oligoynes. Based on these results, we reassess the estimates for both the optical and fundamental gap of carbyne to below 1.6 eV, significantly lower than previously suggested by experimental studies of oligoynes.

Carbyne, a linear sp-hybridized carbon allotrope, is synthetically inaccessible and its properties are extrapolated from those of defined oligomers. Here the authors analyze weak optical bands in two series of oligoynes and reassess the optical and fundamental gap of carbyne to lower values than previously suggested.

Controlled Cyclopolymerization of 1,5-Hexadiynes to Give Narrow Band Gap Conjugated Polyacetylenes Containing Highly Strained Cyclobutenes

For decades, cyclopolymerization of α,ω-diyne derivatives has been an effective method to synthesize various soluble polyacetylenes containing five- to seven-membered rings in the backbone. However, cyclopolymerization to form four-membered carbocycles was considered impossible due to their exceptionally high ring strain (∼30 kcal/mol). Herein, we demonstrate the successful cyclopolymerization of rationally designed 1,5-hexadiyne derivatives to afford various polyacetylenes containing highly strained cyclobutenes in each repeat unit. After screening, Ru catalysts containing bulky diisopropylphenyl groups promoted challenging four-membered ring cyclization efficiently from various monomers, enabling the synthesis of high molecular weight (up to 40 kDa) polyacetylenes in a controlled manner. Furthermore, living polymerization allowed for block copolymer synthesis by combining with ring-opening metathesis polymerization as well as block copolymerization of two different 1,5-hexadiyne monomers to give a fully conjugated polyacetylene. These new polymers unexpectedly showed much narrower band gaps than conventional substituted polyacetylenes by >0.2 eV. Interestingly, computational studies showed much smaller bond length alternation in the conjugated backbone containing cyclobutenes, resulting in highly delocalized π electrons along the polymer chain and lower band gaps.

Ab Initio Study of Low-Lying Excited States of Carotenoid-Derived Polyenes

Knowledge about excited states of carotenoids is essential for understanding photophysical processes underlying photosynthesis. However, due to the presence of a large number of optically dark states, experimental study of the excited-state manifold is limited to a significant extent. In this paper, we apply high-level ab initio quantum chemical methods to study the low-lying excited states of polyenes containing from 8 to 13 conjugated double bonds, which serve as a model for natural carotenoids. Vertical and adiabatic excitation energies from the ground 1A_g^- state to the excited 2A_g^-, 1B_u⁺, and 1B_u^- states were evaluated by means of density matrix renormalization group (DMRG) with NEVPT2 perturbative correction. The energies of all excited states are highly sensitive to nuclear geometry, especially the 2A_g^- state. Thus, the 2A_g^- and 1B_u⁺ states interchange their relative positions upon geometry relaxation, while the vertical excitation energy to the 2A_g^- state is rather high. At the same time, the 1B_u^- state energy is shown to be higher than other studied excited states at any geometry. With relaxed geometries of the excited states, absorption and transient absorption spectra were calculated within the Franck-Condon approximation bridging the gap between experimental spectroscopic data and computational results.

Range-separated hybrid and double-hybrid density functionals: A quest for the determination of the range-separation parameter

We recently derived a new and simple route to the determination of the range-separation parameter in range-separated exchange hybrid and double-hybrid density functionals by imposing an additional constraint to the exchange-correlation energy to recover the total energy of the hydrogen atom [Brémond et al., J. Chem. Phys. 15, 201102 (2019)]. Here, we thoroughly assess this choice by statistically comparing the derived values of the range-separation parameters to the ones obtained using the optimal tuning (OT) approach. We show that both approaches closely agree, thus, confirming the reliability of ours. We demonstrate that it provides very close performances in the computation of properties particularly prone to the one- and many-electron self-interaction errors (i.e., ionization potentials). Our approach arises as an alternative to the OT procedure, conserving the accuracy and efficiency of a standard Kohn-Sham approach to density-functional theory computation.

sp-hybridized carbon allotrope molecular structures: An ongoing challenge for density-functional approximations

The recent synthesis of a C₁₈ monocyclic ring constitutes a major breakthrough as a new all-carbon disclosed form. However, modern density functional theory approaches do not lead to the correct experimental polyynic structure and favor the cumulenic one instead. We demonstrate here that this serious drawback can be solved by recently developed range-separated nonempirical schemes, independently of which kind of functional is being applied (i.e., semilocal, hybrid, or double-hybrid).

Range-separated hybrid density functionals made simple

In this communication, we present a new and simple route to derive range-separated exchange (RSX) hybrid and double hybrid density functionals in a nonempirical fashion. In line with our previous developments [Brémond et al., J. Chem. Theory Comput. 14, 4052 (2018)], we show that by imposing an additional physical constraint to the exchange-correlation energy, i.e., by enforcing to reproduce the total energy of the hydrogen atom, we are able to generalize the nonempirical determination of the range-separation parameter to a family of RSX hybrid density functionals. The success of the resulting models is illustrated by an accurate modeling of several molecular systems and properties, like ionization potentials, particularly prone to the one- and many-electron self-interaction errors.

Reduced bond length alternation and helical molecular orbitals in donor and acceptor substituted linear carbon chains

Increasing chain length and end group substitution of polyynes play a crucial role in molecular electronics and nanomaterials. The studies on linear carbon chains are lesser when compared to other carbon allotropes like graphene, fullerenes, nanotube, etc. This prompted us to study the linear carbon chains of different lengths and substitutions. The electronic and optical properties of X–C[Formula: see text]–X ([Formula: see text]–15 and [Formula: see text], NH2, CN, OH) molecules have been studied by using CAM-B3LYP/6-31G* level of theory of DFT methods. Linear carbon chains with odd values of n show lower bond length alternation (BLA) values similar to that of cumulenes and may have metallic property, but the substitution of donor/acceptor molecules does not decrease the BLA significantly. Molecular orbital analysis of linear carbon chains shows that NH2 or NO2 substituted polyynes have helical molecular orbitals for smaller chain lengths which may make a good candidate for molecular wires in molecular devices. As the chain length increases, the helicity decreases and finally disappears. Also, it is seen that for smaller odd values of [Formula: see text] for donor, substituted polyynes have a singlet ground, whereas all the odd [Formula: see text] values of acceptor substitution have triplet ground state.

A computational study of the vibrationally-resolved electronic circular dichroism spectra of single-chain transoid and cisoid oligothiophenes in chiral conformations

We simulate the vibronic profile of the electronic circular dichroism (ECD) spectra of oligothiophenes in cisoid and transoid chiral arrangements. We consider oligomers of different lengths, from two to fifteen units, and investigate extensively how the ECD spectral shapes depend on the inter-ring torsions. In general, the molecular structures we consider are not stationary points of the ground state potential energy surface. Therefore, in order to perform vibronic calculations, we present a new computational protocol able to define reduced-dimensionality models where the effect of the off-equilibrium modes is removed. This is done adopting a description of the vibrational motions in curvilinear internal coordinates, and vertical harmonic models coupled with an iterative application of projectors to define energy Hessians, and therefore effective normal modes, in the space complementary to the one of the off-equilibrium coordinates. Although we consider both Franck-Condon and Herzberg-Teller contributions, the results show that transoid twisted ribbons always give rise to monosignated ECD spectra, while bi-signated and multi-signated spectra are expected for cisoid helices. These findings are explained on the basis of the different transition strengths of the lowest electronic states imparted by the different spatial arrangement, that is almost linear for transoid structures and more globular for cisoid ones. We predicted the chiroptical response of a large number of possible molecular arrangements. These data are employed to critically discuss the experimental ECD of polythiophenes in different experimental conditions, forming either aggregates or host-guest complexes. The method here proposed to perform vibronic calculations in reduced-dimensionality models is of general applicability and its potential interest goes beyond the practical application presented here.

Ultrafast dynamics of low-energy electron attachment via a non-valence correlation-bound state

The primary electron-attachment process in electron-driven chemistry represents one of the most fundamental chemical transformations with wide-ranging importance in science and technology. However, the mechanistic detail of the seemingly simple reaction of an electron and a neutral molecule to form an anion remains poorly understood, particularly at very low electron energies. Here, time-resolved photoelectron imaging was used to probe the electron-attachment process to a non-polar molecule using time-resolved methods. An initially populated diffuse non-valence state of the anion that is bound by correlation forces evolves coherently in ∼30 fs into a valence state of the anion. The extreme efficiency with which the correlation-bound state serves as a doorway state for low-energy electron attachment explains a number of electron-driven processes, such as anion formation in the interstellar medium and electron attachment to fullerenes.

Quantum-Assisted Metrology of Neutral Vitamins in the Gas Phase

It has recently been shown that matter-wave interferometry can be used to imprint a periodic nanostructure onto a molecular beam, which provides a highly sensitive tool for beam displacement measurements. Herein, we used this feature to measure electronic properties of provitamin A, vitamin E, and vitamin K1 in the gas phase for the first time. The shift of the matter-wave fringes in a static electric field encodes the molecular susceptibility and the time-averaged dynamic electric dipole moment. The dependence of the fringe pattern on the intensity of the central light-wave diffraction grating was used to determine the molecular optical polarizability. Comparison of our experimental findings with molecular dynamics simulations and density functional theory provides a rich picture of the electronic structures and dynamics of these biomolecules in the gas phase with β-carotene as a particularly interesting example.

Performance of Hybrid DFT Compared to MP2 Methods in Calculating Nonlinear Optical Properties of Divinylpyrene Derivative Molecules

A series of divinyl-pyrene derivatives of the form D-vinyl-pyrene-vinyl-A, in which D corresponds to an electron donor group and A to an electron acceptor group, were studied in this work. The first purpose was to determine the optimal HF % exchange as incorporated in a range of hybrid functionals (M06HF, M062X, M06L, CAM-B3LYP, PBE0, BMK, and B3LYP) capable to produce, reliably and as close as possible to those obtained from MP2 calculations, NLO parameters and, in particular, first-order static hyperpolarizabilities. The CAM-B3LYP functional was revealed to be the most suitable one. The pair N(CH₃)₂/NO₂ was then determined as the most efficient pair of groups in producing appreciable NLO responses. The effect of the substitution position on the pyrene moiety was also investigated, whereby aligning the two substituents involving the D and A groups in the direction of the dipole moment as in the (1,6 DVP) derivatives was shown to be most favorable for increasing the NLO parameters.

Experimental and computational evaluation of the barrier to torsional rotation in a butadiyne-linked porphyrin dimer

The barrier to torsional rotation in a butadiyne-linked porphyrin dimer has been determined in solution using variable temperature UV-vis-NIR spectroscopy: ΔH = 5.27 ± 0.03 kJ mol(-1), ΔS = 10.69 ± 0.14 J K(-1) mol(-1). The value of ΔH agrees well with theoretical predictions. Quantum chemical calculations (DFT) were used to predict the torsion angle dependence of the absorption spectrum, and to calculate the vibronic fine structure of the S0 → S1 absorption for the planar dimer, showing that the absorption band of the planar conformer has a vibronic component overlapping with the 〈0|0〉 absorption of the perpendicular conformer. The torsion barrier in the porphyrin dimer is higher than that of 1,4-diphenylbutadiyne (calculated ΔH = 1.1 kJ mol(-1)). Crystallographic bond lengths and IR vibrational frequencies confirm that there is a greater contribution of the cumulenic resonance form in butadiyne-linked porphyrin dimers than in 1,4-diphenylbutadiyne. The DFT frontier orbitals of the twisted conformer of the porphyrin dimer are helical, when calculated in the absence of symmetry. The helical character of these orbitals disappears when D2d symmetry is enforced in the 90° twisted conformer. Helical representations of the frontier orbitals can be generated by linear combinations of the more localised orbitals from a symmetry-constrained calculation but they do not indicate π-conjugation. This work provides insights into the relationship between electronic structure and conformation in alkyne-linked conjugated oligomers.

Excited-State Dipole and Quadrupole Moments: TD-DFT versus CC2

The accuracies of the excited-state dipole and quadrupole moments obtained by TD-DFT are assessed by considering 16 different exchange-correlation functionals and more than 30 medium and large molecules. Except for excited-state presenting a significant charge-transfer character, a relatively limited dependency on the nature of the functional is found. It also turns out that while DFT ground-state dipole moments tend to be too large, the reverse trend is obtained for their excited-state counterparts, at least when hybrid functionals are used. Consequently, the TD-DFT excess dipole moments are often too small, an error that can be fortuitously corrected for charge-transfer transition by selecting a pure or a hybrid functional containing a small share of exact exchange. This error-cancelation phenomena explains the contradictory conclusions obtained in previous investigations. Overall, the largest correlation between CC2 and TD-DFT excess dipoles is obtained with M06-2X, but at the price of a nearly systematic underestimation of this property by ca. 1 D. For the excess quadrupole moments, the average errors are of the order of 0.2–0.6 D·Å for the set of small aromatic systems treated.

Quasiparticle and excitonic gaps of one-dimensional carbon chains

The charge density of a one-dimensional sp-bonded chain composed of 26 carbon atoms terminated by H with alternating single and triple bonds.

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.