Density functional calculations of potential energy surface and charge transfer integrals in molecular triphenylene derivative HAT6

We investigate the effect of structural fluctuations on charge transfer integrals, overlap integrals, and site energies in a system of two stacked molecular 2,3,6,7,10,11-hexakishexyloxytriphenylene (HAT₆), which is a model system for conducting devices in organic photocell applications. A density functional based computational study is reported. Accurate potential energy surface calculations are carried out using an improved meta-hybrid density functional to determine the most stable configuration of the two weakly bound HAT₆ molecules. The equilibrium parameters in terms of the twist angle and co-facial separation are calculated. Adopting the fragment approach within the Kohn–Sham density functional framework, these parameters are combined to a lateral slide, to mimic structural/conformational fluctuations and variations in the columnar phase. The charge transfer and spatial overlap integrals, and site energies, which form the matrix element of the Kohn–Sham Hamiltonian are derived. It is found that these quantities are strongly affected by the conformational variations. The spatial overlap between stacked molecules is found to be of considerable importance since charge transfer integrals obtained using the fragment approach differ significantly from those using the dimer approach.     

Intermolecular potential energy surface of Ar-NO

Rotational spectra of an open-shell complex, Ar-NO, in the electronic ground state have been analyzed by employing an analysis using a free-rotor model, where previously observed data by Mills et al. [J. Phys. Chem. 90, 3331 (1986); 90, 4961 (1986)] and additional transitions observed by Fourier-transform microwave spectroscopy in the present study are simultaneously analyzed with a standard deviation of the least-squares fit to be 27.5 kHz. A two-dimensional intermolecular potential energy surface for Ar-NO has been determined from the analysis. The determined potential energy surface is compared with those of Ar-OH and Ar-SH, which are also complexes containing an open-shell species with the 2Pi ground electronic state.     

The potential energy surface of the triazadiphosphole formation

The formal GaCl₃-assisted [3+2] cycloaddition of two (Me₃Si)₂N–N(SiMe₃)–PCl₂ molecules resulting in the formation of a triazadiphosphole has been studied by means of B3LYP/6-31G(d,p) computations. These calculations revealed a stepwise reaction mechanism starting from the disguised 1,3-dipole and dipolarophile (Me₃Si)₂N–N(SiMe₃)–PCl₂. Comparison of the potential energy surface for the formation of a triazadiphosphole in the presence and without a Lewis acid indicate, that addition of a Lewis acid such as GaCl₃ decreases the activation barriers to Me₃Si–Cl elimination, in accord with experiment. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Nitromethane dimer potential energy surface studies

The interaction energy of several conformations of the nitromethane dimer is investigated at the SCF level. The dispersion energy and counterpoise correction are computed for certain relative orientations of the monomers. Fourth-order many body perturbation theory SDQ-MBPT(4) energies are reported for selected points. Double zeta and double zeta plus polarization basis sets were used. All calculations were done with the monomer fixed at the isolated monomer geometry. Interaction energies as large as 6 kcal/mol are found at minima of hydrogen bonding orientations.     

Transiting the molecular potential energy surface along low energy pathways: The TRREAT algorithm

The Transition Rapidly exploring Random Eigenvector Assisted Tree (TRREAT) algorithm is introduced to perform searches along low curvature pathways on a potential energy surface (PES). The method combines local curvature information about the PES with an iterative Rapidly exploring Random Tree algorithm (LaValle, Computer Science Department, Iowa State University, 1998, TR98–11) that quickly searches high‐dimensional spaces for feasible pathways between local minima. Herein, the method is applied to identifying conformational changes of molecular systems using Cartesian coordinates while avoiding a priori definition of collective variables. We analyze the pathway identification problem for alanine dipeptide, cyclohexane and glycine using nonreactive and reactive forcefields. We show how TRREAT‐identified pathways can be used as valuable input guesses for double‐ended methods such as the Nudged Elastic Band when ascertaining transition state energies. This method can be utilized to improve/extend the reaction databases that lie at the core of automatic chemical reaction mechanism generator software currently developed to build kinetic models of chemical reactions. © 2013 Wiley Periodicals, Inc.     

Fluorobenzene-argon ground-state intermolecular potential energy surface

The ground-state intermolecular potential energy surface for the fluorobenzene-argon van der Waals complex is evaluated using the coupled-cluster singles and doubles including connected triple excitations model, with the augmented correlation consistent polarized valence double-zeta basis set extended with a set of 3s3p2d1f1g midbond functions. In the surface minima the Ar atom is located above and below the fluorobenzene plane at a distance of 3.562 A from the fluorobenzene center of mass and at an angle of 6.33 degrees with respect to the axis perpendicular to the fluorobenzene plane. The corresponding binding energy is 391.1 cm(-1). Both these results and the eigenvalues obtained from the potential compare well with the experimental data available.     

Stationary points on the aminomethanol potential energy surface

Summary In this work we study surface fitting equations for a rigid rotor model of aminomethanol. The energies were obtained from the GAUSSIAN88 package using 3-21G bases and fitted on a least square equation, thus generating a Fourier series expansion of the energy as a function of two dihedral angles. The dihedral angles chosen are those that represent rotation around the C-O and N-C axes in the first case, and rotation around C-O and inversion around the amino group in the second case. Results indicate that the hydroxyl hydrogen is subject to almost free rotation around the C-O axis. Further fully relaxed 6-31G* calculations were performed in order to qualify the results obtained for the rigid rotor model.     

On the OCS2 Singlet potential energy surface

The reaction between atomic oxygen and carbon disulfide is predicted to lead to at least two primary products, which are the dithiiranone ( 1 ) and the oxathiirane-thione ( 2 ) and/or the carbon disulfide S-oxide ( 4 ). The possible intramolecular equilibria 1 ? 2, 1 ? 3, 2 ? 4 , and 2 ? 5 as well as the fragmentations of the possible intermediates 1 – 5 have been studied theoretically within the semiempirical CNDO /B framework as conceivable ground-state reactions. On the basis of MO correlations and potential energy changes along the reaction paths, supplementary with previously reported experimental data, the single molecular transformations and the eventual product formations are discussed.     

Dissociation of toluene cation: a new potential energy surface

The potential energy surface (PES) for the formation of tropylium and benzylium ions from toluene cation (1) has been explored theoretically. Quantum chemical calculations at the B3LYP/6-311++G and G3//B3LYP levels were performed. A pathway to form o-isotoluene (5-methylene-1,3-cyclohexadiene) cation (5) from 1 was found. The isomerization occurs by two consecutive 1,2-H shifts from CH(3) to the ortho position of the aromatic ring via a distonic benzenium cation (2), which is also an intermediate in the well-known isomerization of 1 to cycloheptatriene cation (4). Since the barrier for the formation of 2 is the highest in the two isomerization pathways, 1, 4, and 5 are interconvertible energetically prior to dissociation. The benzylium ion can be produced via 5 as well as from 1 and the tropylium ion via 4. Rice-Ramsperger-Kassel-Marcus model calculations were carried out based on the obtained PES. The result agrees with previous experimental observations. From a theoretical analysis of kinetics of the isomerizations and dissociations, we suggest that 5 plays an important role in the formation of C(7)H(7)(+) from 1.     

Excited-state potential energy surface for the photophysics of adenine

The decay paths on the singlet excited-state surface of 9H-adenine and the associated energy barriers have been calculated at the CAS-PT2//CASSCF level. There are three fundamental paths for the photophysics: two paths for the (1)L(b) state which are virtually barrierless at the present level of theory and correspond to formation of the (n,pi) intermediate and direct decay to the ground state and a third path for ground-state decay of the (n,pi) state with an activation barrier of approximately 0.1 eV. The (1)L(a) state, which has the largest oscillator strength, either decays directly to the ground state or contributes indirectly to the excited-state lifetime by populating the two other states. The results are used to interpret the photophysics in terms of an excited-state plateau for the (1)L(b) state that corresponds to the short-lived excited-state component (approximately 0.1 ps) and a well (i.e., a proper minimum) for the (n,pi) state that gives rise to the long component (1 ps or more). The direct decay to the ground state of the (1)L(b) state is probably the decay channel invoked to explain the experimental wavelength dependence of the relative amplitudes of the two components. In addition to that, the excited-state component in the nanosecond range detected in the time-resolved photoelectron spectrum is proposed to be a triplet (pi,pi) state formed after intersystem crossing from the singlet (n,pi) state.     

Analysis of critical points on the potential energy surface

A simple classification scheme is proposed for critical points, based only on rankr and signatures of the (n,n)-matrixG of harmonic force constants. The determination ofr ands, e.g. by the well-known factorizationG=L ^T gL (L: triangular matrix,g: diagonal matrix), has several theoretical as well as practical (computational) advantages over the inspection of eigenvalues ofG, so far used in quantum chemistry. The eigenvalues are sufficient butnot necessary for a classification whereas rank and signature are the only necessary and sufficient prerequisites for solving the task. For the purpose of presenting a working example, by calculating only a 2×2 torque constant matrix, it is shown that the coplanar ethylbenzene is unstable in the CNDO/2 picture.     

Synthesis and characterization of naphthalene-substituted triphenylene discotic liquid crystals

A series of naphthalene-substituted triphenylene liquid crystals (LCs), viz. triphenylene-2,3,6,7,10,11-hexayl hexakis(6-alkenyloxy-2-naphthoate)s (HTPnN compounds) were synthesized and characterized. Alkenyloxy groups containing three to eleven carbon atoms were used as peripheral spacers in these liquid crystals. The discotic liquid crystals synthesized trap one to three water molecules depending on the length of their peripheral spacers. Differential scanning calorimetry, polarizing optical microscopy and X-ray diffraction measurements confirm both nematic (N_d) and rectangular disordered columnar (Col_rd) phases for most of the synthesized discotic LCs. The clearing temperatures of these discotic LCs increase with increasing peripheral spacer length. A hybrid aligned HTP9N sample was investigated to evaluate its optical performance. Retardation values of this hybrid sample decrease with increasing wavelength and increase with increasing measuring angle.     

Anchoring the potential energy surface of the cyclic water trimer

Six cyclic stationary points on the water trimer potential energy surface have been fully optimized at the MP2 level with the aug-cc-pVQZ basis set. In agreement with previous work, harmonic vibrational frequencies indicate that two structures are minima, three are transition states connecting minima on the surface while the remaining stationary point is a higher-order saddle point. The 1- and n-particle limits of the electronic energies of each of these six structures were estimated by systematically varying both the basis sets and theoretical methods. The former limit was approached with the cc-pVXZ and aug-cc-pVXZ families of basis sets (X=2-7) while MP2, CCSD(T), and BD(TQ) calculations helped examine the latter. Core correlation effects have also been assessed at the MP2 level with the cc-pCVXZ series of basis sets (X=2-5). These data have been combined to provide highly accurate relative energies and dissociation energies for these stationary points.     

Surveying a potential energy surface by eigenvector-following

We have developed a method to search potential energy surfaces which avoids some of the difficulties associated with trapping in local minima. Steps are directly taken between minima using eigenvector-following. Exploration of this space by low temperature Metropolis Monte Carlo is a useful global optimisation tool. This method successfully finds the lowest energy icosahedral minima of Lennard- Jones clusters from random starting configurations, but cannot find the global minimum in a reasonable time for difficult cases such as the 38-atom Lennard-Jones cluster where the face-centred-cubic truncated octahedron is lowest in energy. However, by performing searches at higher temperatures, we have found a pathway between the truncated octahedron and the lowest energy icosahedral minima. Such a pathway may be illustrative of some of the structural transformations that are observed for supported metal clusters by electron microscopy.     

p-Difluorobenzene-argon ground state intermolecular potential energy surface

The ground state intermolecular potential energy surface for the p-difluorobenzene-Ar van der Waals complex is evaluated using the coupled cluster singles and doubles including connected triple excitations [CCSD(T)] model and the augmented correlation consistent polarized valence double-zeta basis set extended with a set of 3s3p2d1f1g midbond functions. The surface minima are characterized by the Ar atom located above and below the difluorobenzene center of mass at a distance of 3.5290 A. The corresponding binding energy is -398.856 cm(-1). The surface is used in the evaluation of the intermolecular level structure of the complex. The results clearly improve previously available data and show the importance of using a good correlation method and basis set when dealing with van der Waals complexes.     

Theoretical study on the potential energy surface and vibrational energy levels of HXeI

The potential energy surface for the electronic ground state of the HXeI molecule is constructed by using the internally contracted multi-reference configuration interaction with the Davidson correction（icMRCI＋Q）method and large basis sets. The stabilities and dissociation barriers are identified from the potential energy surfaces.The three-body dissociation channel is found to be the dominate dissociation channel for HXeI.Based on the obtained potentials,vibrational energy levels of HXeI are calculated using the Lanczos algorithm.Our theoretical results are in excellent agreement with the available observed values.     

The anharmonic potential energy surface of methyl fluoride

A large experimental spectroscopic data set sensitive to the cubic anharmonic potential energy surface (PES) of methyl fluoride has been compiled from the literature for six symmetric and asymmetric top isotopomers of methyl fluoride: 12CH3F, 13CH3FF, 12CD3F, 13CD3F, 12CHD2F and 12CH2DF. This empirical data set has been used to critically assess the best available literature ab initio cubic force field and various 'improved' theoretical force fields. A perturbation-resonance approach to the calculation of the observables from the force constants has been utilized and existing PESs were found to reproduce the data poorly. The careful treatment required for the correct theoretical reproduction of empirical spectroscopic constants arising from highly correlated least-squares fits to the original data is discussed. A new fit to the data has been performed (optimising 19 of the 38 cubic force constants) resulting in a much improved PES. The latter has been used to predict currently unknown spectroscopic constants for the least well-characterised isotopomer 12CH2DF. The prospects for a future empirical determination of the complete cubic force field of methyl fluoride are discussed and new data most likely to yield new information on the PES identified.     

Ab initio intermolecular potential energy surface of He-LiH

The intermolecular potential energy surface of He-LiH complex was studied using the full-electronic complete forth-order Miller-Plesset perturbation (MPPT) method.In ab initio calculations,the bond length of LiH was fixed at 0 159 5 nm.The potential has two local minima of Vm=-179.93 cm for the linear He LiH geormetrv at Rm=0.227 nm and Vm=-10.44 cm-1 for the linear He-HL1 geometry at Rm=0.516 nm The potemal exhibits strong anisotropy The analytic potential function with 31 parameters was determined by fitting to the calculated ab,mtio potentials The influence of variation of LiH bond length on the potential energy surface was also studied