Excitation and relaxation energies of trans-stilbene: Confined singlet, triplet, and charged bipolarons

The π-electronic excitations and excited-state geometries of trans-stilbene (tS) are found by combining exact solutions of the Pariser-Parr-Pople (PPP) model and semiempirical Parametric Method 3 (PM3) calculations. Comprehensive comparisons with tS spectra are obtained and related to the fluorescence and topological alternation of poly(paraphenylenevinylene) (PPV). The one-photon absorption and triplet of tS correspond, respectively, to singlet and triplet bipolarons confined to two phenyls, while the tS2- ground state is a confined charged bipolaron. Independent estimates of the relaxation energy between vertical and adiabatic excitation show the bipolaron binding energy to depend on both charge and spin, as expected for interacting π electrons in correlated or molecular states. Complete configuration interaction within the PPP model of tS accounts for the singlet-triplet gap, for the fine-structure constants and triplet-triplet spectra, for two-photon transitions and intensities, and for one-photon spectra and the radiative lifetime, although the relative position of nearly degenerate covalent and ionic singlets is not resolved. The planar PM3 geometry and low rotational barrier of tS agree with resolved rotational and vibrational spectra in molecular beams. PM3 excitation and relaxation energies for tS bipolarons are consistent with experiment and with PPP results. Instead of the exciton model, we interpret tS excitations in terms of states that are localized on each ring or extended over an alternating chain, as found exactly in Hückel theory, and find nearly degenerate transitions between extended and localized states in the singlet, triplet, and dianion manifolds. The large topological alternation of the extended system increases the ionicity and interchanges the order of the lowest one- and two-photon absorption of PPV relative to polyenes.

1H NMR Study of Molecular Motions and Phase Transitions in Methyl Ammonium Hexabromo Selenate [(CH3NH3)2SeBr6]

The temperature dependence of 1H spin-lattice relaxation time, T1, and that of the second moment, M2, are analysed in the temperature range 390 K to 77 K. A plot of T1 vs inverse temperature shows three phase transitions at 250 K, 167 K and 111 K. At 167 K, T1 displays a large jump while it shows changes in slope at 250 K and 111 K. In the high temperature phase (> 167 K), the correlated motion of CH3 and NH3 groups is found to cause the relaxation while their uncorrelated motion takes over in the low temperature phases (< 167 K). The unusual T1 behaviour in phase II (250 K-167 K) is ascribed to the small angle torsion of the cation. A constant M2 value of ∼ 9.7 G2, throughout the range of temperature studied, indicates the presence of reorientation of CH3 and NH3 groups.

Molecular structure of electron donor-acceptor complexes of metallotetraphenylporphyrins with trinitrobenzene

The metallotetraphenylporphyrins, MTPPs, where M=Co(II), Cu(II) and Ag(II) form one to one molecular complexes in solution with 1,3,5-trinitrobenzene (TNB). The crystal structure of CoTPP.TNB.2CH3OH revealed that the mean separation between the porphyrin and TNB planes is 3·27 Å and the centre of the aryl ring of TNB is displaced by 0·90 Å from the centre of the porphyrin plane. Extended Huckel molecular orbital calculations suggest that acceptor orbitals are predominantly nitro-group based rather than an aryl π framework. The contribution of the metal orbitals in the donation is also seen as predicted by the structural data. Electrochemical redox measurements in solution provide evidence for charge-transfer stabilisation in these complexes.

NDDO-based CI methods for the prediction of electronic spectra and sum-over-states molecular hyperpolarization

NDDO-based (AM1) configuration interaction (CI) calculations have been used to calculate the wavelength and oscillator strengths of electronic absorptions in organic molecules and the results used in a sum-over-states treatment to calculate second-order-hyperpolarizabilities. The results for both spectra and hyperpolarizabilities are of acceptable quality as long as a suitable CI-expansion is used. We have found that using an active space of eight electrons in eight orbitals and including all single and pair-double excitations in the CI leads to results that agree well with experiment and that do not change significantly with increasing active space for most organic molecules. Calculated second-order hyperpolarizabilities using this type of CI within a sum-over-states calculation appear to be of useful accuracy.

Solvation dynamics in a Brownian dipolar lattice. Comparison between computer simulation and various molecular theories of solvation dynamics

Several recent theoretical and computer simulation studies have considered solvation dynamics in a Brownian dipolar lattice which provides a simple model solvent for which detailed calculations can be carried out. In this article a fully microscopic calculation of the solvation dynamics of an ion in a Brownian dipolar lattice is presented. The calculation is based on the non‐Markovian molecular hydrodynamic theory developed recently. The main assumption of the present calculation is that the two‐particle orientational correlation functions of the solid can be replaced by those of the liquid state. It is shown that such a calculation provides an excellent agreement with the computer simulation results. More importantly, the present calculations clearly demonstrate that the frequency‐dependent dielectric friction plays an important role in the long time decay of the solvation time correlation function. We also find that the present calculation provides somewhat better agreement than either the dynamic mean spherical approximation (DMSA) or the Fried–Mukamel theory which use the simulated frequency‐dependent dielectric function. It is found that the dissipative kernels used in the molecular hydrodynamic approach and in the Fried–Mukamel theory are vastly different, especially at short times. However, in spite of this disagreement, the two theories still lead to comparable results in good agreement with computer simulation, which suggests that even a semiquantitatively accurate dissipative kernel may be sufficient to obtain a reliable solvation time correlation function. A new wave vector and frequency‐dependent dissipative kernel (or memory function) is proposed which correctly goes over to the appropriate expressions in both the single particle and the collective limits. This form is expected to lead to better results than all the existing descriptions.

Dynamic coupling between the LID and NMP domain motions in the catalytic conversion of ATP and AMP to ADP by adenylate kinase

The catalytic conversion of adenosine triphosphate (ATP) and adenosine monophosphate (AMP) to adenosine diphosphate (ADP) by adenylate kinase (ADK) involves large amplitude, ligand induced domain motions, involving the opening and the closing of ATP binding domain (LID) and AMP binding domain (NMP) domains, during the repeated catalytic cycle. We discover and analyze an interesting dynamical coupling between the motion of the two domains during the opening, using large scale atomistic molecular dynamics trajectory analysis, covariance analysis, and multidimensional free energy calculations with explicit water. Initially, the LID domain must open by a certain amount before the NMP domain can begin to open. Dynamical correlation map shows interesting cross-peak between LID and NMP domain which suggests the presence of correlated motion between them. This is also reflected in our calculated two-dimensional free energy surface contour diagram which has an interesting elliptic shape, revealing a strong correlation between the opening of the LID domain and that of the NMP domain. Our free energy surface of the LID domain motion is rugged due to interaction with water and the signature of ruggedness is evident in the observed root mean square deviation variation and its fluctuation time correlation functions. We develop a correlated dynamical disorder-type theoretical model to explain the observed dynamic coupling between the motion of the two domains in ADK. Our model correctly reproduces several features of the cross-correlation observed in simulations. (C) 2011 American Institute of Physics. doi:10.1063/1.3516588]

Importance of orientational rearrangement during vitrification of hydrocarbons: dependence on molecular shape

On the basis of Monte Carlo calculations of 2,2-dimethylpropane (neopentane), n-pentane, and 2,2-dimethylbutane (neohexane) at several temperatures, thermodynamic properties and radial distribution functions as well as dimerization and bonding energy distribution functions are reported for both liquid and glassy states. Changes in the radial distribution functions on cooling depend on whether the groups are accessible (peripheral) or inaccessible. Peaks in the radial distribution functions corresponding to peripheral groups do not shift to lower distances on cooling and at times display a large increase in the intensity of the first peak. The peaks due to inaccessible groups, on the other hand, shift to lower distances on cooling. The magnitude of the reorientational contribution in determining the resulting structure of the glass is estimated for the different hydrocarbon molecules investigated. The reorientational contribution is highest for neopentane (26%) followed by isopentane (24%), neohexane (22%), and n-pentane (0%). It appears that molecular geometry has an important role in determining the magnitude of the reorientational contribution to the structure of the glass.

Using Correlated Parameters for Improved Ranking of Protein-Protein Docking Decoys

A successful protein-protein docking study culminates in identification of decoys at top ranks with near-native quaternary structures. However, this task remains enigmatic because no generalized scoring functions exist that effectively infer decoys according to the similarity to near-native quaternary structures. Difficulties arise because of the highly irregular nature of the protein surface and the significant variation of the nonbonding and solvation energies based on the chemical composition of the protein-protein interface. In this work, we describe a novel method combining an interface-size filter, a regression model for geometric compatibility (based on two correlated surface and packing parameters), and normalized interaction energy (calculated from correlated nonbonded and solvation energies), to effectively rank decoys from a set of 10,000 decoys. Tests on 30 unbound binary protein-protein complexes show that in 16 cases we can identify at least one decoy in top three ranks having <= 10 angstrom backbone root mean square deviation from true binding geometry. Comparisons with other state-of-art methods confirm the improved ranking power of our method without the use of any experiment-guided restraints, evolutionary information, statistical propensities, or modified interaction energy equations. Tests on 118 less-difficult bound binary protein-protein complexes with <= 35% sequence redundancy at the interface showed that in 77% cases, at least 1 in 10,000 decoys were identified with <= 5 angstrom backbone root mean square deviation from true geometry at first rank. The work will promote the use of new concepts where correlations among parameters provide more robust scoring models. It will facilitate studies involving molecular interactions, including modeling of large macromolecular assemblies and protein structure prediction. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 32: 787-796, 2011.

Sorbate properties of xenon in cloverite: A molecular dynamics study

Molecular dynamics calculations on xenon adsorbed in the cubic cavity of cloverite, a gallophosphate, is presented. Guest-host radial distribution functions, guest-host energy distribution functions, power spectra, and diffusion coefficients for xenon are reported at 397, 494, and 716 K. Results suggest that xenon is highly mobile at 700 K. A shift in the peak position of the power spectra toward lower frequency is observed with increase in temperature. (C) 1994 Academic Press, Inc.

Correlated electronic states of C-60 fragments

We have studied electronic states of various fragments of C-60 within the Pariser-Parr-Pople (PPP) model and have obtained structural, magnetic and spectral properties of these molecules. The fragments studied include corannulene, fluoranthene and pyracylene. Pyracylene is studied using the exact valence bond (VB) approach while fluoranthene and corannulene are studied using a novel restricted CI technique which employs molecular orbitals for constructing the VB functions. The electronic excitations, bond order and ring currents are calculated for these systems. From these studies, the wide range of absorptions in C-60 can be viewed as those localized on pyracylene units or on the corannulene/fluoranthene units. The bond orders and ring currents show the hexagons to be similar to benzene rings. The pentagon-hexagon bonds are also found to be longer than the hexagon-hexagon bonds.

Electron donor-acceptor complex of ICl with diethyl ether - He I photoelectron spectroscopy and ab initio molecular orbital study

The He I photoelectron spectrum of the diethyl ether-ICl complex has been obtained. The oxygen orbitals are shifted to higher binding energies and that of ICl to lower binding energies owing to complex formation. Ab initio molecular orbital (MO) calculations of the complex molecule showed that the bonding is between the sigma-type lone pair of oxygen and the I atom and that the complex has C-2v symmetry. The binding energy of the complex is computed to be 8.06 kcal mol(-1) at the MP2/3-21G* level. The orbital energies obtained from the photoelectron spectra of the complex are compared and assigned with orbital energies obtained by MO calculations. Natural bond orbital analysis (NBO) shows that charge transfer is from the sigma-type oxygen lone pair to the iodine atom and the magnitude of charge transfer is 0.0744 e.

Solid-state photobehaviour and crystal packing of o-chlorobenzylidene-DL-piperitone: influence of molecular topology on photobehaviour

C17H19ClO, M(r) = 274.7, triclinic, P1BAR, a = 11.154 (3), b = 12.685 (2), c = 12.713 (2) angstrom, alpha = 100.68 (1), beta = 113.58 (1), gamma = 104.50 (2)-degrees, V = 1511.1 (6) angstrom3, Z = 4, D(m) = 1.22, D(x) = 1.215 Mg m-3, Cu K-alpha, lambda = 1.5418 angstrom, mu = 2.16 mm-1, F(000) = 584, T = 293 K, R = 0.057 for 3481 observed reflections. The title compound is photostable in the crystalline state and lattice-energy calculations have been employed to rationalize the photobehaviour. The well-known beta-steering ability of the chloro group is not operative in this system as there are no Cl...Cl interactions in the crystal lattice. All five benzylidene-DL-piperitone structures so far studied are alpha-packed and the molecular topology appears to be a deciding factor even in the presence of steering groups.

Synthesis and characterization of photo-crosslinkable main-chain liquid-crystalline polymers containing bis(benzylidene)cycloalkanone units

A series of new photo-crosslinkable main-chain liquid-crystalline polymers containing bis(benzylidene)cycloalkanone units have been studied. These units in the polymers function as mesogens as well as photoactive centres. Polyesters with three different bis(4-hydroxybenzylidene)cycloalkanones corresponding to three cycloalkanones, namely cyclopentanone, cyclohexanone and cycloheptanone, have been prepared. Three dicarboxylic acids with ether linkages, which were derived from oligoethylene oxides, namely triethylene glycol, tetraethylene glycol and pentaethylene glycol, have been used as spacers in these polymers. Polymerization was carried out by both solution and interfacial polycondensation; the latter method gave high-molecular-weight polymers. Structural characterizations were done by ultra-violet, infra-red and H-1 nuclear magnetic resonance spectroscopy. Liquid-crystalline properties were studied by differential scanning calorimetry and polarized-light optical microscopy. These polymers show a nematic mesophase. Liquid-crystalline transition temperatures were correlated with polymer structure. The decrease in transition temperature with increase in cycloalkanone ring size was explained in terms of the change in geometrical anisotropy of bis(benzylidene)cycloalkanone units. MNDO (modified neglect of differential overlap) calculations were performed on the model compounds, bis(4-acetyloxybenzylidene)cycloalkanone to elucidate the geometrical variation of the mesogenic units with cycloalkanone ring size. Studies of photolysis reveal the two kinds of photoreactions that proceed in these polymer systems, namely photoisomerization and photo-crosslinking. The former reaction disrupts the parallel stacking of the chromophores and is reflected as an increase in the ultra-violet spectral intensity. The favourability of these two reactions depends on the mobility of the polymer chains. When the photolysis was done below T-g, photo-crosslinking dominates over photoisomerization. Above T-g, photoisomerization is followed by photo-crosslinking. The photosensitivity of the polymers decreases with increase in size of the cycloalkanone ring.

UVPES and ab initio molecular orbital studies on the electron donor-acceptor complexes of bromine with methylamines

The He I photoelectron spectra of bromine, methylamine, and their complex have been obtained, and the spectra show that lone-pair orbital energy of nitrogen in methylamine is stabilized by 1.8 eV and the bromine orbital energies are destabilized by about 0.5 eV due to complexation. Ab initio calculations have been performed on the charge-transfer complexes of Br-2 with ammonia and methyl-, dimethyl-, and trimethylamines at the 3-21G*, 6-311G, and 6-311G* levels and also with effective core potentials. Calculations predict donor and acceptor orbital energy shifts upon complexation, and there is a reasonable agreement between the calculated and experimental results. Complexation energies have been corrected for BSSE. Frequency analysis has confirmed that ammonia and trimethylamine form complexes with C-3v symmetry and methylamine and dimethylamine with C-s symmetry. Calculations reveal that the lone-pair orbital of nitrogen in amine and the sigma* orbital of Br-2 are involved in the charge-transfer interaction. LANL1DZ basis seems to be consistent and give a reliable estimate of the complexation energy. The computed complexation energies, orbital energy shifts, and natural bond orbital analysis show that the strength of the complex gradually increases from ammonia to trimethylamine.

Molecular structure and reactions of 1,1 '-bi(acenaphthen-1-ylidene)-2,2 '-dione

The molecular structure of 1,1'-bi(acenaphthen-1-ylidene)-2,2'-dione 1, a potential building-block for the synthesis of fullerene fragments, has been investigated by X-ray crystallography and semi-empirical (AM1 and PM3) calculations. There is a good agreement between the calculated and crystal structure which is essentially planar and has E-configuration. In the solid state, molecules of 1 pack in an interesting manner as corrugated sheets sustained by a network of C-H ... O hydrogen bonds and resulting in the formation of tetrameric loops. While steric factors limit the reactivity of the carbonyl groups in 1, the ene double bond of the ene-dione moiety present in it exhibits propensity toward [4 + 2]-cycloadditions to furnish novel and highly compressed polycycles 8-10.