Compact multipolar representation of the electrostatic potential for flexible molecules

A new method for generating a compact multipolar representation of the electrostatic potential (EP) for flexible molecules is presented. The method is based on a constrained minimization of the difference between the quantum mechanical and the classical EP. The fitting procedure used adopts the least absolute shrinkage and selection operator technique [R. Tibshirani, J. Roy. Stat. Soc. B 58, 267 (1996)] which can be seen as penalized ordinary least squares. The penalty function optimized for the particular molecule of interest effectively removes redundant multipoles. It is shown that the use of multiple conformations is crucial for the predictive ability of the EP model for flexible molecules. The multipole local coordinate systems are chosen in a way that best reflects the key conformational changes. It was demonstrated that such an approach improves the predictive ability of EP models. It also allows to exploit equivalence of atoms in the calculation of multipoles components. In the case of polar flexible molecules, the augmentation of the EP model based on charges by higher multipoles decreases the relative root mean square error by a factor of 1.5-5. The corresponding effect of enlargement of the set of multipoles was significantly reduced.     

Unimolecular photochemistry of 4-thiouracils

Unimolecular phototautomeric reactions in 4-thiouracil, 1-methyl-4-thiouracil and 6-aza-4-thiouracil were studied using the matrix-isolation technique combined with infrared absorption spectroscopy. For monomers of these compounds, isolated in solid argon at 10 K, an intramolecular proton-transfer photoreaction was observed. Upon UV (lambda > 345 nm) irradiation, the initial oxo-thione forms of 4-thiouracils were converted into the corresponding oxo-thiol tautomers. The photogenerated oxo-thiol isomers were identified by comparing their experimental IR spectra with the spectra theoretically calculated at the DFT(B3LYP)/6-311++G(2d,p) level. Good agreement between the observed and predicted pattern of spectral bands allowed a reliable identification. This is the first report on experimental observation of isomeric forms of 4-thiouracils other than the canonical oxo-thione tautomers.     

Direct calculation of lattice sums. A method to account for the crystal field effects

A method of direct calculation of lattice sums in three-dimensional crystals is reported. The method is based on annihilation of some lowest multipole moments of the unit cell by a redefinition of the unit cell content. As a result, properties of the infinite crystal can be calculated as usual by taking a finite cluster of unit cells, but surrounded by an additional surface layer of a charge density (e.g., a layer of point charges). This charge density distribution produces the electric field approximating that one of the rest of the infinite crystal. The method proposed is easily applicable in the SCFLCAO procedure as well as in any method using a cluster representation for an infinite crystal. The validity of the infinite crystal model for a finite crystal is also discussed.     

A study of polyamine complex formation with H+, Cu(II), Zn(II), Pb(II), and Mg(II) in aqueous solution

Summary The composition and stability of the following biogenic amine complexes have been investigated: 1,4-diaminobutane(Put), 4-azaoctane-1,8-diamine(Spd), 4,9-diazadodecan-1, 12-diamine(Spm) as well as homologues such as 1,3-diaminopropane(Put3), 4-azaheptane-1, 7-diamine(Spd3,3) and 4,8-diazaundecan-1,11-diamine(Spm3,3,3) with H⁺, Cu(II), Zn(II), Pb(II) and Mg(II). A potentiometric method was used. The VIS technique enabled the determination of coordination mode in copper/amine systems. It was found that Mg(II) does not form coordination compounds with any of the studied polyamines in solution. An increase in the concentration of ligand and metal was found to result in a stronger tendency towards the formation of protonated compounds accompanied by a decrease in the concentration of hydroxocomplexes. At physiologicalpH (7.4) an increase in the concentration of protonated compounds by approximately 15% was observed within the ligand concentration range from 0.001 mol dm^–3 to 0.0001 mol dm^–3 at a Cu(II) concentration of 0.000177 mol dm^–3.

---

Untersuchungen zur Komplexbildung von Polyaminen mit H⁺, Cu(II), Zn(II), Pb(II) und Mg(II) in wäßriger Lösung

Zusammenfassung Anhand einer Analyse von potentiometrischen Daten wurden Zusammensetzung und Beständigkeit folgender biogener Aminkomplexe untersucht: 1,4-Diaminobutan(Put), 4-Azaoktan-1,8-diamin(Spd), 4,9-Diazadodekan-1,12-diamin(Spm), sowie auch deren Homologen 1,3-Diaminopropan(Put3), 4-Azaheptan-1,7-diamin(Spd3,3) und 4,8-Diazaundekan-1,11-diamin(Spm3,3,3) mit H⁺, Cu(II), Zn(II), Pb(II) und Mg(II). Mit Hilfe der VIS-Technik wurde die Koordinationsweise in Kupfer/Amin-Systemen bestimmt. Es wurde festgestellt, daß Mg(II) keine Koordinationsverbindungen mit den untersuchten Polyaminen bildet. Eine höhere Konzentration von Ligand und Metall führte zu stärkerer Tendenz der Bildung protonierter Verbindungen, wobei die Konzentration von Hydroxokomplexen kleiner wurde. Bei physiologischempH (7.4) wurde im Bereich der Ligand-Konzentration von 0.001 mol dm^–3 bis 0.0001 mol dm^–3 bei einer Cu(II)-Konzentration von 0.000177 mol dm^–3 ein Anstieg der Konzentration protonierter Verbindungen um etwa 15% beobachtet.

     

Molecular tight-binding method. I. Many-electron method for hole bands of molecular crystals

A method for calculating ab initio electronic excitation energies of molecular crystals, based on a many-electron tight-binding approximation, is described. The method follows Frenkel's model for excitons and allows a many-electron treatment of the band-structure problem of molecular crystals. The case of hole bands is studied in detail and various versions of the method are considered. A computational scheme is proposed, in which approximate correlation corrections to the HFR matrix equations of the one-electron LCMO method are calculated. The main effects contributing to these corrections are the effect of relaxation of a molecular ion, the effect of intramolecular electronic-correlation change, and the effect of polarization of the remaining molecules in a crystal. The method developed in the present paper is applied to calculation of the hole bands of the HCP helium crystal.     

Multidimensional gravitational waves. I. Purely radiative spacetimes

The concept of gravitational radiation in multidimensional Einstein gravity is studied. Due to technical difficulties the Petrov classification of the conformal tensor and the peeling off theorem cannot be generalized to spaces of dimensiond > 4. A multidimensionalpp wave is defined by analogy with that ford = 4. If the definition is supplemented by the Lichnerowicz condition for purely radiative spacetimes, the admitted solutions of the dimensionally reduced gravity theory describe a scalar (dilaton) and a gravitationalpp wave moving in the external spacetime in the same direction. The existence of the extra dimensions can be traced back by measuring the frequence of thepp wave in the external space. The condition is equally restrictive in a more general case of the internal space being a flat anisotropic manifold (i.e. an abelian group space).     

The vibrational spectra and normal coordinate analysis of tetracyclopropyllead

All the fundamental frequencies observed in IR and Raman spectra have been assigned to the normal modes of the molecule (C₃H₅)₄Pb (C₃H₅ = cyclopropyl). The calculated force field is compared to that of (C₃H₅)₂Hg and the Cmetal stretching force constants are discussed along with those of Me₄Pb and Me₂Hg; The similarity of the spectra of cyclopropyllead and cyclopropyl-mercury proves that the vibrations of cyclopropyl structural units are isolated in both molecules.     

Photochemical syn-anti isomerization reactions in N2-hydroxyisocytosines--an experimental matrix isolation and theoretical study

N2-hydroxyisocytosine and 1-methyl-N2-hydroxyisocytosine were studied using a matrix isolation technique combined with infrared absorption spectroscopy. For N2-hydroxyisocytosine isolated in an Ar matrix (at 10 K), two imino-oxo isomers, one with the hydroxyimino =N-OH group directed toward the N1-H group (the form called further anti) and the second with the =N-OH group directed toward N3-H (syn), were observed in the ratio 1.4:1. The syn isomer is converted totally to the anti form after UV (lambda > 295 nm) irradiation of the matrix. A small amount of the N(3)H-hydroxy-amino tautomer of N2-hydroxyisocytosine was also detected in the matrix. This form did not react photochemically. For 1-methyl-N2-hydroxyisocytosine, only the syn form of the imino-oxo tautomer was observed after deposition of the matrix. UV (lambda > 295 nm) irradiation induced a photoreaction converting this isomer into the anti form. After 15% of the starting material had been converted into the product, a photostationary state was achieved, and no further progress of the reaction was observed. Subsequent UV irradiation (lambda > 335 nm) caused a back reaction, leading to a disappearance of the anti form and to the recovery of the initial syn isomer. All isomers were identified by comparing their experimental IR spectra with the spectra theoretically calculated at the DFT(B3LYP)/6-31G(d,p) level, where DFT is the density functional theory. Good agreement between the observed and predicted patterns of the spectral lines allowed for reliable identification. The experimental IR spectra were interpreted and discussed. The relative energies of the 12 isomers of N2-hydroxyisocytosine were calculated at the MP2/6-31G(d,p) and MP4//MP2/6-31G(d,p) levels. For six isomers of 1-methyl-N2-hydroxyisocytosine, the calculations were carried out at the MP2/6-31G(d,p) level. The anti form of the imino-oxo tautomer of N-hydroxyisocytosine and the syn form of the imino-oxo tautomer of 1-methyl-N2-hydroxyisocytosine were predicted to be the most stable.     

Thermodynamic properties of the NdBr3–MBr binary systems (M = Na,K)

Thixoforming involves heating different types of alloys to the semisolid state at high heating rates and forming in die-casting machines or conventional presses. At temperatures higher than the solidus and lower than the liquidus, the mush metal behaves like a high-viscosity thixotropic material. Therefore, determining the thermodynamic behavior of the solid-to-liquid transition is paramount to control thixoforming processes. This article describes a simple, novel experimental setup based on differential thermal analysis (DTA) for analyzing the phase transitions in an alloy heated using high heating rates typical of industrial applications. A365 alloy was chosen to demonstrate the effectiveness of the method as the phase transformations for this alloy in semisolid materials (SSM) processing are well understood. Samples were heated to 750 °C using constant linear heating rates of 5, 10, 15, 20, 25, 50, 75, 100 and 125 °C min in a Norax 25 kW 8 kHz induction furnace with an Omron E5CK temperature controller. AISI 316 austenitic stainless steel was used as the inert reference. Comparison of the results of DTA using the proposed method and the results of simulation with Thermo-Calc^® indicates that the proposed in situ DTA device and its method is suitable for analyzing phase transitions when high heating rates are used.