Hypercoordinated carbon in C-doped boron fullerenes: a quantum chemical study

The structures and stability of C-doped boron fullerenes with the three-dimensional arrangement of non-classical pentacoordinated quasi-flat carbon centers were studied using the density functional theory (DFT) B3LYP/6-311+G(d,p) method. The doping with carbon atoms in apical positions above the five-membered rings stabilizes the spherical boron fullerene forms due to multicenter interactions of p_z-orbitals of the carbons and adjacent boron atoms. Increasing in the size of the fullerene cluster is accompanied by change in the bonding pattern and by flattening of the hypercoordinated carbon centers. Endohedral metal atoms significantly affect on the structure and stability of the fullerene systems with hypercoordinated carbon centers.     

Hydrogen storage in magnesium clusters: quantum chemical study

Magnesium hydride is cheap and contains 7.7 wt % hydrogen, making it one of the most attractive hydrogen storage materials. However, thermodynamics dictate that hydrogen desorption from bulk magnesium hydride only takes place at or above 300 degrees C, which is a major impediment for practical application. A few results in the literature, related to disordered materials and very thin layers, indicate that lower desorption temperatures are possible. We systematically investigated the effect of crystal grain size on the thermodynamic stability of magnesium and magnesium hydride, using ab initio Hartree-Fock and density functional theory calculations. Also, the stepwise desorption of hydrogen was followed in detail. As expected, both magnesium and magnesium hydride become less stable with decreasing cluster size, notably for clusters smaller than 20 magnesium atoms. However, magnesium hydride destabilizes more strongly than magnesium. As a result, the hydrogen desorption energy decreases significantly when the crystal grain size becomes smaller than approximately 1.3 nm. For instance, an MgH2 crystallite size of 0.9 nm corresponds to a desorption temperature of only 200 degrees C. This predicted decrease of the hydrogen desorption temperature is an important step toward the application of Mg as a hydrogen storage material.     

Stabilization of platinum clusters by substitutional boron dopants in carbon supports

To stabilize platinum clusters on carbon supports, carbon substituted boron dopants were introduced in graphite and fullerene models, and the binding energies of Pt(1) to Pt(6) clusters on the carbon supports were investigated using first-principles density functional theory (DFT) calculations. The adsorption energies of the Pt clusters increased significantly in the presence of boron. The adsorption energies of Ru(1) and Au(1) were also found to dramatically increase in the presence of boron.     

CO adsorption on pure and binary-alloy gold clusters: a quantum chemical study

We performed density-functional theory analysis of nondissociative CO adsorption on 22 binary Au-alloy (Au(n)M(m)) clusters: n=0-3, m=0-3, and m+n=2 (dimers) or 3 (trimers), M=Cu/Ag/Pd/Pt. We report basis-set superposition error corrections to adsorption energies and include both internal energy of adsorption (DeltaU(ads)) and Gibbs free energy of adsorption (DeltaG(ads)) at standard conditions (298.15 K and 1 atm). We found onefold (atop) CO binding on all the clusters except Pd2 (twofold/bridged), Pt2 (twofold/bridged), and Pd3 (threefold). In agreement with the experimental results, we found that CO adsorption is thermodynamically favorable on pure Au/Cu clusters but not on pure Ag clusters and also observed the following adsorption affinity trend: Pd>Pt>Au>Cu>Ag. For alloy dimers we found the following patterns: Au2>M Au>M2 (M=Ag/Cu) and M2>M Au>Au2 (M=Pd/Pt). Alloying Ag/Cu dimers with (more reactive) Au enhanced adsorption and the opposite effect was observed for PdPt dimers. The Ag-Au, Cu-Au, and Pd-Au trimers followed the trends observed on dimers: Au3>M Au2>M2Au>M3 (M=Ag/Cu) and Pd3>Pd2Au>PdAu2>Au3. Interestingly, Pt-Au trimers reacted differently and alloying with Au systematically increased the adsorption affinity: PtAu2>Pt2Au>Pt3>Au3. A strikingly different behavior of Pt is also manifested by the triplet spin state and onefold (atop) binding in Pt3-CO which is in contradiction with the singlet spin state and threefold binding in Pd3-CO. We found a linear correlation between CO binding energy (BE) and elongation of the CO bond. For Ag-Au and Cu-Au clusters, the increase in CO BE (and elongation of the C-O bond which is probably due to the back donation) is accompanied by the decrease in the cluster-CO distance suggesting that the donation (from 5sigma highest occupied molecular orbital in CO to cluster lowest unoccupied molecular orbital) mechanism also contributes to the BE. For Pd-Au clusters, the cluster-CO distance (and CO bond length) increases with increase in the BE, suggesting that the donation mechanism may not be important for those clusters. No clear trend was observed for Pt-Au clusters.     

Stabilization of silver clusters by matrices of various chemical nature

Silver clusters having unusual optical and chemical properties and characterized by a band at 320–330 nm, are shown to be stabilized by solid matrices of various chemical nature.

---

, , 320–330 , .

     

Stabilization of coordinatively unsaturated Ir4 clusters with bulky ligands: a comparative study of chemical and mechanical effects

The synthesis and characterization of new cluster compounds represented by the series Ir(4)(CO)(12-x)L(x) (L = tert-butyl-calix[4]-arene(OPr)(3)(OCH(2)PPh(2)); x = 2 and 3) is reported using ESI mass spectrometry, NMR spectroscopy, IR spectroscopy and single-crystal X-ray diffraction. Thermally driven decarbonylation of the cluster compound series represented by x = 1-3 according to the formula above is followed via FTIR and NMR spectroscopies, and dynamic light scattering in toluene solution. The propensity of these clusters to decarbonylate in solution is shown to be directly correlated with number density of adsorbed calixarene phosphine ligands and controlled via Pauli repulsion between metal d and CO 5σ orbitals. The tendency for cluster aggregation unintuitively follows a trend that is exactly opposite to the cluster's propensity to decarbonylate. No cluster aggregation is observed for clusters consisting of x = 3, even after extensive decarbonylation via loss of all bridging CO ligands and coordinative unsaturation. Some of the CO lost during thermal treatment via decarbonylation can be rebound to the coordinatively unsaturated cluster consisting of x = 3. In contrast, the clusters consisting of x = 1 and x = 2 both aggregate into large nanoparticles when treated under identical conditions. Clusters in which the calixarene phosphine ligand is replaced with a sterically less demanding PPh(2)Me ligand 6 lead to significantly less coordinative unsaturation upon thermal treatment. Altogether, these data support a mechanical model of accessibility in coordinatively unsaturated metal clusters in solution, which hinges on having at least three sterically bulky organic ligands per Ir(4) core.     

Silver clusters in the cages of zeolites: A quantum chemical study

The electronic structure of zeolite A is developed in a step by step procedure from the simple O_hH₈Si₈O₁₂ molecule, to the _∞ ¹ [(-O)₂H₄[Si₈O₁₂)] chain, to the _∞ ² [(-O₄)(Si₈O₁₂)] layer, and finally to the silica zeolite A framework _∞ ³ [Si₂₄O₄₈]. It is remarkable how well the calculated band structures of both, _∞ ² [(-O)₄(Si₈O₁₂)] and _∞ ³ [Si₂₄O₄₈] correspond to the experimentally determined band structure of α-quartz with a Fermi level of -10.55 eV. The HOMO region consists in each case of nonbonding 2p-oxygen bands which in a localized language can be denoted as oxygen lone pairs ( | O<). We observe in each case the typical behaviour of an insulator with saturated valencies whose electronic structure can be described as being localized and is already present in the starting O_h-H₈Si₈O₁₂ molecule. The double-8-rings D8R of the _∞ ² [(-O)₄(Si₈O₁₂)] layer have a pore diameter of 4.1 Å, the same as the pore opening of zeolite A. It is large enough to accept up to four Ag, forming _∞ ² [(-O)₄(Si₈O₁₂)Ag_n], n = 1, 2, 3, 4, layers, suitable for modelling the electronic interactions between the zeolite cavity embedded silver clusters and between the clusters and the zeolite framework. With one Ag per D8R the band structure is simply a superposition of the 4d, 5s and 5p levels of a layer of nearly noninteracting Ag and the silicon dioxide layer. The Ag-d band lies below the oxygen lone pairs, the Ag-s band lies about 3 eV above the oxygen lone pairs, and the Ag-5p bands are in the antibonding silicon dioxide region. The first electronic transition is of oxygen lone pair to Ag-5s LMCT type. Increasing silver content results in progressive splitting of the 5σ Ag bands and shifts the first (Ag^m+ _n)? ← (| O<) charge transfer transition to lower energies. The filled Ag 4d-bands lie always significantly below the (| O<) HOCOs (highest occupied crystal orbitals) but their band width increases with increasing silver content. In all cases the zeolite environment separates the Ag clusters through antibonding Ag-(← O<) interactions so that the coupling remains weak and it makes sense to describe the Ag clusters in the D8R as quantum dots weakly interacting with each other.     

Photocurrent increase by doping a liquid crystal host with a functionalized fullerene

Doping a photoconductive mesogenic material with a small amount of functionalized fullerene significantly increases the photocurrent. LESR studies show the appearance of a long-lived charge-separated state under light illumination, with a negative charge on the fullerene moiety. This confirms that in the system studied the fullerene unit acts as an electron trap. In the lamellar structure of the liquid crystalline phase the fullerene units and mesogenic cores are separated, which ensures better space separation between negative and positive charges. As a result the charge recombination is slowed and the photocurrent is amplified.     

Molecular inclusion of PCB126 by beta-cyclodextrin: a combined molecular dynamics simulation and quantum chemical study

The effective enrichment and identification of lowly concentrated polychlorinated biphenyls (PCBs) in the environment is attracting much research attention due to human health concerns raised from their emissions. Cyclodextrins (CDs) are known to be capable to form inclusion complexes with a variety of organic molecules. The purpose of this study is to provide theoretical evidences whether CDs can form energetically stable inclusion complexes with PCBs through a host–guest interaction, and if so, whether infrared and Raman techniques are suitable for the detection of CD-modified PCBs. Focusing on a representative PCB molecule, 3,3′,4,4′,5-pentachlorobiphenyl (PCB126), we studied its molecular inclusion by β-CD (BCD) by performing molecular dynamics simulations and density functional theory calculations. Calculated results show that PCB126 and BCD preferentially form the stable 1:1 inclusion complex. The calculated IR spectra of the 1:1 inclusion complexes mainly present the spectra features of BCD and give only a slight indication for bands of the guest molecule. In contrast, the characteristic vibration modes of the guest molecule are remarkably prominent in the Raman spectra of the inclusion complexes. Based on the present results, we propose that BCD can potentially serve as a candidate for including PCB126 to form the stable 1:1 host–guest complex, and that Raman spectroscopy technology is expected to be suitable for the identification of the CD-modified PCBs, whereas IR spectroscopy is not feasible for such an application.     

Making and breaking of small water clusters: A combined quantum chemical and molecular dynamics approach

We present a combined quantum chemical and molecular dynamics study of cyclic and noncyclic water n-mers ([(H₂O]_n, n = 2–6) at four different temperatures and showcase that the dynamics of small water clusters can reproduce the known properties of bulk water reasonably well. We investigate the making and breaking of the water clusters by computing the hydrogen bond strengths, average lifetimes, and relative stabilities, which are important to understand the complex solution dynamics. We compare the behavior of water clusters in the gas phase and in the solution phase as well as the variation in the properties as a function of cluster size and highlight the notably more interesting cluster dynamics of the water trimer when compared to the other water clusters. © 2019 Wiley Periodicals, Inc.     

Dynamical hydrogen atom tunneling in dichlorotropolone: a combined quantum, semiclassical, and classical study

Based on the Cartesian Reaction Surface framework we construct a four-dimensional potential for the tropolone derivative 3,7-dichlorotropolone, a molecule with an intramolecular O-H...O hydrogen bond. The reduced configuration space involves the in-plane hydrogen atom coordinates, a symmetric O-O vibrational mode, and an antisymmetric mode related to deformations of the seven-membered ring. The system is characterized in terms of quantum mechanical computations of the low-lying eigenstates as well as a classical and semiclassical analysis of spectra obtained via Fourier transforming autocorrelation functions. For the semiclassical analysis we utilize the amplitude-free correlation function method [K. Hotta and K. Takatsuka, J. Phys. A 36, 4785 (2003)]. Our results demonstrate substantial anharmonic couplings leading to highly correlated wave functions even at moderate energies. Furthermore, the importance of dynamical tunneling in tropolone is suggested since many low-lying states--including the ground state--lie above the classical saddle point but nevertheless appear as split pairs.     

Environmental effects on the spectroscopic properties of gallic acid: a combined classical and quantum mechanical study

The solvation of gallic acid (in water and acetonitrile) is studied by means of its spectroscopic properties. IR, UV, and NMR spectra are predicted by using various solvation models obtained in terms of both purely classical and density functional approaches. Comparison with experiments is used to validate solvation models. Hydrogen-bond and long-range (or bulk) effects are evaluated by comparing different solvation models. A continuum-only approach, a purely discrete, and a mixed continuum/discrete approach based on quantum-mechanical and classical molecular-dynamics solute-solvent clusters are tested.     

interaction of nitro compounds with radicals: a quantum chemical study

Photoreduction of nitro compounds is accompanied by formation of various radical products that can react with the starting nitro compound, thus causing deviation of the decomposition kinetics from the first-order kinetics with respect to the nitro compound. The results of quantum chemical modeling of the reactions of nitro compounds with radicals and the pathways of further transformations of radical adducts formed in the reactions are presented. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 2, pp. 202–206, February, 2006.     

A quantum chemistry study: a new kind of boron nitrides

At the B3LYP/6-311+G** and the BP86/6-311+G** levels of theory, BNN, H(3)BNN, NNBH(2)-BH(2)NN, (BNN)(2)H(2), NNBBNN, (BNN)(3) (+), (BNN)(4), (BNN)(4) (2+), (BNN)(4) (2-), B(4)(NN)(2), (BNN)(5) (-), (BNN)(6), (BNN)(7) (+), and (BNN)(8) (2+) are investigated. Neutral (BNN)(4) is aromatic with its triplet state but antiaromatic with its singlet state. (BNN)(4) dication favors D(2d) structure, while (BNN)(4) dianion favors a planar D(4h) structure. (BNN)(3) (+), (BNN)(4) (2-), (BNN)(5) (-), (BNN)(6), (BNN)(7) (+), (BNN)(7) (3-), (BNN)(8) (2+), and (BNN)(8) (2-) are all aromatic with planar monocyclic conformation, following the 4n + 2 rule. Moreover, according to the CASSCF and MRCI calculations, the planar B(4)(NN)(2) of D(2h) symmetry prefers to be a sigma-pi diradical in spite of open-shell singlet or triplet and is also aromatic. Akin to the sigma-pi back interaction in compounds containing transition metal, there exists the sigma-pi back interaction between boron and N(2) ligand among some species reported herein, which strengthens B(-)N bond but activates N-N bond, especially in (4)Sigma(-) BNN. The T-shaped structure lies lowest in energy among seven isomers of the (BNN)(2)H(2) dimer, and the parallel-displaced structure is favored between two isomers of the (BNN)(6) dimer.     

一系列新型硼氮化合物的理论研究

借助密度泛函理论,在B3LYP/6-311 G**和BP86/6-311 G**水平上,对一系列包含BNN单元的硼氮化合物BNN,H3BNN,NNBH2-BH2NN,(BNN)2H2,NN-BBNN,(BNN)3 ,(BNN)42 ,(BNN)5-,(BNN)6,(BNN)7 和(BNN)82 进行了理论计算研究.结果表明:在这些硼氮化合物中,尤其是BNN,B原子和N2配体之间存在着σ-π反馈作用.这种作用增强了B-N键,而活化了N-N键.结果还显示,(BNN)3 ,(BNN)42 ,(BNN)5-,(BNN)6,(BNN)7 和(BNN)82 都是具有平面单环结构并且满足4n 2规则的芳香性体系.     

The coordination of uranyl in water: a combined quantum chemical and molecular simulation study

The coordination environment of uranyl in water has been studied using a combined quantum mechanical and molecular dynamics approach. Multiconfigurational wave function calculations have been performed to generate pair potentials between uranyl and water. The quantum chemically determined energies have been used to fit parameters in a polarizable force field with an added charge transfer term. Molecular dynamics simulations have been performed for the uranyl ion and up to 400 water molecules. The results show a uranyl ion with five water molecules coordinated in the equatorial plane. The U-O(H(2)O) distance is 2.40 A, which is close to the experimental estimates. A second coordination shell starts at about 4.7 A from the uranium atom. No hydrogen bonding is found between the uranyl oxygens and water. Exchange of waters between the first and second solvation shell is found to occur through a path intermediate between association and interchange. This is the first fully ab initio determination of the solvation of the uranyl ion in water.     

Elongation of planar boron clusters by hydrogenation: boron analogues of polyenes

Dihydrogenated boron clusters, H(2)B(n)(-) (n = 7-12), were produced and characterized using photoelectron spectroscopy and computational chemistry to have ladderlike structures terminated by a hydrogen atom on each end. The two rows of boron atoms in the dihydrides are bonded by delocalized three-, four-, or five-center σ and π bonds. The π bonding patterns in these boron nanoladders bear similarities to those in conjugated alkenes: H(2)B(7)(-), H(2)B(8), and H(2)B(9)(-), each with two π bonds, are similar to butadiene, while H(2)B(10)(2-), H(2)B(11)(-), and H(2)B(12), each with three π bonds, are analogous to 1,3,5-hexatriene. The boron cluster dihydrides can thus be considered as polyene analogues, or "polyboroenes". Long polyboroenes with conjugated π bonds (analogous to polyacetylenes), which may form a new class of molecular wires, should exist.     

Proton mobility in protonated glycylglycine and N-formylglycylglycinamide: a combined quantum chemical and RKKM study

Theoretical model calculations were performed to investigate the degree of validity of the mobile proton model of protonated peptides. The structures and energies of the most important minima corresponding to different structural isomers of protonated diglycine and their conformers, as well as the barriers separating them, were determined by DFT calculations. The rate coefficients of the proton transfer reactions between the isomers were calculated using the RRKM method in order to obtain a quantitative measure of the time scale of these processes. The proton transfer reactions were found to be very fast already at and above the threshold to the lowest energy decomposition pathway. Two possible mechanisms of b2+-ion formation via water loss from the dipeptide are also discussed. The rate-determining step of the proton migration along a peptide chain is also investigated using the model compound N-formylglycylglycinamide. The investigations revealed that this process very possibly occurs via the protonation of the carbonyl oxygens of the amide bonds, and its rate-determining step is an internal rotation-type transition of the protonated C=O-H group between two adjacent C=O-HellipsisO=C bridges.     

Absorption and fluorescence of PRODAN in phospholipid bilayers: a combined quantum mechanics and classical molecular dynamics study

Absorption and fluorescence spectra of PRODAN (6-propionyl-2-dimethylaminonaphthalene) were studied by means of the time-dependent density functional theory and the algebraic diagrammatic construction method. The influence of environment, a phosphatidylcholine lipid bilayer and water, was taken into account employing a combination of quantum chemical calculations with empirical force-field molecular dynamics simulations. Additionally, experimental absorption and emission spectra of PRODAN were measured in cyclohexane, water, and lipid vesicles. Both planar and twisted configurations of the first excited state of PRODAN were taken into account. The twisted structure is stabilized in both water and a lipid bilayer, and should be considered as an emitting state in polar environments. Orientation of the excited dye in the lipid bilayer significantly depends on configuration. In the bilayer, the fluorescence spectrum can be regarded as a combination of emission from both planar and twisted structures.