An electronic structure theory investigation of the physical chemistry of the intermolecular complexes of cyclopropenylidene with hydrogen halides

The proton accepting and donating abilities of cyclopropenylidene (c‐C₃H₂) on its complexation with hydrogen halides H? X (X = F, Cl, Br) are analyzed using density‐functional theory with three functionals (PBE0, B3LYP, and B3LYP‐D) and benchmarked against second‐order Møller–Plesset (MP2) theory. Standard signatures including, inter alia, dipole moment enhancement, charge transfer from the carbenic lone pair to the antibonding σ*(H? X) orbital, and H? X bond elongation are examined to ascertain the presence of hydrogen bonding in these complexes. The latter property is found to be accompanied with a pronounced red shift in the bond stretching frequency and with a substantial increase in the infrared intensity of the band on complex formation. The MP2/aug‐cc‐pVTZ c‐C₃H₂···H? F complex potential energy surface turns out to be an asymmetric deep single well, while asymmetric double wells are found for the c‐C₃H₂···H? Cl and c‐C₃H₂···H? Br complexes, with an energy barrier of 4.1 kcal mol^?1 for proton transfer along the hydrogen bond in the latter complex. Hydrogen‐bond energy decomposition, with the reduced variational space self‐consistent field approach, indicates that there are large polarization and charge‐transfer interactions between the interacting partners in c‐C₃H₂···H? Br compared to the other two complexes. The C···H bonds are found to be predominantly ionic with partial covalent character, unveiled by the quantum theory of atoms in molecules. The present results reveal that the c‐C₃H₂ carbene divalent carbon can act as a proton acceptor and is responsible for the formation of hydrogen bonds in the complexes investigated. © 2012 Wiley Periodicals, Inc.     

复合物HNO···H2O2内N-H···O蓝移氢键的理论研究

A theoretical study on the blue-shifted H-bond N-H…O and red-shifted H-bond O-H…O in the complexHNO…H_2O_2 was conducted by employment of both standard and counterpoise-corrected methods to calculate thegeometric structures and vibrational frequencies at the MP2/6-31G(d),MP2/6-31 G(d,p),MP2/6-311 q G(d,p),B3LYP/6-31G(d),B3LYP/6-31 G(d,p) and B3LYP/6-311 G(d,p) levels.In the H-bond N-H…O,the calcu-lated blue shift of N-H stretching frequency is in the vicinity of 120 cm~(-1) and this is indeed the largest theoreticalestimate of a blue shift in the X-H…Y H-bond ever reported in the literature.From the natural bond orbital analy-sis,the red-shifted H-bond O-H…O can be explained on the basis of the dominant role of the hyperconjugation.For the blue-shifted H-bond N-H…O,the hyperconjugation was inhibited due to the existence of significant elec-tron density redistribution effect,and the large blue shift of the N-H stretching frequency was prominently due tothe rehybridization of sp~n N-H hybrid orbital.     

Photochemical Generation of Chiral N,B,X‐Heterocycles by Heteroaromatic C−X Bond Scission (X=S,O) and Boron Insertion

Chiral organoboron compounds with a chelate backbone and mesityl/heterocycle substituents (thienyl, furyl, and derivatives thereof) undergo a quantitative phototransformation that yields rare, chiral N,B,X‐containing heterocycles, such as base‐stabilized 1,2‐thiaborinines and 1,2‐oxaborinines. Boriranes were observed as intermediates in some of these transformations. The oxaborinines display further reactivity, generating 4a,12b‐dihydrobenzo[h][1,2]oxaborinino[4,3‐f]quinolines through a sequential conrotatory electrocyclization and a [1,5]‐H shift. The N,B,X‐containing heterocycles display strong blue‐green to orange‐red emission in the solid state. Combined DFT//CASP2T calculations suggest that a common biradical intermediate is responsible for the formation of these compounds as well as their interconversion.     

Photochemical Generation of Chiral N,B,X‐Heterocycles by Heteroaromatic C−X Bond Scission (X=S,O) and Boron Insertion

Chiral organoboron compounds with a chelate backbone and mesityl/heterocycle substituents (thienyl, furyl, and derivatives thereof) undergo a quantitative phototransformation that yields rare, chiral N,B,X‐containing heterocycles, such as base‐stabilized 1,2‐thiaborinines and 1,2‐oxaborinines. Boriranes were observed as intermediates in some of these transformations. The oxaborinines display further reactivity, generating 4a,12b‐dihydrobenzo[h][1,2]oxaborinino[4,3‐f]quinolines through a sequential conrotatory electrocyclization and a [1,5]‐H shift. The N,B,X‐containing heterocycles display strong blue‐green to orange‐red emission in the solid state. Combined DFT//CASP2T calculations suggest that a common biradical intermediate is responsible for the formation of these compounds as well as their interconversion.     

Ab initio and AIM studies on typical π‐type and pseudo‐π‐type halogen bonds: Comparison with hydrogen bonds

Series of typical π‐type and pseudo‐π‐type halogen‐bonded complexes B ··· ClY and B ··· BrY and hydrogen‐bonded complex B ··· HY (B = C₂H₄, C₂H₂, and C₃H₆; Y = F, Cl, and Br) have been investigated using the MP2/aug‐cc‐pVDZ method. A striking parallelism was found in the geometries, vibrational frequencies, binding energies, and topological properties between B ··· XY and B ··· HY (X = Cl and Br). It has been found that the lengths of the weak bond d(X ··· π)/d(H ··· π), the frequencies of the weak bond ν(X ··· π)/ν(H ··· π), the frequency shifts Δν(X? Y)/Δν(H? Y), the electron densities at the bond critical point of the weak bonds ρ_c(X ··· π)/ρ_c(H ··· π), and the electron density changes Δρ_c(X? Y)/Δρ_c(H? Y) could be used as measures of the strengths of typical π‐type and pseudo‐π‐type halogen/hydrogen bonds. The typical π‐type and pseudo‐π‐type halogen bond and hydrogen bond are noncovalent interactions. For the same Y, the halogen bond strengths are in the order B ··· ClY < B ··· BrY. For the same X, the halogen bond strength decreases according to the sequence F > Cl > Br that is in agreement with the hydrogen bond strengths B ··· HF > B ··· HCl > B ··· HBr. All of these typical π‐type and pseudo‐π‐type hydrogen‐bonded and halogen‐bonded complexes have the “conflict‐type” structure. Contour maps of the Laplacian of π electron density indicate that the formation of B ··· XY halogen‐bonded complex and B ··· HY hydrogen‐bonded complex is very similar. Charge transfer is observed from B to XY/HY and both the dipolar polarization and the volume of the halogen atom or hydrogen atom decrease on B ··· XY/B ··· HY complex formation. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

The synthesis of six‐membered P‐heterocycles with sterically demanding substituent on the phosphorus atom

The ring enlargement of 1‐(2,4,6‐trialkylphenyl)2,5‐dihydro‐1H‐phosphole oxides ( 1 ) via 6,6‐dichloro‐3‐Phosphabicyclo[3.1.0]hexanes ( 2 ) afforded the double‐bond isomers of 1,2‐dihydrophosphinine oxides ( 3 ). Catalytic hydrogenation of the isomeric 1‐(di‐tert‐butyltolyl)‐1,2‐dihydrophosphinine oxides ( 3a ) gave the diastereomers of phosphinane oxide ( 4 ), while that of the 1‐(tri‐isopropylphenyl) isomers ( 5 ) led predominantly to phospholane oxides ( 6 ) formed by ring contraction. © 2001 John Wiley & Sons, Inc. Heteroatom Chem 12:528–533, 2001     

Substituted 1,3,2,4‐benzodithiadiazines: Novel derivatives,by‐products,and intermediates*

The synthesis of the title compounds 1 by 1 : 1 condensation of Ar NSNSiMe₃ 2 with SCl₂ followed by intramolecular ortho‐cyclization of each [Ar NSN S Cl] intermediate is complicated by further reaction of 1 with SCl₂ to give Herz salts 3 . With the 2 :SCl₂ ratio of 2:1, the formation of by‐products 3 is reduced and novel compounds 1 are accessible. With ortho‐I containing starting material 2j , the parent compound 1s is obtained as the result of an unexpected I, not H, substitution. The rate of the 1 + SCl₂ reaction depends upon a substituent's position, and the minor 8‐R isomers 1l,p (R = Br, I) are isolated for the first time from mixtures with the major 6‐R isomers due to reduced reactivity toward SCl₂. The synthesized compounds 1–3 are characterized by multinuclear (including nitrogen) NMR and X‐ray crystallography. According to the X‐ray diffraction data, 1j (6‐Br) and 1k (7‐Br) derivatives are planar, whereas 1i (5‐Br) and 1l (8‐Br) are bent along the S1···N4 line by ∼5° and ∼4°, respectively, and the 1r (7‐OCH₃) derivative is planar in contrast to the known 5‐OCH₃ isomer, which possesses a significantly folded heterocycle. The distortion of the planar geometry of some compounds 1 is interpreted in terms of a pseudo‐Jahn‐Teller effect as the result of π‐highest occupied molecular orbital (HOMO)  σ*‐(LUMO) lowest unoccupied molecular orbital + 1 mixing in a planar conformation. The 2p compound is the first structurally defined Ar–N = S = N–SiMe₃ azathiene. The compound Ar–N = S = N–S–NH‐Ar 6 modeling the aforementioned intermediate has been isolated and structurally characterized. We describe the attempts to synthesize compounds 1 from 2‐aminobenzenethiols and (SN)₄ and from salts 3 and Me₃SiN₃, and we discuss the reaction pathways. © 2001 John Wiley & Sons, Inc. Heteroatom Chem 12:563–576, 2001     

Theoretical Studies on the Halogen Bond between Propargyl Radical and Dihalogen Molecules

MP2/aug‐cc‐pVDZ calculations are carried out on the geometries, vibrational frequencies, interaction energies and topological properties for the π‐type halogen‐bonded complexes between propargyl radical and dihalogen molecules ClF, BrF and BrCl. There are two kinds of geometries: complex ( a ) involves the interaction between the X (X=Cl, Br) atom and the midpoint of C(1) –C(2) bond, complex ( b ) involves the interaction between the X atom and C(3) atom. The lengths of the halogen bond, the frequencies of the halogen bond, the elongation extent of the X–Y (XY=ClF, BrF, BrCl) bond, topological parameters at the BCPs of the halogen bond and X–Y bond are all well consistent with the interaction energies. The interaction of complex ( a ) is stronger than that of complex ( b ); the interaction of propargyl···BrF is stronger than that of propargyl···ClF and propargyl···BrCl. For the complexes ( a ) and ( b ), the charge transfer is observed from propargyl radical to XY, the atomic energy, the dipolar polarization, and the volume of the halogen atom X decrease upon complex formation.     

Theoretical study of aromaticity in small hydrogen and metal cation clusters X (X=H,Li, Na,K, and Cu)

Five cation clusters X (X=H, Li, Na, K, and Cu) with two possible isomers, i.e., regular trigonal structure (D_3h) and linear structure (D_∞h), have been investigated using four methods: B3LYP, B3PW91, MP2, CCSD(T) and basis set 6‐311+G(3df). The calculations show that only the regular trigonal structure (D_3h) is stable. The related neutral clusters X₃Cl (X=H, Li, Na, K, and Cu) are also investigated using two methods: B3LYP, MP2, and basis set 6‐311+G(3df). For H₃Cl species, there is no a stable structure to be found. For other four X₃Cl (X=Li, Na, K, and Cu) species, there are two stable isomers, for which the bidentate structures (C_2v‐1) [see Fig. 1 (d)] are global minima. According to the general criteria for aromaticity including resonance energy (RE) and nucleus‐independent chemical shift (NICS), the five trigonal isomers exhibit a higher degree of aromaticity. Molecular orbital analysis reveals that the five trigonal X(X=H, Li, Na, K, and Cu) isomers possess only σ‐aromaticity originating from s orbitals. For the Cu ring the d orbitals do not play a significant role in the electron delocalization effects. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Density Functional Theory Study of Hydrogen Bonding Dimers with 4-Pyridinecarboxylic Acid Hydrazine

The H‐bonding dimers of 4‐pyridinecarboxylic acid hydrazine were studied using density functional theory (DFT) at B3LYP/6‐311++G** level. The results showed that the most stable dimer D1 had two same linear N H···O hydrogen bonds, and the interaction energy between them was 51.038 kJ·mol⁻¹ which was corrected by the basis set superposition error and zero‐point. The stretching vibration frequency of N H bond had a red shift because of the hydrogen bonds. The natural bond orbital analysis showed that each N H···O hydrogen bond in D1 had the biggest interaction stabilization energy of 69.078 kJ·mol⁻¹. Thermodynamic analysis indicated that the formation process of D1 was exothermic and spontaneous at low and room temperatures.     

Tieftemperaturuntersuchung von Wasserstoffbrückenbindungen in Lithiumtetrahydroxoborat mit Raman‐Spektroskopie,Röntgen‐ und Neutronenbeugung (Li11B(OD)4)

Low Temperature Investigation of Hydrogen Bridge Bonds in Lithium Tetrahydroxoborate by Raman Spectroscopy, X‐Ray and Neutron Diffraction (Li¹¹B(OD)₄) Low temperature Raman spectroscopic measurements on isotopically diluted Li¹¹B(OH)₄ with 8 % D and Li¹¹B(OD)₄ with 8 % H reveal four crystallographically different hydrogen bridge bonds. With decreasing temperatures beginning at ～50 K measured down to ～10 K the stretching modes of the hydroxide ions shift to higher wave numbers. For the strongest bond O–D···O the frequency shift is 16 cm^?1and for the weakest 7 cm^?1. For O–H···O the maximum in the frequency shift is 22 cm^?1. X‐ray single crystal (LiB(OH)₄) and neutron powder diffraction (Li¹¹B(OD)₄) data result in bond lengths for the four hydroxide ions in the range of 0.943 (3) Å ≤ d(O–D) ≤ 0.974 (3) Å. The value of the effect of inversion of the stretching mode frequencies seems to correlate with the strength of the hydrogen bridge bonds and is found to be different for the two isotopes H and D in this compound.     

Electron localization function and electron localizability indicator applied to study the bonding in the peroxynitrous acid HOONO

The ground‐state electronic structure of peroxynitrous acid (HOONO) and its singlet biradicaloid form (HO ··· ONO) have been studied using topological analysis of the electron localization function (ELF), together with the electron localizability indicator (ELI‐D), at the DFT (B3LYP, M05, M052X, and M06), CCSD, and CASSCF levels. Three isomers of HOONO (cis‐cis, cis‐perp, and trans‐perp) have been considered. The results show that from all functionals applied, only B3LYP yields the correct geometrical structure. The ELF and ELI‐D‐topology of the O? O and central N? O bonds strongly depends on the wave function used for analysis. Calculations carried out at CAS (14,12)/aug‐cc‐pVTZ//CCSD(T)/aug‐cc‐pVTZ level reveal two bonds of the charge‐shift type: a protocovalent N? O bond with a basin population of 0.82–1.08e, and a more electron depleted O? O bond with a population of 0.66–0.71e. The most favorable dissociation channel (HOONO → HO + ONO) corresponds to breaking of the most electron‐deficient bond (O? O). In the case of cis‐ and trans‐HO ··· ONO, the ELF, ELI‐D, and electron density fields results demonstrate a closed‐shell O ··· O interaction. The α‐spin electrons are found mainly (0.64e) in the lone pairs of oxygen V_{i = 1,2} (O) from the OH group. The β‐spin electrons are delocalized over the ONO group, with the largest concentration (0.34e) on the lone pair of nitrogen V(N). © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Theoretical Study on the Structures and Properties of Hydrogen Bonding Complexes between Diazines and Water

Density functional theory B3LYP method and second-order Moller-Plesset perturbation theory MP2 method were employed to obtain the optimized geometries of the ground state and interaction energy for diazines and water complexes. The results show that the ground state complexes have strong hydrogen bonding interaction with -20.99, -16.73 and -15.31 kJ/mol after basis set superposition error and zero-point vibration energy correction for pyridazine-water, pyrimidine-water and pyrazine-water, respectively, and large red-shift for the symmetric H-O stretching vibration frequencies due to the formation of N…H-O hydrogen bond in the diazine-water complexes. The NBO analysis indicates that intermolecular charge transfer are 0.0316, 0.0255 and 0.0265 e respectively. In addition, the first singlet （n,n＊） vertical excitation energy of the monomer and the hydrogen bonding complexes between diazines and water was investigated by time-dependent density functional theory.     

Competitive interaction between halogen and hydrogen bonds in NH2Br‐HOX (X = F,Cl, and Br) complex

The NH₂Br‐HOX (X = F, Cl, and Br) complexes have been investigated with quantum chemical calculations at the MP2/aug‐cc‐pVTZ level. Five isomers are observed for the Cl and Br complexes, whereas only two isomers are found for the F complex. The geometrical, energetic, and spectroscopic parameters have been analyzed for these complexes. The hydrogen‐bonded complexes are more stable than the halogen‐bonded ones. In most complexes, the associated O? H and O? X bonds are elongated and show a red shift, whereas the distant bonds are contracted and exhibit a blue shift. The complexes have been analyzed with natural bond orbital and atoms in molecules. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

A computational study on the ring stretching modes of halogen‐substituted pyridine involved in H‐bonding

Density functional method B3LYP plus the AUG‐cc‐pVDZ and AUG‐cc‐pVTZ basis sets is used to investigate ring normal modes of halogen‐substituted pyridines involved in the N ··· H? X H‐bonds with HX (X = F, Cl). The results demonstrated that the formation of hydrogen bond leads to an increase in the frequencies of the ring breathing mode v₁, the N‐para‐C stretching mode v_6a and the meta‐CC stretching mode v_8a, whereas there is no change in the triangle mode v₁₂ for free pyridine and a smaller blue shift for substituted pyridines. There is a strong coupling between the C? Y stretching vibration and the triangle mode (ortho‐ and para‐substituted) or the breathing mode (meta‐substituted) in substituted pyridines, which leads to the frequency decrease in the triangle or breathing modes. The natural bond orbital analysis suggests that electrostatic interaction and charge transfer caused by the intermolecular and intramolecular hyperconjugations are the origin of the frequency blue shift in the ring stretching modes. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Hydrogen bonding characterization of XH2NH2···HNO(X = B,Al, Ga) complexes: A theoretical investigation

The complexes of XH₂NH₂···HNO(X = B, Al, Ga) are characterized as head to tail with hydrogen bonding interactions. The structural characteristics can be confirmed by atoms in molecules (AIM) analysis, which also provide comparisons of hydrogen bonds strengths. The calculated interaction energies at G2MP2 level show that stability of complexes decrease as BH₂NH₂···HNO > AlH₂NH₂···HNO > GaH₂NH₂···HNO. On the basis of the vibrational frequencies calculations, there are red‐shifts for ν(X1? H) and blue‐shifts for ν(N? H) in the complexes on dihydrogen bonding formations (X1? H···H? N). On hydrogen bonding formations (N? H···O), there are red‐shifts for ν(N? H) compared to the monomers. Natural bond orbital (NBO) analysis is used to discuss the reasons for the ν(X1? H) and ν(N? H) stretching vibrational shifts by hyperconjugation, electron density redistribution, and rehybridization. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Polarity and structure of P(X)-modified (X = O,S) arylcarbamoylmethylphosphine oxides and sulfides

As shown by the dipole moment method and quantum chemical calculations, (arylcarbamoylmethyl) diphenylphosphine oxides and sulfides exist as equilibrium mixtures of several rotational isomers stabilized by H · · · X intramolecular hydrogen bonds (X = O, S). The most energetically favorable rotamer and its fraction have been determined for each compound.     

1H NMR analysis of O‐methyl‐inositol isomers: a joint experimental and theoretical study

Density functional theory (DFT) calculations of ¹H NMR chemical shifts for l ‐quebrachitol isomers were performed using the B3LYP functional employing the 6‐31G(d,p) and 6‐311 + G(2d,p) basis sets. The effect of the solvent on the B3LYP‐calculated NMR spectrum was accounted for using the polarizable continuum model. Comparison is made with experimental ¹H NMR spectroscopic data, which shed light on the average uncertainty present in DFT calculations of chemical shifts and showed that the best match between experimental and theoretical B3LYP ¹H NMR profiles is a good strategy to assign the molecular structure present in the sample handled in the experimental measurements. Among four plausible O‐methyl‐inositol isomers, the l ‐quebrachitol 2a structure was unambiguously assigned based only on the comparative analysis of experimental and theoretical ¹H NMR chemical shift data. The B3LYP infrared (IR) spectrum was also calculated for the four isomers and compared with the experimental data, with analysis of the theoretical IR profiles corroborating assignment of the 2a structure. Therefore, it is confirmed in this study that a combined experimental/DFT spectroscopic investigation is a powerful tool in structural/conformational analysis studies. Copyright © 2012 John Wiley & Sons, Ltd.     

Multiple‐Well,multiple‐path unimolecular reaction systems. II. 2‐methylhexyl free radicals

Vibrationally excited 2‐methylhexyl radicals formed by shock wave activation or by chemical activation can isomerize by multiple pathways to form any of six stable isomers, can fragment by multiple C? H and C? C bond fission pathways, and can be collisionally stabilized. Master equation simulations of chemical activation and of shock wave activation are used to explore the generic behavior of this complicated coupled system. Selecting the argon pressure in chemical activation systems that produce the 2‐methyl‐1‐hexyl radical isomer ( 1 ) can control the yield of specific isomers. Shock heating of 1 also shows a highly regular sequence of isomer formation. This regular behavior is because the first isomerization steps are faster than subsequent steps. Other radical isomers, such as 2‐methyl‐3‐hexyl ( 3 ), do not show such regular behavior, because the first isomerization step is slower than subsequent steps. Incubation and unimolecular rate‐constant fall‐off are observed in the shock wave simulations. The unimolecular rate‐constant fall‐off for the coupled system produces low‐pressure limiting rate constants proportional to [M]ⁿ, where n can be greater than unity. The fact that n can be greater than unity is a natural feature of multichannel coupled unimolecular reaction systems, but detection of the effect in experiments may be very demanding. © 2001 John Wiley & Sons, Inc. Int J Chem Kinet 33: 246–261, 2001