Singlet and triplet molecular interaction in excimers

A perturbative approach is developed to determine, term by term, the contributions of the various forces to the excimer potentials of the singlet and triplet excimers. The results show that the singlet excimer of naphthalene is more stable than the corresponding triplet excimer primarily due to large contributions of the exciton-resonance and the dispersion energy terms. The variation of the various terms with the conformations of the excimers suggests that the singlet and triplet excimers of naphthalene cannot have identical structure.

Study of the dynamics of protein folding through minimalistic models

Starting with the Levinthal paradox, a brief introduction to the protein folding problem is presented. The existing theories of protein folding, including the folding funnel scenario, are discussed. After briefly discussing different simulation studies of model proteins, we discuss our recent work on the dynamics of folding of the model HP-36 (the chicken villin headpiece) protein by using a simplified hydropathy scale. Special attention has been paid to the statics and dynamics of contact formation among the hydrophobic residues. The results obtained from this simple model appear to be surprisingly similar to several features observed in the folding of real proteins. The account concludes with a discussion of future problems.

Partition function for hindered, one-and multi-dimensional rotors

Writing the hindered rotor (hr) partition function as the trace of (rho) over cap = e(-beta(H) over cap hr), we approximate it by the sum of contributions from a set of points in position space. The contribution of the density matrix from each point is approximated by performing a local harmonic expansion around it. The highlight of this method is that it can be easily extended to multidimensional systems. Local harmonic expansion leads to a breakdown of the method a low temperatures. In order to calculate the partition function at low temperatures, we suggest a matrix multiplication procedure. The results obtained using these methods closely agree with the exact partition function at all temperature ranges. Our method bypasses the evaluation of eigenvalues and eigenfunctions and evaluates the density matrix for internal rotation directly. We also suggest a procedure to account for the antisymmetry of the total wavefunction in the same. (C) 2012 Elsevier B.V. All rights reserved.

Study of structures, energies and vibrational frequencies of (O-2)(n)(+) (n=2-5) clusters by GGA and meta-GGA density functional methods

Using Generalized Gradient Approximation (GGA) and meta-GGA density functional methods, structures, binding energies and harmonic vibrational frequencies for the clusters O-4(+), O-6(+), O-8(+) and O-10(+) have been calculated. The stable structures of O-4(+), O-6(+), O-8(+) and O-10(+) have point groups D-2h, D-3h, D-4h, and D-5h optimized on the quartet, sextet, octet and dectet potential energy surfaces, respectively. Rectangular (D-2h) O-4(+) has been found to be more stable compared to trans-planar (C-2h) on the quartet potential energy surface. Cyclic structure (D-3h) of CA cluster ion has been calculated to be more stable than other structures. Binding energy (B.E.) of the cyclic O-6(+) is in good agreement with experimental measurement. The zero-point corrected B.E. of O-8(+) with D4h symmetry on the octet potential energy surface and zero-point corrected B.E. of O-10(+) with D-5h symmetry on the dectet potential energy surface are also in good agreement with experimental values. The B.E. value for O-4(+) is close to the experimental value when single point energy is calculated by Brueckner coupled-cluster method, BD(T). (C) 2014 Elsevier B.V. All rights reserved.

Estimating successive pK(a) values of polyprotic acids from ab initio molecular dynamics using metadynamics: the dissociation of phthalic acid and its isomers

Estimation of the dissociation constant, or pK(a), of weak acids continues to be a central goal in theoretical chemistry. Here we show that ab initio Car-Parrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free energy profile of the dissociation reaction can provide reasonable estimates of the successive pK(a) values of polyprotic acids. We use the distance-dependent coordination number of the protons bound to the hydroxyl oxygen of the carboxylic group as the collective variable to explore the free energy profile of the dissociation process. Water molecules, sufficient to complete three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. Two distinct minima corresponding to the dissociated and un-dissociated states of the acid are observed and the difference in their free energy values provides the estimate for pK(a), the acid dissociation constant. We show that the method predicts the pK(a) value of benzoic acid in good agreement with experiment and then show using phthalic acid (benzene dicarboxylic acid) as a test system that both the first and second pK(a) values as well, as the subtle difference in their values for different isomers can be predicted in reasonable agreement with experimental data.

Spin populations in one-dimensional alternant spin systems: A theoretical approach

Spin-density maps, deduced from polarized neutron diffraction experiments, for both the pair and chain compounds of the system Mn2+Cu2+ have been reported recently. These results have motivated us to investigate theoretically the spin populations in such alternant mixed-spin systems. In this paper, we report our studies on the one-dimensional ferrimagnetic systems (S-A,S-B)(N) where hi is the number of AB pairs. We have considered all cases in which the spin Sri takes on allowed values in the range I to 7/2 while the spin S-B is held fixed at 1/2. The theoretical studies have been carried out on the isotropic Heisenberg model, using the density matrix renormalization group method. The effect of the magnitude of the larger spin SA On the quantum fluctuations in both A and B sublattices has been studied as a function of the system size N. We have investigated systems with both periodic and open boundary conditions, the latter with a view to understanding end-of-chain effects. The spin populations have been followed as a function of temperature as well as an applied magnetic field. High-magnetic fields are found to lead to interesting re-entrant behavior. The ratio of spin populations P-A-P-B is not sensitive to temperature at low temperatures.

Theoretical evidence and chemical origin of the magnetism-dependent electrostructural coupling in La2NiMnO6

Using first-principles density-functional calculations, we determine and analyze the Born effective charges Z(*) that describe the coupling between electric field and atomic displacements for ferromagnetic double-perovskite compound, La2NiMnO6. We find that th Born effective charge matrix of Ni in La2NiMnO6, has an anomalously large antisymmetric component, whose magnitude reduces substantially upon change in the magnetic ordering between Ni and Mn, showing it to be a magnetism-dependent electrostructural coupling. We use a local picture of the electronic structure obtained with Wannier functions, along with its band-by-band decomposition to determine its electronic origin.

Large carrier mobilities in octathio[8]circulene crystals: a theoretical study

The electron and hole mobilities of octathio[8]circulene (sulflower) crystal have been calculated using quantum chemical methods, with accurate determination of reorganization energies and the rate of charge transfer, the key parameters controlling the charge carriers conductance. We find this molecular crystal to be an excellent conductor with large mobilities for both the charge carriers. Moreover, the hole mobility is found to be slightly larger than the electron mobility. Such an ambipolar organic crystal with substantial carrier mobilities shows possibilities of sophisticated device fabrication in advanced electronics.

A Theoretical Investigation of the Ni(II)-Catalyzed Hydrovinylation of Styrene

We report a detailed and full computational investigation on the hydrovinylation reaction of styrene with the Ni(II)-phospholane catalytic system, which was originally presumed to proceed through a cationic mechanism involving a nickel hydride intermediate. The following general features emerge from this study on a specific catalyst complex that was found to give quantitative yield and moderate selectivity: (a) the activation barrier for the initiation (18.8 kcal/mol) is higher than that for the reaction due to a low-lying square-planar pentenyl chelate intermediate originating from a Ni(II)-allyl catalyst precursor. Consequently there is an induction period for the catalysis; (b) the exit of product from the catalyst is via a β-H-transfer step instead of the usual β-H elimination pathway, which has a very high activation energy due to a trans effect of the phospholane ligand; (c) the turnover-limiting and enantio- determining transition state is also the β-H-transfer; (d) because of the absence of a hydride intermediate, the unwanted isomerization of the product is prevented; (e) since the enantio-discrimination is decided at the H-transfer stage itself, the configuration of the product in a catalytic cycle influences the enantioselectivity in the subsequent cycle; (f) the trans effect of the sole strong ligand in the d8 square-planar Ni(II), the stability of the η3-benzyl intermediate, and the availability of three coordination sites enable regioselective hydrovinylation over the possible oligomerization/polymerization of the olefin substrates and linear hydrovinylation. This work has also confirmed the previously recognized role of the hemilabile group at various stages in the mechanism.

Theoretical studies of lithium bonding in lithium chloride/aliphatic amine complexes

In the systematic study of amine … LiCl [amines = NH3, CH3NH2, (CH3)2NH] complexes the possibility of an ion-pair structure and the effect of methylation on the stabilization energy is investigated. ΔEis evaluated by the SCF/4-31G method and augmented by the approximate dispersion energy calculated perturbationally. The interaction energy decreases with the increasing number of methyl groups in the amine. The dispersion energy plays a negligible role in the stabilization of complexes. None of the systems studied are ion pairs; their Li bonds are of a so-called molecular type. Due to the divergence of the multipole expansion, the attempt to correct the 4-31G stabilization energies via the electrostatic energy fails. The relative order of the ΔE in the series of complexes is verified instead in the extended basis set calculation. The lithium bonds are compared with their H-bonded analogues.

Which One Is Preferred: Myers-Saito Cyclization of Ene-Yne-Allene or Garratt-Braverman Cyclization of Conjugated Bisallenic Sulfone? A Theoretical and Experimental Study

A competitive scenario between Myers-Saito (MS) and Garraff-Braverman (GB) cyclization has been created in a molecule. High-level computations indicate a preference for GB over MS cyclization. The activation energies for the rate-determining steps of the GB and MS cyclizations were found to be the same (24.4 kcal/mol) at the B3LYP/6-31G* level of theory; thus, from the kinetic point of view, both reactions are feasible. However, the main biradical intermediate GB2 of the GB reaction is 6.2 kcal/mol lower in energy than the biradical MS2, which is the main intermediate of MS reaction, so GB cyclization is thermodynamically favored over MS cyclization. To verify the prediction by computational techniques, bisenediynyl sulfones 1-4 and bisenediynyl sulfoxide 17 were synthesized. Under basic conditions, these molecules isomerized to a system possessing both the ene-yne-allene and the bisallenic sulfone. The isolation of only one product, identified as the corresponding naphthalene- or benzene-fused sulfone 8-11, indicated the occurrence of GB cyclization as the sole reaction pathway. No product corresponding to the MS cyclization pathway could be isolated. Though the theoretical prediction showed a preference for the GB pathway over the MS pathway, the exclusive preference for GB over MS cyclization is very striking. Further analysis showed that the intramolecular self-quenching nature of the GB pathway may play an important role in the complete preference for this reaction. Apart from the mechanistic studies, these sulfones showed DNA cleavage activity that had an inverse relation with the reactivity order. Our findings are important for the design of artificial DNA-cleaving agents.

Theoretical Investigation on the Effect of Different Nitrogen Donors on Intramolecular Se center dot center dot center dot N Interactions

The effect of different donor nitrogen atoms on the strength and nature of intramolecular Se center dot center dot center dot N interactions is evaluated for organoselenium compounds having N,N-dimethylaminomethyl (dime), oxazoline (oxa) and pyridyl (py) substituents. Quantum chemical calculations on three series of compounds [2-(dime)C6H4SeX (1a-g), 2-(oxa)C6H4SeX (2a-g), 2- (py)C6H4SeX (3-ag); X=Cl, Br, OH, CN, SPh, SePh, CH3] at the B3LYP/6-31G(d) level show that the stability of different conformers depends on the strength of intramolecular nonbonded Se center dot center dot center dot N interactions. Natural bond orbital (NBO), NBO deletion and atoms in molecules (AIM) analyses suggest that the nature of the Se center dot center dot center dot N interaction is predominantly covalent and involves nN ->sigma*(Se-X) orbital interaction. In the three series of compounds, the strength of the Se center dot center dot center dot N interaction decreases in the order 3>2>1 for a particular X, and it decreases in the order Cl > Br > OH>SPh approximate to CN approximate to SePh>CH3 for all the three series 1-3. However, further analyses suggest that the differences in strength of Se center dot center dot center dot N interaction in 1-3 is predominantly determined by the distance between the Se and N atoms, which in turn is an outcome of specific structures of 1, 2 and 3, and the nature of the donor nitrogen atoms involved has very little effect on the strength of Se center dot center dot center dot N interaction. It is also observed that Se center dot center dot center dot N interaction becomes stronger in polar solvents such as CHCl3, as indicated by the shorter r(Se center dot center dot center dot N) and higher E-Se center dot center dot center dot N values in CHCl3 compared to those observed in the gas phase.

Investigation of inter-ion interactions in N,N,N',N'-tetramethylethylenediammonium dithiocyanate via experimental and theoretical charge density studies

The crystal structure of the N,N,N',N'-tetramethylethylenediammonium dithiocyanate salt has been examined by experimental charge density studies from high-resolution X-ray diffraction data. The corresponding results are compared with multipole refinements, using theoretical structure factors obtained from a periodic density functional theory calculation at the B3LYP level with a 6-31G** basis set. The salt crystallizes in space group P (1) over bar and contains only a single ion pair with an inversion center in the cation. The salt has thus one unique classical N+-H center dot center dot center dot(NCS)(-) hydrogen bond but also has six other weaker interactions: four C-H center dot center dot center dot S, one C-H center dot center dot center dot N, and one C-H center dot center dot center dot C-pi. The nature of all these interactions has been examined topologically using Bader's quantum theory of "atoms in molecules" and all eight of the Koch-Popelier criteria. The experimental and theoretical approaches agree well and both show that the inter-ion interactions, even in this simplest of systems, play an integrated and complex role in the packing of the ions in the crystal. Electrostatic potential maps are derived from experimental charge densities. This is the first time such a system has been examined in detail by these methods.