A pseudopotential study of the hydrogen bond in H2O·H2S,H2S·H2S and H2O·H2Se systems

Intermolecular potential energy curves for the hydrogen bonded systems H₂O·H₂S, H₂O·H₂Se and H₂S·H₂S were calculated with nonempirical pseudopotentials using optimized-in-molecules basis sets augmented by polarization functions. The H₂O·H₂O interaction energy curve has been also considered as a test case. The present results for H₂O·H₂S and H₂S·H₂S indicate much weaker intermolecular interactions than those found in previous ab initio calculations. The H₂O·H₂Se interaction was found to be quite similar to H₂O·H₂S.This work was partly supported by the Polish Academy of Sciences within the Project PAN-09, 7.1.1.1On leave from Quantum Chemistry Laboratory, Dept. of Chemistry, University of Warsaw, Pasteura 1, 02-093. Warsaw, Poland     

On the effect of the BSSE on intermolecular potential energy surfaces. Comparison of a priori and a posteriori BSSE correction schemes

A comparative study of geometrical parameters is performed on the complexes HF–HF, H₂O–H₂O, and HF–H₂O using 12 different basis sets at the RHF, MP2, and DFT (BLYP and B3LYP) levels of theory. The equilibrium geometries were obtained from uncorrected, a posteriori (counterpoise, CP) and a priori (Chemical Hamiltonian Approach, CHA) BSSE‐corrected potential energy surfaces. The calculation of equilibrium geometries using the CP and CHA schemes is described in details. The effect of the BSSE on various intermolecular parameters is discussed and the performance of the applied theoretical models is critically evaluated from the BSSE point of view. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 765–786, 2001     

Rotation–vibration energy levels from recent potential energy surfaces for the ground electronic states of NO2 and H2O

Using recently published potential energy surfaces, rovibrational energy levels are computed for the ground electronic states of H₂O and NO₂, using three‐dimensional discrete variable representation (DVR) algorithms. Calculations are presented for H₂O to demonstrate the accuracy of these algorithms. In the case of NO₂, comparisons with previous calculations are made as a test for the accuracy of the new potential energy surface. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 535–540, 2000     

Calculation of the Gibbs Energy of the Reaction OH-·(H2O)n-1 + H2O = OH-·(H2O)n by the Monte-Carlo Method

The energy, the Gibbs energy of the reaction OH^-·(H₂O) _{n- 1} + H₂O = OH^-·(H₂O) _n are calculated by the Monte-Carlo method with a large canonical ensemble for n = 1, ..., 20. The ion-waternonpair interaction potential was obtained by numerical fitting of calculated Gibbs energy and entropy of (H₂O)_n clusters (n = 1, ..., 5) to experimental ones. A good fit to experiment both of the internal energy and the Gibbs energy can be obtained in terms of a model allowing for nonpair interaction. It is shown that constructing an ion-water interaction potential without allowance for the entropy factor can lead to considerable errors in the Gibbs energy of cluster formation and in the nucleation rate.     

Intra‐ and intermolecular Se...X (X = Se,O) interactions in selenium‐containing heterocycles: 3‐benzoylimino‐5‐(morpholin‐4‐yl)‐1,2,4‐diselenazole

In the selenium‐containing heterocyclic title compound {systematic name: N‐[5‐(morpholin‐4‐yl)‐3H‐1,2,4‐diselenazol‐3‐ylidene]benzamide}, C₁₃H₁₃N₃O₂Se₂, the five‐membered 1,2,4‐diselenazole ring and the amide group form a planar unit, but the phenyl ring plane is twisted by 22.12 (19)° relative to this plane. The five consecutive N—C bond lengths are all of similar lengths [1.316 (6)–1.358 (6) Å], indicating substantial delocalization along these bonds. The Se...O distance of 2.302 (3) Å, combined with a longer than usual amide C=O bond of 2.252 (5) Å, suggest a significant interaction between the amide O atom and its adjacent Se atom. An analysis of related structures containing an Se—Se...X unit (X = Se, S, O) shows a strong correlation between the Se—Se bond length and the strength of the Se...X interaction. When X = O, the strength of the Se...O interaction also correlates with the carbonyl C=O bond length. Weak intermolecular Se...Se, Se...O, C—H...O, C—H...π and π–π interactions each serve to link the molecules into ribbons or chains, with the C—H...O motif being a double helix, while the combination of all interactions generates the overall three‐dimensional supramolecular framework.     

The investigation on the relative stability of different clusters of thiourea dioxide in water using gas phase quantum chemical calculations

The relative stability of different clusters of thiourea dioxide (TDO) in water is examined using gas phase quantum chemical calculations at the MP2 and B3LYP level with 6‐311++G(d,p) basis set. The possible equilibrium structures and other energetic and geometrical data of the thiourea dioxide clusters, TDO‐(H₂O)_n (n is the number of water molecules), are obtained. The calculation results show that a strong interaction exists between thiourea dioxide and water molecules, as indicated by the binding energies of the TDO clusters progressively increased by adding water molecules. PCM model is used to investigate solvent effect of TDO. We obtained a negative hydration energy of ?20.6 kcal mol^?1 and free‐energy change of ?21.0 kcal mol^?1 in hydration process. On the basis of increasing binding energies with adding water molecules and a negative hydration energy by PCM calculation, we conclude thiourea dioxide can dissolve in water molecules. Furthermore, the increases of the C? S bond distance by the addition of water molecules show that the strength of the C? S bonds is attenuated. We find that when the number of water molecules was up to 5, the C? S bonds of the clusters, TDO‐(H₂O)₅ and TDO‐(H₂O)₆ were ruptured. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Potential energy surface and thermochemistry for the direct gas phase reaction of germane and water

Gas phase reaction between germane GeH₄ and water H₂O was investigated at CCSD(T)/[aug-cc-pVTZ-pp for Ge + Lanl2dz for H and O]//MP2/6-31G(d,p) level. Only the hydrogen elimination channels are monitored. Within the energy range of 100 kcal/mol, we located nine equilibrium and six transition states on the potential energy surface (PES) of the Ge–O–H systems. GeH₄ reacts with H₂O exothermically (by 2.37 kcal/mol) without a barrier to form a non-planar complex GeH₄·H₂O which isomerizes to GeH₃OH·H₂ and H₂GeOH₂·H₂ with a barrier of 44.34 kcal/mol and 53.75 kcal/mol respectively. The first step of hydrogen elimination gives two non-planar species, GeH₃OH and H₂GeOH₂ but germinol GeH₃OH is found to be more stable. Further thermal decomposition reactions of GeH₃OH involving hydrogen elimination have been studied extensively using the same method. The final hydrogen elimination step gives HGeOH which can exist in cis and trans forms. As the trans form is more stable, only the trans form is considered on the potential energy surface (PES) of the reaction. The important thermochemical parameters (∆_rE_tot + ZPE), ∆_rH and ∆_rG for the H₂ elimination pathways are predicted accurately.     

Synthesis and Crystal Structures of Octahedral Metal Complexes containing the New Dianion [PhP(Se,O)Se‐Se(O,Se)PPh]2−

The synthesis and structures of three new compounds are reported. [Mg₂{PhP(Se,O)Se‐Se(O,Se)PPh}₂(thf)₇(H₂O)₃] ( 1 ), [Mg{PhP(Se,O)Se‐Se(O,Se)PPh}(thf)₃(H₂O)] ( 2 ), and [Mn{PhP(Se,O)Se‐Se(O,Se)PPh}(thf)₃(H₂O)] ( 3 ) were prepared by treatment of Woollins' reagent [PhP(Se)(μ‐Se)]₂ with the corresponding hydrated metal acetates.     

A model study of the intermolecular interactions of amino acids in aqueous solution: The glycine-water system

We have performed calculations of the glycine zwitterion surrounded by water molecules with the help of the mutually consistent field (MCF) method and perturbation theoretical expressions. Two different models for the hydration shell have been chosen, the glycine·6H₂O and glycine·12H₂O complexes, representing the most probable first and second solvation shell, respectively. To calculate the exchange and charge transfer energy contributions we have applied approximative expressions derived from perturbation theory for weakly overlapping subunits. For the sake of comparison we also calculated the interaction energy in the supermolecule approach for the smaller of the two solvation complexes. Furthermore, we have investigated the part of the potential energy surface which is determined by varying the lengths of the hydrogen bonds between glycine and water in the complex glycine·12H₂O using the electrostatic approach. The exchange energy contribution to the interaction energy for different points on the surface was approximated with the help of an analytical expression fitted to three directly calculated points. For the charge transfer energy a polynomial expansion of second order was established on the basis of five values, computed with the aid of the perturbation theoretical expression. To get a more detailed insight in the relatively strong hydrogen bonds between the water molecules and the ionic hydrophilic parts of glycineab initio model studies on NH ₄ ⁺ ·3H₂O and HCOO⁻·3H₂O systems are reported.     

Aromaticity and Chemical Bonding of Chalcogen-Bonded Capsules Featuring Enhanced Magnetic Anisotropy

We present a theoretical study of chalcogen bonded container capsules ( A_X+A_X ) where X=O, S, Se, and Te, and their encapsulation complexes with n-C₉H₂₀ (n-C₉H₂₀@ A_X+A_X ). Both Se and Te encapsulation complexes have significant experimental and computed binding energies, analogous to the hydrogen bonded counterparts, while the S and O capsules and their encapsulation complexes show only weak binding energies, which are attributed to different types of bonding: chalcogen S⋅⋅⋅N bonds for S-capsules and π–π stacking and weak hydrogen bonds for the O case. All A_X+A_X and C₉H₂₀@ A_X+A_X present unusually high magnetic anisotropies in their interiors. The ¹H NMR spectra of the encapsulation complexes display the proton signals of the encapsulated n-nonane highly upfield shifted, in agreement with the available experimental data for the Se capsule. We found that different factors contribute to the observed magnetic anisotropy of the capsule's interior: for the Te capsule the most important factor is Te's large polarizability; for the O analogue the inductive effects produced by the electronegative nature of the O and N heteroatoms; and for the S and Se capsules, the polarizability of the heteroatoms combines with electric field effects.     

New Polydentate Trimethylsilyl Chalcogenide Reagents for the Assembly of Polyferrocenyl Architectures

A series of polychalcogenotrimethylsilane complexes Ar(CH₂ESiMe₃)_n, (Ar=aryl; E=S, Se; n=2, 3, and 4) can be prepared from the corresponding polyorganobromide and M[ESiMe₃] (M=Na, Li). These represent the first examples of the incorporation of such a large number of reactive ?ESiMe₃ moieties onto an organic molecular framework. They are shown to be convenient reagents for the preparation of the polyferrocenylseleno‐ and thioesters from ferrocenoyl chloride. The synthesis, structures, and spectroscopic properties of the new silyl chalcogen complexes 1,4‐(Me₃SiECH₂)₂(C₆Me₄) (E=S, 1 ; E=Se, 2 ), 1,3,5‐(Me₃SiECH₂)₃(C₆Me₃) (E=S, 3 ; E=Se, 4 ) and 1,2,4,5‐(Me₃SiECH₂)₄(C₆H₂) (E=S, 5 ; E=Se, 6 ) and the polyferrocenyl chalcogenoesters [1,4‐{FcC(O)ECH₂}₂(C₆Me₄)] (E=S, 7 ; E=Se, 8 ), [1,3,5‐{FcC(O)ECH₂}₃(C₆Me₃)] (E=S, 9 ; E=Se, 10 ) and [1,2,4,5‐{FcC(O)ECH₂}₄(C₆H₂)] (E=S, 11 illustrated; E=Se, 12 ) are reported. The new polysilylated reagents and polyferrocenyl chalcogenoesters have been characterized by multinuclear NMR spectroscopy (¹H, ¹³C, ⁷⁷Se), electrospray ionization mass spectrometry and, for complexes 1 , 2 , 3 , 4 , 7 , 8 , and 11 , single‐crystal X‐ray diffraction. The cyclic voltammograms of complexes 7 – 11 are presented.     

3,3′‐Seleno­bis­(propionic acid)

In contrast to Se[CH₂C(O)OH]₂versus S[CH₂C(O)OH]₂, the title compound, Se[CH₂CH₂C(O)OH]₂ or C₆H₁₀O₄Se, is structurally quite similar to its sulfur analogue. The mol­ecule has twofold symmetry. The C—Se—C bond angle is 96.48 (8)° and the Se—C bond lengths are 1.9610 (14) Å. The shortest Se?O intermolecular distance is 3.5410 (11) Å. The O?O distances in the carboxyl­ic acid dimers are 2.684 (2) Å. The temperature dependence of the IR spectrum suggests tautomerism in the solid state.     

Zustandsbarogramme — Zustandsdiagramme durch Gesamtdruckmessungen

From phasebarograms it is possible to construct phase diagrams. We show on example of the system Bi/Se/O how follow ternary sections from the total pressure measurements in equilibrium and how we can attribute the barogram of the ternary region Bi₂Se₃/Bi₂O₂Se/Se to the binary system Bi₂Se₃/Se. The knowledge of the ternary system Bi/Se/I and its coexistent pressure courses allow to followT-p-x conditions for the chemical transport of phases from this system.

     

Electron energy loss spectroscopy (EELS) of H2O,H2S,H2Se and H2Te

Electronic excitation in H₂O, H₂S, H₂Se and H₂Te molecules has been studied by the EELS technique. Spectra of H₂S and H₂Se are remarkably similar with the 1b₁-nd transition most intense. The intensity of the first transition 1b₁-nsa₁ decreases through H₂O to H₂Se and this transition is absent in H₂Te. Transitions observed by EELS have been compared with optical absorption studies. A correlation diagram of the occupied and the excited states has been provided for these four molecules by making use of UVPES and EELS.     

SCF,MP2, and CEPA-1 calculations on the OH ‥ O hydrogen bonded complexes (H2O)2 and (H2O-H2CO)

Counterpoise corrected ab initio calculations are reported for (H₂O)₂ and H₂O-H₂CO. Geometry searches were done in the moment-optimized basis DZP' at the SCF, MP2, and CEPA-1 levels of theory, followed by more accurate single-point calculations in basis ESPB, which includes bondfunctions to saturate the dispersion energy. The final equilibrium binding energies obtained are ?4.7 ±0.3 kcal/mol for a near-linear (H₂O)₂ structure and ?4.6 ±0.3 kcal/mol for a strongly bent HOH ‥ OCH₂ structure. The energy difference between these systems is much smaller than in all previous ab initio work. Cyclic (C_2h) and bifurcated (C_2v) transition structures for (H₂O)₂ are located at 1.0 ±0.1 kcal/mol and 1.9 ±0.3 kcal/mol above the global minimum, respectively. A new partitioning scheme is presented that rigorously partitions the MP2 correlation interaction energy in intra and intermolecular (dispersion) contributions. These terms are large (up to 2 kcal/mol) but of opposite sign for most geometries studied and hence their overall effect upon the final structures is relatively small. The relative merits of the MP2 and CEPA-1 approaches are discussed are discussed and it is concluded that for economical reasons MP2 is to be preferred, especially for larger systems.     

SCF-CI studies of correlation effects on hydrogen bonding and ion hydration

Previous single-determinant Hartree-Fock studies on the equilibrium structures and stabilities of H₂ O, H₃ O⁺ as well as of the monohydrated ionic systems Li⁺ · H₂O, F^? · H₂O and the hydrogen bonded water dimer, H₂ O · HOH, are extended by large scale configuration interaction calculations including all the possible single and double excitations arising from the canonical set of Hartree-Fock molecular orbitals. The correlation energy effects on the equilibrium geometrical parameters of the systems under consideration are found to be quite small. The contributions of the correlation energy to the total binding energies of the weakly interacting composed systems are obtained to be of the order of 1 kcal/mole, leading to a considerable increase of the hydrogen bond strength in F^? · H₂O and H₂O · HOH and to a small decrease of the binding energy in Li⁺ · H₂ O. The observed strengthening of the hydrogen bonding interaction due to correlation is shown to be partly compensated by the change in the vibrational zero-point energy of the composed systems compared to the non-interacting subsystems. Approximate force constants corresponding to the intersystem vibrations in Li⁺ · H₂O, F^? · H₂ O, and H₂O · HOH are deduced from the calculated potential curve data on the SCF and the CI level of accuracy.     

SCF LCAO MO studies on the hydration of ions: The system F− · 2H2O

The energy surface of the dihydrated fluoride anion (F·2H₂O)^–1 is studied for a number of different geometry points near the equilibrium structure within the SCF LCAO MO framework, using an extended gaussian basis set to approximate the molecular wavefunctions. For the first and second hydration step of the fluoride anion the corresponding hydration energies are calculated to beB ₁ ^scf =24.1 kcal/mole andB ₂ ^SCF =20.8 kcal/mole (experimental measurements: 23.3 kcal/mole and 16.6 kcal/mole, respectively). The hydration energies and equilibrium bond distances obtained for the dihydrated fluoride anion (F·2H₂O)^– are compared with those found for the monohydrate (FHOH)^– and with corresponding results of the dihydrated lithium cation (Li · 2H₂O)⁺. The system (F·2H₂O)^– is taken as a very simple model to discuss some basic features of the hydration process of small ions and to study the influence of a negative ion on an adjacent hydrogen bond.We would like to thank our technical staff for valuable help in carrying out these calculations.     

An ab initio-based global potential energy surface for the SH3 system and full-dimensional state-to-state quantum dynamics study for the H2 + HS → H2S + H reaction

An accurate potential energy surface for the ground electronic state of SH₃ system has been constructed with 41,882 high level ab initio energy points and the neural network fitting method. The time-dependent wave packet method has been used to calculate the first state-to-state differential cross sections for the title reaction up to 1.2 eV in full dimensions, based on the reactant–product decoupling scheme. It is found that the majority of H₂S are produced in the ground vibrational state, with a large fraction of available energy for the reaction ending up as product translational motion. The differential cross sections at the threshold energy are dominated by a very narrow peak in the backward direction. With the increase of collision energy, the width of the angular distribution increases considerably, which is a typical feature of a direct reaction via abstract mechanism, similar to the H₂ + OH → H₂O + H reaction. © 2018 Wiley Periodicals, Inc.