Researcher profile

Mark R. Pederson

Mark R. Pederson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

physics.chem-ph cond-mat.mtrl-sci physics.comp-ph physics.atm-clus cond-mat.other cond-mat.str-el cond-mat.supr-con physics.atom-ph quant-ph

Trust snapshot

Quick read

Trust 21 - Emerging

7works

0followers

9topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2015arXiv

Fermi orbital self-interaction corrected electronic structure of molecules beyond local density approximation

The correction of the self-interaction error (SIE) that is inherent to all standard density functional theory (DFT) calculations is an object of increasing interest. In this article we apply the very recently developed Fermi-orbital based approach for the self-interaction correction (FOSIC) to a set of different molecular systems. Our study covers systems ranging from simple diatomic to large organic molecules. We focus our analysis on the direct estimation of the ionization potential from orbital eigenvalues.

preprint2015arXiv

Practical and Rigorous reduction of the Many-Electron Quantum-Mechanical Coulomb Problem to O(N^(2/3) Storage

It is tacitly accepted that, for practical basis sets consisting of N functions, solution of the two-electron Coulomb problem in quantum mechanics requires storage of O(N^4) integrals in the small N limit. For localized functions, in the large N limit, or for planewaves, due to closure, the storage can be reduced to O(N^2) integrals. Here, it is shown that the storage can be further reduced to O(N^{2/3}) for separable basis functions. A practical algorithm, that uses standard one-dimensional Gaussian-quadrature sums, is demonstrated. The resulting algorithm allows for the simultaneous storage, or fast reconstruction, of any two-electron coulomb integral required for a many-electron calculation, on each and every processor of massively parallel computers even if such processors have very limited memory and disk space. For example, for calculations involving a basis of 9171 planewaves, the memory required to effectively store all coulomb integrals decreases from 2.8Gbytes to less than 2.4 Mbytes.

preprint2014arXiv

Fermi Orbital Derivatives in Self-Interaction Corrected Density Functional Theory: Applications to Closed Shell Atoms

A recent modification of the Perdew-Zunger self-interaction-correction (SIC) to the density-functional formalism (Pederson, Ruzsinszky, Perdew) has provided a framework for explicitly restoring unitary invariance to the expression for the total energy. The formalism depends upon construction of Lowdin orthonormalized Fermi-orbitals (Luken et al) which parametrically depend on variational quasi-classical electronic positions. Derivatives of these quasi-classical electronic positions, required for efficient minimization of the self-interaction corrected energy, are derived and tested here on atoms. Total energies and ionization energies in closed-shell atoms, where correlation is less important, using the PW92 LDA functional are in very good to excellent agreement with experiment and non-relativistic Quantum-Monte-Carlo (QMC) results.

preprint2009arXiv

Tuning Molecule-Mediated Spin Coupling in Bottom-Up Fabricated Vanadium-TCNE Nanostructures

We have fabricated hybrid magnetic complexes from V atoms and tetracyanoethylene (TCNE) ligands via atomic manipulation with a cryogenic scanning tunneling microscope. Using tunneling spectroscopy we observe spin-polarized molecular orbitals as well as Kondo behavior. For complexes having two V atoms, the Kondo behavior can be quenched for different molecular arrangements, even as the spin-polarized orbitals remain unchanged. This is explained by variable spin-spin (i.e., V-V) ferromagnetic coupling through a single TCNE molecule, as supported by density functional calculations.

preprint2007arXiv

Comparative study of unscreened and screened molecular static linear polarizability in the Hartree-Fock, hybrid-density functional, and density functional models

The sum-over-states (SOS) polarizabilities are calculated within the approximate independent electron theories such as the Hartree-Fock approximation and density functional models using the eigenvalues and orbitals obtained from the self-consistent solution of single-particle equations. The SOS polarizabilities are then compared with those calculated using the finite-field method within three widely used single particle models: (1) the Hartree-Fock (HF) method, (2) the three parameter hybrid generalized gradient approximation (B3LYP), and (3) the parameter free generalized gradient approximation due to Perdew-Burke-Ernzerhof (PBE). The comparison is carried out for polarizabilities of 142 molecules calculated using the 6-311++G(d,p) orbital basis at the geometries optimized at the B3LYP/6-311G** level. The results show that the SOS method almost always overestimate the FF polarizabilities in the PBE and B3LYP models. This trend is reversed in the HF method however exception to the trend do exist. The mean absolute errors (MAE) in the screened (FF) and unscreened (SOS) polarizability are 0.78 Å$^3$ for the HF method, 1.87 Å$^3$ for the B3LYP, and 3.44 Å$^3$ for the PBE-GGA. The Hartree-Fock approximation predicts finite-field polarizabilities that are smaller than those predicted by the PBE and B3LYP models. Finally, a simple scheme is devised to obtain FF quality estimate of polarizability from those obtained using the SOS method by exploiting the observed trends.

preprint2005arXiv

Electron-phonon coupling as an order-one problem

The coupling between electrons and phonons plays important roles in physics, chemistry and biology. However, the accurate calculation of the electron-phonon coupling constants is computationally expensive as it involves solving the Schrodinger equation for O(N) nuclear configurations, where N is the number of nuclei. Herein we show that by considering the forces on the nuclei caused by the addition or subtraction of an arbitrarily small electronic charge one may calculate the electron-phonon coupling constants from O(1) solutions of the Schrodinger equation. We show that Janak's theorem means that this procedure is exact within the density functional formalism. We demonstrate that the O(1) approach produces numerically accurate results by calculating the electron-phonon coupling constants for a series of molecules ranging in size from H_2 to C_60. We use our approach to introduce a computationally fast approximation for the adiabatic ionisation potentials and electron affinities which is shown to be accurate for large molecules. We also show that our approach allows for the calculation isotope effects in O(0) time and discuss the deuteration driven Mott transition in kappa-(BEDT-TTF)_2Cu[N(CN)_2]Br.

preprint2004arXiv

Electronic structure, vibrational stability, infra-red, and Raman spectra of B24N24 cages

We examine the vibrational stability of three candidate structures for the B24N24 cage and report their infra-red (IR) and Raman spectra. The candidate structures considered are a round cage with octahedral O symmetry, a cage with S_4 symmetry that satisfies the isolated square rule, and a cage of S_8 symmetry, which combines the caps of the (4,4) nanotube, and contains two extra squares and octagons. The calculations are performed within density functional theory, at the all electron level, with large basis sets, and within the generalized gradient approximation. The vertical ionization potential (VIP) and static dipole polarizability are also reported. The S_4 and S_8 cages are energetically nearly degenerate and are favored over the O cage which has six extra octagons and squares. The IR and Raman spectra of the three clusters show notable differences providing thereby a way to identify and possibly synthesize the cages.