DFT Modeling of Polyethylene Chains Cross-linked by Selected ns1 and ns2 Metal Atoms.

We analyze the structures, stabilities, and thermochemical properties of polyethylene (PE) oligomer chains cross-linked by metal (M) atoms through C-M-C bonds. Representative PEn-Mm-PEn complexes contain between 7 and 15 carbon atoms in each oligomer and one to three Li, Be, Mg, Zn, Ag, or Au cross-linking metal elements. PEn-Mm-PEn complexes are quasiplanar with nearly parallel PE chains. Their stability is determined by covalent C-M-C bonds accompanied by noncovalent dispersion interactions between PEn chains. Using the CAM-B3LYP+D3BJ+ABC functional, the binding energies of PE15-M-PE15 with respect to two PE15 radicals and metal fragments are -225, -230, -322, -551, -289, and -303 kJ/mol for Li, Ag, Au, Be, Mg, and Zn atoms, respectively. Entropy contributions (109 to 121 kJ/mol at 298.15 K) destabilize all complexes significantly. With two cross-linking metal elements in PE15-M2-PE15 complexes, binding energies are about double. Complexes with several open-shell Li, Ag, or Au doublet atoms have spins located on separated C-M-C bonds. High-spin PE15-Mm1-PE15-Mm2-PE15 complexes of three PE oligomers cross-linked by up to five doublet metal atoms create parallel PE tubes, which are suggested as elementary cells for modeling magnetic polymer tubes.

Accurate Nonlinear Optical Properties of Solvated para-Nitroaniline Predicted by an Electrostatic Discrete Local Field Approach.

A general computational protocol for accurate predictions of nonlinear optical (NLO) properties of solvated molecules based on the rigorous local field (RLF) approach taking all relevant effects into account is presented. para-Nitroaniline (pNA) was taken as a model NLO system dissolved in cyclohexane, tetrahydrofuran, and 1,4-dioxane. Molecular dynamics (MD) simulations employing either non-polarizable or polarizable force fields were used to generate representative sets of structures of the solutions. The static NLO properties of the solute were calculated at the MP2/aug-cc-pVDZ level of theory with the multiplicative scaling method used to account for the frequency dispersion of the properties. Focusing on the electric field-induced second harmonic generation (EFISH) and hyper-Rayleigh scattering (HRS), a good agreement between calculated results and experimental measurements was achieved with a polarizable force field. While the solvent effects on the vibrational contributions to the static molecular properties are significant, they remain small for both EFISH and HRS. Our results show that the proposed approach yields reliable predictions of dynamic NLO properties of solvated chromophores, which paves a route to further applications of the RLF approach to study a wide range of NLO phenomena in heterogeneous environments.

MS-CASPT2 study of the ground and low lying states of CsH.

Correlated ab initio methods [CASPT2 and R-CCSD(T)] in conjunction with the ANO-RCC basis sets in large contraction were used to calculate potential energy curves (PECs) of the ground and excited electronic states of CsH+ (doublets and quartets) with the inclusion of the scalar relativistic effects and spin-orbit interaction. The ground X2Σ+ state is a rather fragile van der Waals molecular ion. The binding energy of this X2Σ+ state provided by both computational methods is estimated to be 0.02-0.04 eV, and is compared with the reported experimental binding energy (0.51-0.77 eV). This large binding energy can be attributed to the A2Σ+ state, and can thus explain the apparent disagreement between theory and experiment. The spectroscopic constants of all bound states were calculated from the PECs and compared with previous published data for X2Σ+ and A2Σ+ states. Graphical abstract Low-lying Ω states of cesium hydride cation.

Caesium hydride: MS-CASPT2 potential energy curves and A1Σ+→X1Σ+ absorption/emission spectroscopy.

Correlated ab initio methods (CASPT2 and CCSD(T)) in conjunction with the ANO-RCC basis sets were used to calculate potential energy curves (PECs) of the ground, valence, and Rydberg electronic states of CsH with the inclusion of the scalar relativistic effects. The spectroscopic constants of bound states were calculated from the PECs and compared with previous theoretical and/or available experimental data. Absorption and emission spectra arising from the transition between X1Σ+ and A1Σ+ states were modelled using vibrational and rotational energy levels and corresponding nuclear wave functions obtained via the direct numerical integration of one-dimensional rovibrational Schrödinger equation in the CASPT2/ANO-RCC electronic potentials. The anharmonic shape of the A1Σ+ potential and the shape of the pertinent vibrational wave functions have an interesting impact on the final shape of the spectrum and result in the complicated fine structure of individual emission bands.

Electrodiffusion kinetics of ionic transport in a simple membrane channel.

We employ numerical techniques for solving time-dependent full Poisson-Nernst-Planck (PNP) equations in 2D to analyze transient behavior of a simple ion channel subject to a sudden electric potential jump across the membrane (voltage clamp). Calculated spatiotemporal profiles of the ionic concentrations and electric potential show that two principal exponential processes can be distinguished in the electrodiffusion kinetics, in agreement with original Planck's predictions. The initial fast process corresponds to the dielectric relaxation, while the steady state is approached in a second slower exponential process attributed to the nonlinear ionic redistribution. Effects of the model parameters such as the channel length, height of the potential step, boundary concentrations, permittivity of the channel interior, and ionic mobilities on electrodiffusion kinetics are studied. Numerical solutions are used to determine spatiotemporal profiles of the electric field, ionic fluxes, and both the conductive and displacement currents. We demonstrate that the displacement current is a significant transient component of the total electric current through the channel. The presented results provide additional information about the classical voltage-clamp problem and offer further physical insights into the mechanism of electrodiffusion. The used numerical approach can be readily extended to multi-ionic models with a more structured domain geometry in 2D or 3D, and it is directly applicable to other systems, such as synthetic nanopores, nanofluidic channels, and nanopipettes.

Atmospheric reactivity of CH2ICl with OH radicals: high-level OVOS CCSD(T) calculations for the X-abstraction pathways (X = H, Cl, or I).

Reactants, weak molecular complexes, transition states, and products for the H-, Cl-, and I-abstraction channels in the reaction of OH radicals with chloroiodomethane CH(2)ICl as well as the energy profiles at 0 K have been determined using high-level all-electron ab initio methods. The results showed that all-electron DK-CCSD(T)/ANO-RCC approach performed very well in predicting the reactivity of iodine. In terms of activation enthalpy at 0 K, the energy profile for the Cl-abstraction showed that this reaction pathway is not energetically favorable in contrast to the two other channels (H- and I-abstractions), which are competitive. The H-abstraction was strongly exothermic (-87 kJ mol(-1)), while the I-abstraction was modestly endothermic (11.8 kJ mol(-1)). On the basis of our calculations including the following corrections to the potential energies: basis set saturation, valence and core-valence electron correlation, relativistic effects, spin-adaptation, vibration contributions, and tunneling corrections, rate constants were predicted using canonical transition state theory over the temperature range 250-500 K for each abstraction pathway. The overall rate constant at 298 K was estimated to be 4.29 × 10(-14) and 5.44 × 10(-14) cm(3) molecule(-1) s(-1) for complex and direct abstraction mechanisms, respectively. In addition, the overall rate constant computed at 277 K was used in the estimation of the atmospheric lifetime for CH(2)ICl. On the basis of our theoretical calculations, the atmospheric lifetime for the OH removal process is predicted to be close to 1 year. In terms of atmospheric lifetime, the OH reaction is not competitive with the Cl reaction and photolysis processes.

Off-Center Gaussian Functions, an Alternative Atomic Orbital Basis Set for Accurate Noncovalent Interaction Calculations of Large Systems.

Proper description of noncovalent interactions requires, among other things, the use of diffuse atomic orbital (AO) basis sets. However, the presence of diffuse functions, especially in extended molecular systems, can lead to linear dependent AO basis sets. This in turn results, for example, in molecular orbital optimization problems or, when dependencies are removed in unpredictable and possibly geometry-dependent accuracy fluctuations. In this work, an alternative approach is proposed which suffers no linear dependence problems and delivers comparably accurate noncovalent interaction energies. An algorithm is proposed and implemented to construct a grid of off-center s-type Gaussian functions surrounding the molecule; substituting the presence of atom-centered diffuse basis functions. While the number of basis functions in the grid is comparable to the number of diffuse basis functions in aug-cc-pVXZ (for each cardinality number "X") basis sets for small molecular systems, the ratio becomes more favorable with increasing system size. The grid is constructed in a way that it is unique for a molecule (monomer) and, thus, independent of noncovalent complex/cluster geometry. The grid parameters, such as the density of grid points and s-function exponents, are obtained via optimization toward the S22 data set on the MP2 level. The quality, transferability, and versatility of the grid is tested on the S66 data set as well as on several cuts through the potential energy surface for noncovalent complexes, such as methyl-guanine···methyl-cytosine conversion from stacked to hydrogen-bonded structure.

Electron affinities of uracil: microsolvation effects and polarizable continuum model.

We present adiabatic electron affinities (AEAs) and the vertical detachment energies (VDEs) of the uracil molecule interacting with one to five water molecules. Credibility of MP2 and DFT/B3LYP calculations is supported by comparison with available benchmark CCSD(T) data. AEAs and VDEs obtained by MP2 and DFT/B3LYP methods copy trends of benchmark CCSD(T) results for the free uracil and uracil-water complexes in the gas phase being by 0.20 - 0.28 eV higher than CCSD(T) values depending on the particular structure of the complex. AEAs and VDEs from MP2 are underestimated by 0.09-0.15 eV. For the free uracil and uracil-(H(2)O)(n) (n = 1,2,3,5) complexes, we also consider the polarizable continuum model (PCM) and discuss the importance of the microsolvation when combined with PCM. AEAs and VDEs of uracil and uracil-water complexes enhance rapidly with increasing relative dielectric constant (ε) of the solvent. Highest AEAs and VDEs of the U(H(2)O)(5) complexes from B3LYP with ε = 78.4 are 2.03 and 2.81 eV, respectively, utilizing the correction from CCSD(T). Specific structural features of the microsolvated uracil-(H(2)O)(n) complexes and their anions are preserved also upon considering PCM in calculations of AEAs and VDEs.

Coupled Cluster and Møller-Plesset Perturbation Theory Calculations of Noncovalent Intermolecular Interactions using Density Fitting with Auxiliary Basis Sets from Cholesky Decompositions.

We compute noncovalent intermolecular interaction energies for the S22 test set [Phys. Chem. Chem. Phys.2006, 8, 1985-1993] of molecules at the Møller-Plesset and coupled cluster levels of supermolecular theory using density fitting (DF) to approximate all two-electron integrals. The error due to the DF approximation is analyzed for a range of auxiliary basis sets derived from Cholesky decomposition (CD) in conjunction with correlation consistent and atomic natural orbital valence basis sets. A Cholesky decomposition threshold of 10(-4)Eh for full molecular CD and its one-center approximation (1C-CD) generally yields errors below 0.03 kcal/mol, whereas 10(-3)Eh is sufficient to obtain the same level of accuracy or better with the atomic CD (aCD) and atomic compact CD (acCD) auxiliary basis sets. Comparing to commonly used predefined auxiliary basis sets, we find that while the aCD and acCD sets are larger by a factor of 2-4 with triple-ζ AO basis sets, they provide results 1-2 orders of magnitude more accurate.

Electron affinities of small uracil-water complexes: a comparison of benchmark CCSD(T) calculations with DFT.

We present benchmark CCSD(T) calculations of the adiabatic electron affinities (AEA) and the vertical detachment energies (VDE) of the uracil molecule interacting with one to three water molecules. Calculations with rather large aug-cc-pVTZ basis set were only tractable when the space of virtual orbitals was reduced to about 60% of the full space employing the OVOS (Optimized Virtual Orbital Space) technique. Because of the microhydration, the valence-bound uracil anion is stabilized leading to gradually more positive values of both AEA and VDE with increasing number of participating water molecules. This agrees with experimental findings. Upon hydration by three water molecules, the electron affinity of uracil increased in comparison with AEA of the isolated molecule by about 250 up to 570 meV, depending on the geometry of the complex. CCSD(T) results confirm trends determined by DFT calculations of the microhydrated uracil and its anion, even if electron affinities of the free and hydrated uracil molecule are overestimated by DFT by up to 300 meV.

MOLCAS 7: the next generation.

Some of the new unique features of the MOLCAS quantum chemistry package version 7 are presented in this report. In particular, the Cholesky decomposition method applied to some quantum chemical methods is described. This approach is used both in the context of a straight forward approximation of the two-electron integrals and in the generation of so-called auxiliary basis sets. The article describes how the method is implemented for most known wave functions models: self-consistent field, density functional theory, 2nd order perturbation theory, complete-active space self-consistent field multiconfigurational reference 2nd order perturbation theory, and coupled-cluster methods. The report further elaborates on the implementation of a restricted-active space self-consistent field reference function in conjunction with 2nd order perturbation theory. The average atomic natural orbital basis for relativistic calculations, covering the whole periodic table, are described and associated unique properties are demonstrated. Furthermore, the use of the arbitrary order Douglas-Kroll-Hess transformation for one-component relativistic calculations and its implementation are discussed. This section especially focuses on the implementation of the so-called picture-change-free atomic orbital property integrals. Moreover, the ElectroStatic Potential Fitted scheme, a version of a quantum mechanics/molecular mechanics hybrid method implemented in MOLCAS, is described and discussed. Finally, the report discusses the use of the MOLCAS package for advanced studies of photo chemical phenomena and the usefulness of the algorithms for constrained geometry optimization in MOLCAS in association with such studies.

Scaled MP3 non-covalent interaction energies agree closely with accurate CCSD(T) benchmark data.

Scaled MP3 interaction energies calculated as a sum of MP2/CBS (complete basis set limit) interaction energies and scaled third-order energy contributions obtained in small or medium size basis sets agree very closely with the estimated CCSD(T)/CBS interaction energies for the 22 H-bonded, dispersion-controlled and mixed non-covalent complexes from the S22 data set. Performance of this so-called MP2.5 (third-order scaling factor of 0.5) method has also been tested for 33 nucleic acid base pairs and two stacked conformers of porphine dimer. In all the test cases, performance of the MP2.5 method was shown to be superior to the scaled spin-component MP2 based methods, e.g. SCS-MP2, SCSN-MP2 and SCS(MI)-MP2. In particular, a very balanced treatment of hydrogen-bonded compared to stacked complexes is achieved with MP2.5. The main advantage of the approach is that it employs only a single empirical parameter and is thus biased by two rigorously defined, asymptotically correct ab-initio methods, MP2 and MP3. The method is proposed as an accurate but computationally feasible alternative to CCSD(T) for the computation of the properties of various kinds of non-covalently bound systems.

Toward more efficient CCSD(T) calculations of intermolecular interactions in model hydrogen-bonded and stacked dimers.

Interaction energies of the model H-bonded complexes, the formamide and formamidine dimers, as well as the stacked formaldehyde and ethylene dimers are calculated by the coupled cluster CCSD(T) method. These systems serve as a model for H-bonded and stacking interactions, typical in molecules participating in biological systems. We use the optimized virtual orbital space (OVOS) technique, by which the dimension of the space of virtual orbitals in coupled cluster CCSD(T) calculations can be significantly reduced. We demonstrate that when the space of virtual orbitals is reduced to 50% of the full space, which means reducing computational demands by 1 order of magnitude, the interaction energies for both H-bonded and stacked dimers are affected by no more than 0.1 kcal/mol. This error is much smaller than the error when interaction energies are calculated using limited basis sets.

Highly accurate CCSD(T) and DFT-SAPT stabilization energies of H-bonded and stacked structures of the uracil dimer.

The CCSD(T) interaction energies for the H-bonded and stacked structures of the uracil dimer are determined at the aug-cc-pVDZ and aug-cc-pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug-cc-pVDZ and aug-cc-pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug-cc-pVTZ and aug-cc-pVQZ basis sets) and extrapolated DeltaCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H-bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS-MP2 (SCS: spin-component-scaled), SCS(MI)-MP2 (MI: molecular interaction), MP3, dispersion-augmented DFT (DFT-D), M06-2X, and DFT-SAPT (SAPT: symmetry-adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)-MP2 method. Remarkably good binding energies are also obtained with the DFT-SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.