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.

Diagrammatic valence bond studies on hemocyanin

The Diagrammatic Valence Bond studies on the active sites of hemocyanin, consisting of two Cu(I) ions and an oxygen molecule, are performed to find out the stable geometrical pattern and electronic structure. Different parameters used in this theoretical approach are taken from existing literature on high T-c superconductors. Attempts have been made to find out the differences in electronic structure of [Cu2O2](+2) and [Cu2O2N4](+2) as it is observed that coordination of nitrogen ligand do affect electronic structure i.e. spin excitation gaps and charge and spin density distribution. A comparison of our results with earlier theoretical results are also presented.

Theoretical analysis of x-ray-absorption near-edge fine structure at the O and metal K edges of LaFeO3 and LaCoO3

We present experimental x-ray-absorption spectra at the oxygen and 3d transition-metal K edges of LaFeO3 and LaCoO3. We interpret the experimental results in terms of detailed theoretical calculations based on multiple-scattering theory. Along with providing an understanding of the origin of various experimental features, we investigate the effects of structural distortions and the core-hole potential in determining the experimental spectral shape. The results indicate that the core-hole potential as well as many-body effects within the valence electrons do not have any strong effect on the spectra suggesting that the spectral features can be directly interpreted in terms of the electronic structure of such compounds.

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.

He II spectra of La, Ce and Yb: Novel features in the valence band region

He II photoelectron spectra of La, Ce and Yb show features which cannot be explained in terms of single electron excitations. It is proposed that these are due to formation of electron-hole paris.

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.

3d and 4d core level X-ray photoelectron spectroscopy of mixed-valence CePd3

3d and 4d core-level XPS spectra for CePd3, a mixed-valence system, have been measured. Each spectrum exhibits two sets of structures, each corresponding to one of the valence states of cerium. Thus the usefulness of XPS, which has so far not been used extensively to investigate the mixed-valence cerium systems, is pointed out.

Valence bond analysis of the Kondo chain Hamiltonian

Accurate extrapolations for the ground state energy per site of the one - dimensional Kondo chain system is obtained from exact finite system calculations carried out employing a valence bond scheme. An analysis of the ground state wave function indicates that the localized spin is quenched for all nonzero values of the Kondo exchange constant in one dimension.