Relativistic molecular calculations of superheavy molecules

Relativistic molecular calculations within the Dirac-Slater scheme have been used in a study of the electronic structure of 6d-metal superheavy hexafluorides. The theoretical results are compared with calculations and measurements of the homolog 4d- and 5d-metal hexafluorides. Large spin-orbit splitting dominates the electronic structure and even has the same order of magnitude as the crystal-field splitting for the valence electrons for the superheavy molecules. Ionization energies have been calculated using a transition state procedure.

Interpretation of noncharacteristic M X-rays in heavy colliding systems by selfconsistent relativistic molecular calculations

The result of the first calculation of a self-consistent relativistic many electron correlation diagram ever done (for the system Au - I) leads to a good agreement of the spectral shape and position of the observed noncharacteristic X-rays within the quasi adiabatic model.

Self-consistent relativistic molecular calculations of superheavy molecules: (_110 X)F_6

Using new relativistic molecular calculations within the Dirac-Slater scheme it is now feasible to study theoretically molecules containing superheavy elements. This opens a new era for the prediction of the physics and chemistry of superheavy elements. As an example we present the results for (_110 X) F_6, where it is shown that relativistic effects are nearly of the same order of magnitude as the crystal-field splitting.

Non-relativistic and relativistic molecular calculations for the chalcogen hexafluorides: SF_6, SeF_6, TeF_6, PoF_6

Non-relativistic Hartree-Fock-Slater and relativistic Dirac-Slater self-consistent orbital models are applied for the analysis of the electronic structure of the chalcogen hexafluorides: SF_6, SeF_6, TeF_6 and PoF_6. The molecular eigenfunctions and eigenvalues are generated using the discrete variational method (DVM) with numerical basis functions. The results obtained for SF_6 are compared with other ab initio calculations. Information about relativistic level shifts and spin-orbit splitting has been obtained by comparison between the non-relativistic and relativistic results.

Ab Initio correlated all electron dirac-fock calculations for Eka-Francium Fluoride (E119F)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Relativistic LCAO with Minimax Principle and New Balanced Basis Sets

Relativistic density functional theory is widely applied in molecular calculations with heavy atoms, where relativistic and correlation effects are on the same footing. Variational stability of the Dirac Hamiltonian is a very important field of research from the beginning of relativistic molecular calculations on, among efforts for accuracy, efficiency, and density functional formulation, etc. Approximations of one- or two-component methods and searching for suitable basis sets are two major means for good projection power against the negative continuum. The minimax two-component spinor linear combination of atomic orbitals (LCAO) is applied in the present work for both light and super-heavy one-electron systems, providing good approximations in the whole energy spectrum, being close to the benchmark minimax finite element method (FEM) values and without spurious and contaminated states, in contrast to the presence of these artifacts in the traditional four-component spinor LCAO. The variational stability assures that minimax LCAO is bounded from below. New balanced basis sets, kinetic and potential defect balanced (TVDB), following the minimax idea, are applied with the Dirac Hamiltonian. Its performance in the same super-heavy one-electron quasi-molecules shows also very good projection capability against variational collapse, as the minimax LCAO is taken as the best projection to compare with. The TVDB method has twice as many basis coefficients as four-component spinor LCAO, which becomes now linear and overcomes the disadvantage of great time-consumption in the minimax method. The calculation with both the TVDB method and the traditional LCAO method for the dimers with elements in group 11 of the periodic table investigates their difference. New bigger basis sets are constructed than in previous research, achieving high accuracy within the functionals involved. Their difference in total energy is much smaller than the basis incompleteness error, showing that the traditional four-spinor LCAO keeps enough projection power from the numerical atomic orbitals and is suitable in research on relativistic quantum chemistry. In scattering investigations for the same comparison purpose, the failure of the traditional LCAO method of providing a stable spectrum with increasing size of basis sets is contrasted to the TVDB method, which contains no spurious states already without pre-orthogonalization of basis sets. Keeping the same conditions including the accuracy of matrix elements shows that the variational instability prevails over the linear dependence of the basis sets. The success of the TVDB method manifests its capability not only in relativistic quantum chemistry but also for scattering and under the influence of strong external electronic and magnetic fields. The good accuracy in total energy with large basis sets and the good projection property encourage wider research on different molecules, with better functionals, and on small effects.

Unsaturated σ-hydrocarbyl transition-metal complexes. Part 4. Crystal and molecular structure of trans-chlorobis(diethylphenylphosphine)-(phenylethynyl)platinum(II) and comments on the relative trans influence of various carbon ligands

The molecular structure of trans-[PtCl(CCPh)(PEt2Ph)2] has been determined by X-ray diffraction methods. The crystals are monoclinic, space group P21, with a= 12.359(3), b= 13.015(3), c= 9.031(2)Å, β= 101.65(2)°, and Z= 2. The structure has been solved by the heavy-atom method and refined by full-matrix least squares to R 0.046 for 1 877 diffractometric intensity data. The crystals contain discrete molecules in which the platinum coordination is square planar. The phenylethynyl group is non-linear, with a Pt–CC angle of 163(2)°. Selected bond lengths are Pt–Cl 2.407(5) and Pt–C 1.98(2)Å. The structural trans influences of CCPh, CHCH2, and CH2SiMe3 ligands in platinum(II) complexes are compared; there is only a small dependence on hybridization at the ligating carbon atom.

The role of carbon impurities on the Si(001)-c(4 x 4) surface reconstruction: Theoretical calculations

In this work we employ the state-of-the-art pseudopotential method, within a generalized gradient approximation to the density functional theory, combined with a recently developed method for the calculation of HREELS spectra to study a series of different proposed models for carbon incorporation on the silicon (001) surface. A fully discussion on the geometry, energetics and specially the comparison between experimental and theoretical STM images and electron energy loss spectra indicate that the Si(100)-c(4 x 4) is probably induced by Si-C surface dinners, in agreement with recent experimental findings. (C) 2009 Elsevier B.V. All rights reserved.

Ab Initio correlated all electron Dirac-Fock calculations for eka-francium fluoride (E119F)

Resultados obtidos por cálculos Dirac-Fock correlacionados de 4 componentes para o fluoreto do elemento E119 (Eka-Frâncio) com base estável e precisa, livre de prolapso variacional, são reportados neste trabalho. No nível CCSD(T), a distância de equilíbrio R_e, frequência harmônica ω_e e energia de dissociação D_e são 2,432 Å, 354,97 cm_-1 e 116,92 kcal mol_-1, respectivamente. Também são reportados base livre de prolapso variacional de 4 componentes para o elemento 119, uma curva analítica de energia potencial precisa e o espectro vibracional a partir dos dados obtidos no nível CCSD(T). Nossos resultados sugerem que a molécula E119F deva ser menos iônica que seus fluoretos alcalinos homólogos mais leves, em contraste com o senso químico comum baseado nas propriedades periódicas - era de se esperar nesta molécula a ligação química mais iônica possível. Também encontramos que a correção do tipo modelo de carga para negligenciar as integrais do tipo SS resulta em erros insignificantes e acelera os cálculos cerca de 3 vezes no nível CCSD(T) e cerca de 4 vezes no nível DFT/B3LYP.

Beschreibung von Oberflächenphänomenen durch relativistische Clusterrechnungen unter Verwendung eines Einbettungsverfahrens

For the theoretical investigation of local phenomena (adsorption at surfaces, defects or impurities within a crystal, etc.) one can assume that the effects caused by the local disturbance are only limited to the neighbouring particles. With this model, that is well-known as cluster-approximation, an infinite system can be simulated by a much smaller segment of the surface (Cluster). The size of this segment varies strongly for different systems. Calculations to the convergence of bond distance and binding energy of an adsorbed aluminum atom on an Al(100)-surface showed that more than 100 atoms are necessary to get a sufficient description of surface properties. However with a full-quantummechanical approach these system sizes cannot be calculated because of the effort in computer memory and processor speed. Therefore we developed an embedding procedure for the simulation of surfaces and solids, where the whole system is partitioned in several parts which itsself are treated differently: the internal part (cluster), which is located near the place of the adsorbate, is calculated completely self-consistently and is embedded into an environment, whereas the influence of the environment on the cluster enters as an additional, external potential to the relativistic Kohn-Sham-equations. The basis of the procedure represents the density functional theory. However this means that the choice of the electronic density of the environment constitutes the quality of the embedding procedure. The environment density was modelled in three different ways: atomic densities; of a large prepended calculation without embedding transferred densities; bulk-densities (copied). The embedding procedure was tested on the atomic adsorptions of 'Al on Al(100) and Cu on Cu(100). The result was that if the environment is choices appropriately for the Al-system one needs only 9 embedded atoms to reproduce the results of exact slab-calculations. For the Cu-system first calculations without embedding procedures were accomplished, with the result that already 60 atoms are sufficient as a surface-cluster. Using the embedding procedure the same values with only 25 atoms were obtained. This means a substantial improvement if one takes into consideration that the calculation time increased cubically with the number of atoms. With the embedding method Infinite systems can be treated by molecular methods. Additionally the program code was extended by the possibility to make molecular-dynamic simulations. Now it is possible apart from the past calculations of fixed cores to investigate also structures of small clusters and surfaces. A first application we made with the adsorption of Cu on Cu(100). We calculated the relaxed positions of the atoms that were located close to the adsorption site and afterwards made the full-quantummechanical calculation of this system. We did that procedure for different distances to the surface. Thus a realistic adsorption process could be examined for the first time. It should be remarked that when doing the Cu reference-calculations (without embedding) we begun to parallelize the entire program code. Only because of this aspect the investigations for the 100 atomic Cu surface-clusters were possible. Due to the good efficiency of both the parallelization and the developed embedding procedure we will be able to apply the combination in future. This will help to work on more these areas it will be possible to bring in results of full-relativistic molecular calculations, what will be very interesting especially for the regime of heavy systems.

Fully relativistic calculations on the potential energy surfaces of the lowest 23 states of molecular chlorine

The electronic structure and spectroscopic properties (R(e), omega(e), omega(e)x(e), beta(e), and T(e)) of the ground state and the 22 lowest excited states of chlorine molecule were studied within a four-component relativistic framework using the MOLFDIR program package. The potential energy curves of all possible 23 covalent states were calculated using relativistic complete open shell configuration interaction approach. In addition, four component multireference configuration interaction with single and double excitation calculations were performed in order to infer the effects due to dynamical correlation in vertical excitations. The calculated properties are in good agreement with the available experimental data.

Fully Relativistic 4-Components DFT Investigation on Bonding and Dissociation Energy of HgO

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Many-electron relativistic calculation and interpretation of atomic processes in time dependent heavy-ion scattering

The time dependence of a heavy-ion-atom collision system is solved via a set of coupled channel equations using energy eigenvalues and matrix elements from a self-consistent field relativistic molecular many-electron Dirac-Fock-Slater calculation. Within this independent particle model we give a full many-particle interpretation by performing a small number of single-particle calculations. First results for the P(b) curves for the Ne K-hole excitation for the systems F{^8+} - Ne and F{^6+} - Ne as examples are discussed.

The electronic states of SeF: A reinterpretation of the chemiluminescent emission of the reaction of selenium with fluorine

The low-lying doublet and quartet electronic states of the species SeF correlating with the first dissociation channel are investigated theoretically at a high-level of electronic correlation treatment, namely, the complete active space self-consistent field/multireference single and double excitations configuration interaction (CASSCF/MRSDCI) using a quintuple-zeta quality basis set including a relativistic effective core potential for the selenium atom. Potential energy curves for (Lambda+S) states and the corresponding spectroscopic properties are derived that allows for an unambiguous assignment of the only spectrum known experimentally as due to a spin-forbidden X (2)Pi-a (4)Sigma(-) transition, and not a A (2)Pi-X (2)Pi transition as assumed so far. For the bound excited doublets, yet unknown experimentally, this study is the first theoretical characterization of their spectroscopic properties. Also the spin-orbit coupling constant function for the X (2)Pi state is derived as well as the spin-orbit coupling matrix element between the X (2)Pi and a (4)Sigma(-) states. Dipole moment functions and vibrationally averaged dipole moments show SeF to be a very polar species. An overview of the lowest-lying spin-orbit (Omega) states completes this description. (C) 2010 American Institute of Physics. [doi: 10.1063/1.3426315]