Absolute calculations of acidity of C-substituted tetrazoles in solution

The CBS-QB3 method was used to calculate the gas-phase free energy difference between nine tetrazole derivatives and their anions, and the DPCM and CPCM continuum solvation methods were applied to calculate the free energy differences of solvation. The calculations were performed on both gas-phase and solvent-phase optimized structures. Absolute pKa calculations using the CPCM method and the gas-phase optimized structures yielded mean unsigned error of 0.4 pKa unit. The calculations were made with the routine settings implemented in Gaussian 98. The study is as accurate as the best reported so far for six carboxylic acids and phenols and, to our knowledge, the best reported for the acidities of heterocyclic compounds in solution.     

C-substituted diethyleneimides of acids of phosphorus

Conclusions We have obtained the hitherto undescribed 2,2-dimethyl- and 2-ethylethyleneimides of alkylphosphorous acids and, by addition of sulfur and selenium to these compounds, the corresponding ethyleneimldes of alkylthio-and selenophosphoric acids.Translated from Izvestiya Akademii Nauk SSSR, Seriya Khimicheskaya, No. 5, pp. 1143–1144, May, 1968.     

Application of natural bond orbital analysis to delocalization and aromaticity in C-substituted tetrazoles.

Energies of two tautomeric forms of 10 tetrazole derivatives substituted at C5 were established by DFT/B3LYP calculations carried out at the 6-311++G level. In each case the calculated energy of the 2H-tautomer was lower than that of the 1H. Furthermore, three geometric aromaticity indices of both forms were calculated, as were the values of nuclear independent nuclear shift and aromatic stabilization energy. The electronic properties were evaluated with the help of the natural bonding orbital theory. Following this a new pi-delocalization parameter, the root-mean square of pi-electron density localized on the atoms of the five-membered tetrazole ring, SDn, was introduced. It was concluded that the electronic delocalization can be described equally well by three different parameters: SDn, the extent of the transfer of electron density from the p(z) orbital of one nitrogen to the rest of the pi electron system, and population of two antibonding pi orbitals. Arguably, the information provided by the electronic parameters is similar to that contained in the geometric (structural) aromaticity indices except for tetrazole substituted by -BH(2).     

C-substituted macrocycles as candidates for radioimmunotherapy

The reaction between aryl aldehydes, the macrocyclic ligand 6-methyl-1,4,8,11-tetraazacyclotetradecane-6-amine (L1), and NaBH3CN produces the corresponding benzyl-substituted ligands in good yield. Copper(II) complexes of the ligands derived from salicylaldehyde (L2), p-hydroxybenzaldehyde (L4), and p-carboxybenzaldehyde (L5) were structurally characterized: [CuL2](ClO4)2.3H2O (monoclinic, P2(1)/c, a = 11.915(6) A, b = 13.861(2) A, c = 17.065(8) A, beta = 102.14(2) degrees, Z = 4); [CuL4](ClO4)2 (monoclinic, P2(1)/n, a = 9.550(3) A, b = 17.977(2) A, c = 14.612(4) A, beta 96.76(1) degrees, Z = 4), and [CuL4](ClO4)2 (monoclinic, P2(1)/n, a = 9.286(2) A, b = 11.294(1) A, c = 23.609(8) A, beta 93.68(1) degrees, Z = 4). Conjugation of several CuII complexes to a protein (bovine serum albumin) has been pursued with a view to the application of these macrocycles as bifunctional chelating agents in radioimmunotherapy.     

Aldehyde separation by polymer-supported oligo(ethyleneimines)

Chlorosulfonated styrene (10%) divinylbenzene resin beads reacted with an excess of ethylenediamine (EDA), diethylenetriamine (DETA), and triethylenetetramine (TETA) to give the corresponding sulfonamides with pendant oligo(ethyleneimines). The resulting modified resins are useful in the separation of aldehydes from hydrocarbon mixtures. Sorption of aldehydes occurs through formation of both Schiff base and five-membered (imidazoline) rings. Sorbed aldehydes can readily be stripped from the resins by treating with dilute acid solutions. Since the sulfamide bond has a reasonable stability toward acid-base hydrolysis, the loaded resins can be regenerated and recycled by simple acid-base washings, without losing their activity. In the present study, sorption and desorption kinetics of acetaldehyde, benzaldehyde, and salicylaldehyde have been investigated under different conditions. The aldehyde sorption obeys second-order kinetics. The method presented is applicable for all aromatic aldehydes. However, in the case of aliphatic aldehydes carrying an α-hydrogen, aldol condensation products form in solution. So aliphatic aldehydes and their aldol products are sorbed together by the resins. This limits the recovery of aliphatic aldehydes. Consequently, the resins described are cost effective sorbents for the removal and recovery of aromatic aldehydes from various mixtures. © 1997 John Wiley & Sons, Inc. J Polym Sci A: Polym Chem 35 : 2857–2864, 1997     

Structural and photophysical studies of trans-AB(2)C-substituted porphyrin ligands and their zinc and copper complexes

trans-AB(2)C porphyrins with A = C(6)H(4)-COOR, C = C(6)H(4)-NX(2) and B = C(6)H(5) (R = CH(3), H; X = O, H) have been synthesised by a rational high-yield procedure (1a-1d) and their zinc(ii) and copper(ii) complexes have been prepared (2a-2d, 3a-3d ).1a, 2a .THF and 3a display different distortions of the porphyrin core as shown by single crystal X-ray crystallography and NSD analyses. The Soret and Q bands of free-base and metalated porphyrins with mixed electron donating and withdrawing substituents (NH(2)/COOR) are red-shifted as are the corresponding emission bands of free-base and zinc porphyrins. The electronic asymmetry revealed by spectrocopy is rationalised by DFT calculations.     

Syntheses of C-substituted icosahedral dicarbaboranes bearing the 8-dioxane-cobalt bisdicarbollide moiety

The treatment of 1,2-, 1,7- and 1,12-carbaborane lithiated isomers with [3,3′-Co-8-(CH₂CH₂O)₂-(1,2-C₂B₉H₁₀)-(1′,2′-C₂B₉H₁₁)] (1) at molar ratios 1:1 or 1:2 at room temperature in THF leads generally to the formation of a series of orange double-cluster mono and dianions. These were characterized by NMR and MS methods as [1′′-X-1′′,2′′-closo-C₂B₁₀H₁₁]⁻, [2]⁻; [1′′-X-1′′,7′′-closo-C₂B₁₀H₁₁]⁻, [3]⁻ and [1′′-X-1′′,12′′-closo-C₂B₁₀H₁₁], [4]⁻ for the monoanions, whereas [1′′,2′′-X₂-1′′,2′′-closo-C₂B₁₀H₁₀]²⁻, [2]²⁻; [1′′,7′′-X₂-1′′,7′′-closo-C₂B₁₀H₁₀]²⁻, [3]²⁻; and [1′′,12′′-X₂-1′′,12′′-closo-C₂B₁₀H₁₀]²⁻, [4]²⁻ for the dianions (where X = 3,3′-Co-8-(CH₂CH₂O)₂-(1,2-C₂B₉H₁₀)-1′,2′-(C₂B₉H₁₁)). Moreover, these borane-cage subunits can be easily modified via attaching variable substituents onto cage carbon and boron vertices, which makes these compounds structurally flexible potential candidates for BNCT of cancer and HIV-PR inhibition.     

Ab initio study of exchange coupling for the consistent understanding of the magnetic ordering at room temperature in V[TCNE] x

We report quantum chemical calculations providing the exchange coupling constants of the V[TCNE]₂ model system, describing the amorphous room temperature molecular magnet V[TCNE] _x (TCNE = tetracyanoethylene, x ~ 2). The geometry is optimized for the ideal lattice using density functional theory (DFT) calculations with periodic boundary conditions. Broken-symmetry DFT calculations indicate antiparallel spin alignment resulting in ferrimagnetic ordering, but heavily overestimate the value of the exchange coupling. Better estimates of the exchange coupling parameters between the V(II) ion and the [TCNE]^? anionic radical are obtained by means of multiconfigurational calculations performed on smaller molecular models cut from the optimized crystal lattice. Complete active space self-consistent field and multireference second-order perturbation theory calculations provide the sign and the strength of the nearest-neighbor as well as next-nearest-neighbor interactions along all three crystallographic directions. We are able to explain also intuitively the mechanism for antiferromagnetic spin coupling in terms of the superexchange pathways, discussing the role of the main four types of contributions to superexchange. Moreover, we clarify the influence of the transition metal ion on the strength of the exchange interaction and on the critical temperature for long-range ferrimagnetic ordering.     

Thermodynamical and structural study of protactinium(V) oxalate complexes in solution

The complexation of protactinium(V) by oxalate was studied by X-ray absorption spectroscopy (XAS), density functional theory (DFT) calculations, capillary electrophoresis coupled with inductively coupled plasma mass spectrometry (CE-ICP-MS) and solvent extraction. XAS measurements showed unambiguously the presence of a short single oxo-bond, and the deduced structure agrees with theoretical calculations. CE-ICP-MS results indicated the formation of a highly charged anionic complex. The formation constants of PaO(C(2)O(4))(+), PaO(C(2)O(4))(2)(-), and PaO(C(2)O(4))(3)(3-) were determined from solvent extraction data by using protactinium at tracer scale (C(Pa) < 10(-10) M). Complexation reactions of Pa(V) with oxalate were found to be exothermic with relatively high positive entropic variation.     

Verification of DFT-predicted hydrogen storage capacity of VC3H3 complex using molecular dynamics simulations

Density functional theory (DFT) and Fourth‐order Möller–Plesset (MP4) perturbation theory calculations are performed to examine the possibility of hydrogen storage in V‐capped VC₃H₃ complex. Stability of bare and H₂ molecules adsorbed V‐capped VC₃H₃ complex is verified using DFT and MP4 method. Thermo‐chemistry calculations are carried out to estimate the Gibbs free corrected averaged H₂ adsorption energy which reveals whether H₂ adsorption on V‐capped VC₃H₃ complex is energetically favorable, at different temperatures. We use different exchange and correlation functionals employed in DFT to see their effect on H₂ adsorption energy. Molecular dynamic (MD) simulations are performed to confirm whether this complex adsorbs H₂ molecules at a finite temperature. We elucidate the correlation between H₂ adsorption energy obtained from density functional calculations and retaining number of H₂ molecules on VC₃H₃ complex during MDs simulations at various temperatures. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Computational study on the reaction mechanism of hydrosilylation of carbonyls catalyzed by high-valent rhenium(V)-di-oxo complexes

Density functional theory (DFT) calculations have been performed to elucidate the reaction mechanism of hydrosilylation of carbonyl compounds catalyzed by high-valent rhenium(V)-di-oxo complexes ReO2I(PR3)2 (R = Me, Ph). The calculations suggest that the most favorable mechanism involves the rate-determining dissociative [2 + 2] addition of the Si-H bond across a Re=O bond to form a Re(V) hydrido siloxy intermediate; this is followed by carbonyl coordination, reduction of the carbonyl, rearrangement, and final intramolecular nucleophilic attack from the alkoxy group to the silyl center (dissociative retro-[2 + 2] addition). It was also found that the additional oxo ligand in the ReO2I(PR3)2 complexes promotes the [2 + 2] addition across the rhenium-oxo bond both kinetically and thermodynamically, as compared to the neutral rhenium(V)-mono-oxo complex ReOCl3(PMe3)2. The effect of different silanes on the [2 + 2] addition barriers is also studied.     

Computation of multiphase equilibrium phenomena with an equation of state

This paper presents a general approach to multiphase equilibrium computations using an equation of state. Various types of calculations are described: flash calculations, saturation-pressure and saturation-temperature calculations, and calculations of the feed composition corresponding to a given saturation pressure and temperature.A stagewise procedure for flash calculations using the QNSS method allows efficient computation of phase equilibria for systems with any number of phases. The initial guess of the composition of an extra phase which is required to start the iteration is obtained from an analysis of the Gibbs-energy surface.Newton's method is used to construct multiphase boundaries on various types of diagrams (pressure-temperature, pressure-composition and temperature-composition). Multiphase critical points on these boundaries are also estimated by interpolation.The algorithms developed are tested for a typical reservoir-oil-CO₂ mixture which exhibits three-phase liquid 1-liquid 2-vapor (L₁L₂V) separation. The predictions obtained using the Peng-Robinson of state are satisfactory.     

51V solid-state NMR and density functional theory studies of vanadium environments in V(V)O2 dipicolinic acid complexes

(51)V solid-state NMR and density functional theory (DFT) investigations are reported for a series of pentacoordinate dioxovanadium(V)-dipicolinate [V(V)O(2)-dipicolinate] and heptacoordinate aquahydroxylamidooxovanadium(V)-dipicolinate [V(V)O-dipicolinate] complexes. These compounds are of interest because of their potency as phosphatase inhibitors as well as their insulin enhancing properties and potential for the treatment of diabetes. Experimental solid-state NMR results show that the electric field gradient tensors in the V(V)O(2)-dipicolinate derivatives are affected significantly by substitution on the dipicolinate ring and range from 5.8 to 8.3 MHz. The chemical shift anisotropies show less dramatic variations with respect to the ligand changes and range between -550 and -600 ppm. To gain insights on the origins of the NMR parameters, DFT calculations were conducted for an extensive series of the V(V)O(2)- and V(V)O-dipicolinate complexes. To assess the level of theory required for the accurate calculation of the (51)V NMR parameters, different functionals, basis sets, and structural models were explored in the DFT study. It is shown that the original x-ray crystallographic geometries, including all counterions and solvation water molecules within 5 A of the vanadium, lead to the most accurate results. The choice of the functional and the basis set at a high level of theory has a relatively minor impact on the outcome of the chemical shift anisotropy calculations; however, the use of large basis sets is necessary for accurate calculations of the quadrupole coupling constants for several compounds of the V(V)O(2) series. These studies demonstrate that even though the vanadium compounds under investigations exhibit distorted trigonal bipyramidal coordination geometry, they have a "perfect" trigonal bipyramidal electronic environment. This observation could potentially explain why vanadate and vanadium(V) adducts are often recognized as potent transition state analogs.     

Determination of the spin-lattice relevant for the quaternary magnetic oxide Bi4Cu3V2O14 on the basis of tight-binding and density functional calculations

The quaternary magnetic oxide Bi4Cu3V2O14 consists of Cu4O8 triple chains made up of corner-sharing CuO4 square planes. To determine its spin-lattice, the spin exchange interactions of Bi4Cu3V2O14 were evaluated by performing a spin dimer analysis based on tight-binding calculations and a mapping analysis based on first principles density functional theory calculations. Both calculations show that the spin-lattice of Bi4Cu3V2O14 is not an antiferromagnetically coupled diamond chain, which results from an idealized view of the structure of the Cu4O8 triple chain and a neglect of super-superexchange interactions. The correct spin-lattice is an antiferromagnetic chain made up of antiferromagnetic linear trimers coupled through their midpoints via super-superexchange interaction, which predicts that Bi4Cu3V2O14 has an antiferromagnetic spin ground state and has no spin frustration, both in agreement with experiment.     

Liquid—liquid equilibria of ternary systems of 1-hexene/hexane and extraction solvents

Olefins and paraffins are important basic chemical raw materials with so similar molecular structure and volatility that their separation is a difficult and energy-consuming process. Liquid—liquid equilibrium (LLE) data were determined for ternary systems of: 1-hexene + hexane + solvent at 25°C, 30°C, and 35°C, and N-formylmorpholine (NFM), N-methyl-2-pyrrolidone (NMP), γ-butyrolactone (GBL), and 1-methylimidazole (1-MEI) as the solvents studied. Liquid—liquid equilibrium data for each system were correlated to the NRTL (Non-Random Two Liquid) equation and the interaction parameters presented. The calculated results agree well with the experimental data.     

Analysis of approximations and errors in equations of motion method calculations

We analyze a number of fundamental questions associated with the use of a finite one-particle orbital basis in equations of motion (EOM) method calculations of excitation energies etc., of atomic and molecular systems. This approximation yields an approximate n_e-electron ground state and say, N excited states, while there are (N + 1)² different possible basis operators for EOM calculations. We show that sets of at most 2N basis operators can contribute to the EOM calculations. Any set of 2N basis operators, satisfying certain conditions, provides the exact EOM energies which are equivalent to complete configuration interaction results within the same orbital basis. We investigate the use of particle-particle shifting operators which are not employed in EOM calculations in model calculations on He with operator bases smaller than the complete 2V to consider the convergence of the expansion. The dependence of EOM calculations on the quality of the approximate ground state wavefunction is studied through calculations for Be where additional support is provided for the frequent need for multiconfigurational zeroth order reference functions (as corrected perturbatively). Excited state EOM wavefunctions from EOM calculations are shown to not necessarily be orthogonal to either the exact or approximate ground state wavefunction, suggesting implications in the use of EOM methods to evaluate excited state properties. The He and Be examples and a simple two-level problem are also utilized to illustrate questions concerning the use of the EOM equations to obtain an iteratively improved ground state wavefunction.     

Has AuF7 been made?

Quantum chemical calculations at DFT (BP86, B3LYP, BHLYP), MP2, CCSD, and CCSD(T) levels have been carried out on various fluoro complexes of gold in oxidation states +V through +VII to evaluate the previously claimed existence of AuF7. The calculations indicate clearly that elimination of F2 from AuF7 is a strongly exothermic reaction with a low activation barrier. This is inconsistent with the reported stability of AuF7 up to room temperature. A reported experimental vibrational frequency at 734 cm(-1) for AuF7 could not be verified computationally. It is concluded that the reported observation of AuF7 was probably erroneous. As the calculations indicate also an extremely large electron affinity and little stability for AuF6, Au(V) remains the highest well-established gold oxidation state.     

Application of the MC SCF method to the π → π* excitation energies of ethylene

The MC SCF method is employed to calculate the N → T and N → V π → π^* vertical excitation energies of ethylene. To obtain accurate excitation energies it is found to be necessary to utilize an expanded valence space containing two π and two π^* orbitals. Relatively small MC SCF calculations, allowing at most one-electron excitations from the sigma space, are found to yield excitation energies and spatial extents of the excited states in excellent agreement with the predictions of large multi-reference or iterative-natural-orbital CI calculations. These results show that within an MC SCF framework σ-σ correlation is unimportant for describing the π → π^* processes. We also conclude that the neglect of the effects of unlinked cluster terms in some of the CI calculations may have introduced small, but important, errors in the excitation energies and predictions of the spatial extent of the V state.     

Evidence for quantum interference in SAMs of arylethynylene thiolates in tunneling junctions with eutectic Ga-In (EGaIn) top-contacts

This paper compares the current density (J) versus applied bias (V) of self-assembled monolayers (SAMs) of three different ethynylthiophenol-functionalized anthracene derivatives of approximately the same thickness with linear-conjugation (AC), cross-conjugation (AQ), and broken-conjugation (AH) using liquid eutectic Ga-In (EGaIn) supporting a native skin (~1 nm thick) of Ga(2)O(3) as a nondamaging, conformal top-contact. This skin imparts non-Newtonian rheological properties that distinguish EGaIn from other top-contacts; however, it may also have limited the maximum values of J observed for AC. The measured values of J for AH and AQ are not significantly different (J ≈ 10(-1)A/cm(2) at V = 0.4 V). For AC, however, J is 1 (using log averages) or 2 (using Gaussian fits) orders of magnitude higher than for AH and AQ. These values are in good qualitative agreement with gDFTB calculations on single AC, AQ, and AH molecules chemisorbed between Au contacts that predict currents, I, that are 2 orders of magnitude higher for AC than for AH at 0 < |V| < 0.4 V. The calculations predict a higher value of I for AQ than for AH; however, the magnitude is highly dependent on the position of the Fermi energy, which cannot be calculated precisely. In this sense, the theoretical predictions and experimental conclusions agree that linearly conjugated AC is significantly more conductive than either cross-conjugated AQ or broken conjugate AH and that AQ and AH cannot necessarily be easily differentiated from each other. These observations are ascribed to quantum interference effects. The agreement between the theoretical predictions on single molecules and the measurements on SAMs suggest that molecule-molecule interactions do not play a significant role in the transport properties of AC, AQ, and AH.     

First‐Principles Calculations on Electronic Structures of N/V‐Doped and N‐V‐Dodoped Anatase TiO2 (101) Surfaces

The energetic and electronic properties of N/V‐doped and N‐V‐codoped anatase TiO₂ (101) surfaces are investigated by first‐principles calculations, with the aim to elucidate the relationship between the electronic structure and the photocatalytic performance of N‐V‐codoped TiO₂. Several substitutional and interstitial configurations for the N and/or V impurities in the bulk phase and on the surface are studied, and the relative stability of different doping configurations is compared by the impurity formation energy. Systematic calculations reveal that N and V impurities can be encapsulated by TiO₂ to form stable structures as a result of strong N‐V interactions both in the bulk and the surface model. Through analyzing and comparing the electronic structures of different doping systems, the synergistic doping effects are discussed in detail. Based on these discussions, we suggest that N_OV_Ti codoping cannot only narrow the band gap of anatase TiO₂, but also forms impurity states, which are propitious for the separation of photoexcited electron–hole pairs. In the case of N_OV_Ti‐codoped TiO₂ (101) surfaces, this phenomenon is especially prominent. Finally, a feasible synthesis route for N_OV_Ti codoping into anatase TiO₂ is proposed.