A Heterogeneous CPU + GPU Algorithm for Variational Two-Electron Reduced-Density Matrix Driven Complete Active Space Self-Consistent Field Theory

<div><div><div><p>We present a heterogeneous CPU+GPU algorithm for the direct variational optimization of the two-electron reduced-density matrix (2RDM) under two-particle N-representability conditions. This variational 2RDM (v2RDM) approach is the driver for a polynomially-scaling approximation to configuration-interaction-driven complete active space self-consistent field (CASSCF) theory. For v2RDM-based CASSCF com- putations involving an active space consisting of 50 electrons in 50 orbitals [denoted (50e,50o)], we observe a speedup of a factor of 3.7 when the code is executed on a combination of an NVIDIA TITAN V GPU and an Intel Core i7-6850k CPU, relative to the case when the code is executed on the CPU alone. We use this GPU-accelerated v2RDM-CASSCF algorithm to explore the electronic structure of the 3,k-circumacene and 3,k-periacene series (k=2–7) and compare indicators of polyradical character in the lowest-energy singlet states to those observed for oligoacene molecules. The singlet states in larger circumacene and periacene molecules display the same polyradical characteristics observed in oligoacenes, with the onset of this behavior occuring at smallest k for periacenes, followed by the circumacenes and then the oligoacenes. However, the unpaired electron density that accumulates along the zig-zag edge of the circumacenes is slightly less than that which accumulates in the oligoacenes, while periacenes clearly exhibit the greatest build-up of unpaired electron density in this region.</p></div></div></div>

A Heterogeneous CPU + GPU Algorithm for Variational Two-Electron Reduced-Density Matrix Driven Complete Active Space Self-Consistent Field Theory

<div><div><div><p>We present a heterogeneous CPU+GPU algorithm for the direct variational optimization of the two-electron reduced-density matrix (2RDM) under two-particle N-representability conditions. This variational 2RDM (v2RDM) approach is the driver for a polynomially-scaling approximation to configuration-interaction-driven complete active space self-consistent field (CASSCF) theory. For v2RDM-based CASSCF com- putations involving an active space consisting of 50 electrons in 50 orbitals [denoted (50e,50o)], we observe a speedup of a factor of 3.7 when the code is executed on a combination of an NVIDIA TITAN V GPU and an Intel Core i7-6850k CPU, relative to the case when the code is executed on the CPU alone. We use this GPU-accelerated v2RDM-CASSCF algorithm to explore the electronic structure of the 3,k-circumacene and 3,k-periacene series (k=2–7) and compare indicators of polyradical character in the lowest-energy singlet states to those observed for oligoacene molecules. The singlet states in larger circumacene and periacene molecules display the same polyradical characteristics observed in oligoacenes, with the onset of this behavior occuring at smallest k for periacenes, followed by the circumacenes and then the oligoacenes. However, the unpaired electron density that accumulates along the zig-zag edge of the circumacenes is slightly less than that which accumulates in the oligoacenes, while periacenes clearly exhibit the greatest build-up of unpaired electron density in this region.</p></div></div></div>

Electronic and Spatial Structures of Water-Soluble Dinitrosyl Iron Complexes with Thiol-Containing Ligands Underlying Their Ability to Act as Nitric Oxide and Nitrosonium Ion Donors

The ability of mononuclear dinitrosyl iron commplexes (M-DNICs) with thiolate ligands to act as NO donors and to trigger S-nitrosation of thiols can be explain only in the paradigm of the model of the [Fe+(NO+)2] core ({Fe(NO)2}7 according to the Enemark-Feltham classification). Similarly, the {(RS−)2Fe+(NO+)2}+ structure describing the distribution of unpaired electron density in M-DNIC corresponds to the low-spin (S=1/2) state with a d7 electron configuration of the iron atom and predominant localization of the unpaired electron on MO(dz2) and the square planar structure of M-DNIC. On the other side, the formation of molecular orbitals of M-DNIC including orbitals of the iron atom, thiolate and nitrosyl ligands results in a transfer of electron density from sulfur atoms to the iron atom and nitrosyl ligands. Under these conditions, the positive charge on the nitrosyl ligands diminishes appreciably, the interaction of the ligands with hydroxyl ions or with thiols slows down and the hydrolysis of nitrosyl ligands and the S-nitrosating effect of the latter are not manifested. Most probably, the S-nitrosating effect of nitrosyl ligands is a result of weak binding of thiolate ligands to the iron atom under conditions favoring destabilization of M-DNIC.