Errors in "Wavefunctions for the methane molecule"

Two errors in my paper “Wave functions for the methane molecule” [1] are corrected. They concern my f-harmonic approximation to the wave-function in the equilibrium configuration, for which the final expression for the wave function, the energy lowering, and the density function were all in error.

The determination of molecular properties from MULTIMODE with an application to the calculation of Franck-Condon factors for photoionization of CF3 to CF3+

Extensions to the code MULTIMODE to obtain rovibrational wave functions and properties are described. An application of these new capabilities is made to a calculation of the Franck-Condon factors for photoionization of CF3 to CF3+. These calculations make use of a new, full-dimensional ab initio potential energy surface, which is also described here.

Anharmonic vibrational levels of the two cyclic isomers of SiC3

Using coupled-cluster approach full six-dimensional analytic potential energy surfaces for two cyclic SiC3 isomers [C-C transannular bond (I) and Si-C transannular bond (II)] have been generated and used to calculate anharmonic vibrational wave functions. Several strong low-lying anharmonic resonances have been found. In both isomers already some of the fundamental transitions cannot be described within the harmonic approximation. Adiabatic electron affinities and ionization energies have been calculated as well. The Franck-Condon factors for the photodetachment processes c-SiC3-(I)-> c-SiC3(I) and c-SiC3-(II)-> c-SiC3(II) are reported. (c) 2006 American Institute of Physics.

Error analysis of Trefftz-discontinuous Galerkin methods for the time-harmonic Maxwell equations

In this paper, we extend to the time-harmonic Maxwell equations the p-version analysis technique developed in [R. Hiptmair, A. Moiola and I. Perugia, Plane wave discontinuous Galerkin methods for the 2D Helmholtz equation: analysis of the p-version, SIAM J. Numer. Anal., 49 (2011), 264-284] for Trefftz-discontinuous Galerkin approximations of the Helmholtz problem. While error estimates in a mesh-skeleton norm are derived parallel to the Helmholtz case, the derivation of estimates in a mesh-independent norm requires new twists in the duality argument. The particular case where the local Trefftz approximation spaces are built of vector-valued plane wave functions is considered, and convergence rates are derived.

Why the Heyd-Scuseria-Ernzerhof hybrid functional description of VO2 phases is not correct

In contrast with recent claims that the Heyd-Scuseria-Ernzerhof (HSE) screened hybrid functional can provide a good description of the electronic and magnetic structures of VO2 phases [Eyert, Phys. Rev. Lett.PRLTAO0031-900710.1103/PhysRevLett.107.016401 107, 016401 (2011)], we show here that the HSE lowest-energy solutions for both the low-temperature monoclinic (M1) phase and the high-temperature rutile (R) phase, which are obtained upon inclusion of spin polarization, are at odds with experimental observations. For the M1 phase the ground state is (but should not be) magnetic, while the ground state of the R phase, which is also spin polarized, is not (but should be) metallic. The energy difference between the low-temperature and high-temperature phases has strong discrepancies with the experimental latent heat.

CASPT2 Study of the Potential Energy Surface of the HSO(2) System

The importance of the HSO(2) system in atmospheric and combustion chemistry has motivated several works dedicated to the study of associated structures and chemical reactions. Nevertheless controversy still exists in connection with the reaction SH + O(2) -> H + SO(2) and also related to the role of the HSOO isomers in the potential energy surface (PES). Here we report high-level ab initio calculation for the electronic ground state of the HSO(2) system. Energetic, geometric, and frequency properties for the major stationary states of the PES are reported at the same level of calculations:,CASPT2/aug-cc-pV(T+d)Z. This study introduces three new stationary points (two saddle points and one minimum). These structures allow the connection of the skewed HSOOs and the HSO(2) minima defining new reaction paths for SH + O(2) -> H + SO(2) and SH + O(2) -> OH + SO. In addition, the location of the HSOO isomers in the reaction pathways have been clarified.

Surface and Quantum Confinement Effects in ZnO Nanocrystals

ZnO nanocrystals are studied using theoretical calculations based on the density functional theory. The two main effects related to the reduced size of the nanocrystals are investigated: quantum confinement and a large surface:volume ratio. The effects of quantum confinement are studied by saturating the surface dangling bonds of the nanocrystals with hypothetical H atoms. To understand the effects of the surfaces of the nanocrystals, all saturation is removed and the system is relaxed to its minimum energy position. Several different surface motifs are reported, which should be observed experimentally. Spin-polarized calculations are performed in the nonsaturated nanocrystals, leading to different magnetic moments. We propose that this magnetic moment can be responsible for the intrinsic magnetism observed in ZnO nanostructures.

Hyperfine and magnetic properties of Fe-Cu clusters and Fe precipitates embedded in a Cu matrix

Using the first-principles real-space linear muffin-tin orbital method within the atomic sphere approximation (RS-LMTO-ASA) we study hyperfine and local magnetic properties of substituted pure Fe and Fe-Cu clusters in an fcc Cu matrix. Spin and orbital contributions to magnetic moments, hyperfine fields and the Mossbauer isomer shifts at the Fe sites in Fe precipitates and Fe-Cu alloy clusters of sizes up to 60 Fe atoms embedded in the Cu matrix are calculated and the influence of the local environment on these properties is discussed.

The effect of additional Si and SiGe layers on the confinement potential of GeMn diluted ferromagnetic semiconductor

In this work we analyze the spin-polarized charge density distribution in the GeMn diluted ferromagnetic semiconductors (DFS). The calculations are performed within a self-consistent k.p method, in which the exchange correlation effects in the local density approximation, as well as the strain effects due to the lattice mismatch, are taken into account. Our findings show that the extra confinement potential provided by the barriers and the variation of the Mn content in the DFS are responsible for a separation between the different spin charge densities, giving rise to higher mobility spin-polarized currents or high ferromagnetism transition temperatures systems. (c) 2008 Elsevier B.V. All rights reserved.

Transfer reactions in the investigation of light-nuclei nucleosynthesis

Cross sections for the (6)Li(p,gamma)(7)Be, (7)Li(n,gamma)(8)Li (8)Li(n,gamma)(9)Li and (8)Li(p,gamma)(9)Be capture reactions have been investigated in the framework of the potential model. The main ingredients of the potential model are the potentials used to generate the continuum and bound-state wave functions and spectroscopic factors of the corresponding bound systems. The spectroscopic factors for the (7)Li circle times n=(8)Li(gs), (8)Li circle times n=(9)Li(gs) bound systems were obtained from a FR-DWBA analysis of neutron transfer reactions induced by (8)Li radioactive beam on a (9)Be target, while spetroscopic factor for the (8)Li circle times n=(9)Be(gs) bound system were obained from a proton transfer reaction. From the obtained capture reaction cross section, reaction rate for the (8)Li(n,gamma)(9)Li and (8)Li(p,gamma)(9)Be direct neutron and proton capture were determined and compared with other experimental and calculated values.

Empirical analysis of the Lieb-Oxford bound in ions and molecules

Universal properties of the Coulomb interaction energy apply to all many-electron systems. Bounds on the exchange-correlation energy, in particular, are important for the construction of improved density functionals. Here we investigate one such universal property-the Lieb-Oxford lower bound-for ionic and molecular systems. In recent work [J Chem Phys 127, 054106 (2007)], we observed that for atoms and electron liquids this bound may be substantially tightened. Calculations for a few ions and molecules suggested the same tendency, but were not conclusive due to the small number of systems considered. Here we extend that analysis to many different families of ions and molecules, and find that for these, too, the bound can be empirically tightened by a similar margin as for atoms and electron liquids. Tightening the Lieb-Oxford bound will have consequences for the performance of various approximate exchange-correlation functionals. (C) 2008 Wiley Periodicals Inc.

Predicting and characterizing new molecular species at a high-level theoretical approach: HSeCl, HClSe, and SeCl

New molecular species HSeCl, HClSe, and SeCl were investigated at a high level of theory, CCSD(T), with a series of correlation consistent basis sets with extrapolation to the CBS limit. Account has been taken for valence-only and core-valence correlation effects, and of anharmonic effects on the vibrational frequencies. HSeCl is 43.25 kcal mol (1) more stable than HClSe. A barrier (Delta G(#)) of 47.20 kcal mol (1) separates these species. Internuclear distances are generally overestimated by 0.008 angstrom in the valence-only correlation calculations. Inclusion of anharmonicity leads to much improved vibrational frequencies. For SeCl, we estimate Delta H(f) (0 K) = 23.96 and Delta H(f) (298.15 K) = 24.64 kcal mol (1); for HSeCl, we had 4.20 and 4.97 kcal mol (1), respectively. (C) 2010 Elsevier B.V. All rights reserved.

Structural and spectroscopic properties of the diazocarbene radical (CNN) and its ions CNN(+) and CNN(-): a high-level theoretical investigation

The diazocarbene radical, CNN, and the ions CNN(+) and CNN(-) were investigated at a high level of theory. Very accurate structural parameters for the states X(3)Sigma(-) and A(3)Pi of CNN, and X(2)Pi of both CNN(+) and CNN(-) were obtained with the UCCSD(T) method using correlated-consistent basis functions with extrapolations to the complete basis set limit, with valence only and also with all electrons correlated. Harmonic and anharmonic frequencies were obtained for all species and the Renner parameter and average frequencies evaluated for the Pi states. At the UCCSD(T)/CBS(T-5) level of theory, Delta(f)H(0 K) = 138.89 kcal/mol and Delta(f)H(298 K) = 139.65 kcal/mol were obtained for diazocarbene; for the ionization potential and the electron affinity of CNN, 10.969 eV (252.95 kcal/mol), and 1.743 eV (40.19 kcal/mol), respectively, are predicted. Geometry optimization was also carried out with the CASSCF/MRCI/CBS(T-5) approach for the states X(3)Sigma(-) A(3)Pi, and a(1)Delta of CNN, and with the CASSCF/MRSDCI/aug-cc-pVTZ approach for the states b(1)Sigma(+), c(1)Pi, d(1)Sigma(-), and B(3)Sigma(-), and excitation energies (T(e)) evaluated. Vertical energies were calculated for 15 electronic states, thus improving on the accuracy of the five transitions already described, and allowing for a reliable overview of a manifold of other states, which is expected to guide future spectroscopic experiments. This study corroborates the experimental assignment for the vertical transition X (3)Sigma(-) <- E (3)Pi.

Exploring new molecular species on the (1)[H, Se, F] singlet potential energy surface: Energetics, structures, IR spectra, and heats of formation

Structural, energetic, and vibrational properties of new molecular species, HSeF and HFSe, the associated transition state, and dissociation fragments are investigated using a state-of-the-art theoretical approach, CCSD(T)/CBS. HSeF is a normal covalently bonded molecule 38.98 kcal mol (1) more stable than the complex HF-Se, which shows an unusual structure with a central fluorine atom and a bond angle of 101.8 degrees.A barrier (Delta G(#)) of 49.01 kcal mol (1) separates the two species. Vibrational frequencies are also quite distinct. Heats of formation are evaluated for the diatomic fragments and HSeF. Final Delta(f)H values depend on the experimental accuracy of those of Se(g) and H(2)Se. (c) 2009 Elsevier B.V. All rights reserved.

New Relativistic Atomic Natural Orbital Basis Sets for Lanthanide Atoms with Applications to the Ce Diatom and LuF3

New basis sets of the atomic natural orbital (ANO) type have been developed for the lanthanide atoms La-Lu. The ANOs have been obtained from the average density matrix of the ground and lowest excited states of the atom, the positive ions, and the atom in an electric field. Scalar relativistic effects are included through the use of a Douglas-Kroll-Hess Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second-order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of ionization energies and some excitation energies. Computed ionization energies have an accuracy better than 0.1 eV in most cases. Two molecular applications are inluded as illustration: the cerium diatom and the LuF3 molecule. In both cases it is shown that 4f orbitals are not involved in the chemical bond in contrast to an earlier claim for the latter molecule.