Source author record

Paul E. Grabowski

Paul E. Grabowski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

quant-ph physics.atom-ph physics.chem-ph physics.comp-ph astro-ph.SR cond-mat.quant-gas nucl-th physics.plasm-ph

Catalog footprint

What is connected

6works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

preprint2016arXiv

Exact conditions on the temperature dependence of density functionals

Universal exact conditions guided the construction of most ground-state density functional approximations in use today. We derive the relation between the entropy and Mermin free energy density functionals for thermal density functional theory. Both the entropy and sum of kinetic and electron-electron repulsion functionals are shown to be monotonically increasing with temperature, while the Mermin functional is concave downwards. Analogous relations are found for both exchange and correlation. The importance of these conditions is illustrated in two extremes: the Hubbard dimer and the uniform gas.

preprint2016arXiv

Thermal density functional theory: Time-dependent linear response and approximate functionals from the fluctuation-dissipation theorem

The van Leeuwen proof of linear-response time-dependent density functional theory (TDDFT) is generalized to thermal ensembles. This allows generalization to finite temperatures of the Gross-Kohn relation, the exchange-correlation kernel of TDDFT, and fluctuation dissipation theorem for DFT. This produces a natural method for generating new thermal exchange-correlation (XC) approximations.

preprint2014arXiv

A Review of Wave Packet Molecular Dynamics

Warm dense matter systems created in the laboratory are highly dynamical. In such cases electron dynamics is often needed to accurately simulate the evolution and properties of the system. Large systems force one to make simple approximations enabling computationally feasibility. Wave packet molecular dynamics (WPMD) provides a simple framework for simulating time-dependent quantum plasmas. Here, this method is reviewed. The different variants of WPMD are shown and compared and their validity is discussed.

preprint2014arXiv

Quantum critical benchmark for density functional theory

Two electrons at the threshold of ionization represent a severe test case for electronic structure theory. A pseudospectral method yields a very accurate density of the two-electron ion with nuclear charge close to the critical value. Highly accurate energy components and potentials of Kohn-Sham density functional theory are given, as well as a useful parametrization of the critical density. The challenges for density functional approximations and the strength of correlation are also discussed.

preprint2011arXiv

Pseudospectral Calculation of Helium Wave Functions, Expectation Values, and Oscillator Strength

The pseudospectral method is a powerful tool for finding highly precise solutions of Schrödinger's equation for few-electron problems. We extend the method's scope to wave functions with non-zero angular momentum and test it on several challenging problems. One group of tests involves the determination of the nonrelativistic electric dipole oscillator strength for the helium $1^1$S $\to 2^1$P transition. The result achieved, $0.27616499(27)$, is comparable to the best in the literature. Another group of test applications is comprised of well-studied leading order finite nuclear mass and relativistic corrections for the helium ground state. A straightforward computation reaches near state-of-the-art accuracy without requiring the implementation of any special-purpose numerics. All the relevant quantities tested in this paper -- energy eigenvalues, S-state expectation values and bound-bound dipole transitions for S and P states -- converge exponentially with increasing resolution and do so at roughly the same rate. Each individual calculation samples and weights the configuration space wave function uniquely but all behave in a qualitatively similar manner. Quantum mechanical matrix elements are directly and reliably calculable with pseudospectral methods. The technical discussion includes a prescription for choosing coordinates and subdomains to achieve exponential convergence when two-particle Coulomb singularities are present. The prescription does not account for the wave function's non-analytic behavior near the three-particle coalescence which should eventually hinder the rate of the convergence. Nonetheless the effect is small in the sense that ignoring the higher-order coalescence does not appear to affect adversely the accuracy of any of the quantities reported nor the rate at which errors diminish.

preprint2010arXiv

Pseudospectral Calculation of the Wavefunction of Helium and the Negative Hydrogen Ion

We study the numerical solution of the non-relativistic Schrödinger equation for two-electron atoms in ground and excited S-states using pseudospectral (PS) methods of calculation. The calculation achieves convergence rates for the energy, Cauchy error in the wavefunction, and variance in local energy that are exponentially fast for all practical purposes. The method requires three separate subdomains to handle the wavefunction's cusp-like behavior near the two-particle coalescences. The use of three subdomains is essential to maintaining exponential convergence. A comparison of several different treatments of the cusps and the semi-infinite domain suggest that the simplest prescription is sufficient. For many purposes it proves unnecessary to handle the logarithmic behavior near the three-particle coalescence in a special way. The PS method has many virtues: no explicit assumptions need be made about the asymptotic behavior of the wavefunction near cusps or at large distances, the local energy is exactly equal to the calculated global energy at all collocation points, local errors go down everywhere with increasing resolution, the effective basis using Chebyshev polynomials is complete and simple, and the method is easily extensible to other bound states. This study serves as a proof-of-principle of the method for more general two- and possibly three-electron applications.