NiTi过渡金属团簇的结构和芳香性的理论研究

在B3PW91/6-311+G(d)计算水平上, 计算并讨论了Ni₄Ti₂, [Ni₄Ti₂]²⁺, [Ni₄Ti₂]^2-与Ni₄Ti₄, [Ni₄Ti₄]²⁺, [Ni₄Ti₄]^2-团簇的几何结构和芳香性. 在构型优化过程中得到了Ni₄Ti₂(D_4h), [Ni₄Ti₂]²⁺(D_4h), [Ni₄Ti₂]^2-(D_4h)和Ni₄Ti₄(D_2h)4个稳定构型, 发现当引入上下2个Ti原子后, Ni₄环成为了平面正方形构型. 核无关化学位移(NICS)计算结果表明, Ni₄Ti₂(D_4h)与Ni₄Ti₄(D_2h)的NICS值为正, 而[Ni₄Ti₂]²⁺(D_4h)和[Ni₄Ti₂]^2-(D_4h)的NICS值为负, 且[Ni₄Ti₂]^2-(D_4h)的NICS值更负. 同时还发现, 由s与d轨道参与形成的反磁性环流是引起[Ni₄Ti₂]²⁺(D_4h)和[Ni₄Ti₂]^2-(D_4h)具有较大芳香性的主要原因; 其中Ti原子主要提供d_z²与s轨道, 而Ni原子主要利用其d_z²与d_x²-y²轨道形成正方形环, 它们之间构成了球状的d轨道环流, 且[Ni₄Ti₂]²⁺(D_4h)和[Ni₄Ti₂]^2-(D_4h)中还有非常明显的π轨道环流.     

Computational studies on the dimers and the thermal dimerization of norbornadiene

In this study, 14 norbornadiene (NBD) dimers and the thermal dimerization mechanism were studied using the hybrid density functional theory (B3LYP) and the second-order multiconfigurational perturbation theory (CASPT2). In the process of dimerization, the biradical stationary points were located using the unrestricted, broken-spin, symmetry approach. The pathways were divided into eight parts to aid the analysis of their mechanisms. Our results indicated that the process for the formation of the cage-like heptacyclo[6.6.0.0.260.3130.411059.01014] tetradecane (HCTD, D14) is highly exothermic (92.15 kcal/mol), indicating that D14 is the most stable NBD dimer. However, the formation of D14 is very difficult to achieve kinetically because of a higher barrier in the thermal dimerization. On the contrary, the formation of exo-cis-exo (D5) is kinetically favorable, but thermodynamically unfavorable at higher temperature. Therefore, the combination of both thermodynamic and kinetic factors indicated that the formation of exo-exo (D9), which resembles the product of the pseudo-Diels-Alder addition, is most likely in the NBD dimerization.     

Massively parallel implementation of 3D‐RISM calculation with volumetric 3D‐FFT

A new three‐dimensional reference interaction site model (3D‐RISM) program for massively parallel machines combined with the volumetric 3D fast Fourier transform (3D‐FFT) was developed, and tested on the RIKEN K supercomputer. The ordinary parallel 3D‐RISM program has a limitation on the number of parallelizations because of the limitations of the slab‐type 3D‐FFT. The volumetric 3D‐FFT relieves this limitation drastically. We tested the 3D‐RISM calculation on the large and fine calculation cell (2048³ grid points) on 16,384 nodes, each having eight CPU cores. The new 3D‐RISM program achieved excellent scalability to the parallelization, running on the RIKEN K supercomputer. As a benchmark application, we employed the program, combined with molecular dynamics simulation, to analyze the oligomerization process of chymotrypsin Inhibitor 2 mutant. The results demonstrate that the massive parallel 3D‐RISM program is effective to analyze the hydration properties of the large biomolecular systems. © 2014 Wiley Periodicals, Inc.     

Geometries, Electronic Structures, and Electron Detachment Energies of Small Boron Sulfide Anions： A Density Functional Theory Investigation

A density functional theory investigation on the geometries, electronic structures, and electron detachment energies of BS, BS2, B（BS）2 and B（BS）3 has been performed in this work. The linear ground-state structures of BS （C∞v, ^1∑^＋） and BS2^- （O∞h, ^1∑g^＋） prove to be similar to the previously reported BO and BO2 with systematically lower electron detachment energies. Small boron sulfide clusters are found to favor the formation of -B=S groups which function basically as a-radicals and dominate the ground-state structures of the systems. The perfect linear B（BS）2^-（D∞h, ^3∑g） and beautiful equilateral triangle B（BS）3^- （D3h,^2A1”） turn out to be analogous to the well-known C2v BH2 and O3h BH3, respectively. The electron affinities of BS, BS2, B（BS）2 and B（BS）3 are predicted to be 2.3, 3.69, 3.00 and 3.45 eV, respectively. The electron detachment energies calculated for BS^-, BS2^-, B（BS）2^-, and B（BS）3^- may facilitate future photoelectron spectroscopy measurements to characterize the geometrical and electronic structures of these anions.     

From trifluoroacetate complex precursors to monodisperse rare-earth fluoride and oxyfluoride nanocrystals with diverse shapes through controlled fluorination in solution phase

We report the first systematic synthesis of monodisperse rare-earth (RE=La to Lu, Y) fluoride and oxyfluoride nanocrystals with diverse shapes (trigonal REF3 triangular, truncated-triangular, hexagonal, and polygonal nanoplates; orthorhombic REF3 quadrilateral and zigzag-shaped nanoplates; cubic REOF nanopolyhedra and nanorods) from single-source precursors (SSP) of [RE(CF(3)COO)(3)] through controlled fluorination in oleic acid (OA)/oleylamine (OM)/1-octadecene (ODE). To selectively obtain REF3 or REOF nanocrystals, the fluorination of the RE-O bond to the RE-F bond at the nucleation stage was controlled by finely tuning the ratio of OA/ODE or OA/OM, and the reaction temperature. For phase-pure REF3 or REOF naocrystals, their shape-selective syntheses could be realized by further modifying the reaction conditions. The two-dimensional growth of the REF3 nanoplates and the one-dimensional growth of the REOF nanorods were likely due to the selective adsorption of the capping ligands on specific crystal planes of the nanocrystals. Those well-shaped nanocrystals with diverse geometric symmetries (such as D(3h), D(6h), C(2h), O(h), and D(nh)) displayed a remarkable capability to form self-assembled superlattices. By manipulating the solvent-substrate combination, the plate-shaped REF3 nanocrystals could form highly ordered nanoarrays by means of either the face-to-face formation or the edge-to-edge formation. By using this SSP strategy, we also obtained high-quality LaF3:Eu and LaF3:Eu/LaF3 triangular nanoplates that showed photoluminescent red emissions of Eu3+ ions sensitive to the surface effect.     

Sc3N@C80: computations on the two-isomer equilibrium at high temperatures.

Reported herein are computations on the relative concentrations of the two experimentally known isomers of Sc3N@C80 , that is, those produced by encapsulation of Sc3N in two particular C80 cages that obey the isolated-pentagon rule, namely, with I(h) and D(5h) symmetries. The calculations are based on density functional methods and have been carried out using the Gibbs energy over a broad temperature interval. It has been computed that, if a relatively free motion of the encapsulate inside the cages is allowed, the observed populations of 10 and 17 % for the D(5h) Sc3N@C80 species are reached at temperatures of 2100 and 2450 K, respectively. The inclusion of the entropy term is essential as, if it is neglected, the D(5h) Sc3N@C80 population at a temperature of 2100 K would be a mere 1 %, owing to the relatively large interisomeric separation potential energy of 19 kcal mol(-1).     

Coordination field analysis of rare earth complexes with triangular symmetry

The calculation of the complex matrixes in odd triangular symmetry was accomplished.The configurations of the coordination unit with various triangular symmetries and different ligand numbers were discussed.On the basis of the double-sphere coordination point-charge (DSCPCF) model,the detailed forms of the DSCPCF parameters Bmk and the expressions of the perturbation matrix elements in triangular field (D3,D3h,D3d) were derived.Thereby,the calculation scheme of coordination field perturbation energy of the rare earth complexes with triangular symmetry was constructed After the calculation scheme was programmed,the Stark energies of the crystalline TbAl3(BO3)4 were calculated The results were considerably close to the experimental values     

(ClAlNH)n簇合物的结构与稳定性

用量子化学从头算方法（HF／6－31G）和密度泛函理论（DFT）的B3LYP方法，以6－31G标准基组加一个极化函数，对（ClAlNH)n(n=1-10)簇合物的几何构型，电子结构和红外光谱进行了优化，并讨论了聚合反应（ClAlNH)m→(ClAlNH)n的热力学效应，结果表明，（ClAlNH)n系列簇合物的基态稳定结构为Cs(n=1),D2h(n=2),D3h(n=3),Td(n=4),Cs(n=5),D3d(n=6),Cs(n=7),S4(n=8),D3h(n=9),C2h(n=10,n=2,4,6,8,10等偶数对应的（ClAlNH)n簇化合物的结构比n等于奇数量更稳定。     

A Density Functional Study on the Geometries of Compounds Fe(HCN)_n~+ (n = 1～6)

1 INTRODUCTION The solvation of metal ion by different types ofsolvents is of great interest for a wide variety of app-lications[1]. In the experimental and theoretical inve-stigations, most of such studies are focused on ion-ligand systems complexed by…     

High‐resolution mass spectrometry and hydrogen/deuterium exchange study of mitorubrin azaphilones and nitrogenized analogues

Azaphilones represent numerous groups of wild fungal secondary metabolites that exhibit exceptional tendency to bind to nitrogen atoms in various molecules, especially those containing the amine group. Nitrogenized analogues of mitorubrin azaphilones, natural secondary metabolites of Hypoxylon fragiforme fungus, have been detected in the fungal methanol extract in very low concentrations. Positive electrospray ionization interfaced with high‐resolution mass spectrometry was applied for confirmation of the elemental composition of protonated species. Collision‐induced dissociation (CID) experiments have been performed, and fragmentation mechanisms have been proposed. Additional information regarding both secondary metabolite analogue families has been reached by application of gas‐phase proton/deuterium (H/D) exchanges performed in the collision cell of a triple quadrupole mass spectrometer. An incomplete H/D exchange with one proton less than expected was observed for both protonated mitorubrin azaphilones and their nitrogenized analogues. By means of the density functional theory, an appropriate explanation of this behavior was provided, and it revealed some information concerning gas‐phase H/D exchange mechanism and protonation sites. Copyright © 2012 John Wiley & Sons, Ltd.     

Evolution of lone pair of excess electrons inside molecular cages with the deformation of the cage in e2@C60F60 systems

For unusual e(2)@C(60)F(60)(I(h), D(6h), and D(5d)) cage structures with two excess electrons, it is reported that not only the lone pair in singlet state but also two single excess electrons in triplet state can be encapsulated inside the C(60)F(60) cages to form single molecular solvated dielectrons. The interesting relationship between the shape of the cage and the spin state of the system has revealed that ground states are singlet state for spherical shaped e(2)@C(60)F(60)(I(h)) and triplet states for short capsular shaped e(2)@C(60)F(60)(D(6h)) and long capsular shaped e(2)@C(60)F(60)(D(5d)), which shows a spin evolution from the singlet to triplet state with the deformation of the cage from spherical to capsular shape. For these excess electron systems, the three ground state structures have large vertical electron detachment energies (VDEs (I) of 1.720-2.283 eV and VDEs (II) of 3.959-5.288 eV), which shows their stabilities and suggests that the large C(60)F(60) cage is the efficient container of excess electrons.     

A solvation‐free‐energy functional: A reference‐modified density functional formulation

The three‐dimensional reference interaction site model (3D‐RISM) theory, which is one of the most applicable integral equation theories for molecular liquids, overestimates the absolute values of solvation‐free‐energy (SFE) for large solute molecules in water. To improve the free‐energy density functional for the SFE of solute molecules, we propose a reference‐modified density functional theory (RMDFT) that is a general theoretical approach to construct the free‐energy density functional systematically. In the RMDFT formulation, hard‐sphere (HS) fluids are introduced as the reference system instead of an ideal polyatomic molecular gas, which has been regarded as the appropriate reference system of the interaction‐site‐model density functional theory for polyatomic molecular fluids. We show that using RMDFT with a reference HS system can significantly improve the absolute values of the SFE for a set of neutral amino acid side‐chain analogues as well as for 504 small organic molecules. © 2015 Wiley Periodicals, Inc.     

A Density Functional Study on the Geometries of Compounds Fe(HCN)n+(n = 1～6)

The possible geometries of Fe(HCN)n+ (n = 1～6) compounds were studied by using DFT/UB3LYP/6-31G(2df) method. The structure and ground state of each fragmental ion are C∞v (4∑+ or 6∑+), D∞h(4∑g+), D3h (4A1'), C2v or Td or C3v (4A1), and D3h (4A1') or C4v( (2A1') sequentially with n = 1～5. For the compound Fe(HCN)6+, the possible geometry was not obtained. The sequential incremental interaction energy (-△(△E)), dissociation energy (△D0), enthalpy (-△(△H)) and Gibbs free energy (-△(△G)), and frequencies for HCN-Fe(HCN)n-1+ were also calculated, and the results are all in good agreement with the experiments. The bond length of Fe-N is lengthened with the increase of cluster size, and the strength of Fe+-N coordination bond varies nonmon- tonically as increasing the number of ligands. The Fe+-N bond of Fe(HCN)2+ is the strongest in all compounds.     

Role of doping and sheet size in tailoring optoelectronic properties of germanene: A TDDFT study

Germanene is a novel 2D material with promising optoelectronic properties, tuning of which is to be explored. This work demonstrates that doping and increasing the sheet size can alter optical and electronic properties of germanene via perturbation of the band structure. This feature has also been observed in other nanostructures, notably, silicon nanostructures, and may be attributed to quantum confinement effects. Our main findings on H‐terminated germanene are, (i) band gap can be reduced by 30%, (ii) exciton binding energy can be reduced by 60%, and (iii) absorption spectra can be tuned from UV to visible range. We employ time‐dependent density functional theory to investigate the role of dopants, boron (B), phosphorus (P), carbon (C), silicon (Si), and zirconium (Zr). Width of the germanene sheet is varied from 0.78 nm to 2.78 nm. Frequency and energy calculations are carried out to analyze the infrared (IR) and ultra‐violet (UV)‐visible (VIS) spectra.     

Accurate reaction barrier heights of pericyclic reactions: Surprisingly large deviations for the CBS‐QB3 composite method and their consequences in DFT benchmark studies

Accurate barrier heights are obtained for the 26 pericyclic reactions in the BHPERI dataset by means of the high‐level Wn‐F12 thermochemical protocols. Very often, the complete basis set (CBS)‐type composite methods are used in similar situations, but herein it is shown that they in fact result in surprisingly large errors with root mean square deviations (RMSDs) of about 2.5 kcal mol^?1. In comparison, other composite methods, particularly G4‐type and estimated coupled cluster with singles, doubles, and quasiperturbative triple excitations [CCSD(T)/CBS] approaches, show deviations well below the chemical‐accuracy threshold of 1 kcal mol^?1. With the exception of SCS‐MP2 and the herein newly introduced MP3.5 approach, all other tested Møller‐Plesset perturbative procedures give poor performance with RMSDs of up to 8.0 kcal mol^?1. The finding that CBS‐type methods fail for barrier heights of these reactions is unexpected and it is particularly troublesome given that they are often used to obtain reference values for benchmark studies. Significant differences are identified in the interpretation and final ranking of density functional theory (DFT) methods when using the original CBS‐QB3 rather than the new Wn‐F12 reference values for BHPERI. In particular, it is observed that the more accurate Wn‐F12 benchmark results in lower statistical errors for those methods that are generally considered to be robust and accurate. Two examples are the PW6B95‐D3(BJ) hybrid‐meta‐general‐gradient approximation and the PWPB95‐D3(BJ) double‐hybrid functionals, which result in the lowest RMSDs of the entire DFT study (1.3 and 1.0 kcal mol^?1, respectively). These results indicate that CBS‐QB3 should be applied with caution in computational modeling and benchmark studies involving related systems. © 2015 Wiley Periodicals, Inc.     

Interaction mechanisms between poly(amido‐amine) and nano‐silicon dioxide

To gain insight into the attachment of ?Si⁺ (SC) and ?SiO^? (SOA) ions (regarded as guests) to the lowest generation, ? NH₂‐terminated poly(amidoamine) (PAMAM) dendrimers (regarded as host) in the gas phase, density functional theory is used to investigate the structures and energetics of the complexes with B3LYP/6‐31+G (d) and HF/6‐31G basis sets. The initial parameters are obtained through the initial optimizations at the HF level using the most basic STO‐3G basis set. Various initial configurations of the ions bound to PAMAM are tested, and four stable conformers are found, i.e., types A to D. Types 1A and 2C are the most stable due to the chemical bond formations of Si? N° and Si? O, respectively. For type B, SC coordination to amide O sites occurs via electrostatic induction. For type D, SOA coordination to amide hydrogen and amine hydrogen sites occurs via hydrogen bond interaction. Spatial hindrance, electrostatic induction force, and hydrogen‐bond interaction play important roles in the complexation process. © 2013 Wiley Periodicals, Inc. This article was published online on 5 July 2012. An error was subsequently identified. This notice is included in the online and print version to indicate that both have been corrected on 3 August 2012.     

Calculations of solid‐state 43Ca NMR parameters: A comparison of periodic and cluster approaches and an evaluation of DFT functionals

We present a computational study of magnetic‐shielding and quadrupolar‐coupling tensors of ⁴³Ca sites in crystalline solids. A comparison between periodic and cluster‐based approaches for modeling solid‐state interactions demonstrates that cluster‐based approaches are suitable for predicting ⁴³Ca NMR parameters. Several model chemistries, including Hartree–Fock theory and 17 DFT approximations (SVWN, CA‐PZ, PBE, PBE0, PW91, B3PW91, rPBE, PBEsol, WC, PKZB, BMK, M06‐L, M06, M06‐2X, M06‐HF, TPSS, and TPSSh), are evaluated for the prediction of ⁴³Ca NMR parameters. Convergence of NMR parameters with respect to basis sets of the form cc‐pVXZ (X = D, T, Q) is also evaluated. All DFT methods lead to substantial, and frequently systematic, overestimations of experimental chemical shifts. Hartree–Fock calculations outperform all DFT methods for the prediction of ⁴³Ca chemical‐shift tensors. © 2017 Wiley Periodicals, Inc.     

Efficient treatment of solvation shells in 3D molecular theory of solvation

We developed a technique to decrease memory requirements when solving the integral equations of three‐dimensional (3D) molecular theory of solvation, a.k.a. 3D reference interaction site model (3D‐RISM), using the modified direct inversion in the iterative subspace (MDIIS) numerical method of generalized minimal residual type. The latter provides robust convergence, in particular, for charged systems and electrolyte solutions with strong associative effects for which damped iterations do not converge. The MDIIS solver (typically, with 2 × 10 iterative vectors of argument and residual for fast convergence) treats the solute excluded volume (core), while handling the solvation shells in the 3D box with two vectors coupled with MDIIS iteratively and incorporating the electrostatic asymptotics outside the box analytically. For solvated systems from small to large macromolecules and solid–liquid interfaces, this results in 6‐ to 16‐fold memory reduction and corresponding CPU load decrease in MDIIS. We illustrated the new technique on solvated systems of chemical and biomolecular relevance with different dimensionality, both in ambient water and aqueous electrolyte solution, by solving the 3D‐RISM equations with the Kovalenko–Hirata (KH) closure, and the hypernetted chain (HNC) closure where convergent. This core–shell‐asymptotics technique coupling MDIIS for the excluded volume core with iteration of the solvation shells converges as efficiently as MDIIS for the whole 3D box and yields the solvation structure and thermodynamics without loss of accuracy. Although being of benefit for solutes of any size, this memory reduction becomes critical in 3D‐RISM calculations for large solvated systems, such as macromolecules in solution with ions, ligands, and other cofactors. © 2012 Wiley Periodicals, Inc.     

Theoretical study of metal‐free organic dyes based on different configurations for efficient dye‐sensitized solar cells

Considering different solar dyes configuration, four novel metal‐free organic dyes based on phenoxazine as electron donor, thiophene and cyanovinylene linkers as the ‐conjugation bridge and cyanoacrylic acid as electron acceptor were designed to optimize open circuit voltage and short circuit current parameters and theoretically inspected. Density functional theory and time‐dependent density functional theory calculations were used to study frontier molecular orbital energy states of the dyes and their optical absorption spectra. The results indicated that D2‐4 dyes can be suitable candidates as sensitizers for application in dye sensitized solar cells and among these three dyes, D3 showed a broader and more bathochromically shifted absorption band compared to the others. The dye also showed the highest molar extinction coefficient. This work suggests optimizing the configuration of metal‐free organic dyes based on simple D‐ ‐A configuration containing alkyl chain as substitution, starburst conformation, and symmetric double D‐ ‐A chains would produce good photovoltaic properties.     

Xanthone O‐Glycosides from the Roots of Pteris Multifida

Two new xanthone glycosides and six known compounds were isolated from the roots of Pteris multifida. Based on spectroscopic and chemical methods, the structures of the new compounds were elucidated as 1‐hydroxy‐4,7‐dimethoxy‐8‐(3‐methyl‐2‐butenyl)‐6‐O‐α‐L‐rhamnopyranosyl‐(1→2)‐[β‐D‐glucopyranosyl‐(1→3)]‐β‐D‐glucopyranosylxanthone ( 1 ), and 1,3‐dihydroxy‐7‐methoxy‐8‐(3‐methyl‐2‐butenyl)‐6‐O‐α‐L‐rhamnopyranosyl‐(1 →2)‐[β‐D‐glucopyranosyl‐(1→3)]‐β‐D‐glucopyranosylxanthone ( 2 ), respectively.