C D : a computer program to generate 1D, 2D and 3D grids of functions dependent on the molecular ab initio electron density

A program to compute many functions dependent on the electron density ρ(r) from the results of ab initio molecular calculations is presented. The program allows the generation of different one-, two-, and three-dimensional grids for further graphical representation or numerical analysis. Other options like extracting separate atom contributions to the function computed or locating maximum and minimum values are also implemented. A number of illustrative applications regarding different ρ(r)-dependent functions are presented and the performance and portability of the program is discussed.     

MgO films grown on yttria-stabilized zirconia by molecular beam epitaxy

MgO films were grown on (0 0 1) yttria-stabilized zirconia (YSZ) substrates by molecular beam epitaxy (MBE). The crystalline structures of these films were investigated using X-ray diffraction and transmission electron microscopy. Growth temperature was varied from 350 to 550 °C, with crystalline quality being improved at higher temperatures. The MgO films had a domain structure: (1 1 1)[1 1 2¯]_MgO(0 0 1)[1 0 0]_YSZ with four twin variants related by a 90° in-plane rotation about the [1 1 1]_MgO axis. The observed epitaxial orientation was compared to previous reports of films grown by pulsed laser deposition and sputtering and explained as resulting in the lowest interface energy.     

Atomistic phenomena in dense fluid shock waves

The shock structure problem is one of the classical problems of fluid mechanics and at least for non-reacting dilute gases it has been considered essentially solved. Here we present a few recent findings, to show that this is not the case. There are still new physical effects to be discovered provided that the numerical technique is general enough to not rule them out a priori. While the results have been obtained for dense fluids, some of the effects might also be observable for shocks in dilute gases.     

卤代对苝醌类光敏剂光敏活性影响的量子化学研究

以一种新型的醌类光敏剂—菌生素 (HMB)为模型化合物 ,利用量子化学从头算HF/ 6 31G和含时密度泛函TD B3LYP/ 6 31G方法计算研究了卤代作用对醌类光敏剂分子性质和光敏活性的影响 .结果发现 ,卤代作用降低了HMB母体的HOMO和LUMO能级 (EHOMO和ELUMO)及其差值△E ,随卤族元素从上到下 ,EHOMO和ELUMO呈增大趋势 ,而△E呈减小趋势 ,使得分子激发光谱有不同程度的红移 ,在增加PQP三重态量子产率的同时降低了分子最低三重激发态的能量 ,两种作用相互抵消 ,使得卤代物的 1O2 量子产率与母体相似 ;增加了分子的绝热电子亲合势 ,使分子产生O2 -的能力下降 ;卤原子的引入 ,降低了HMB母体分子内氢键的强度 ,同时使顺式构型的分子内质子传递势垒增大 ,而使反式构型势垒减小 .     

Kinetics and equilibrium of small metallic clusters: Ab initio confinement molecular dynamics study of 4

The ab initio molecular dynamics (AIMD) [1] is combined with the heuristic, successive confinement method of surveying a potential energy surface (PES) [2], thereby offering a framework for the simulation study of kinetics and equilibrium properties of metallic clusters. This approach is applied to the study of Au₄, a cluster possessing a simple but specific PES, which consists of very shallow and deep basins and due to this presents a challenge to the conventional AIMD methods. Among other things, the probabilities of the transitions between isomers have been found, and on this basis, both the time-dependent and equilibrium populations of the isomers have been calculated for the conditions typical of the NeNePo experiments [3] in the femtosecond pump-probe spectroscopy.     

Rigid nitronyl-nitroxide-labelled anchoring molecules: syntheses, structural and magnetic investigations

The syntheses of two rigid organic molecular rods bearing a nitronyl-nitroxide radical and a terminal nitrogen-based functionality like a pyridine or a cyano group are reported. Both new paramagnetic molecules are fully characterized, including crystal structure analysis. Furthermore their magnetic behaviours in the crystalline state are investigated and their spin concentration corroborate their excellent purity. While the pyridine functionalized rod is synthesized by converting the corresponding benzaldehyde to the phenyl-nitronyl-nitroxide radical, the synthesis of the cyano functionalized rod demonstrates the accessibility of highly sophisticated spin-labelled molecules via cross-coupling reaction with a meta-iodo-phenyl-nitronyl-nitroxide moiety.     

Dispersion analysis of the least‐squares finite‐element shallow‐water system

The frequency or dispersion relation for the least‐squares mixed formulation of the shallow‐water equations is analysed. We consider the use of different approximation spaces corresponding to co‐located and staggered meshes, respectively. The study includes the effect of Coriolis, and the dispersion properties are compared analytically and graphically with those of the mixed Galerkin formulation. Numerical solutions of a test problem to simulate slow Rossby modes illustrate the theoretical results. Copyright © 2003 John Wiley & Sons, Ltd.     

Low-energy ion-assisted control of interfacial structures in metallic multilayers

A molecular dynamics method has been used to simulate the argon ion-assisted deposition of Cu/Co/Cu multilayers and to explore ion beam assistance strategies that can be used during or after the growth of each layer to control interfacial structures. A low-argon ion energy of 5–10 eV was found to minimize a combination of interfacial roughness and interlayer mixing (alloying) during the ion-assisted deposition of multilayers. However, complete flattening with simultaneous ion assistance could not be achieved without some mixing between the layers when a constant ion energy approach was used. It was found that multilayers with lower interfacial roughness and intermixing could be grown either by modulating the ion energy during the growth of each metal layer or by utilizing ion assistance only after the completion of each layers deposition. In these latter approaches, relatively high-energy ions could be used since the interface is buried and less susceptible to intermixing. The interlayer mixing dependence upon the thickness of the over layer has been determined as a function of ion energy.     

A new approach to potential flow around axisymmetric bodies

A new method is introduced to solve potential flow problems around axisymmetric bodies. The approach relies on expressing the infinite series expansion of the Laplace equation solution in terms of a finite sum which preserves the Laplace solution for the potential function under a Neumann-type boundary condition. Then the coefficients of the finite sum are calculated in a least squares approximation sense using the Gram-Schmidt orthonormalization method. Sample benchmark problems are presented and discussed in some detail. The solutions are accurate and converged faster when a rather small number of terms were used. The method is simple and can be easily programmed.     

Least‐squares finite element method on adaptive grid for PDEs with shocks

We apply the least‐squares finite element method with adaptive grid to nonlinear time‐dependent PDEs with shocks. The least‐squares finite element method is also used in applying the deformation method to generate the adaptive moving grids. The effectiveness of this method is demonstrated by solving a Burgers' equation with shocks. Computational results on uniform grids and adaptive grids are compared for the purpose of evaluation. The results show that the adaptive grids can capture the shock more sharply with significantly less computational time. For moving shock, the adaptive grid moves correctly with the shock. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006