Breakdown of the mirror image symmetry in the optical absorption/emission spectra of oligo(para-phenylene)s

The absorption and emission spectra of most luminescent, pi-conjugated, organic molecules are the mirror image of each other. In some cases, however, this symmetry is severely broken. In the present work, the asymmetry between the absorption and fluorescence spectra in molecular systems consisting of para-linked phenyl rings is studied. The vibronic structure of the emission and absorption bands is calculated from ab-initio quantum chemical methods and a subsequent, rigorous Franck-Condon treatment. Good agreement with experiment is achieved. A clear relation can be established between the strongly anharmonic double-well potential for the phenylene ring librations around the long molecular axis and the observed deviation from the mirror image symmetry. Consequences for related compounds and temperature dependent optical measurements are also discussed.Comment: 12 pages, 13 Figure

Adsorption of Cu, Ag, and Au atoms on graphene including van der Waals interactions

We performed a systematic density functional study of the adsorption of copper, silver, and gold adatoms on graphene, especially accounting for van der Waals interactions by the vdW-DF and the PBE+D2 methods. In particular, we analyze the preferred adsorption site (among top, bridge, and hollow positions) together with the corresponding distortion of the graphene sheet and identify diffusion paths. Both vdW schemes show that the coinage metal atoms do bind to the graphene sheet and that in some cases the buckling of the graphene can be significant. The results for silver are at variance with those obtained with GGA, which gives no binding in this case. However, we observe some quantitative differences between the vdW-DF and the PBE+D2 methods. For instance the adsorption energies calculated with the PBE+D2 method are systematically higher than the ones obtained with vdW-DF. Moreover, the equilibrium distances computed with PBE+D2 are shorter than those calculated with the vdW-DF method

Reversible Metal-Free Carbon Dioxide Binding by Frustrated Lewis Pairs

Frustrated Lewis pairs comprising phosphine and borane react to reversibly bind and release CO2, offering a rare example of metal-free CO2 sequestration. The mechanism of formation of CO2 derivatives by almost simultaneous P-C and O-B bond formation was characterized by quantum chemical calculations

Performance of ab initio and density functional methods for conformational equilibria of CnH2n+2 alkane isomers (n=2-8)

Conformational energies of n-butane, n-pentane, and n-hexane have been calculated at the CCSD(T) level and at or near the basis set limit. Post-CCSD(T) contribution were considered and found to be unimportant. The data thus obtained were used to assess the performance of a variety of density functional methods. Double-hybrid functionals like B2GP-PLYP and B2K-PLYP, especially with a small Grimme-type empirical dispersion correction, are capable of rendering conformational energies of CCSD(T) quality. These were then used as a `secondary standard' for a larger sample of alkanes, including isopentane and the branched hexanes as well as key isomers of heptane and octane. Popular DFT functionals like B3LYP, B3PW91, BLYP, PBE, and PBE0 tend to overestimate conformer energies without dispersion correction, while the M06 family severely underestimates GG interaction energies. Grimme-type dispersion corrections for these overcorrect and lead to qualitatively wrong conformer orderings. All of these functionals also exhibit deficiencies in the conformer geometries, particularly the backbone torsion angles. The PW6B95 and, to a lesser extent, BMK functionals are relatively free of these deficiencies. Performance of these methods is further investigated to derive conformer ensemble corrections to the enthalpy function, $H_{298}-H_0$, and the Gibbs energy function, ${\rm gef}(T)\equiv - [G(T)-H_0]/T$, of these alkanes. While $H_{298}-H_0$ is only moderately sensitive to the level of theory, ${\rm gef}(T)$ exhibits more pronounced sensitivity. Once again, double hybrids acquit themselves very well.Comment: J. Phys. Chem. A, revised [Walter Thiel festschrift

Elastic and vibrational properties of alpha and beta-PbO

The structure, electronic and dynamic properties of the two layered alpha (litharge) and beta (massicot) phases of PbO have been studied by density functional methods. The role of London dispersion interactions as leading component of the total interaction energy between layers has been addressed by using the Grimme's approach, in which new parameters for Pb and O atoms have been developed. Both gradient corrected and hybrid functionals have been adopted using Gaussian-type basis sets of polarized triple zeta quality for O atoms and small core pseudo-potential for the Pb atoms. Basis set superposition error (BSSE) has been accounted for by the Boys-Bernardi correction to compute the interlayer separation. Cross check with calculations adopting plane waves that are BSSE free have also been performed for both structures and vibrational frequencies. With the new set of proposed Grimme's type parameters structures and dynamical parameters for both PbO phases are in good agreement with experimental data.Comment: 8 pages, 5 figure

Accurate Treatment of Large Supramolecular Complexes by Double-Hybrid Density Functionals Coupled with Nonlocal van der Waals Corrections

In this work, we present a thorough assessment of the performance of some representative double-hybrid density functionals (revPBE0-DH-NL and B2PLYP-NL) as well as their parent hybrid and GGA counterparts, in combination with the most modern version of the nonlocal (NL) van der Waals correction to describe very large weakly interacting molecular systems dominated by noncovalent interactions. Prior to the assessment, an accurate and homogeneous set of reference interaction energies was computed for the supramolecular complexes constituting the L7 and S12L data sets by using the novel, precise, and efficient DLPNO-CCSD(T) method at the complete basis set limit (CBS). The correction of the basis set superposition error and the inclusion of the deformation energies (for the S12L set) have been crucial for obtaining precise DLPNO-CCSD(T)/CBS interaction energies. Among the density functionals evaluated, the double-hybrid revPBE0-DH-NL and B2PLYP-NL with the three-body dispersion correction provide remarkably accurate association energies very close to the chemical accuracy. Overall, the NL van der Waals approach combined with proper density functionals can be seen as an accurate and affordable computational tool for the modeling of large weakly bonded supramolecular systems.Financial support by the “Ministerio de Economía y Competitividad” (MINECO) of Spain and European FEDER funds through projects CTQ2011-27253 and CTQ2012-31914 is acknowledged. The support of the Generalitat Valenciana (Prometeo/2012/053) is also acknowledged. J.A. thanks the EU for the FP7-PEOPLE-2012-IEF-329513 grant. J.C. acknowledges the “Ministerio de Educación, Cultura y Deporte” (MECD) of Spain for a predoctoral FPU grant

Spin states of zigzag-edged Mobius graphene nanoribbons from first principles

Mobius graphene nanoribbons have only one edge topologically. How the magnetic structures, previously associated with the two edges of zigzag-edged flat nanoribbons or cyclic nanorings, would change for their Mobius counterparts is an intriguing question. Using spin-polarized density functional theory, we shed light on this question. We examine spin states of zigzag-edged Mobius graphene nanoribbons (ZMGNRs) with different widths and lengths. We find a triplet ground state for a Mobius cyclacene, while the corresponding two-edged cyclacene has an open-shell singlet ground state. For wider ZMGNRs, the total magnetization of the ground state is found to increase with the ribbon length. For example, a quintet ground state is found for a ZMGNR. Local magnetic moments on the edge carbon atoms form domains of majority and minor spins along the edge. Spins at the domain boundaries are found to be frustrated. Our findings show that the Mobius topology (i.e., only one edge) causes ZMGNRs to favor one spin over the other, leading to a ground state with non-zero total magnetization.Comment: 17 pages, 4 figure