1
|
Tehrani A, Richer M, Heidar-Zadeh F. CuGBasis: High-performance CUDA/Python library for efficient computation of quantum chemistry density-based descriptors for larger systems. J Chem Phys 2024; 161:072501. [PMID: 39158048 DOI: 10.1063/5.0216781] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2024] [Accepted: 06/17/2024] [Indexed: 08/20/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Alireza Tehrani
- Department of Chemistry, Queen's University, Kingston, Ontario K7L-3N6, Canada
| | - Michelle Richer
- Department of Chemistry, Queen's University, Kingston, Ontario K7L-3N6, Canada
| | - Farnaz Heidar-Zadeh
- Department of Chemistry, Queen's University, Kingston, Ontario K7L-3N6, Canada
| |
Collapse
|
2
|
Barrera Y, Anderson JSM. Does the radical GPRI strongly depend on the population scheme? A comparative study to predict radical attack on unsaturated molecules with the radical general-purpose reactivity indicator. J Comput Chem 2024; 45:1152-1159. [PMID: 38299704 DOI: 10.1002/jcc.27314] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2023] [Revised: 07/15/2023] [Accepted: 01/09/2024] [Indexed: 02/02/2024]
Abstract
The reactivity of 22 unsaturated molecules undergoing attack by a methyl radical (⋅CH3) have been elucidated using the condensed radical general-purpose reactivity indicator (condensed radical GPRI) appropriate for relatively nucleophilic or electrophilic molecules. Using the appropriate radical GPRI equation for electrophilic attack or nucleophilic radical attack, seven different population schemes were used to assign the most reactive atoms in each of the 22 molecules. The results show that the condensed radical GPRI is sensitive to the population scheme chosen, but less sensitive than the radical Fukui function. Therefore, the reliability of these methods depends on the population scheme. Our investigation indicates that the condensed radical GPRI is most accurate in predicting the dominant products of the methyl radical addition reactions on a variety of unsaturated molecules when the Hirshfeld, Merz-Singh-Kollman, or Voronoi deformation density population schemes are used. Furthermore, for all populations schemes in the majority of instances where the radical Fukui function failed the radical GPRI was able to identify the most reactive atom under certain reactivity conditions.
Collapse
Affiliation(s)
- Yoshio Barrera
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Mexico City, Mexico
| | - James S M Anderson
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Mexico City, Mexico
| |
Collapse
|
3
|
He X, Li M, Rong C, Zhao D, Liu W, Ayers PW, Liu S. Some Recent Advances in Density-Based Reactivity Theory. J Phys Chem A 2024; 128:1183-1196. [PMID: 38329898 DOI: 10.1021/acs.jpca.3c07997] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
Establishing a chemical reactivity theory in density functional theory (DFT) language has been our intense research interest in the past two decades, exemplified by the determination of steric effect and stereoselectivity, evaluation of electrophilicity and nucleophilicity, identification of strong and weak interactions, and formulation of cooperativity, frustration, and principle of chirality hierarchy. In this Featured Article, we first overview the four density-based frameworks in DFT to appreciate chemical understanding, including conceptual DFT, use of density associated quantities, information-theoretic approach, and orbital-free DFT, and then present a few recent advances of these frameworks as well as new applications from our studies. To that end, we will introduce the relationship among these frameworks, determining the entire spectrum of interactions with Pauli energy derivatives, performing topological analyses with information-theoretic quantities, and extending the density-based frameworks to excited states. Applications to examine physiochemical properties in external electric fields and to evaluate polarizability for proteins and crystals are discussed. A few possible directions for future development are followed, with the special emphasis on its merger with machine learning.
Collapse
Affiliation(s)
- Xin He
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| | - Meng Li
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, China
| | - Chunying Rong
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China), Hunan Normal University, Changsha, Hunan 410081, China
| | - Dongbo Zhao
- Institute of Biomedical Research, Yunnan University, Kunming 650500, China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton ONL8S, Canada
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, North Carolina 27599-3420, United States
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599-3290, United States
| |
Collapse
|
4
|
Miranda-Quintana RA, Vela A, De Proft F, Martínez González M, Gázquez JL. Can we predict ambident regioselectivity using the chemical hardness? Phys Chem Chem Phys 2023; 25:13611-13622. [PMID: 37144347 DOI: 10.1039/d3cp00876b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
The hard/soft acid/base (HSAB) principle is a cornerstone in our understanding of chemical reactivity preferences. Motivated by the success of the original ("global") version of this rule, a "local" counterpart was readily proposed to account for regioselectivity preferences, in particular, in ambident reactions. However, ample experimental evidence indicates that the local HSAB principle often fails to provide meaningful predictions. Here we examine the assumptions behind the standard proof of the local HSAB rule, showing that it is based on a flawed premise. By solving this issue, we show that it is critical to consider not only the charge transferred between the different reacting centers but also the charge reorganization within the non-reacting parts of the molecule. We propose different reorganization models and derive the corresponding regioselectivity rules for each.
Collapse
Affiliation(s)
| | - Alberto Vela
- Departamento de Química, Centro de Investigacion y Estudios Avanzados, Av. Instituto Politécnico Nacional, 2508, Ciudad de, México 07360, Mexico
| | - Frank De Proft
- Eenheid Algemene Chemie (ALGC) Vrije Universiteit Brussel Pleinlaan 2, 1050 Brussels, Belgium
| | - Marco Martínez González
- Department of Chemistry & Chemical Biology, McMaster University, Hamilton, Ontario, L8S 4M1, Canada
| | - José L Gázquez
- Departamento de Química, Universidad Autónoma Metropolitana-Iztapalapa, Ave. San Rafael Atlixco 186, Ciudad de México 09340, Mexico
| |
Collapse
|
5
|
Miranda-Quintana RA, Heidar-Zadeh F, Fias S, Chapman AEA, Liu S, Morell C, Gómez T, Cárdenas C, Ayers PW. Molecular interactions from the density functional theory for chemical reactivity: Interaction chemical potential, hardness, and reactivity principles. Front Chem 2022; 10:929464. [PMID: 35936089 PMCID: PMC9352952 DOI: 10.3389/fchem.2022.929464] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2022] [Accepted: 06/27/2022] [Indexed: 11/13/2022] Open
Abstract
In the first paper of this series, the authors derived an expression for the interaction energy between two reagents in terms of the chemical reactivity indicators that can be derived from density functional perturbation theory. While negative interaction energies can explain reactivity, reactivity is often more simply explained using the “|dμ| big is good” rule or the maximum hardness principle. Expressions for the change in chemical potential (μ) and hardness when two reagents interact are derived. A partial justification for the maximum hardness principle is that the terms that appear in the interaction energy expression often reappear in the expression for the interaction hardness, but with opposite sign.
Collapse
Affiliation(s)
- Ramón Alain Miranda-Quintana
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, FL, United States
- *Correspondence: Ramón Alain Miranda-Quintana, ; Tatiana Gómez, Carlos Cárdenas, ; Paul W. Ayers,
| | | | - Stijn Fias
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON, Canada
| | - Allison E. A. Chapman
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON, Canada
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, NC, United states
| | - Christophe Morell
- Université de Lyon, Universit́e Claude Bernard Lyon 1, Institut des Sciences Analytiques, UMR CNRS 5280, Villeurbanne Cedex, France
| | - Tatiana Gómez
- Theoretical and Computational Chemistry Center, Institute of Applied Chemical Sciences, Faculty of Engineering, Universidad Autonoma de Chile, Santiago, Chile
- *Correspondence: Ramón Alain Miranda-Quintana, ; Tatiana Gómez, Carlos Cárdenas, ; Paul W. Ayers,
| | - Carlos Cárdenas
- Departamento de Fisica, Facultad de Ciencias, Universidad de Chile, Santiago, Chile
- Centro para el desarrollo de la Nanociencias y Nanotecnologia, CEDENNA, Santiago, Chile
- *Correspondence: Ramón Alain Miranda-Quintana, ; Tatiana Gómez, Carlos Cárdenas, ; Paul W. Ayers,
| | - Paul W. Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON, Canada
- *Correspondence: Ramón Alain Miranda-Quintana, ; Tatiana Gómez, Carlos Cárdenas, ; Paul W. Ayers,
| |
Collapse
|