Core Binding Energy Calculations: A Scalable Approach with the Quantum Embedding-Based Equation-of-Motion Coupled-Cluster Method.

We investigated the use of density matrix embedding theory to facilitate the computation of core ionization energies (IPs) of large molecules at the equation-of-motion coupled-cluster singles doubles with perturbative triples (EOM-CCSD*) level in combination with the core-valence separation (CVS) approximation. The unembedded IP-CVS-EOM-CCSD* method with a triple-ζ basis set produced ionization energies within 1 eV of experiment with a standard deviation of ∼0.2 eV for the core65 data set. The embedded variant contributed very little systematic error relative to the unembedded method, with a mean unsigned error of 0.07 eV and a standard deviation of ∼0.1 eV, in exchange for accelerating the calculations by many orders of magnitude. By employing embedded EOM-CC methods, we computed the core ionization energies of the uracil hexamer, doped fullerene, and chlorophyll molecule, utilizing up to ∼4000 basis functions within 1 eV from experimental values. Such calculations are not currently possible with the unembedded EOM-CC method.

Aliovalent Substitution Tunes Physical Properties in a Conductive Bis(dithiolene) Two-Dimensional Metal-Organic Framework.

Two-dimensional conductive metal-organic frameworks have emerged as promising electronic materials for applications in (opto)electronic, thermoelectric, magnetic, electrocatalytic, and energy storage devices. Many bottom-up or postsynthetic protocols have been developed to isolate these materials or further modulate their electronic properties. However, some methodologies commonly used in classic semiconductors, notably, aliovalent substitution, are conspicuously absent. Here, we demonstrate how aliovalent Fe(III) to Ni(II) substitution enables the isolation of a Ni bis(dithiolene) material from a previously reported Fe analogue. Detailed characterization supports the idea that aliovalent substitution of Fe(III) to Ni(II) results in an in situ oxidation of the organic dithiolene linker. This substitution-induced redox tuning modulates the electronic properties in the system, leading to higher electrical conductivity and Hall mobility but slightly lower carrier densities and weaker antiferromagnetic interactions. Moreover, this aliovalent substitution improves the material's electrochemical stability and thus enables pseudocapacitive behavior in the Ni material. These results demonstrate how classic aliovalent substitution strategies in semiconductors can also be leveraged in conductive MOFs and add further support to this class of compounds as emerging electronic materials.

State-specific frozen natural orbital for reduced-cost algebraic diagrammatic construction calculations: The application to ionization problem.

We have developed a reduced-cost algebraic diagrammatic construction (ADC) method based on state-specific frozen natural orbital and natural auxiliary functions. The newly developed method has been benchmarked on the GW100 test set for the ionization problem. The use of state-specific natural orbitals drastically reduces the size of the virtual space with a systematically controllable accuracy and offers a significant speedup over the standard ionization potential (IP)-ADC(3) method. The accuracy of the method can be controlled by two thresholds and nearly a black box to use. The inclusion of the perturbative correction significantly improves the accuracy of the calculated IP values, and the efficiency of the method has been demonstrated by calculating the IP of a molecule with 60 atoms and more than 2216 basis functions.

A reduced cost four-component relativistic coupled cluster method based on natural spinors.

We present the theory, implementation, and benchmark results for a frozen natural spinors based reduced cost four-component relativistic coupled cluster method. The natural spinors are obtained by diagonalizing the one-body reduced density matrix from a relativistic second-order Møller-Plesset calculation based on a four-component Dirac-Coulomb Hamiltonian. The correlation energy in the coupled cluster method converges more rapidly with respect to the size of the virtual space in the frozen natural spinor basis than that observed in the standard canonical spinors obtained from the Dirac-Hartree-Fock calculation. The convergence of properties is not smooth in the frozen natural spinor basis. However, the inclusion of the perturbative correction smoothens the convergence of the properties with respect to the size of the virtual space in the frozen natural spinor basis and greatly reduces the truncation errors in both energy and property calculations. The accuracy of the frozen natural spinor based coupled cluster methods can be controlled by a single threshold and is a black box to use.