Ultrafast non-adiabatic dynamics of methyl substituted ethylenes: the π3s Rydberg state

Excited state unimolecular reactions of some polyenes exhibit localization of their dynamics at a single ethylenic double bond. Here we present studies of the fundamental photophysical processes in the ethylene unit itself. Combined femtosecond time-resolved photoelectron spectroscopy (TRPES) and ab initio quantum chemical calculations was applied to the study of excited state dynamics in cis-butene, trans-butene, trimethylethylene, and tetramethylethylene, following initial excitation to their respective π3s Rydberg states. The wavelength dependence of the π3s Rydberg state dynamics of tetramethylethylene was investigated in more detail. The π3s Rydberg to ππ(?) valence state decay rate varies greatly with substituent: the 1,2-di- and tri-methyl substituted ethylenes (cis-butene, trans-butene, and trimethylethylene) show an ultrafast decay (～20 fs), whereas the fully methylated tetramethylethylene shows a decay rate of 2 to 4 orders of magnitude slower. These observations are rationalized in terms of topographical trends in the relevant potential energy surfaces, as found from ab initio calculations: (1) the barrier between the π3s state and the ππ? state increases with increasing methylation, and (2) the π3s∕ππ? minimum energy conical intersection displaces monotonically away from the π3s Franck-Condon region with increasing methylation. The use of systematic methylation in combination with TRPES and ab initio computation is emerging as an important tool in discerning the excited state dynamics of unsaturated hydrocarbons.     

Chromophores in spectroscopy: Ab initio studies of localized descriptions of molecular electronic excitations

We review a computationally efficient approach, based jointly on the Random Phase Approximation (RPA) and on localized molecular orbitals, for calculating and analyzing electronic excitations in terms of the nature of the chromophore and its interaction with its molecular surroundings. The method is applied to two typical chromophoric systems using ab initio extended-basis calculations: the non-conjugated but electronically coupled ethylenic double bonds in norbornadiene (NBD, bicyclo[2.2.1]hepta-2,5-diene) and the chirally perturbed carbonyl chromophore inequatorial 4-methyladamantanone (EMAO). The analyses are a posteriori in nature but provide insights into the spectroscopic properties of medium-sized molecules.     

A Quantitative Assessment of „Through-space”︁ and „Through-bond”︁ Interactions. Application to Semi-empirical SCF Models

The scheme of ‘through-space’ and ‘through-bond’ interaction of (semi)localized orbitals, originally proposed by Hoffmann, is reexamined in terms of SCF many-electron treatments. It is shown that the two types of interaction can be characterized by examining the corresponding off-diagonal matrix elements of the Hartree-Fock matrices of the localized or the symmetry adapted localized orbitals and of the partially diagonalized Hartree-Fock matrices referring to ‘precanonical orbitals’. The procedure outlined is applied to three practical examples using the semiempirical many-electron treatments SPINDO, MINDO/2 and CNDO/2:

• a A reassessment of ‘through-space’ and ‘through-bond’ interaction in norbornadiene indicates, that the latter type of interaction is also of importance for the orbital based mainly on the antisymmetric combination of the localized x-orbitals. The differences in the predictions derived from the three models are critically examined.
• b The competition between ‘through-space’ and ‘through-bond’ interaction in the series of bicyclic dienes from norbornadiene to bicyclo[4.2.2]-dcca-7,9-diene and in cyclohexa-1,4-diene, i. e. their dependence on the dihedral angle UI is reexamined. It is found that the rationalization for the orbital crossing near ω = 130° deduccd from PE. spectroscopic data can not be as simple as originally suggested and that the relay’ orbitals responsible for ‘through-bond interaction affecting both the symmetric and the antisymmetric combination of the π-orbitals extend over the whole CC-σ-system of the six membered ring.
• c ‘Through-bond’ interaction of the two lone pair orbitals in 1,4-diazabicyclo[2.2.2]octane is found to be large for their symmetric and the antisymmetric linear combination.

The analysis quoted, draws attention to some of the dangers involved in using semiempirical treatments for the interpretation of PE. data in terms of Koopmans′ theorem, without due caution.     

Gradients of the polarization energy in the effective fragment potential method

The effective fragment potential (EFP) method is an ab initio based polarizable classical method in which the intermolecular interaction parameters are obtained from preparative ab initio calculations on isolated molecules. The polarization energy in the EFP method is modeled with asymmetric anisotropic dipole polarizability tensors located at the centroids of localized bond and lone pair orbitals of the molecules. Analytic expressions for the translational and rotational gradients (forces and torques) of the EFP polarization energy have been derived and implemented. Periodic boundary conditions (the minimum image convention) and switching functions have also been implemented for the polarization energy, as well as for other EFP interaction terms. With these improvements, molecular dynamics simulations can be performed with the EFP method for various chemical systems.     

The concept of bond order

A method for obtaining precise charge densities in defined regions of space from ab initio molecular wavefunctions is employed to place the concept of bond order on a firm theoretical footing. The bond orders obtained for carboncarbon bonds in a range of organic compounds are assessed: those for buta-1,3-diene confirm that it consists of essentially localised double and single bonds.     

Probing weak molecular orbital interactions in non-conjugated diene molecules by means of near-edge X-ray absorption spectroscopy

Carbon and oxygen near-edge X-ray absorption fine structure (NEXAFS) spectra of 1,4-cyclohexadiene, p-benzoquinone, norbornadiene, norbornadienone, and cis-cis-[4,4,2]propella-3,8-diene-11,12-dione were calculated by means of Hartree-Fock and hybrid density functional theory using the static-exchange (STEX) approximation. The NEXAFS spectra are used as a probe to identify weak molecular interactions between the two non-conjugated ethylenic pi* orbitals present in these molecules. We show that the X-ray absorption spectrum of 1,4-cyclohexadiene exhibits some particular spectral structures in the discrete energy region that evidence diene through-bond orbital interaction, whereas absorption peaks are identified in the norbornadiene and norbornadienone spectra that indicate effective through-space orbital interactions. The molecular structure of the cis-cis-[4,4,2]propella-3,8-diene-11,12-dione isomer is such that the indirect through-bond or through-space diene orbital interactions are too weak to be assigned by its C1s NEXAFS spectrum.     

Optical and electrical properties of three-dimensional interlinked gold nanoparticle assemblies

The optical and electrical properties of 11-20 nm thick films composed of approximately 4 nm gold nanoparticles (Au-NPs) interlinked by six organic dithiol or bis-dithiocarbamate derivatives were compared to investigate how these properties depend on the core of the linker molecule (benzene or cyclohexane) and its metal-binding substituents (thiol or dithiocarbamate). Films prepared with the thiol-terminated linker molecules, (1,4-bis(mercaptomethyl)benzene, 1,4-bis(mercaptomethyl)cyclohexane, 1,4-bis(mercaptoacetamido)benzene, and 1,4-bis(mercaptoacetamido)cyclohexane), exhibit thermally activated charge transport. The activation energies lie between 59 and 71 meV. These films show distinct plasmon absorption bands with maxima between 554 and 589 nm. In contrast, the film prepared with 1,4-cyclohexane-bis(dithiocarbamate) has a significantly red-shifted plasmon band ( approximately 626 nm) and a pronounced absorbance in the near infrared. The activation energy for charge transport is only 14 meV. These differences are explained in terms of the formation of a resonant state at the interface due to overlap of the molecular orbital and metal wave function, leading to an apparent increase in NP diameter. The film prepared with 1,4-phenylene-bis(dithiocarbamate) exhibits metallic properties, indicating the full extension of the electron wave function between interlinked NPs. In all cases, the replacement of the benzene ring with a cyclohexane ring in the center of the linker molecule leads to a 1 order of magnitude decrease in conductivity. A linear relationship is obtained when the logarithm of conductivity is plotted as a function of the number of nonconjugated bonds in the linker molecules. This suggests that nonresonant tunneling along the nonconjugated parts of the molecule governs the electron tunneling decay constant (beta(N)(-)(CON)), while the contribution from the conjugated parts of the molecule is weak (corresponding to resonant tunneling). The obtained value for beta(N)(-)(CON) is approximately 1.0 (per non-conjugated bond) and independent of the nanoparticle-binding group. Hence, the molecules can be viewed as consisting of serial connections of electrically insulating (nonconjugated) and conductive (conjugated) parts.     

Ab initio centroid path integral molecular dynamics: application to vibrational dynamics of diatomic molecular systems

An ab initio centroid molecular dynamics (CMD) method is developed by combining the CMD method with the ab initio molecular orbital method. The ab initio CMD method is applied to vibrational dynamics of diatomic molecules, H2 and HF. For the H2 molecule, the temperature dependence of the peak frequency of the vibrational spectral density is investigated. The results are compared with those obtained by the ab initio classical molecular dynamics method and exact quantum mechanical treatment. It is shown that the vibrational frequency obtained from the ab initio CMD approaches the exact first excitation frequency as the temperature lowers. For the HF molecule, the position autocorrelation function is also analyzed in detail. The present CMD method is shown to well reproduce the exact quantum result for the information on the vibrational properties of the system.     

Probing the ruthenium-cumulene bonding interaction: synthesis and spectroscopic studies of vinylidene- and allenylidene-ruthenium complexes supported by tetradentate macrocyclic tertiary amine and comparisons with diphosphine analogues of ruthenium and osmium

The synthesis and spectroscopic properties of trans-[Cl(16-TMC)Ru[double bond]C[double bond]CHR]PF(6) (16-TMC = 1,5,9,13-tetramethyl-1,5,9,13-tetraazacyclohexadecane, R = C(6)H(4)X-4, X = H (1), Cl (2), Me (3), OMe (4); R = CHPh(2) (5)), trans-[Cl(16-TMC)Ru[double bond]C[double bond]C[double bond]C(C(6)H(4)X-4)(2)]PF(6) (X = H (6), Cl (7), Me (8), OMe (9)), and trans-[Cl(dppm)(2)M[double bond]C[double bond]C[double bond]C(C(6)H(4)X-4)(2)]PF(6) (M = Ru, X = H (10), Cl (11), Me (12); M = Os, X = H (13), Cl (14), Me (15)) are described. The crystal structures of 1, 5, 6, and 8 show that the Ru-C(alpha) and C(alpha)-C(beta) distances of the allenylidene complexes fall between those of the vinylidene and acetylide relatives. Two reversible redox couples are observed by cyclic voltammetry for 6-9, with E(1/2) values ranging from -1.19 to -1.42 and 0.49 to 0.70 V vs Cp(2)Fe(+/0), and they are both 0.2-0.3 and 0.1-0.2 V more reducing than those for 10-12 and 13-15, respectively. The UV-vis spectra of the vinylidene complexes 1-4 are dominated by intense high-energy bands at lambda(max) < or = 310 nm (epsilon(max) > or = 10(4) dm(3) mol(-1) cm(-1)), while weak absorptions at lambda(max) > or = 400 nm (epsilon(max) < or = 10(2) dm(3) mol(-1) cm(-1)) are tentatively assigned to d-d transitions. The resonance Raman spectrum of 5 contains a nominal nu(C[double bond]C) stretch mode of the vinylidene ligand at 1629 cm(-1). The electronic absorption spectra of the allenylidene complexes 6-9 exhibit an intense absorption at lambda(max) = 479-513 nm (epsilon(max) = (2-3) x 10(4) dm(3) mol(-1) cm(-1)). Similar electronic absorption bands have been found for 10-12, but the lowest energy dipole-allowed transition is blue-shifted by 1530-1830 cm(-1) for the Os analogues 13-15. Ab initio calculations have been performed on the ground state of trans-[Cl(NH(3))(4)Ru[double bond]C[double bond]C[double bond]CPh(2)](+) at the MP2 level, and imply that the HOMO is not localized purely on the metal center or allenylidene ligand. The absorption band of 6 at lambda(max) = 479 nm has been probed by resonance Raman spectroscopy. Simulations of the absorption band and the resonance Raman intensities show that the nominal nu(C[double bond]C[double bond]C) stretch mode accounts for ca. 50% of the total vibrational reorganization energy, indicating that this absorption band is strongly coupled to the allenylidene moiety. The excited-state reorganization of the allenylidene ligand is accompanied by rearrangement of the Ru[double bond]C and Ru[bond]N (of 16-TMC) fragments, which supports the existence of bonding interaction between the metal and C[double bond]C[double bond]C unit in the electronic excited state.     

Molecular dynamics and metadynamics simulations of [2 + 2] photocycloaddition

Triplet state mechanism of [2 + 2] photocycloaddition forming a cyclobutane ring from two ethylenes is investigated in the context of photocatalysis. High‐level ab initio calculations are combined with ab initio adiabatic molecular dynamics and ab initio metadynamics for rare events modeling. In a photocatalytic scheme, a reactant reaches the triplet state either via intersystem crossing (ISC) or triplet sensitization. The model system adopts a biradical structure, which represents energy intersection with the ground state. The system either completes cyclization or undergoes fragmentation into two olefinic units. The potential and free energy surfaces of the cyclobutane/ethylenes system are mapped with multireference approaches describing possible reaction pathways. To obtain a full picture of a double bond photoreactivity, ab initio adiabatic dynamical calculations were used to estimate reaction yields and to model the effects of excess energy. The potential use of density functional theory based approaches for [2 + 2] photocycloaddition was investigated for future simulations and design of realistic photocatalytic systems. Dynamical aspects of [2 + 2] photocycloaddition via a triplet state manifold are investigated by combining ab initio multireference methods and ab initio molecular dynamics and metadynamics. The reaction pathways are studied for a model system of two ethylenes forming a cyclobutane ring to provide a basis for further studies on design of photocatalytic systems.     

Magnitudes and chemical consequences of R(3)N(+)-C-H...O[double bond]C hydrogen bonding

The magnitude of the stabilizing interaction between an aliphatic C[bond]H bond attached to an ammonium nitrogen and a carbonyl oxygen was evaluated by ab initio calculations at the MP2/6-311++G** level of theory. Attractive R(3)N(+)-C-H...O[double bond]C interactions play an important role in supramolecular recognition and various types of stereoselective catalysis. Our calculations show that R(3)N(+)-C-H...O[double bond]C is the strongest hydrogen bond of the C-H...O type known to date. Such hydrogen bonds remain as stabilizing interactions even in water for amide acceptors.     

Ab initio molecular dynamics of excited-state intramolecular proton transfer around a three-state conical intersection in malonaldehyde

Excited-state potential energy surface (PES) characterization is carried out at the CASSCF and MRSDCI levels, followed by ab initio dynamics simulation of excited-state intramolecular proton transfer (ESIPT) on the S2(pipi*) state in malonaldehyde. The proton-transfer transition state lies close to an S2/S1 conical intersection, leading to substantial coupling of proton transfer with electronic relaxation. Proton exchange proceeds freely on S2, but its duration is limited by competition with twisting out of the molecular plane. This rotamerization pathway leads to an intersection of the three lowest singlet states, providing the first detailed report of ab initio dynamics around a three-state intersection (3SI). There is a significant energy barrier to ESIPT on S1, and further pyramidalization of the twisted structure leads to the minimal energy S1/S0 intersection and energetic terminal point of excited-state dynamics. Kinetics and additional mechanistic details of these pathways are discussed. Significant depletion of the spectroscopic state and recovery of the ground state is seen within the first 250 fs after photoexcitation.     

Auger spectrum of a water molecule after single and double core ionization

The high intensity of free electron lasers opens up the possibility to perform single-shot molecule scattering experiments. However, even for small molecules, radiation damage induced by absorption of high intense x-ray radiation is not yet fully understood. One of the striking effects which occurs under intense x-ray illumination is the creation of double core ionized molecules in considerable quantity. To provide insight into this process, we have studied the dynamics of water molecules in single and double core ionized states by means of electronic transition rate calculations and ab initio molecular dynamics (MD) simulations. From the MD trajectories, photoionization and Auger transition rates were computed based on electronic continuum wavefunctions obtained by explicit integration of the coupled radial Schro?dinger equations. These rates served to solve the master equations for the populations of the relevant electronic states. To account for the nuclear dynamics during the core hole lifetime, the calculated electron emission spectra for different molecular geometries were incoherently accumulated according to the obtained time-dependent populations, thus neglecting possible interference effects between different decay pathways. We find that, in contrast to the single core ionized water molecule, the nuclear dynamics for the double core ionized water molecule during the core hole lifetime leaves a clear fingerprint in the resulting electron emission spectra. The lifetime of the double core ionized water was found to be significantly shorter than half of the single core hole lifetime.     

NONCONJUGATED CONDUCTIVE POLYMERS

Conjugation is not a prerequisite for a polymer to be conductive. A polymer must have at least one double bond in the repeat to become conductive. Interaction with a dopant (e.g., electron acceptor) causes transfer of an electron from the double bond to the dopant creating a hole at the double bond site. Electrical conduction occurs via intersite hopping of holes. Various spectroscopic methods (FTIR, optical absorption, solid-state ¹³C NMR, etc.) along with electrical measurements have been used to elucidate the mechanism of conduction in specific nonconjugated conductive polymers. Examples of these polymers include 1,4-polyisoprene which has one double bond and three single bonds in the repeat. The conductivity of polyisoprene increases 100 billion times upon doping with iodine to a maximum value of 10 S/m. Polyisoprene (natural rubber) is used nonconjugated conductive polymers have a wide range of applications in antistatics, various sensors and optoelectronics.     

The intermolecular potential in NO-N2 and (NO-N2)+ systems: implications for the neutralization of ionic molecular aggregates

The characterization of the non covalent interaction potential, responsible for the intermolecular bond in NO-N(2) and (NO-N(2))(+) molecular aggregates, has been achieved by coupling the predictions of a semiempirical method with the results of a scattering experiment and ab initio calculations. The potential wells for the most stable configurations of the neutral and ionic state, having approximatively a T shape in both cases, fall in the same intermolecular distance range. In addition, in the ionic state, the charge is completely localized on the NO partner. Important implications on the dynamics of the neutralization process, occurring as a vertical transition from ionic to neutral state, are obtained by exploiting the analytical formulation of the interaction and calculating energy spacings and relevant Franck-Condon factors for both intramolecular and intermolecular vibration modes.     

Carbon chains and the (5,5) single-walled nanotube: structure and energetics versus length

Reliable thermochemistry is computed for infinite stretches of pure-carbon materials including acetylenic and cumulenic carbon chains, graphene sheet, and single-walled carbon nanotubes (SWCNTs) by connection to the properties of finite size molecules that grow into the infinitely long systems. Using ab initio G3 theory, the infinite cumulenic chain (:C[double bond]C[double bond]C[double bond]C:) is found to be 1.9+/-0.4 kcal/mol per carbon less stable in free energy at room temperature than the acetylenic chain (.C[triple bond]C-C[triple bond]C.) which is 24.0 kcal/mol less stable than graphite. The difference between carbon-carbon triple, double, and single bond lengths (1.257, 1.279, and 1.333 A, respectively) in infinite chains is evident but much less than with small hydrocarbon molecules. These results are used to evaluate the efficacy of similar calculations with the less rigorous PM3 semiempirical method on the (5,5) SWCNT, which is too large to be studied with high-level ab initio methods. The equilibrium electronic energy change for C(g)-->C[infinite (5,5) SWCNT] is -166.7 kcal/mol, while the corresponding free energy change at room temperature is -153.3 kcal/mol (6.7 kcal/mol less stable than graphite). A threefold alternation (6.866, 6.866, and 6.823 A) in the ring diameter of the equilibrium structure of infinitely long (5,5) SWCNT is apparent, although the stability of this structure over the constant diameter structure is small compared to the zero point energy of the nanotube. In general, different (n,m) SWCNTs have different infinite tube energetics, as well as very different energetic trends that vary significantly with length, diameter, and capping.     

Ab initio structural investigation of methyl and ethyl carbamate, and carbamyl choline

The molecular structures of four conformations of methylcarbamate, three forms of ethylcarbamate, four forms of ethylacetate, and of the trans-form of carbamylcholine, were determined by ab initio gradient geometry refinement on the 4-21G level, and the results are compared with the geometries of homologous systems. Significant changes in bond distances and angles are observed with torsional changes, but, barring long-range non-bonded interactions, they are to a large extent localized in that part of the system which is directly involved with the torsional transition; i.e., through-bond effects in a bond distance chain begin to be insignificant after a sequence of three bonds.     

Molecular electrostatic fields from bond fragments

We apply completely transferable, strictly localized molecular orbitals for the calculation of molecular electrostatic fields. This approach, derived from our previous bond fragment method for the calculation of molecular electrostatic potentials, reduces computational efforts drastically. The fields around small molecules containing first- and second-row atoms are systematically overestimated as compared with ab initio calculations with a minimal STO -3G basis set. However, deviations can be corrected by a simple multiplicative factor, which means that the overall shape of the potential and field around the molecule is correctly reproduced. Our approximate field can be used to determine possible hydration sites around molecules as proposed earlier by Peinel and coworkers. Application of the method is illustrated on the formamide molecule.