Source author record

John Parkhill

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

8 published item(s)

preprint2018arXiv

Compressing physical properties of atomic species for improving predictive chemistry

The answers to many unsolved problems lie in the intractable chemical space of molecules and materials. Machine learning techniques are rapidly growing in popularity as a way to compress and explore chemical space efficiently. One of the most important aspects of machine learning techniques is representation through the feature vector, which should contain the most important descriptors necessary to make accurate predictions, not least of which is the atomic species in the molecule or material. In this work we introduce a compressed representation of physical properties for atomic species we call the elemental modes. The elemental modes provide an excellent representation by capturing many of the nuances of the periodic table and the similarity of atomic species. We apply the elemental modes to several different tasks for machine learning algorithms and show that they enable us to make improvements to these tasks even beyond simply achieving higher accuracy predictions.

preprint2016arXiv

Cost-effective description of strong correlation: efficient implementations of the perfect quadruples and perfect hextuples models

Novel implementations based on dense tensor storage are presented for the singlet-reference perfect quadruples (PQ) [Parkhill, Lawler, and Head-Gordon, J. Chem. Phys. 130, 084101 (2009)] and perfect hextuples (PH) [Parkhill and Head-Gordon, J. Chem. Phys. 133, 024103 (2010)] models. The methods are obtained as block decompositions of conventional coupled-cluster theory that are exact for four electrons in four orbitals (PQ) and six electrons in six orbitals (PH), but that can also be applied to much larger systems. PQ and PH have storage requirements that scale as the square, and as the cube of the number of active electrons, respectively, and exhibit quartic scaling of the computational effort for large systems. Applications of the new implementations are presented for full-valence calculations on linear polyenes (C n H n+2 ), which highlight the excellent computational scaling of the present implementations that can routinely handle active spaces of hundreds of electrons. The accuracy of the models is studied in the π space of the polyenes, in hydrogen chains (H 50 ), and in the π space of polyacene molecules. In all cases, the results compare favorably to density matrix renormalization group values. With the novel implementation of PQ, active spaces of 140 electrons in 140 orbitals can be solved in a matter of minutes on a single core workstation, and the relatively low polynomial scaling means that very large systems are also accessible using parallel computing.

preprint2015arXiv

How Electronic Dynamics with Pauli Exclusion Produces Fermi-Dirac Statistics

It is important that any dynamics method approaches the correct population distribution at long times. In this paper, we derive a one-body reduced density matrix dynamics for electrons in energetic contact with a bath. We obtain a remarkable equation of motion which shows that in order to reach equilibrium properly, rates of electron transitions depend on the density matrix. Even though the bath drives the electrons towards a Boltzmann distribution, hole blocking factors in our equation of motion cause the electronic populations to relax to a Fermi-Dirac distribution. These factors are an old concept, but we show how they can be derived with a combination of time-dependent perturbation theory, and the extended normal ordering of Mukherjee and Kutzelnigg. The resulting non-equilibrium kinetic equations generalize the usual Redfield theory to many-electron systems, while ensuring that the orbital occupations remain between zero and one. In numerical applications of our equations, we show that relaxation rates of molecules are not constant because of the blocking effect. Other applications to model atomic chains are also presented which highlight the importance of treating both dephasing and relaxation. Finally we show how the bath localizes the electron density matrix.

preprint2015arXiv

The Kinetic Energy of Hydrocarbons as a Function of Electron Density and Convolutional Neural Networks

We demonstrate a convolutional neural network trained to reproduce the Kohn-Sham kinetic energy of hydrocarbons from electron density. The output of the network is used as a non-local correction to the conventional local and semi-local kinetic functionals. We show that this approximation qualitatively reproduces Kohn-Sham potential energy surfaces when used with conventional exchange correlation functionals. Numerical noise inherited from the non-linearity of the neural network is identified as the major challenge for the model. Finally we examine the features in the density learned by the neural network to anticipate the prospects of generalizing these models.

preprint2014arXiv

Relaxation Between Bright Optical Wannier Excitons in Perovskite Solar Absorber CH$_3$NH$_3$PbI$_3$

We study the light-absorbing states of the mixed-halide perovskite CH$_{3}$NH$_{3}$PbI$_2$Cl and tri-iodide perovskite CH$_{3}$NH$_{3}$PbI$_3$ with density functional and many-body calculations to explain the desirable photovolatic features of these materials. The short-lived electron-hole bound states produced in this photovoltaic material are of halide to lead electron transfer character, with a Wannier-type exciton. Bethe-Salpeter (GW+BSE) calculations of the absorption cross section reveal strong screening of the electron-hole interaction. The atomic character of the exciton retains ligand-to-metal character within the visible spectrum, with differing degrees of localization outside the unit cell. The average electron-hole separation in the lowest exciton is found to be about 5$A^{\circ}$, slightly larger than the Pb-I bond length. Finally, we determine the role of methylammonium's dipole in the ultrafast relaxation by preparing an atomistic model of the picosecond electronic dynamics in the tri-iodide, PbI$_3$. Our model allows us to identify phonon modes which couple strongly to the electronic excitations, and explain the picosecond timescale intra-band relaxation dynamics seen in recent transient absorption experiments. We largely substantiate the conjectured three band model for the dynamics, but also identify other possible relaxation channels in the tri-iodide.

preprint2013arXiv

Force-Field Functor Theory: Classical Force-Fields which Reproduce Equilibrium Quantum Distributions

Feynman and Hibbs were the first to variationally determine an effective potential whose associated classical canonical ensemble approximates the exact quantum partition function. We examine the existence of a map between the local potential and an effective classical potential which matches the exact quantum equilibrium density and partition function. The usefulness of such a mapping rests in its ability to readily improve Born-Oppenheimer potentials for use with classical sampling. We show that such a map is unique and must exist. To explore the feasibility of using this result to improve classical molecular mechanics, we numerically produce a map from a library of randomly generated one-dimensional potential/effective potential pairs then evaluate its performance on independent test problems. We also apply the map to simulate liquid para-hydrogen, finding that the resulting radial pair distribution functions agree well with path integral Monte Carlo simulations. The surprising accessibility and transferability of the technique suggest a quantitative route to adapting Born-Oppenheimer potentials, with a motivation similar in spirit to the powerful ideas and approximations of density functional theory.

preprint2013arXiv

More accurate and efficient bath spectral densities from super-resolution

Quantum transport and other phenomena are typically modeled by coupling the system of interest to an environment, or bath, held at thermal equilibrium. Realistic bath models are at least as challenging to construct as models for the quantum systems themselves, since they must incorporate many degrees of freedom that interact with the system on a wide range of timescales. Owing to computational limitations, the environment is often modeled with simple functional forms, with a few parameters fit to experiment to yield semi-quantitative results. Growing computational resources have enabled the construction of more realistic bath models from molecular dynamics (MD) simulations. In this paper, we develop a numerical technique to construct these atomistic bath models with better accuracy and decreased cost. We apply a novel signal processing technique, known as super-resolution, combined with a dictionary of physically-motivated bath modes to derive spectral densities from MD simulations. Our approach reduces the required simulation time and provides a more accurate spectral density than can be obtained via standard Fourier transform methods. Moreover, the spectral density is provided as a convenient closed-form expression which yields an analytic time-dependent bath kernel. Exciton dynamics of the Fenna-Matthews-Olsen light-harvesting complex are simulated with a second order time-convolutionless master equation, and spectral densities constructed via super-resolution are shown to reproduce the dynamics using only a quarter of the amount of MD data.

preprint2011arXiv

Exciton coherence lifetimes from electronic structure

We model the coherent energy transfer of an electronic excitation within covalently linked aromatic homodimers from first-principles, to answer whether the usual models of the bath calculated via detailed electronic structure calculations can reproduce the key dynamics. For these systems the timescales of coherent transport are experimentally known from time-dependent polarization anisotropy measurements, and so we can directly assess the whether current techniques might be predictive for this phenomenon. Two choices of electronic basis states are investigated, and their relative merits discussed regarding the predictions of the perturbative model. The coupling of the electronic degrees of freedom to the nuclear degrees of freedom is calculated rather than assumed, and the fluorescence anisotropy decay is directly reproduced. Surprisingly we find that although TDDFT absolute energies are routinely in error by orders of magnitude more than the coupling energy, the coherent transport properties of these dimers can be semi-quantitatively reproduced from first-principles. The directions which must be pursued to yield predictive and reliable prediction of coherent transport are suggested.