Ab initio absorption spectrum of NiO combining molecular dynamics with the embedded cluster approach in a discrete reaction field

We developed a procedure that combines three complementary computational methodologies to improve the theoretical description of the electronic structure of nickel oxide. The starting point is a Car-Parrinello molecular dynamics simulation to incorporate vibrorotational degrees of freedom into the material model. By means ofcomplete active space self-consistent field second-order perturbation theory (CASPT2) calculations on embedded clusters extracted from the resulting trajectory, we describe localized spectroscopic phenomena on NiO with an efficient treatment of electron correlation. The inclusion of thermal motion into the theoretical description allowsus to study electronic transitions that, otherwise, would be dipole forbidden in the ideal structure and results in a natural reproduction of the band broadening. Moreover, we improved the embedded cluster model by incorporating self-consistently at the complete active space self-consistent field (CASSCF) level a discrete (or direct) reaction field (DRF) in the cluster surroundings. The DRF approach offers an efficient treatment ofelectric response effects of the crystalline embedding to the electronic transitions localized in the cluster. We offer accurate theoretical estimates of the absorption spectrum and the density of states around the Fermi level of NiO, and a comprehensive explanation of the source of the broadening and the relaxation of the charge transferstates due to the adaptation of the environment

Calculation of Franck-Condon factors including anharmonicity: simulation of the C2 H4+ X̃2B 3u←C2 H4X̃1 Ag band in the photoelectron spectrum of ethylene

Our new simple method for calculating accurate Franck-Condon factors including nondiagonal (i.e., mode-mode) anharmonic coupling is used to simulate the C2H4+X2B 3u←C2H4X̃1 Ag band in the photoelectron spectrum. An improved vibrational basis set truncation algorithm, which permits very efficient computations, is employed. Because the torsional mode is highly anharmonic it is separated from the other modes and treated exactly. All other modes are treated through the second-order perturbation theory. The perturbation-theory corrections are significant and lead to a good agreement with experiment, although the separability assumption for torsion causes the C2 D4 results to be not as good as those for C2 H4. A variational formulation to overcome this circumstance, and deal with large anharmonicities in general, is suggested

The effect of counterpoise correction and relaxation energy term to the internal rotation barriers: application to the BF3 ...NH3 and C2H4 ...SO2 dimers

The relevance of the fragment relaxation energy term and the effect of the basis set superposition error on the geometry of the BF3⋯NH3 and C2H4⋯SO2 van der Waals dimers have been analyzed. Second-order Møller-Plesset perturbation theory calculations with the d95(d,p) basis set have been used to calculate the counterpoise-corrected barrier height for the internal rotations. These barriers have been obtained by relocating the stationary points on the counterpoise-corrected potential energy surface of the processes involved. The fragment relaxation energy can have a large influence on both the intermolecular parameters and barrier height. The counterpoise correction has proved to be important for these systems

Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

Variational calculation of vibrational linear and nonlinear optical properties

A variational approach for reliably calculating vibrational linear and nonlinear optical properties of molecules with large electrical and/or mechanical anharmonicity is introduced. This approach utilizes a self-consistent solution of the vibrational Schrödinger equation for the complete field-dependent potential-energy surface and, then, adds higher-level vibrational correlation corrections as desired. An initial application is made to static properties for three molecules of widely varying anharmonicity using the lowest-level vibrational correlation treatment (i.e., vibrational Møller-Plesset perturbation theory). Our results indicate when the conventional Bishop-Kirtman perturbation method can be expected to break down and when high-level vibrational correlation methods are likely to be required. Future improvements and extensions are discussed

Magnetic moments of the Lambda(1405) and Lambda(1670) resonances

By using techniques of unitarized chiral perturbation theory, where the Lamda(1405) and Lamda(1670) resonances are dynamically generated, we evaluate the magnetic moments of these resonances and their transition magnetic moment. The results obtained here differ appreciably from those obtained with existing quark models. The width for the Lamda(1670)->Lamda(1405)gamma transition is also evaluated, leading to a branching ratio of the order of 210-6.

The infrared behaviour of the static potential in perturbative QCD

The definition of the quark-antiquark static potential is given within an effective field theory framework. The leading infrared divergences of the static singlet potential in perturbation theory are explicitly calculated.

The QCD Potential at O(1/m2): Complete spin-dependent and spin-independent result

Within an effective field theory framework, we obtain an expression, with O(1/m2) accuracy, for the energies of the gluonic excitations between heavy quarks, which holds beyond perturbation theory. For the singlet heavy-quarkantiquark energy, in particular, we also obtain an expression in terms of Wilson loops. This provides, twenty years after the seminal work of Eichten and Feinberg, the first complete expression for the heavy quarkonium potential up to O(1/m2) for pure gluodynamics. Several errors present in the previous literature (also in the work of Eichten and Feinberg) have been corrected. We also briefly discuss the power counting of NRQCD in the nonperturbative regime.

Temperature-induced resonances and Landau damping of collective modes in Bose-Einstein condensed gases in spherical traps

Interaction between collective monopole oscillations of a trapped Bose-Einstein condensate and thermal excitations is investigated by means of perturbation theory. We assume spherical symmetry to calculate the matrix elements by solving the linearized Gross-Pitaevskii equations. We use them to study the resonances of the condensate induced by temperature when an external perturbation of the trapping frequency is applied and to calculate the Landau damping of the oscillations.

Landau damping of transverse quadrupole oscillations of an elongated Bose-Einstein condensate

We have studied the interaction between the low-lying transverse collective oscillations and the thermal excitations of an elongated Bose-Einstein condensate by means of perturbation theory. We consider a cylindrical trapped condensate and calculate the transverse elementary excitations at zero temperature by solving the linearized Gross-Pitaevskii equations in two dimensions (2D). We use them to calculate the matrix elements between the thermal excited states and the quasi-2D collective modes. The Landau damping of transverse collective modes is studied as a function of temperature. At low temperatures, the corresponding damping rate is in agreement with the experimental data for the decay of the transverse quadrupole mode, but it is too small to explain the observed slow decay of the transverse breathing mode. The reason for this discrepancy is discussed.

Microscopic calculations of the excitation spectrum of one 3He impurity in liquid 4He

We calculate the chemical potential ¿0 and the effective mass m*/m3 of one 3He impurity in liquid 4He. First a variational wave function including two- and three-particle dynamical correlations is adopted. Triplet correlations bring the computed values of ¿0 very close to the experimental results. The variational estimate of m*/m3 includes also backflow correlations between the 3He atom and the particles in the medium. Different approximations for the three-particle distribution function give almost the same values for m*/m3. The variational approach underestimates m*/m3 by ~10% at all of the considered densities. Correlated-basis perturbation theory is then used to improve the wave function to include backflow around the particles of the medium. The perturbative series built up with one-phonon states only is summed up to infinite order and gives results very close to the variational ones. All the perturbative diagrams with two independent phonons have then been summed to compute m*/m3. Their contribution depends to some extent on the form used for the three-particle distribution function. When the scaling approximation is adopted, a reasonable agreement with the experimental results is achieved.

Interaccion among systems of finite size in predictive relativistic mechanics I. General framework

We derive a Hamiltonian formulation for the three-dimensional formalism of predictive relativistic mechanics. This Hamiltonian structure is used to derive a set of dynamical equations describing the interaction among systems in perturbation theory.

Local size segregation in polydisperse hard sphere fluids

The structure of polydisperse hard sphere fluids, in the presence of a wall, is studied by the Rosenfeld density functional theory. Within this approach, the local excess free energy depends on only four combinations of the full set of density fields. The case of continuous polydispersity thereby becomes tractable. We predict, generically, an oscillatory size segregation close to the wall, and connect this, by a perturbation theory for narrow distributions, with the reversible work for changing the size of one particle in a monodisperse reference fluid.

General Interaction Picture from Action Principle for Mechanics

In this paper we consider a general action principle for mechanics written by means of the elements of a Lie algebra. We study the physical reasons why we have to choose precisely a Lie algebra to write the action principle. By means of such an action principle we work out the equations of motion and a technique to evaluate perturbations in a general mechanics that is equivalent to a general interaction picture. Classical or quantum mechanics come out as particular cases when we make realizations of the Lie algebra by derivations into the algebra of products of functions or operators, respectively. Later on we develop in particular the applications of the action principle to classical and quantum mechanics, seeing that in this last case it agrees with Schwinger's action principle. The main contribution of this paper is to introduce a perturbation theory and an interaction picture of classical mechanics on the same footing as in quantum mechanics.

Ground- and excited-state properties of M-center oxygen vacancy aggregates in the bulk and the surface of MgO

Aggregates of oxygen vacancies (F centers) represent a particular form of point defects in ionic crystals. In this study we have considered the combination of two oxygen vacancies, the M center, in the bulk and on the surface of MgO by means of cluster model calculations. Both neutral and charged forms of the defect M and M+ have been taken into account. The ground state of the M center is characterized by the presence of two doubly occupied impurity levels in the gap of the material; in M+ centers the highest level is singly occupied. For the ground-state properties we used a gradient corrected density functional theory approach. The dipole-allowed singlet-to-singlet and doublet-to-doublet electronic transitions have been determined by means of explicitly correlated multireference second-order perturbation theory calculations. These have been compared with optical transitions determined with the time-dependent density functional theory formalism. The results show that bulk M and M+ centers give rise to intense absorptions at about 4.4 and 4.0 eV, respectively. Another less intense transition at 1.3 eV has also been found for the M+ center. On the surface the transitions occur at 1.6 eV (M+) and 2 eV (M). The results are compared with recently reported electron energy loss spectroscopy spectra on MgO thin films.