Relative energy of the high-(5T2g) and low-(1A1g) spin states of the ferrous complexes [Fe(L)(NHS4)]: CASPT2 versus density functional theory

High-level ab initio calculations using multiconfigurational perturbation theory [complete active space with second-order perturbation theory (CASPT2)] were performed on the transition energy between the lowest high-spin (corresponding to (5T2g) in Oh) and low-spin (corresponding to 1A1g in Oh) states in the series of six-coordinated Fe(II) molecules [Fe(L)(NHS4)], where NHS4 is 2,2'-bis(2-mercaptophenylthio)diethylamine dianion and L=NH3, N2H4, PMe3, CO, and NO+. The results are compared to (previous and presently obtained) results from density functional theory (DFT) calculations with four functionals, which were already shown previously by Casida and co-workers [Fouqueau et al., J. Chem. Phys. 120, 9473 (2004); Ganzenmuller et al., ibid. 122, 234321 (2005); Fouqueau et al., ibid. 122, 044110 (2005); Lawson Daku et al., ChemPhysChem 6, 1393 (2005)] to perform well for the spin-pairing problem in these and other Fe(II) complexes, i.e., OLYP, PBE0, B3LYP, and B3LYP*. Very extended basis sets were used both for the DFT and CASPT2 calculations and were shown to be necessary to obtain quantitative results with both types of method. This work presents a sequel to a previous DFT/CASPT2 study of the same property in the complexes [Fe(H2O)6]2+, [Fe(NH3)6]2+, and [Fe(bpy)3]2+ [Pierloot et al., J. Chem. Phys. 125, 124303 (2006)]. The latter work was extended with new results obtained with larger basis sets and including the OLYP functional. For all considered complexes, the CASPT2 method predicts the correct ground state spin multiplicity. Since experimental data for the actual quintet-singlet (free) energy differences are not available, the performance of the different DFT functionals was judged based on the comparison between the DFT and CASPT2 results. From this, it was concluded that the generalized gradient OLYP functional performs remarkably well for the present series of ferrous compounds, whereas the success of the three hybrid functionals varies from case to case.     

Sliding friction between polymer surfaces: a molecular interpretation

For two contacting rigid bodies, the friction force F is proportional to the normal load and independent of the macroscopic contact area and relative velocity V (Amonton's law). With two mutually sliding polymer samples, the surface irregularities transmit deformation to the underlying material. Energy loss along the deformation cycles is responsible for the friction force, which now appears to depend strongly on V [see, e.g., N. Maeda et al., Science 297, 379 (2002)]. We base our theoretical interpretation on the assumption that polymer chains are mainly subjected to oscillatory "reptation" along their "tubes." At high deformation frequencies-i.e., with a large sliding velocity V-the internal viscosity due to the rotational energy barriers around chain bonds hinders intramolecular mobility. As a result, energy dissipation and the correlated friction force strongly diminish at large V. Derived from a linear differential equation for chain dynamics, our results are basically consistent with the experimental data by Maeda et al. [Science 297, 379 (2002)] on modified polystyrene. Although the bulk polymer is below T(g), we regard the first few chain layers below the surface to be in the liquid state. In particular, the observed maximum of F vs V is consistent with physically reasonable values of the molecular parameters. As a general result, the ratio FV is a steadily decreasing function of V, tending to V(-2) for large velocities. We evaluate a much smaller friction for a cross-linked polymer under the assumption that the junctions are effectively immobile, also in agreement with the experimental results of Maeda et al. [Science 297, 379 (2002)].     

Pseudo-Jahn-Teller origin of geometry and pseudorotations in second row tetra-atomic clusters X4 (X=Na, Mg, Al, Si, P, S)

Experimentally determined or ab initio calculated molecular geometries carry no information about their origin. Employing the Jahn-Teller (JT) vibronic coupling effects as the only source of instability and consequent distortions of high-symmetry molecular configurations, we have worked out a procedure that allows us to trace the origin of particular geometries and determine the detailed electronic mechanism of their formation. This procedure is illustrated by considering a series of X(4) clusters with X=Na, Mg, Al, Si, P, and S. It shows explicitly why Na(4), Si(4), and Al(4) have a rhombic geometry in the ground state, while Mg(4) and P(4) are tetrahedral, whereas S(4) is a trapezium. Even when the minimum-energy geometries are the same (as in the case of rhombic Na(4), Si(4), and Al(4)), the electronic mechanism of their formation is quite different. In particular, in Na(4) and Si(4) the rhombic minima are produced by a strong pseudo JT coupling between two excited states in the square-planar configuration (different in the two cases) that stabilizes one of them and makes it the ground state by rhombic distortions. The rhombic configuration of Al(4) is due to the pseudo JT effect in its ground-state square-planar configuration, and the trapezium in S(4) is formed by two pseudo JT couplings essentially involving excited states. In several cases this analysis shows also the tunneling paths between equivalent configurations.     

Space-time contours to treat intense field-dressed molecular states. I. Theory

A molecular system exposed to an intense external field is considered. The strength of the field is measured by the number L of electronic states that become populated during this process. In the present article the authors discuss a rigorous way, based on the recently introduced space-time contours [R. Baer, et al., J. Chem. Phys. 119, 6998 (2003)], to form N coupled Schrodinger equations where N相似文献   

Equation of state of nitrogen (N2) at high pressures and high temperatures: molecular dynamics simulation

Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures above 800 K was studied by molecular dynamics (MD) simulations. The dynamics of the N(2) molecules is treated in hard rotor approximation, i.e., it accounts both translational and rotational degrees of freedom. The rotational motion of the N(2) molecule is treated assuming constant moment of inertia of the nitrogen molecule. The new MD program fully accounts anisotropic molecular nitrogen interaction. The N(2)-N(2) interaction potential has been derived by van der Avoird et al. [J. Chem. Phys. 84, 1629 (1986)] using the results of high precision Hartree-Fock ab initio quantum mechanical calculations. The potential, fully accounts rotational symmetry of the N(2)-N(2) system, by employing 6-j Wigner symbols, i.e., preserving full rotational symmetry of the system. Various numerical algorithms were tested, in order to achieve the energy preservation during the simulation. It has been demonstrated that the standard Verlet algorithm was not preserving the energy for the standard MD time step, equal to 5x10(-16) s. Runge-Kutta fourth order method was able to preserve the energy within 10(-4) relative error, but it requires calculation of the force four times for each time step and therefore it is highly inefficient. A predictor-corrector method of the fifth order (PC5) was found to be efficient and precise and was therefore adopted for the simulation of the molecular nitrogen properties at high pressure. Singer and Fincham algorithms were tested and were found to be as precise as PC5 algorithm and they were also used in the simulation of the equation of state. Results of MD simulations are in very good agreement with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars). For higher pressures, up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 85, 1262 (1984)], the obtained numerical results provide new data of the experimentally unexplored region. These data were formulated in the analytical form of pressure-density-temperature equation of state.     

On the equilibrium structures of the complexes H2C3H+ · Ar and c-C3H3(+) · Ar: results of explicitly correlated coupled cluster calculations

Explicitly correlated coupled cluster theory at the CCSD(T)-F12x (x = a, b) level [T. B. Adler et al., J. Chem. Phys. 127, 221106 (2007)] has been employed in a study of the potential energy surfaces for the complexes H(2)C(3)H(+) · Ar and c-C(3)H(3)(+) · Ar. For the former complex, a pronounced minimum with C(s) symmetry was found (D(e) ≈ 780 cm(-1)), well below the local "H-bound" minimum with C(2v) symmetry (D(e) ≈ 585 cm(-1)). The absorption at 3238 cm(-1) found in the recent infrared photodissociation spectra [A. M. Ricks et al., J. Chem. Phys. 132, 051101 (2010)] is, thus, interpreted as an essentially free acetylenic CH stretching vibration of the propargyl cation. A global minimum of C(s) symmetry was also obtained for c-C(3)H(3)(+) (D(e) ≈ 580 cm(-1)), but the energy difference with respect to the local C(2v) minimum is only 54 cm(-1).     

The adiabatic‐to‐diabatic transformation angle and the berry phase for coupled jahn–teller/renner–teller systems: The F + H2 as a case study

The approach to calculate improved, two‐state, adiabatic‐to‐diabatic transformation angles (also known as mixing angles), presented before (see Das et al., J Chem Phys 2010, 133, 084107), was used here while studying the F + H₂ system. However, this study is characterized by two new features: (a) it is the first of its kind in which is studied the interplay between Renner–Teller (RT) and Jahn–Teller (JT) nonadiabatic coupling terms (NACT); (b) it is the first of its kind in which is reported the effect of an upper singular RT‐NACT on a lower two‐state (JT) mixing angle. The fact that the upper NACT is singular (it is shown to be a quasi‐Dirac δ‐function) enables a semi‐analytical solution for this perturbed mixing angle. The present treatment, performed for the F + H₂ system, revealed the existence of a novel parameter, η, the Jahn–Renner coupling parameter (JRCP), which yields, in an unambiguous way, the right intensity of the RT coupling (as resembled, in this case, by the quasi‐Dirac δ‐function) responsible for the fact that the final end‐of‐the contour angle (identified with the Berry phase) is properly quantized. This study implies that the numerical value of this parameter is a pure number (independent of the molecular system): η = $ 2\sqrt 2 /\pi $ (= 0.9003) and that there is a good possibility that this value is a novel characteristic molecular constant for a certain class of tri‐atomic systems. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Collision-induced absorption by H2 pairs: from hundreds to thousands of kelvin

An interaction-induced dipole surface (IDS) and a potential energy surface (PES) of collisionally interacting molecular hydrogen pairs H(2)-H(2) was recently obtained using quantum chemical methods (Li, X.; et al. Computational Methods in Science and Engineering, ICCMSE. AIP Conf. Proc. 2009, ; see also Li, X.; et al. Int. J. Spectrosc. 2010, ID 371201). The data account for substantial rotovibrational excitations of the H(2) molecules, as encountered at temperatures of thousands of kelvin (e.g., in the atmospheres of "cool" stars). In this work we use these results to compute the binary collision-induced absorption (CIA) spectra of dense hydrogen gas in the infrared at temperatures up to several thousand kelvin. The principal interest of the work is in the spectra at such higher temperatures, but we also compare our computations with existing laboratory measurements of CIA spectra of dense hydrogen gas and find agreement.     

Equilibrium geometries of low-lying isomers of some Li clusters, within Hartree-Fock theory plus bond order or MP2 correlation corrections

In a recent study by Kornath et al. [J. Chem. Phys. 118, 6957 (2003)], the Li(n) clusters with n=2, 4, and 8 have been isolated in argon matrices at 15 K and characterized by Raman spectroscopy. This has prompted us to carry out a theoretical study on such clusters up to n=10, using Hartree-Fock theory, plus low-order M?ller-Plesset perturbation corrections. To check against the above study of Kornath et al., as a by-product we have made the same approximations for n=6 and 8 as we have for n=10. This has led us to emphasize trends with n through the Li(n) clusters for (i) ground-state energy, (ii) HOMO-LUMO energy gap, (iii) dissociation energy, and (iv) Hartree-Fock eigenvalue sum. The role of electron correlation in distinguishing between low-lying isomers is plainly crucial, and will need a combination of experiment and theory to obtain decisive results such as that of Kornath et al. for Li(8). In particular, it is shown that Hartree-Fock theory plus bond order correlations does account for the experimentally observed symmetry T(d) symmetry for Li(8).     

Theoretical analysis of the Jahn-Teller distortions in tetrathiolato iron(II) complexes

Crystallographic studies of [Fe(SR)(4)](2-) (R is an alkyl or aryl residue) have shown that the Fe(II)S(4) cores of these complexes have (pseudo) D2d symmetry. Here we analyze the possibility that these structures result from a Jahn-Teller (JT) distortion that arises from the e(3z(2) - r(2), x(2) - y(2)) orbital ground state of Fe(II) in T(d)symmetry. Special attention is paid to the influence of the second-nearest neighbors of Fe, which lowers the symmetry and reduces the full JT effect to a smaller, pseudo JT effect (PJT). To estimate the size of the PJT distortion, we have determined the vibronic parameters and orbital state energies for a number of [Fe(SR)(4)](2-) models using density functional theory (DFT). Subsequently, this information is used for evaluating the adiabatic potential surfaces in the space of the JT-active coordinates of the FeS(4) moiety. The surfaces reveal that the JT effect of Fe(II) is completely quenched by the tetrathiolate coordination.     

Protein-matrix coupling/uncoupling in "dry" systems of photosynthetic reaction center embedded in trehalose/sucrose: the origin of trehalose peculiarity

Trehalose is a nonreducing disaccharide of glucose found in organisms, which can survive adverse conditions such as extreme drought and high temperatures. Furthermore, isolated structures, as enzymes or liposomes, embedded in trehalose are preserved against stressing conditions [see, e.g., Crowe, L. M. Comp. Biochem. Physiol. A 2002, 131, 505-513]. Among other hypotheses, such protective effect has been suggested to stem, in the case of proteins, from the formation of a water-mediated, hydrogen bond network, which anchors the protein surface to the water-sugar matrix, thus coupling the internal degrees of freedom of the biomolecule to those of the surroundings [Giuffrida, S.; et al. J. Phys. Chem. B 2003, 107, 13211-13217]. Analogous protective effect is also accomplished by other saccharides, although with a lower efficiency. Here, we studied the recombination kinetics of the primary, light-induced charge separated state (P(+)Q(A)(-)) and the thermal stability of the photosynthetic reaction center (RC) of Rhodobacter sphaeroides in trehalose-water and in sucrose-water matrixes of decreasing water content. Our data show that, in sucrose, at variance with trehalose, the system undergoes a "nanophase separation" when the water/sugar mole fraction is lower than the threshold level approximately 0.8. We rationalize this result assuming that the hydrogen bond network, which anchors the RC surface to its surrounding, is formed in trehalose but not in sucrose. We suggest that both the couplings, in the case of trehalose, and the nanophase separation, in the case of sucrose, start at low water content when the components of the system enter in competition for the residual water.     

Dithiane- and trithiane-based photolabile scaffolds for molecular recognition

[see structure]. A modular synthetic approach to novel dithiane- and trithiane-based photolabile molecular hosts equipped with elements of molecular recognition is developed. The approach provides ready access to a family of amino-derivatized photocleavable molecular systems capable of hydrogen-bonding-based recognition of biologically relevant molecules, e.g., ureas, barbiturates etc. These systems undergo efficient photofragmentation in the presence of external (e.g., benzophenone) or internal (e.g., nitropyridine) electron-transfer sensitizers.     

Diffusion in linear porous media with periodic entropy barriers: A tube formed by contacting spheres

The problem of transport in quasi-one-dimensional periodic structures has been studied recently by several groups [D. Reguera et al., Phys. Rev. Lett.96, 130603 (2006); P. S. Burada et al., Phys. Rev. E75, 051111 (2007); B. Q. Ai and L. G. Liu, ibid.74, 051114 (2006); B. Q. Ai et al., ibid.75, 061126 (2007); B. Q. Ai and L. G. Liu, J. Chem. Phys.126, 204706 (2007); 128, 024706 (2008); E. Yariv and K. D. Dorfman, Phys. Fluids19, 037101 (2007); N. Laachi et al., Europhys. Lett.80, 50009 (2007); A. M. Berezhkovskii et al., J. Chem. Phys.118, 7146 (2003); 119, 6991 (2003)]. Using the concept of "entropy barrier" [R. Zwanzig, J. Phys. Chem.96, 3926 (1992)] one can classify such structures based on the height of the entropy barrier. Structures with high barriers are formed by chambers, which are weakly connected with each other because they are connected by small apertures. To escape from such a chamber a diffusing particle has to climb a high entropy barrier to find an exit that takes a lot of time [I. V. Grigoriev et al., J. Chem. Phys.116, 9574 (2002)]. As a consequence, the particle intrachamber lifetime tau(esc) is much larger than its intrachamber equilibration time, tau(rel), tau(esc)>tau(rel). When the aperture is not small enough, the intrachamber escape and relaxation times are of the same order and the hierarchy fails. This is the case of low entropy barriers. Transport in this case is analyzed in the works of Schmid and co-workers, Liu and co-workers, and Dorfman and co-workers, while the work of Berezhkovskii et al. is devoted to diffusion in the case of high entropy barriers.     

Combined experimental-theoretical characterization of the hydrido-cobaloxime [HCo(dmgH)2(PnBu3)

A combined theoretical and experimental approach has been employed to characterize the hydrido-cobaloxime [HCo(dmgH)(2)(PnBu(3))] compound. This complex was originally investigated by Schrauzer et al. [Schrauzer et al., J. Am. Chem. Soc. 1971, 93,1505] and has since been referred to as a key, stable analogue of the hydride intermediate involved in hydrogen evolution catalyzed by cobaloxime compounds [Artero, V. et al. Angew. Chem., Int. Ed. 2011, 50, 7238-7266]. We employed quantum chemical calculations, using density functional theory and correlated RI-SCS-MP2 methods, to characterize the structural and electronic properties of the compound and observed important differences between the calculated (1)H NMR spectrum and that reported in the original study by Schrauzer and Holland. To calibrate the theoretical model, the stable hydrido tetraamine cobalt(III) complex [HCo(tmen)(2)(OH(2))](2+) (tmen = 2,3-dimethyl-butane-2,3-diamine) [Rahman, A. F. M. M. et al. Chem. Commun. 2003, 2748-2749] was subjected to a similar analysis, and, in this case, the calculated results agreed well with those obtained experimentally. As a follow-up to the computational work, the title hydrido-cobaloxime compound was synthesized and recharacterized experimentally, together with the Co(I) derivative, giving results that were in agreement with the theoretical predictions.     

The role of CH-π interaction in the charge transfer properties in tris(8-hydroxyquinolinato)aluminium(III)

The charge mobility is a key property in many electro-optical materials, with charge transfer (CT) taking place in a solid matrix of molecules. Large intermolecular electronic interaction is one of the key factors for a good CT rate, which is dependent on both intra- and intermolecular structures. The connection of the molecular structure with the intermolecular CT property would facilitate the search for a new material with desirable CT property, but currently it is still quite limited by the lack of knowledge for intermolecular configurations. In the present work, we study factors influencing the intermolecular configurations, and subsequently the CT property, in tris(8-hydroxyquinolinato) aluminium(III) (AlQ(3)) from all currently available crystal structures. We found that there exists a pair of CH-π interactions in a good majority of the π-π stacked bimolecular configurations. Such CH-π and π-π interacting structures are also seen in the crystal structures of many other similar molecules. With both experimental and simulated structures, we show that the CH-π interaction stabilizes the bimolecular configurations, and drives the structure towards a region with a higher electron transfer coupling and lower hole transfer coupling. This effect likely affects the electron transport property of AlQ(3), since it is consistent with recent experimental results, where AlQ(3) analogs with their CH-π interaction blocked either require a higher operating voltage in light-emitting devices [Sapochak et al., J. Am. Chem. Soc., 2001, 123, 6300], or become bipolar in their charge mobilities [Liao et al., J. Am. Chem. Soc., 2009, 131, 763]. CH-π interaction is commonly seen in aromatic molecules, which are frequently used as building blocks in molecules for electro-optical applications. Our work points out a possible way to enhance the desired CT property in the design of new materials.     

H2 storage in isostructural UiO-67 and UiO-66 MOFs

The recently discovered UiO-66/67/68 class of isostructural metallorganic frameworks (MOFs) [J. H. Cavka et al. J. Am. Chem. Soc., 2008, 130, 13850] has attracted great interest because of its remarkable stability at high temperatures, high pressures and in the presence of different solvents, acids and bases [L. Valenzano et al. Chem. Mater., 2011, 23, 1700]. UiO-66 is obtained by connecting Zr(6)O(4)(OH)(4) inorganic cornerstones with 1,4-benzene-dicarboxylate (BDC) as linker resulting in a cubic MOF, which has already been successfully reproduced in several laboratories. Here we report the first complete structural, vibrational and electronic characterization of the isostructural UiO-67 material, obtained using the longer 4,4'-biphenyl-dicarboxylate (BPDC) linker, by combining laboratory XRPD, Zr K-edge EXAFS, TGA, FTIR, and UV-Vis studies. Comparison between experimental and periodic calculations performed at the B3LYP level of theory allows a full understanding of the structural, vibrational and electronic properties of the material. Both materials have been tested for molecular hydrogen storage at high pressures and at liquid nitrogen temperature. In this regard, the use of a longer ligand has a double benefit: (i) it reduces the density of the material and (ii) it increases the Langmuir surface area from 1281 to 2483 m(2) g(-1) and the micropore volume from 0.43 to 0.85 cm(3) g(-1). As a consequence, the H(2) uptake at 38 bar and 77 K increases from 2.4 mass% for UiO-66 up to 4.6 mass% for the new UiO-67 material. This value is among the highest values reported so far but is lower than those reported for MIL-101, IRMOF-20 and MOF-177 under similar pressure and temperature conditions (6.1, 6.2 and 7.0 mass%, respectively) [A. G. Wong-Foy et al. J. Am. Chem. Soc., 2006, 128, 3494; M. Dinca and J. R. Long. Angew. Chem., Int. Ed., 2008, 47, 6766]. Nevertheless the remarkable chemical and thermal stability of UiO-67 and the absence of Cr in its structure would make this material competitive.     

Bound states of the positron with nitrile species with a configuration interaction multi-component molecular orbital approach

Characteristic features of the positron binding structure of some nitrile (-CN functional group) species such as acetonitrile, cyanoacetylene, acrylonitrile, and propionitrile are discussed with the configuration interaction scheme of multi-component molecular orbital calculations. This method can take the electron-positron correlation contribution into account through single electronic-single positronic excitation configurations. Our PA value of acetonitrile with the electronic 6-31++G(2df,2pd) and positronic [15s15p3d2f1g] basis set is calculated as 4.96 mhartree, which agrees to within 25% with the recent experimental value of 6.6 mhartree by Danielson et al. [Phys. Rev. Lett., 2010, 104, 233201]. Our PA values of acrylonitrile and propionitrile (5.70 and 6.04 mhartree) are the largest among these species, which is consistent with the relatively large dipole moments of the latter two systems.