Solution of the Kohn–Sham equations for many-electron atoms confined by penetrable walls

CuGBasis: High-performance CUDA/Python library for efficient computation of quantum chemistry density-based descriptors for larger systems

CuGBasis is a free and open-source CUDA®/Python library for efficient computation of scalar, vector, and matrix quantities crucial for the post-processing of electronic structure calculations. CuGBasis integrates high-performance Graphical Processing Unit (GPU) computing with the ease and flexibility of Python programming, making it compatible with a vast ecosystem of libraries. We showcase its utility as a Python library and demonstrate its seamless interoperability with existing Python software to gain chemical insight from quantum chemistry calculations. Leveraging GPU-accelerated code, cuGBasis exhibits remarkable performance, making it highly applicable to larger systems or large databases. Our benchmarks reveal a 100-fold performance gain compared to alternative software packages, including serial/multi-threaded Central Processing Unit and GPU implementations. This paper outlines various features and computational strategies that lead to cuGBasis's enhanced performance, guiding developers of GPU-accelerated code.

Computational tools to study non-covalent interactions and confinement effects in chemical systems

Confinement is a very common phenomenon in chemistry, for example, when molecules are located inside cavities. In these conditions, the electronic structure of atoms and molecules is modified. These changes could be mapped through the interaction with other molecules since non-covalent interactions between molecules are also influenced by confinement. In this work we address both topics, non-covalent interactions, and confined systems, using quantum chemistry tools with new software, emphasizing the importance of analyzing both fields simultaneously.

Finite element method as an alternative to study the electronic structure of confined atoms

The finite element method (FEM) based on a nonregular mesh is used to solve Hartree-Fock and Kohn-Sham equations for three atoms (hydrogen, helium, and beryllium) confined by finite and infinite potentials, defined in terms of piecewise functions or functions with a well-defined first derivative. This approach's reliability is shown when contrasted with Roothaan's approach, which depends on a basis set. Therefore, its exponents must be optimized for each confinement imposed over each atom, which is a monumental task. The comparison between our numerical approach and Roothaan's approach is made by using total and orbitals energies from the Hartree-Fock method, where there are several comparison sources. Regarding the Kohn-Sham method, there are few published data and consequently the results reported here can be used as a benchmark for future comparisons. The way to solve Hartree-Fock or Kohn-Sham equations by the FEM is entirely appropriate to study confined atoms with any form of confinement potential. This article represents a step toward developing a fully numerical quantum chemistry code free of basis sets to obtain the electronic structure of many-electron atoms confined by arbitrary confinement.

Electronegativity under Confinement

The electronegativity concept was first formulated by Pauling in the first half of the 20th century to explain quantitatively the properties of chemical bonds between different types of atoms. Today, it is widely known that, in high-pressure regimes, the reactivity properties of atoms can change, and, thus, the bond patterns in molecules and solids are affected. In this work, we studied the effects of high pressure modeled by a confining potential on different definitions of electronegativity and, additionally, tested the accuracy of first-order perturbation theory in the context of density functional theory for confined atoms of the second row at the Hartree–Fock level. As expected, the electronegativity of atoms at high confinement is very different than that of their free counterparts since it depends on the electronic configuration of the atom, and, thus, its periodicity is modified at higher pressures.

GPU-Accelerated Large-Scale Excited-State Simulation Based on Divide-and-Conquer Time-Dependent Density-Functional Tight-Binding

The present study implemented the divide-and-conquer time-dependent density-functional tight-binding (DC-TDDFTB) code on a graphical processing unit (GPU). The DC method, which is a linear-scaling scheme, divides a total system into several fragments. By separately solving local equations in individual fragments, the DC method could reduce slow central processing unit (CPU)-GPU memory access, as well as computational cost, and avoid shortfalls of GPU memory. Numerical applications confirmed that the present code on GPU significantly accelerated the TDDFTB calculations, while maintaining accuracy. Furthermore, the DC-TDDFTB simulation of 2-acetylindan-1,3-dione displays excited-state intramolecular proton transfer and provides reasonable absorption and fluorescence energies with the corresponding experimental values. © 2019 Wiley Periodicals, Inc.

Squeezing All Elements in the Periodic Table: Electron Configuration and Electronegativity of the Atoms under Compression

We present a quantum mechanical model capable of describing isotropic compression of single atoms in a non-reactive neon-like environment. Studies of 93 atoms predict drastic changes to ground-state electronic configurations and electronegativity in the pressure range of 0-300 GPa. This extension of atomic reference data assists in the working of chemical intuition at extreme pressure and can act as a guide to both experiments and computational efforts. For example, we can speculate on the existence of pressure-induced polarity (red-ox) inversions in various alloys. Our study confirms that the filling of energy levels in compressed atoms more closely follows the hydrogenic aufbau principle, where the ordering is determined by the principal quantum number. In contrast, the Madelung energy ordering rule is not predictive for atoms under compression. Magnetism may increase or decrease with pressure, depending on which atom is considered. However, Hund's rule is never violated for single atoms in the considered pressure range. Important (and understandable) electron shifts, s→p, s→d, s→f, and d→f are essential chemical and physical consequences of compression. Among the specific intriguing changes predicted are an increase in the range between the most and least electronegative elements with compression; a rearrangement of electronegativities of the alkali metals with pressure, with Na becoming the most electropositive s¹ element (while Li becomes a p group element and K and heavier become transition metals); phase transitions in Ca, Sr, and Ba correlating well with s→d transitions; spin-reduction in all d-block atoms for which the valence d-shell occupation is d ⁿ (4 ≤ n ≤ 8); d→f transitions in Ce, Dy, and Cm causing Ce to become the most electropositive element of the f-block; f→d transitions in Ho, Dy, and Tb and a s→f transition in Pu. At high pressure Sc and Ti become the most electropositive elements, while Ne, He, and F remain the most electronegative ones.