Co-reporter:Ling Lin;Xin Wang;Peter Lievens;G. Gopakumar;Jorg De Haeck
The Journal of Physical Chemistry C June 25, 2009 Volume 113(Issue 25) pp:10858-10867
Publication Date(Web):2017-2-22
DOI:10.1021/jp900950k
We report a combined experimental and theoretical study of lithium-doped germanium clusters. The doped germanium clusters are produced using a dual-laser dual-target vaporization source and subsequently analyzed using a reflectron time-of-flight (RTOF) mass spectrometer. The RTOF mass spectrum showing the relative abundance of the GenLim clusters, photoionized by an ArF laser, is subjected to a detailed theoretical investigation within the framework of the most popular density functional theory formalism. The periodic appearance of the peaks corresponding to Ge9Li4−5, Ge18Li5−7, Ge27Li9−12, Ge36Li13−16, and Ge45Li16−21 in the mass spectrum clearly suggests the existence of the Ge9 building blocks, which is well-known for the group 14 elements. On the basis of our complementary theoretical investigation, we were able to propose a theoretical model which rationalizes the experimental spectrum. Our systematic theoretical investigation supports the fact that doping with lithium atoms enhances the stability of the Ge cluster, where the Li atom mainly acts as charge balancer. We report for the first time the presence of oligomers of [Ge9]x— units in a gas-phase experiment.
Co-reporter:Vu Thi Thu Huong, Truong Ba Tai and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2016 vol. 18(Issue 8) pp:6259-6267
Publication Date(Web):27 Jan 2016
DOI:10.1039/C5CP07277H
Organic semiconducting materials play an important role in the fabrication of high performance organic electronic devices. In the present work, we theoretically designed a series of organic semiconductors based on nickel complexes. Their characteristics of charge transport were investigated using DFT computational approaches. Based on the computed results, all compounds designed are found to be excellent candidates for ambipolar organic semiconductors with low reorganization energies for both holes and electrons. The (I–V) characteristics and transmission spectra of materials show that the replacement of benzene rings by thiophene rings results in an increase of their HOMO and LUMO energy levels. HOMOs of compounds containing thiophene end-groups are likely dominant for their conductance, while LUMOs of compounds containing benzene end-groups mainly affect their conductance. The electron distribution in these frontier MOs is identified as the main reason which makes the conductance of the compounds in the first series higher than those in the later series.
Co-reporter:Truong Ba Tai, Sang Uck Lee and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2016 vol. 18(Issue 17) pp:11620-11623
Publication Date(Web):29 Feb 2016
DOI:10.1039/C5CP07342A
The cage-like structures containing octagonal holes are located as the lowest-lying isomers for the B0/+42. The presence of octagonal holes, which have been found for the first time, not only gives us new insight into the bonding motif, but also marks a breakthrough in the structural characteristics of boron clusters since they were never expected to be stable units for elemental clusters. These cages are composed of both delocalized σ and π electron systems that consequently make them aromatic and thermodynamically stable.
Co-reporter:Hung Tan Pham, Kie Zen Lim, Remco W. A. Havenith and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2016 vol. 18(Issue 17) pp:11919-11931
Publication Date(Web):25 Feb 2016
DOI:10.1039/C5CP07391J
The planarity of small boron-based clusters is the result of an interplay between geometry, electron delocalization, covalent bonding and stability. These compounds contain two different bonding patterns involving both σ and π delocalized bonds, and up to now, their aromaticity has been assigned mainly using the classical (4N + 2) electron count for both types of electrons. In the present study, we reexplored the aromatic feature of different types of planar boron-based clusters making use of the ring current approach. B3+/−, B42−, B5+/−, B6, B7−, B82−, B9−, B102−, B11−, B12, B13+, B142− and B162− are characterized by magnetic responses to be doubly σ and π aromatic species in which the π aromaticity can be predicted using the (4N + 2) electron count. The triply aromatic character of B12 and B13+ is confirmed. The π electrons of B182−, B19− and B202− obey the disk aromaticity rule with an electronic configuration of [1σ21π41δ42σ2] rather than the (4N + 2) count. The double aromaticity feature is observed for boron hydride cycles including B@B5H5+, Li7B5H5 and M@BnHnq clusters from both the (4N + 2) rule and ring current maps. The double π and σ aromaticity in carbon-boron planar cycles B7C−, B8C, B6C2, B9C−, B8C2 and B7C3− is in conflict with the Hückel electron count. This is also the case for the ions B11C5+/− whose ring current indicators suggest that they belong to the class of double aromaticity, in which the π electrons obey the disk aromaticity characteristics. In many clusters, the classical electron count cannot be applied, and the magnetic responses of the electron density expressed in terms of the ring current provide us with a more consistent criterion for determining their aromatic character.
Co-reporter:Yassin Aweis Jeilani, Phoenix N. Williams, Sofia Walton and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2016 vol. 18(Issue 30) pp:20177-20188
Publication Date(Web):10 May 2016
DOI:10.1039/C6CP02686A
The reaction pathways for the prebiotic formation of nucleobases are complex and lead to the formation of a mixture of products. In the past 50 years, there has been a concerted effort for identifying a unified mechanism for the abiotic origin of the biomolecules but with little success. In the present theoretical study, we identified two prominent precursors for the building up of RNA and DNA nucleobases under prebiotic conditions: (a) 1,2-diaminomaleonitrile (DAMN), which is a tetramer of hydrogen cyanide (HCN), and (b) formamide, a hydrolysis product of HCN; it is important to emphasize that HCN is the source of both precursors. We find that free radical pathways are potentially appropriate to account for the origin of nucleobases from HCN. The current study unites the formamide pathways with the DAMN pathways. The mechanisms for the formation of the RNA and DNA nucleobases (uracil, adenine, purine, cytosine) were studied by quantum chemical computations using density functional theory at the B3LYP/6-311G(d,p) level. All the routes involved proceed with relatively low energy barriers (within the error margin of DFT methods). We showed that the radical mechanisms for the formation of nucleobases could be unified through common precursors. The results demonstrated that 4-aminoimidazole-5-carbonitrile (AICN), which is a known precursor for nucleobases, is a product of DAMN. The overall mechanisms are internally consistent with the abiotic formation of the nucleobases, namely (a) under a meteoritic impact scenario on the early Earth's surface that generated high internal energy, and/or (b) in the (gas phase) interstellar regions without the presence of catalysts.
In a recent paper (Nanoscale, 2014, 6, 11692), based on the results computed using DFT and MP2 methods, the all-boron fullerene I was reported to be the global minimum of the cluster B38 and was much more stable than the quasi-planar II. In this comment, we have shown that at higher level of theory CCSD(T), both structure I and quasi-planar II are almost degenerate in energy and the B38 can be considered to be of a transition size between 2D and 3D boron clusters. While the MP2 method favours the 3D structure I, the CCSD method tends to overestimate the relative stability of the 2D structure II.
Chemical Communications 2015 vol. 51(Issue 36) pp:7677-7680
Publication Date(Web):24 Mar 2015
DOI:10.1039/C5CC01252J
The neutral B32 exhibits an aromatic bowl structure containing one heptagonal hole, while two anionic species, one having a bowl structure and the other a quasi-planar structure, are almost degenerate in energy. These findings not only give more insight into the structural features of boron clusters, but also present a key result explaining the presence of heptagonal holes in the fullerene B40.
Co-reporter:Nguyen Minh Tam, Hung Tan Pham, Long Van Duong, My Phuong Pham-Ho and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2015 vol. 17(Issue 5) pp:3000-3003
Publication Date(Web):16 Dec 2014
DOI:10.1039/C4CP04279D
Stabilized fullerene and tubular forms can be produced in boron clusters Bn in small sizes from n ∼ 14 to 20 upon doping by transition metal atoms. B14Fe and B16Fe are stable tubes whereas B18Fe and B20Fe are stable fullerenes. Their formation and stability suggest the use of dopants to induce different growth paths leading to larger cages, fullerenes and tubes of boron.
Physical Chemistry Chemical Physics 2015 vol. 17(Issue 20) pp:13672-13679
Publication Date(Web):22 Apr 2015
DOI:10.1039/C5CP01851J
Boron clusters have been of great interest over the last few decades due to their unique chemical and physical properties. In the present work, we performed a theoretical study of geometrical and electronic structures of boron clusters Bn with n = 26–29 in both neutral and anionic states using DFT and MO computational methods. The photoelectron spectra of anionic species were simulated using TDDFT methods. Our results predict that in the neutral state both the B26 and B27 clusters exhibit tubular forms, whereas the larger species B28 and B29 are quasi-planar structures. The anionic species Bn− are more favourable for 2D shapes. More importantly, based on known geometrical characteristics, we now establish a general growth mechanism of boron clusters, which gives us more insight into the formation and existence of boron based nanomaterials.
Physical Chemistry Chemical Physics 2015 vol. 17(Issue 26) pp:16927-16936
Publication Date(Web):01 Jun 2015
DOI:10.1039/C5CP01456E
Thermally feasible decomposition pathways of formamide (FM) in the presence of vanadium VO(X4Σ−) and titanium TiO(X3Δ) monoxides are determined using density functional theory (BP86 functional) and coupled-cluster theory (CCSD(T)) computations with large basis sets. These diatomic metal oxides have been shown to be present in the prebiotic conditions. The dehydration, decarbonylation and dehydrogenation reactions of the molecular and dissociative complexes of FM and MO (M = V, Ti) turn out to be more favourable than those of the ground state isolated FM. The effect of addition of one or two water molecules on energy barriers is also probed for these reaction pathways. In some cases, a combined catalytic effect when adding water is observed. This enhanced catalytic effect was not observed in previously reported cases of FM transformation, for example, when adding water molecules into the mineral-catalyzed isomerizations of FM. The dehydration process of MO–FM complexes without the presence of water is found to be more feasible than the decarbonylation and dehydrogenation. The overall energy barrier for the non-water VO–FM dehydration is ∼3 kcal mol−1 lower than the reference energy of the separated systems, whereas those of the two latter reactions are higher than the reference. Although the TiO–FM dehydration has a larger overall barrier of 14 kcal mol−1 as compared to the VO–FM counterpart, the two other decomposition pathways still have much higher energy barriers. Direct formation of urea and H2CO from a FM dimer and indirect formation of urea from FM via the intermediate HNCO are also established. Urea formation in an indirect pathway is preferred. These low-energy-barrier pathways leading to the formation of important prebiotic molecules suggest that metal monoxides MO could play an important catalytic role in the prebiotic reactions of FM.
Physical Chemistry Chemical Physics 2015 vol. 17(Issue 26) pp:17335-17345
Publication Date(Web):08 Jun 2015
DOI:10.1039/C5CP01650A
Using density functional theory with the TPSSh functional and the 6-311+G(d) basis set, we extensively searched for the global minima of two metallic atoms doped boron clusters B6M2, B7M2, B12M2 and B14M2 with transition metal element M being Co and Fe. Structural identifications reveal that B7Co2, B7Fe2 and B7CoFe clusters have global minima in a B-cyclic motif, in which a perfectly planar B7 is coordinated with two metallic atoms placed along the C7 axis. The B6 cluster is too small to form a cycle with the presence of two metals. Similarly, the B12 cluster is not large enough to stabilize the metallic dimer within a double ring 2 × B6 tube. The doped B14M2 clusters including B14Co2, B14Fe2 and B14CoFe have a double ring 2 × B7 tubular shape in which one metal atom is encapsulated by the B14 tube and the other is located at an exposed position. Dissociation energies demonstrate that while bimetallic cyclic cluster B7M2 prefers a fragmentation channel that generates the B7 global minimum plus metallic dimer, the tubular structure B14M2 tends to dissociate giving a bimetallic cyclic structure B7M2 and a B@B6 cluster. The enhanced stability of the bimetallic doped boron clusters considered can be understood from the stabilizing interactions between the anti-bonding MOs of metal–metal dimers and the levels of a disk aromatic configuration (for bimetallic cyclic structures), or the eigenstates of the B14 tubular form (in case of bimetallic tubular structure).
Co-reporter:Hung Tan Pham, Thu-Thuy Phan, Nguyen Minh Tam, Long Van Duong, My Phuong Pham-Ho and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2015 vol. 17(Issue 27) pp:17566-17570
Publication Date(Web):09 Jun 2015
DOI:10.1039/C5CP02257F
The smallest triple ring tubular silicon cluster Mn2@Si15 is reported for the first time. Theoretical structural identification shows that the Mn2@Si15 tubular structure whose triple ring is composed by three five-membered Si rings in anti-prism motif, is stable in high symmetry (D5h) and singlet ground state (1A1′). The dimer Mn2 is placed inside the tubular along the C5 axis, and the Mn dopant form single Si–Mn bonds with Si skeleton, whereas the Mn–Mn is characterized as a triple bond. The effect of Mn2 on the stability of the Si15 triple ring structure arises from strong orbital overlap of Mn2 with Si15.
Co-reporter:Vu Thi Thu Huong, Truong Ba Tai and Minh Tho Nguyen
RSC Advances 2015 vol. 5(Issue 31) pp:24167-24174
Publication Date(Web):24 Feb 2015
DOI:10.1039/C4RA16485G
In the present work, a series of new tetra-hetero[8]circulenes were theoretically designed in which heteroatoms include O, S, Se and N. Their electronic structure and characteristics of charge transport were investigated using DFT based computational methods. Except for the compounds containing Se-atoms (3a and 3b), all remaining compounds exhibit planar and highly symmetrical structures featuring novel aromatic features: the inner eight-membered ring is anti-aromatic and the outer fused rings are aromatic. The predicted UV spectrum and reduction potential of 1a agree well with available experimental values. Following replacement of H-atoms by F-atoms, the energy levels of the frontier orbitals of 1b–4b are consistently decreased as compared to those of 1a–4a. Based on the calculated properties of electrochemistry and charge transport, the molecules 1a–4a and 4b are suggested to be good candidates for p-type semiconductors. More importantly, the molecules 1b–3b are revealed to be potential ambipolar organic semiconductors.
The Journal of Physical Chemistry C 2015 Volume 119(Issue 9) pp:4524-4539
Publication Date(Web):February 11, 2015
DOI:10.1021/jp511668z
In previous papers (Nguyen et al. J. Phys. Chem. C2008, 112, 5662–5667 and J. Phys. Chem. C2009, 113, 18914–18926), formation of H2 molecules from ammonia alane monomer (AAl) and dimers (AAl)2 was shown to be facilitated by the addition of one or more alane or ammonia molecules that can play the role of efficient bifunctional catalyst. Ammonia alane emerges as a good starting compound for building up materials for chemical hydrogen storage (CHS). Further exploration of the products based on the H2 release from ammonia alane were carried out using coupled-cluster theory computations together with the aug-cc-pVTZ basis set (based on MP2/aug-cc-pVDZ optimized geometries). Our ab initio MO calculations for the first time led to the identification of cyclotrialazane [(H2AlNH2)3], alazine [(HAlNH)3], and its oligomer [H2Al(HNAlH)2NH2] that are produced along the multistep dehydrogenation processes from the reactions of ammonia alane and AlH3NH2AlH2NH3. The formation of alazine (homologue of borazine) as the final product in our H2 elimination reactions is an important feature because of its long-debated existence. The present reaction path analyses show that the formation of this compound is an important phenomenon for explaining the entire dehydrogenation process.
Co-reporter:Huyen Thi Nguyen, Yassin A. Jeilani, Huynh Minh Hung, and Minh Tho Nguyen
The Journal of Physical Chemistry A 2015 Volume 119(Issue 33) pp:8871-8883
Publication Date(Web):July 21, 2015
DOI:10.1021/acs.jpca.5b03625
The prebiotic formation of nucleobases, the building blocks of RNA/DNA, is of current interest. Highly reactive radical species present in the atmosphere under irradiation have been suggested to be involved in the prebiotic synthesis of nucleobases from formamide (FM). We studied several free radical reaction pathways for the synthesis of pyrimidine bases (cytosine, uracil, and thymine) from FM under cold conditions. These pathways are theoretically determined using density functional theory (DFT) computations to examine their kinetic and thermodynamic feasibilities. These free radical reaction pathways share some common reaction types such as H-rearrangement, •H/•OH/•NH2 radical loss, and intramolecular radical cyclization. The rate-determining steps in these pathways are characterized with low energy barriers. The energy barriers of the ring formation steps are in the range of 3–7 kcal/mol. Although DFT methods are known to significantly underestimate the barriers for addition of •H radical to neutral species, many of these reactions are highly exergonic with energy release of −15 to −52 kcal/mol and are thus favorable. Among the suggested pathways for formation of cytosine (main route, routes 7a and 1a), uracil (main route, routes 7b and 1b), and thymine (main route and route 26a), the main routes are in general thermodynamically more exergonic and more kinetically favored than other alternative routes with lower overall energy barriers. The reaction energies released following formation of cytosine, uracil, and thymine from FM via the main radical routes amount to −59, −81, and −104 kcal/mol, respectively. Increasing temperature induces unfavorable changes in both kinetic and thermodynamic aspects of the suggested routes. However, the main routes are still more favored than the alternative pathways at the temperature up to the boiling point of FM.
The Journal of Physical Chemistry A 2015 Volume 119(Issue 24) pp:6493-6503
Publication Date(Web):May 20, 2015
DOI:10.1021/acs.jpca.5b02492
The ground state geometries, electronic structures, and thermochemical properties of binary alkaline-earth-metal silicon clusters Si3M with M = Be, Mg, Ca in neutral, cationic, and anionic states were investigated using quantum chemical computations. Lowest-lying isomers of the clusters were determined on the basis of the composite G4 energies. Along with total atomization energies, thermochemical parameters were determined for the first time by means of the G4 and coupled-cluster theory with complete basis set CCSD(T)/CBS approaches. The most favored equilibrium formation sequences for Si3M clusters emerge as follows: all Si3M+/0/– clusters are formed by attaching the M atom into the corresponding cation, neutral and anion silicon trimer Si3+/0/–, except for the Si3Mg+ and Si3Ca+ where the metal cations are bound to the neutral Si3. The resulting mixed tetramers exhibit geometrical and electronic features similar to those of the pure silicon tetramer Si4+/0/–. Electron localization function (ELF) and ring current analyses point out that the σ-aromatic character of silicon tetramer remains unchanged upon substituting one Si atom by one alkaline-earth-metal atom.
Co-reporter:Truong Ba Tai, Long Van Duong, Hung Tan Pham, Dang Thi Tuyet Mai and Minh Tho Nguyen
Chemical Communications 2014 vol. 50(Issue 13) pp:1558-1560
Publication Date(Web):03 Dec 2013
DOI:10.1039/C3CC48392D
The B30 boron cluster has a bowl rather than a double-ring or a triple-ring tubular structure. This bowl isomer exhibits disk-aromaticity similar to that found for B202− and B19− clusters. We confirmed that the concept of disk-aromaticity can be applied to both planar and non-planar systems.
Co-reporter:Long Van Duong, Hung Tan Pham, Nguyen Minh Tam and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2014 vol. 16(Issue 36) pp:19470-19478
Publication Date(Web):29 Jul 2014
DOI:10.1039/C4CP01996B
We determined the geometries and chemical bonding phenomena of the B27 system in its dicationic, cationic, neutral, anionic and dianionic states using DFT computations. In both cationic and neutral states, the triple ring tubular forms correspond to the lowest-energy isomers, especially in B27+. The cation B27+ represents the first stable hollow cylinder having a triple ring among the pure boron clusters. In the anionic and dianionic states, the quasi-planar structures are favoured due to a charge effect. In the triple ring tube B27+, strong diatropic responses to external magnetic field occur in both radial and tangential types of electrons, and thus confer it a characteristic tubular aromaticity. The presence of a consistent aromatic character contributes to its high thermodynamic stability. The shapes of calculated MOs of B27+ TR can be predicted by the eigenstates of a simple model of a particle on a hollow cylinder. The number of electrons in a hollow cylinder should attain a number of (4N + 2M) with M = 0 and 1 for both radial and tangential electrons, depending on the number of non-degenerate MOs occupied, in order to properly fulfill the closed electron shells. In the case of B27+, M = 0 for radial electrons and M = 1 for tangential electrons.
Co-reporter:Yassin A. Jeilani, Huyen Thi Nguyen, Beatriz H. Cardelino, Minh Tho Nguyen
Chemical Physics Letters 2014 Volume 598() pp:58-64
Publication Date(Web):8 April 2014
DOI:10.1016/j.cplett.2014.02.053
•A unified mechanism is proposed for the formation of xanthine and isoguanine.•The mechanisms explain formation of biomolecules and promote prebiotic biodiversity.•The mechanisms are appropriate for non-aqueous scenarios of prebiotic synthesis.Free radical pathways for the synthesis of xanthine and isoguanine from formamide were studied using density functional theory (B3LYP/6-311G(d,p)). The proposed mechanisms are complex and appropriate for the non-aqueous scenario of prebiotic reactions. Formation of the carbonyl bond in the nucleobases proceeds through enol-keto tautomerization since the direct formation of the CO bond is a highly endothermic step. The mechanisms show 2-amino-imidazole as a precursor for nucleobases and polyazaporphyrin. The proposed mechanisms contribute to a further understanding of the origin of biomolecules.
Co-reporter:Nguyen Minh Tam, Vu Thi Ngan, Minh Tho Nguyen
Chemical Physics Letters 2014 Volumes 595–596() pp:272-276
Publication Date(Web):18 March 2014
DOI:10.1016/j.cplett.2014.02.015
•CSi42+andCGe42+ dications contain a planar tetracoordinate carbon atom (ptC).•The neutral CSi9 and CGe9 clusters also exhibit a ptC.•The driving force for C-planarization is electron delocalization along C–X bonds.•The X5 group stabilizes the cage by large electron transfer in maintaining a ptC.•A ptC can be stabilized in a cluster using both electronic and mechanical effects.Using quantum chemical computations and analysis of electron distribution (MO, DOS, ELF) we showed that in some carbon-doped silicon and germanium clusters, it is possible to achieve a planar tetracoordinate carbon with enhanced stability. While the driving force for C-planarization in the square dications CX42+ (with X = Si, Ge) is electron delocalization on X4 frame together with single bonds along C–X bonds, the larger neutral CSi9 and CGe9 clusters enjoy combined stabilization from both electronic effect and geometrical constraint of the X9 cages. In CX9, an additional electrostatic interaction reinforces stabilization within the CX4 moiety in maintaining the ptC configuration.A ptC in CSi9 and CGe9.
The Journal of Physical Chemistry A 2014 Volume 118(Issue 34) pp:7017-7023
Publication Date(Web):July 29, 2014
DOI:10.1021/jp5053216
The effects of kaolinite mineral surfaces on the unimolecular rearrangements of formamide (FM) were investigated using periodic density functional theory in conjunction with pseudopotential plane-wave approach. Surface hydroxyl groups covering the octahedral surface of kaolinite were found to play the role of catalysts in the transformations of FM. They induce a reduction of 31 kcal/mol on the energy barrier for formation of its isomer aminohydroxymethylene (AHM), which is close to the reduction amount calculated for water-catalyzed reactions. This suggests that the kaolinite octahedral surface exerts a catalytic effect similar to that of the water molecule. As the tetrahedral surface does not contain catalytic surface hydroxyl groups, only water-assisted FM transformation was therefore studied on this surface whose energy barrier amounts to ∼17 kcal/mol. The combined effect of both water and kaolinite on FM rearrangements via triple hydrogen transfer reactions does not significantly lower the energy barriers, as compared to those of double hydrogen transfer reactions. The triple hydrogen transfer energy barriers amount to ∼20 and ∼36 kcal/mol, and the double ones are ∼21 and ∼40 kcal/mol for formation of formimic acid and AHM isomers, respectively. However, the energies of the systems in water-catalyzed channels lie below the available energies of the original reactants, and thus these channels are more favored than the water-free ones. With its multiple functions as both a supporting plate-form and a catalyst for FM reactions, kaolinite can thus be regarded as an important natural catalyst for prebiotic synthesis.
The Journal of Physical Chemistry A 2014 Volume 118(Issue 23) pp:4079-4086
Publication Date(Web):May 15, 2014
DOI:10.1021/jp5013945
The efficient formation of HCN/HNC from formamide (FM) combining the advantages of water-assistance, self-catalyzed reactions, and the mineral surfaces was investigated. Periodic density functional theory calculations with plane-wave pseudopotential basis sets were performed to study the interaction of FM with pyrite (100) ideal and defect surfaces. Effects of sulfur vacancy defect and water on tautomerization and rearrangement barriers of FM on the (100) surface were evaluated. Calculated results show that FM adsorbs more strongly on the defect surface than on the ideal surface, with the lowest adsorption energy on the defect surface being −22 kcal/mol. The energy barriers for rearrangements of FM on these two surfaces being close to each other suggests that the adsorptions on the surfaces have small effects on the energy barriers. The energy barriers for formimic acid isomer formations are 44.5 and 46.0 kcal/mol, and those of aminohydroxymethylene formations are 72.6 and 71.9 kcal/mol on the ideal and defect surfaces, respectively. A reduction of ∼30 kcal/mol in tautomerization energy barriers is observed in water-assisted process on the defect surface. Because this reduction is close to that of the gas-phase reactions, the catalytic effect is clearly due to the presence of water molecule instead of the interaction with the surface. In this case, the pyrite surfaces with the ability to accumulate reactive species only play the role of connecting bridges between the two steps of the proposed reaction mechanism: the water-assisted rearrangement and the self-catalyzed dehydration.
Co-reporter:Vu Thi Thu Huong, Truong Ba Tai, and Minh Tho Nguyen
The Journal of Physical Chemistry A 2014 Volume 118(Issue 18) pp:3335-3343
Publication Date(Web):April 10, 2014
DOI:10.1021/jp500899k
The characteristics of molecular structure and charge transport of some new n-type organic semiconductors containing thiazole 1a–6a and oxazole 1b–6b frameworks and trifluoromethylphenyl as terminal groups were predicted using density functional theory (DFT) methods. The energy levels of HOMO and LUMO of these compounds are decreased when thiophene and furan units are replaced by thiazole and oxazole units, respectively. The same trend was observed when benzo[1,2-d:4,5-d′]bisthiazole groups were replaced with benzo[1,2-d:4,5-d′]bisthiazole-4,8-diones. The reorganization energies for electron of compounds are computed in a range of 0.21–0.37 eV, which is comparable to the value of 0.25 eV of well-known n-type semiconductors such as perfluoropentacene. Some important trends can be pointed out as follows: (i) replacing the core thiazolothiazole unit of compounds 1a and 2a by the larger core benzo[1,2-d:4,5-d′]bisthiazole units of 3a and 4a decreases both reorganization energies for electron (λe); (ii) the λe values of compounds containing thiazole 2a, 4a, and 6a are smaller than those of compounds containing thiophene 1a, 3a, and 5a, respectively; (iii) there is no clear trend when replacing benzene rings of compounds 3a and 4a by quinone rings of 5a and 6a. The λe values of 5 and 6 are only somewhat larger. The same trend is also found for compounds containing oxazole 1b–6b. The intermolecular charge transports in solid state of these compounds mainly occur among molecules in the same molecular layer, whereas intermolecular interactions between molecules in different molecular layers are very small. Generally, beside some experimentally reported molecules 1a–4a, the remaining molecules designed here are good candidates for n-type organic semiconducting materials with small reorganization energies for electron and low energy levels of LUMO.
The Journal of Physical Chemistry C 2014 Volume 118(Issue 41) pp:24181-24187
Publication Date(Web):September 26, 2014
DOI:10.1021/jp507901n
We analyzed the chemical bonding phenomena of the boron B2n tubes (n = 10–14). The B2n tubes represent stable hollow cylinders each having a double ring (DR) among pure boron clusters. The shapes of the molecular orbitals of a B2n DR can be predicted by the eigenstates of a simple model of a particle on a hollow cylinder, which shows both radial and tangential components. In a DR tube, strong diatropic responses to external magnetic field occur in both radial and tangential types of electrons, and thus confer it a characteristic tubular aromaticity. The presence of a consistent aromatic character contributes to the high thermodynamic stability of a DR. The number of electrons in a hollow cylinder should attain (4N + 2M) with M = 0 and 1 for both series of radial and tangential electrons, depending on the number of nondegenerate MOs occupied, to properly fill the electron shells. In the case of B20, M = 1 for both radial and tangential electrons, and the classical Hückel counting rule is thus recovered.
Co-reporter:Truong Ba Tai, Vu Thi Thu Huong and Minh Tho Nguyen
Chemical Communications 2013 vol. 49(Issue 98) pp:11548-11550
Publication Date(Web):17 Oct 2013
DOI:10.1039/C3CC47573E
The heteropolycyclic compounds containing borole units were theoretically designed. The presence of electron deficient boron atoms results in full electron delocalization and remarkably affects their aromaticity. While molecules 1 and 2a exhibit antiaromaticity for inner rings and non-aromaticity for outer rings, 2b and 2c are completely aromatic.
Chemical Communications 2013 vol. 49(Issue 9) pp:913-915
Publication Date(Web):10 Dec 2012
DOI:10.1039/C2CC38038B
The B6Li8 cluster is a symmetrical 3D complex whose high stability can be understood through the Wade rule and aromaticity. A new mechanism of B–Li chemical bonding is proposed. Importantly, B6Li8 is predicted to be a promising candidate for hydrogen storage material with gravimetric density reaching up to a theoretical limit of 24%.
The B202– cluster is predicted to exhibit a planar sheet-like structure with a circular circumference. Orbital plots and energy correlations demonstrate the close correspondence between the electronic structure of B202– and the Bessel functions describing the waves of a quantum mechanical particle confined to a disk. The π-band of B202–, and its B19– congener, contains 12 π-electrons, forming a (1σ)2(1π)4(1δ)4(2σ)2 configuration, which corresponds to a “disk aromaticity” electron count. The analogy not only applies to the π-band, but also extends to the 50 valence σ-electrons. The occupied σ-orbitals are assigned on the basis of radial and angular nodes of the scalar disk waves. The magnetic response of the cluster was examined by Nucleus Independent Chemical Shift (NICS) values and current density calculations based on the ipsocentric model. B202– is found to exhibit a remarkable inner paratropic current in the σ-channel and an outer diatropic current in the π-channel. The orbital excitations responsible for the antiaromaticity in σ and the disk-aromaticity in π are identified.
Physical Chemistry Chemical Physics 2013 vol. 15(Issue 37) pp:15404-15415
Publication Date(Web):10 Jul 2013
DOI:10.1039/C3CP51017D
Optical properties of silver Agn nanoclusters are demonstrated to be dependent on their size, structure and charge state. It is found that when being contained in the sodalite cavity of LTA zeolite the tetradecanuclear hexacation silver cluster Ag146+ is stable. Its lower-lying states and optical spectrum are theoretically determined using the quantum chemical TD-DFT method. Its ground state possesses an outer-shell electron configuration of A1g2T2g6 mimicking the s2p6 valence of noble gas atoms. These frontier orbitals are constructed from 5s,5p(Ag)-AOs with contributions from framework oxygen atoms. Light absorption of Ag146+ embedded in the sodalite cage which is characterized by strong peaks centered at 331 and 476 nm (transitions 5s,p(Ag) → 5s,p(Ag)) leads to much longer wavelength emission. The sodalite cage, as a container, stabilizes the central Ag146+ cluster by electrostatic attraction. The absorption spectrum of the isovalent neutral Ag8 cluster embedded inside the same sodalite cavity is also simulated using TD-DFT and CASPT2 methods. This absorption spectrum which is similar to that of the Ag146+ cluster has two absorption bands in the near UV and visible regions.
Co-reporter:Yassin A. Jeilani, Huyen Thi Nguyen, Domnique Newallo, Jean-Marie D. Dimandja and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2013 vol. 15(Issue 48) pp:21084-21093
Publication Date(Web):25 Oct 2013
DOI:10.1039/C3CP53108B
Modeling the complicated chemical reactions in the interstellar medium and surface materials of Titan is nontrivial. Since both the atmosphere and the surface are rich in organic molecules, the chemistry may have important implications for the origin of biomolecules. Prebiotic synthesis of DNA nucleobases from simple molecules such as formamide has been known for more than half a century. In this study, new free radical pathways leading to the synthesis of guanine, hypoxanthine, purine, and adenine have been studied using density functional theory (B3LYP with the 6-311G(d,p) basis set). The pathways of the selected nucleobases demonstrate the importance of free radicals in the production of useful biomolecules under conditions appropriate for the interstellar medium or on Titan. The pathways may be universal in nature and proceed without solvent requirements. Calculations indicate that radical pathways yield lower reaction barriers as compared to previously reported pathways. Overall, these results suggest that the chemistry on Titan's surface and/or the growth of organic particulates in the haze layers in Titan's atmosphere likely involve free radicals. The mechanisms demonstrate that important prebiotic precursors can be predicted. The reaction sequences reported here may lead to the production and build-up of molecules with prebiotic relevance.
Co-reporter:Hung Tan Pham, Long Van Duong, Buu Quoc Pham, Minh Tho Nguyen
Chemical Physics Letters 2013 Volume 577() pp:32-37
Publication Date(Web):9 July 2013
DOI:10.1016/j.cplett.2013.05.041
•The 2D–3D relative energy of small boron cluster (Bn, n = 20–24) linearly varies with respect to the actual net charge.•Electron addition favors 2D planar structures.•An early 2D-to-3D geometry hopping occurs for cationic boron clusters, and reversely for anions.•Planarization of molecular skeleton is a driving force facilitating electron delocalization in 2D structures.DFT TPSSh/6-311+G(d) calculations are carried out on a series of 2D and 3D forms of Bn, n = 20, 22 and 24 in different charge states. For a certain size, the relative energy within a pair of two-dimensional quasi-planar (2D) and three-dimensional staggered double-ring (3D) boron cluster isomers may shift the sign as they reach a certain charge state. Specifically, electron addition tends to enhance the stability of the 2D over the corresponding 3D isomer irrespective of the available electrons. Linear correlations between 2D–3D relative energy and net charge are established. Along with 2D-to-3D geometry hopping at critical size, our results suggest a local 2D–3D geometry hopping via critical charge.
Chemical Physics Letters 2013 Volume 584() pp:147-154
Publication Date(Web):1 October 2013
DOI:10.1016/j.cplett.2013.08.051
•Total atomization energies and heats of formation of Sin clusters are determined.•Patterns of EA, IE, binding and dissociation energies of Sin clusters are followed.•Deviations between computed and available experimental TEAs are large.Total atomization energies and heats of formation of small silicon clusters Sin and their ions are calculated using G4 (n = 2–13) and CCSD(T)/CBS (aug-cc-pV(n+d)Z for n = 2–6) methods. Experimental data for Sin were available with large uncertainties. A new ground state structure for Si11- was located. Using ΔfH∘(Si,298K)=451.5 kJ/mol, ΔfH∘(Sin,298K) are computed as: Si2: 588/588 kJ/mol (G4/CBS), Si3: 625/632, Si4: 633/639, Si5: 669/692, Si6: 675/701, Si7: 698, Si8: 866, Si9: 872, Si10: 833, Si11: 996, Si12: 1051 and Si13: 1158. Adiabatic electron affinities, ionization, binding and dissociation energies of Sin are determined.
Co-reporter:Vinh Son Nguyen, D. Majumdar, Jerzy Leszczynski, Minh Tho Nguyen
Chemical Physics Letters 2013 Volume 584() pp:30-36
Publication Date(Web):1 October 2013
DOI:10.1016/j.cplett.2013.08.014
•Energy barriers for H2 release from XH3PH3, X = B, Al are higher than bond energies.•Small hydrides (borane, alane, galane) act as bifunctional catalysts for H2 release.•Alane reduces the energy barrier for H2 release in BH3PH3 down to 20 kcal/mol.•In AlH3PH3 + AlH3, the TS for H2 release is located even below the starting reactants.The H2 release mechanism from phosphine borane and phosphine alane was investigated using quantum chemical methods (MP2/aug-cc-pVTZ geometry optimization and coupled-cluster energies were obtained through complete basis set extrapolation, CCSD(T)/CBS). The effect of catalysts borane, alane and galane on the processes was also explored. As the energy barriers for the release of H2 from BH3PH3 and AlH3PH3 are much higher than the B–P and Al–P bond energies, the presence of inherent catalysts can reduce substantially such energy barriers (using BH3 for BH3PH3, while AlH3 and GaH3 for AlH3PH3), and these systems could be useful as probable hydrogen source.
Co-reporter:Pham Cam Nam, Asit K. Chandra, Minh Tho Nguyen
Chemical Physics Letters 2013 Volume 555() pp:44-50
Publication Date(Web):3 January 2013
DOI:10.1016/j.cplett.2012.10.072
Abstract
Integration of the (RO)B3LYP/6-311++G(2df,2p) with the PM6 method into a two-layer ONIOM is found to produce reasonably accurate BDE(O–H)s of phenolic compounds. The chosen ONIOM model contains only two atoms of the breaking bond as the core zone and is able to provide reliable evaluation for BDE(O–H) for phenols and tocopherol. Deviation of calculated values from experiment is ±(1–2) kcal/mol. BDE(O–H) of several curcuminoids and flavanoids extracted from ginger and tea are computed using the proposed model. The BDE(O–H) values of enol curcumin and epigallocatechin gallate are predicted to be 83.3 ± 2.0 and 76.0 ± 2.0 kcal/mol, respectively.
Co-reporter:Vinh Son Nguyen, Thomas M. Orlando, Jerzy Leszczynski, and Minh Tho Nguyen
The Journal of Physical Chemistry A 2013 Volume 117(Issue 12) pp:2543-2555
Publication Date(Web):March 5, 2013
DOI:10.1021/jp312853j
Formamide (NH2CHO, FM) has been considered an active key precursor in prebiotic chemistry on early Earth. Under certain conditions such as dry lagoons, FM can decompose to produce reactants that lead to formation of more complex biomolecules. Specifically, FM decomposition follows many reactive channels producing small molecules such as H2, CO, H2O, HCN, HNC, NH3, and HNCO with comparable energy barriers in the range of 73–82 kcal/mol. Due to the likely presence of water on prebiotic Earth and the intrinsic presence of water following FM decomposition, we explore the effects of water oligomers, (H2O)n with n = 1–3, on its dehydration, dehydrogenation, and decarbonylation reactions using quantum chemical computations. Geometries are optimized using MP2/aug-cc-pVxZ calculations (x = D,T), and relative energies are evaluated using coupled-cluster theory CCSD(T) with the aug-cc-pVxZ basis sets (x = D, T, Q). Where possible the coupled-cluster energies are extrapolated to the complete basis set limit (CBS). Water classically acts as an efficient bifunctional catalyst for decomposition. With the presence of one water molecule, the dehydration pathway leading to HCN is favored. When two and three water molecules are involved, dehydration remains energetically favored over other channels and attains an energy barrier of ∼30 kcal/mol.
Co-reporter:Huyen Thi Nguyen, Vinh Son Nguyen, Nguyen Tien Trung, Remco W. A. Havenith, and Minh Tho Nguyen
The Journal of Physical Chemistry A 2013 Volume 117(Issue 33) pp:7904-7917
Publication Date(Web):July 26, 2013
DOI:10.1021/jp405657y
Unimolecular decompositions of neutral (NH2CHO) and protonated (NH3CHO+) formamide, an active precursor of biomolecules in prebiotic chemistry, are investigated in the ground (S0) and first triplet (T1) and singlet (S1) excited states. Different decomposition channels including the homolytic bond dissociations, dehydration, decarbonylation, dehydrogenation, etc., are explored using coupled-cluster theory (CCSD(T)/CBS method) for both S0 and T1 states and RASPT2(18,15)/6-31G(d,p) computations for the S1 state. On S1 and T1 energy surfaces, formamide preferentially follows C–N homolytic bond cleavages forming NH2 + HCO radical pairs. Formation of HCN and HNC from dehydration of neutral and protonated formamide via formimic acid and aminohydroxymethylene isomers has higher energy barriers. A strong stabilization upon triplet excitation of the two latter isomers significantly facilitates the interconversions between isomers, and thus considerably reduces the energy barriers for dehydration pathways. The most probable pathways for HCN and HNC generation are found to be dehydration of formamide in the T1 state. Dehydration pathways from the neutral S1 and protonated T1 forms lead to stable complexes of HCN and HNC with water but are associated with large energy barriers. Overall, in the lower-lying excited states of either neutral or protonated formamide, dehydration is not competitive with homolytic C–N bond cleavages, which finally lead to formation of CO.
Co-reporter:Nguyen Minh Tam, Truong Ba Tai, Vu Thi Ngan, and Minh Tho Nguyen
The Journal of Physical Chemistry A 2013 Volume 117(Issue 31) pp:6867-6882
Publication Date(Web):July 9, 2013
DOI:10.1021/jp405280c
A systematic examination of the aluminum doped silicon clusters, SinAlm with n = 1–11 and m = 1–2, in both neutral and anionic states, is carried out using quantum chemical calculations. Lowest-energy equilibrium structures of the clusters considered are identified on the basis of G4 energies. High accuracy total atomization energies and thermochemical properties are determined for the first time using the G4 and CCSD(T)/CBS (coupled-cluster theory with complete basis set up to n = 3) methods. In each size, substitution of Si atoms at different positions of a corresponding pure silicon clusters by Al dopants invariably leads to a spectrum of distinct binary structures but having similar shape and comparable energy content. Such an energetic degeneracy persists in the larger cluster sizes, in particular for the anions. The equilibrium growth sequences for Al-doped Si clusters emerge as follows: (i) neutral singly doped SinAl clusters favor Al atom substitution into a Si position in the structure of the corresponding cation Sin+1+, whereas the anionic SinAl– has one Si atom of the isoelectronic neutral Sin+1 being substituted by the Al impurity; and (ii) for doubly doped SinAl20/– clusters, the neutrals have the shape of Sin+1 counterparts in which one Al atom substitutes a Si atom and the other Al adds on an edge or a face of it, whereas the anions have both Al atoms substitute two Si atoms in the Sin+2+ frameworks. The Al dopant also tends to avoid high coordination position.
Co-reporter:Truong Ba Tai, Vu Thi Thu Huong, and Minh Tho Nguyen
The Journal of Physical Chemistry C 2013 Volume 117(Issue 29) pp:14999-15008
Publication Date(Web):July 3, 2013
DOI:10.1021/jp4049154
Heteropolycyclic compounds have been emerging as promising materials for electronic organic devices. In this context, structural features and charge transport properties of triarylborane 1, which was experimentally synthesized, and its derivatives 2–5 are investigated using ab initio calculations. The FMO energy levels predicted by using the SMD/IEF-PCM solvation model reveal that all compounds 1–5 have high oxidation stability. Triarylborane 1 possesses good hole and electron transport characters with hole mobility of 2.97 × 10–2 cm2 V–1 s–1 and electron mobility of 1.96 cm2 V–1 s–1. Substitution of −C(CH3)2 groups of 1 by functional groups of −NH, −O and −S reduces reorganization energies for hole and increases rates of hole hopping and lowers energy barriers for gold electrode of derivatives 3–5. On the basis of the predicted results, we suggest that these compounds 1–5 are good candidates for efficient p-type semiconductors.
Co-reporter:Vu Thi Thu Huong, Huyen Thi Nguyen, Truong Ba Tai, and Minh Tho Nguyen
The Journal of Physical Chemistry C 2013 Volume 117(Issue 19) pp:10175-10184
Publication Date(Web):April 17, 2013
DOI:10.1021/jp401191a
We performed a theoretical investigation on a series of π-conjugated organic molecules containing naphtho[2,3-b]thiophene and their derivatives using density functional theory calculations. All molecules considered exhibit planar structures and aromaticity. Energy levels of frontier orbitals and reduction and oxidation potentials of these compounds predicted by our solvation model reveal good agreement with available experimental values. The UV absorption spectra point out a clear trend that maximum peaks corresponding HOMO–LUMO transitions are red-shifted: (i) from compounds containing O to those containing Se, (ii) from dimers 1a–3a and 1b–3b to trimers 4a–6a and 4b–6b, and (iii) from parent compounds 1a–6a to perfluorinated derivatives 1b–6b. Parent compounds 1a–6a can be considered as p-type semiconducting materials with low reorganization energies, high transfer integrals, and hole mobility. Perfluorinated compounds 1b–6b are suggested to be very good candidates for ambipolar semiconducting materials. Introduction of fused-ring core molecules considerably improves the charge transport characteristics of the co-oligomers 4a–6a and 4b–6b as compared to those of corresponding molecules 1a–3a and 1b–3b. Accordingly, the former have lower reorganization energies, higher electron transfer integrals, and higher rates of charge hopping.
Co-reporter:Vinh Son Nguyen, Rehab M. Ibrahim Elsamra, Jozef Peeters, Shaun A. Carl and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2012 vol. 14(Issue 20) pp:7456-7470
Publication Date(Web):29 Mar 2012
DOI:10.1039/C2CP40367F
We investigated the rate constants and reaction mechanism of the gas phase reaction between the ethynyl radical and nitrous oxide (C2H + N2O) using both experimental methods and electronic structure calculations. A pulsed-laser photolysis/chemiluminescence technique was used to determine the absolute rate coefficient over the temperature range 570 K to 836 K. In this experimental temperature range, the measured temperature dependence of the overall rate constants can be expressed as: k(T) (C2H + N2O) = 2.93 × 10−11 exp((−4000 ± 1100) K/T) cm3 s−1 (95% statistical confidence). Portions of the C2H + N2O potential energy surface (PES), containing low-energy pathways, were constructed using the composite G3B3 method. A multi-step reaction route leading to the products HCCO + N2 is clearly preferred. The high selectivity between product channels favouring N2 formation occurs very early. The pathway corresponds to the addition of the terminal C atom of C2H to the terminal N atom of N2O. Refined calculations using the coupled-cluster theory whose electronic energies were extrapolated to the complete basis set limit CCSD(T)/CBS led to an energy barrier of 6.0 kcal mol−1 for the entrance channel. The overall rate constant was also determined by application of transition-state theory and Rice–Ramsperger–Kassel–Marcus (RRKM) statistical analyses to the PES. The computed rate constants have similar temperature dependence to the experimental values, though were somewhat lower.
Co-reporter:Vu Thi Thu Huong, Truong Ba Tai and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2012 vol. 14(Issue 43) pp:14832-14841
Publication Date(Web):02 Aug 2012
DOI:10.1039/C2CP42474F
Following the theme of this special issue, two new compounds, the P-flowers C16(PH)8 and C16(PF)8, are designed by us and subsequently characterized by quantum chemical computations. Their geometries and infrared signatures are analyzed and compared to those of the well-known sulflower C16S8. Their electronic structure and aromaticity are examined using the electron localization function (ELF) and also by the total and partial densities of state (DOS). Both C16(PF)6 and C16(PH)8 molecules exhibit small energy barrier of electron injection (Φe = 0.33 eV for the gold electrode for the former, and Φe = 0.1 eV for the calcium electrode for the latter), remarkably low reorganization energy and high rate of electron hopping. Thus, both theoretically designed P-flower molecules are predicted to be excellent candidates for organic n-type semiconductors.
Properties of a series of MnXn with X = H, F, Cl, Br and n = 1–4 are investigated using DFT, CCSD(T) and CASPT2 computations. The B3P86/6-311++G(3df,2d) method appears to be suitable for predicting their structures whose geometries and IR spectra are dependent on the charge state. While MnX2 are linear, MnX3 and MnX4 are characterized by high symmetry shape. The π-bonding type is observed for MnH30/+ and MnH40/+. In halides, a different type of bonds is formed as p-orbitals of halogens can overlap with empty metal d-orbitals allowing a more effective electron transfer and high spin ground electronic states. Vibrational frequencies and basic energetic quantities are computed and compared with available experiments. Several previous thermochemical quantities are re-evaluated, and the heats of formation of Mn-compounds can be determined with reasonable accuracy using the B3LYP functional. However, while calculated ionization energies are in agreement with experiment, electron affinities are obtained with large deviations.Graphical abstractHighlights► The B3P86 functional is found to be reliable in predictions of molecular structures and vibrational spectra. ► The hybrid B3LYP is more reliable for energetic parameters such as heats of formation. ► We also propose several new assignments for heats of formation and ionization energies of a number of species considered.
Co-reporter:Nguyen Minh Tam, Truong Ba Tai, and Minh Tho Nguyen
The Journal of Physical Chemistry C 2012 Volume 116(Issue 37) pp:20086-20098
Publication Date(Web):August 27, 2012
DOI:10.1021/jp306037q
A systematic investigation of the boron-doped silicon clusters SinB with n ranging from 1 to 10 in the neutral, anionic, and cationic states is performed using quantum chemical calculations. Lowest-energy minima of the clusters considered are identified on the basis of the B3LYP, G4, and CCSD(T) energies. Total atomization energies and thermochemical properties such as ionization energy, electron affinity, and dissociation energies are obtained using the high accuracy G4 (B3LYP-MP4-CCSD(T)) and CCSD(T)/CBS (complete basis set up to n = 4) methods. Theoretical heats of formation are close to each other and used to assess the available experimental values. The growth mechanism for boron-doped silicon clusters SinB with n = 1–10 emerges as follows: (i) each SinB cluster is formed by adding one excess Si-atom into the smaller sized Sin–1B, rather than by adding B into Sin, (ii) a competition between the exposed (exohedral) and enclosed (endohedral) structures occurs at the size Si8B where both structures become close in energy, and (iii) the larger size clusters Si9B and Si10B exhibit endohedral structures where the B-impurity is located at the center of the corresponding Sin cages. The species Si9B–, Si9B, and Si10B+ are identified as enhanced stability systems with larger average binding energies and embedded energies. The higher stability of the closed shells Si9B– and Si10B+ can be rationalized in terms of the jellium electron shell model and spherical aromaticity.
The Journal of Physical Chemistry A 2012 Volume 116(Issue 27) pp:7405-7418
Publication Date(Web):June 11, 2012
DOI:10.1021/jp302279j
A comprehensive theoretical investigation on structures and properties of niobium clusters in the range from 13 to 20 atoms, in three different charged states, is performed by using the BPW91 and M06 functionals and the cc-pVDZ-PP basis set. These species are predicted to prefer low spin ground state, i.e., singlet (for even electron) and doublet (for odd electron) systems. In terms of growth mechanism, a compact structure with one Nb encapsulated by a cage formed from five and six triangles is found to be favored over an icosahedral evolution. Unlike many 3d metals, whose volumes are much smaller, 13 and 19 Nb atoms clusters do not exist as icosahedra and double-icosahedra. A distinct case is Nb15 as it bears a slightly distorted bcc structure. For some systems, several lower lying isomers are computed to be so close in energy that DFT computations cannot clearly establish their ground electronic states. The existence of structural isomers with comparable energy content is established for Nbn species with n = 13, 18, 19, and 20 in both neutral and charged states. The vibrational (IR) spectra are also calculated. While the spectra of smaller systems are strongly dependent on addition or removal of an electron from the neutral, the spectra of the larger size clusters are mostly independent of the charged state. The neutrals and their corresponding ions usually have a quite similar IR pattern. Electron affinities (EA), ionization energies (IE), average binding energies, dissociation energies, and frontier orbital energy gaps are evaluated. The computed EAs and IEs are generally in fair agreement with experiment. The Nb15 system is observed to be stable and it can form a highly symmetric structure in all charged states with both open and closed electron shells.
Journal of Chemical Theory and Computation 2011 Volume 7(Issue 4) pp:1119-1130
Publication Date(Web):March 7, 2011
DOI:10.1021/ct1006482
Investigations on germanium clusters in the neutral, anionic, and dianion states Genx (n = 2−12 and x = 0, −1, −2) are performed using quantum chemical calculations with the B3LYP functional and the coupled-cluster singles and doubles [CCSD(T)] methods, in conjunction with the 6-311+G(d) basis set. An improved stochastic method is implemented for searching the low-lying isomers of clusters. Comparison of our results with previous reports on germanium clusters shows the efficiency of the search method. The Ge8 system is presented in detail. The anionic clusters Gen−/2− are studied theoretically and systematically for the first time, and their energetics are in good agreement with available experiments. The clusters Ge10, Ge102−, and Ge122− are, in their ground state, characterized by large highest occupied molecular orbital−lowest unoccupied molecular orbital gaps, high vertical and adiabatic detachment energies, and substantial average binding energies. The enhanced stability of these magic clusters can consistently be rationalized using the jellium electron shell model and the spherical aromatic character.
Electronic structure calculations suggest that hydrazine bisalane (AlH3NH2NH2AlH3, alhyzal) is a promising compound for chemical hydrogen storage (CHS). Calculations are carried out using the coupled-cluster theory CCSD(T) with the aug-cc-pVTZ basis set. Potential energy surfaces are constructed to probe the formation of, and hydrogen release from, hydrazine bisalane which is initially formed from the reaction of hydrazine with dialane. Molecular and electronic characteristics of both gauche and transalhyzal are determined for the first time. The gauche hydrazine bisalane is formed from starting reactants hydrazine + dialane following a movement of an AlH3 group from AlH3AlH3NH2NH2 rather than by a direct attachment of a separate AlH3 group, generated by predissociation of dialane, to AlH3NH2NH2. The energy barriers for dehydrogenation processes from gauche and transalhyzal are in the range of 21–28 kcal mol−1, which are substantially smaller than those of ca. 40 kcal mol−1 previously determined for the isovalent hydrazine bisborane (bhyzb) system. H2 release from hydrazine bisalane is thus more favored over that from hydrazine bisborane, making the Al derivative an alternative candidate for CHS.
Co-reporter:Vinh Son Nguyen, Saartje Swinnen, Jerzy Leszczynski and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2011 vol. 13(Issue 14) pp:6649-6656
Publication Date(Web):07 Mar 2011
DOI:10.1039/C0CP02484H
The reactivity of hydrazine in the presence of diborane has been investigated using ab initio quantum chemical computations (MP2 and CCSD(T) methods with the aug-cc-pVTZ basis set). Portions of the relevant potential energy surface were constructed to probe the formation mechanism of the hydrazine diborane (BH3BH3NH2NH2) and hydrazine bisborane (BH3NH2NH2BH3). The differences between both adducts are established. The release of hydrogen molecules from hydrazine bisborane adducts has also been characterized. Our results suggest that the BH3NH2NH2BH3 adduct, which has been prepared experimentally, is formed from the starting reactants hydrazine + diborane. The observed adduct is produced by a transfer of a BH3 group from BH3BH3NH2NH2 rather than by the direct attachment of a separate BH3 group, generated by predissociation of diborane, to BH3NH2NH2.
Co-reporter:Nguyen Tien Trung, Nguyen Phi Hung, Tran Thanh Hue and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2011 vol. 13(Issue 31) pp:14033-14042
Publication Date(Web):16 Jun 2011
DOI:10.1039/C1CP20533A
In this study, 16 gas phase complexes of the pairs of XCHZ and CO2 (X = F, Cl, Br; Z = O, S) have been identified. Interaction energies calculated at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ level including both BSSE and ZPE corrections range from −5.6 to −10.5 kJ mol−1 for XCHO⋯CO2 and from −5.7 to −9.1 kJ mol−1 for XCHS⋯CO2. Substitution of one H atom by one halogen in formaldehyde and thioformaldehyde reduces the interaction energy of XCHZ⋯CO2, while a CH3 substitution increases the interaction energy of both CH3CHO⋯CO2 and CH3CHS⋯CO2. NBO and AIM analyses also point out that the strength of Lewis acid–base interactions decreases going from >C1S3⋯C6 to >C1O3⋯C6 and to >C1–X4⋯C6. This result suggests the higher capacity of solubility of thiocarbonyl compounds in scCO2, providing an enormous potential application for designing CO2-philic materials based on the >CS functional group in competition with >CO. The Lewis acid–base interaction of the types >CS⋯C, >C–Cl⋯C and >C–Br⋯C is demonstrated for the first time. The contribution of the hydrogen bonding interaction to the total interaction energy is larger for XCHS⋯CO2 than for XCHO⋯CO2. Upon complexation, a contraction of the C1–H2 bond length and a blue shift of its stretching frequency have been observed, as compared to the isolated monomer, indicating the existence of a blue-shifting hydrogen bond in all complexes examined. Calculated results also lend further support for the viewpoint that when acting as proton donor, a C–H bond having a weaker polarization will induce a stronger distance contraction and frequency blue shift upon complexation, and vice versa.
Physical Chemistry Chemical Physics 2011 vol. 13(Issue 36) pp:16254-16264
Publication Date(Web):11 Aug 2011
DOI:10.1039/C1CP22078K
A systematic quantum chemical investigation on the electronic, geometric and energetic properties of AunV clusters with n = 1–14 in both neutral and anionic states is performed using BP86/cc-pVTZ-PP calculations. Most clusters having an even number of electrons prefer a high spin state. For odd-electron systems, a quartet state is consistently favoured as the ground state up to Au8V. The larger sized Au10V, Au12V and Au14V prefer a doublet state. The clusters prefer 2D geometries up to Au8V involving a weak charge transfer. The larger systems bear 3D conformations with a more effective electron transfer from Au to V. The lowest-energy structure of a size AunV is built upon the most stable form of Aun−1V. During the growth, V is endohedrally doped in order to maximize its coordination numbers and augment the charge transfer. Energetic properties, including the binding energies, embedding energies and second-order energy differences, show that the presence of a V atom enhances considerably the thermodynamic stability of odd-numbered gold clusters but reduces that of even-numbered systems. The atomic shape has an apparently more important effect on the clusters stability than the electronic structure. Especially, if both atomic shape and electronic condition are satisfied, the resulting cluster becomes particularly stable such as the anion Au12V−, which can thus combine with the cation Au+ to form a superatomic molecule of the type [Au12V]Au. Numerous lower-lying electronic states of these clusters are very close in energy, in such a way that DFT computations cannot clearly establish their ground electronic states. Calculated results demonstrate the existence of structural isomers with comparable energy content for several species including Au9V, Au10V, Au13V and Au14V.
Co-reporter:Saartje Swinnen, Rehab M.I. Elsamra, Vinh Son Nguyen, Jozef Peeters, Shaun A. Carl, Minh Tho Nguyen
Chemical Physics Letters 2011 Volume 513(4–6) pp:201-207
Publication Date(Web):15 September 2011
DOI:10.1016/j.cplett.2011.07.098
Abstract
Absolute rate constants for the gas-phase reaction C2H + SO2 are experimentally determined over the temperature range 295–800 K. C2H radicals are generated by pulsed 193-nm photolysis of C2H2 in the presence of SO2 and He or N2 buffer gas. The temperature dependence of the rate constants is established as kSO2 (T) = (0.86 ± 0.08) T3.80±0.13 exp[−(1222 ± 79)/T] cm3 s−1. The rate constants are moderately high at these temperatures and show negative temperature dependence. CCSD(T)/6-311++G(3df,2p) calculations of the potential energy surface confirm experimental findings and show that the reaction products are HCCO + SO.
Co-reporter:Saartje Swinnen, Vinh Son Nguyen, Minh Tho Nguyen
Chemical Physics Letters 2011 Volume 513(4–6) pp:195-200
Publication Date(Web):15 September 2011
DOI:10.1016/j.cplett.2011.07.097
Abstract
A ruthenium complex bearing cooperative PNP ligands showed unprecedented activities in homogenous catalysis of the reversible dehydrogenation of ammonia borane at both N-functionality and ethylene backbone. Quantum chemical calculations on a simple model of Ru-complex are carried out to probe the H2 release mechanism. The energy barrier for H2 formation from NH3BH3 is found to be strongly reduced with the presence of Ru-catalyst, and one additional H2 molecule can be produced. The Ru-complex can also be used as a catalyst for H2 release from other potential materials for chemical hydrogen storage such as hydrazine and ammonia alane.
Co-reporter:Saartje Swinnen, Vinh Son Nguyen, Minh Tho Nguyen
Chemical Physics Letters 2011 Volume 517(1–3) pp:22-28
Publication Date(Web):28 November 2011
DOI:10.1016/j.cplett.2011.10.010
Abstract
Ammonia borane and lithium amidoborane proved to have potential as efficient hydrogen storage materials. Although a combination of both compounds, LiNH2BH3–NH3BH3 has been shown to have even better hydrogen release characteristics, the inherent process was not clear. We determined the reaction mechanism of H2 release from LiNH2BH3–NH3BH3 using computational quantum chemistry methods (MP2 and CCSD(T) methods and aug-cc-pVnZ basis sets). The energy barrier for H2 release becomes lower than in pure LiNH2BH3 or NH3BH3 in agreement with experiment. Relevant transition structures are analyzed in detail and the initial formation of a cyclic adduct is likely responsible for the strong reduction of energy barriers.
Co-reporter:Truong Ba Tai, Hue Minh Thi Nguyen, Minh Tho Nguyen
Chemical Physics Letters 2011 Volume 502(4–6) pp:187-193
Publication Date(Web):27 January 2011
DOI:10.1016/j.cplett.2010.12.021
The coinage metal encapsulated group 14 cationic clusters X10M+ (X = Ge, Sn, Pb and M = Cu, Ag, Au) are investigated using quantum chemical calculations. While Cu- and Au-doped clusters have similar stability, Ag-doped counterparts are less stable. The D4d global minima have large frontier orbital gaps and binding energies, and are magic clusters of 40 valence electrons that satisfy the jellium shell model. The concept of doubly spherical aromaticity, based on the number of 2(N + 1)2 π and σ valence electrons, is proposed to account for the enhanced stability. Predictions are in good agreement with experiment for available X10Cu+ clusters.Graphical abstractResearch highlights► The coinage metal encapsulated group 14 cationic clusters X10M+ with X = Ge, Sn, Pb and M = Cu, Ag, Au are magic cluster with high stability and aromaticity. ► The concept of doubly spherical aromaticity, based on the number of 2(N + 1)2 π and σ valence electrons, is proposed to account for the enhanced stability. ► Predictions are in good agreement with experiment for available X10Cu+ clusters.
Co-reporter:Truong Ba Tai, Nguyen Minh Tam, Minh Tho Nguyen
Chemical Physics 2011 Volume 388(1–3) pp:1-8
Publication Date(Web):22 September 2011
DOI:10.1016/j.chemphys.2011.06.041
Abstract
Small zinc-doped tin clusters SnnZn, n = 1–12, are studied using DFT and CCSD(T) methods. The isomers are located using a stochastic search algorithm. The growth mechanism can be formulated as follows: (i) small clusters SnnZn with n ⩽ 8 are formed by capping Zn on a surface of Snn, (ii) a competition between exohedral and endohedral structures appears at n = 9 and 10, from which the endohedral structures become predominant, (iii) for n = 11 and 12, the clusters are formed by encapsulating Zn into the empty cages Snn. Icosahedral Sn12Zn (Ih) and Sn10Zn (D4d) are magic clusters with large HOMO–LUMO gaps, high binding energies and embedding energies. While both Sn12Zn and Sn10Zn clusters can be considered to be spherically aromatic with 8 valence π-electrons that satisfy the electron count rule of 2(N + 1)2, the enhanced stability of Sn12Zn (Ih) can further be rationalized in terms of its closed crystal field splitting shell.
The Journal of Physical Chemistry A 2011 Volume 115(Issue 35) pp:9993-9999
Publication Date(Web):July 12, 2011
DOI:10.1021/jp111324n
The group 14 clusters encapsulated by coinage metals in neutral and anionic states X10M0/– (X = Ge, Sn, Pb and M = Cu, Ag, Au) are investigated using quantum chemical calculations with the DFT/B3LYP functional and coupled-cluster CCSD(T) theory. Addition of transition metals into the empty cages forms high symmetry endohedral structures, except for Ge10Ag0/–. In agreement with experiments available for X10Cu, the D4d global minima of the anions are calculated to be magic clusters with large frontier orbital gaps, high vertical and adiabatic detachment energies, and large embedding energies and binding energies as compared to those of the empty cages X102–. The enhanced stability of these magic clusters can be rationalized by the three-dimensional aromaticity.
Co-reporter:Truong Ba Tai, Pham Vu Nhat, Minh Tho Nguyen, Shenggang Li, and David A. Dixon
The Journal of Physical Chemistry A 2011 Volume 115(Issue 26) pp:7673-7686
Publication Date(Web):June 14, 2011
DOI:10.1021/jp200992u
The stability, electronic structure, and thermochemical properties of the pure Lin and boron-doped LinB (n = 1–8) clusters in both neutral and cationic states are studied using electronic structure methods. The global equilibrium structures are established, and their heats of formation are evaluated using the G3B3 and CCSD(T)/CBS methods based on the density functional theory geometries. Theoretical adiabatic ionization energies (IEa) for the Lin clusters are in good agreement with experiment: Li2 (G3B3, 5.21 eV; CCSD(T), 5.14 eV; expt, 5.1127 ± 0.0003 eV), Li3 (4.16, 4.11, 4.08 ± 0.10), Li4 (4.76, 4.68, 4.70 ± 0.05), Li5 (4.11, 4.06, 4.02 ± 0.10), Li6 (4.46, 4.32, 4.20 ± 0.10), Li7 (4.07, 3.99, 3.94 ± 0.10), and Li8 (4.49, 4.31, 4.16 ± 0.10). The Li4 experimental IEa has been revised on the basis of the Franck–Condon simulations. Species Li5B, Li6B+, Li7B, and Li8B+ exhibit high stability as compared to their neighbors, which can be understood by considering the magic numbers of the phenomenological shell model (PSM).
Co-reporter:Pham Vu Nhat, Vu Thi Ngan, Truong Ba Tai, and Minh Tho Nguyen
The Journal of Physical Chemistry A 2011 Volume 115(Issue 15) pp:3523-3535
Publication Date(Web):March 29, 2011
DOI:10.1021/jp110758p
Geometric and electronic structures, vibrational properties, and relative stabilities of niobium clusters Nbn, n = 7−12, are studied using both DFT (BPW91 and M06 functionals) and CCSD(T) calculations with the cc-pVnZ-PP basis set. In each cluster, various lower-lying states are very close in energy in such a way that the ground state cannot be unambiguously established by DFT computations. Nb clusters tend to prefer the lowest possible spin state as the ground state, except for Nb12 (3Ag). The optimal structure of the cluster at a certain size does not simply grow from that of the smaller one by adding an atom randomly. Instead, the Nb clusters prefer a close-packed growth behavior. Nb10 has a spherically aromatic character, high chemical hardness and large HOMO−LUMO gap. Electron affinities, ionization energies, binding energy per atom, and the stepwise dissociation energies are evaluated. Energetic properties exhibit odd−even oscillations. Comparison with experimental values shows that both BPW91 and M06 functionals are reliable in predicting the EA and IE values, but the BPW91 is deficient in predicting the binding and dissociation energies. We re-examine in particular the experimental far IR spectra previously recorded using the IR-MPD and free electron laser spectrometric techniques and propose novel assignments for Nb7 and Nb9 systems. The IR spectra of the anions are also predicted.
We perform a systematic investigation on small silicon-doped boron clusters BnSi (n=1–7) in both neutral and anionic states using density functional (DFT) and coupled-cluster (CCSD(T)) theories. The global minima of these BnSi0/− clusters are characterized together with their growth mechanisms. The planar structures are dominant for small BnSi clusters with n≤5. The B6Si molecule represents a geometrical transition with a quasi-planar geometry, and the first 3D global minimum is found for the B7Si cluster. The small neutral BnSi clusters can be formed by substituting the single boron atom of Bn+1 by silicon. The Si atom prefers the external position of the skeleton and tends to form bonds with its two neighboring B atoms. The larger B7Si cluster is constructed by doping Si-atoms on the symmetry axis of the Bn host, which leads to the bonding of the silicon to the ring boron atoms through a number of hyper-coordination. Calculations of the thermochemical properties of BnSi0/− clusters, such as binding energies (BE), heats of formation at 0 K (ΔHf0) and 298 K (ΔHf[298]), adiabatic (ADE) and vertical (VDE) detachment energies, and dissociation energies (De), are performed using the high accuracy G4 and complete basis-set extrapolation (CCSD(T)/CBS) approaches. The differences of heats of formation (at 0 K) between the G4 and CBS approaches for the BnSi clusters vary in the range of 0.0–4.6 kcal mol−1. The largest difference between two approaches for ADE values is 0.15 eV. Our theoretical predictions also indicate that the species B2Si, B4Si, B3Si− and B7Si− are systems with enhanced stability, exhibiting each a double (σ and π) aromaticity. B5Si− and B6Si are doubly antiaromatic (σ and π) with lower stability.
Co-reporter:Vinh Son Nguyen, Heather L. Abbott, M. Michele Dawley, Thomas M. Orlando, Jerzy Leszczynski, and Minh Tho Nguyen
The Journal of Physical Chemistry A 2011 Volume 115(Issue 5) pp:841-851
Publication Date(Web):January 13, 2011
DOI:10.1021/jp109143j
The chemical transformations of formamide (NH2CHO), a molecule of prebiotic interest as a precursor for biomolecules, are investigated using methods of electronic structure computations and Rice−Rampserger−Kassel−Marcus (RRKM) theory. Specifically, quantum chemical calculations applying the coupled-cluster theory CCSD(T), whose energies are extrapolated to the complete basis set limit (CBS), are carried out to construct the [CH3NO] potential energy surface. RRKM theory is then used to systematically examine decomposition channels leading to the formation of small molecules including CO, NH3, H2O, HCN, HNC, H2, HNCO, and HOCN. The energy barriers for the decarboxylation, dehydrogenation, and dehydration processes are found to be in the range of 73−78 kcal/mol. H2 loss is predicted to be a one-step process although a two-step process is competitive. CO elimination is found to prefer a two-step pathway involving the carbene isomer NH2CHO (aminohydroxymethylene) as an intermediate. This CO-elimination channel is also favored over the one-step H2 loss, in agreement with experiment. The H2O loss is a multistep process passing through a formimic acid conformer, which subsequently undergoes a rate-limiting dehydration. The dehydration appears to be particularly favored in the low-temperature regime. The new feature identifies aminohydroxymethylene as a transient but crucial intermediate in the decarboxylation of formamide.
Co-reporter:Vu Thi Ngan ; Philipp Gruene ; Pieterjan Claes ; Ewald Janssens ; André Fielicke ; Minh Tho Nguyen ;Peter Lievens
Journal of the American Chemical Society 2010 Volume 132(Issue 44) pp:15589-15602
Publication Date(Web):October 20, 2010
DOI:10.1021/ja105099u
The growth mechanisms of small cationic silicon clusters containing up to 11 Si atoms, exohedrally doped by V and Cu atoms, are described. We find that as dopants, V and Cu follow two different paths: while V prefers substitution of a silicon atom in a highly coordinated position of the cationic bare silicon clusters, Cu favors adsorption to the neutral or cationic bare clusters in a lower coordination site. The different behavior of the two transition metals becomes evident in the structures of SinM+ (n = 4−11 for M = V, and n = 6−11 for M = Cu), which are investigated by density functional theory and, for several sizes, confirmed by comparison with their experimental vibrational spectra. The spectra are measured on the corresponding SinM+·Ar complexes, which can be formed for the exohedrally doped silicon clusters. The comparison between experimental and calculated spectra indicates that the BP86 functional is suitable to predict far-infrared spectra of these clusters. In most cases, the calculated infrared spectrum of the lowest-lying isomer fits well with the experiment, even when various isomers and different electronic states are close in energy. However, in a few cases, namely Si9Cu+, Si11Cu+, and Si10V+, the experimentally verified isomers are not the lowest in energy according to the density functional theory calculations, but their structures still follow the described growth mechanism. The different growth patterns of the two series of doped Si clusters reflect the role of the transition metal’s 3d orbitals in the binding of the dopant atoms.
Co-reporter:Truong Ba Tai, Pham Vu Nhat and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2010 vol. 12(Issue 37) pp:11477-11486
Publication Date(Web):02 Aug 2010
DOI:10.1039/C003401K
Quantum chemical calculations are performed on the aluminium doped lithium clusters LinAl at both neutral and cationic states using the DFT/B3LYP and CCSD(T) methods in conjugation with the aug-cc-pVaZ (a = D,T,Q) basis sets. The global minima are located and the growth mechanism is established. The electronic structure, geometrical parameters and energetic properties, such as average binding energy Eb, second difference of energy Δ2E, adiabatic and vertical ionization energy, and dissociated enthalpy, are evaluated using the coupled-cluster CCSD(T) method, whose energies are extrapolated to the complete basis set limit (CBS). The high stability of Li5Al, Li7Al, Li6Al+ and Li8Al+ that have the “magic numbers” of valence electrons, can be understood using the phenomenological shell model.
Co-reporter:Ling Lin, Pieterjan Claes, Tibor Höltzl, Ewald Janssens, Torsten Wende, Risshu Bergmann, Gabriele Santambrogio, Gerard Meijer, Knut R. Asmis, Minh Tho Nguyen and Peter Lievens
Physical Chemistry Chemical Physics 2010 vol. 12(Issue 42) pp:13907-13913
Publication Date(Web):20 Sep 2010
DOI:10.1039/C0CP00911C
The geometric and electronic structure of the Au6Y+ cation is studied by gas phase vibrational spectroscopy combined with density functional theory calculations. The infrared photodissociation spectrum of Au6Y+·Ne is measured in the 95–225 cm−1 energy range and exhibits two characteristic absorption bands at 181 cm−1 and 121 cm−1. Based on DFT/BP86 quantum chemical calculations, the infrared spectrum is assigned to the lowest energy species found, an eclipsed C3v geometry. The 3D structure of Au6Y+ is considerably different from those previously found for both the neutral Au6Y (quasi-planar circular geometry) and the anionic Au6Y− (planar D6h symmetry). The different geometries are related to different electronic structures in agreement with 2D and 3D phenomenological shell models for metal clusters.
Co-reporter:Saartje Swinnen, Vinh Son Nguyen, Minh Tho Nguyen
Chemical Physics Letters 2010 Volume 496(1–3) pp:25-31
Publication Date(Web):20 August 2010
DOI:10.1016/j.cplett.2010.07.024
Abstract
The catalytic effect of LiH, 2LiH, BeH2 and 2BeH2 on the hydrogen release from hydrazine is studied using quantum chemical calculations. Potential energy surfaces of the relevant reactions are constructed using MP2 and CCSD(T) methods and the aug-cc-pVnZ basis functions and extrapolated to the complete basis sets. No extra energy is needed for hydrogen release, and regeneration of the catalyst in the case of 2LiH. A small amount of energy is needed for the regeneration of 2BeH2. The results suggest the capacities of small hydride molecules to act as efficient catalysts for H2 production.
Co-reporter:Tibor Höltzl, Tamás Veszprémi, Minh Tho Nguyen
Chemical Physics Letters 2010 Volume 499(1–3) pp:26-30
Publication Date(Web):20 October 2010
DOI:10.1016/j.cplett.2010.08.085
Abstract
Substitution of C–H groups in 1,3,5-trimethylenebenzene (1,3,5-TMB) triradical by one, two or three N, P and As atoms is studied using CASSCF/CASPT2 computations with a ANO-RCC triple-ζ basis set. The computations showed that consistently with the previous results, the unpaired electrons in 1,3,5-TMB are highly delocalized, but our computations show that they become localized on the P and As atoms, with a more pronounced localization effect of As. All studied compounds have a quartet ground state, but the results showed in this Letter indicate that the quartet–doublet energy gap is reduced upon substitution. Localization of the unpaired electrons on the ring reduces the average distances between them, which induces a destabilization of the quartet states compared to the doublet. Therefore the position of the unpaired electrons and the doublet–quartet gap can be tuned by chemical methods, which opens the way to fine-tune the properties of the building blocks of organic magnets.
Co-reporter:Ling Lin, Peter Lievens, Minh Tho Nguyen
Chemical Physics Letters 2010 Volume 498(4–6) pp:296-301
Publication Date(Web):8 October 2010
DOI:10.1016/j.cplett.2010.08.068
Abstract
CO adsorption on AunY (n = 1–9) is investigated using density functional theory, and compared with pure gold clusters Aun+1. CO is confirmed to prefer on-top binding to the least coordinated Au atom for the pure gold clusters, whereas for Y-doped gold clusters, the most favorable sites are dependent on the cluster size; the Au on-top sites are favored from n = 6. CO behaves as an electron acceptor, and the C–O vibrational frequency of the metal–CO complexes correlates with the degree of charge transfer to CO. Reactivity with CO is decreased after one gold atom is replaced by one yttrium for most sizes.
Chemical Physics Letters 2010 Volume 498(1–3) pp:120-124
Publication Date(Web):30 September 2010
DOI:10.1016/j.cplett.2010.08.030
CO oxidation on Au2n- (n = 1 – 4) cluster anions are investigated using DFT methods. O2 prefers to bind on-top to a gold atom. Binding energies of O2 on Au2n- are in the range of 0.8–1.2 eV. CO molecule can then insert into the O–O bond of Au2nO2- to form carbonate species, that are the most stable Au2nCO3- isomers, and thus expected to play an important role in CO oxidation. The negative charges of the gold carbonates are mainly distributed on the CO3 moieties. The gold carbonate anions can yield CO2 by direct dissociation or react with another CO to produce 2(CO2).The gold carbonates are found to be the most stable structures among the possible Au2nCO3- isomers; they are thus expected to play an important role in the CO oxidation on Au2n- following an Eley–Rideal mechanism.
Co-reporter:Saartje Swinnen, Vinh Son Nguyen, Minh Tho Nguyen
Chemical Physics Letters 2010 Volume 489(4–6) pp:148-153
Publication Date(Web):9 April 2010
DOI:10.1016/j.cplett.2010.02.060
Abstract
Lithium amidoboranes show potential as hydrogen storage materials. The mechanism for production of LiNH2BH3 from LiH and NH3BH3 is studied using quantum chemical methods (MP2. CCSD(T)). The reaction is exothermic with a low energy barrier. Pathways of H2 production from LiNH2BH3, NaNH2NH3 and (LiNH2BH3)2 show that a H-transfer from B to Li or Na constitutes the rate-determining step. Substitution of one borane hydrogen by a more electron-donating group (OH and NH2) tends to lower the energy barriers for H2 release. Calculated results point out that LiNH2BH2OH and LiNH2BH2NH2 can be considered as potentially promising materials for hydrogen storage.
Chemical Physics Letters 2010 Volume 492(4–6) pp:290-296
Publication Date(Web):7 June 2010
DOI:10.1016/j.cplett.2010.04.072
Abstract
Structure and stability of the Ge12Mx clusters with M = Li, Na, Be, Mg, B, Al, and x from −1 to +1, each containing 50 valence electrons are investigated using DFT calculations. The global minima turn to be the high symmetry icosahedral structures, with large HOMO–LUMO gaps and high detachment or ionization energies. In particular, Li is found for the first time to be located at the center of an icosahedron Ge12Li−. The high thermodynamic stability of the icosahedra arises from a combination of their closed crystal field shells, spherical aromaticity and electrostatic attraction force.
Chemical Physics 2010 Volume 375(Issue 1) pp:35-45
Publication Date(Web):14 September 2010
DOI:10.1016/j.chemphys.2010.07.015
Abstract
The electronic structure and molecular properties of a series of small mixed lithium boron clusters BnLi (n = 1–8) in both neutral and anionic states are investigated using quantum chemical methods. The lowest-energy structures are identified. Standard heats of formation, adiabatic electron affinities (EAs) and vertical detachment energies (VDEs) are predicted using coupled-cluster theory CCSD(T)/6-311+G(d) calculations. Addition of Li to Bn to form BnLi clusters marginally distorts the geometries of the parent Bn. The chemical bond in BnLi has a highly ionic character in which the positive charges are located on Li. MO analysis confirms a similar degree of aromaticity between BnLi and . The relative stabilities of the clusters are evaluated using the average binding energy (Eb), second order difference in energy (Δ2E), resonance energy (RE) and normalized resonance energy (NRE). Both REs and NREs quantities of the Li-doped boron clusters are consistently larger than those of the pure Bn clusters. The B3Li species exhibits a remarkably high stability within the BnLi series.
Co-reporter:Ngo Tuan Cuong, Truong Ba Tai, Vu Thi Thu Ha, Minh Tho Nguyen
The Journal of Chemical Thermodynamics 2010 Volume 42(Issue 4) pp:437-440
Publication Date(Web):April 2010
DOI:10.1016/j.jct.2009.10.006
Thermochemical parameters of caffeine 1, theophylline 2, xanthine 3, uracil, and imidazole derivatives are determined by quantum chemical calculations. Using the composite G3B3 method, the standard heat of formation of caffeine in the gaseous phase amounts to ΔfHg∘(1)=-243±8kJ·mol-1, which lends a support for the recent experimental value of −237.0 ± 2.5 kcal · mol−1. We also obtain ΔfHg∘(2)=-232±8kJ·mol-1andΔfHg∘(3)=-209±8kJ·mol-1. The adiabatic ionization energies are IEa(1) = 7.9 eV, IEa(2) = 8.1 eV, and IEa(3) = 8.5 eV using B3LYP calculations. The enhanced ability of caffeine to eject electron, as compared to the parent compounds and cyclic components, is of interest with regard to its potential use as a corrosion inhibitor.
Co-reporter:Ling Lin, Peter Lievens, Minh Tho Nguyen
Journal of Molecular Structure: THEOCHEM 2010 Volume 943(1–3) pp:23-31
Publication Date(Web):15 March 2010
DOI:10.1016/j.theochem.2009.09.004
A theoretical search for aromatic seven-membered rings has been carried out using density functional theory calculations with the B3LYP functional. The rings considered include the C7H7+ cation and its hetero-derivatives by replacing C, CH or CC units by B, Al, Ga, Si, Ge, N, P, As and BN group, the B82− dianion, the CB7− anion, the neutral S3N4 ring and its derivatives by substituting one or two N atoms by CH group, the S4N3+ cation and the all-metallic cycles M73− and M7T with M = Cu, Ag, Au and T = Y, Sc. Most molecules studied have planar structure in their electronic ground state. Vibrational spectra of some derivatives are plotted. Nucleus independent chemical shift (NICS) indices show that the parent molecule C7H7+ has a higher degree of aromaticity than its derivatives. Substitution of N in S3N4 by CH marginally influences the aromaticity in such a way that S3N4, S3N3(CH) and S3N2(CH)2 are similarly aromatic. The B82− dianion and both D7h and C2v isomers of CB7− possess comparable aromatic character. As for the all-metal clusters Cu73−, Ag73−, Au73−, Cu7Sc, Cu7Y, Ag7Y and Au7Y, the binary clusters become more aromatic than the pure metal anions, and are thus characterized by σ-aromaticity.
The geometric, spectroscopic, and electronic properties of neutral yttrium-doped gold clusters AunY (n=1–9) are studied by far-infrared multiple photon dissociation (FIR-MPD) spectroscopy and quantum chemical calculations. Comparison of the observed and calculated vibrational spectra allows the structures of the isomers present in the molecular beam to be determined. Most of the isomers for which the IR spectra agree best with experiment are calculated to be the energetically most stable ones. Attachment of xenon to the AunY cluster can cause changes in the IR spectra, which involve band shifts and band splittings. In some cases symmetry changes, as a result of the attachment of xenon atoms, were also observed. All the AunY clusters considered prefer a low spin state. In contrast to pure gold clusters, which exhibit exclusively planar lowest-energy structures for small sizes, several of the studied species are three-dimensional. This is particularly the case for Au4Y and Au9Y, while for some other sizes (n=5, 8) the 3D structures have an energy similar to that of their 2D counterparts. Several of the lowest-energy structures are quasi-2D, that is, slightly distorted from planar shapes. For all the studied species the Y atom prefers high coordination, which is different from other metal dopants in gold clusters.
Co-reporter:Pham Vu Nhat, Vu Thi Ngan and Minh Tho Nguyen
The Journal of Physical Chemistry C 2010 Volume 114(Issue 31) pp:13210-13218
Publication Date(Web):July 16, 2010
DOI:10.1021/jp103484k
Geometries, electronic structures, and vibrational spectra of a series of small niobium clusters Nbn (n = 2−6) and their cations and anions are reinvestigated with DFT calculations with cc-pVaZ-PP basis sets. CCSD(T) calculations are also carried out for relative energies. In each cluster, different lower lying states are close in energy or quasidegenerate. Stable Nbn clusters prefer high coordination state and 3D shape. Clusters with an odd number of electrons usually have a doublet ground state (except for Nb2+). Neutral and cationic clusters with an even number of electrons tend to possess a triplet ground state, with an exception for Nb4 (Td, 1A1). For anions, a competition between low and high spin manifolds is observed. Due to the closed electronic shells, the systems possessing 10 (Nb) and 20 valence electrons (Nb4) are observed to be more thermodynamically stable than their neighbors. Electron affinities and ionization energies are reevaluated. In particular, new assignments for the vibrational spectra of the Nb5 and Nb6 systems are proposed on the basis of a comparison of calculated results with available experimental data.
Co-reporter:Truong Ba Tai, Minh Tho Nguyen and David A. Dixon
The Journal of Physical Chemistry A 2010 Volume 114(Issue 8) pp:2893-2912
Publication Date(Web):January 29, 2010
DOI:10.1021/jp909512m
The molecular and electronic structures of a series of small boron monoxide and dioxide clusters BnOm (n = 5−10, m = 1, 2) plus their anions were predicted. The enthalpies of formation (ΔHf’s), electron affinities (EAs), vertical detachment energies, and energies of different fragmentation processes are predicted using the G3B3 method. The G3B3 results were benchmarked with respect to more accurate CCSD(T)/CBS values for n = 1−4 with average deviations of ±1.5 kcal/mol for ΔHf’s and ±0.03 eV for EAs. The results extend previous observations on the growth mechanism for boron oxide clusters: (i) The low spin electronic state is consistently favored. (ii) The most stable structure of a neutral boron monoxide BnO is obtained either by condensing O on a BB edge of a Bn cycle, or by binding one BO group to a Bn−1 ring. The balance between both factors is dependent on the inherent stability of the boron cycles. (iii) A boron dioxide is formed by incorporating the second O atom into the corresponding monoxide to form BO bonds. (iv) A BnOm− anion is constructed with BO groups bound to the Bn−1 or Bn−2 rings (yielding the Bn−2(BO)2− species). This becomes the preferred geometry for the larger boron dioxides, even in the neutral state. The boronyl group mainly behaves as an electron-withdrawing substituent reducing the binding energy and resonance energy of the oxides. (v) The boron oxides conserve some of the properties of the parent boron clusters such as the planarity and multiple aromaticity.
The Journal of Physical Chemistry A 2010 Volume 114(Issue 28) pp:7609-7615
Publication Date(Web):June 28, 2010
DOI:10.1021/jp103180y
The geometrical and electronic structures of the Si82− dianion and isovalent silicon clusters doped by main second-row elements including Li@Si8−, Be@Si8, B@Si8+, C@Si82+, N@Si83+, and O@Si84+, are investigated using quantum chemical methods. The analyses of phenomenological shell model (PSM) combined with partial electron localizability indicators (pELI-D) rationalize the existence of cubic silicon clusters. A cubic cluster can be formed, in the cases of Be@Si8, B@Si8+, and C@Si82+, when three conditions are satisfied, namely, a full occupancy of electronic shells (34 electrons), a presence of positive charge at the center, and a type of spherical aromaticity. A chemical bonding picture for the cubic cage of the doped silicon clusters is illustrated. Each Si atom has four lobes of sp3 hybridization in which three lobes make three covalent σ bonds with other Si atoms, and the fourth lobe makes a chemical bond with the dopant. The eight delocalized electrons distributed on the fourth lobes describing the bonding between dopant and Si cage follow the Hirsch rule. We demonstrate that a way of applying electron counting rule is to take into account delocalized electrons on the shell orbitals with N > 1 (2S and 2P shell orbitals).
Co-reporter:Vitaly G. Kiselev, Saartje Swinnen, Vinh Son Nguyen, Nina P. Gritsan and Minh Tho Nguyen
The Journal of Physical Chemistry A 2010 Volume 114(Issue 17) pp:5573-5579
Publication Date(Web):April 9, 2010
DOI:10.1021/jp911655a
Different uni- and bimolecular reactions of hydroxymethylene, an important intermediate in the photochemistry of formaldehyde, as well as its halogenated derivatives (XCOH, X = H, F, Cl, Br), have been considered using high-level CCSD(T)/CBS quantum chemical methods. The Wentzel−Kramers−Brillouin (WKB) and Eckart approximations were applied to estimate the tunneling rate constant of isomerization of trans-HCOH to H2CO, and the WKB procedure was found to perform better in this case. In agreement with recent calculations and experimental observations [Schreiner et al., Nature 2008, 453, 906], the half-life of HCOH at the low temperature limit in the absence of bimolecular processes was found to be very long (∼2.1 h). The corresponding half-life at room temperature was also noticeable (∼35 min). Bimolecular reactions of trans-hydroxymethylene with parent formaldehyde yield primarily more thermodynamically favorable glycolaldehyde via the specific mechanism involving 5-center transition state. The most preferable reaction of cis-hydroxymethylene with formaldehyde yields carbon monoxide and methanol. Due to very low activation barriers, both processes occur with nearly a collision rate. If the concentration of HCOH (and its halogenated analogues XCOH as well) is high enough, the bimolecular reactions of this species with itself become important, and H2CO (or X(H)CO) is then formed with a collision rate. The singlet−triplet energy separation of trans-HCOH is confirmed to be ∼−25 kcal/mol.
Co-reporter:Truong Ba Tai, Daniel J. Grant, Minh Tho Nguyen and David A. Dixon
The Journal of Physical Chemistry A 2010 Volume 114(Issue 2) pp:994-1007
Publication Date(Web):December 11, 2009
DOI:10.1021/jp9085848
Thermochemical parameters of a set of small-sized neutral (Bn) and anionic (Bn−) boron clusters, with n = 5−13, were determined using coupled-cluster theory CCSD(T) calculations with the aug-cc-pVnZ (n = D, T, and Q) basis sets extrapolated to the complete basis set limit (CBS) plus addition corrections and/or G3B3 calculations. Enthalpies of formation, adiabatic electron affinities (EA), vertical (VDE), and adiabatic (ADE) detachment energies were evaluated. Our calculated EAs are in good agreement with recent experiments (values in eV): B5 (CBS, 2.29; G3B3, 2.48; exptl., 2.33 ± 0.02), B6 (CBS, 2.59; G3B3, 3.23; exptl., 3.01 ± 0.04), B7 (CBS, 2.62; G3B3, 2.67; exptl., 2.55 ± 0.05), B8 (CBS, 3.02; G3B3, 3.11; exptl., 3.02 ± 0.02), B9 (G3B3, 3.03; exptl., 3.39 ± 0.06), B10 (G3B3, 2.85; exptl., 2.88 ± 0.09), B11 (G3B4, 3.48;, exptl., 3.43 ± 0.01), B12 (G3B3, 2.33; exptl., 2.21 ± 0.04), and B13 (G3B3, 3.62; exptl., 3.78 ± 0.02). The difference between the calculated adiabatic electron affinity and the adiabatic detachment energy for B6 is due to the fact that the geometry of the anion is not that of the ground-state neutral. The calculated adiabatic detachment energies to the 3Au, C2h and 1Ag, D2h excited states of B6, which have geometries similar to the 1Ag, D2h state of B6−, are 2.93 and 3.06 eV, in excellent agreement with experiment. The VDEs were also well reproduced by the calculations. Partitioning of the electron localization functions into π and σ components allows probing of the partial and local delocalization in global nonaromatic systems. The larger clusters appear to exhibit multiple aromaticity. The binding energies per atom vary in a parallel manner for both neutral and anionic series and approach the experimental value for the heat of atomization of B. The resonance energies and the normalized resonance energies are convenient indices to quantify the stabilization of a cluster of elements.
Co-reporter:Nguyen Tien Trung, Tran Thanh Hue and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2009 vol. 11(Issue 6) pp:926-933
Publication Date(Web):15 Dec 2008
DOI:10.1039/B816112G
The hydrogen-bonded complexes formed from interaction of trihalomethanes CHX3 (X = F, Cl, Br) with nitrosyl hydride HNO were studied using ab initio MO calculations at the second-order perturbation theory (MP2/6-311++G(d,p)). Each interaction contains at least five separate equilibrium structures. Calculated binding energies range from 4 to 8 kJ mol−1 with both ZPE and BSSE corrections. While CHBr3 leads to the most stable complexes with HNO, CHF3 forms the least stable counterparts. The strength of complexes thus tends to increase from F to Cl to Br, which is consistent with a decrease of deprotonation enthalpy of the corresponding C–H bonds. It is remarkable that all the C–H and N–H bonds are shortened upon complexation, giving rise to an increase of their stretching frequencies. A blue shift is thus observed for the N–H bonds of the type N–H⋯X (X = F, Cl, Br); such a contraction of the covalent N–H bond is extremely rare.
Co-reporter:Saartje Swinnen, Vinh Son Nguyen, Shogo Sakai, Minh Tho Nguyen
Chemical Physics Letters 2009 Volume 472(4–6) pp:175-180
Publication Date(Web):20 April 2009
DOI:10.1016/j.cplett.2009.02.078
Abstract
Producing H2 from water is a very challenging task. Using quantum chemical calculations with the MP2 and CCSD(T) methods and the aug-cc-pVnZ basis sets (extrapolated to CBS), we investigated the possibilities of B2H6, AlH3, Al2H6 and AlH3BH3 to act as catalysts in the reactions that split water. Hydrogen production from H2O is greatly accelerated in the presence of alane, dialane or borane–alane in such a way that reaction H2O + Al2H6 is a nearly spontaneous process. The main catalytic effect of AlH3 arises from the occurrence of a strong dihydrogen bond of the type Al–Hδ−–δ+H–O within a cyclic transition structure.
Co-reporter:Tibor Höltzl, Nele Veldeman, Tamás Veszprémi, Peter Lievens, Minh Tho Nguyen
Chemical Physics Letters 2009 Volume 469(4–6) pp:304-307
Publication Date(Web):17 February 2009
DOI:10.1016/j.cplett.2009.01.001
Abstract
The Cu6Sc+ cation was detected experimentally as a stable bimetallic cluster [N. Veldeman, T. Höltzl, S. Neukermans, T. Veszprémi, M.T. Nguyen, Phys. Rev. A. 76 (2007) 011201]. DFT calculations using the BP86 functional with various basis sets yield a D3h tricapped tetrahedron shaped structure as the most stable isomer. This cluster has a S2P6 electron configuration with a closed electronic shell, and its electrons tend to localize around Sc thereby fully occupying electron shell orbitals. Stability, high symmetry, and the NICS show that this cluster is aromatic. The neutral Cu5Sc cluster has the same electron configuration as Cu6Sc+. In agreement with qualitative predictions, DFT computations show that the C4v octahedron-shaped Cu5Sc cluster features a closed electronic structure and an aromatic character.
Chemical Physics Letters 2009 Volume 476(4–6) pp:236-239
Publication Date(Web):16 July 2009
DOI:10.1016/j.cplett.2009.06.008
Abstract
Electronic structure of Al2Na was reinvestigated using MP2, CCSD(T) and MRMP2/CASSCF(7,7) calculations with aug-cc-pVnZ (n = D, T and Q) basis sets. In contrast to a recent report, the triangular Al2Na has a 2B1 ground state with a 4B1–2B1 gap of 2.5 kcal/mol, rather than a 2A1 state (13.7 kcal/mol higher). In addition, Al2Na+ exhibits a linear 1Σ+ ground state with a singlet–triplet gap of −5.2 kcal/mol rather than a cyclic 1A1 state (6.0 kcal/mol above the linear). Al2Na− has a triangular 3B1 ground state and a singlet–triplet gap of 6.9 kcal/mol. We obtained IEa(Al2Na) = 5.3 ± 0.2 eV and EAa(Al2Na) = 1.3 ± 0.2 eV.
Chemical Physics Letters 2009 Volume 483(1–3) pp:35-42
Publication Date(Web):24 November 2009
DOI:10.1016/j.cplett.2009.10.054
The isomers of B6O3, B6O3- and B6O32- were studied using DFT, MP2 and G3B3 methods. The dianion B3(BO)32- possesses a perfect triangular D3h form, that undergoes to lower symmetry in the anionic (C2v) and neutral (Cs) states. The adiabatic electron affinity of B6O3 and vertical detachment energy of B6O3- amount to 3.97 and 4.07 eV, respectively. NBO and ELF analyses reveal that B6O32- possesses a global π-aromaticity, a local σ-electron delocalization, and a three-center-two-electron bond. The high stability of the dianion illustrates the growth mechanism of boron oxide clusters that is governed by formation of strong BO bonds and small boron cycles.The B6O32- dianion has a global π-aromaticity, local σ-electron delocalization and three-center-two-electron bond, and its structure illustrates the growth pattern of boron oxides governed by formation of BO bonds and small boron cycles.
Co-reporter:Vu Thi Ngan, Jorg De Haeck, Hai Thuy Le, G. Gopakumar, Peter Lievens and Minh Tho Nguyen
The Journal of Physical Chemistry A 2009 Volume 113(Issue 32) pp:9080-9091
Publication Date(Web):July 21, 2009
DOI:10.1021/jp9056913
We report a combined experimental and quantum chemical study of the small neutral and cationic germanium-doped lithium clusters LinGe0,+ (n = 1−7). The clusters were detected by time-of-flight mass spectrometry after laser vaporization and ionization. The molecular geometries and electronic structures of the clusters were investigated using quantum chemical calculations at the DFT/B3LYP and CCSD(T) levels with the aug-cc-pVnZ basis sets. While Li3Ge0,+ and Li4Ge+ prefer planar structures, the clusters from Li4Ge to Li7Ge and the corresponding cations (except Li4Ge+) exhibit nonplanar forms. Clusters having from 4 to 6 valence electrons prefer high spin structures, and low spin ground states are derived for the others because valence electron configurations are formed by filling the electron shells 1s/1p/2s/2p based on Pauli’s and Hund’s rules. Odd−even alternation is observed for both neutral and cationic clusters. Because of the closed electronic shells, the 8- and 10-electron systems are more stable than the others, and the 8-electron species (Li4Ge, Li5Ge+) are more favored than the 10-electron ones (Li6Ge, Li7Ge+). This behavior for Ge is different from C in their doped Li clusters, which can be attributed to the difference in atomic radii. The averaged binding energy plot for neutrals tends to increase slowly with the increasing number of Li atoms, while the same plot for cations shows a maximum at Li5Ge+, which is in good agreement with the mass spectrometry experiment. Atom-in-molecules (AIM) analysis suggests that Li atoms do not bond to one another but through Ge or pseudoatoms, and an essentially ionic character can be attributed to the cluster chemical bonds. An interesting finding is that the larger clusters have the smallest adiabatic ionization energies known so far (IEa ≈ 3.5 eV).
Co-reporter:Vinh Son Nguyen, Saartje Swinnen, Minh Tho Nguyen and David A. Dixon
The Journal of Physical Chemistry C 2009 Volume 113(Issue 43) pp:18914-18926
Publication Date(Web):October 5, 2009
DOI:10.1021/jp904344p
Electronic structure calculations at the CCSD(T) level with the aug-cc-pVnZ and aug-cc-pV(n+d)Z basis sets (n = D, T, and Q) were employed to construct the potential energy surfaces for H2 release from a series of derivatives of ammonia alane. H2 production from AlH3NH3 is facilitated by the addition of alane, which plays the role of a Lewis acid bifunctional catalyst. H2 release is not competitive from the combination of two AlH3NH3 monomers in the form of a stable dimer. An alternative route involves the reactions of one AlH3NH3 monomer with two separated molecules of AlH3 and NH3, in two successive steps: (i) an initial condensation of AlH3NH3 with NH3 leads to a stable linear NH3AlH3NH3 species, and (ii) a subsequent combination of the latter with AlH3 gives rise to a bicycle framed by dihydrogen bonds. From the new bicyclic form, H2 production becomes a facile process. Condensation of dialane (Al2H6) with two separated ammonia molecules (2NH3) is also a potential route for H2 release. The reaction enthalpies and free energies at 298 K were also evaluated for the pathways considered. Although the effects of temperature and entropy are important, the ΔG298 preferred pathways remain similar to the energy (ΔH(0 K)) profiles.
Co-reporter:Nguyen Tien Trung, Tran Thanh Hue and Minh Tho Nguyen
The Journal of Physical Chemistry A 2009 Volume 113(Issue 13) pp:3245-3253
Publication Date(Web):March 3, 2009
DOI:10.1021/jp810826z
Weak interactions of monosubstituted formaldehydes and thioformaldehydes with nitrosyl hydride were investigated by using ab initio MO calculations at the MP2/aug-cc-pVTZ level. Thirty two equilibrium structures having different complex forms were located on the corresponding potential energy surfaces (all having Cs symmetry). Obtained binding energies, which include both ZPE and BSSE corrections, range from 7 to 14 kJ·mol−1 and 6 to 12 kJ·mol−1 for complexes of substituted formaldehydes and thioformaldehydes, respectively. In each geometrical structure, the (XCHO,HNO) complex is consistently more stable than the (XCHS,HNO) complex. The H-bond strength significantly increases when one H atom is replaced by a methyl group in both formaldehyde and thioformaldehyde. When replacing H by a halogen atom, the binding energy tends to decrease. It is remarkable that all the C−H and N−H bonds are shortened upon complexation, resulting in an increase of their stretching frequencies. Furthermore, the blue shifts are consistently observed for the interacting N−H bonds in N−H···X, Z, with X = F, Cl, Br, and Z = O, S; such contraction of a covalent N−H bond is extremely rare. In addition, the N−H bond length contraction and its frequency blue shift in the N−H···S complex have been revealed for the first time.
Chemical Physics Letters 2008 Volume 456(4–6) pp:141-145
Publication Date(Web):5 May 2008
DOI:10.1016/j.cplett.2008.03.036
Abstract
Quantum chemical calculations (G3B3 and B3LYP/6-311++G(d,p)) and tandem mass spectrometric experiments demonstrate the higher stability of 1,3-cyclohexadienylidenemethanone radical cation compared to ionized benzaldehyde. Characterized by a heat of formation of 833 kJ mol−1 (at 298 K), this ketene C7H6O+ isomer is found 43 kJ mol−1 more stable. It has been generated by ion/molecule reaction between ionized benzaldehyde and neutral methanol, a new example of proton transport catalysis tautomerization. The greater stability of the ketene ion is due to the low IE of the neutral (7.49 eV).
Co-reporter:Tibor Höltzl;Ewald Janssens Dr.;Nele Veldeman Dr.;Tamás Veszprémi Dr.;Peter Lievens Dr. Dr.
ChemPhysChem 2008 Volume 9( Issue 6) pp:833-838
Publication Date(Web):
DOI:10.1002/cphc.200700752
Abstract
Density functional theory calculations demonstrate that the global minimum of the Cu7Sc potential energy surface is a seven-membered ring of copper atoms with scandium in its center, yielding a planar D7h structure. Nucleus-independent chemical shifts [NICS(1)zz and NICS(2)zz] show that this cluster has aromatic character, which is consistent with the number of 4s electrons of copper and scandium plus the 3d electrons of scandium satisfying Hückel’s rule. According to a canonical MO decomposition of NICS(1)zz and NICS(2)zz, the MOs consisting of the 4s atomic orbitals are mainly responsible for the aromatic behavior of the cluster. The electron localizability indicator (ELI–D) and its canonical MO decomposition (partial ELI–D) suggest that a localized basin is formed in Cu7Sc by the copper atoms whereas the two circular localized domains are situated below and above the ring. The planar Cu7Sc cluster can thus be considered as a σ-aromatic species. These findings agree with the phenomenological shell model.
Co-reporter:G. Gopakumar, Vu Thi Ngan, Peter Lievens and Minh Tho Nguyen
The Journal of Physical Chemistry A 2008 Volume 112(Issue 47) pp:12187-12195
Publication Date(Web):November 5, 2008
DOI:10.1021/jp805173n
Quantum chemical calculations were applied to investigate the electronic structure of germanium hydrides, GenH (n = 1, 2, 3), their cations, and anions. Computations using a multiconfigurational quasi-degenerate perturbation approach (MCQDPT2) based on complete active space wave functions (CASSCF), multireference perturbation theory (MRMP2), and density functional theory reveal that Ge2H has a 2B1 ground state with a doublet-quartet gap of ∼39 kcal/mol. A quasidegenerate 2A1 state has been derived to be 2 kcal/mol above the ground state (MCQDPT2/aug-cc-pVTZ). In the case of the cation Ge3H+ and anion Ge3H−, singlet low-lying electronic states are derived, that is, 1A′ and 1A1, respectively. The singlet−triplet energy gap is estimated to 6 kcal/mol for the cation. An “Atoms in Molecules” (AIM) analysis shows a certain positive charge on the Gen (n = 1, 2, 3) unit in its hydrides, in accordance with the NBO analysis. The topologies of the electron density of the germanium hydrides are different from that of the lithium-doped counterparts. On the basis of our electron localization function (ELF) analysis, the Ge−H bond in Ge2H is characterized as a three-center-two-electron bond. Some key thermochemical parameters of GenH have also been derived.
Chemical Physics Letters 2007 Volume 448(4–6) pp:183-188
Publication Date(Web):14 November 2007
DOI:10.1016/j.cplett.2007.10.033
CCSD(T)/CBS calculations plus corrections predict the heats of formation of formaldehyde O-oxide 1 and dioxirane 2: ΔHf0(CH2OO,1)=28.1 and 26.4 kcal/mol, and ΔHf0(CH2O2,2)=3.0 and 1.2 kcal/mol at 0 K and 298 K, respectively. The adiabatic ionization energies are IEa(1) = 9.98 eV and IEa(2) = 10.82 eV. Protonation of carbonyl oxide takes place at terminal oxygen with a proton affinity of PA(1) = 203.5 ± 1.0 kcal/mol. The vertical triplet state 3A′ of 1 is located at 1.82 eV above the ground state (with errors of ±1.0 kcal/mol or ±0.05 eV). The parent Criegee intermediate is unstable with respect to O–O bond cleavage but becomes more stable upon ionization.CCSD(T)/CBS calculations plus corrections predict the heats of formation of formaldehyde O-oxide 1 (Criegee intermediate) and dioxirane 2 to be ΔHf0(1)=26.4 kcal/mol and ΔHf0(2)=1.2 kcal/mol at 298 K.
Chemical Physics Letters 2007 Volume 449(1–3) pp:11-17
Publication Date(Web):26 November 2007
DOI:10.1016/j.cplett.2007.10.015
The electronic structure of indigo was determined using B3LYP functional and different basis sets. The first triplet 3Au (C2h) state of indigo is characterized by the singlet–triplet gaps of ΔET–S(vertical) = 1.23 eV and ΔET–S(adiabatic) = 0.95 eV (recent experiment: 1.04 ± 0.1 eV). The vertical S1–T1 gap amounts to 1.0 eV (exptl: 0.91 ± 0.1 eV). The electron localization function (ELF) and spin density analysis show that the singlet–triplet excitation is accompanied by a migration of electrons from the CC bond and N-atoms to the adjacent C–C and C–N bonds. A low ionization energy is confirmed for indigo IEa = 6.9 eV.The first excited state of indigo was characterized with a singlet–triplet energy gap of around 1.0 eV.
Co-reporter:G. Gopakumar, Vinh Son Nguyen, Minh Tho Nguyen
Journal of Molecular Structure: THEOCHEM 2007 Volume 811(1–3) pp:77-89
Publication Date(Web):1 June 2007
DOI:10.1016/j.theochem.2007.01.028
There is significant interest in the development of new materials for chemical hydrogen storage for the transportation sector. Electronic structure calculations using various methods, reveal that the borane molecule (BH3) could act as an efficient bifunctional acid–base catalyst in the H2-elimination reactions of XHnYHn systems (X, Y = B, C, N). Such a catalyst is needed as the elimination of H2 from isoelectronic ethane and borane amine compounds, proceeds with an energy barrier much higher than that of the A–B bond energy. The catalytic effect of BH3 has been probed by an analysis of the electronic densities of the transition structures using the atoms-in-molecules (AIM) and electron localization function (ELF) approaches.
Chemical Physics Letters 2006 Volume 419(1–3) pp:139-143
Publication Date(Web):15 February 2006
DOI:10.1016/j.cplett.2005.11.063
Abstract
MS/MS/MS experiments were carried out using a large-scale tandem mass spectrometer featuring a sector-quadrupole-sector geometry to probe the and radical cations. These ions (m/z 92 and 66) were produced from metastable ionized precursors or charge exchange reactions in the quadrupole cell; subsequent fragmentations and double ionization were performed after mass-selection and collisional activation with oxygen. Charge stripping spectra show significant variations in the relative intensities of the same doubly charged species generated in different ways. These variations are likely to be associated with the inherent differences in internal energy content.
Chemistry - A European Journal 2006 Volume 12(Issue 31) pp:
Publication Date(Web):26 JUL 2006
DOI:10.1002/chem.200600228
High-level ab initio (CCSD(T), CBS-QB3 and CASSCF, CASPT2, MR-ACPF, MR-ACPF-2) and density functional theory (B3LYP) calculations were carried out to study the dimerization of phosphaacetylene or phosphaethyne (HCP). Seventeen low energy closed-shell and five open-shell phosphaacetylene dimers were found on the potential energy surface. Two head-to-head, one head-to-tail and three other dimerization reaction pathways were determined, all with high activation barriers, suggesting that closed-shell minima are usually kinetically stable. An open-shell head-to-head reaction pathway has also been found with moderate initial barrier (95.0 kJ mol−1) leading to 1,2- and 1,3-diphosphacyclobutadiene, suggesting that polymerization of HCP and oligomerization of its derivatives have open-shell mechanism. Formation of 1,2-diphosphacyclobutadiene is both thermodynamically and kinetically favored over 1,3-diphosphacyclobutadiene. A head-to-head reaction involving LiBr as a catalyst was also studied. It has been pointed out that LiBr catalyze the closed-shell mechanism. All the four possible reaction channels of this reaction yield 1,4-diphosphatriafulvene with a fairly low activation Gibbs-free energy (44.8 kJ mol−1), suggesting that this compound could be synthesized. This finding fully supports the experimental results.
Co-reporter:Hue Minh Thi Nguyen, Tibor Höltzl, G. Gopakumar, Tamás Veszprémi, Jozef Peeters, Minh Tho Nguyen
Chemical Physics 2005 Volume 316(1–3) pp:125-140
Publication Date(Web):19 September 2005
DOI:10.1016/j.chemphys.2005.05.006
Abstract
Quantum chemical calculations were applied to investigate the electronic structure of the parent 1,3,5-tridehydrobenzene triradical (C6H3, TDB) and its anion , cation and protonated form . Our results obtained using the state-averaged complete active space self-consistent-field (CASSCF) followed by second-order multi-state multi-configuration perturbation theory, MS-CASPT2, and MRMP2 in conjunction with the large ANO-L and 6-311++G(3df,2p) basis set, confirm and reveal the followings: (i) TDB has a doublet 2A1 ground state with a 4B2–2A1 energy gap of 29 kcal/mol, (ii) the ground state of the anion in the triplet 3B2 being 4 kcal/mol below the 1A1 state. (iii) the electron affinity (EA), ionization energy (IE) and proton affinity (PA) are computed to be: EA = 1.6 eV, IE = 7.2 eV, PA = 227 kcal/mol using UB3LYP/6-311++G(3df,2p) + ZPE; standard heat of formation ΔHf(298 K, 1 atm)(TDB) = 179 ± 2 kcal/mol was calculated with CBS-QB3 method. An atoms-in-molecules (AIM) analysis of the structure reveals that the topology of the electron density is similar in all compounds: hydrogens connect to a six-membered ring, except for the case of the 2A2 state of (MBZ+) which is bicyclic with fused five- and three-membered rings. Properties of the chemical bonds were characterized with Electron Localization Function (ELF) analysis, as well as Wiberg indices, Laplacian and spin density maps. We found that the radicals form separate monosynaptic basins on the ELF space, however its pair character remains high. In the 2A1 state of TDB, the radical center is mainly localized on the C1 atom, while in the 2B2 state it is equally distributed between the C3 and C5 atoms and, due to the symmetry, in the 4B2 state the C1, C2 and C3 atoms have the same radical character. There is no C3–C5 bond in the 2A1 state of TDB, but the interaction between these atoms is strong. The ground state of cation (DHP), 1A1, is not a diradical and has a doubly aromatic character. Aromaticity of the different compounds was studied within the ELF framework and the standard deviation of the bond lengths and bond orders. The Jahn–Teller distorted 2A1 and 2B2(C2v) states of TDB were found to exhibit an aromaticity comparable to that of benzene. Overall protonation of the TDB reinforces the stability of the low-spin doublet states, the classical Hund’s rule is not obeyed. In a view, these species could better be regarded as radicals than triradicals.
Chemical Physics Letters 2005 Volume 411(4–6) pp:450-456
Publication Date(Web):15 August 2005
DOI:10.1016/j.cplett.2005.06.063
Abstract
1,3,5-Trimethylenebenzene (TMB, C9H9) is a prototypical hydrocarbon triradical. CASPT2/ANO-L calculations confirm that TMB has a quartet ground state with a gap of 12.4 ± 2.0 kcal/mol; both 2B1 and 2A2 states are quasi-degenerate. has nearly degenerate 3B2, 3A1 and 1A1 states and nearly degenerate 3A1, 1A1 and 3B2 states. Exocyclic protonation is favoured over ring-protonation; both processes induce a low-spin doublet ground state. Ionization energy, electron affinity, proton affinity are: IEa(TMB) = 7.5 ± 0.2 eV, EA(TMB) = 22 ± 2 kcal/mol, PA(C9H9) = 235 ± 2 kcal/mol.
Co-reporter:Truong Ba Tai, Long Van Duong, Hung Tan Pham, Dang Thi Tuyet Mai and Minh Tho Nguyen
Chemical Communications 2014 - vol. 50(Issue 13) pp:NaN1560-1560
Publication Date(Web):2013/12/03
DOI:10.1039/C3CC48392D
The B30 boron cluster has a bowl rather than a double-ring or a triple-ring tubular structure. This bowl isomer exhibits disk-aromaticity similar to that found for B202− and B19− clusters. We confirmed that the concept of disk-aromaticity can be applied to both planar and non-planar systems.
Chemical Communications 2015 - vol. 51(Issue 36) pp:NaN7680-7680
Publication Date(Web):2015/03/24
DOI:10.1039/C5CC01252J
The neutral B32 exhibits an aromatic bowl structure containing one heptagonal hole, while two anionic species, one having a bowl structure and the other a quasi-planar structure, are almost degenerate in energy. These findings not only give more insight into the structural features of boron clusters, but also present a key result explaining the presence of heptagonal holes in the fullerene B40.
Chemical Communications 2013 - vol. 49(Issue 9) pp:NaN915-915
Publication Date(Web):2012/12/10
DOI:10.1039/C2CC38038B
The B6Li8 cluster is a symmetrical 3D complex whose high stability can be understood through the Wade rule and aromaticity. A new mechanism of B–Li chemical bonding is proposed. Importantly, B6Li8 is predicted to be a promising candidate for hydrogen storage material with gravimetric density reaching up to a theoretical limit of 24%.
Co-reporter:Truong Ba Tai, Vu Thi Thu Huong and Minh Tho Nguyen
Chemical Communications 2013 - vol. 49(Issue 98) pp:NaN11550-11550
Publication Date(Web):2013/10/17
DOI:10.1039/C3CC47573E
The heteropolycyclic compounds containing borole units were theoretically designed. The presence of electron deficient boron atoms results in full electron delocalization and remarkably affects their aromaticity. While molecules 1 and 2a exhibit antiaromaticity for inner rings and non-aromaticity for outer rings, 2b and 2c are completely aromatic.
Electronic structure calculations suggest that hydrazine bisalane (AlH3NH2NH2AlH3, alhyzal) is a promising compound for chemical hydrogen storage (CHS). Calculations are carried out using the coupled-cluster theory CCSD(T) with the aug-cc-pVTZ basis set. Potential energy surfaces are constructed to probe the formation of, and hydrogen release from, hydrazine bisalane which is initially formed from the reaction of hydrazine with dialane. Molecular and electronic characteristics of both gauche and transalhyzal are determined for the first time. The gauche hydrazine bisalane is formed from starting reactants hydrazine + dialane following a movement of an AlH3 group from AlH3AlH3NH2NH2 rather than by a direct attachment of a separate AlH3 group, generated by predissociation of dialane, to AlH3NH2NH2. The energy barriers for dehydrogenation processes from gauche and transalhyzal are in the range of 21–28 kcal mol−1, which are substantially smaller than those of ca. 40 kcal mol−1 previously determined for the isovalent hydrazine bisborane (bhyzb) system. H2 release from hydrazine bisalane is thus more favored over that from hydrazine bisborane, making the Al derivative an alternative candidate for CHS.
Co-reporter:Truong Ba Tai, Sang Uck Lee and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 17) pp:NaN11623-11623
Publication Date(Web):2016/02/29
DOI:10.1039/C5CP07342A
The cage-like structures containing octagonal holes are located as the lowest-lying isomers for the B0/+42. The presence of octagonal holes, which have been found for the first time, not only gives us new insight into the bonding motif, but also marks a breakthrough in the structural characteristics of boron clusters since they were never expected to be stable units for elemental clusters. These cages are composed of both delocalized σ and π electron systems that consequently make them aromatic and thermodynamically stable.
Co-reporter:Hung Tan Pham, Kie Zen Lim, Remco W. A. Havenith and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 17) pp:NaN11931-11931
Publication Date(Web):2016/02/25
DOI:10.1039/C5CP07391J
The planarity of small boron-based clusters is the result of an interplay between geometry, electron delocalization, covalent bonding and stability. These compounds contain two different bonding patterns involving both σ and π delocalized bonds, and up to now, their aromaticity has been assigned mainly using the classical (4N + 2) electron count for both types of electrons. In the present study, we reexplored the aromatic feature of different types of planar boron-based clusters making use of the ring current approach. B3+/−, B42−, B5+/−, B6, B7−, B82−, B9−, B102−, B11−, B12, B13+, B142− and B162− are characterized by magnetic responses to be doubly σ and π aromatic species in which the π aromaticity can be predicted using the (4N + 2) electron count. The triply aromatic character of B12 and B13+ is confirmed. The π electrons of B182−, B19− and B202− obey the disk aromaticity rule with an electronic configuration of [1σ21π41δ42σ2] rather than the (4N + 2) count. The double aromaticity feature is observed for boron hydride cycles including B@B5H5+, Li7B5H5 and M@BnHnq clusters from both the (4N + 2) rule and ring current maps. The double π and σ aromaticity in carbon-boron planar cycles B7C−, B8C, B6C2, B9C−, B8C2 and B7C3− is in conflict with the Hückel electron count. This is also the case for the ions B11C5+/− whose ring current indicators suggest that they belong to the class of double aromaticity, in which the π electrons obey the disk aromaticity characteristics. In many clusters, the classical electron count cannot be applied, and the magnetic responses of the electron density expressed in terms of the ring current provide us with a more consistent criterion for determining their aromatic character.
Co-reporter:Vu Thi Thu Huong, Truong Ba Tai and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 8) pp:NaN6267-6267
Publication Date(Web):2016/01/27
DOI:10.1039/C5CP07277H
Organic semiconducting materials play an important role in the fabrication of high performance organic electronic devices. In the present work, we theoretically designed a series of organic semiconductors based on nickel complexes. Their characteristics of charge transport were investigated using DFT computational approaches. Based on the computed results, all compounds designed are found to be excellent candidates for ambipolar organic semiconductors with low reorganization energies for both holes and electrons. The (I–V) characteristics and transmission spectra of materials show that the replacement of benzene rings by thiophene rings results in an increase of their HOMO and LUMO energy levels. HOMOs of compounds containing thiophene end-groups are likely dominant for their conductance, while LUMOs of compounds containing benzene end-groups mainly affect their conductance. The electron distribution in these frontier MOs is identified as the main reason which makes the conductance of the compounds in the first series higher than those in the later series.
Co-reporter:Yassin Aweis Jeilani, Phoenix N. Williams, Sofia Walton and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2016 - vol. 18(Issue 30) pp:NaN20188-20188
Publication Date(Web):2016/05/10
DOI:10.1039/C6CP02686A
The reaction pathways for the prebiotic formation of nucleobases are complex and lead to the formation of a mixture of products. In the past 50 years, there has been a concerted effort for identifying a unified mechanism for the abiotic origin of the biomolecules but with little success. In the present theoretical study, we identified two prominent precursors for the building up of RNA and DNA nucleobases under prebiotic conditions: (a) 1,2-diaminomaleonitrile (DAMN), which is a tetramer of hydrogen cyanide (HCN), and (b) formamide, a hydrolysis product of HCN; it is important to emphasize that HCN is the source of both precursors. We find that free radical pathways are potentially appropriate to account for the origin of nucleobases from HCN. The current study unites the formamide pathways with the DAMN pathways. The mechanisms for the formation of the RNA and DNA nucleobases (uracil, adenine, purine, cytosine) were studied by quantum chemical computations using density functional theory at the B3LYP/6-311G(d,p) level. All the routes involved proceed with relatively low energy barriers (within the error margin of DFT methods). We showed that the radical mechanisms for the formation of nucleobases could be unified through common precursors. The results demonstrated that 4-aminoimidazole-5-carbonitrile (AICN), which is a known precursor for nucleobases, is a product of DAMN. The overall mechanisms are internally consistent with the abiotic formation of the nucleobases, namely (a) under a meteoritic impact scenario on the early Earth's surface that generated high internal energy, and/or (b) in the (gas phase) interstellar regions without the presence of catalysts.
Co-reporter:Nguyen Minh Tam, Hung Tan Pham, Long Van Duong, My Phuong Pham-Ho and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2015 - vol. 17(Issue 5) pp:NaN3003-3003
Publication Date(Web):2014/12/16
DOI:10.1039/C4CP04279D
Stabilized fullerene and tubular forms can be produced in boron clusters Bn in small sizes from n ∼ 14 to 20 upon doping by transition metal atoms. B14Fe and B16Fe are stable tubes whereas B18Fe and B20Fe are stable fullerenes. Their formation and stability suggest the use of dopants to induce different growth paths leading to larger cages, fullerenes and tubes of boron.
Physical Chemistry Chemical Physics 2015 - vol. 17(Issue 20) pp:NaN13679-13679
Publication Date(Web):2015/04/22
DOI:10.1039/C5CP01851J
Boron clusters have been of great interest over the last few decades due to their unique chemical and physical properties. In the present work, we performed a theoretical study of geometrical and electronic structures of boron clusters Bn with n = 26–29 in both neutral and anionic states using DFT and MO computational methods. The photoelectron spectra of anionic species were simulated using TDDFT methods. Our results predict that in the neutral state both the B26 and B27 clusters exhibit tubular forms, whereas the larger species B28 and B29 are quasi-planar structures. The anionic species Bn− are more favourable for 2D shapes. More importantly, based on known geometrical characteristics, we now establish a general growth mechanism of boron clusters, which gives us more insight into the formation and existence of boron based nanomaterials.
Physical Chemistry Chemical Physics 2015 - vol. 17(Issue 26) pp:NaN16936-16936
Publication Date(Web):2015/06/01
DOI:10.1039/C5CP01456E
Thermally feasible decomposition pathways of formamide (FM) in the presence of vanadium VO(X4Σ−) and titanium TiO(X3Δ) monoxides are determined using density functional theory (BP86 functional) and coupled-cluster theory (CCSD(T)) computations with large basis sets. These diatomic metal oxides have been shown to be present in the prebiotic conditions. The dehydration, decarbonylation and dehydrogenation reactions of the molecular and dissociative complexes of FM and MO (M = V, Ti) turn out to be more favourable than those of the ground state isolated FM. The effect of addition of one or two water molecules on energy barriers is also probed for these reaction pathways. In some cases, a combined catalytic effect when adding water is observed. This enhanced catalytic effect was not observed in previously reported cases of FM transformation, for example, when adding water molecules into the mineral-catalyzed isomerizations of FM. The dehydration process of MO–FM complexes without the presence of water is found to be more feasible than the decarbonylation and dehydrogenation. The overall energy barrier for the non-water VO–FM dehydration is ∼3 kcal mol−1 lower than the reference energy of the separated systems, whereas those of the two latter reactions are higher than the reference. Although the TiO–FM dehydration has a larger overall barrier of 14 kcal mol−1 as compared to the VO–FM counterpart, the two other decomposition pathways still have much higher energy barriers. Direct formation of urea and H2CO from a FM dimer and indirect formation of urea from FM via the intermediate HNCO are also established. Urea formation in an indirect pathway is preferred. These low-energy-barrier pathways leading to the formation of important prebiotic molecules suggest that metal monoxides MO could play an important catalytic role in the prebiotic reactions of FM.
Physical Chemistry Chemical Physics 2015 - vol. 17(Issue 26) pp:NaN17345-17345
Publication Date(Web):2015/06/08
DOI:10.1039/C5CP01650A
Using density functional theory with the TPSSh functional and the 6-311+G(d) basis set, we extensively searched for the global minima of two metallic atoms doped boron clusters B6M2, B7M2, B12M2 and B14M2 with transition metal element M being Co and Fe. Structural identifications reveal that B7Co2, B7Fe2 and B7CoFe clusters have global minima in a B-cyclic motif, in which a perfectly planar B7 is coordinated with two metallic atoms placed along the C7 axis. The B6 cluster is too small to form a cycle with the presence of two metals. Similarly, the B12 cluster is not large enough to stabilize the metallic dimer within a double ring 2 × B6 tube. The doped B14M2 clusters including B14Co2, B14Fe2 and B14CoFe have a double ring 2 × B7 tubular shape in which one metal atom is encapsulated by the B14 tube and the other is located at an exposed position. Dissociation energies demonstrate that while bimetallic cyclic cluster B7M2 prefers a fragmentation channel that generates the B7 global minimum plus metallic dimer, the tubular structure B14M2 tends to dissociate giving a bimetallic cyclic structure B7M2 and a B@B6 cluster. The enhanced stability of the bimetallic doped boron clusters considered can be understood from the stabilizing interactions between the anti-bonding MOs of metal–metal dimers and the levels of a disk aromatic configuration (for bimetallic cyclic structures), or the eigenstates of the B14 tubular form (in case of bimetallic tubular structure).
Co-reporter:Hung Tan Pham, Thu-Thuy Phan, Nguyen Minh Tam, Long Van Duong, My Phuong Pham-Ho and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2015 - vol. 17(Issue 27) pp:NaN17570-17570
Publication Date(Web):2015/06/09
DOI:10.1039/C5CP02257F
The smallest triple ring tubular silicon cluster Mn2@Si15 is reported for the first time. Theoretical structural identification shows that the Mn2@Si15 tubular structure whose triple ring is composed by three five-membered Si rings in anti-prism motif, is stable in high symmetry (D5h) and singlet ground state (1A1′). The dimer Mn2 is placed inside the tubular along the C5 axis, and the Mn dopant form single Si–Mn bonds with Si skeleton, whereas the Mn–Mn is characterized as a triple bond. The effect of Mn2 on the stability of the Si15 triple ring structure arises from strong orbital overlap of Mn2 with Si15.
Co-reporter:Long Van Duong, Hung Tan Pham, Nguyen Minh Tam and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2014 - vol. 16(Issue 36) pp:NaN19478-19478
Publication Date(Web):2014/07/29
DOI:10.1039/C4CP01996B
We determined the geometries and chemical bonding phenomena of the B27 system in its dicationic, cationic, neutral, anionic and dianionic states using DFT computations. In both cationic and neutral states, the triple ring tubular forms correspond to the lowest-energy isomers, especially in B27+. The cation B27+ represents the first stable hollow cylinder having a triple ring among the pure boron clusters. In the anionic and dianionic states, the quasi-planar structures are favoured due to a charge effect. In the triple ring tube B27+, strong diatropic responses to external magnetic field occur in both radial and tangential types of electrons, and thus confer it a characteristic tubular aromaticity. The presence of a consistent aromatic character contributes to its high thermodynamic stability. The shapes of calculated MOs of B27+ TR can be predicted by the eigenstates of a simple model of a particle on a hollow cylinder. The number of electrons in a hollow cylinder should attain a number of (4N + 2M) with M = 0 and 1 for both radial and tangential electrons, depending on the number of non-degenerate MOs occupied, in order to properly fulfill the closed electron shells. In the case of B27+, M = 0 for radial electrons and M = 1 for tangential electrons.
Physical Chemistry Chemical Physics 2013 - vol. 15(Issue 37) pp:NaN15415-15415
Publication Date(Web):2013/07/10
DOI:10.1039/C3CP51017D
Optical properties of silver Agn nanoclusters are demonstrated to be dependent on their size, structure and charge state. It is found that when being contained in the sodalite cavity of LTA zeolite the tetradecanuclear hexacation silver cluster Ag146+ is stable. Its lower-lying states and optical spectrum are theoretically determined using the quantum chemical TD-DFT method. Its ground state possesses an outer-shell electron configuration of A1g2T2g6 mimicking the s2p6 valence of noble gas atoms. These frontier orbitals are constructed from 5s,5p(Ag)-AOs with contributions from framework oxygen atoms. Light absorption of Ag146+ embedded in the sodalite cage which is characterized by strong peaks centered at 331 and 476 nm (transitions 5s,p(Ag) → 5s,p(Ag)) leads to much longer wavelength emission. The sodalite cage, as a container, stabilizes the central Ag146+ cluster by electrostatic attraction. The absorption spectrum of the isovalent neutral Ag8 cluster embedded inside the same sodalite cavity is also simulated using TD-DFT and CASPT2 methods. This absorption spectrum which is similar to that of the Ag146+ cluster has two absorption bands in the near UV and visible regions.
Co-reporter:Yassin A. Jeilani, Huyen Thi Nguyen, Domnique Newallo, Jean-Marie D. Dimandja and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2013 - vol. 15(Issue 48) pp:NaN21093-21093
Publication Date(Web):2013/10/25
DOI:10.1039/C3CP53108B
Modeling the complicated chemical reactions in the interstellar medium and surface materials of Titan is nontrivial. Since both the atmosphere and the surface are rich in organic molecules, the chemistry may have important implications for the origin of biomolecules. Prebiotic synthesis of DNA nucleobases from simple molecules such as formamide has been known for more than half a century. In this study, new free radical pathways leading to the synthesis of guanine, hypoxanthine, purine, and adenine have been studied using density functional theory (B3LYP with the 6-311G(d,p) basis set). The pathways of the selected nucleobases demonstrate the importance of free radicals in the production of useful biomolecules under conditions appropriate for the interstellar medium or on Titan. The pathways may be universal in nature and proceed without solvent requirements. Calculations indicate that radical pathways yield lower reaction barriers as compared to previously reported pathways. Overall, these results suggest that the chemistry on Titan's surface and/or the growth of organic particulates in the haze layers in Titan's atmosphere likely involve free radicals. The mechanisms demonstrate that important prebiotic precursors can be predicted. The reaction sequences reported here may lead to the production and build-up of molecules with prebiotic relevance.
Co-reporter:Vu Thi Thu Huong, Truong Ba Tai and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2012 - vol. 14(Issue 43) pp:NaN14841-14841
Publication Date(Web):2012/08/02
DOI:10.1039/C2CP42474F
Following the theme of this special issue, two new compounds, the P-flowers C16(PH)8 and C16(PF)8, are designed by us and subsequently characterized by quantum chemical computations. Their geometries and infrared signatures are analyzed and compared to those of the well-known sulflower C16S8. Their electronic structure and aromaticity are examined using the electron localization function (ELF) and also by the total and partial densities of state (DOS). Both C16(PF)6 and C16(PH)8 molecules exhibit small energy barrier of electron injection (Φe = 0.33 eV for the gold electrode for the former, and Φe = 0.1 eV for the calcium electrode for the latter), remarkably low reorganization energy and high rate of electron hopping. Thus, both theoretically designed P-flower molecules are predicted to be excellent candidates for organic n-type semiconductors.
Co-reporter:Vinh Son Nguyen, Rehab M. Ibrahim Elsamra, Jozef Peeters, Shaun A. Carl and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2012 - vol. 14(Issue 20) pp:NaN7470-7470
Publication Date(Web):2012/03/29
DOI:10.1039/C2CP40367F
We investigated the rate constants and reaction mechanism of the gas phase reaction between the ethynyl radical and nitrous oxide (C2H + N2O) using both experimental methods and electronic structure calculations. A pulsed-laser photolysis/chemiluminescence technique was used to determine the absolute rate coefficient over the temperature range 570 K to 836 K. In this experimental temperature range, the measured temperature dependence of the overall rate constants can be expressed as: k(T) (C2H + N2O) = 2.93 × 10−11 exp((−4000 ± 1100) K/T) cm3 s−1 (95% statistical confidence). Portions of the C2H + N2O potential energy surface (PES), containing low-energy pathways, were constructed using the composite G3B3 method. A multi-step reaction route leading to the products HCCO + N2 is clearly preferred. The high selectivity between product channels favouring N2 formation occurs very early. The pathway corresponds to the addition of the terminal C atom of C2H to the terminal N atom of N2O. Refined calculations using the coupled-cluster theory whose electronic energies were extrapolated to the complete basis set limit CCSD(T)/CBS led to an energy barrier of 6.0 kcal mol−1 for the entrance channel. The overall rate constant was also determined by application of transition-state theory and Rice–Ramsperger–Kassel–Marcus (RRKM) statistical analyses to the PES. The computed rate constants have similar temperature dependence to the experimental values, though were somewhat lower.
Co-reporter:Vinh Son Nguyen, Saartje Swinnen, Jerzy Leszczynski and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2011 - vol. 13(Issue 14) pp:NaN6656-6656
Publication Date(Web):2011/03/07
DOI:10.1039/C0CP02484H
The reactivity of hydrazine in the presence of diborane has been investigated using ab initio quantum chemical computations (MP2 and CCSD(T) methods with the aug-cc-pVTZ basis set). Portions of the relevant potential energy surface were constructed to probe the formation mechanism of the hydrazine diborane (BH3BH3NH2NH2) and hydrazine bisborane (BH3NH2NH2BH3). The differences between both adducts are established. The release of hydrogen molecules from hydrazine bisborane adducts has also been characterized. Our results suggest that the BH3NH2NH2BH3 adduct, which has been prepared experimentally, is formed from the starting reactants hydrazine + diborane. The observed adduct is produced by a transfer of a BH3 group from BH3BH3NH2NH2 rather than by the direct attachment of a separate BH3 group, generated by predissociation of diborane, to BH3NH2NH2.
Co-reporter:Nguyen Tien Trung, Nguyen Phi Hung, Tran Thanh Hue and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2011 - vol. 13(Issue 31) pp:NaN14042-14042
Publication Date(Web):2011/06/16
DOI:10.1039/C1CP20533A
In this study, 16 gas phase complexes of the pairs of XCHZ and CO2 (X = F, Cl, Br; Z = O, S) have been identified. Interaction energies calculated at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ level including both BSSE and ZPE corrections range from −5.6 to −10.5 kJ mol−1 for XCHO⋯CO2 and from −5.7 to −9.1 kJ mol−1 for XCHS⋯CO2. Substitution of one H atom by one halogen in formaldehyde and thioformaldehyde reduces the interaction energy of XCHZ⋯CO2, while a CH3 substitution increases the interaction energy of both CH3CHO⋯CO2 and CH3CHS⋯CO2. NBO and AIM analyses also point out that the strength of Lewis acid–base interactions decreases going from >C1S3⋯C6 to >C1O3⋯C6 and to >C1–X4⋯C6. This result suggests the higher capacity of solubility of thiocarbonyl compounds in scCO2, providing an enormous potential application for designing CO2-philic materials based on the >CS functional group in competition with >CO. The Lewis acid–base interaction of the types >CS⋯C, >C–Cl⋯C and >C–Br⋯C is demonstrated for the first time. The contribution of the hydrogen bonding interaction to the total interaction energy is larger for XCHS⋯CO2 than for XCHO⋯CO2. Upon complexation, a contraction of the C1–H2 bond length and a blue shift of its stretching frequency have been observed, as compared to the isolated monomer, indicating the existence of a blue-shifting hydrogen bond in all complexes examined. Calculated results also lend further support for the viewpoint that when acting as proton donor, a C–H bond having a weaker polarization will induce a stronger distance contraction and frequency blue shift upon complexation, and vice versa.
Physical Chemistry Chemical Physics 2011 - vol. 13(Issue 36) pp:NaN16264-16264
Publication Date(Web):2011/08/11
DOI:10.1039/C1CP22078K
A systematic quantum chemical investigation on the electronic, geometric and energetic properties of AunV clusters with n = 1–14 in both neutral and anionic states is performed using BP86/cc-pVTZ-PP calculations. Most clusters having an even number of electrons prefer a high spin state. For odd-electron systems, a quartet state is consistently favoured as the ground state up to Au8V. The larger sized Au10V, Au12V and Au14V prefer a doublet state. The clusters prefer 2D geometries up to Au8V involving a weak charge transfer. The larger systems bear 3D conformations with a more effective electron transfer from Au to V. The lowest-energy structure of a size AunV is built upon the most stable form of Aun−1V. During the growth, V is endohedrally doped in order to maximize its coordination numbers and augment the charge transfer. Energetic properties, including the binding energies, embedding energies and second-order energy differences, show that the presence of a V atom enhances considerably the thermodynamic stability of odd-numbered gold clusters but reduces that of even-numbered systems. The atomic shape has an apparently more important effect on the clusters stability than the electronic structure. Especially, if both atomic shape and electronic condition are satisfied, the resulting cluster becomes particularly stable such as the anion Au12V−, which can thus combine with the cation Au+ to form a superatomic molecule of the type [Au12V]Au. Numerous lower-lying electronic states of these clusters are very close in energy, in such a way that DFT computations cannot clearly establish their ground electronic states. Calculated results demonstrate the existence of structural isomers with comparable energy content for several species including Au9V, Au10V, Au13V and Au14V.
Co-reporter:Truong Ba Tai, Pham Vu Nhat and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2010 - vol. 12(Issue 37) pp:NaN11486-11486
Publication Date(Web):2010/08/02
DOI:10.1039/C003401K
Quantum chemical calculations are performed on the aluminium doped lithium clusters LinAl at both neutral and cationic states using the DFT/B3LYP and CCSD(T) methods in conjugation with the aug-cc-pVaZ (a = D,T,Q) basis sets. The global minima are located and the growth mechanism is established. The electronic structure, geometrical parameters and energetic properties, such as average binding energy Eb, second difference of energy Δ2E, adiabatic and vertical ionization energy, and dissociated enthalpy, are evaluated using the coupled-cluster CCSD(T) method, whose energies are extrapolated to the complete basis set limit (CBS). The high stability of Li5Al, Li7Al, Li6Al+ and Li8Al+ that have the “magic numbers” of valence electrons, can be understood using the phenomenological shell model.
Co-reporter:Ling Lin, Pieterjan Claes, Tibor Höltzl, Ewald Janssens, Torsten Wende, Risshu Bergmann, Gabriele Santambrogio, Gerard Meijer, Knut R. Asmis, Minh Tho Nguyen and Peter Lievens
Physical Chemistry Chemical Physics 2010 - vol. 12(Issue 42) pp:NaN13913-13913
Publication Date(Web):2010/09/20
DOI:10.1039/C0CP00911C
The geometric and electronic structure of the Au6Y+ cation is studied by gas phase vibrational spectroscopy combined with density functional theory calculations. The infrared photodissociation spectrum of Au6Y+·Ne is measured in the 95–225 cm−1 energy range and exhibits two characteristic absorption bands at 181 cm−1 and 121 cm−1. Based on DFT/BP86 quantum chemical calculations, the infrared spectrum is assigned to the lowest energy species found, an eclipsed C3v geometry. The 3D structure of Au6Y+ is considerably different from those previously found for both the neutral Au6Y (quasi-planar circular geometry) and the anionic Au6Y− (planar D6h symmetry). The different geometries are related to different electronic structures in agreement with 2D and 3D phenomenological shell models for metal clusters.
Co-reporter:Nguyen Tien Trung, Tran Thanh Hue and Minh Tho Nguyen
Physical Chemistry Chemical Physics 2009 - vol. 11(Issue 6) pp:NaN933-933
Publication Date(Web):2008/12/15
DOI:10.1039/B816112G
The hydrogen-bonded complexes formed from interaction of trihalomethanes CHX3 (X = F, Cl, Br) with nitrosyl hydride HNO were studied using ab initio MO calculations at the second-order perturbation theory (MP2/6-311++G(d,p)). Each interaction contains at least five separate equilibrium structures. Calculated binding energies range from 4 to 8 kJ mol−1 with both ZPE and BSSE corrections. While CHBr3 leads to the most stable complexes with HNO, CHF3 forms the least stable counterparts. The strength of complexes thus tends to increase from F to Cl to Br, which is consistent with a decrease of deprotonation enthalpy of the corresponding C–H bonds. It is remarkable that all the C–H and N–H bonds are shortened upon complexation, giving rise to an increase of their stretching frequencies. A blue shift is thus observed for the N–H bonds of the type N–H⋯X (X = F, Cl, Br); such a contraction of the covalent N–H bond is extremely rare.