Excitation energies from density functional perturbation theory

We consider two perturbative schemes to calculate excitation energies, each employing the Kohn-Sham Hamiltonian as the unperturbed system. Using accurate exchange-correlation potentials generated from essentially exact densities and their exchange components determined by a recently proposed method, we evaluate energy differences between the ground state and excited states in first-order perturbation theory for the Helium, ionized Lithium and Beryllium atoms. It was recently observed that the zeroth-order excitations energies, simply given by the difference of the Kohn-Sham eigenvalues, almost always lie between the singlet and triplet experimental excitations energies, corrected for relativistic and finite nuclear mass effects. The first-order corrections provide about a factor of two improvement in one of the perturbative schemes but not in the other. The excitation energies within perturbation theory are compared to the excitations obtained within $\Delta$SCF and time-dependent density functional theory. We also calculate the excitation energies in perturbation theory using approximate functionals such as the local density approximation and the optimized effective potential method with and without the Colle-Salvetti correlation contribution

On optimality of kernels for approximate Bayesian computation using sequential Monte Carlo

Approximate Bayesian computation (ABC) has gained popularity over the past few years for the analysis of complex models arising in population genetics, epidemiology and system biology. Sequential Monte Carlo (SMC) approaches have become work-horses in ABC. Here we discuss how to construct the perturbation kernels that are required in ABC SMC approaches, in order to construct a sequence of distributions that start out from a suitably defined prior and converge towards the unknown posterior. We derive optimality criteria for different kernels, which are based on the Kullback-Leibler divergence between a distribution and the distribution of the perturbed particles. We will show that for many complicated posterior distributions, locally adapted kernels tend to show the best performance. We find that the added moderate cost of adapting kernel functions is easily regained in terms of the higher acceptance rate. We demonstrate the computational efficiency gains in a range of toy examples which illustrate some of the challenges faced in real-world applications of ABC, before turning to two demanding parameter inference problems in molecular biology, which highlight the huge increases in efficiency that can be gained from choice of optimal kernels. We conclude with a general discussion of the rational choice of perturbation kernels in ABC SMC settings

Impact of Electron-Electron Cusp on Configuration Interaction Energies

The effect of the electron-electron cusp on the convergence of configuration interaction (CI) wave functions is examined. By analogy with the pseudopotential approach for electron-ion interactions, an effective electron-electron interaction is developed which closely reproduces the scattering of the Coulomb interaction but is smooth and finite at zero electron-electron separation. The exact many-electron wave function for this smooth effective interaction has no cusp at zero electron-electron separation. We perform CI and quantum Monte Carlo calculations for He and Be atoms, both with the Coulomb electron-electron interaction and with the smooth effective electron-electron interaction. We find that convergence of the CI expansion of the wave function for the smooth electron-electron interaction is not significantly improved compared with that for the divergent Coulomb interaction for energy differences on the order of 1 mHartree. This shows that, contrary to popular belief, description of the electron-electron cusp is not a limiting factor, to within chemical accuracy, for CI calculations.Comment: 11 pages, 6 figures, 3 tables, LaTeX209, submitted to The Journal of Chemical Physic

Alleviation of the Fermion-sign problem by optimization of many-body wave functions

We present a simple, robust and highly efficient method for optimizing all parameters of many-body wave functions in quantum Monte Carlo calculations, applicable to continuum systems and lattice models. Based on a strong zero-variance principle, diagonalization of the Hamiltonian matrix in the space spanned by the wav e function and its derivatives determines the optimal parameters. It systematically reduces the fixed-node error, as demonstrated by the calculation of the binding energy of the small but challenging C$_2$ molecule to the experimental accuracy of 0.02 eV

Simultaneous Contralateral Vestibular Schwannoma and Middle Ear Paraganglioma Tumor

To the best of our knowledge, only 2 cases of a simultaneous contralateral vestibular schwannoma (VS) and middle ear paraganglioma (MEP) have previously been reported in literature. We report the third case observed in a 43-year-old male, who presented with an 11-year history of right-sided hearing loss and a 1-year history of left-sided pulsatile tinnitus. A magnetic resonance imaging (MRI) showed a VS on the right side and computer tomography (CT) identified a Fisch type A1 paraganglioma on the left side. The VS was treated using a translabyrinthine approach and the MEP was kept under radiological observation for 1 year. Due to the growth of the MEP (Fisch type A2), it was treated with excision via a retroauricular approach. Our case was very challenging because there was a different and important pathology on each side, both carrying a risk of deafness as a consequence of the disease and/or the treatments

Early assessment of vestibular function after unilateral cochlear implant surgery

Introduction : Cochlear implantation (CI) has been reported to negatively effect on the vestibular function. The study of the vestibular function has variably been conducted by different types of diagnostic tools. The combined use of modern, rapidly performable diagnostic tools could reveal useful for standardizing the evaluation protocol. Methods: In a group of 28 subjects undergoing CI, the video Head Impulse Test (vHIT), the cervical Vestibular Evoked Myogenic Potentials (cVEMPS) and the short-form of Dizziness Handicap Inventory (DHI) questionnaire were investigated pre-operatively and post-operatively (implant on and off) in both the implanted and the contralateral, non-implanted ear. All surgeries were performed with a round window approach (RWA), except for three otosclerosis cases were the extended RWA (eRWA) was used. Results: The vHIT of the lateral semicircular canal showed a pre-operative vestibular involvement in nearly 50% of the cases, whilst the three canals were contemporarily affected in only 14% of them. In all the hypo-functional subjects, cVEMPs were absent. A low VOR gain in all the investigated SSCC was found in 4 subjects (14%). In those subjects, (21.7%) in whom cVEMPs were pre-operatively present and normal in the operated side, absence of response was post-operatives recorded. Discussion/Conclusion: The vestibular protocol applied for the study showed to be appropriate for distinguishing between the CI operated and the non-operated ear. In this regard, cVEMPs showed to be more sensitive than vHIT for revealing a vestibular sufferance after CI, although without statistical significance. Finally, the use of the RWA surgery was apparently not avoiding signs of vestibular impairment to occur

A heuristic framework for the bi-objective enhanced index tracking problem

The index tracking problem is the problem of determining a portfolio of assets whose performance replicates, as closely as possible, that of a financial market index chosen as benchmark. In the enhanced index tracking problem the portfolio is expected to outperform the benchmark with minimal additional risk. In this paper, we study the bi-objective enhanced index tracking problem where two competing objectives, i.e., the expected excess return of the portfolio over the benchmark and the tracking error, are taken into consideration. A bi-objective Mixed Integer Linear Programming formulation for the problem is proposed. Computational results on a set of benchmark instances are given, along with a detailed out-of-sample analysis of the performance of the optimal portfolios selected by the proposed model. Then, a heuristic procedure is designed to build an approximation of the set of Pareto optimal solutions. We test the proposed procedure on a reference set of Pareto optimal solutions. Computational results show that the procedure is significantly faster than the exact computation and provides an extremely accurate approximation

Energy and variance optimization of many body wave functions

We present a simple, robust and efficient method for varying the parameters in a many-body wave function to optimize the expectation value of the energy. The effectiveness of the method is demonstrated by optimizing the parameters in flexible Jastrow factors, that include 3-body electron-electron-nucleus correlation terms, for the NO$_2$ and decapentaene (C$_{10}$H$_{12}$) molecules. The basic idea is to add terms to the straightforward expression for the Hessian that are zero when the integrals are performed exactly, but that cancel much of the statistical fluctuations for a finite Monte Carlo sample. The method is compared to what is currently the most popular method for optimizing many-body wave functions, namely minimization of the variance of the local energy. The most efficient wave function is obtained by optimizing a linear combination of the energy and the variance.Comment: 4 pages, 4 figures, minor corrections of inexact statements, missing

Sc substitution for Mg in MgB2: effects on Tc and Kohn anomaly

Here we report synthesis and characterization of Mg_{1-x}Sc_{x}B_{2} (0.12T_{c}>6 K. We find that the Sc doping moves the chemical potential through the 2D/3D electronic topological transition (ETT) in the sigma band where the ``shape resonance" of interband pairing occurs. In the 3D regime beyond the ETT we observe a hardening of the E_{2g} Raman mode with a significant line-width narrowing due to suppression of the Kohn anomaly over the range 0<q<2k_{F}.Comment: 8 pages, 4 EPS figures, to be published in Phys. Rev.

Correlated sampling in quantum Monte Carlo: a route to forces

In order to find the equilibrium geometries of molecules and solids and to perform ab initio molecular dynamics, it is necessary to calculate the forces on the nuclei. We present a correlated sampling method to efficiently calculate numerical forces and potential energy surfaces in diffusion Monte Carlo. It employs a novel coordinate transformation, earlier used in variational Monte Carlo, to greatly reduce the statistical error. Results are presented for first-row diatomic molecules.Comment: 5 pages, 2 postscript figure