BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Paul W. Ayers

Paul W. Ayers contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phcond-mat.str-elquant-phphysics.comp-phastro-ph.EPphysics.ao-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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.

Published work

12 published item(s)

preprint2022arXiv

Bivariational Principle for an Antisymmetrized Product of Nonorthogonal Geminals Appropriate for Strong Electron Correlation

We develop a bivariational principle for an antisymmetric product of nonorthogonal geminals. Special cases reduce to the antisymmetric product of strongly-orthogonal geminals (APSG), the generalized valence bond-perfect pairing (GVB-PP), and the antisymmetrized geminal power (AGP) wavefunctions. The presented method employs wavefunctions of the same type as Richardson-Gaudin (RG) states, but which are not eigenvectors of a model Hamiltonian which would allow for more freedom in the mean-field. The general idea is to work with the same state in a primal picture in terms of pairs, and in a dual picture in terms of pair-holes. This leads to an asymmetric energy expression which may be optimized bivariationally, and is strictly variational when the two representations are consistent. The general approach may be useful in other contexts, such as for computationally feasible variational coupled-cluster methods.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fanpy: A Python Library for Prototyping Multideterminant Methods in Ab Initio Quantum Chemistry

Fanpy is a free and open-source Python library for developing and testing multideterminant wavefunctions and related ab initio methods in electronic structure theory. The main use of Fanpy is to quickly prototype new methods by making it easier to transfer the mathematical conception of a new wavefunction ansätze to a working implementation. Fanpy uses the framework of our recently introduced Flexible Ansatz for N-electron Configuration Interaction (FANCI), where multideterminant wavefunctions are represented by their overlaps with Slater determinants of orthonormal spin-orbitals. In the simplest case, a new wavefunction ansatz can be implemented by simply writing a function for evaluating its overlap with an arbitrary Slater determinant. Fanpy is modular in both implementation and theory: the wavefunction model, the system's Hamiltonian, and the choice of objective function are all independent modules. This modular structure makes it easy for users to mix and match different methods and for developers to quickly try new ideas. Fanpy is written purely in Python with standard dependencies, making it accessible for most operating systems; it adheres to principles of modern software development, including comprehensive documentation, extensive testing, and continuous integration and delivery protocols. This article is considered to be the official release notes for the Fanpy library.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Flexible ansatz for N-body Configuration Interaction

We present a Flexible Ansatz for N-body Configuration Interaction (FANCI) that includes any multideterminant wavefunction. This ansatz is a generalization of the Configuration Interaction (CI) wavefunction, where the coefficients are replaced by a specified function of certain parameters. By making an appropriate choice for this function, we can reproduce popular wavefunction structures like CI, Coupled-Cluster, Tensor Network States, and geminal-product wavefunctions. The universality of this framework suggests a programming structure that allows for the easy construction and optimization of arbitrary wavefunctions. Here, we will discuss the structures of the FANCI framework and its implications for wavefunction properties, particularly accuracy, cost, and size-consistency. We demonstrate the flexibility of this framework by reconstructing popular wavefunction ansätze and modifying them to construct novel wavefunction forms. FANCI provides a powerful framework for exploring, developing, and testing new wavefunction forms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

CRAHCN-O: A Consistent Reduced Atmospheric Hybrid Chemical Network Oxygen Extension for Hydrogen Cyanide and Formaldehyde Chemistry in CO2-, N2-, H2O-, CH4-, and H2-Dominated Atmospheres

Hydrogen cyanide (HCN) and formaldehyde (H2CO) are key precursors to biomolecules such as nucleobases and amino acids in planetary atmospheres; However, many reactions which produce and destroy these species in atmospheres containing CO2 and H2O are still missing from the literature. We use a quantum chemistry approach to find these missing reactions and calculate their rate coefficients using canonical variational transition state theory and Rice-Ramsperger-Kassel-Marcus/master equation theory at the BHandHLYP/aug-cc-pVDZ level of theory. We calculate the rate coefficients for 126 total reactions, and validate our calculations by comparing with experimental data in the 39% of available cases. Our calculated rate coefficients are most frequently within an factor of 2 of experimental values, and generally always within an order of magnitude of these values. We discover 45 previously unknown reactions, and identify 6 from this list that are most likely to dominate H2CO and HCN production and destruction in planetary atmospheres. We highlight $^1$O + CH3 $\rightarrow$ H2CO + H as a new key source, and H2CO + $^1$O $\rightarrow$ HCO + OH as a new key sink, for H2CO in upper planetary atmospheres. In this effort, we develop an oxygen extension to our consistent reduced atmospheric hybrid chemical network (CRAHCN-O), building off our previously developed network for HCN production in N2-, CH4-, and H2-dominated atmospheres (CRAHCN). This extension can be used to simulate both HCN and H2CO production in atmospheres dominated by any of CO2, N2, H2O, CH4, and H2.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Multi-Reference Nature of Plutonium Oxides: PuO$_2^{2+}$, PuO$_2$, PuO$_3$ and PuO$_2$(OH)$_2$

Actinide-containing complexes present formidable challenges for electronic structure methods due to the large number of degenerate or quasi-degenerate electronic states arising from partially occupied 5f and 6d shells. Conventional multi-reference methods can treat active spaces that are often at the upper limit of what is required for a proper treatment of species with complex electronic structures, leaving no room for verifying their suitability. In this work we address the issue of properly defining the active spaces in such calculations, and introduce a protocol to determine optimal active spaces based on the use of the Density Matrix Renormalization Group algorithm and concepts of quantum information theory. We apply the protocol to elucidate the electronic structure and bonding mechanism of volatile plutonium oxides (PuO$_3$ and PuO$_2$(OH)$_2$), species associated with nuclear safety issues for which little is known about the electronic structure and energetics. We show how, within a scalar relativistic framework, orbital-pair correlations can be used to guide the definition of optimal active spaces which provide an accurate description of static/non-dynamic electron correlation, as well as to analyse the chemical bonding beyond a simple orbital model. From this bonding analysis we are able to show that the addition of oxo- or hydroxo-groups to the plutonium dioxide species considerably changes the pi-bonding mechanism with respect to the bare triatomics, resulting in bent structures with considerable multi-reference character.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Longitudinal static optical properties of hydrogen chains: finite field extrapolations of matrix product state calculations

We have implemented the sweep algorithm for the variational optimization of SU(2) x U(1) (spin and particle number) invariant matrix product states (MPS) for general spin and particle number invariant fermionic Hamiltonians. This class includes non-relativistic quantum chemical systems within the Born-Oppenheimer approximation. High-accuracy ab-initio finite field results of the longitudinal static polarizabilities and second hyperpolarizabilities of one-dimensional hydrogen chains are presented. This allows to assess the performance of other quantum chemical methods. For small basis sets, MPS calculations in the saturation regime of the optical response properties can be performed. These results are extrapolated to the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Exact Ionization Potentials from Wavefunction Asymptotics: The Extended Koopmans' Theorem, Revisited

A simple explanation is given for the exactness of the extended Koopmans' theorem for computing the removal energy of any many-electron system to the lowest-energy ground state ion of a given symmetry. In particular, by removing the electron from a "removal orbital" of appropriate symmetry that is concentrated in the asymptotic region, one obtains the exact ionization potential and the exact Dyson orbital for the corresponding state of the ion. It is argued that the EKT is not restricted to many-electron systems but holds for any finite many-body system, provided the interaction vanishes for increasing interparticle distance. A necessary and sufficient condition for the validity of the extended Koopmans' theorem for any state (not just the lowest-energy states of a given symmetry) in terms of the third-order reduced density matrix is stated and derived.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fast density-matrix based partitioning of the energy over the atoms in a molecule consistent with the Hirshfeld-I partitioning of the electron density

For the Hirshfeld-I atom-in-molecule model, associated single-atom energies and interaction energies at the Hartree-Fock level are determined efficiently in one-electron Hilbert space. In contrast to most other approaches, the energy terms are fully consistent with the partitioning of the underlying one-electron density matrix. Starting from the Hirshfeld-I atom-in-molecule model for the electron density, the molecular one-electron density matrix is partitioned with a previously introduced double-atom scheme [Vanfleteren D. et al., J Chem Phys 2010, 132, 164111]. Single-atom density matrices are constructed from the atomic and bond contributions of the double-atom scheme. Since the Hartree-Fock energy can be expressed solely in terms of the one-electron density matrix, the partitioning of the latter over the atoms in the molecule leads naturally to a corresponding partitioning of the Hartree-Fock energy. When the size of the molecule or the molecular basis set does not grow too large, the method shows considerable computational advantages compared to other approaches that require cumbersome numerical integration of the molecular energy integrals weighted by atomic weight functions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High-accuracy ab-initio quantum chemistry by means of an SU(2) x U(1) invariant matrix product state Ansatz: the static second hyperpolarizability

We have implemented the single-site density matrix renormalization group algorithm for the variational optimization of SU(2) \times U(1) (spin and particle number) invariant matrix product states for general spin and particle number symmetric fermionic Hamiltonians. This class also includes non-relativistic quantum chemical systems within the Born-Oppenheimer approximation. High-accuracy ab-initio finite field results of the longitudinal static polarizabilities and second hyperpolarizabilities of one-dimensional hydrogen chains are obtained with the algorithm. A comparison with other methods is made.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hilbert-space partitioning of the molecular one-electron density matrix with orthogonal projectors

A double-atom partitioning of the molecular one-electron density matrix is used to describe atoms and bonds. All calculations are performed in Hilbert space. The concept of atomic weight functions (familiar from Hirshfeld analysis of the electron density) is extended to atomic weight matrices. These are constructed to be orthogonal projection operators on atomic subspaces, which has significant advantages in the interpretation of the bond contributions. In close analogy to the iterative Hirshfeld procedure, self-consistency is built in at the level of atomic charges and occupancies. The method is applied to a test set of about 67 molecules, representing various types of chemical binding. A close correlation is observed between the atomic charges and the Hirshfeld-I atomic charges.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Partitioning of the molecular density matrix over atoms and bonds

A double-index atomic partitioning of the molecular first-order density matrix is proposed. Contributions diagonal in the atomic indices correspond to atomic density matrices, whereas off-diagonal contributions carry information about the bonds. The resulting matrices have good localization properties, in contrast to single-index atomic partitioning schemes of the molecular density matrix. It is shown that the electron density assigned to individual atoms, when derived from the density matrix partitioning, can be made con- sistent with well-known partitions of the electron density over AIM basins, either with sharp or with fuzzy boundaries. The method is applied to a test set of about 50 molecules, representative for various types of chemical binding. A close correlation is observed between the trace of the bond matrices and the SEDI (shared electron density index) bond index.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Subsystem constraints in variational second order density matrix optimization: curing the dissociative behavior

A previous study of diatomic molecules revealed that variational second-order density matrix theory has serious problems in the dissociation limit when the N-representability is imposed at the level of the usual two-index (P, Q, G) or even three-index (T1, T2) conditions [H. van Aggelen et al., Phys. Chem. Chem. Phys. 11, 5558 (2009)]. Heteronuclear molecules tend to dissociate into fractionally charged atoms. In this paper we introduce a general class of N-representability conditions, called subsystem constraints, and show that they cure the dissociation problem at little additional computational cost. As a numerical example the singlet potential energy surface of BeB+ is studied. The extension to polyatomic molecules, where more subsystem choices can be identified, is also discussed.

0 citations0 reviewsNo code linkedOpen access