Calculations of bond dissociation energies and dipole moments in energetic materials using density-functional methods

We compare the effectiveness of six exchange/correlation functional combinations (Becke/Lee, Yang and Parr; Becke-3/Lee, Yang and Parr; Becke/Perdew-Wang 91; Becke-3/Perdew-Wang 91; Becke/Perdew 86; Becke-3/Perdew 86) for computing C-N, O-O and N-NO2 dissociation energies and dipole moments of five compounds. The studied compounds are hexabydro-1,3,5-trinitro-1,3,5-triazine (RDX), dimethylnitramine, cyanogen, nitromethane and ozone. The Becke-3/Perdew 86 in conjunction with 6-31G

Density-functional study of lanthanum, ytterbium, and lutetium dimers

La-2, Yb-2, and Lu-2 have been studied by use of the density-functional methods B3LYP, BLYP, B3PW91, BHLYP, BP86, B3P86, MPW1PW91, and PBE1PBE. In these density-functional methods, the exchange functional is from either Becke's three-parameter HF-DFT hybrid exchange functional (B3), pure DFT exchange functional of 1988 (B), a modification of the half-and-half HF/DFT hybrid method (BH), Perdew-Wang 1991 (PW91), or Barone's modified PW91 (MPW1), while the correlation functional is from either Lee, Yang, and Parr (LYP), Perdew-Wang 1991 (PW91), or Perdew 86 (P86). PBE1PBE is the generalized-gradient-approximation exchange-correlation functional of Perdew, Burke, and Ernzerhof. For La-2, the calculated bond distance is in reasonable agreement with the experiment, but the calculated vibrational frequency is underestimated significantly compared with the experiment. Only BP86 and B3P86 have the best performance in reproducing the experimental dissociation energy for La-2. For the van der Waals dimer Yb-2, three functionals, B3LYP, BLYP, and BHLYP have excellent performance in reproducing the spectroscopic constants compared with both the experiment and previous theoretical studies.

Structures and vibrational spectra of C-2 and LaC2+ clusters

C-2 and LaC2+ were studied using Hartree-Fock(HF), B3LYP (Becke 3-paremeter-Lee-Yang-Parr) density functional method, second-order Moller-Plesset perturbation (MP2) and coupled cluster singles and doubles with non-iterative triples(CCSD(T)) methods. The basis set employed was LANL1DZ. Geometries, vibrational frequencies and other quantities were reported. The results showed that for C-2, all the methods performed well for low spin state (singlet), while only HF and B3LYP remained so for high spin state (triplet). For LaC2+, four isomers were presented and fully optimized. The results suggested that linear isomers with C-infinity v and D-infinity h symmetries were predicted to be saddle points on the energy surface for all the methods, while for isomers with C-2 upsilon and C-s symmetries, they were local minima except C-2 upsilon at B3LYP level, and were isoenergetic at HF, MP2 and CCSD(T) levels, near isoenergetic at B3LYP level. From the differences between HOMO and LUMO, it is also known that the isomers with C-2 upsilon and C-s symmetries offer the largest values and therefore correspond to the most stable structure. For La-C bond lengths, B3LYP gives the shortest, the order is B3LYP

Theoretical investigation of LaC2 and LaC2+ clusters

LaC2 (with doublet and quartet states) and LaC2+ (with singlet and triplet states) cluster have been studied by using the B3LYP (Becke three-parameter/Lee-Yang-Parr) density functional method and the HF (Hartree-Fock) method with LANLIDZ basis set. For each cluster, four possible isomers in C-2v, C-s, C-proportional to v and D-proportional to h symmetries have been investigated. The results indicate that structures in C-s symmetry are local minima in all cases and, in most cases (particularly for high spin states), our initial guess in C-s symmetry converges to structures in C-2v symmetry. For the isomers in C-2v, C-proportional to v and D-proportional to h symmetries, local minima were found to be dependent on the method and spin state. The two clusters may also exist as linear chains. The ordering of the binding energies for the isomers in all spin states is C-s similar to C-2v < C-proportional to v < D-proportional to h. The ionization potential of LaC2 is reported as well. (C) 1998 Elsevier Science B.V.

Theoretical investigation of LaC3n+ (n = 0, 1, 2) clusters by density functional theory

LaC3n+ (n = 0, 1, 2) clusters have been studied using B3LYP (Becke 3-parameter-Lee-Yang-Parr) density functional method. The basis set is Dunning/ Huzinaga valence double zeta for carbon and [2s2p2d] for lanthanum, denoted LANL1DZ. Four isomers are presented for each cluster; two of them are edge binding isomers with C-2 upsilon symmetry, the other two are Linear chains with C-infinity upsilon symmetry. Meanwhile, two spin states for each isomer, that is, singlet and triplet for LaC3+, doublet and quartet for LaC3 and LaC32+, respectively, are also considered. Geometries, vibrational frequencies, infrared intensities, and other quantities are reported and discussed. The results indicate that at some spin states; the C-2 upsilon symmetry isomers are the dominant structures, while for the other spin states, linear isomers are energetically favored. (C) 1998 John Wiley & Sons, Inc.

Modification of Bertrand's theorem and extended Runge-Lenz vector

It is shown that for a particle with suitable angular moments in the screened Coulomb potential or isotropic harmonic potential, there still exist closed orbits rather than ellipse, characterized by the conserved aphelion and perihelion vectors, i.e. extended Runge-Lenz vector, which implies a higher dynamical symmetry than the geometrical symmetry O-3. The closeness of a planar orbit implies the radial and angular motional frequencies are commensurable.

The Complete Proof of the Virial Theorem in the Refined TFD Theory for all Electrons of an Atom in a Solid

The complete proof of the virial theorem in refined Thomas-Fermi-Dirac theory for all electrons of an atom in a solid is given.

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.