Publications

Order by: publication date number of citations

• A benchmark study of dioxygen complexes based on coupled cluster and density functional theory

  Marcel Swart, Marc Reimann

  SciPost Chem. 3, 001 (2024) · published 13 August 2024 | Toggle abstract · pdf

  A set of five compounds containing peroxo, superoxo or bis-μ-oxo moieties has been studied in the gas phase using CCSD(T)/aug-cc-pVQZ, also in combination with Goodson’s continued fraction approach. The corresponding analytical frequencies corroborate assignments of bands from experiments, and thus provide a consistent set of reference data that can be used for benchmarking a range of density functional approximations. A total of 100 density functionals have been checked for the general bond lengths, the specific peroxo/superoxo bond lengths, angles, and vibrational frequencies. There is not one density functional that performs equally well for all of these properties, not even within one class of density functionals.

• Do machines dream of atoms? Crippen's logP as a quantitative molecular benchmark for explainable AI heatmaps

  Maria H. Rasmussen, Diana S. Christensen, Jan H. Jensen

  SciPost Chem. 2, 002 (2023) · published 27 June 2023 | Toggle abstract · pdf

  While there is a great deal of interest in methods aimed at explaining machine learning predictions of chemical properties, it is difficult to quantitatively benchmark such methods, especially for regression tasks. We show that the Crippen logP model (J. Chem. Inf. Comput. Sci. 1999, 39, 868) provides an excellent benchmark for atomic attribution/heatmap approaches, especially if the ground truth heatmaps can be adjusted to reflect the molecular representation. The "atom attribution from finger prints"-method developed by Riniker and Landrum (J. Chem. Inf. Comput. Sci. 2013, 5, 43) gives atomic attribution heatmaps that are in reasonable agreement with the atomic contribution heatmaps of the Crippen logP model for most molecules, with average heatmap overlaps of up to 0.54. The agreement is increased significantly (to 0.75) when the atomic contributions are adjusted to match the fact that the molecular representation is fragment-based rather than atom-based (the finger print-adapted (FPA) ground truth vector). Most heatmaps and the corresponding FPA overlaps are relatively insensitive to the training set size and the results are close to converged for a training set size of 1000 molecules, although for molecules with low overlap some heatmaps change significantly. Using the "remove an atom" approach for graph convolutional neural networks (GCNNs) suggested by Matveieva and Polishchuk (J. Cheminform. 2021, 13, 41) we find an average heatmap overlap of 0.47 for the atomic contribution heatmaps of the Crippen logP model. Like the simpler attribution benchmarks for classification tasks that have come before it, this work sets the bar for regression tasks.

  4 citations

• Ether cleavage and chemical removal of SU-8

  José M. Ripalda, Raquel Álvaro, María Luisa Dotor

  SciPost Chem. 2, 001 (2023) · published 13 April 2023 | Toggle abstract · pdf

  The high chemical stability of SU-8 makes it irreplaceable for a wide range of applications, most notably as a lithography photoresist for micro and nanotechnology. This advantage becomes a problem when there is a need to remove SU-8 from the fabricated devices. Researchers have been struggling for two decades with this problem, and although a number of partial solutions have been found, this difficulty has limited the applications of SU-8. Here we demonstrate a fast, reproducible, and comparatively gentle method to chemically remove SU-8 photoresist. An ether cleavage mechanism for the observed reaction is proposed, and the hypothesis is tested with ab initio quantum chemical calculations. We also describe a complementary removal method based on atomic hydrogen inductively coupled plasma.

  1 citation

• CrystalNets.jl: Identification of Crystal Topologies

  Lionel Zoubritzky, François-Xavier Coudert

  SciPost Chem. 1, 005 (2022) · published 17 June 2022 | Toggle abstract · pdf

  We present here an open-source Julia library for the topological identification of crystalline materials, with algorithmic and computational improvements over the previously available software in the field, resulting in a speed increase of one order of magnitude. This new algorithm and implementation can therefore be used at large scale in high-throughput screening methodologies. We have validated and benchmarked CrystalNets.jl against a diverse set of crystal databases, covering in particular metal–organic frameworks, aluminophosphates, zeolites, and other inorganic compounds.

  11 citations

• Relativistic reduced density matrix functional theory

  Mauricio Rodríguez-Mayorga, Klaas J.H. Giesbertz, Lucas Visscher

  SciPost Chem. 1, 004 (2022) · published 10 June 2022 | Toggle abstract · pdf

  As a new approach to efficiently describe correlation effects in the relativistic quantum world we propose to consider reduced density matrix functional theory, where the key quantity is the first-order reduced density matrix (1-RDM). In this work, we first introduce the theoretical foundations to extend the applicability of this theory to the relativistic domain. Then, using the so-called no-pair (np) approximation, we arrive at an approximate treatment of the relativistic effects by focusing on electronic wavefunctions and neglecting explicit contributions from positrons. Within the np approximation the theory becomes similar to the nonrelativistic case, with as unknown only the functional that describes the electron-electron interactions in terms of the 1-RDM. This requires the construction of functional approximations, and we therefore also present the relativistic versions of some common RDMFT approximations that are used in the nonrelativistic context and discuss their properties.

  12 citations

• Fast and automated identification of reactions with low barriers: the decomposition of 3-hydroperoxypropanal

  Mads Koerstz, Maria H. Rasmussen, Jan H. Jensen

  SciPost Chem. 1, 003 (2021) · published 1 October 2021 | Toggle abstract · pdf

  We show how fast semiempirical QM methods can be used to significantly decrease the CPU requirements for automated reaction mechanism discovery, using two different method for generating reaction products: graph-based systematic enumeration of all possible products and the meta-dynamics approach by Grimme (J. Chem. Theory. Comput. 2019, 15, 2847). We test the two approaches on the low-barrier reactions of 3-hydroperoxypropanal, which have been studied by a large variety of reaction discovery approaches and therefore provides a good benchmark. By using PM3 and GFN2-xTB for reaction energy and barrier screening the systematic approach identifies 64 reactions (out of 27,577 possible reactions) for DFT refinement, which in turn identifies the three reactions with lowest barriers plus a previously undiscovered reaction. With optimised hyperparameters meta-dynamics followed by PM3/GFN2-xTB-based screening identifies 15 reactions for DFT refinement, which in turn identifies the three reactions with lowest barrier. The number of DFT refinements can be further reduced to as little as six for both approaches by first verifying the transition states with GFN1-xTB. The main conclusion is that the semiempirical methods are accurate and fast enough to automatically identify promising candidates for DFT refinement for the low barrier reactions of 3-hydroperoxypropanal in about 15-30 minutes using relatively modest computational resources.

  8 citations

• Relativistic density-functional theory based on effective quantum electrodynamics

  Julien Toulouse

  SciPost Chem. 1, 002 (2021) · published 18 May 2021 | Toggle abstract · pdf

  A relativistic density-functional theory based on a Fock-space effective quantum-electrodynamics (QED) Hamiltonian using the Coulomb or Coulomb-Breit two-particle interaction is developed. This effective QED theory properly includes the effects of vacuum polarization through the creation of electron-positron pairs but does not include explicitly the photon degrees of freedom. It is thus a more tractable alternative to full QED for atomic and molecular calculations. Using the constrained-search formalism, a Kohn-Sham scheme is formulated in a quite similar way to non-relativistic density-functional theory, and some exact properties of the involved density functionals are studied, namely charge-conjugation symmetry and uniform coordinate scaling. The usual no-pair Kohn-Sham scheme is obtained as a well-defined approximation to this relativistic density-functional theory.

  5 citations

• The seniority quantum number in Tensor Network States

  Klaas Gunst, Dimitri Van Neck, Peter Andreas Limacher, Stijn De Baerdemacker

  SciPost Chem. 1, 001 (2021) · published 15 January 2021 | Toggle abstract · pdf

  We employ tensor network methods for the study of the seniority quantum number -- defined as the number of unpaired electrons in a many-body wave function -- in molecular systems. Seniority-zero methods recently emerged as promising candidates to treat strong static correlations in molecular systems, but are prone to deficiencies related to dynamical correlation and dispersion. We systematically resolve these deficiencies by increasing the allowed seniority number using tensor network methods. In particular, we investigate the number of unpaired electrons needed to correctly describe the binding of the neon and nitrogen dimer and the $D_{6h}$ symmetry of benzene.

  6 citations