Q-Chem: an engine for innovation

A Vision for the Future of Multiscale Modeling

The rise of modern computer science enabled physical chemistry to make enormous progresses in understanding and harnessing natural and artificial phenomena. Nevertheless, despite the advances achieved over past decades, computational resources are still insufficient to thoroughly simulate extended systems from first principles. Indeed, countless biological, catalytic and photophysical processes require ab initio treatments to be properly described, but the breadth of length and time scales involved makes it practically unfeasible. A way to address these issues is to couple theories and algorithms working at different scales by dividing the system into domains treated at different levels of approximation, ranging from quantum mechanics to classical molecular dynamics, even including continuum electrodynamics. This approach is known as multiscale modeling and its use over the past 60 years has led to remarkable results. Considering the rapid advances in theory, algorithm design, and computing power, we believe multiscale modeling will massively grow into a dominant research methodology in the forthcoming years. Hereby we describe the main approaches developed within its realm, highlighting their achievements and current drawbacks, eventually proposing a plausible direction for future developments considering also the emergence of new computational techniques such as machine learning and quantum computing. We then discuss how advanced multiscale modeling methods could be exploited to address critical scientific challenges, focusing on the simulation of complex light-harvesting processes, such as natural photosynthesis. While doing so, we suggest a cutting-edge computational paradigm consisting in performing simultaneous multiscale calculations on a system allowing the various domains, treated with appropriate accuracy, to move and extend while they properly interact with each other. Although this vision is very ambitious, we believe the quick development of computer science will lead to both massive improvements and widespread use of these techniques, resulting in enormous progresses in physical chemistry and, eventually, in our society.

Two- and one-photon absorption spectra of aqueous thiocyanate anion highlight the role of symmetry in the condensed phase

We present the two-photon absorption (2PA) spectrum of aqueous thiocyanate calculated using high-level quantum-chemistry methods. The 2PA spectrum is compared to the one-photon absorption (1PA) spectrum computed using the same computational protocol. Although the two spectra probe the same set of electronic states, the intensity patterns are different, leading to an apparent red-shift of the 2PA spectrum relative to the 1PA spectrum. The presented analysis explains the intensity patterns and attributes the differences between the 1PA and 2PA spectra to the native symmetry of isolated SCN  - , which influences the spectra in the low-symmetry solvated environment. The native symmetry also manifests itself in variations of the polarization ratio (e.g., parallel vs. perpendicular cross sections) across the spectrum. The presented results highlight the potential of 2PA spectroscopy and high-level quantum-chemistry methods in studies of condensed-phase phenomena.

Revisiting the benzene excimer using [2,2] paracyclophane model system: Experiment and theory

We report high-level calculations of the excited states of [2,2]-paracyclophane (PCP), which was recently investigated experimentally by ultrafast pump-probe experiments on oriented single crystals [Haggag et al., ChemPhotoChem 6 e202200181 (2022)]. PCP, in which the orientation of the two benzene rings and their range of motion are constrained, serves as a model for studying benzene excimer formation. The character of the excimer state and the state responsible for the brightest transition are similar to those of the benzene dimer. The constrained structure of PCP allows one to focus on the most important degree of freedom, the inter-ring distance. The calculations explain the main features of the transient absorption spectral evolution. This brightest transition of the excimer is polarized along the inter-fragment axis. The absorption of the light polarized in the plane of the rings reveals the presence of other absorbing states of Rydberg character, with much weaker intensities. We also report new transient absorption data obtained by a broadband 8 fs pump, which time-resolve strong modulations of the excimer absorption. The combination of theory and experiment provides a detailed picture of the evolution of the electronic structure of the PCP excimer in the course of a single molecular vibration.

Membrane Permeation of Psychedelic Tryptamines by Dynamic Simulations

Renewed scientific interest in psychedelic compounds represents one of the most promising avenues for addressing the current burden of mental health disorders. Classic psychedelics are a group of compounds that exhibit structural similarities to the naturally occurring neurotransmitter serotonin (5-HT). Acting on the 5-HT type 2A receptors (HT2ARs), psychedelics induce enduring neurophysiological changes that parallel their therapeutic psychological and behavioral effects. Recent preclinical evidence suggests that the ability of psychedelics to exert their action is determined by their ability to permeate the neuronal membrane to target a pool of intracellular 5-HT2ARs. In this computational study, we employ classical molecular dynamics simulations and umbrella sampling techniques to investigate the permeation behavior of 12 selected tryptamines and to characterize the interactions that drive the process. We aim at elucidating the impact of N-alkylation, indole ring substitution and positional modifications, and protonation on their membrane permeability. Dimethylation of the primary amine group and the introduction of a methoxy group at position 5 exhibited an increase in permeability. Moreover, there is a significant influence of positional substitutions on the indole groups, and the protonation of the molecules substantially increases the energy barrier at the center of the bilayer, making the compounds highly impermeable. All the information extracted from the trends predicted by the simulations can be applied in future drug design projects to develop psychedelics with enhanced activity.

State-Interaction Approach for Evaluating g-Tensors within EOM-CC and RAS-CI Frameworks: Theory and Benchmarks

Among various techniques designed for studying open-shell species, electron paramagnetic resonance (EPR) spectroscopy plays an important role. The key quantity measured by EPR is the g-tensor, describing the coupling between an external magnetic field and molecular electronic spin. One theoretical framework for quantum chemistry calculations of g-tensors is based on response theory, which involves substantial developments that are specific to the underlying electronic structure models. A simplified and easier-to-implement approach is based on the state-interaction scheme, in which perturbation is included by considering a small number of states. We describe and benchmark the state-interaction approach using equation-of-motion coupled-cluster and restricted-active-space configuration interaction wave functions. The analysis confirms that this approach can deliver accurate results and highlights caveats of applying it, such as a choice of the reference state, convergence with respect to the number of states used in calculations, etc. The analysis also contributes toward a better understanding of challenges in calculations of higher-order properties using approximate wave functions.

A Theoretical Study on the Underlying Factors of the Difference in Performance of Organic Solar Cells Based on ITIC and Its Isomers

Recently, non-fullerene-based organic solar cells (OSCs) have made great breakthroughs, and small structural differences can have dramatic impacts on the power conversion efficiency (PCE). We take ITIC and its isomers as examples to study their effects on the performance of OSCs. ITIC and NFBDT only differed in the side chain position, and they were used as models with the same donor molecule, PBDB-T, to investigate the main reasons for the difference in their performance in terms of theoretical methods. In this work, a detailed comparative analysis of the electronic structure, absorption spectra, open circuit voltage and interfacial parameters of the ITIC and NFBDT systems was performed mainly by combining the density functional theory/time-dependent density functional theory and molecular dynamics simulations. The results showed that the lowest excited state of the ITIC molecule possessed a larger ∆q and more hybrid FE/CT states, and PBDB-T/ITIC had more charge separation paths as well as a larger kCS and smaller kCR. The reason for the performance difference between PBDB-T/ITIC and PBDB-T/NFBDT was elucidated, suggesting that ITIC is a superior acceptor based on a slight modulation of the side chain and providing a guiding direction for the design of superior-performing small molecule acceptor materials.

Molecular-Orbital Framework of Two-Electron Processes: Application to Auger and Intermolecular Coulomb Decay

States with core- or inner-shell vacancies, which are commonly created by absorption of high-energy photons, can decay by a two-electron process in which one electron fills the core hole and the second one is ejected. These processes accompany many X-ray spectroscopies. Depending on the nature of the initial core- or inner-shell-hole state and the decay valence-hole state, these processes are called Auger decay, intermolecular Coulomb decay, or electron-transfer-mediated decay. To connect many-body wave functions of the initial and final states with the molecular orbital picture of the decay, we introduce the concept of natural Auger orbitals (NAOs). NAOs are obtained by a two-step singular value decomposition of the two-body Dyson orbitals, reduced quantities that enter the expression of the decay rate in the Feshbach-Fano treatment. NAOs afford chemical insight and interpretation of the high-level ab initio calculations of Auger decay and related two-electron relaxation processes.

Multiscale Simulations of the Covalent Inhibition of the SARS-CoV-2 Main Protease: Four Compounds and Three Reaction Mechanisms

We report the results of computational modeling of the reactions of the SARS-CoV-2 main protease (MPro) with four potential covalent inhibitors. Two of them, carmofur and nirmatrelvir, have shown experimentally the ability to inhibit MPro. Two other compounds, X77A and X77C, were designed computationally in this work. They were derived from the structure of X77, a non-covalent inhibitor forming a tight surface complex with MPro. We modified the X77 structure by introducing warheads capable of reacting with the catalytic cysteine residue in the MPro active site. The reaction mechanisms of the four molecules with MPro were investigated by quantum mechanics/molecular mechanics (QM/MM) simulations. The results show that all four compounds form covalent adducts with the catalytic cysteine Cys 145 of MPro. From the chemical perspective, the reactions of these four molecules with MPro follow three distinct mechanisms. The reactions are initiated by a nucleophilic attack of the thiolate group of the deprotonated cysteine residue from the catalytic dyad Cys145-His41 of MPro. In the case of carmofur and X77A, the covalent binding of the thiolate to the ligand is accompanied by the formation of the fluoro-uracil leaving group. The reaction with X77C follows the nucleophilic aromatic substitution SNAr mechanism. The reaction of MPro with nirmatrelvir (which has a reactive nitrile group) leads to the formation of a covalent thioimidate adduct with the thiolate of the Cys145 residue in the enzyme active site. Our results contribute to the ongoing search for efficient inhibitors of the SARS-CoV-2 enzymes.

Quantum Mechanical Prediction of Dissociation Constants for Thiazol-2-imine Derivatives

As weak acids or bases, in solution, drug molecules are in either their ionized or nonionized states. A high degree of ionization is essential for good water solubility of a drug molecule and is required for drug-receptor interactions, whereas the nonionized form improves a drug's lipophilicity, allowing the ligand to cross the cell membrane. The penetration of a drug ligand through cell membranes is mainly governed by the pKa of the drug molecule and the membrane environment. In this study, with the aim of predicting the acetonitrile pKa's (pKa(MeCN)) of eight drug-like thiazol-2-imine derivatives, we propose a very accurate and computationally affordable protocol by using several quantum mechanical approaches. Benchmark studies were conducted on a set of training molecules, which were selected from the literature with known pKa(water) and pKa(MeCN). Highly well-correlated pKa values were obtained when the calculations were performed with the isodesmic method at the M062X/6-31G** level of theory in conjunction with SMD solvation model for nitrogen-containing heterocycles. Finally, experimentally unknown pKa(MeCN) values of eight thiazol-2-imine structures, which were previously synthesized by some of us, are proposed.

Influence of Transition Metal Electron Configuration on the Structure of Metal-EDTA Complexes

The vibrational spectra of cold complexes of ethylenediaminetetraacetic acid (EDTA) with transition metal dications in vacuo show how the electronic structure of the metal provides a geometric template for interaction with the functional groups of the binding pocket. The OCO stretching modes of the carboxylate groups of EDTA serve as structural probes, informing on the spin state of the ion as well as the coordination number in the complex. The results highlight the flexibility of EDTA in accepting a large range of metal cations in its binding site.

The Auger spectrum of benzene

We present an ab initio computational study of the Auger electron spectrum of benzene. Auger electron spectroscopy exploits the Auger-Meitner effect, and although it is established as an analytic technique, the theoretical modeling of molecular Auger spectra from first principles remains challenging. Here, we use coupled-cluster theory and equation-of-motion coupled-cluster theory combined with two approaches to describe the decaying nature of core-ionized states: (i) Feshbach-Fano resonance theory and (ii) the method of complex basis functions. The spectra computed with these two approaches are in excellent agreement with each other and also agree well with experimental Auger spectra of benzene. The Auger spectrum of benzene features two well-resolved peaks at Auger electron energies above 260 eV, which correspond to final states with two electrons removed from the 1e_1g and 3e_2g highest occupied molecular orbitals. At lower Auger electron energies, the spectrum is less well resolved, and the peaks comprise multiple final states of the benzene dication. In line with theoretical considerations, singlet decay channels contribute more to the total Auger intensity than the corresponding triplet decay channels.

XPS and quantum chemical analysis of 4Me-BODIPY derivatives

The results of a X-ray photoelectron spectroscopy (XPS) and steady-state absorption spectroscopy study of the electronic structure, and cationic and excited states of a series of 1,3,5,7-tetramethyl-substituted BODIPYs (4Me,2R-BODIPYs) are presented. The experimental data were interpreted using high-level ab initio quantum chemical computations, including the algebraic diagrammatic construction method for the polarization propagator of the second order (ADC(2)), the outer-valence Green's function (OVGF) method, the density functional (DFT) approach, and the time-dependent DFT (TD-DFT) approach. Substitution effects on the XPS and absorption spectra were determined for 2,6-positions of 4Me,2R-BODIPY pyrrole nuclei (R = H, Br, Bu, benzyl). A very satisfactory performance of the DFT Koopmans theorem analogue was demonstrated with respect to the energy intervals between the electronic levels of 4Me,2R-BODIPY above 13 eV (BHHLYP functional) and the values of the HOMO-LUMO energy gap (ωB97X functional).

MoBioTools: A toolkit to setup quantum mechanics/molecular mechanics calculations

We present a toolkit that allows for the preparation of QM/MM input files from a conformational ensemble of molecular geometries. The package is currently compatible with trajectory and topology files in Amber, CHARMM, GROMACS and NAMD formats, and has the possibility to generate QM/MM input files for Gaussian (09 and 16), Orca (≥4.0), NWChem and (Open)Molcas. The toolkit can be used in command line, so that no programming experience is required, although it presents some features that can also be employed as a python application programming interface. We apply the toolkit in four situations in which different electronic-structure properties of organic molecules in the presence of a solvent or a complex biological environment are computed: the reduction potential of the nucleobases in acetonitrile, an energy decomposition analysis of tyrosine interacting with water, the absorption spectrum of an azobenzene derivative integrated into a voltage-gated ion channel, and the absorption and emission spectra of the luciferine/luciferase complex. These examples show that the toolkit can be employed in a manifold of situations for both the electronic ground state and electronically excited states. It also allows for the automatic correction of the active space in the case of CASSCF calculations on an ensemble of geometries, as it is shown for the azobenzene derivative photoswitch case.

Magnetic exchange interactions in binuclear and tetranuclear iron(III) complexes described by spin-flip DFT and Heisenberg effective Hamiltonians

Low-energy spectra of single-molecule magnets (SMMs) are often described by Heisenberg Hamiltonians. Within this formalism, exchange interactions between magnetic centers determine the ground-state multiplicity and energy separation between the ground and excited states. In this contribution, we extract exchange coupling constants (J) for a set of iron (III) binuclear and tetranuclear complexes from all-electron calculations using non-collinear spin-flip time-dependent density functional theory (NC-SF-TDDFT). For 12 binuclear complexes with J-values ranging from -6 to -132 cm^-1 , our benchmark calculations using the short-range hybrid ωPBEh functional and 6-31G(d,p) basis set agree well with the experimentally derived values (mean absolute error of 4.7 cm^-1 ). For the tetranuclear SMMs, the computed J constants are within 6 cm^-1 from the experimentally derived values. We explore the range of applicability of the Heisenberg model by analyzing bonding patterns in these Fe(III) complexes using natural orbitals (NO), their occupations, and the number of effectively unpaired electrons. The results illustrate the efficiency of the spin-flip protocol for computing the exchange couplings and the utility of the NO analysis in assessing the validity of effective spin Hamiltonians.

Improving Lurasidone Hydrochloride's Solubility and Stability by Higher-Order Complex Formation with Hydroxypropyl-β-cyclodextrin

The biopharmaceutical classification system groups low-solubility drugs into two groups: II and IV, with high and low permeability, respectively. Most of the new drugs developed for common pathologies present solubility issues. This is the case of lurasidone hydrochloride-a drug used for the treatment of schizophrenia and bipolar depression. Likewise, the stability problems of some drugs limit the possibility of preparing them in liquid pharmaceutical forms where hydrolysis and oxidation reactions can be favored. Lurasidone hydrochloride presents the isoindole-1,3-dione ring, which is highly susceptible to alkaline hydrolysis, and the benzisothiazole ring, which is susceptible to a lesser extent to oxidation. Herein, we propose to study the increase in the solubility and stability of lurasidone hydrochloride by the formation of higher-order inclusion complexes with hydroxypropyl-β-cyclodextrin. Several stoichiometric relationships were studied at between 0.5 and 3 hydroxypropyl-β-cyclodextrin molecules per drug molecule. The obtained products were characterized, and their solubility and stability were assessed. According to the obtained results, the formation of inclusion complexes dramatically increased the solubility of the drug, and this increased with the increase in the inclusion ratio. This was associated with the loss of crystalline state of the drug, which was in an amorphous state according to infrared spectroscopy, calorimetry, and X-ray analysis. This was also correlated with the stabilization of lurasidone by the cyclodextrin inhibiting its recrystallization. Phase solubility,¹H-NMR, and docking computational characterization suggested that the main stoichiometric ratio was 1:1; however, we cannot rule out a 1:2 ratio, where a second cyclodextrin molecule could bind through the isoindole-1,3-dione ring, improving its stability as well. Finally, we can conclude that the formation of higher-order inclusion complexes of lurasidone with hydroxypropyl-β-cyclodextrin is a successful strategy to increase the solubility and stability of the drug.

Modeling Light-Induced Chromophore Hydration in the Reversibly Photoswitchable Fluorescent Protein Dreiklang

We report the results of a computational study of the mechanism of the light-induced chemical reaction of chromophore hydration in the fluorescent protein Dreiklang, responsible for its switching from the fluorescent ON-state to the dark OFF-state. We explore the relief of the charge-transfer excited-state potential energy surface in the ON-state to locate minimum energy conical intersection points with the ground-state energy surface. Simulations of the further evolution of model systems allow us to characterize the ground-state reaction intermediate tentatively suggested in the femtosecond studies of the light-induced dynamics in Dreiklang and finally to arrive at the reaction product. The obtained results clarify the details of the photoswitching mechanism in Dreiklang, which is governed by the chemical modification of its chromophore.

Atomic Group Decomposition of Charge Transfer Excitation Global Indexes

A model for decomposing the Le Bahers, Adamo, and Ciofini Charge Transfer (CT) Excitations global indexes ( J. Chem. Theory Comput. 2011, 7, 2498-2506) into molecular subdomains contributions is presented and a software, DOCTRINE (atomic group Decomposition Of the Charge TRansfer INdExes) for the implementation of this novel model has been coded. Although our method applies to any fuzzy or to any disjoint exhaustive partitioning of the real space, it is here applied using a definition of chemically relevant molecular subdomains based on the Atoms in Molecules Bader basins. This choice has the relevant advantage of associating intra or inter subdomain contributions to rigorously defined quantum objects, yet bearing a clear chemical meaning. Our method allows for a quantitative evaluation of the subdomain contributions to the charge transfer, the charge transfer excitation length and the dipole moment change upon excitation. All these global indexes may be obtained either from the electron density increment or the electron density depletion upon excitation. However, the subdomain contributions obtained from the two distributions generally differ, therefore allowing to distinguish whether the contribution to a given property of a given subdomain is dominated by one of the two distributions or if both are playing a significant role. As a toy system for the first application of our model, a typical [D-π-A, π = conjugated bridge] compound belonging to the merocyanine dyes family is selected, and the first four excited states of this compound in a strongly polar protic solvent and in a weakly polar solvent are thoroughly investigated.

How Reproducible Are QM/MM Simulations? Lessons from Computational Studies of the Covalent Inhibition of the SARS-CoV-2 Main Protease by Carmofur

This work explores the level of transparency in reporting the details of computational protocols that is required for practical reproducibility of quantum mechanics/molecular mechanics (QM/MM) simulations. Using the reaction of an essential SARS-CoV-2 enzyme (the main protease) with a covalent inhibitor (carmofur) as a test case of chemical reactions in biomolecules, we carried out QM/MM calculations to determine the structures and energies of the reactants, the product, and the transition state/intermediate using analogous QM/MM models implemented in two software packages, NWChem and Q-Chem. Our main benchmarking goal was to reproduce the key energetics computed with the two packages. Our results indicate that quantitative agreement (within the numerical thresholds used in calculations) is difficult to achieve. We show that rather minor details of QM/MM simulations must be reported in order to ensure the reproducibility of the results and offer suggestions toward developing practical guidelines for reporting the results of biosimulations.

[2.2.2.2]Paracyclophanetetraenes (PCTs): cyclic structural analogues of poly( p‑phenylene vinylene)s (PPVs)

Background: Poly( p-phenylene vinylene)s ( PPVs) and [2.2.2.2]paracyclophanetetraene ( PCT) are both composed of alternating π-conjugated para-phenylene and vinylene units. However, while the former constitute a class of π-conjugated polymers that has been used in organic electronics for decades, the latter is a macrocycle that only recently revealed its potential for applications such as organic battery electrodes. The cyclic structure endows PCT with unusual properties, and further tuning of these may be required for specific applications. Methods: In this article, we adopt an approach often used for tuning the properties of PPVs, the introduction of alkoxy (or alkylthio) substituents at the phenylene units, for tuning the optoelectronic properties of PCT. The resulting methoxy- and methylthio-substituted PCTs, obtained by Wittig cyclisation reactions, are studied by UV-vis absorption, photoluminescence, and cyclic voltammetry measurements, and investigated computationally using the visualisation of chemical shielding tensors (VIST) method. Results: The measurements show that substitution leads to slight changes in terms of absorption/emission energies and redox potentials while having a pronounced effect on the photoluminescence intensity. The computations show the effect of the substituents on the ring currents and chemical shielding and on the associated local and global (anti)aromaticity of the macrocycles, highlighting the interplay of local and global aromaticity in various electronic states. Conclusions: The study offers interesting insights into the tuneability of the properties of this versatile class of π-conjugated macrocycles.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Functional group introduction and aromatic unit variation in a set of π-conjugated macrocycles: revealing the central role of local and global aromaticity

π-Conjugated macrocycles are molecules with unique properties that are increasingly exploited for applications and the question of whether they can sustain global aromatic or antiaromatic ring currents is particularly intriguing. However, there are only a small number of experimental studies that investigate how the properties of π-conjugated macrocycles evolve with systematic structural changes. Here, we present such a systematic experimental study of a set of [2.2.2.2]cyclophanetetraenes, all with formally Hückel antiaromatic ground states, and combine it with an in-depth computational analysis. The study reveals the central role of local and global aromaticity for rationalizing the observed optoelectronic properties, ranging from extremely large Stokes shifts of up to 1.6 eV to reversible fourfold reduction, a highly useful feature for charge storage/accumulation applications. A recently developed method for the visualization of chemical shielding tensors (VIST) is applied to provide unique insight into local and global ring currents occurring in different planes along the macrocycle. Conformational changes as a result of the structural variations can further explain some of the observations. The study contributes to the development of structure-property relationships and molecular design guidelines and will help to understand, rationalize, and predict the properties of other π-conjugated macrocycles. It will also assist in the design of macrocycle-based supramolecular elements with defined properties.

Equation-of-Motion Coupled-Cluster Protocol for Calculating Magnetic Properties: Theory and Applications to Single-Molecule Magnets

We present a new computational protocol for computing macroscopic magnetic properties of transition-metal complexes using the equation-of-motion coupled-cluster (EOM-CC) framework. The approach follows a two-step state-interaction scheme: we first compute zero-order states using nonrelativistic EOM-CC and then use these states to evaluate matrix elements of the spin-orbit and Zeeman operators. Diagonalization of the resulting Hamiltonian yields spin-orbit- and field-perturbed eigenstates. Temperature- and field-dependent magnetization and susceptibility are computed by numerical differentiation of the partition function. To compare with powder-sample experiments, these quantities are numerically averaged over field orientations. We applied this protocol to several single-molecule magnets (SMMs) with Fe(II) and Fe(III) in trigonal pyramidal, linear, and trigonal bipyramidal coordination environments. We described the underlying electronic structure by the electron-attachment (EOM-EA) and spin-flip (EOM-SF) variants of EOM-CC. The computed energy barriers for spin inversion, and macroscopic magnetization and susceptibility agree well with experimental data. Trends in magnetic anisotropy and spin-reversal energy barriers are explained in terms of a molecular orbital picture rigorously distilled from spinless transition density matrices between many-body states. The results illustrate excellent performances of EOM-CC in describing magnetic behavior of mononuclear transition-metal SMMs.

Modified Indulines: From Dyestuffs to In Vivo Theranostic Agents

The efficient, versatile, and straightforward synthesis of the first N-alkyl analogues of induline 3B (8a and 8b) is reported. Thanks to the introduction of lipophilic substituents and their attractive photophysical properties (far-red emission and production of singlet oxygen), phenazinium 8b can be used as a theranostic agent and shows, at very low concentrations (100 nM), a remarkable ability to (i) image cells and zebrafish embryos with high quality under both mono- (514 nm) and biphotonic (790 and 810 nm) excitations, (ii) efficiently and quickly penetrate cancer cells rather than healthy fibroblasts, and (iii) induce a total or almost total cancer cell death in vitro and in vivo after illumination (λ_exc = 540-560 nm). The molecular structure of 8b is based on a triamino-phenazinium core only, with no need for additional components, highlighting the emergence of a minimalistic and versatile class of fluorescent probes for targeted photodynamic cancer therapy.

The orbital picture of the first dipole hyperpolarizability from many-body response theory

We present an approach for obtaining a molecular orbital picture of the first dipole hyperpolarizability (β) from correlated many-body electronic structure methods. Ab initio calculations of β rely on quadratic response theory, which recasts the sum-over-all-states expression of β into a closed-form expression by calculating a handful of first- and second-order response states; for resonantly enhanced β, damped response theory is used. These response states are then used to construct second-order response reduced one-particle density matrices (1PDMs), which, upon visualization in terms of natural orbitals (NOs), facilitate a rigorous and black-box mapping of the underlying electronic structure with β. We explain the interpretation of different components of the response 1PDMs and the corresponding NOs within both the undamped and damped response theory framework. We illustrate the utility of this new tool by deconstructing β for cis-difluoroethene, para-nitroaniline, and hemibonded OH^· + H₂O complex, computed within the framework of coupled-cluster singles and doubles response theory, in terms of the underlying response 1PDMs and NOs for a range of frequencies.

Exploitation of Baird Aromaticity and Clar’s Rule for Tuning the Triplet Energies of Polycyclic Aromatic Hydrocarbons

Polycyclic aromatic hydrocarbons (PAH) are a prominent substance class with a variety of applications in molecular materials science. Their electronic properties crucially depend on the bond topology in ways that are often highly non-intuitive. Here, we study, using density functional theory, the triplet states of four biphenylene-derived PAHs finding dramatically different triplet excitation energies for closely related isomeric structures. These differences are rationalised using a qualitative description of Clar sextets and Baird quartets, quantified in terms of nucleus independent chemical shifts, and represented graphically through a recently developed method for visualising chemical shielding tensors (VIST). The results are further interpreted in terms of a 2D rigid rotor model of aromaticity and through an analysis of the natural transition orbitals involved in the triplet excited states showing good consistency between the different viewpoints. We believe that this work constitutes an important step in consolidating these varying viewpoints of electronically excited states.

Equation-of-motion coupled-cluster method with double electron-attaching operators: Theory, implementation, and benchmarks

We report a production-level implementation of the equation-of-motion (EOM) coupled-cluster (CC) method with double electron-attaching (DEA) EOM operators of 2p and 3p1h types, EOM-DEA-CCSD. This ansatz, suitable for treating electronic structure patterns that can be described as two-electrons-in-many orbitals, represents a useful addition to the EOM-CC family of methods. We analyze the performance of EOM-DEA-CCSD for energy differences and molecular properties. By considering reduced quantities, such as state and transition one-particle density matrices, we compare EOM-DEA-CCSD wave functions with wave functions computed by other EOM-CCSD methods. The benchmarks illustrate that EOM-DEA-CCSD is capable of treating diradicals, bond-breaking, and some types of conical intersections.

Feshbach-Fano approach for calculation of Auger decay rates using equation-of-motion coupled-cluster wave functions. I. Theory and implementation

X-ray absorption creates electron vacancies in the core shell. These highly excited states often relax by Auger decay-an autoionization process in which one valence electron fills the core hole and another valence electron is ejected into the ionization continuum. Despite the important role of Auger processes in many experimental settings, their first-principles modeling is challenging, even for small systems. The difficulty stems from the need to describe many-electron continuum (unbound) states, which cannot be tackled with standard quantum-chemistry methods. We present a novel approach to calculate Auger decay rates by combining Feshbach-Fano resonance theory with the equation-of-motion coupled-cluster single double (EOM-CCSD) framework. We use the core-valence separation scheme to define projectors into the bound (square-integrable) and unbound (continuum) subspaces of the full function space. The continuum many-body decay states are represented by products of an appropriate EOM-CCSD state and a free-electron state, described by a continuum orbital. The Auger rates are expressed in terms of reduced quantities, two-body Dyson amplitudes (objects analogous to the two-particle transition density matrix), contracted with two-electron bound-continuum integrals. Here, we consider two approximate treatments of the free electron: a plane wave and a Coulomb wave with an effective charge, which allow us to evaluate all requisite integrals analytically; however, the theory can be extended to incorporate a more sophisticated description of the continuum orbital.

Frontiers in Multiscale Modeling of Photoreceptor Proteins

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

Spectroscopy and Scattering Studies Using Interpolated Ab Initio Potentials

The Born-Oppenheimer potential energy surface (PES) has come a long way since its introduction in the 1920s, both conceptually and in predictive power for practical applications. Nevertheless, nearly 100 years later-despite astonishing advances in computational power-the state-of-the-art first-principles prediction of observables related to spectroscopy and scattering dynamics is surprisingly limited. For example, the water dimer, (H₂O)₂, with only six nuclei and 20 electrons, still presents a formidable challenge for full-dimensional variational calculations of bound states and is considered out of reach for rigorous scattering calculations. The extremely poor scaling of the most rigorous quantum methods is fundamental; however, recent progress in development of approximate methodologies has opened the door to fairly routine high-quality predictions, unthinkable 20 years ago. In this review, in relation to the workflow of spectroscopy and/or scattering studies, we summarize progress and challenges in the component areas of electronic structure calculations, PES fitting, and quantum dynamical calculations.

Interplay between Locally Excited and Charge Transfer States Governs the Photoswitching Mechanism in the Fluorescent Protein Dreiklang

We present the results of high-level electronic structure and dynamics simulations of the photoactive protein Dreiklang. With the goal of understanding the details of the Dreiklang photocycle, we carefully characterize the excited states of the ON- and OFF-forms of Dreiklang. The key finding of our study is the existence of a low-lying excited state of a charge-transfer character in the neutral ON form and that population of this state, which is nearly isoenergetic with the locally excited bright state, initiates a series of steps that ultimately lead to the formation of the hydrated dark chromophore (OFF state). These results allow us to refine the mechanistic picture of Dreiklang's photocycle and photoactivation.

From Y6 to BTPT-4F: a theoretical insight into the influence of the individual change of fused-ring skeleton length or side alkyl chains on molecular arrangements and electron mobility

Organic solar cells (OSCs) based on non-fullerene acceptor (NFA) Y6 have drawn tremendous attention due to the great progress in their power conversion efficiencies (PCEs).

Data of interaction of supported ionic liquids phases onto copper nanoparticles: A density functional theory study

This work contains data on the computational, structural, and electronic characterization of supported ionic liquids phases anchored to copper nanoparticles using Density Functional theory calculations. The data supplement the paper "Interaction of supported ionic liquids phases onto copper nanoparticles: A Density Functional Theory study" [1], based on the adsorption of ionic liquid onto a Cu nanoparticle is analyzed from a chemical and physical point of view. The chemical analysis is based on Atoms in Molecule theory (AIM) and allows us to differentiate the chemical binding nature between ionic liquid and copper nanoparticle. On the other hand, the energy decomposition analysis based on absolutely localized molecular orbital (ALMO-EDA) describes the physical contributions that govern the interaction between ionic liquid and the copper nanoparticles. Herein, detailed and extended information in the synthesis and computational characterization are presented.

Bridge Resonance Effects in Singlet Fission

A major benefit of intramolecular singlet fission (iSF) materials, in which through-bond interactions mediate triplet pair formation, is the ability to control the triplet formation dynamics through molecular engineering. One common design strategy is the use of molecular bridges to mediate interchromophore interactions, decreasing electronic coupling by increasing chromophore-chromophore separation. Here, we report how the judicious choice of aromatic bridges can enhance chromophore-chromophore electronic coupling. This molecular engineering strategy takes advantage of "bridge resonance", in which the frontier orbital energies are nearly degenerate with those of the covalently linked singlet fission chromophores, resulting in fast iSF even at large interchromophore separations. Using transient absorption spectroscopy, we investigate this bridge resonance effect in a series of pentacene and tetracene-bridged dimers, and we find that the rate of triplet formation is enhanced as the bridge orbitals approach resonance. This work highlights the important role of molecular connectivity in controlling the rate of iSF through chemical bonds and establishes critical design principles for future use of iSF materials in optoelectronic devices.

Interplay of Open-Shell Spin-Coupling and Jahn-Teller Distortion in Benzene Radical Cation Probed by X-ray Spectroscopy

We report a theoretical investigation and elucidation of the X-ray absorption spectra of neutral benzene and of the benzene cation. The generation of the cation by multiphoton ultraviolet (UV) ionization and the measurement of the carbon K-edge spectra of both species using a table-top high-harmonic generation source are described in the companion experimental paper [Epshtein, M.; et al. J. Phys. Chem. A http://dx.doi.org/10.1021/acs.jpca.0c08736]. We show that the 1s_C → π transition serves as a sensitive signature of the transient cation formation, as it occurs outside of the spectral window of the parent neutral species. Moreover, the presence of the unpaired (spectator) electron in the π-subshell of the cation and the high symmetry of the system result in significant differences relative to neutral benzene in the spectral features associated with the 1s_C → π* transitions. High-level calculations using equation-of-motion coupled-cluster theory provide the interpretation of the experimental spectra and insight into the electronic structure of benzene and its cation. The prominent split structure of the 1s_C → π* band of the cation is attributed to the interplay between the coupling of the core → π* excitation with the unpaired electron in the π-subshell and the Jahn-Teller distortion. The calculations attribute most of the splitting (∼1-1.2 eV) to the spin coupling, which is visible already at the Franck-Condon structure, and we estimate the additional splitting due to structural relaxation to be around ∼0.1-0.2 eV. These results suggest that X-ray absorption with increased resolution might be able to disentangle electronic and structural aspects of the Jahn-Teller effect in the benzene cation.

Equation-of-Motion Coupled-Cluster Theory to Model L-Edge X-ray Absorption and Photoelectron Spectra

We present an extension of the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) theory for computing X-ray L-edge spectra, both in the absorption (XAS) and in the photoelectron (XPS) regimes. The approach is based on the perturbative evaluation of spin-orbit couplings using the Breit-Pauli Hamiltonian and nonrelativistic wave functions described by the fc-CVS-EOM-CCSD ansatz (EOM-CCSD within the frozen-core core-valence separated (fc-CVS) scheme). The formalism is based on spinless one-particle density matrices. The approach is illustrated by modeling XAS and XPS of several model systems ranging from Ar to small molecules containing sulfur and silicon.

A Series of Discrete Dichloride Dihydrates: Characterisation and Symmetry Effects

A series of three discrete dichloride dihydrates [Cl₂ (H₂ O)₂ ]^2- have been isolated with different triaminocyclopropenium (TAC) cations and with different crystallographic symmetries. The cluster exhibits D_2h symmetry with the tris(dimethylamino)cyclopropenium cation [C₃ (NMe₂ )₃ ]⁺ , C_2h symmetry with the fluorinated cation [C₃ (N(CH₂ CF₃ )₂ )(NBu₂ )₂ ]⁺ (containing two 2,2,2-trifluoroethyl substituents) and C_2v symmetry with the more fluorinated [C₃ (N(CH₂ CF₃ )₂ )₂ (NBu₂ )]⁺ cation. The effect of symmetry on the infrared spectra of the dichloride ion-pair clusters, as well as deuterated analogues, has been investigated. The D_2h - and C_2h -symmetric clusters each exhibit two stretching bands in the infrared at 3427 and 3368 cm^-1 for D_2h symmetry and 3444 and 3392 cm^-1 for C_2h symmetry, whereas the C_2v -symmetric cluster exhibits three bands at 3475, 3426 and 3373 cm^-1 . Computational studies were carried out on a [Cl₂ (H₂ O)₂ ]^2- cluster with C_2v symmetry to aid the infrared band assignments.