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

Martin Head-Gordon

Martin Head-Gordon contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phquant-phcond-mat.otherphysics.comp-phcond-mat.mtrl-scicond-mat.str-elnucl-thphysics.atm-clusphysics.atom-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
9topics
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

17 published item(s)

preprint2023arXiv

Generalization of one-center non orthogonal configuration interaction singles to open shell singlet reference states: Theory and application to valence-core pump-probe states in acetylacetone

We formulate a one-center non-orthogonal configuration interaction singles (1C-NOCIS) theory for the computation of core excited states of an initial singlet state with two unpaired electrons. This model, which we refer to as 1C-NOCIS two-electron open-shell (2eOS), is appropriate for computing the K-edge near-edge X-ray absorption spectra (NEXAS) of the valence excited states of closed-shell molecules relevant to pump-probe time-resolved (TR) NEXAS experiments. With inclusion of core hole relaxation effects and explicit spin adaptation, 1C-NOCIS 2eOS requires mild shifts to match experiment, is free of artifacts due to spin contamination, and can capture the high-energy region of the spectrum beyond the transitions into the singly occupied molecular orbitals (SOMO). Calculations on water and thymine illustrate the different key features of excited-state NEXAS, namely the core-to-SOMO transition as well as shifts and spin-splittings in the transitions analogous to those of the ground state. Finally, simulations of the TR-NEXAS of acetylacetone after excitation onto its pi to pi-star singlet excited state at the carbon K-edge - an experiment carried out recently - showcases the ability of 1C-NOCIS 2eOS to efficiently simulate NEXAS based on non-adiabatic molecular dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach

We investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without any use of response theory or equation-of-motion formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for Delta-CCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed Delta-CCSD schemes seem to fare better than the widely used equation-of-motion CCSD (EOM-CCSD) with core-valence separation protocol, with statistical errors being reduced by more than a factor of two when compared to FC-CVS-EOM-CCSD.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Computing X-ray absorption spectra from linear-response particles atop optimized holes

State specific orbital optimized density functional theory (OO-DFT) methods like restricted open-shell Kohn-Sham (ROKS) can attain semiquantitative accuracy for predicting X-ray absorption spectra of closed-shell molecules. OO-DFT methods however require that each state be individually optimized. In this work, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, that is capable of giving reasonable accuracy without requiring state-specific optimization. This is achieved by fully optimizing the core-hole through the core-ionized state, followed by use of electron-addition configuration interaction singles (EA-CIS) to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method, and can attain $\sim 0.6$ eV RMS error for the K-edges of C-F through the use of local functionals like PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment, and can thus help reduce the computational cost for obtaining OO-DFT quality spectra. ROKS(STEX) therefore appears to be a useful technique for efficient prediction of X-ray absorption spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Jahn-Teller Distortion and Dissociation of CCl$_4^+$ by Transient X-ray Spectroscopy Simultaneously at the Carbon K- and Chlorine L-Edge

X-ray Transient Absorption Spectroscopy (XTAS) and theoretical calculations are used to study CCl$_4^+$ prepared by 800 nm strong-field ionization. XTAS simultaneously probes atoms at the carbon K-edge (280-300 eV) and chlorine L-edge (195-220 eV). Comparison of experiment to X-ray spectra computed by orbital-optimized density functional theory (OO-DFT) indicates that after ionization, CCl$_4^+$ undergoes symmetry breaking driven by Jahn-Teller distortion away from the initial tetrahedral structure (T$_d$) in 6$\pm$2 fs. The resultant symmetry-broken covalently bonded form subsequently separates to a noncovalently bound complex between CCl$_3^+$ and Cl over 90$\pm$10 fs, which is again predicted by theory. Finally, after more than 800 fs, L-edge signals for atomic Cl are observed, indicating dissociation to free CCl$_3^+$ and Cl. The results for Jahn-Teller distortion to the symmetry-broken form of CCl$_4^+$ and formation of the Cl -- CCl$_3^+$ complex characterize previously unobserved new species along the route to dissociation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Relativistic Orbital Optimized Density Functional Theory for Accurate Core-Level Spectroscopy

Core-level spectra of 1s electrons of elements heavier than Ne show significant relativistic effects. We combine advances in orbital optimized DFT (OO-DFT) with the spin-free exact two-component (X2C) model for scalar relativistic effects, to study K-edge spectra of third period elements. OO-DFT/X2C is found to be quite accurate at predicting energies, yielding $\sim 0.5$ eV RMS error vs experiment with the modern SCAN (and related) functionals. This marks a significant improvement over the $>50$ eV deviations that are typical for the popular time-dependent DFT (TDDFT) approach. Consequently, experimental spectra are quite well reproduced by OO-DFT/X2C, sans empirical shifts for alignment. OO-DFT/X2C combines high accuracy with ground state DFT cost and is thus a promising route for computing core-level spectra of third period elements. We also explored K and L edges of 3d transition metals to identify limitations of the OO-DFT/X2C approach in modeling the spectra of heavier atoms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Revisiting the performance of time-dependent density functional theory for electronic excitations: Assessment of 43 popular and recently developed functionals from rungs one to four

In this paper, the performance of more than 40 popular or recently developed density functionals is assessed for the calculation of 463 vertical excitation energies against the large and accurate QuestDB benchmark set. For this purpose, the Tamm-Dancoff approximation offers a good balance between performance and accuracy. The functionals $ω$B97X-D and BMK are found to offer the best performance overall with a Root-Mean Square Error (RMSE) of 0.28 eV, better than the computationally more demanding CIS(D) wavefunction method with a RMSE of 0.36 eV. The results also suggest that Jacob's ladder still holds for TDDFT excitation energies, though hybrid meta-GGAs are not generally better than hybrid GGAs. Effects of basis set convergence, gauge invariance correction to meta-GGAs, and nonlocal correlation (VV10) are also studied, and practical basis set recommendations are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Say NO to Optimization: A Non-Orthogonal Quantum Eigensolver

A balanced description of both static and dynamic correlations in electronic systems with nearly degenerate low-lying states presents a challenge for multi-configurational methods on classical computers. We present here a quantum algorithm utilizing the action of correlating cluster operators to provide high-quality wavefunction ansätze employing a non-orthogonal multireference basis that captures a significant portion of the exact wavefunction in a highly compact manner, and that allows computation of the resulting energies and wavefunctions at polynomial cost with a quantum computer. This enables a significant improvement over the corresponding classical non-orthogonal solver, which incurs an exponential cost when evaluating off-diagonal matrix elements between the ansatz states, and is therefore intractable. We implement the non-orthogonal quantum eigensolver (NOQE) here with an efficient ansatz parameterization inspired by classical quantum chemistry methods that succeed in capturing significant amounts of electronic correlation accurately. By taking advantage of classical methods for chemistry, NOQE provides a flexible, compact, and rigorous description of both static and dynamic electronic correlation, making it an attractive method for the calculation of electronic states of a wide range of molecular systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ultrafast X-ray Spectroscopy of Intersystem Crossing in Hexafluoroacetylacetone: Chromophore Photophysics and Spectral Changes in the Face of Electron Withdrawing Groups

Intersystem crossings between singlet and triplet states represent a crucial relaxation pathway in photochemical processes. Herein, we probe the intersystem crossing in hexafluoro-acetylacetone with ultrafast X-ray transient absorption spectroscopy at the carbon K-edge. We observe the excited state dynamics following excitation with 266 nm UV light to the $^1ππ^{*}$ (S$_2$) state with element and site-specificity using a broadband soft X-ray pulse produced by high harmonic generation. These results are compared to X-ray spectra computed from orbital optimized density functional theory methods. It is found that the electron withdrawing fluorine atoms decongest the X-ray absorption spectrum by enhancing separation between features originating from different carbon atoms. This facilitates the elucidation of structural and electronic dynamics at the chromophore. The evolution of the core-to-valence resonances at the carbon K-edge reveals an ultrafast population transfer between the $^1nπ^{*}$ (S$_1$) and $^3ππ^{*}$ (T$_1$) states on a $1.6\pm0.4$ ps timescale, which is similar to the 1.5 ps timescale earlier observed for acetylacetone [J. Am. Chem. Soc. 139, 16576 (2017)]. It therefore appears that terminal fluorination has little influence on the intersystem crossing rate of the acetylacetone chromophore. In addition, the significant role of hydrogen-bond opened and twisted rotational isomers is elucidated in the excited state dynamics by comparison of the experimental transient X-ray spectra with theory.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Assessment of the Performance of Density Functionals for Predicting Potential Energy Curves in Hydrogen Storage Applications

The availability of accurate computational tools for modeling and simulation is vital to accelerate the discovery of materials capable of storing hydrogen (H2) under given parameters of pressure swing and temperature. Previously, we compiled the H2Bind275 dataset consisting of equilibrium geometries and assessed the performance of 55 density functionals over this dataset (Veccham, S. P.; Head-Gordon, M. J. Chem. Theory Comput., 2020, 16, 4963--4982). As it is crucial for computational tools to accurately model the entire potential energy curve (PEC), in addition to the equilibrium geometry, we have extended this dataset with 389 new data points to include two compressed and three elongated geometries along 78 PECs for H2 binding, forming the H2Bind78x7 dataset. Assessing the performance of 55 density functionals on this significantly larger and more comprehensive H2Bind78x7 dataset, we have identified the best performing density functionals for H2 binding applications: PBE0-DH, $ω$B97X-V, $ω$B97M-V, and DSD-PBEPBE-D3(BJ). Addition of Hartree Fock exchange improves the performance of density functionals, albeit not uniformly throughout the PEC. We recommend the usage of wB97X-V and wB97M-V density functionals as they give good performance for both geometries and energies. In addition, we have also identified B97M-V and B97M-rV as the best semi-local density functionals for predicting H2 binding energy at its equilibrium geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Changes in polarization dictate necessary approximations for modeling electronic de-excitation intensity: an application to X-ray emission

We systematically investigate the underlying relations among different levels of approximation for simulating electronic de-excitations, with a focus on modeling X-ray emission spectroscopy (XES). Using Fermi's golden rule and explicit modeling of the initial, core-excited state and the final, valence-hole state, we show that XES can be accurately modeled by using orbital optimization for the various final states within a Slater-determinant framework. However, in this paper, we introduce a much cheaper approach reliant only on a single self-consistent field for all the final states, and show that it is typically sufficient. Further approximations reveal that these fundamentally many-body transitions can be reasonably approximated by projections of ground state orbitals, but that the ground state alone is insufficient. Furthermore, except in cases where the core-ionization induces negligible changes in polarization, linear-response approaches within the adiabatic approximation will have difficulty in accurately modeling de-excitation to the core level. Therefore, change in the net dipole moment of the valence electrons can serve as a metric for the validity of the linear-response approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Third order Møller-Plesset theory made more useful? The role of density functional theory orbitals

The practical utility of Møller-Plesset (MP) perturbation theory is severely constrained by the use of Hartree-Fock (HF) orbitals. It has recently been shown that use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error (J. Phys. Chem. Lett. 10, 4170, 2019). In this work we examine whether density functional theory (DFT) optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, non-covalent interactions, and dipole moments relative to standard HF based MP theory. Indeed MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled cluster singles and doubles (CCSD) for several datasets. We also found that the results are not particularly functional sensitive in most cases, (although range-separated hybrid functionals with low delocalization error perform the best). MP3 based on DFT orbitals thus appears to be an efficient, non-iterative $O(N^6)$ scaling wave function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Too big, too small or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory (DFT) prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher order moments origin-dependent. Therefore, instead of using the $3 \times 3$ quadrupole moment matrix, we utilize the translationally invariant $3 \times 3$ matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by $\mathcal{K}$). The principal components of ${\mathcal{K}}$ are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the prinicpal components of ${\mathcal{K}}$ for 100 small molecules at the coupled cluster singles and doubles with perturbative triples (CCSD(T)) at the complete basis set (CBS) limit is developed, resulting in 213 independent ${\mathcal{K}}$ components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 yield reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H and Be atoms in particular are challenging for nearly all methods, indicating that future functional development could benefit from inclusion of their density information in training or testing protocols.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a dataset (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this dataset using coupled cluster theory. Secondly, we have assessed the performance of 55 density functional approximations for predicting H$_2$ interaction energies and have identified two hybrid density functionals ($ω$B97X-V and $ω$B97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semi-local density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals like revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1kJ/mol.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Highly Accurate Prediction of Core Spectra of Molecules at Density Functional Theory Cost: Attaining sub eV Error from a Restricted Open-Shell Kohn-Sham Approach

We present the use of the recently developed Square Gradient Minimization (SGM) algorithm for excited state orbital optimization, to obtain spin-pure Restricted Open-Shell Kohn-Sham (ROKS) energies for core excited states of molecules. The SGM algorithm is robust against variational collapse, and offers a reliable route to converging orbitals for target excited states at only 2-3 times the cost of ground state orbital optimization (per iteration). ROKS/SGM with the modern SCAN/$ω$B97X-V functionals is found to predict the K edge of C,N,O and F to a root mean squared error of $\sim$0.3 eV. ROKS/SGM is equally effective at predicting L edge spectra of third period elements, provided a perturbative spin-orbit correction is employed. This high accuracy can be contrasted with traditional TDDFT, which typically has greater than 10 eV error and requires translation of computed spectra to align with experiment. ROKS is computationally affordable (having the same scaling as ground state DFT, and a slightly larger prefactor) and can be applied to geometry optimizations/ab-initio molecular dynamics of core excited states, as well as condensed phase simulations. ROKS can also model doubly excited/ionized states with one broken electron pair, which are beyond the ability of linear response based methods.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Excited state orbital optimization via minimizing the square of the gradient: General approach and application to singly and doubly excited states via density functional theory

We present a general approach to converge excited state solutions to any quantum chemistry orbital optimization process, without the risk of variational collapse. The resulting Square Gradient Minimization (SGM) approach only requires analytic energy/Lagrangian orbital gradients and merely costs 3 times as much as ground state orbital optimization (per iteration), when implemented via a finite difference approach. SGM is applied to both single determinant $Δ$SCF and spin-purified Restricted Open-Shell Kohn-Sham (ROKS) approaches to study the accuracy of orbital optimized DFT excited states. It is found that SGM can converge challenging states where the Maximum Overlap Method (MOM) or analogues either collapse to the ground state or fail to converge. We also report that $Δ$SCF/ROKS predict highly accurate excitation energies for doubly excited states (which are inaccessible via TDDFT). Singly excited states obtained via ROKS are also found to be quite accurate, especially for Rydberg states that frustrate (semi)local TDDFT. Our results suggest that orbital optimized excited state DFT methods can be used to push past the limitations of TDDFT to doubly excited, charge-transfer or Rydberg states, making them a useful tool for the practical quantum chemist's toolbox for studying excited states in large systems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller-Plesset Perturbation Theory

We present a systematically improvable tensor hypercontraction (THC) factorization based on interpolative separable density fitting (ISDF). We illustrate algorithmic details to achieve this within the framework of Becke's atom-centered quadrature grid. A single ISDF parameter $c_\text{ISDF}$ controls the tradeoff between accuracy and cost. In particular, $c_\text{ISDF}$ sets the number of interpolation points used in THC, $N_\text{IP} = c_\text{ISDF}\times N_\text{X}$ with $N_\text{X}$ being the number of auxiliary basis functions. In conjunction with the resolution-of-the-identity (RI) technique, we develop and investigate the THC-RI algorithms for cubic-scaling exact exchange for Hartree-Fock and range-separated hybrids (e.g., $ω$B97X-V) and quartic-scaling second- and third-order Møller-Plesset theory (MP2 and MP3). These algorithms were evaluated over the W4-11 thermochemistry (atomization energy) set and A24 non-covalent interaction benchmark set with standard Dunning basis sets (cc-pVDZ, cc-pVTZ, aug-cc-pVDZ, and aug-cc-pVTZ). We demonstrate the convergence of THC-RI algorithms to numerically exact RI results using ISDF points. Based on these, we make recommendations on $c_\text{ISDF}$ for each basis set and method. We also demonstrate the utility of THC-RI exact exchange and MP2 for larger systems such as water clusters and $\text{C}_{20}$. We stress that more challenges await in obtaining accurate and numerically stable THC factorization for wavefunction amplitudes as well as the space spanned by virtual orbitals in large basis sets and implementing sparsity-aware THC-RI algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Poisson-Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients

We present the theory and implementation of a Poisson-Boltzmann implicit solvation model for electrolyte solutions. This model can be combined with arbitrary electronic structure methods that provide an accurate charge density of the solute. A hierarchy of approximations for this model includes a linear approximation for weak electrostatic potentials, finite size of the mobile electrolyte ions and a Stern-layer correction. Recasting the Poisson-Boltzmann equations into Euler-Lagrange equations then significantly simplifies the derivation of the free energy of solvation for these approximate models. The parameters of the model are then either fit directly to experimental observables, e.g. the finite ion size, or optimized for agreement with experimental results. Experimental data for this optimization is available in the form of Sechenov coefficients that describe the linear dependence of the salting-out effect of solutes with respect to the electrolyte concentration. In the final part we rationalize the qualitative disagreement of the finite ion size modification to the Poisson-Boltzmann model with experimental observations by taking into account the electrolyte concentration dependence of the Stern layer. A route towards a revised model that captures the experimental observations while including the finite ion size effects is then outlined. This implementation paves the way for the study of electrochemical and electrocatalytic processes of molecules and cluster models with accurate electronic structure methods.

0 citations0 reviewsNo code linkedOpen access