Effects of trimethylaluminium and tetrakis(ethylmethylamino) hafnium in the early stages of the atomic-layer-deposition of aluminum oxide and hafnium oxide on hydroxylated GaN nanoclusters

We calculate the interactions of two atomic layer deposition (ALD) reactants, trimethylaluminium (TMA) and tetrakis(ethylmethylamino) hafnium (TEMAH) with the hydroxylated Ga-face of GaN clusters when aluminum oxide and hafnium oxide, respectively, are being deposited. The GaN clusters are suitable as testbeds for the actual Ga-face on practical GaN nanocrystals of importance not only in electronics but for several other applications in nanotechnology. We find that TMA spontaneously interacts with hydroxylated GaN; however it does not follow the atomic layer deposition reaction path unless there is an excess in potential energy introduced in the clusters at the beginning of the optimization, for instance, using larger bond lengths of various bonds in the initial structures. TEMAH also does not interact with hydroxylated GaN, unless there is an excess in potential energy. The formation of a Ga—N(CH₃)(CH₂CH₃) bond during the ALD of HfO₂ using TEMAH as the reactant without breaking the Hf—N bond could be the key part of the mechanism behind the formation of an interface layer at the HfO₂/GaN interface.

Penta- and heteropentadienyl ligands coordinated to beryllium

In this work we have performed a systematic study of new organometallic complexes containing penta- and heteropentadienyl (CH₂CHCHCHX, X?=?CH₂, O, NH, S) ligands coordinated to beryllium. Calculated complexes were studied using the density functional theory (PBE) in combination with the 6-311++G(3d,2p) basis set. The coordination number on the beryllium atom varies according to the type of ligand. Pentadienyl ligand shows hapticities η¹ and η⁵, while heteropentadienyl ligands display η¹ and η² hapticities. A Wiberg bond indices study was performed in order to get information about their bond orders.

Computational investigation on microsolvation of the osmolyte glycine betaine [GB (H2O)1-7]

The preferential interactions of glycine betaine (GB) with solvent components and the effect of solvent on its stability have been examined. In particular, the microsolvation of organic osmolyte and widely important osmoprotectant in nature as glycine betaine has been reported by using M06 method. A number of configurations (b_X (a-z)) of the clusters for one to seven water molecules (×?=?1-7) have been considered for the microsolvation. Structures of stable conformers are obtained and denoted as b1a, b2a, b3a, b4a, b5a, b6a and b7a. It is observed from the interaction energy difference (?E) that only seven water molecules can be accommodated in the first solvation shell to stabilize GB. It is also observed that the calculated relative energy using M06 is in close agreement with calculations at the MP2 level of theory.

Theoretical study on cooperative effects between X⋯N and X⋯Carbene halogen bonds (X = F,Cl,Br and I)

Quantum chemical calculations are performed to study the interplay between halogen?nitrogen and halogen?carbene interactions in NCX?NCX?CH₂ complexes, where X?=?F, Cl, Br and I. Molecular geometries and interaction energies of dyads and triads are investigated at the MP2/aug-cc-pVTZ level of theory. It is found that the X?N and X?C_carbene interaction energies in the triads are larger than those in the dyads, indicating that both the halogen bonding interactions are enhanced. The estimated values of cooperative energy E _coop are all negative with much larger E _coop in absolute value for the systems including iodine. The nature of halogen bond interactions of the complexes is analyzed using parameters derived from the quantum theory atoms in molecules methodology and energy decomposition analysis.

Increase of fungitoxicity of mercuric chloride by methionine,ethionine and S-methylcysteine

The fungitoxicity of mercuric chloride to Aspergillus niger was increased in the presence of d-, l-, dl-methionine, dl-ethionine, dl-S-methylcysteine or sodium methylmercaptide. The same effect was observed with methionine for two other fungi investigated: Cladosporium cucumerinum and Scopulariopsis brevicaulis. It is suggested that this effect can be ascribed to the formation of CH₃SHg⁺ or (CH₃S)₂Hg, or the corresponding ethyl compounds. CH₃SHgCl and (CH₃S)₂Hg were synthetically prepared and proved indeed far more fungitoxic than HgCl₂. The hypothesis was further substantiated by the observation that A. niger rapidly converts dl-methionine into CH₃SH, which undoubtedly reacts with Hg²⁺ to give the above mentioned methylthiomercury compounds.     

An assessment of DFT methods for predicting the thermochemistry of ion-molecule reactions of group 14 elements (Si, Ge, Sn)

Experimental mass-spectrometry data on thermochemistry of methide transfer reactions (CH₃)₃M⁺ + M'(CH₃)₄ ? M(CH₃)₄?+?(CH₃)₃M'⁺ (M, M'?=?Si, Ge or Sn) and the formation energy of the [(CH₃)₃Si-CH₃-Si(CH₃)₃]⁺ complex are used as benchmarks for DFT methods (B3LYP, BMK, M06L, and ωB97XD). G2 and G3 theory methods are also used for the prediction of thermochemical data. BMK, M06L, and ωB97XD methods give the best fit to experimental data (close to chemical accuracy) as well as to G2 and G3 results, while B3LYP demonstrates poor performance. From the first three methods M06L gives the best overall result. Structures and formation energies of intermediate “mixed” [(CH₃)₃M-CH₃- M′(CH₃)₃] complexes not observed in experiment are predicted. Their structures, better described as M(CH₃)₄?[M′(CH₃)₃]⁺ complexes, explain their fast decompositions.

DFT study on the reactions of ClO–/BrO– with RCl (R = CH3, C2H5, and C3H7) in gas phase

Gas-phase reactions of ClO^–/BrO^– with RCl (R = CH₃, C₂H₅, and C₃H₇) have been investigated in detail using the popular DFT functional BHandHLYP/aug-cc-pVDZ level of theory. As a result, our findings strongly suggest that the type of reaction is firstly initiated by a typical S_N2 fashion. Subsequently, two competitive substitution steps, named as S_N2-induced substitution and S_N2-induced elimination, respectively, would proceed before the initial S_N2 product ion-dipole complex separates, in which the former exhibits less reactivity than the latter. Those are consistent with relevant experimental results. Moreover, we have also explored reactivity difference for the title reactions in term of some factors derived from methyl group, p-π electronic conjugation, ionization energy (IE), as well as molecular orbital (MO) analysis.

Theoretical investigation into optical and electronic properties of 1,8-naphthalimide derivatives

A series of 1,8-naphthalimide derivatives has been designed to explore their optical, electronic, and charge transport properties as charge transport and/or luminescent materials for organic light-emitting diodes (OLEDs). The frontier molecular orbitals (FMOs) analysis have shown that the vertical electronic transitions of absorption and emission are characterized as intramolecular charge transfer (ICT) for electron-donating and aromatic groups substituted derivatives. However, the ICT character of the electron-withdrawing substituted derivatives is not significant. The calculated results show that their optical and electronic properties are affected by the substituent groups in 4-position of 1,8-naphthalimide. Our results suggest that 1,8-naphthalimide derivatives with electron-donating ?OCH₃ and ?N(CH₃)₂ (1 and 2), electron-withdrawing ?CN and?COCH₃ (3 and 4), 2-(thiophen-2-yl)thiophene (5), 2,3-dihydrothieno[3,4-b][1, 4]dioxine (6), 2-phenyl-1,3,4-oxadiazole (7), and benzo[c][1,2,5]thiadiazole (8) fragments are expected to be promising candidates for luminescent materials for OLEDs, particularly for 5 and 7. In addition, 3 and 7 can be used as promising hole transport materials for OLEDs. This study should be helpful in further theoretical investigations on such kind of systems and also to the experimental study for charge transport and/or luminescent materials for OLEDs.

A new interaction mechanism of LiNH2 with MgH2: magnesium bond

Quantum chemical calculations were performed for LiNH₂–HMgX (X?=?H, F, Cl, Br, CH₃, OH, and NH₂) complexes to propose a new interaction mechanism between them. This theoretical survey showed that the complexes are stabilized through the combinative interaction of magnesium and lithium bonds. The binding energies are in the range of 63.2–66.5 kcal mol^?1, i.e., much larger than that of the lithium bond. Upon complexation, both Mg–H and Li–N bonds are lengthened. Substituents increase Mg-H bond elongation and at the same time decrease Li-N bond elongation. These cyclic complexes were characterized with the presence of a ring critical point and natural population analysis charges.

Effect of surface hydroxyls on dimethyl ether synthesis over the γ-Al2O3 in liquid paraffin: a computational study

In a recent paper (Zuo et al., Appl Catal A 408:130–136, 2011), the mechanism of dimethyl ether (DME) synthesis from methanol dehydration over γ-Al₂O₃ (110) was studied using density functional theory (DFT). Using the same method, the effect of surface hydroxyls on γ-Al₂O₃ in liquid paraffin during DME synthesis from methanol dehydration is investigated. It is found that DME is mainly formed from two adsorbed CH₃O groups via methanol dehydrogenation on both dehydrated and hydrated γ-Al₂O₃ in liquid paraffin. No close correlation between catalytic activity and acid intensity was found. Before and after water adsorption at typical catalytic conditions (e.g., 553 K), the reaction rate is not obviously changed on γ-Al₂O₃(100) surface in liquid paraffin, but the reaction rate decreases by about 11 times on the (110) in liquid paraffin. Considering the area of the (110) and (100) surfaces under actual conditions, the catalytic activity decreased mainly because the Al3 sites on the (110) surface gradually become inactive. Catalytic activity decreased mainly due to surface hydrophilicity. The calculated results were consistent with the experiment.

Theoretical investigations on the synthesis mechanism of cyanuric acid from NH3 and CO2

In the synthesis of cyanuric acid from NH₃ and CO₂, urea and isocyanic acid OCNH are two pivotal intermediates. Based on density functional theory (DFT) calculations, the synthesis mechanism of cyanuric acid from NH₃ + CO₂ was investigated systematically. Urea can be synthesized from NH₃ and CO₂, and cyanuric acid can be obtained from urea or NH₃ + CO₂. In the stepwise mechanism of cyanuric acid from urea or NH₃ + CO₂, the energy barriers are relatively high, and the condition of high pressure and temperature does not decrease the energy barriers. Our theoretical model shows that cyanuric acid is actually acquired from OCNH via a one-step cycloaddition reaction.

Insights into the strength and nature of carbene···halogen bond interactions: a theoretical perspective

Halogen-bonding, a noncovalent interaction between a halogen atom X in one molecule and a negative site in another, plays critical roles in fields as diverse as molecular biology, drug design and material engineering. In this work, we have examined the strength and origin of halogen bonds between carbene CH₂ and XCCY molecules, where X?=?Cl, Br, I, and Y?=?H, F, COF, COOH, CF₃, NO₂, CN, NH₂, CH₃, OH. These calculations have been carried out using M06-2X, MP2 and CCSD(T) methods, through analyses of surface electrostatic potentials V _S(r) and intermolecular interaction energies. Not surprisingly, the strength of the halogen bonds in the CH₂···XCCY complexes depend on the polarizability of the halogen X and the electron-withdrawing power of the Y group. It is revealed that for a given carbene···X interaction, the electrostatic term is slightly larger (i.e., more negative) than the dispersion term. Comparing the data for the chlorine, bromine and iodine substituted CH₂···XCCY systems, it can be seen that both the polarization and dispersion components of the interaction energy increase with increasing halogen size. One can see that increasing the size and positive nature of a halogen’s σ-hole markedly enhances the electrostatic contribution of the halogen-bonding interaction.

Spin-flip reactions of Zr + C2H6 researched by relativistic density functional theory

Density functional theory (DFT) with relativistic corrections of zero-order regular approximation (ZORA) has been applied to explore the reaction mechanisms of ethane dehydrogenation by Zr atom with triplet and singlet spin-states. Among the complicated minimum energy reaction path, the available states involves three transition states (TS), and four stationary states (1) to (4) and one intersystem crossing with spin-flip (marked by ?): ³ Zr + C ₂ H ₆ → ³ Zr-CH ₃ -CH ₃ ( ³ 1) → ³ TS _1/2 → ³ ZrH-CH ₂ -CH ₃ ( ³ 2) → ³ TS _2/3 ? ¹ ZrH₂-CH₂ = CH₂ ( ¹ 3) → ¹ TS _3/4 → ¹ ZrH ₃ -CH = CH ₂ ( ¹ 4). The minimum energy crossing point is determined with the help of the DFT fractional-occupation-number (FON) approach. The spin inversion leads the reaction pathway transferring from the triplet potential energy surface (PES) to the singlet’s accompanying with the activation of the second C-H bond. The overall reaction is calculated to be exothermic by about 231 kJ mol^?1. Frequency and NBO analysis are also applied to confirm with the experimental observed data.

Theoretical study of the triangular bonding complex formed by carbon tetrabromide, a halide, and a solvent molecule in the gas phase

MP2(full)/aug-cc-pVDZ(-PP) computations predict that new triangular bonding complexes (where X^? is a halide and H–C refers to a protic solvent molecule) consist of one halogen bond and two hydrogen bonds in the gas phase. Carbon tetrabromide acts as the donor in the halogen bond, while it acts as an acceptor in the hydrogen bond. The halide (which commonly acts as an acceptor) can interact with both carbon tetrabromide and solvent molecule (CH₃CN, CH₂Cl₂, CHCl₃) to form a halogen bond and a hydrogen bond, respectively. The strength of the halogen bond obeys the order CBr₄???Cl^? > CBr₄???Br^? > CBr₄???I^?. For the hydrogen bonds formed between various halides and the same solvent molecule, the strength of the hydrogen bond obeys the order C-H???Cl^? > C-H???Br^? > C-H???I^?. For the hydrogen bonds formed between the same halide and various solvent molecules, the interaction strength is proportional to the acidity of the hydrogen in the solvent molecule. The diminutive effect is present between the hydrogen bonds and the halogen bond in chlorine and bromine triangular bonding complexes. Complexes containing iodide ion show weak cooperative effects.

Computationally designed prodrugs of statins based on Kirby’s enzyme model

DFT calculations at B3LYP/6-31G(d,p) for intramolecular proton transfer in Kirby’s enzyme models 1–7 demonstrated that the reaction rate is dependent on the distance between the two reacting centers, r_GM, and the hydrogen bonding angle, α, and the rate of the reaction is linearly correlated with r_GM and α. Based on these calculation results three simvastatin prodrugs were designed with the potential to provide simvastatin with higher bioavailability. For example, based on the calculated log EM for the three proposed prodrugs, the interconversion of simvastatin prodrug ProD 3 to simvastatin is predicted to be about 10 times faster than that of either simvastatin prodrug ProD 1 or simvastatin ProD 2. Hence, the rate by which the prodrug releases the statin drug can be determined according to the structural features of the promoiety (Kirby’s enzyme model).

Revealing substitution effects on the strength and nature of halogen-hydride interactions: a theoretical study

A quantum chemistry study was carried out to investigate the strength and nature of halogen bond interactions in HXeH···XCCY complexes, where X = Cl, Br and Y = H, F, Cl, Br, CN, NC, C₂H, CH₃, OH, SH, NH₂. Examination of the electrostatic potentials V(r) of the XCCY molecules reveals that the addition of substituents has a significant effect upon the most positive electrostatic potential on the surface of the interacting halogen atom. We found that the magnitude of atomic charges and multipole moments depends upon the halogen atom X and is rather sensitive to the electron-withdrawing/donating power of the remainder of the molecule. An excellent correlation was found between the most positive electrostatic potentials on the halogen atom and the interaction energies. For either HXeH···ClCCY or HXeH···BrCCY complexes, an approximate linear correlation between the interaction energies and halogens multipole moments are established, indicating that the electrostatic and polarization interactions are responsible for the stability of the complexes. According to energy decomposition analysis, it is revealed that the electrostatic interactions are the major source of the attraction in the HXeH···XCCY complexes. Furthermore, the changes in the electrostatic term are mainly responsible for the dependence of interaction energy on the halogen atom.

Theoretical study about the 5-azido-1H-tetrazole and its ion salts

Periodic DFT method has been firstly used to calculate the bulk structure, electronic structure, electrical transferring and thermodynamic properties of crystalline 5-azido-1H-tetrazole (HCN₇) and its four different salts. The anion CN₇ ^? was included in all of the salts such as ammonium 5-azidotetrazolate ([NH₄]⁺[CN₇]^?), hydrazinium 5-azidotetrazolate ([N₂H₅]⁺[CN₇]^?), guanidinium 5-azidotetrazolate ([CH₆N₃]⁺[CN₇]^??·?H₂O) and 1-aminoguanidinium 5- azidotetrazolate ([CH₇N₄]⁺[CN₇]^?). The simulation is in reasonable agreement with the experimental results. It is found the salts of HCN₇ are more stable than itself because the band gap of the salts is larger. The density of state shows the p states of them (including HCN₇ and its four salts) have played a very significant role in the reaction.

The average local ionization energy as a tool for identifying reactive sites on defect-containing model graphene systems

In a continuing effort to further explore the use of the average local ionization energy $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ as a computational tool, we have investigated how well $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ computed on molecular surfaces serves as a predictive tool for identifying the sites of the more reactive electrons in several nonplanar defect-containing model graphene systems, each containing one or more pentagons. They include corannulene (C₂₀H₁₀), two inverse Stone-Thrower-Wales defect-containing structures C₂₆H₁₂ and C₄₂H₁₆, and a nanotube cap model C₂₂H₆, whose end is formed by three fused pentagons. Coronene (C₂₄H₁₂) has been included as a reference planar defect-free graphene model. We have optimized the structures of these systems as well as several monohydrogenated derivatives at the B3PW91/6-31G* level, and have computed their $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ on molecular surfaces corresponding to the 0.001 au, 0.003 au and 0.005 au contours of the electronic density. We find that (1) the convex sides of the interior carbons of the nonplanar models are more reactive than the concave sides, and (2) the magnitudes of the lowest $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ surface minima (the $ {{\overline{\mathrm{I}}}_{{\mathrm{S}\text{,}\min }}} $ ) correlate well with the interaction energies for hydrogenation at these sites. These $ {{\overline{\mathrm{I}}}_{{\mathrm{S}\text{,}\min }}} $ values decrease in magnitude as the nonplanarity of the site increases, consistent with earlier studies. A practical benefit of the use of $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ is that a single calculation suffices to characterize the numerous sites on a large molecular system, such as graphene and defect-containing graphene models.