Architectures,electronic structures,and stabilities of Cu-doped Ge<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript> clusters: density functional modeling

The present study reports the geometries, electronic structures, growth behavior, and stabilities of neutral and ionized copper-doped germanium clusters containing 1–20 Ge atoms within the framework of linear combination of atomic orbitals density functional theory (DFT) under the spin-polarized generalized gradient approximation. It was found that Cu-capped Ge _n (or Cu-substituted Ge _n+1) and Cu-encapsulated Ge _n clusters mostly occur in the ground state at a particular cluster size (n). In order to explain the relative stabilities of the ground-state clusters, parameters such as the average binding energy per atom (BE), the embedding energy (EE), and the fragmentation energy (FE) of the clusters were calculated, and the resulting values are discussed. To explain the chemical stabilities of the clusters, parameters such as the energy gap between the highest occupied and the lowest unoccupied molecular orbitals (the HOMO–LUMO gap), the ionization energy (IP), the electron affinity (EA), the chemical potential (μ), the chemical hardness (η), and the polarizability were calculated, and the resulting values are also discussed. Natural atomic orbital (NAO) and natural bond orbital (NBO) analyses were also used to determine the electron-counting rule that should be applied to the most stable Ge₁₀Cu cluster. Finally, the relevance of the calculated results to the design of Ge-based superatoms is discussed.

Open image in new window Open image in new window

Figure Contributions of the valance orbitals of the Ge and Cu atom(s) to the HOMO of the ground-state icosahedral Ge₁₀Cu cluster obtained from NBO analysis. The numbers below the clusters represent the occupancies of the HOMO orbitals

     

CO2 adsorption on polar surfaces of ZnO

Physical and chemical adsorption of CO₂ on ZnO surfaces were studied by means of two different implementations of periodic density functional theory. Adsorption energies were computed and compared to values in the literature. In particular, it was found that the calculated equilibrium structure and internuclear distances are in agreement with previous work. CO₂ adsorption was analyzed by inspection of the density of states and electron localization function. Valence bands, band gap and final states of adsorbed CO₂ were investigated and the effect of atomic displacements analyzed. The partial density of states (PDOS) of chemical adsorption of CO₂ on the ZnO(0001) surface show that the p orbitals of CO₂ were mixed with the ZnO valence band state appearing at the top of the valence band and in regions of low-energy conduction band.

In pursuit of negative Fukui functions: examples where the highest occupied molecular orbital fails to dominate the chemical reactivity

In our quest to explore molecules with chemically significant regions where the Fukui function is negative, we explored reactions where the frontier orbital that indicates the sites for electrophilic attack is not the highest occupied molecular orbital. The highest occupied molecular orbital (HOMO) controls the location of the regions where the Fukui function is negative, supporting the postulate that negative values of the Fukui function are associated with orbital relaxation effects and nodal surfaces of the frontier orbitals. Significant negative values for the condensed Fukui function, however, were not observed.

Non-covalent interactions – QTAIM and NBO analysis

MP2(full)/6-311++G(3df,3pd) calculations were carried out on complexes linked through various non-covalent Lewis acid – Lewis base interactions. These are: hydrogen bond, dihydrogen bond, hydride bond and halogen bond. The quantum theory of ´atoms in molecules´ (QTAIM) as well as the natural bond orbitals (NBO) method were applied to analyze properties of these interactions. It was found that for the A-H…B hydrogen bond as well as for the A-X…B halogen bond (X designates halogen) the complex formation leads to the increase of s-character in the A-atom hybrid orbital aimed toward the H or X atom. In opposite, for the A…H-B hydride bond, where the H-atom possesses negative charge, the decrease of s-character in the B-atom orbital is observed. All these changes connected with the redistribution of the electron charge being the effect of the complex formation are in line with Bent´s rule. The numerous correlations between energetic, geometrical, NBO and QTAIM parameters were also found.

Spectroscopic investigations and hydrogen bond interactions of 8-aza analogues of xanthine, theophylline and caffeine: a theoretical study

The structure, spectral properties and the hydrogen bond interactions of 8-aza analogues of xanthine, theophylline and caffeine have been studied by using quantum chemical methods. The time-dependent density functional theory (TD-DFT) and the singly excited configuration interaction (CIS) methods are employed to optimize the excited state geometries of isolated 8-azaxanthine, 8-azatheophylline tautomers and 8-azacaffeine in both the gas and solvent phases. The solvent phase calculations are performed using the polarizable continuum model (PCM). The absorption and emission spectra are calculated using the time-dependent density functional theory (TD-DFT) method. The results from the TD-DFT calculations reveal that the excitation spectra are red shifted relative to absorption in aqueous medium. These changes in the transition energies are qualitatively comparable to the experimental data. The examination of molecular orbital reveals that the molecules with a small H→L energy gap possess maximum absorption and emission wavelength. The relative stability and hydrogen bonded interactions of mono and heptahydrated 8-azaxanthine, 8-azatheophylline tautomers and 8-azacaffeine have been studied using the density functional theory (DFT) and Møller Plesset perturbation theory (MP2) implementing the 6-311++G(d,p) basis set. The formation of strong N-H…O bond has resulted in the highest interaction energy among the monohydrates. Hydration does not show any significant impact on the stability of heptahydrated complexes. The atoms in molecule (AIM) and natural bonding orbital (NBO) analyses have been performed to elucidate the nature of the hydrogen bond interactions in these complexes.

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.

Natures of benzene-water and pyrrole-water interactions in the forms of σ and π types: theoretical studies from clusters to liquid mixture

A combined and sequential use of quantum mechanical (QM) calculations and classical molecular dynamics (MD) simulations was made to investigate the σ and π types of hydrogen bond (HB) in benzene-water and pyrrole-water as clusters and as their liquid mixture, respectively. This paper aims at analyzing similarities and differences of these HBs resulted from QM and MD on an equal footing. Based on the optimized geometry at ωb97xD/aug-cc-pVTZ level of theory, the nature and property of σ and π types of HBs are unveiled by means of atoms in molecules (AIM), natural bond orbital (NBO) and energy decomposition analysis (EDA). In light of the above findings, MD simulation with OPLS-AA and SPC model was applied to study the liquid mixture at different temperatures. The MD results further characterize the behavior and structural properties of σ and π types HBs, which are somewhat different but reasonable for the clusters by QM. Finally, we provide a reasonable explanation for the different solubility between benzene/water and pyrrole/water.

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.

Mechanisms on electrical breakdown strength increment of polyethylene by acetophenone and its analogues addition: a theoretical study

A theoretical investigation is completed on the mechanism of electrical breakdown strength increment of polyethylene. It is shown that it is one of the most important factors for increasing electrical breakdown strength of polyethylene through keto-enol isomerization of acetophenone and its analogues at the ground state S₀ and the lowest triplet state T₁. The minimum structures and transition states of the keto- and the enol-tautomer of acetophenone and its analogues at the S₀ and T₁ states are obtained at the B3LYP/6-311+G(d,p) level, as well as the harmonic vibration frequencies of the equilibrium geometries and the minimum energy path (MEP) by the intrinsic reaction coordinate (IRC) theory at the same level. The two C–C bond cleavage reaction channels have been identified in acetophenone. The calculated results show that the energy barriers of keto-enol isomerization of acetophenone and its analogues at S₀ and T₁ states are much smaller than the average C-C bond energy of polyethylene, and the acetophenone doping or bond linked into polyethylene can increase the electrical breakdown strength and inhibit polyethylene electrical tree initiation and aging.

Density functional studies on photophysical properties and chemical reactivities of the triarylboranes: effect of the constraint of planarity

The geometric and electronic structures, absorption spectra, transporting properties, chemical reactivity indices and electrostatic potentials of the planar three-coordinate organoboron compounds 1-2 and twisted reference compound Mes ₃ B, have been investigated by employing density functional theory (DFT) and conceptual DFT methods to shed light on the planarity effects on the photophysical properties and the chemical reactivity. The results show that the planar compounds 1-2 exhibit significantly lower HOMO level than Mes ₃ B, owing to the stronger electronic induction effect of boron centers. This feature conspicuously induces a blue shifted absorption for 1, although 1 seemingly possesses more extended conjugation framework than Mes ₃ B. Importantly, the reactivity strength of the boron atoms in 1-2 is much lower than that in Mes ₃ B, despite the fact that the tri-coordinate boron centers of 1-2 are completely naked. The interesting and abnormal phenomenon is caused by the strong p-π electronic interactions, that is, the empty p-orbital of boron center is partly filled by π-electron of the neighbor carbon atoms in 1-2, which are confirmed by the analysis of Laplacian of the electron density and natural bond orbitals. Furthermore, the negative electrostatic potentials of the boron centers in 1-2 also interpret that they are not the most preferred sites for incoming nucleophiles. Moreover, it is also found that the planar compounds 1-2 can act as promising electron transporting materials since the internal reorganization energies for electron are really small.

Theoretical study of geometrical and nonlinear optical properties of pyridinum N-phenolate betaine dyes

This paper presents an ab initio quantum chemical investigation of the geometrical structures and the non-linear optical properties (NLO) of three structural isomers of pyridinium N-phenolate betaine dye. The ground state geometrical parameters and the first-order hyperpolarizabilities were calculated using the Hartree-Fock (HF) as well as the second-order perturbation Møller-Pleset (MP2) method with the 6–31G, 6–31G(d), 6–31G(d,p), 6–31+G(d), 6–31++G(d,p), 6–311+G(d), aug-cc-PVDZ and the recently developed Z3PolX basis sets. Moreover, the first-order hyperpolarizability was calculated at the coupled cluster singles and doubles (CCSD/6–31+G(d)) level of theory. The analysis of the results of calculations for the investigated isomers indicates that there are important differences in their NLO activities. Additionally, it was shown that Z3PolX basis set works reasonable well for betaine dyes.

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.

Structures, spectroscopic and thermodynamic properties of U2On (n = 0 ∼ 2, 4) molecules: a density functional theory study

The equilibrium structures, spectroscopic and thermodynamic parameters [entropy (S), internal energy (E), heat capacity (C _p)] of U₂, U₂O, U₂O₂ and U₂O₄ uranium oxide molecules were investigated systematically using density functional theory (DFT). Our computations indicated that the ground electronic state of U₂ is the septet state and the equilibrium bond length is 2.194 Å; the ground electronic state of U₂O and U₂O₂ were found to be $ {\tilde{X}}^3\varPhi $ and $ {\tilde{X}}^3{\sum}_{\mathrm{g}} $ with stable C _∞v and D _∞h linear structures, respectively. The bridge-bonded structure with D _2h symmetry and $ {\tilde{X}}^3{\mathrm{B}}_{1\mathrm{g}} $ state is the most stable configuration for the U₂O₄ molecule. Mulliken population analyses show that U atoms always lose electrons to become the donor and O atoms always obtain electrons as the acceptor. Molecular orbital analyses demonstrated that the frontier orbitals of the title molecules were contributed mostly by 5f atomic orbitals of U atoms. Vibrational frequencies analyses indicate that the maximum absorption peaks stem from the stretching mode of U–O bonds in U₂O, U₂O₂ and U₂O₄. In addition, thermodynamic data of U₂O_n (n?=?0?～?4) molecules at elevated temperatures of 293.0 K to 393.0 K was predicted.     

Variation of the electronic dipole polarizability on the reaction path

The reaction force and the electronic flux, first proposed by Toro-Labbé et al. (J Phys Chem A 103:4398, 1999) have been expressed by the existing conceptual DFT apparatus. The critical points (extremes) of the chemical potential, global hardness and softness have been identified by means of the existing and computable energy derivatives: the Hellman-Feynman force, nuclear reactivity and nuclear stiffness. Specific role of atoms at the reaction center has been unveiled by indicating an alternative method of calculation of the reaction force and the reaction electronic flux. The electron dipole polarizability on the IRC has been analyzed for the model reaction HF + CO→HCOF. The electron polarizability determined on the IRC α _e (ξ) was found to be reasonably parallel to the global softness curve S(ξ). The softest state on the IRC (not TS) coincides with zero electronic flux.

B30H8, B39H9 2−, B42H10, B48H10, and B72H12: polycyclic aromatic snub hydroboron clusters analogous to polycyclic aromatic hydrocarbons

Calculations performed at the ab initio level using the recently reported planar concentric π-aromatic B₁₈H₆ ²⁺(1) [Chen Q et al. (2011) Phys Chem Chem Phys 13:20620] as a building block suggest the possible existence of a new class of B_3n H _m polycyclic aromatic hydroboron (PAHB) clusters—B₃₀H₈(2), B₃₉H₉ ^2?(3), B₄₂H₁₀(4/5), B₄₈H₁₀(6), and B₇₂H₁₂(7)—which appear to be the inorganic analogs of the corresponding C _n H _m polycyclic aromatic hydrocarbon (PAHC) molecules naphthalene C₁₀H₈, phenalenyl anion C₁₃H₉ ^?, phenanthrene/anthracene C₁₄H₁₀, pyrene C₁₆H₁₀, and coronene C₂₄H₁₂, respectively, in a universal atomic ratio of B:C?=?3:1. Detailed canonical molecular orbital (CMO), adaptive natural density partitioning (AdNDP), and electron localization function (ELF) analyses indicate that, as they are hydrogenated fragments of a boron snub sheet [Zope RR, Baruah T (2010) Chem Phys Lett 501:193], these PAHB clusters are aromatic in nature, and exhibit the formation of islands of both σ- and π-aromaticity. The predicted ionization potentials of PAHB neutrals and electron detachment energies of small PAHB monoanions should permit them to be characterized experimentally in the future. The results obtained in this work expand the domain of planar boron-based clusters to a region well beyond B₂₀, and experimental syntheses of these snub B_3n H _m clusters through partial hydrogenation of the corresponding bare B_3n may open up a new area of boron chemistry parallel to that of PAHCs in carbon chemistry.

Structural transitions in mixed ternary noble gas clusters

The properties of noble gas systems can be greatly extended by heterogeneous mixtures of elements. The geometrical structures and energies of mixed Ar–Kr–Xe clusters were investigated using ternary Lennard-Jones (TLJ) potential. For the Ar₁₉Kr _n Xe₁₉, Ar₁₉Kr₁₉Xe _n , and Ar _n Kr₁₉Xe₁₉ (n?=?0–17) clusters investigated, the results show that only two minimum energy configurations exist, i.e., polytetrahedron and six-fold pancake. The inner core of all these clusters is composed mainly of Ar atoms, and Kr and Xe atoms are distributed on the surface with well mixed pattern for polytetrahedral and segregate pattern for six-fold pancake configurations. The relative stability property of Ar–Kr–Xe clusters with a certain composition is discussed. Moreover, the role of heterogeneity on the strain was investigated, and reduced strain energies in Ar–Kr–Xe clusters were studied to find possible ways of reducing strain. The results showed that the strain energies were affected mainly by Ar–Ar, Ar–Kr, and Xe–Xe bonds.

Relativistic theoretical studies on hydrogen bonds and the electronic structure of aqueous solvated bis(uranyl) complex: an insight into explicit and/or implicit solvent effects

To understand the chemical behavior of uranyl complexes in water, a bis-uranyl [(phen)(UO₂)(μ₂–F)(F)]₂ (A; phen?=?phenanthroline, μ₂?=?doubly bridged) and its hydrated form A?·?(H₂O)_n (n?=?2, 4 and 6) were examined using scalar relativistic density functional theory. The addition of water caused the phen ligands to deviate slightly from the U₂(μ₂–F)₂ plane, and red-shifts the U–F-terminal and U?=?O stretching vibrations. Four types of hydrogen bonds are present in the optimized hydrated A?·?(H₂O)_n complexes; their energies were calculated to fall within the range 4.37–6.77 kcal mol^-1, comparable to the typical values of 5.0 kcal mol^-1 reported for hydrogen bonds. An aqueous environment simulated by explicit and/or implicit models lowers and re-arranges the orbitals of the bis-uranyl complex.

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.

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.