Computational study on the donor-adducts of phosphinidene complexes

The bonding of adducts between electron donors (D), such as NH₃ and PH₃, with phosphinidenes RP: (R = H, CH₃, NH₂, PH₂, PF₂) was examined using density functional theory calculations. The structures were optimized at the B3LYP/6-311++G⁷ level of theory. The main interest of this work is to find the appropriate donor that can stabilize the phosphinidene. Accordingly, the magnitude of bond strength correlates well with the difference in ionization potential (ΔIP) of the orbitals involved in the D?P interaction.     

Computational study of the thermodynamic stabilities of hydrogen-bonded complexes in solution

Understanding of the multiple H-bonding arrays of heterocyclic compounds is essential to design effective building blocks of supramolecular polymers. We have carried out a comprehensive computational study on the thermodynamic stabilities of thirty-six H-bonded complexes with all possible H-bonding arrays in the gas phase and chloroform solvent by using M06-2X, SMD calculations and cc-pVDZ basis set. The multiple H-bonding arrays include donor acceptor–acceptor donor (DA–AD), DD–AA for the doubly H-bonded pairs, and DAD–ADA, DDA–AAD and DDD–AAA for the triply H-bonded pairs. The computational results have provided insights into the geometrical, energetic and solvation effects on the stabilities of these H-bonded complexes. The calculated free energies of association for the DD–AA (8–9) and the DDD–AAA (33–35, 36–35) H-bonded complexes are found to be inconsistent with the experimental measurements and observations that these complexes are the most strongly doubly and triply H-bonded pairs in solution, respectively, while the calculated binding free energies for all other H-bonding arrays are in good agreement with experimental values. The computational protocol can be used by practical chemists and undergraduate researchers as an efficient and state-of-the-art tool to study H-bonding interactions in supramolecular chemistry.     

Computational study of substituent effect in para substituted platinabenzene complexes

The electronic structure and properties of the platinabenzene and para substituted platinabenzenes have been investigated using the hybrid density functional mpw1pw91 theory. The substituent effect in structure parameters, frontier orbital energies, aromaticity indexes, and hyperpolarizability has been studied. The calculations show that, in all molecules HOMO → LUMO transition makes the major contribution in the most intense electronic transition.     

Computational study of noncovalent complexes between formamide and formic acid

The geometries and binding energies of 1:1, 1:2, and 1:4 formic acid-formamide complexes (FA-FMA) are calculated by quantum chemical procedures. Vibrational spectra and intermolecular distances of the most stable FA-FMA dimers as well as the influence of the basis set superposition error (BSSE) on the geometries and energies of the dimers are also discussed. All FA-FMA dimers are optimized at the B3LYP/cc-pVTZ, the MP2/cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ levels of theory to study the influence of the level of theory on the calculated geometries and energies. CCSD(T)/cc-pVTZ single-point calculations at the MP2/aug-cc-pVTZ-optimized geometries were performed as reference for estimating the quality of lower level calculations. These calculations allow us to qualitatively describe the competition between different types of hydrogen-bonding interactions in FA-FMA complexes. FA-FMA dimers are compared to other formamide complexes and to the FA-FMA crystal structure.     

HFI and DFT study of the bonding in complexes of halogen and interhalogen diatomics with Lewis base

We have analyzed by means of DFT calculations with use of the pseudo-potential the nuclear quadrupole coupling constants of a range of XYB complexes (n aσ type in Mulliken notation) formed between diatomic interhalogen molecules XY and Lewis bases B. The geometrical parameters, rotational and halogen nuclear quadrupole coupling constants obtained by these calculations substantially corresponded to the data of microwave spectroscopy in the gas phase. An analysis of the quality of the calculations that employ the pseudo-potential and the expanded basis set for the halogen compounds was carried out. The ZORA model is shown to be a viable alternative to the computationally demanding BH and HLYP model for the calculation of halogen and nitrogen coupling constants in molecules. In addition, the ZORA model, in contrast to the pseudo-potential model, leads to realistic values of iodine nuclear quadrupole coupling constants. From electron partitioning analyses and Klopman's approach it follows that for the IClB complexes the electrostatic bonding is predominant relative to covalent bonding.     

Computational study of formamide–water complexes using the SAPT and AIM methods

In this work, the complexes formed between formamide and water were studied by means of the SAPT and AIM methods. Complexation leads to significant alterations in the geometries and electronic structure of formamide. Intermolecular interactions in the complexes are intense, especially in the cases where the solvent interacts with the carbonyl and amide groups simultaneously. In the transition states, the interaction between the water molecule and the lone pair on the amide nitrogen is also important. In all the complexes studied herein, the electrostatic interactions between formamide and water are the main attractive force, and their contribution may be five times as large as the corresponding contribution from dispersion, and twice as large as the contribution from induction. However, an increase in the resonance of planar formamide with the successive addition of water molecules may suggest that the hydrogen bonds taking place between formamide and water have some covalent character.     

Computational study of methane activation by mercury(II) complexes

A computational investigation of methane activation by Hg^II complexes is reported. Calculated geometries and energetics of Hg(II)-containing reactants and products are consistent with available experimental data for a wide range of diverse ligand types. Calculated reaction enthalpies and activation barriers for Hg^II complexes cover a wide range of values for different ligands. This diversity suggests that the kinetics and thermodynamics of methane activation by Hg(II) and related medium-valent complexes can be tailored through rational modification of the ligand environment. Calculated activation barriers and reaction enthalpies for methane activation by Hg(II) complexes indicates that harder, more electronegative, ligands are kinetically and thermodynamically preferred. Potential donor groups on the activating ligand can stabilize the transition state versus the ground state reactants and hence result in substantially lower methane activation barriers. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 902–911, 1998     

Computational studies on stable triplet states of heteroacetylenes and the effects of halogen substituents

This paper describes theoretical studies of halogen-substituted heteroacetylenes (XCMY, M = Si and Ge; X, Y = H, Cl and F) performed at the QCISD(T)/6-311G//QCISD/6-31G level of theory. The electronegative halogen substituents destabilize the singlet state such that the triplet state tends to become favorable. The triplet state has the bifunctional electronic structure of a triplet carbene joined to a heavy singlet carbene. We found that the substituents effectively reduce the energy of the donor-acceptor interactions (E(D-A)) between the two in-plane lone pairs of electrons of the singlet state; therefore, the remaining pi bond is less favorable energetically than the triplet state with a sigma bond. A related phenomenon occurs for the homonuclear heavy acetylenes in singlets in which the lead compound RPbPbR switches to a Pb-Pb sigma bond from the pi bonds observed for the lighter acetylenes.     

Energetics of halogen bonding of group 10 metal fluoride complexes

A study is presented of the thermodynamics of the halogen-bonding interaction of C(6)F(5)I with a series of structurally similar group 10 metal fluoride complexes trans-[Ni(F)(2-C(5)NF(4))(PCy(3))(2)] (2), trans-[Pd(F)(4-C(5)NF(4))(PCy(3))(2)] (3), trans-[Pt(F){2-C(5)NF(2)H(CF(3))}(PR(3))(2)] (4a, R = Cy; 4bR = iPr) and trans-[Ni(F){2-C(5)NF(2)H(CF(3))}(PCy(3))(2)] (5a) in toluene solution. (19)F NMR titration experiments are used to determine binding constants, enthalpies and entropies of these interactions (2.4 ≤ K(300) ≤ 5.2; -25 ≤ ΔH(o) ≤ -16 kJ mol(-1); -73 ≤ ΔS(o) ≤ -49 J K(-1) mol(-1)). The data for -ΔH(o) for the halogen bonding follow a trend Ni < Pd < Pt. The fluoropyridyl ligand is shown to have a negligible influence on the thermodynamic data, but the influence of the phosphine ligand is significant. We also show that the value of the spin-spin coupling constant J(PtF) increases substantially with adduct formation. X-ray crystallographic data for Ni complexes 5a and 5c are compared to previously published data for a platinum analogue. We show by experiment and computation that the difference between Pt-X and Ni-X (X = F, C, P) bond lengths is greatest for X = F, consistent with F(2pπ)-Pt(5dπ) repulsive interactions. DFT calculations on the metal fluoride complexes show the very negative electrostatic potential around the fluoride. Calculations of the enthalpy of adduct formation show energies of -18.8 and -22.8 kJ mol(-1) for Ni and Pt complexes of types 5 and 4, respectively, in excellent agreement with experiment.     

Computational study of pH‐responsive di‐lanthanide complexes

Lanthanide complexes have found extensive use as luminescent probes for biological and medical investigations. Recently, a di‐europium complex that exhibits pH‐dependent luminescence‐decay was reported, and the ligand in that complex includes a large number of ionizable sites. To better understand the pH‐dependence of luminescence‐decay of this complex, the pK _a's of all tautomers of the di‐Lu³⁺ version of this complex were calculated computationally. The calculated Boltzmann‐averaged pK _a's of the complex are 5.85, ?0.21, and ?1.47 for the di‐Lu³⁺ complex in its first, second, and third protonation states, respectively. These pK _a values across protonation states indicate that changes in luminescence‐decay rate at physiologically relevant pH may be related to first protonation event of the complex exclusively.     

Computational study of alkynes insertion into metal-hydride bonds catalyzed by bimetallic complexes

Density Functional Theory investigations on the insertion mechanism of phenylacetylene into metal-hydride bonds in bimetallic (Pt,Os) catalysts have been carried out. The results obtained have been also compared with the non-reactive monometallic (Os-based) system, to elucidate the cooperative effects and to explain the observed absence of reactivity. The identified reaction path involves phenylacetylene coordination followed by the insertion into the metal-hydride bond, leading to the formation of the experimentally observed products. Both steps do not require large energies compatible with the experimental conditions. The comparison with the reaction path for the monometallic species gives some hints on the cooperative effects due to the presence of the second metal which is related to its role in the CO release for creating a coordination site for phenylacetylene and not in the insertion energetics. The calculations provide a detailed analysis of the reaction complexity and provide a rationale for the efficiency of the process.     

Computational study of the lowest triplet state of ruthenium polypyridyl complexes used in artificial photosynthesis

The potential energy surfaces of the first excited triplet state of some ruthenium polypyridyl complexes were investigated by means of density functional theory. Focus was placed on the interaction between the geometrical changes accompanying the photoactivity of these complexes when used as antenna complexes in artificial photosynthesis and dye-sensitized solar cells and the accompanying changes in electronic structure. The loss process (3)MLCT --> (3)MC can be understood by means of ligand-field splitting, traced down to the coordination of the central ruthenium atom.     

Theoretical investigation of the halogen bonded complexes between carbonyl bases and molecular chlorine

The halogen bonded complexes between six carbonyl bases and molecular chlorine are investigated theoretically. The interaction energies calculated at the CCSD(T)/aug‐cc‐pVTZ level range between ?1.61 and ?3.50 kcal mol^?1. These energies are related to the ionization potential, proton affinity, and also to the most negative values (V_s,min) on the electrostatic potential surface of the carbonyl bases. A symmetry adapted perturbation theory decomposition of the energies has been performed. The interaction results in an elongation of the Cl? Cl bond and a contraction of the CF and CH bonds accompanied by a blue shift of the ν(CH) vibrations. The properties of the Cl₂ molecules are discussed as a function of the σ*(Cl? Cl) occupation, the hybridization, and the occupation of the Rydberg orbitals of the two chlorine atoms. Our calculations predict a large enhancement of the infrared and Raman intensities of the ν(Cl? Cl) vibration on going from isolated to complexed Cl₂. © 2015 Wiley Periodicals, Inc.     

Computational study of hydrogen-bonded complexes between the most stable tautomers of L-leucine and bases of RNA

A total of 16 hydrogen-bonded complexes between the lowest energy tautomers of L-leucine and each base of RNA have been characterized at the DFT level of theory. The most stable complexes are formed by L-leucine and guanine. Considering backbone of RNA, the affinity order between the L-leucine and bases is guanine> cytosine> adenine≈uracil in perfect accordance with the experimentation. The interplay between the transformed H–X bonds’ structural parameters from two monomers to the dimer accompanying with the shifts of the frequency for the H–X stretching mode and interaction energies has been discussed.     

A post‐HF study on the halogen bonding interaction of pyrene with diatomic halogen molecules

We present a detailed SCS‐MP2 study on the potential energy curves (PEC) for interactions between diatomic halogen molecules and pyrene. BSSE corrected CCSD[T] energies at equilibrium distances are computed and compared to CCSD(T) energies. The most stable conformation of these weakly bound van der Waals complexes is almost linear in the perpendicular direction to the pyrene plane. The complexes of highly polarizable bromine and iodine molecules with pyrene are very stable and they carry rather large number of vibrational states. Despite its small size, F₂ also forms strong halogen bonding similar to Br₂ and I₂. The interaction between Cl₂ and pyrene is the weakest and it is attributed to the highest polarizability / molar mass ratio of chlorine among the others. I₂‐pyrene is found to be the most stable complex due to the strongest mutual polarization effects and is carrying more than 60 vibrational states. Due to the rather large number of electrons in some complexes, the relativistic corrections are also considered. © 2016 Wiley Periodicals, Inc.     

Substituent effects in halogen bonding complexes between aromatic donors and acceptors: a comprehensive ab initio study

Substituent effects in halogen bonding complexes involving aromatic rings are investigated. We have analyzed how the interaction energy (the RI-MP2/aug-cc-pVDZ level of theory) is affected by the substitution in both halogen bond donor and acceptor aromatic moieties. In addition, we have used two different aromatic electron donor molecules pyridine and cyanobenzene, which allow us to study the effect of having the electron donor nitrogen atom forming part of the ring or outside the ring (-CN). Interestingly, the effect of the substituents on the interaction energies is similar in both cases. We have obtained the Hammett's plots for four combinations of aromatic donors and acceptors and in all cases we have obtained good regression plots (interaction energies vs. Hammett's σ parameter). We have also studied and compared bifurcated halogen bonds using both possible combinations, that is two donors and one acceptor and vice versa. In addition, we have analyzed the effect of the solvent on the interaction energies using COSMO. Finally, we have used Bader's theory of "atoms-in-molecules" to demonstrate that the electron density computed at the bond critical point that emerges upon complexation can be used as a measure of bond order in this noncovalent interaction.     

Platinum complexes of diazo ligands. Studies of regioselective aromatic ring amination, oxidative halogen addition and reductive halogen elimination reactions

2-(Arylazo)pyridine ligands, L1a-1c react with the salt K2[PtCl4] to give the mononuclear complexes [PtCl2(L1)](1), which readily react with ArNH2 to yield the monochloro complexes of type [PtCl(L2)](HL2= 2-[(2-(arylamino)phenyl)azo]pyridine)(2) via regioselective ortho-amine fusion at the pendent aryl ring of coordinated L1. Oxidative addition of the electrophiles Y2(Y = Cl, Br, I) to the square-planar platinum(II) complex, has led to syntheses of the corresponding octahedral platinum(IV) complexes, [PtY3(L2)](3) in high yields. Ascorbate ion reductions of the platinum(IV) complexes, , resulted in reductive halogen elimination to revert to the platinum(II) complexes almost quantitatively. Isolation of products and X-ray structure determination of the representative complexes followed all these chemical reactions. In crystal packing, the compound [PtCl2(L1c)](1c) forms dimeric units with a Pt...Pt distance of 3.699(1) A. In contrast, the crystal packing of 2b revealed that the molecules are arranged in an antiparallel fashion to form a noncovalent 1D chain to accommodate pi(aryl)-pi(pyridyl) and Pt-pi(aryl) interactions. Notably, the oxidation of [Pt(II)Cl(L2a)](2a) by I2 produced a mixed halide complex [Pt(IV)ClI2(L2a)](5), which, in turn, is reduced by ascorbate ion to produce [Pt(II)I(L2a)] with the elimination of ClI. All the platinum(II) complexes are brown, the platinum(IV) complexes, on the other hand, are green. Low-energy visible range transitions in the complexes of the extended ligand [L2]- are ascribed to ligand basedpi-pi* transitions. Cyclic voltammetric behaviour of the complexes is reported.     

FCl:PCX complexes: old and new types of halogen bonds

MP2/aug'-cc-pVTZ calculations have been performed to investigate the halogen-bonded complexes FCl:PCX, for X = NC, CN, F, H, CCH, CCF, CH(3), Li, and Na. Although stable complexes with a F-Cl···P halogen bond exist that form through the lone pair at P (configuration I), except for FCl:PCCN, the more stable complexes are those in which FCl interacts with the C≡P triple bond through a perturbed π system (configuration II). In complexes I, the nature of the halogen bond changes from traditional to chlorine-shared and the interaction energies increase, as the electron-donating ability of X increases. The anionic complex FCl:PC(-) has a chlorine-transferred halogen bond. SAPT analyses indicate that configuration I complexes with traditional halogen bonds are stabilized primarily by the dispersion interaction. The electrostatic interaction is the most important for configuration I complexes with chlorine-shared halogen bonds and for configuration II complexes except for FCl:PCNa for which the induction term is most important. The F-Cl stretching frequency is red-shifted upon complexation. EOM-CCSD/(qzp,qz2p) spin-spin coupling constants have been obtained for all FCl:PCX complexes with configuration I. (1)J(F-Cl) decreases upon complexation. (2X)J(F-P) values are quadratically dependent upon the F-P distance and are very sensitive to halogen-bond type. (1X)J(Cl-P) tends to increase as the Cl-P distance decreases but then decreases dramatically in the chlorine-transferred complex FCl:PC(-) as the Cl-P interaction approaches that of a covalent Cl-P bond. Values of (1)J(F-Cl) for configuration II are reduced relative to configuration I, reflecting the longer F-Cl distances in II compared to those of the neutral complexes of I. Although the F-P and Cl-P distances in configuration II complexes are shorter than these distances in the corresponding configuration I complexes, (2X)J(F-P) and (1X)J(Cl-P) values are significantly reduced, indicating that coupling through the perturbed C-P π bond is less efficient. The nature of F-P coupling for configuration II is also significantly different, as evidenced by the relative importance of PSO, FC, and SD components.     

Benchmarking of model core potentials: application to the halogen complexes of group 4 metals

The reliability of the model core potential (MCP) method was probed in a systematic RHF and MP2 study of the geometries of the group 4 metal halogen complexes (MX4; M = Ti, Zr, Hf and X = F, Cl, Br, I). The computed bond lengths were compared with experimental values, as well as those predicted using effective core potentials. Provided that electrons from the outermost core shell of the metal atom are treated explicitly in the calculation, both the MCP and ECP methods predict M-X bond lengths within 0.02-0.03 A of experiment. The reaction energies for a simple set of halogen substitution reactions of the MX4 complexes leading to the mixed halogen complex, MX2Y2, were also studied. Although no experimental values are available for these reactions, comparison was made with the values computed using effective core potentials. The predictability of the different pseudopotential techniques and the importance of the metal atom valence basis set contraction scheme and polarization space are discussed.