Overview. Post-excitation dynamics in biomimetic photo-catalysts. Strain-induced metal insulator transitions in strongly correlated transition metal oxide materials. Tracking ultrafast charge migration at a molecule-material interface with x-ray spectroscopy. Due in large part to the number of and complicated interactions between the electrons in systems like these, current theoretical methods are often unable to offer reliable predictions in support of the design and analysis of new experiments and technology. Our group specializes in developing new theoretical models, algorithms, and software in order to achieve reliable predictive power over ever larger and more complicated collections of electrons in molecules and materials.
Treating Excited States on an Equal Footing. Due to both the historical priority of ground states and the inherent limitations of the traditional variational principle, electronic structure methods for modeling excited states tend to be more approximate than those for ground states. This imbalance is present both in chemistry, where excitations are most often modeled by linear response theory, and in solid state physics, where band gaps and spectra are most commonly derived from DFT and the many-body perturbation theory of Greens functions. While substantially more robust methods exist in both fields for modeling ground states (e.g. DMRG, QMC, fully relaxed coupled cluster, and so on), applying these methods to excited states is typically not possible without additional approximations, such as the linear response approximation or the use of state-averaged orbital shapes. A major direction of investigation in our group is into methods that exploit alternative variational principles that can treat excited states on the same footing as ground states, thus avoiding the need for different levels of approximation for different states. These approaches are especially important for excitations that frustrate traditional approaches, such as charge transfer excitations, band gaps near a Mott transition, and double excitations.
Developing New Wave Function Approximations. Whether dealing with ground or excited states, current methods for modeling large groups of electrons are either too expensive or too unreliable to provide the desired degree of predictive power in a wide range of molecular and materials systems. These include the multi-metal catalytic cores of many enzymes, various strongly correlated functional materials, and the charge-transfer relays responsible for moving electrons in and out of reaction centers in both natural and artificial photosynthesis. Our group is pursuing new directions in wave function approximation to help meet these challenges, including variational analogues of the coupled cluster expansion, functional forms that can relax themselves in the presence of their own excited state linear response, and wave function stenciling, in which one wave function component detects and deletes the unwanted portions of another. As in all areas of our work, this research requires careful consideration of the intersections between physical approximation and algorithmic efficacy.
Improving Underlying Algorithms. Much of our group’s work is made possible by developing our own optimization algorithms so as to overcome the numerical challenges that arise from the physical approximations we employ. In some cases, we gain access to generalized variational principles though statistical sampling techniques whose random nature invalidates much of standard optimization theory. In other cases, we are faced with optimization problems that require high-order derivatives of energies and wave functions, which we have found can be evaluated efficiently for our purposes via recursive formulations of algorithmic differentiation. Regardless of the specific application, designing our algorithms to take proper advantage of modern supercomputers adds an additional layer of consideration and entertainment.
QMCPACK. As partners in the Center for Predictive Simulation of Functional Materials (CPSFM) we contribute to the QMCPACK software package, an open-source and highly-parallel suite of methods for quantum Monte Carlo simulations in molecules and materials. Two major goals of CPSFM are to make QMC methods more black-box and to reduce their dependence on DFT-based molecular orbitals. We assist with both of these efforts through our work on wave function optimization methods.