Ab initio molecular orbital calculations on ion-molecule and ion pair-molecule complexes of the water-lithium cyanide system

Ab initio molecular orbital (MO) calculations with the 3-21G and 6-31G basis sets were performed on a series of ion-molecule and ion pair-molecule complexes for the H2O + LiCN system. Stabilisation energies (with counter-poise corrections), geometrical parameters, internal force constants and harmonic vibrational frequencies were evaluated for 16 structures of interest. Although the interaction energies are smaller, the geometries and relative stabilities of the monohydrated contact ion pair are reminiscent of those computed for the complexes of the individual ions. Thus, interaction of the oxygen lone pair with lithium leads to a highly stabilised C2v structure, while the coordination of water to the cyanide ion involves a slightly non-linear hydrogen bond. Symmetrical bifurcated structures are computed to be saddle points on the potential energy surface, and to have an imaginary frequency for the rocking mode of the water molecule. On optimisation the geometries of the solvent shared ion pair structures (e.g. Li+cdots, three dots, centered OH2cdots, three dots, centered CN−) revealed a proton transfer from the water molecule leading to hydrogen bonded forms such as Li-O-Hcdots, three dots, centered HCN. The variation in the force constants and harmonic frequencies in the various structures considered are discussed in terms of ion-molecular and ion pair-molecule interactions.

Band to correlated crossover in alternating Hubbard and Pariser-Parr-Pople chains: Nature of the lowest singlet excitation of conjugated polymers

The evolution with increasing Coulomb correlations of a semiconductor to a magnetic insulator is related to an excited-state crossover in pi-electron models for conjugated polymers. We associate strong fluorescence with a lowest singlet excitation S1 that is dipole allowed, on the band side, while S1 becomes two-photon allowed on the correlated side. S1/S2 crossovers in Hubbard, Pariser-Parr-Pople, or other chains with electron-hole symmetry and alternating transfer integral t(1 +/- delta) are based on exact results at delta=0 and 1, on molecular exciton theory at large delta, and on oligomer calculations up to twelve sites.

Chiral polyhydroxylated tetrahydrothiophene derivatives: novel synthesis and structural elucidation by X-ray crystallography, NMR spectroscopy and molecular mechanics calculations

The dideoxygenation reaction of 1,3;4,6-di-O-alkylidene-2,5-di-S-methylthiocarbonyl-D-mannitol derivatives under Barton-McCombie reaction conditions gave the hexahydrodipyranothiophenes 4 and 7 instead of the expected 2,5-dideoxy products. Structural and conformational information on these novel derivatives has been obtained by NMR spectroscopy, single-crystal X-ray crystallography and molecular mechanics calculations.

Ab initio molecular orbital calculations on ion pair-water complexes of metal halides and oxides

Ab initio MO calculations are performed on a series of ion-molecular and ion pair-molecular complexes of H2O + MX (MX = LiF, LiCl, NaCl, BeO and MgO) systems. BSSE-corrected stabilization energies, optimized geometrical parameters, internal force constants and harmonic vibrational frequencies have been evaluated for all the structures of interest. The trends observed in the geometrical parameters and other properties calculated for the mono-hydrated contact ion pair complexes parallel those computed for the complexes of the individual ions. The bifurcated structures are found to be saddle points with an imaginary frequency corresponding to the rocking mode of water molecules. The solvent-shared ion pair complexes have high interaction energies. Trends in the internal force constant and harmonic frequency values are discussed in terms of ion-molecular and ion-pair molecular interactions.

Extended Huckel molecular orbital calculations on hydroxobridged heterotrinuclear cobalt(III) complexes

EHT calculations on heterotrinuclear cobalt(III) complexes of the type [Cu{(OH)(2)Co(L(4))}(2)](4+) where L(4) denotes (en)(2) or (NH3)(4), en = ethylenediamine and their component species have been carried out. The results regarding bonding and structure for the trinuclear complexes are compared with those for the monomer components such as [Co(en)(2)(OH)(2)](+), [Co(NH3)(4)(OH)(2)](+) and [Cu(OH)(4)](2-) are discussed.

H2O-CH4 and H2S-CH4 complexes: a direct comparison through molecular beam experiments and ab initio calculations

New molecular beam scattering experiments have been performed to measure the total ( elastic plus inelastic) cross sections as a function of the velocity in collisions between water and hydrogen sulfide projectile molecules and the methane target. Measured data have been exploited to characterize the range and strength of the intermolecular interaction in such systems, which are of relevance as they drive the gas phase molecular dynamics and the clathrate formation. Complementary information has been obtained by rotational spectra, recorded for the hydrogen sulfide-methane complex, with a pulsed nozzle Fourier transform microwave spectrometer. Extensive ab initio calculations have been performed to rationalize all the experimental findings. The combination of experimental and theoretical information has established the ground for the understanding of the nature of the interaction and allows for its basic components to be modelled, including charge transfer, in these weakly bound systems. The intermolecular potential for H2S-CH4 is significantly less anisotropic than for H2O-CH4, although both of them have potential minima that can be characterized as `hydrogen bonded'.

Density functional theory embedding for correlated wavefunctions

Methods that exploit the intrinsic locality of molecular interactions show significant promise in making tractable the electronic structure calculation of large-scale systems. In particular, embedded density functional theory (e-DFT) offers a formally exact approach to electronic structure calculations in which the interactions between subsystems are evaluated in terms of their electronic density. In the following dissertation, methodological advances of embedded density functional theory are described, numerically tested, and applied to real chemical systems.

First, we describe an e-DFT protocol in which the non-additive kinetic energy component of the embedding potential is treated exactly. Then, we present a general implementation of the exact calculation of the non-additive kinetic potential (NAKP) and apply it to molecular systems. We demonstrate that the implementation using the exact NAKP is in excellent agreement with reference Kohn-Sham calculations, whereas the approximate functionals lead to qualitative failures in the calculated energies and equilibrium structures.

Next, we introduce density-embedding techniques to enable the accurate and stable calculation of correlated wavefunction (CW) in complex environments. Embedding potentials calculated using e-DFT introduce the effect of the environment on a subsystem for CW calculations (WFT-in-DFT). We demonstrate that WFT-in-DFT calculations are in good agreement with CW calculations performed on the full complex.

We significantly improve the numerics of the algorithm by enforcing orthogonality between subsystems by introduction of a projection operator. Utilizing the projection-based embedding scheme, we rigorously analyze the sources of error in quantum embedding calculations in which an active subsystem is treated using CWs, and the remainder using density functional theory. We show that the embedding potential felt by the electrons in the active subsystem makes only a small contribution to the error of the method, whereas the error in the nonadditive exchange-correlation energy dominates. We develop an algorithm which corrects this term and demonstrate the accuracy of this corrected embedding scheme.

X-ray crystal structure and molecular mechanics calculations of (2,6-iso-dipropyl-phenylamide) dimethyl (tetra-methylcyclopenta-dienyl) silane titanium dichloride

Crystal and molecular structure of (2.6-dipropylphenylamide) dimethyl (tetra-methyl cyclopentadienyl) silane titanium dichloride (I) was fully characterized by X-ray diffraction. The crystal is obtained from a mixture of ether/hexane as orthorhombic. with a = 12.658 (3) Angstrom. b = 16.62 (3) Angstrom. c = 11.760 (2) Angstrom. V = 2474.2 (9) Angstrom(3). Z = 4, space group Pnma. R = 0.0399; Componud I compose of the pi-bounded ring with its dimethylsilyl-dipropyl phenyl amido group and the two terminal chloride atoms coordinated to central metal to form a so-called constrained geometry catalyst (CGC) structure. The result of molecular mechanics (MM) calculations on compound I shows that bond lengths and bond angles from the MM calculation are comparable to the data obtained from the X-ray diffraction study. The relation of the structure of CGCs and their catalytic activity by MM calculations is also discussed.

Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free-energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behavior above the transition temperature.

Correlated electron transport in molecular electronics

Theoretical and experimental values to date for the resistances of single molecules commonly disagree by orders of magnitude. By reformulating the transport problem using boundary conditions suitable for correlated many-electron systems, we approach electron transport across molecules from a new standpoint. Application of our correlated formalism to benzene-dithiol gives current-voltage characteristics close to experimental observations. The method can solve the open system quantum many-body problem accurately, treats spin exactly, and is valid beyond the linear response regime.

Beyond Ehrenfest: correlated non-adiabatic molecular dynamics

A method for introducing correlations between electrons and ions that is computationally affordable is described. The central assumption is that the ionic wavefunctions are narrow, which makes possible a moment expansion for the full density matrix. To make the problem tractable we reduce the remaining many-electron problem to a single-electron problem by performing a trace over all electronic degrees of freedom except one. This introduces both one- and two-electron quantities into the equations of motion. Quantities depending on more than one electron are removed by making a Hartree-Fock approximation. Using the first-moment approximation, we perform a number of tight binding simulations of the effect of an electric current on a mobile atom. The classical contribution to the ionic kinetic energy exhibits cooling and is independent of the bias. The quantum contribution exhibits strong heating, with the heating rate proportional to the bias. However, increased scattering of electrons with increasing ionic kinetic energy is not observed. This effect requires the introduction of the second moment.

Reduced-size polarized basis sets for calculations of molecular electric properties. III. Second-row atoms

Reduced-size polarized (ZmPolX) basis sets are developed for the second-row atoms X = Si, P, S, and Cl. The generation of these basis sets follows from a simple physical model of the polarization effect of the external electric field which leads to highly compact polarization functions to be added to the chosen initial basis set. The performance of the ZmPolX sets has been investigated in calculations of molecular dipole moments and polarizabilities. Only a small deterioration of the quality of the calculated molecular electric properties has been found. Simultaneously the size of the present reduced-size ZmPolX basis sets is about one-third smaller than that of the usual polarized (PolX) sets. This reduction considerably widens the range of applications of the ZmPolX sets in calculations of molecular dipole moments, dipole polarizabilities, and related properties.

On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies

It is well known that shape corrections have to be applied to the local-density (LDA) and generalized gradient (GGA) approximations to the Kohn-Sham exchange-correlation potential in order to obtain reliable response properties in time dependent density functional theory calculations. Here we demonstrate that it is an oversimplified view that these shape corrections concern primarily the asymptotic part of the potential, and that they affect only Rydberg type transitions. The performance is assessed of two shape-corrected Kohn-Sham potentials, the gradient-regulated asymptotic connection procedure applied to the Becke-Perdew potential (BP-GRAC) and the statistical averaging of (model) orbital potentials (SAOP), versus LDA and GGA potentials, in molecular response calculations of the static average polarizability alpha, the Cauchy coefficient S-4, and the static average hyperpolarizability beta. The nature of the distortions of the LDA/GGA potentials is highlighted and it is shown that they introduce many spurious excited states at too low energy which may mix with valence excited states, resulting in wrong excited state compositions. They also lead to wrong oscillator strengths and thus to a wrong spectral structure of properties like the polarizability. LDA, Becke-Lee-Yang-Parr (BLYP), and Becke-Perdew (BP) characteristically underestimate contributions to alpha and S-4 from bound Rydberg-type states and overestimate those from the continuum. Cancellation of the errors in these contributions occasionally produces fortuitously good results. The distortions of the LDA, BLYP, and BP spectra are related to the deficiencies of the LDA/GGA potentials in both the bulk and outer molecular regions. In contrast, both SAOP and BP-GRAC potentials produce high quality polarizabilities for 21 molecules and also reliable Cauchy moments and hyperpolarizabilities for the selected molecules. The analysis for the N-2 molecule shows, that both SAOP and BP-GRAC yield reliable energies omega(i) and oscillator strengths f(i) of individual excitations, so that they reproduce well the spectral structure of alpha and S-4.(C) 2002 American Institute of Physics.

PAMOP: Petascale atomic, molecular and optical collision calculations

Petaflop architectures are currently being utilized efficiently to perform large scale computations in Atomic, Molecular and Optical Collisions. We solve the Schr\"odinger or Dirac equation for the appropriate collision problem using the R-matrix or R-matrix with pseudo-states approach. We briefly outline the parallel methodology used and implemented for the current suite of Breit-Pauli and DARC codes. In this report, various examples are shown from our theoretical results compared with experimental results obtained from Synchrotron Radiation facilities where the Cray architecture at HLRS is playing an integral part in our computational projects.