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

Susi Lehtola

Susi Lehtola contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phphysics.comp-phphysics.atom-phphysics.atm-clusMachine Learning

Trust snapshot

Quick read

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

9 published item(s)

preprint2021arXiv

DQC: a Python program package for Differentiable Quantum Chemistry

Automatic differentiation represents a paradigm shift in scientific programming, where evaluating both functions and their derivatives is required for most applications. By removing the need to explicitly derive expressions for gradients, development times can be be shortened, and calculations simplified. For these reasons, automatic differentiation has fueled the rapid growth of a variety of sophisticated machine learning techniques over the past decade, but is now also increasingly showing its value to support {\it ab initio} simulations of quantum systems, and enhance computational quantum chemistry. Here we present an open-source differentiable quantum chemistry simulation code, DQC, and explore applications facilitated by automatic differentiation: (1) calculating molecular perturbation properties; (2) reoptimizing a basis set for hydrocarbons; (3) checking the stability of self-consistent field wave functions; and (4) predicting molecular properties via alchemical perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Meta-local density functionals: a new rung on Jacob's ladder

The homogeneous electron gas (HEG) is a key ingredient in the construction of most exchange-correlation functionals of density-functional theory. Often, the energy of the HEG is parameterized as a function of its spin density $n$, leading to the local density approximation (LDA) for inhomogeneous systems. However, the connection between the electron density and kinetic energy density of the HEG can be used to generalize the LDA by evaluating it on a weighted geometric average of the local spin density and the spin density of a HEG that has the local kinetic energy density of the inhomogeneous system, with a mixing ratio $x$. This leads to a new family of functionals that we term meta-local density approximations (meta-LDAs), which are still exact for the HEG, which are derived only from properties of the HEG, and which form a new rung of Jacob's ladder of density functionals. The first functional of this ladder, the local $τ$ approximation (LTA) of Ernzerhof and Scuseria that corresponds to $x=1$ is unfortunately not stable enough to be used in self-consistent field calculations, because it leads to divergent potentials as we show in this work. However, a geometric averaging of the LDA and LTA densities with smaller values of $x$ not only leads to numerical stability of the resulting functional, but also yields more accurate exchange energies in atomic calculations than the LDA, the LTA, or the tLDA functional ($x=1/4$) of Eich and Hellgren. We choose $x=0.50$ as it gives the best total energy in self-consistent exchange-only calculations for the argon atom. Atomization energy benchmarks confirm that the choice $x=0.50$ also yields improved energetics in combination with correlation functionals in molecules, almost eliminating the well-known overbinding of the LDA and reducing its error by two thirds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An overview of self-consistent field calculations within finite basis sets

A uniform derivation is presented of the self-consistent field equations in a finite basis set. Both restricted and unrestricted Hartree-Fock (HF) theory as well as various density functional (DF) approximations are considered. The unitary invariance of the HF and DF models is discussed, paving the way for the use of localized molecular orbitals. The self-consistent field equations are derived in a non-orthogonal basis set, and their solution is discussed in the presence of linear dependencies in the basis set. It is argued why iterative diagonalization of the Kohn-Sham-Fock matrix leads to the minimization of the total energy. Alternative methods for the solution of the self-consistent field equations via direct minimization as well as stability analysis are also briefly discussed. Explicit expressions are given for the contributions to the Kohn-Sham-Fock matrix up to meta-GGA functionals. Range-separated hybrids and non-local correlation functionals are also briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets

The superposition of atomic potentials (SAP) approach has recently been shown to be a simple and efficient way to initialize electronic structure calculations [S. Lehtola, J. Chem. Theory Comput. 15, 1593 (2019)]. Here, we study the differences between effective potentials from fully numerical density functional and optimized effective potential calculations for fixed configurations. We find that the differences are small, overall, and choose exchange-only potentials at the local density approximation level of theory computed on top of Hartree-Fock densities as a good compromise. The differences between potentials arising from different atomic configurations are also found to be small at this level of theory. Furthermore, we discuss the efficient Gaussian-basis implementation of SAP via error function fits to fully numerical atomic radial potentials. The guess obtained from the fitted potentials can be easily implemented in any Gaussian-basis quantum chemistry code in terms of two-electron integrals. Fits covering the whole periodic table from H to Og are reported for non-relativistic as well as fully relativistic four-component calculations that have been carried out with fully numerical approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the accurate reproduction of strongly repulsive interatomic potentials

Knowledge of the repulsive behavior of potential energy curves $V(R)$ at $R\to0$ is necessary for understanding and modeling irradiation processes of practical interest. $V(R)$ is in principle straightforward to obtain from electronic structure calculations; however, commonly-used numerical approaches for electronic structure calculations break down in the strongly repulsive region due to the closeness of the nuclei. In the present work, we show by comparison to fully numerical reference values that a recently developed procedure [S. Lehtola, J. Chem. Phys. 151, 241102 (2019)] can be employed to enable accurate linear combination of atomic orbitals calculations of $V(R)$ even at small $R$ by a study of the seven nuclear reactions He2 <=> Be, HeNe <=> Mg, Ne2 <=> Ca, HeAr <=> Ca, MgAr <=> Zn, Ar2 <=> Kr, and NeCa <=> Zn.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Polarized Gaussian basis sets from one-electron ions

We demonstrate that basis sets suitable for electronic structure calculations can be obtained from simple accuracy considerations for the hydrogenic one-electron ions $Y^{(Y-1)+}$ for $Y\in[1,Z]$, necessitating no self-consistent field calculations at all. It is shown that even-tempered basis sets with parameters from the commonly-used universal Gaussian basis set (UGBS) [E. V. R. de Castro and F. E. Jorge, J. Chem. Phys. 108, 5225 (1998)] reproduce non-relativistic spin-restricted spherical Hartree-Fock total energies from fully numerical calculations to better accuracy than UGBS, which is shown to exhibit huge errors for some elements, e.g. 0.19 $E_{h}$ for Th$^+$ and 0.13 $E_{h}$ for Lu as it has been parametrized for a single atomic configuration. Having shown the feasibility of the one-electron approach, partially energy-optimized basis sets are formed for all atoms in the periodic table, $1\leq Z\leq118$, by optimizing the even-tempered parameters for $Z^{(Z-1)+}$. As the hydrogenic Gaussian basis sets suggested in the present work are built strictly from first principles, also polarization shells can be obtained in the same fashion in contrast to previous approaches. The accuracy of the polarized basis sets is demonstrated by calculations on a small set of molecules by comparison to fully numerical reference values, which show that chemical accuracy can be reached even for challenging cases like SF$_6$. The present approach is straightforward to extend to relativistic calculations, and could facilitate studies beyond the established periodic table.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

PyFLOSIC: Python-based Fermi-Löwdin orbital self-interaction correction

We present PyFLOSIC, an open-source, general-purpose Python implementation of the Fermi-Löwdin orbital self-interaction correction (FLO-SIC), which is based on the Python simulation of chemistry frame-work (PySCF) electronic structure and quantum chemistry code. Thanks to PySCF, PyFLOSIC can be used with any kind of Gaussian-type basis set, various kinds of radial and angular quadrature grids, and all exchange-correlation functionals within the local density approximation (LDA), generalized-gradient approximation (GGA), and meta-GGA provided in the Libxc and XCFun libraries. A central aspect of FLO-SIC are Fermi-orbital descriptors, which are used to estimate the self-interaction correction. Importantly, they can be initialized automatically within PyFLOSIC and optimized with an interface to the atomic simulation environment, a Python library which provides a variety of powerful gradient-based algorithms for geometry optimization. Although PyFLOSIC has already facilitated applications of FLO-SIC to chemical studies, it offers an excellent starting point for further developments in FLO-SIC approaches, thanks to its use of a high-level programming language and pronounced modularity.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Fully numerical calculations on atoms with fractional occupations. Range-separated exchange functionals

A recently developed finite element approach for fully numerical atomic structure calculations [S. Lehtola, Int. J. Quantum Chem. 119, e25945 (2019)] is extended to the description of atoms with spherically symmetric densities via fractionally occupied orbitals. Specialized versions of Hartree-Fock as well as local density and generalized gradient approximation density functionals are developed, allowing extremely rapid calculations at the basis set limit on the ground and low-lying excited states even for heavy atoms. The implementation of range-separation based on the Yukawa or complementary error function (erfc) kernels is also described, allowing complete basis set benchmarks of modern range-separated hybrid functionals with either integer or fractional occupation numbers. Finally, computation of atomic effective potentials at the local density or generalized gradient approximation levels for the superposition of atomic potentials (SAP) approach [S. Lehtola, J. Chem. Theory Comput. 15, 1593 (2019)] that has been shown to be a simple and efficient way to initialize electronic structure calculations is described. The present numerical approach is shown to afford beyond microhartree accuracy with a small number of numerical basis functions, and to reproduce literature results for the ground states of atoms and their cations for $1 \leq Z \leq 86 $. Our results indicate that the literature values deviate by up to 10 μEh from the complete basis set limit. The numerical scheme for the erfc kernel is shown to work by comparison to results from large Gaussian basis set calculations from the literature. Spin-restricted ground states are reported for Hartree-Fock and Hartree-Fock-Slater calculations with fractional occupations for $1 \leq Z \leq 118$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Fully numerical electronic structure calculations on diatomic molecules in weak to strong magnetic fields

We present fully numerical electronic structure calculations on diatomic molecules exposed to an external magnetic field at the unrestricted Hartree-Fock limit, using a modified version of a recently developed finite element program, HelFEM. We have performed benchmark calculations on a few low-lying states of H2, HeH+, LiH, BeH+, BH, and CH+ as a function of the strength of an external magnetic field parallel to the molecular axis. The employed magnetic fields are in the range of $B=[0,10]~B_0$ atomic units, where $B_0 \approx 2.35 \times 10^5$ T. We have compared the results of the fully numerical calculations to ones obtained with the LONDON code using a large uncontracted gauge-including Cartesian Gaussian (GICG) basis set with exponents adopted from the Dunning aug-cc-pVTZ basis set. By comparison to the fully numerical results, we find that the basis set truncation error in the gauge-including Gaussian basis set is of the order of 1 kcal/mol at zero field, that the truncation error grows rapidly when the strength of the magnetic field increases, and that the largest basis set truncation error at $B=10~B_0$ exceeds 1000 kcal/mol. Studies in larger Gaussian basis sets suggest that reliable results can be obtained in GICG basis sets at fields stronger than $B=B_0$, provided that a sufficient coverage of higher-angular-momentum functions is included in the basis set.

0 citations0 reviewsNo code linkedOpen access