采用优化的DFTB参数对铜(111)表面碳二聚化的分子动力学研究（英文）

针对铜表面化学反应,我们发展了一套铜-碳体系的密度泛函紧束缚(DFTB)参数。测试结果表明这套参数可以很好的描述吸附铜或碳原子前后铜表面的几何结构和能量。基于这套参数,我们对Cu(111)表面的碳二聚化过程进行了分子模拟研究。即使在高温下,直接的分子动力学模拟也很难观察到碳二聚体的形成。这是因为高温下铜表面显著的结构弛豫一定程度上阻止了二聚化。为了研究高温下铜表面碳二聚化的机理,我们进行了赝动力学模拟。发现在二聚化的过程中,碳原子形成C-Cu-C桥状结构以后,会绕中间Cu原子转动,最后形成碳二聚体。1300 K下碳二聚化的自由能垒约0.9 eV。     

氧化锌有序结构在Au(111)和Cu(111)上的生长（英文）

利用NO_2或O_2作为氧化剂,研究了氧化锌在Au(111)和Cu(111)上的生长和结构。NO_2表现了更好的氧化性能,有利于有序氧化锌纳米结构或薄膜的生长。在Au(111)和Cu(111)这两个表面上,化学计量比氧化锌都形成非极性的平面化ZnO(0001)的表面结构。在Au(111)上,NO_2气氛下室温沉积锌倾向于形成双层氧化锌纳米结构;而在更高的沉积温度下,在NO_2气氛中沉积锌则可同时观测到单层和双层氧化锌纳米结构。O_2作为氧化剂时可导致形成亚化学计量比的ZnOx结构。由于铜和锌之间的强相互作用会促进锌的体相扩散,并且铜表面可以被氧化形成表面氧化物,整层氧化锌在Cu(111)上的生长相当困难。我们通过使用NO_2作为氧化剂解决了这个问题,生长出了覆盖Cu(111)表面的满层有序氧化锌薄膜。这些有序氧化锌薄膜表面显示出莫尔条纹,表明存在一个ZnO和Cu(111)之间的莫尔超晶格。实验上观察到的超晶格结构与最近理论计算提出的Cu(111)上的氧化锌薄膜结构相符,具有最小应力。我们的研究表明,氧化锌薄膜的表界面结构可能会随氧化程度或氧化剂的不同而变化,而Cu(111)的表面氧化也可能影响氧化锌的生长。当Cu(111)表面被预氧化成铜表面氧化物时,ZnOx的生长模式会发生变化,锌原子会受到铜氧化物晶格的限域形成单位点锌。我们的研究表明了氧化锌的生长需要抑制锌向金属基底的扩散,并阻止亚化学计量比ZnOx的形成。因此,使用原子氧源有利于在Au(111)和Cu(111)表面上生长有序氧化锌薄膜。     

氧化锌有序结构在Au(111)和Cu(111)上的生长

利用NO₂或O₂作为氧化剂,研究了氧化锌在Au(111)和Cu(111)上的生长和结构。NO₂表现了更好的氧化性能,有利于有序氧化锌纳米结构或薄膜的生长。在Au(111)和Cu(111)这两个表面上,化学计量比氧化锌都形成非极性的平面化ZnO(0001)的表面结构。在Au(111)上,NO₂气氛下室温沉积锌倾向于形成双层氧化锌纳米结构;而在更高的沉积温度下,在NO₂气氛中沉积锌则可同时观测到单层和双层氧化锌纳米结构。O₂作为氧化剂时可导致形成亚化学计量比的ZnO_x结构。由于铜和锌之间的强相互作用会促进锌的体相扩散,并且铜表面可以被氧化形成表面氧化物,整层氧化锌在Cu(111)上的生长相当困难。我们通过使用NO₂作为氧化剂解决了这个问题,生长出了覆盖Cu(111)表面的满层有序氧化锌薄膜。这些有序氧化锌薄膜表面显示出莫尔条纹,表明存在一个ZnO和Cu(111)之间的莫尔超晶格。实验上观察到的超晶格结构与最近理论计算提出的Cu(111)上的氧化锌薄膜结构相符,具有最小应力。我们的研究表明,氧化锌薄膜的表界面结构可能会随氧化程度或氧化剂的不同而变化,而Cu(111)的表面氧化也可能影响氧化锌的生长。当Cu(111)表面被预氧化成铜表面氧化物时,ZnO_x的生长模式会发生变化,锌原子会受到铜氧化物晶格的限域形成单位点锌。我们的研究表明了氧化锌的生长需要抑制锌向金属基底的扩散,并阻止亚化学计量比ZnO_x的形成。因此,使用原子氧源有利于在Au(111)和Cu(111)表面上生长有序氧化锌薄膜。     

NO双分子和二聚体与Cu2作用的理论计算

采用密度泛函理论(DFT)中的B3LYP方法，在Lanl2DZ基组下，对NO双分子和二聚体与铜原子簇相互作用的结构进行了研究. 结果表明，NO可以在铜表面相邻的两个铜原子上形成稳定的双分子吸附和二聚体吸附，而在双分子吸附形式中NO以氮原子吸附在铜上的构型最稳定，且顶点吸附的稳定性不如非顶点吸附形式.在二聚体吸附形式中， N－N键被加强，而N－O键被削弱的程度大于双分子吸附形式，说明二聚体的形成有利于NO在金属铜表面的直接分解.同时电荷布居分析表明，单重态的二聚体与铜作用时，铜原子上的平均电荷达到0.66 e,说明在这种吸附形式中铜被离子化的倾向较大，而且这种吸附形式最有利于NO的分解.这些结果说明NO经二聚体形式在铜表面直接催化分解是可行的.     

油页岩的~(13)C-NMR特征及FLASHCHAIN热解模拟研究

采用固体13C-NMR核磁共振技术表征了甘肃窑街矿区油页岩的碳骨架结构,分析并计算了油母质团簇化学结构参数,包括团簇平均碳原子数、芳碳原子数、脂碳原子数及芳环数。在热重红外分析仪(TG-FTIR)上进行了油页岩的热解实验,得到了热解产物的生成规律。结合样品的团簇化学结构参数,采用基于油页岩结构的FLASHCHAIN模型模拟其热解产物的生成过程;模拟结果与TG-FITR实验数据符合较好,印证了模型预测的合理性。     

CO32-对羟基乙叉二膦酸体系铜电结晶的影响

采用循环伏安和计时安培法研究了羟基乙叉二膦酸(HEDPA)镀铜液中铜在玻碳电极上电结晶的初期行为。结果表明：羟基乙叉二膦酸(HEDPA)镀铜体系中,铜的电沉积过程经历了晶核形成过程;当溶液中不含CO32-时,其电结晶按连续三维成核方式进行,而CO32-的加入,使得铜电结晶按瞬时三维成核方式进行;成核数密度都随着电位的提高而增加。这可能是CO32- -以第二配体形式进入HEDPA和Cu2+构成的络合结构,从而形成更稳定的络合物吸附在电极表面所致。     

固体碳的一种新形态——富勒烯

一富勒烯(Fullerene)的发现人类对碳元素的研究和应用亘古至今,源远流长。人们熟知碳有两种同素异形体——石墨和金刚石。直到80年代中期,发现了富勒烯碳原子簇,尤其是近年来,对富勒烯的结构、性质深入广泛地研究,确认碳元素还存在着第三种晶体形态,已是顺理成章了。碳原子簇的研究始于天体物理学家对宇宙尘埃形成的兴趣,为了模拟星际空间及恒星附近碳原子簇的形成过程,早在1942年O Hahn等用质谱法证实了原子簇C_n(n<15)的存在。1984年Rohlfing等用质谱仪研究在超声氦气流中以激光蒸发石墨所得产物时,发现碳可以形成n<200的C_n原子簇。当n>40时,簇中碳原子教仅为偶数,并且还发现C_(60)的质谱峰明显高     

全氟并五苯在半金属Ga表面的分子二聚化及自组装单层

利用超高真空扫描隧道显微镜(UHV-STM)和有机分子束沉积(OMBD)方法研究了全氟并五苯(perfluoropentancene,PFP)分子在半金属Ga表面的吸附和两维自组装. 在低覆盖度下单个PFP分子在Ga表面上表现出很高的迁移性. 在1分子单层(monolayer, ML)时PFP分子发生二聚化并在 Ga 表面上无序排列. 轻度热退火可导致PFP两维自组装: 二聚体排列为高度有序的一维分子带阵列, 带中 PFP二聚体排列为砖墙(brick wall)结构. 在高分辨 STM图中, PFP分子两端出现亮暗相反的圆形突起, 并且相邻分子的亮暗极性相反, 表明PFP分子带有电偶极矩, PFP二聚体带有电四极矩. 因此, PFP分子二聚体的形成机制可唯像解释为反向电偶极矩之间的静电吸引作用; 二聚体的砖墙排列结构可归结为同向电四极矩之间的静电排斥作用.     

羊毛铬青R-联吡啶光度法测定痕量铜

铜联吡啶羊毛铬青R在非离子表面活性剂存在下可形成两种组成不同的三元络合物。利用1:1:2络合物测定铜已有过报导。我们发现改用阴离子表面活性剂,并适当增大它的浓度后,可抑制1:1:2络合物的形成。据此,我们利用形成1:1:1络合物的反应,拟订了一个测定铜的新方法。羊毛铬青R(ECR)-联吡啶(BPY)-铜可形成两种组成不同的三元络合物,西田宏等研究了它的组成,并利用1:1:2络合物(Cu:BPY:ECR)测定铜。我们在研究1:1:2络合     

CuNaY分子筛的有效吸附位与其脱硫性能的关联性研究

以液相离子交换法制备了一系列不同Cu负载量的Cu Na Y分子筛;采用XRD及N2吸附-脱附表征分子筛的微观结构和织构性质,采用动态吸附法考察其对噻吩模拟油的吸附脱硫性能,结合NH_3-TPD和Py-FTIR方法对CuNaY分子筛的酸量和有效Cu~+物种进行定量分析,研究了CuNaY分子筛的表面酸性和铜物种形态结构对其吸附脱硫性能的影响机制。结果表明,通过改变铜负载量可有效调控改性Y分子筛的表面酸性以及铜物种化学形态;适量铜物种的引入可以最大限度的形成有效吸附位,从而获得最优吸附脱硫性能,而过量的Cu物种会在Y分子筛笼内形成多核铜物种结构,导致有效吸附位点的减少,影响其对噻吩的吸附能力。     

Nanocatalyst structure as a template to define chirality of nascent single-walled carbon nanotubes

Chirality is a crucial factor in a single-walled carbon nanotube (SWCNT) because it determines its optical and electronic properties. A chiral angle spanning from 0° to 30° results from twisting of the graphene sheet conforming the nanotube wall and is equivalently expressed by chiral indexes (n,m). However, lack of chirality control during SWCNT synthesis is an obstacle for a widespread use of these materials. Here we use first-principles density functional theory (DFT) and classical molecular dynamics (MD) simulations to propose and illustrate basic concepts supporting that the nanocatalyst structure may act as a template to control the chirality during nanotube synthesis. DFT optimizations of metal cluster (Co and Cu)∕cap systems for caps of various chiralities are used to show that an inverse template effect from the nascent carbon nanostructure over the catalyst may exist in floating catalysts; such effect determines a negligible chirality control. Classical MD simulations are used to investigate the influence of a strongly interacting substrate on the structure of a metal nanocatalyst and illustrate how such interaction may help preserve catalyst crystallinity. Finally, DFT optimizations of carbon structures on stepped (211) and (321) cobalt surfaces are used to demonstrate the template effect imparted by the nanocatalyst surface on the growing carbon structure at early stages of nucleation. It is found that depending on the step structure and type of building block (short chains, single atoms, or hexagonal rings), thermodynamics favor armchair or zigzag termination, which provides guidelines for a chirality controlled process based on tuning the catalyst structure and the type of precursor gas.     

Cu(II)-alkyl chlorocomplexes: stable compounds or transients? DFT prediction of their structure and EPR parameters

DFT calculations were used for studying the structure and reactivity of organocuprates(II) usually considered as intermediates with very weak Cu-C bond. It was found that calculated principal g-tensor values of model compounds RCu(II)Cl(2(-)) are similar to the experimentally found values for organocopper product of photolysis of quaternary ammonium tetrachlorocuprates. The calculations confirm that the most of organocuprates(II) could be stable at ambient conditions, and short lifetimes of organocuprates(II) in solutions or soft matrices are caused by their high reactivity in various bimolecular processes; the rate of those may be close to the rate of diffusion controlled reactions. The charges, spin densities, and d-orbital populations of the Cu atom in them are typical for bivalent copper complexes. Natural bond orbital analysis of organochlorocuprates(II) confirms the formation of polar σ-bond between copper and carbon atoms.     

Lattice mismatch induced nonlinear growth of graphene

As a two-dimensional material, graphene can be obtained via epitaxial growth on a suitable substrate. Recently, an interesting nonlinear behavior of graphene growth has been observed on some metal surfaces, but the underlying mechanism is still elusive. Taking the Ir(111) surface as an example, we perform a mechanistic study on graphene growth using a combined approach of first-principles calculations and kinetic Monte Carlo (kMC) simulations. Small carbon clusters on the terrace or at step sites are studied first. Then, we investigate how these small carbon species are attached to graphene edges. Generally, attachment of carbon atoms is thermodynamically favorable. However, due to substrate effect, there are also some edge sites where graphene growth must proceed via cluster attachment. The overall growth rate is determined by these cluster attachment processes, which have a much lower chance of happening compared to the monomer attachment. On the basis of such an inhomogeneous growth picture, kMC simulations are performed by separating different time scales, and the experimentally found quintic-like behavior is well reproduced. Different nonlinear growth behaviors are predicted for different graphene orientations, which is consistent with previous experiments. Inhomogeneity induced by lattice mismatch revealed in this study is expected to be a universal phenomenon and will play an important role in the growth of many other heteroepitaxial systems.     

Density-functional tight-binding for phosphine-stabilized nanoscale gold clusters

We report a parameterization of the second-order density-functional tight-binding (DFTB2) method for the quantum chemical simulation of phosphine-ligated nanoscale gold clusters, metalloids, and gold surfaces. Our parameterization extends the previously released DFTB2 “auorg” parameter set by connecting it to the electronic parameter of phosphorus in the “mio” parameter set. Although this connection could technically simply be accomplished by creating only the required additional Au–P repulsive potential, we found that the Au 6p and P 3d virtual atomic orbital energy levels exert a strong influence on the overall performance of the combined parameter set. Our optimized parameters are validated against density functional theory (DFT) geometries, ligand binding and cluster isomerization energies, ligand dissociation potential energy curves, and molecular orbital energies for relevant phosphine-ligated Au_n clusters (n = 2–70), as well as selected experimental X-ray structures from the Cambridge Structural Database. In addition, we validate DFTB simulated far-IR spectra for several phosphine- and thiolate-ligated gold clusters against experimental and DFT spectra. The transferability of the parameter set is evaluated using DFT and DFTB potential energy surfaces resulting from the chemisorption of a PH₃ molecule on the gold (111) surface. To demonstrate the potential of the DFTB method for quantum chemical simulations of metalloid gold clusters that are challenging for traditional DFT calculations, we report the predicted molecular geometry, electronic structure, ligand binding energy, and IR spectrum of Au₁₀₈S₂₄(PPh₃)₁₆.

We report a parameterization of the density-functional tight-binding (DFTB) method for the accurate prediction of molecular, electronic and vibrational structure of phosphine-ligated nanoscale gold clusters, metalloids, and gold surfaces.     

<Emphasis Type="Italic">Ab initio</Emphasis> calculations of chemical bond parameters and the band structure of a two-dimensional system: Graphene/MnO(001)

The results of the DFT studies of the band structure, Fermi surface, and chemical bond in ultrathin graphene/MnO(001) and MnO(001) films are presented, and the features of the interatomic interactions at the initial stage of the growth of graphene islands on the manganese oxide surface are considered. The features of spin state in the valence band and at the Fermi level in these systems are discussed. The magnetic moment on the Mn atom is estimated, and the effect of spin polarization on oxygen and carbon atoms is found. Their nature is discussedBased on the structural energy calculations of 2D graphene/MnO(001) and 2D MnO(001), the stability of the systems is established, and the chemical bond energy is determined.     

Monolayer structure of tetracene on Cu (100) surface: parallel geometry

The geometrical arrangement of tetracene on Cu (100) surface at monolayer coverage is studied by using scanning tunneling microscopy measurement and density functional theory (DFT) calculations. Tetracene molecule is found to be oriented with its molecular plane parallel to the substrate surface, and no perpendicular geometry is observed at this coverage. The molecule is aligned either in the [011] or [011] direction due to the fourfold symmetry of the Cu (100) surface. DFT calculations show that the molecule with the "flat-lying" mode has larger adsorption energy than that with the "upright standing" mode, indicating that the former is the more stable structure. With the flat-lying geometry, the carbon atoms prefer to be placed between surface Cu atoms. The molecular center prefers to be located at the bridge site between two nearest surface Cu atoms.     

Magic carbon clusters in the chemical vapor deposition growth of graphene

Ground-state structures of supported C clusters, C(N) (N = 16, ..., 26), on four selected transition metal surfaces [Rh(111), Ru(0001), Ni(111), and Cu(111)] are systematically explored by ab initio calculations. It is found that the core-shell structured C(21), which is a fraction of C(60) possessing three isolated pentagons and C(3v) symmetry, is a very stable magic cluster on all these metal surfaces. Comparison with experimental scanning tunneling microscopy images, dI/dV curves, and cluster heights proves that C(21) is the experimentally observed dominating C precursor in graphene chemical vapor deposition (CVD) growth. The exceptional stability of the C(21) cluster is attributed to its high symmetry, core-shell geometry, and strong binding between edge C atoms and the metal surfaces. Besides, the high barrier of two C(21) clusters' dimerization explains its temperature-dependent behavior in graphene CVD growth.     

Ab initio calculations of the reaction pathways for methane decomposition over the Cu (111) surface

Growth of large-area, few-layer graphene has been reported recently through the catalytic decomposition of methane (CH(4)) over a Cu surface at high temperature. In this study, we used ab initio calculations to investigate the minimum energy pathways of successive dehydrogenation reactions of CH(4) over the Cu (111) surface. The geometries and energies of all the reaction intermediates and transition states were identified using the climbing image nudged elastic band method. The activation barriers for CH(4) decomposition over this Cu surface are much lower than those in the gas phase; furthermore, analysis of electron density differences revealed significant degrees of charge transfer between the adsorbates and the Cu atoms along the reaction path; these features reveal the role of Cu as the catalytic material for graphene growth. All the dehydrogenation reactions are endothermic, except for carbon dimer (C(2)) formation, which is, therefore, the most critical step for subsequent graphene growth, in particular, on Cu (111) surface.     

Communication: Stable carbon nanoarches in the initial stages of epitaxial growth of graphene on Cu(111)

To fully exploit the device potential of graphene, reliable production of large-area, high-quality samples is required. Epitaxial growth on metal substrates have shown promise in this regard, but further improvement would be facilitated by a more complete understanding of the atomistic processes involved in the early growth stages. Using first-principles calculations within density functional theory, we have investigated the energetics and kinetics of graphene nucleation and growth on a Cu(111) surface. Our calculations have revealed an energetic preference for the formation of stable one-dimensional carbon nanoarches consisting of 3-13 atoms when compared to two-dimensional compact islands of equal sizes. We also estimate the critical cluster size that marks the transition from nanoarch dominance to island dominance in the growth sequence. Our findings may provide the structural link between nucleated carbon dimers and larger carbon nanodomes, and are expected to stimulate future experimental efforts.     

The Structure and Stability of Magic Carbon Clusters Observed in Graphene Chemical Vapor Deposition Growth on Ru(0001) and Rh(111) Surfaces

To improve the atomically controlled growth of graphene by chemical vapor deposition (CVD), understanding the evolution from various carbon species to a graphene nucleus on various catalyst surfaces is essential. Experimentally, an ultrastable carbon cluster on Ru(0001) and Rh(111) surfaces was observed, while its structure and formation process were still under debate. Using ab initio calculations and kinetic analyses, we disclose a specific type of carbon cluster, composed of a C₂₁ core and a few dangling C atoms, which is exceptional stable in a size range from 21 to 27 C atoms. The most stable one of them, an isomer of C₂₄ characterized by three dangling C atoms attached to the C₂₁ core (denoted as C₂₁‐3C), is the most promising candidate of the experimental observation. The ultrastability of C₂₁‐3C originates from both the stable core and the appropriate passivation of the dangling carbon atoms by the catalyst surface.