The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

Discriminating Clockwise and Counterclockwise Photoisomerization Paths in Achiral Photoswitches by Excited-State Electronic Circular Dichroism

Despite the numerous investigations of photoisomerization reactions from both the computational and experimental points of view, even in complex environments, to date there is no direct demonstration of the direction of rotation of the retinal chromophore, initiating the vision process in several organisms, occurring upon light irradiation. In the literature, many proposals have been formulated to shed light on the details of this process, most of which are extracted from semiclassical simulations. Although high hopes are held in the development of time-resolved X-ray spectroscopy, I argue in this work that simpler but less known techniques can be used to unravel the details of this fascinating photochemical process. In fact, chiroptical spectroscopy would unambiguously prove the direction of the rotatory motion of the chromophore during the photoisomerization process by probing excited state chirality, a piece of information that, so far, has been exclusively extracted from atomistic simulations. I demonstrate this statement by computing the expected chiroptical response along photoisomerization pathways for several models of the retinal chromophores that are found in nature bound to rhodopsins, including nuclear ensemble spectra from semiclassical dynamics simulations, that can be compared with time-resolved experiments.

Designing Quantum Spaces of Higher Dimensionality from a Tetranuclear Erbium-Based Single-Molecule Magnet

The spin relaxation of an Er3+ tetranuclear single-molecule magnet, [Er(hdcCOT)I]4, (hdcCOT = hexahydrodicyclopentacyclooctatetraenide dianion), is modeled as a near-tetrahedral arrangement of Ising-type spins. Combining evidence from single-crystal X-ray diffraction, magnetometry, and computational techniques, the slow spin relaxation is interpreted as a consequence of symmetry restrictions imposed on quantum tunneling within the cluster core. The union of spin and spatial symmetries describe a ground state spin-spin coupled manifold wherein 16 eigenvectors generate the 3D quantum spin-space described by the vertices of a rhombic dodecahedron. Analysis of the experimental findings in this context reveals a correlation between the magnetic transitions and edges connecting cubic and octahedral subsets of the eigenspace convex hull. Additionally, the model is shown to map to a theoretically proposed quantum Cayley network, indicating an underexplored synergy between mathematical descriptions of molecular spin interactions and quantum computing configuration spaces.

Permutation symmetry in spin-adapted many-body wave functions

In the domain of exchange-coupled polynuclear transition-metal (PNTM) clusters, local emergent symmetries exist which can be exploited to greatly increase the sparsity of the configuration interaction (CI) eigensolutions of such systems. Sparsity of the CI secular problem is revealed by exploring the site permutation space within spin-adapted many-body bases, and highly compressed wave functions may arise by finding optimal site orderings. However, the factorial cost of searching through the permutation space remains a bottleneck for clusters with a large number of metal centers. In this work, we explore ways to reduce the factorial scaling, by combining permutation and point group symmetry arguments, and using commutation relations between cumulative partial spin and the Hamiltonian operators, . Certain site orderings lead to commuting operators, from which more sparse wave functions arise. Two graphical strategies will be discussed, one to rapidly evaluate the commutators of interest, and one in the form of a tree search algorithm to predict how many and which distinct site permutations are to be analyzed, eliminating redundancies in the permutation space. Particularly interesting is the case of the singlet spin states for which an additional reversal symmetry can be utilized to further reduce the number of distinct site permutations.

Identification of a Weak Interaction between the Spin-Only 4f-Electronic System of Trivalent Gd Ion and a Photoexcited Cyclic π Electronic System

An electronic interaction between the spin-only 4f-electronic system and the photoexcited π cyclic system has been identified in the anionic bis(phthalocyaninato)gadolinium(III) complex [Pc2Gd]- by using variable-temperature variable-magnetic-field magnetic circular dichroism (VTVH MCD). Two positive MCD A-term patterns, corresponding to the QH and QL absorption bands, were observed to increase in intensity as the temperature decreased, indicating a ferromagnetic-type interaction between spin angular momentum S of the 4f system and orbital angular momentum L of the photoexcited π system. A theoretical model incorporating the S-L interaction constant CSL was fitted to the VTVH MCD spectra, confirming the presence of a weak interaction, whose intensity was determined as CSL = 0.27 cm-1 for QL and CSL = 0.02 cm-1 for QH. Ab initio calculations using the RASSCF and RASSI methods revealed that the inclusion of the excited state calculations with the spin multiplicity (2S + 1) of 6 was crucial for assigning ΔSL and CSL values. This underscores the essential role of the 6PJ term in accurately modeling the S-L interaction in [Pc2Gd]-.

Forte: A suite of advanced multireference quantum chemistry methods

Forte is an open-source library specialized in multireference electronic structure theories for molecular systems and the rapid prototyping of new methods. This paper gives an overview of the capabilities of Forte, its software architecture, and examples of applications enabled by the methods it implements.

The Use of Effective Core Potentials with Multiconfiguration Pair-Density Functional Theory

The reliable and accurate prediction of chemical properties is a key goal in quantum chemistry. Transition-metal-containing complexes can often pose difficulties to quantum mechanical methods for multiple reasons, including many electron configurations contributing to the overall electronic description of the system and the large number of electrons significantly increasing the amount of computational resources required. Often, multiconfigurational electronic structure methods are employed for such systems, and the cost of these calculations can be reduced by the use of an effective core potential (ECP). In this work, we explore both theoretical considerations and performances of ECPs applied in the context of multiconfiguration pair-density functional theory (MC-PDFT). A mixed-basis set approach is used, using ECP basis sets for transition metals and all-electron basis sets for nonmetal atoms. We illustrate the effects that an ECP has on the key parameters used in the computation of MC-PDFT energies, and we explore how ECPs affect the prediction of physical observables for chemical systems. The dissociation curve for a metal dimer was explored, and ionization energies for transition metal-containing diatomic systems were computed and compared to experimental values. In general, we find that ECP approaches employed with MC-PDFT are able to predict ionization energies with improved accuracy compared to traditional Kohn-Sham density functional theory approaches.

The Complex (Organic) Puzzle of the Formation of Hydrogen Cyanide and Isocyanide on Interstellar Ice Analogues

Aiming to constrain the surface formation of HCN and HNC in the dense interstellar medium on ice-covered dust grains, we investigate the interaction of CN radicals with H2O and CO ices and their subsequent reactivity with H and H2. CN radicals can physisorb on both ices. However, on H2O ice, a hemibond formation is the most common binding mode, while on CO ice, the CN-CO van der Waals complex can form NCCO with a small energy barrier. We show low barrier or barrierless pathways to the formation of HCN and HNC for the reaction H + CN on both ices. Reactivity with H2 involves activation energy barriers to form HCN, which may be overcome by quantum tunneling, while HNC formation is unlikely. The formation of HCN and HNC competes with the formation of NH2CHO on H2O and HCOCN on CO.

Triplet-Singlet Emission of d-Block Metal Complexes Characterized by Spin-Orbit Natural Transition Orbitals

Spin-orbit natural transition orbital (SO-NTO) methodology, recently developed in our group for complete and restricted active space (CAS/RAS) wavefunction calculations, is applied to analyze triplet-to-singlet emission in transition metal complexes. The lowest-energy (longest-wavelength) spin-forbidden transitionT 1 → S 0 ${{{\rm T}}_{1}\to {{\rm S}}_{0}}$ is studied for for [Ir(pbt)2(acac)] and [Re(CO)4(pbt)] and the complexes [W(CO)4(bpy)] and [Mo(CO)4(bpy)]. For the latter complexes, spin-forbidden transitions from higher spin-triplet levels are additionally analyzed. SO-NTOs are compared with spin-free NTOs for the transitions under consideration. The major assignment of a spin-forbidden transition is obtained from the spin-free NTO analysis, while the source of intensity of the electronic transition is revealed by the SO-NTOs.

Simulating ultrafast transient absorption spectra from first principles using a time-dependent configuration interaction probe

Transient absorption spectroscopy (TAS) is among the most common ultrafast photochemical experiments, but its interpretation remains challenging. In this work, we present an efficient and robust method for simulating TAS signals from first principles. Excited-state absorption and stimulated emission (SE) signals are computed using time-dependent complete active space configuration interaction (TD-CASCI) simulations, leveraging the robustness of time-domain simulation to minimize electronic structure failure. We demonstrate our approach by simulating the TAS signal of 1'-hydroxy-2'-acetonapthone (HAN) from ab initio multiple spawning nonadiabatic molecular dynamics simulations. Our results are compared to gas-phase TAS data recorded from both jet-cooled (T ∼ 40 K) and hot (∼403 K) molecules via cavity-enhanced TAS (CE-TAS). Decomposition of the computed spectrum allows us to assign a rise in the SE signal to excited-state proton transfer and the ultimate decay of the signal to relaxation through a twisted conical intersection. The total cost of computing the observable signal (∼1700 graphics processing unit hours for ∼4 ns of electron dynamics) was markedly less than that of performing the ab initio multiple spawning calculations used to compute the underlying nonadiabatic dynamics.

Near-infrared luminescent open-shell π-conjugated systems with a bright lowest-energy zwitterionic singlet excited state

Open-shell systems with extensive π-conjugation have fascinating properties due to their narrow bandgaps and spin interactions. In this work, we report neutral open-shell di- and polyradical conjugated materials exhibiting intriguing optical and magnetic properties. Our key design advance is the planarized geometry allowing for greater interaction between adjacent spins. This results in absorption and emission in the near infrared at 803 and 1050 nanometers, respectively, and we demonstrate a unique electronic structure where a bright zwitterionic excited state is the lowest-accessible electronic transition. Electron paramagnetic resonance spectroscopy and superconducting quantum interference device measurements reveal that our materials are open-shell singlets with different degrees of spin interactions, dynamics, and antiferromagnetic properties, which likely contributed to the formation of their emissive zwitterionic singlet excited state and near-infrared emission. In addition, our materials show reversible and stable electrochromic switching with more than 500 cycles, indicating their potential for optoelectronic and electrochemical energy storage applications.

Tris-Silanide f-Block Complexes: Insights into Paramagnetic Influence on NMR Chemical Shifts

The paramagnetism of f-block ions has been exploited in chiral shift reagents and magnetic resonance imaging, but these applications tend to focus on 1H NMR shifts as paramagnetic broadening makes less sensitive nuclei more difficult to study. Here we report a solution and solid-state (ss) 29Si NMR study of an isostructural series of locally D 3h -symmetric early f-block metal(III) tris-hypersilanide complexes, [M{Si(SiMe3)3}3(THF)2] (1-M; M = La, Ce, Pr, Nd, U); 1-M were also characterized by single crystal and powder X-ray diffraction, EPR, ATR-IR, and UV-vis-NIR spectroscopies, SQUID magnetometry, and elemental analysis. Only one SiMe3 signal was observed in the 29Si ssNMR spectra of 1-M, while two SiMe3 signals were seen in solution 29Si NMR spectra of 1-La and 1-Ce. This is attributed to dynamic averaging of the SiMe3 groups in 1-M in the solid state due to free rotation of the M-Si bonds and dissociation of THF from 1-M in solution to give the locally C 3v -symmetric complexes [M{Si(SiMe3)3}3(THF) n ] (n = 0 or 1), which show restricted rotation of M-Si bonds on the NMR time scale. Density functional theory and complete active space self-consistent field spin-orbit calculations were performed on 1-M and desolvated solution species to model paramagnetic NMR shifts. We find excellent agreement of experimental 29Si NMR data for diamagnetic 1-La, suggesting n = 1 in solution and reasonable agreement of calculated paramagnetic shifts of SiMe3 groups for 1-M (M = Pr and Nd); the NMR shifts for metal-bound 29Si nuclei could only be reproduced for diamagnetic 1-La, showing the current limitations of pNMR calculations for larger nuclei.

Non-adiabatic dynamics of photoexcited cyclobutanone: Predicting structural measurements from trajectory surface hopping with XMS-CASPT2 simulations

Over the years, theoretical calculations and scalable computer simulations have complemented ultrafast experiments, as they offer the advantage of overcoming experimental restrictions and having access to the whole dynamics. This synergy between theory and experiment promises to yield a deeper understanding of photochemical processes, offering valuable insights into the behavior of complex systems at the molecular level. However, the ability of theoretical models to predict ultrafast experimental outcomes has remained largely unexplored. In this work, we aim to predict the electron diffraction signals of an upcoming ultrafast photochemical experiment using high-level electronic structure calculations and non-adiabatic dynamics simulations. In particular, we perform trajectory surface hopping with extended multi-state complete active space with second order perturbation simulations for understanding the photodissociation of cyclobutanone (CB) upon excitation at 200 nm. Spin-orbit couplings are considered for investigating the role of triplet states. Our simulations capture the bond cleavage after ultrafast relaxation from the 3s Rydberg state, leading to the formation of the previously observed primary photoproducts: CO + cyclopropane/propene (C3 products), ketene, and ethene (C2 products). The ratio of the C3:C2 products is found to be about 1:1. Within 700 fs, the majority of trajectories transition to their electronic ground state, with a small fraction conserving the initial cyclobutanone ring structure. We found a minimal influence of triplet states during the early stages of the dynamics, with their significance increasing at later times. We simulate MeV-ultrafast electron diffraction (UED) patterns from our trajectory results, linking the observed features with specific photoproducts and the underlying structural dynamics. Our analysis reveals highly intense features in the UED signals corresponding to the photochemical processes of CB. These features offer valuable insights into the experimental monitoring of ring opening dynamics and the formation of C3 and C2 photoproducts.

Distinguishing homolytic vs heterolytic bond dissociation of phenylsulfonium cations with localized active space methods

Modeling chemical reactions with quantum chemical methods is challenging when the electronic structure varies significantly throughout the reaction and when electronic excited states are involved. Multireference methods, such as complete active space self-consistent field (CASSCF), can handle these multiconfigurational situations. However, even if the size of the needed active space is affordable, in many cases, the active space does not change consistently from reactant to product, causing discontinuities in the potential energy surface. The localized active space SCF (LASSCF) is a cheaper alternative to CASSCF for strongly correlated systems with weakly correlated fragments. The method is used for the first time to study a chemical reaction, namely the bond dissociation of a mono-, di-, and triphenylsulfonium cation. LASSCF calculations generate smooth potential energy scans more easily than the corresponding, more computationally expensive CASSCF calculations while predicting similar bond dissociation energies. Our calculations suggest a homolytic bond cleavage for di- and triphenylsulfonium and a heterolytic pathway for monophenylsulfonium.

XMS-CASPT2//XMS-CASPT2 and XMS-CASPT2//CASSCF at comparison: The impact of dynamic correlation in the excited state optimization of nitronaphthalene

Here, analytical extended multi-state complete active space second-order perturbation method (XMS-CASPT2) gradients are used to rationalize the decreasing triplet quantum yield trend in 2-nitronaphthalene, 1-nitronaphthalene, and 2-methyl-1-nitronaphthalene, a series of nitro-substituted aromatic compounds. Comparison with the XMS-CASPT2//CASSCF (where CASSCF stands for complete active space self-consistent field method) results highlights the importance of dynamic correlation in geometry optimization and challenges the validity of an XMS-CASPT2//CASSCF approach: XMS-CASPT2 S1 minima leads to planar structures, while CASSCF optimizations trigger a pyramidalization of the nitro group. The XMS-CASPT2 results correlate the reported decreasing triplet quantum yield trend in these species to a decrease in S1 to T2 population transfer and an increase in S1-S0 decay, while no such correlation is observed when using XMS-CASPT2//CASSCF data.

Intramolecular bridging strategies to suppress two-phonon Raman spin relaxation in dysprosocenium single-molecule magnets

Dy(III) bis-cyclopentadienyl (Cp) sandwich compounds exhibit extremely strong single-ion magnetic anisotropy which imbues them with magnetic memory effects such as magnetic hysteresis, and has put them at the forefront of high-performance single-molecule magnets (SMMs). Owing to the great success of design principles focused on maximising the anisotropy barrier, ever higher Ueff values have been reported leading to significant slow down of single-phonon Orbach spin relaxation. However, anisotropy-based SMM design has largely ignored two-phonon Raman spin relaxation, which is still limiting the temperatures at which a memory effect can be observed. In this work, we study the suppression of Raman relaxation through covalent bridging of the Cp ligands by alkyl chains, testing the hypothesis that increasing the rigidity of the ligand framework results in a blue shift of low frequency vibrations in the first coordination sphere of the Dy(III) ion. This reshaping of the vibrational low-energy density of states (DOS) results in lower occupation of pseudo-acoustic phonons available to drive Raman relaxation at low temperatures. We simulate Orbach and Raman spin relaxation in a series of zero-, mono-, di- and tri-bridged [Dy(Cpttt)2]+ analogues fully ab initio, using a quantum mechanics (QM)/molecular mechanics (MM) condensed phase embedding protocol in a periodic solvent matrix as a generic and experimentally testable environment model that can include (pseudo-)acoustic phononic degrees of freedom. We show that this approach can simulate magnetic relaxation dynamics in the condensed phase for the existing non-bridged [Dy(Cpttt)2]+ compound with quantitative experimental accuracy. Subsequently, we find a significant slowing down of Raman relaxation can be achieved for the singly-bridged SMM, while the introduction of further bridges leads to faster relaxation. A key result being that we find the two-phonon Raman rates correlate with the purity of the first-excited Kramers doublet in terms of its mJ = ±13/2 content. Even though the bridging design principle is successful at progressively reshaping the low-energy DOS, the introduction of linker atoms in the equatorial plane successively degrades magnetic anisotropy, suggesting the importance of refined design of the linker chemistry. The accuracy of our results emphasises the value of a generic periodic solvent embedding model, such that it permits the modelling of molecular spin dynamics in the condensed phase without knowledge of a crystal structure. This allows the study of hypothetical molecules or aggregates under real-world conditions, which we expect to have utility beyond the field of molecular magnetism.

SCINE-Software for chemical interaction networks

The software for chemical interaction networks (SCINE) project aims at pushing the frontier of quantum chemical calculations on molecular structures to a new level. While calculations on individual structures as well as on simple relations between them have become routine in chemistry, new developments have pushed the frontier in the field to high-throughput calculations. Chemical relations may be created by a search for specific molecular properties in a molecular design attempt, or they can be defined by a set of elementary reaction steps that form a chemical reaction network. The software modules of SCINE have been designed to facilitate such studies. The features of the modules are (i) general applicability of the applied methodologies ranging from electronic structure (no restriction to specific elements of the periodic table) to microkinetic modeling (with little restrictions on molecularity), full modularity so that SCINE modules can also be applied as stand-alone programs or be exchanged for external software packages that fulfill a similar purpose (to increase options for computational campaigns and to provide alternatives in case of tasks that are hard or impossible to accomplish with certain programs), (ii) high stability and autonomous operations so that control and steering by an operator are as easy as possible, and (iii) easy embedding into complex heterogeneous environments for molecular structures taken individually or in the context of a reaction network. A graphical user interface unites all modules and ensures interoperability. All components of the software have been made available as open source and free of charge.

Field-Induced Slow Magnetic Relaxation in Mononuclear Cobalt(II) Complexes Decorated by Macrocyclic Pentaaza Ligands

Two cobalt(II) complexes [CoL1](OTf)2 (1, L1 = 6,6''-di(anilino)-4'-phenyl-2,2':6',2''-terpyridine) and [CoL2](OTf)2·MeOH (2, L2 = 6,6''-di(N,N-dimethylamino)-4'-phenyl-2,2':6',2''-terpyridine) were synthesized and characterized. Crystal structure analyses showed that the spin carries were coordinated by five N atoms from the neutral pentaaza ligands, forming distorted trigonal bipyramidal coordination environments. Ab initio calculations revealed large easy-axial anisotropy in complexes 1 and 2. Magnetic measurements suggest that complexes 1 and 2 are field-induced single-molecule magnets, whose relaxations are mainly predominated by Raman and direct processes.

WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation

We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.

Multiconfigurational actinide nitrides assisted by double Möbius aromaticity

Understanding the bonding nature between actinides and main-group elements remains a key challenge in actinide chemistry due to the involvement of f orbitals. Herein, we propose a unique "aromaticity-assisted multiconfiguration" (AAM) model to elucidate the bonding nature in actinide nitrides (An2N2, An = Ac, Th, Pa, U). Each planar four-membered An2N2 with equivalent An-N bonds possesses four delocalized π electrons and four delocalized σ electrons, forming a new family of double Möbius aromaticity that contributes to the molecular stability. The unprecedented aromaticity further supports actinide nitrides to exhibit multiconfigurational characters, where the unpaired electrons (2, 4 or 6 in naked Th2N2, Pa2N2 or U2N2, respectively) either are spin-free and localized on metal centres or form metal-ligand bonds. High-level multiconfigurational computations confirm an open-shell singlet ground state for actinide nitrides, with small energy gaps to high spin states. This is consistent with the antiferromagnetic nature observed experimentally in uranium nitrides. The novel AAM bonding model can be authenticated in both experimentally identified compounds containing a U2N2 motif and other theoretically modelled An2N2 clusters and is thus expected to be a general chemical bonding pattern between actinides and main-group elements.

Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone

Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α-CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.

Toward Real Chemical Accuracy on Current Quantum Hardware Through the Transcorrelated Method

Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and connectivity) hamper the implementation of most quantum algorithms and call for more noise-resilient solutions. We propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach to target these major roadblocks directly. This method transfers, without any approximation, correlations from the wave function directly into the Hamiltonian, thus reducing the resources needed to achieve accurate results with noisy quantum devices. We show that the TC approach allows for shallower circuits and improves the convergence toward the complete basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, thus, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using two and four qubits, respectively. To demonstrate our approach's current and near-term potential, we perform hardware experiments, where our results confirm that the TC method paves the way toward accurate quantum chemistry calculations already on today's quantum hardware.

Substituent Positioning Effects on the Magnetic Properties of Sandwich-Type Erbium(III) Complexes with Bis(trimethylsilyl)-Substituted Cyclooctatetraenyl Ligands

Lanthanide complexes with judiciously designed ligands have been extensively studied for their potential applications as single-molecule magnets. With the influence of ligands on their magnetic properties generally established, recent research has unearthed certain effects inherent to site differentiation due to the different types and varying numbers of substituents on the same ligand platform. Using two new sandwich-type Er(III) complexes with cyclooctatetraenyl (COT) ligands featuring two differently positioned trimethylsilyl (TMS) substituents, namely, [Li(DME)Er(COT1,5-TMS2)2]n (Er1) and [Na(DME)3][Er(COT1,3-TMS2)2] (Er2) [COT1,3-TMS2 and COT1,5-TMS2 donate 1,3- and 1,5-bis(trimethylsilyl)-substituted cyclooctatetraenyl ligands, respectively; DME = 1,2-dimethoxyethane], and with reference to previously reported [Li(DME)3][Er(COT1,4-TMS2)2] (A) and [K(DME)2][Er(COT1,4-TMS2)2] (B), any possible substituent position effects have been explored for the first time. The rearrangement of the TMS substituents from the starting COT1,4-TMS2 to COT1,3-TMS2 and COT1,5-TMS2, by way of formal migration of the TMS group, was thermally induced in the case of Er1, while for the formation of Er2, the use of Na+ in the placement of its Li+ and K+ congeners is essential. Both Er1 and Er2 display single-molecule magnetic behaviors with energy barriers of 170(3) and 172(6) K, respectively. Magnetic hysteresis loops, butterfly-shaped for Er1 and wide open for Er2, were observed up to 12 K for Er1 and 13 K for Er2. Studies of magnetic dynamics reveal the different pathways for relaxation of magnetization below 10 K, mainly by the Raman process for Er1 and by quantum tunneling of magnetization for Er2, leading to the order of magnitude difference in magnetic relaxation times and sharply different magnetic hysteresis loops.

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear-electronic coupling in addition to accurately predicting Franck-Condon vertical excitations. In this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

Modulated spin dynamics of [Co2] coordination helicates via differential strand composition

Coordination supramolecular chemistry provides a versatile entry into materials with functionalities of technological relevance at the nanoscale. Here, we describe how two different bis-pyrazolylpyridine ligands (L1 and L2) assemble with Co(II) ions into dinuclear triple-stranded helicates, in turn, encapsulating different anionic guests. These constructs are described as (Cl@[Co2(L1)3])3+, (SiF6@[Co2(L1)(L2)3])2+ and (ClO4@[Co2(L2)3])3+, as established by single-crystal X-ray diffraction. Extensive magnetic and calorimetric measurements, numerical treatments and theoretical calculations reveal that the individual Co(II) centers of these supramolecular entities exhibit field-induced slow relaxation of magnetization, dominated by direct and Raman mechanisms. While the small variations in the spin dynamics are not easily correlated with the evident structural differences among the three species, the specific heat measurements suggest two vibronic pathways of magnetic relaxation: one that would be associated with the host lattice and another linked with the guest.

Interoperable workflows by exchanging grid-based data between quantum-chemical program packages

Quantum-chemical subsystem and embedding methods require complex workflows that may involve multiple quantum-chemical program packages. Moreover, such workflows require the exchange of voluminous data that go beyond simple quantities, such as molecular structures and energies. Here, we describe our approach for addressing this interoperability challenge by exchanging electron densities and embedding potentials as grid-based data. We describe the approach that we have implemented to this end in a dedicated code, PyEmbed, currently part of a Python scripting framework. We discuss how it has facilitated the development of quantum-chemical subsystem and embedding methods and highlight several applications that have been enabled by PyEmbed, including wave-function theory (WFT) in density-functional theory (DFT) embedding schemes mixing non-relativistic and relativistic electronic structure methods, real-time time-dependent DFT-in-DFT approaches, the density-based many-body expansion, and workflows including real-space data analysis and visualization. Our approach demonstrates, in particular, the merits of exchanging (complex) grid-based data and, in general, the potential of modular software development in quantum chemistry, which hinges upon libraries that facilitate interoperability.

Ultrafast photochemistry and electron-diffraction spectra in n → (3s) Rydberg excited cyclobutanone resolved at the multireference perturbative level

We study the ultrafast time evolution of cyclobutanone excited to the singlet n → Rydberg state through non-adiabatic surface-hopping simulationsperformed at extended multi-state complete active space second-order perturbation (XMS-CASPT2) level of theory. These dynamics predict relaxation to the ground-state with a timescale of 822 ± 45 fs with minimal involvement of the triplets. The major relaxation path to the ground-state involves a three-state degeneracy region and leads to a variety of fragmented photoproducts. We simulate the resulting time-resolved electron-diffraction spectra, which track the relaxation of the excited state and the formation of various photoproducts in the ground state.

Approaching the Complete Basis Set Limit for Spin-State Energetics of Mononuclear First-Row Transition Metal Complexes

Convergence to the complete basis set (CBS) limit is analyzed for the problem of spin-state energetics in mononuclear first-row transition metal (TM) complexes by taking under scrutiny a benchmark set of 18 energy differences between spin states for 13 chemically diverse TM complexes. The performance of conventional CCSD(T) and explicitly correlated CCSD(T)-F12a/b calculations in approaching the CCSD(T)/CBS limits is systematically studied. An economic computational protocol is developed based on the CCSD-F12a approximation and (here proposed) modified scaling of the perturbative triples term (T#). This computational protocol recovers the relative spin-state energetics of the benchmark set in excellent agreement with the reference CCSD(T)/CBS limits (mean absolute deviation of 0.4, mean signed deviation of 0.2, and maximum deviation of 0.8 kcal/mol) and enables performing canonical CCSD(T) calculations for mononuclear TM complexes sized up to ca. 50 atoms, which is illustrated by application to heme-related metalloporphyrins. Furthermore, a good transferability of the basis set incompleteness error (BSIE) is demonstrated for spin-state energetics computed using CCSD(T) and other wave function methods (MP2, CASPT2, CASPT2/CC, NEVPT2, and MRCI + Q), which justifies efficient focal-point approximations and simplifies the construction of multimethod benchmark studies.

Effect of diamagnetic Zn(II) ions on the SMM properties of a series of trinuclear ZnDy2 and tetranuclear Zn2Dy2 (LnIII = Dy, Tb, Gd) complexes: combined experimental and theoretical studies

To study the effect of diamagnetic ions on magnetic interactions, utilizing a compartmental ligand (Z)-2-(hydroxymethyl)-4-methyl-6-((quinolin-8-ylimino)methyl)phenol (LH2), two different series of ZnII-LnIII complexes, namely the trinuclear series of [DyZn2(L)2(μ2-OAc)2(CH3OH)2]·NO3·MeOH (1), [TbZn2(L)2(μ2-OAc)2(CH3OH)2]·NO3·5MeOH·H2O (2), and [GdZn2(L)2(μ2-OAc)2(CH3OH)2]·NO3·MeOH·CHCl3 (3) and the tetranuclear series of [Dy2Zn2(LH)4(NO3)4(μ2OAc)]·NO3·MeOH·H2O (4), [Tb2Zn2(LH)4(NO3)4(μ2-OAc)]·NO3·MeOH·2H2O (5), and [Gd2Zn2(LH)4(NO3)4(μ2-OAc)]·NO3·MeOH·2H2O (6), were synthesized. Trinuclear ZnII-LnIII complexes 1-3 consist of one LnIII ion sandwiched between two peripheral ZnII ions forming a bent type ZnII-DyIII-ZnII array with an angle of 110.64°. Tetranuclear ZnII-LnIII complexes 4-6 are basically a combination of two dinuclear moieties of [LnZn(LH)2(NO3)2]+ connected by one bidentate bridging acetate ion in μ2-OAc coordination mode. The detailed magnetic analysis reveals that complexes 1 and 4 are single molecule magnets having energy barriers of 34.98 K and 46.71 K with relaxation times (τ0) of 5.05 × 10-4 s and 5.24 × 10-4 s, respectively. Ab initio calculations were employed to analyze the magnetic anisotropy and magnetic exchange interaction between the ZnII and LnIII centers with the aim of gaining better insights into the magnetic dynamics of complexes 1-6.

Unraveling the Ultrafast Photochemical Dynamics of Nitrobenzene in Aqueous Solution

Nitroaromatic compounds are major constituents of the brown carbon aerosol particles in the troposphere that absorb near-ultraviolet (UV) and visible solar radiation and have a profound effect on the Earth's climate. The primary sources of brown carbon include biomass burning, forest fires, and residential burning of biofuels, and an important secondary source is photochemistry in aqueous cloud and fog droplets. Nitrobenzene is the smallest nitroaromatic molecule and a model for the photochemical behavior of larger nitroaromatic compounds. Despite the obvious importance of its droplet photochemistry to the atmospheric environment, there have not been any detailed studies of the ultrafast photochemical dynamics of nitrobenzene in aqueous solution. Here, we combine femtosecond transient absorption spectroscopy, time-resolved infrared spectroscopy, and quantum chemistry calculations to investigate the primary steps following the near-UV (λ ≥ 340 nm) photoexcitation of aqueous nitrobenzene. To understand the role of the surrounding water molecules in the photochemical dynamics of nitrobenzene, we compare the results of these investigations with analogous measurements in solutions of methanol, acetonitrile, and cyclohexane. We find that vibrational energy transfer to the aqueous environment quenches internal excitation, and therefore, unlike the gas phase, we do not observe any evidence for formation of photoproducts on timescales up to 500 ns. We also find that hydrogen bonding between nitrobenzene and surrounding water molecules slows the S1/S0 internal conversion process.

Galván I, Lindh R. A Story of Three Levels of Sophistication in SCF/KS-DFT Orbital Optimization Procedures

In this work, three versions of self-consistent field/Kohn-Sham density functional theory (SCF/KS-DFT) orbital optimization are described and benchmarked. The methods are a modified version of the geometry version of the direct inversion in the iterative subspace approach (which we call r-GDIIS), the modified restricted step rational function optimization method (RS-RFO), and the novel subspace gradient-enhanced Kriging method combined with restricted variance optimization (S-GEK/RVO). The modifications introduced are aimed at improving the robustness and computational scaling of the procedures. In particular, the subspace approach in S-GEK/RVO allows the application to SCF/KS-DFT optimization of a machine learning technique that has proven to be successful in geometry optimizations. The performance of the three methods is benchmarked for a large number of small- to medium-sized organic molecules, at equilibrium structures and close to a transition state, and a second set of molecules containing closed- and open-shell transition metals. The results indicate the importance of the resetting technique in boosting the performance of the r-GDIIS procedure. Moreover, it is demonstrated that already at the inception of the subspace version of GEK to optimize SCF wave functions, it displays superior and robust convergence properties as compared to those of the standard state-of-the-art SCF/KS-DFT optimization methods.

Variational Pair-Density Functional Theory: Dealing with Strong Correlation at the Protein Scale

Multiconfigurational pair-density functional theory (MC-PDFT) offers a promising solution to the challenges faced by traditional density functional theory (DFT) in addressing molecular systems containing transition metals, open-shells, or strong correlations in general. By utilizing both the density and on-top pair-density, MC-PDFT can make use of a more flexible multiconfigurational wave function to capture the necessary static correlation, while the pair-density functional also includes the effect of dynamic correlation. So far, MC-PDFT has been used after a multiconfigurational self-consistent field (MCSCF) step, using the orbitals and configuration interaction coefficients from the converged MCSCF wave function to compute PDFT energies and properties. Here, instead, we propose to perform a direct optimization of the wave function using the pair-density functionals, resulting in a variational formulation of MC-PDFT. We derive the expressions for the wave function gradient and illustrate their similarity to standard MCSCF equations. Furthermore, we illustrate the accuracy on a set of singlet-triplet gaps as well as dissociation curves. Our findings highlight one of MC-PDFT's standout features: a reduced dependency on the active space size compared to conventional multiconfigurational wave function methodologies. Additionally, we show that the computational cost of MC-PDFT is potentially lower than MCSCF and often on-par with standard Kohn-Sham DFT, which is demonstrated by performing a MC-PDFT calculation of the entire ferredoxin protein with 1447 atoms and nearly 12 000 basis functions.

Comparison of Variational and Perturbative Spin-Orbit Coupling within Two-Component CASSCF

The modeling of spin-orbit coupling (SOC) remains a challenge in computational chemistry due to the high computational cost. With the rising popularity of spin-driven processes and f-block metals in chemistry and materials science, it is incumbent on the community to develop accurate multiconfigurational SOC methods that scale to large systems and understand the limits of different treatments of SOC. Herein, we introduce an implementation of perturbative SOC in scalar-relativistic two-component CASSCF (srX2C-CASSCF-SO). Perspectives on the limitations and accuracy of srX2C-CASSCF-SO are presented via benchmark calculations.

Renormalized-Residue-Based Multireference Configuration Interaction Method for Strongly Correlated Systems

The implementation of multireference configuration interaction (MRCI) methods in quantum systems with large active spaces is hindered by the expansion of configuration bases or the intricate handling of reduced density matrices (RDMs). In this work, we present a spin-adapted renormalized-residue-based MRCI (RR-MRCI) approach that leverages renormalized residues to effectively capture the entanglement between active and inactive orbitals. This approach is reinforced by a novel efficient algorithm, which also facilitates an efficient deployment of spin-adapted matrix product state MRCI (MPS-MRCI). The RR-MRCI framework possesses several advantages: (1) It considers the orbital entanglement and utilizes highly compressed MPS structure, improving computational accuracy and efficiency compared with internally contracted (ic) MRCI. (2) Utilizing small-sized buffer environments of a few external orbitals as probes based on quantum information theory, it enhances computational efficiency over MPS-MRCI and offers potential application to large molecular systems. (3) The RR framework can be implemented in conjunction with ic-MRCI, eliminating the need for high-rank RDMs, by using distinct renormalized residues. We evaluated this method across nine diverse molecular systems, including Cu2O22+ with an active space of (24e,24o) and two complexes of lanthanide and actinide with active space (38e,36o), demonstrating the method's versatility and efficacy.

Minimum-Energy Conical Intersections by Compressed Multistate Pair-Density Functional Theory

Compressed multistate pair-density functional theory (CMS-PDFT) is a multistate version of multiconfiguration pair-density functional theory that can capture the correct topology of coupled potential energy surfaces (PESs) around conical intersections. In this work, we develop interstate coupling vectors (ISCs) for CMS-PDFT in the OpenMolcas and PySCF/mrh electronic structure packages. Yet, the main focus of this work is using ISCs to calculate minimum-energy conical intersections (MECIs) by CMS-PDFT. This is performed using the projected constrained optimization method in OpenMolcas, which uses ISCs to restrain the iterations to the conical intersection seam. We optimize the S1/S0 MECIs for ethylene, butadiene, and benzene and show that CMS-PDFT gives smooth PESs in the vicinities of the MECIs. Furthermore, the CMS-PDFT MECIs are in good agreement with the MECI calculated by the more expensive XMS-CASPT2 method.

Reconciling TD-DFT and CASPT2 electronic structure methods for describing the photophysics of DNA

Time-dependent density functional theory (TD-DFT) and multiconfigurational second-order perturbation theory (CASPT2) are two of the most widely used methods to investigate photoinduced dynamics in DNA-based systems. These methods sometimes give diverse dynamics in physiological environments usually modeled by quantum mechanics/molecular mechanics (QM/MM) protocol. In this work, we demonstrate for the uridine test case that the underlying topology of the potential energy surfaces of electronic states involved in photoinduced relaxation is similar in both electronic structure methods. This is verified by analyzing surface-hopping dynamics performed at the QM/MM level on aqueous solvated uridine at TD-DFT and CASPT2 levels. By constraining the dynamics to remain onπ π * state we observe similar fluctuations in energy and relaxation lifetimes in surface-hopping dynamics in both TD-DFT and experimentally validated CASPT2 methods. This finding calls for a systematic comparison of the ES potential energy surfaces of DNA and RNA nucleosides at the single- and multi-reference levels of theory. The anomalous long excited state lifetime at the TD-DFT level is explained byn π * trapping due to the tendency of TD-DFT in QM/MM schemes with electrostatic embedding to underestimate the energy of theπ π * state leading to a wrongπ π * / n π * energetic order. A study of the FC energetics suggests that improving the description of the surrounding environment through polarizable embedding or by the expansion of QM layer with hydrogen-bonded waters helps restore the correct state order at TD-DFT level. Thus by combining TDDFT with an accurate modeling of the environment, TD-DFT is positioned as the standout protocol to model photoinduced dynamics in DNA-based aggregates and multimers.

A First Proposal on the Nitrobenzene Photorelease Mechanism of NO2 and Its Relation to NO Formation through a Roaming Mechanism

Despite the fact that NO2 is considered to be the main photoproduct of nitrobenzene photochemistry, no mechanism has ever been proposed to rationalize its formation. NO photorelease is instead a more studied process, probably due to its application in the drug delivery sector and the study of roaming mechanisms. In this contribution, a photoinduced mechanism accounting for the formation of NO2 in nitrobenzene is theorized based on CASPT2, CASSCF, and DFT electronic structure calculations and CASSCF classical dynamics. A triplet nπ* state is shown to evolve toward C-NO2 dissociation, being, in fact, the only low-lying excited state favoring such a deformation. Along the triplet dissociation path, the possibility to decay to the singlet ground state results in the frustration of the dissociation and in the recombination of the fragments, either back to the nitro or the nitrite isomer. The thermal decomposition of the latter to NO constitutes globally a roaming mechanism of NO formation.

Photoisomerization of (Cyanomethylene)cyclopropane (C5H5N) to 1-Cyano-2-methylenecyclopropane in an Argon Matrix

Broad-band ultraviolet photolysis (λ > 200 nm) of (cyanomethylene)cyclopropane (5) in an argon matrix at 20 K generates 1-cyano-2-methylenecyclopropane (7), a previously unknown compound. This product was initially identified by comparison of its infrared spectrum to that predicted by an anharmonic MP2/6-311+G(2d,p) calculation. This assignment was unambiguously confirmed by the synthesis of 1-cyano-2-methylenecyclopropane (7) and observation of its authentic infrared spectrum, which proved identical to that of the observed photoproduct. We investigated the singlet and triplet potential energy surfaces associated with this isomerization process using density functional theory and multireference calculations. The observed rearrangement of compound 5 to compound 7 is computed to be endothermic (3.3 kcal/mol). We were unable to observe the reverse reaction (7 → 5) under the photochemical conditions.

Dissociation and Isomerization Following Ionization of Ethylene: Insights from Nonadiabatic Dynamics Simulations

Photoionized and electronically excited ethylene C2H4+ can undergo H-loss, H2-loss, and ethylene-ethylidene isomerization, where the latter entails a hydrogen migration. Recent pioneering experiments with few-femtosecond extreme ultraviolet pulses and complementary theoretical studies have shed light on the photodynamics of this prototypical organic cation. However, no theoretical investigation based on dynamics simulations reported to date has described the mechanisms and time scales of dissociation and isomerization. Herein, we simulate the coupled electron-nuclear dynamics of ethylene following vertical ionization and electronic excitation to its four lowest-lying cationic states. The electronic structure is treated at the CASSCF level, with an active space large enough to describe bond breaking and formation. The simulations indicate that dissociation and isomerization take place mainly on the cationic ground state and allow the probing of previous hypotheses concerning the correlation between the photochemical outcome and the traversed conical intersections. The results, moreover, support the long-standing view that H2-loss may occur from the ethylidene form. However, the ethylene-ethylidene isomerization time predicted by the simulations is considerably longer than those previously inferred from indirect experimental measurements.

Dissociation Time, Quantum Yield, and Dynamic Reaction Pathways in the Thermolysis of trans-3,4-Dimethyl-1,2-dioxetane

The thermolysis of trans-3,4-dimethyl-1,2-dioxetane is studied by trajectory surface hopping. The significant difference between long and short dissociation times is rationalized by frustrated dissociations and the time spent in triplet states. If the C-C bond breaks through an excited state channel, then the trajectory passes over a ridge of the potential energy surface of that state. The calculated triplet quantum yields match the experimental results. The dissociation half-times and quantum yields follow the same ascending order as per the product states, justifying the conjecture that the longer dissociation time leads to a higher quantum yield, proposed in the context of the methylation effect. The populations of the molecular Coulomb Hamiltonian and diagonal states reach equilibrium, but the triplet populations with different Sz components fluctuate indefinitely. Certain initial velocities, leading the trajectories to given product states, can be identified as the most characteristic features for sorting trajectories according to their product states.

Theoretical Prediction of Core-Level Binding Energies: Analysis of Unexpected Errors

The analysis of the C(1s) and O(1s) core-level binding energies (CLBEs) of selected molecules computed by means of total energy Hartree-Fock (ΔSCF-HF) differences shows that in some cases, the calculated values for the C(1s) are larger than the experiment, which is unexpected. The origin of these unexpected errors of the Hartree-Fock ΔSCF BEs is shown to arise from static, nondynamical, electron correlation effects which are larger for the ion than for the neutral system. Once these static correlation effects are included by using complete active space self-consistent field (CASSCF) wave functions that include internal correlation terms, the resulting ΔSCF BEs are, as expected, smaller than measured values.

Which Electronic Structure Method to Choose in Trajectory Surface Hopping Dynamics Simulations? Azomethane as a Case Study

Nonadiabatic dynamics simulations have become a standard approach to explore photochemical reactions. Such simulations require underlying potential energy surfaces and couplings between them, calculated at a chosen level of theory, yet this aspect is rarely assessed. Here, in combination with the popular trajectory surface hopping dynamics method, we use a high-accuracy XMS-CASPT2 electronic structure level as a benchmark for assessing the performances of various post-Hartree-Fock methods (namely, CIS, ADC(2), CC2, and CASSCF) and exchange-correlation functionals (PBE, PBE0, and CAM-B3LYP) in a TD-DFT/TDA context, using the isomerization around a double bond as test case. Different relaxation pathways are identified, and the ability of the different methods to reproduce their relative importance and time scale is discussed. The results show that multireference electronic structure methods should be preferred, when studying nonadiabatic decay between excited and ground states. If not affordable, TD-DFT with TDA and hybrid functionals and ADC(2) are efficient alternatives but overestimate the nonradiative decay yield and thus may miss deexcitation pathways.

Toward a Stochastic Complete Active Space Second-Order Perturbation Theory

In this work, an internally contracted stochastic complete active space second-order perturbation theory, stochastic-CASPT2, is reported. The method relies on stochastically sampled reduced density matrices (RDMs) up to rank four and contractions thereof with the generalized Fock matrix. A new protocol for calculating higher-order RDMs in full configuration interaction quantum Monte Carlo (FCIQMC) has been designed based on (1) restricting sampling of the corresponding excitations to a deterministic subspace, (2) averaging the RDMs from independent dynamics and (3) projecting them onto the closest positive semi-definite matrix. Our protocol avoids previously encountered numerical conditioning problems in the orthogonalization of the perturber overlap matrix stemming from numerical noise. The chromium dimer CASSCF(12,12)/CASPT2 binding curve is computed as a proof of concept.

Auger-Meitner and X-ray Absorption Spectra of Ethylene Cation: Insight into Conical Intersection Dynamics

We present a theoretical investigation of the near-edge X-ray absorption fine structure and the Auger-Meitner decay spectra of ethylene and its cation. Herein, we demonstrate that our method, coupled with the nuclear ensemble approach, successfully reproduces the natural bandwidth structure of the experimental resonant Auger-Meitner decay spectra of ethylene, which is not very well reproduced within the Franck-Condon approximation. Furthermore, we analyze the Auger-Meitner decay spectra of the ethylene cation in light of minimum energy conical intersection structures involving the two lowest cationic states (D1 and D0), providing valuable insights into the ultrafast D1/D0 relaxation dynamics. Our results suggest that Auger-Meitner electron spectroscopy can help elucidate the mechanism behind the initial 20 fs of the relaxation dynamics.

Time-Resolved X-ray Absorption Spectroscopy: An MCTDH Quantum Dynamics Protocol

Expressions for linear and nonlinear spectroscopy simulation in the X-ray window in which the time evolution of a photoexcited molecular system is treated via quantum dynamics are derived. By leveraging on the peculiar properties of core-excited/ionized states, first- and third-order response functions are recast in the limit of time-scale separation between the extremely short core-state lifetime and the (comparably longer) electronic-state transfer and nuclear vibrational motion. This work is a natural extension of Segatta et al. (J. Chem. Theory Comput. 2023, 19, 2075-2091), in which some of the present authors coupled MCTDH quantum dynamics to spectroscopy simulation at different levels of sophistication. Full quantum dynamics and approximate expressions are compared by simulating X-ray transient absorption spectroscopy at the carbon K-edge in the pyrene molecule.

Photochemical Transformations of Diverse Biologically Active Resveratrol Analogs in Batch and Flow Reactors

Previous biological tests have shown that some resveratrol analogs exhibited significant antioxidative and cholinesterase inhibitory potential, as evidenced by lower IC50 values compared to the established standards, resveratrol and galantamine, respectively. Photochemical transformations were made in parallel on these compounds in the presence of porphyrin photocatalysts in batch and microreactor, showing the significant advantage of flow photochemistry concerning productivity, selectivity, and yields. In this research, the products of photocatalysis and direct irradiation (photolysis) of resveratrol analogs were compared to elucidate how the types and ratios of the products depend on the excitation energy, to reveal the effects of the substituent on the photoinduced reactions and to rationalize experimentally and computationally the nature and ratio of the obtained products. Thus, two main paths were computed in agreement with the experimental results: isomerization with the participation of triplet state intermediates to yield the experimentally detected cis-isomers and subsequent cyclization following a pathway not available for the trans-isomers. The investigation of five model compounds confirmed the advantages of the flow photoreactor in the photochemical reactions of heterocyclic resveratrol analogs.

Thermally Stable Terbium(II) and Dysprosium(II) Bis-amidinate Complexes

The thermostable four-coordinate divalent lanthanide (Ln) bis-amidinate complexes [Ln(Piso)2] (Ln = Tb, Dy; Piso = {(NDipp)2CtBu}, Dipp = C6H3iPr2-2,6) were prepared by the reduction of parent five-coordinate Ln(III) precursors [Ln(Piso)2I] (Ln = Tb, Dy) with KC8; halide abstraction of [Ln(Piso)2I] with [H(SiEt3)2][B(C6F5)] gave the respective Ln(III) complexes [Ln(Piso)2][B(C6F5)]. All complexes were characterized by X-ray diffraction, ICP-MS, elemental analysis, SQUID magnetometry, UV-vis-NIR, ATR-IR, NMR, and EPR spectroscopy and ab initio CASSCF-SO calculations. These data consistently show that [Ln(Piso)2] formally exhibit Ln(II) centers with 4fn5dz21 (Ln = Tb, n = 8; Dy, n = 9) valence electron configurations. We show that simple assignments of the f-d coupling to either L-S or J-s schemes are an oversimplification, especially in the presence of significant crystal field splitting. The coordination geometry of [Ln(Piso)2] is intermediate between square planar and tetrahedral. Projecting from the quaternary carbon atoms of the CN2 ligand backbones shows near-linear C···Ln···C arrangements. This results in strong axial ligand fields to give effective energy barriers to magnetic reversal of 1920(91) K for the Tb(II) analogue and 1964(48) K for Dy(II), the highest values observed for mononuclear Ln(II) single-molecule magnets, eclipsing 1738 K for [Tb(C5iPr5)2]. We tentatively attribute the fast zero-field magnetic relaxation for these complexes at low temperatures to transverse fields, resulting in considerable mixing of mJ states.

Exploiting Locality in Full Configuration Interaction Quantum Monte Carlo for Fast Excitation Generation

In this paper, we propose an improved excitation generation algorithm for the full configuration interaction quantum Monte Carlo method, which is particularly effective in systems described by localized orbitals. The method is an extension of the precomputed heat-bath strategy of Holmes et al., with more effective sampling of double excitations and a novel approach for nonuniform sampling of single excitations. We demonstrate the effectiveness of the algorithm for a chain of 30 hydrogen atoms with atom-localized orbitals, a stack of benzene molecules, and an Fe(II)-porphyrin model complex, whereby we show an overall efficiency gain by a factor of two to four, as measured by variance reduction per wall-clock time.

Comparison of Computational Strategies for the Calculation of the Electronic Coupling in Intermolecular Energy and Electron Transport Processes

Electronic couplings in intermolecular electron and energy transfer processes calculated by six different existing computational techniques are compared to nonorthogonal configuration interaction for fragments (NOCI-F) results. The paper addresses the calculation of the electronic coupling in diketopyrrolopyrol, tetracene, 5,5'-difluoroindigo, and benzene-Cl for hole and electron transport, as well as the local exciton and singlet fission coupling. NOCI-F provides a rigorous computational scheme to calculate these couplings, but its computational cost is rather elevated. The here-considered ab initio Frenkel-Davydov (AIFD), Dimer projection (DIPRO), transition dipole moment coupling, Michl-Smith, effective Hamiltonian, and Mulliken-Hush approaches are computationally less demanding, and the comparison with the NOCI-F results shows that the NOCI-F results in the couplings for hole and electron transport are rather accurately predicted by the more approximate schemes but that the NOCI-F exciton transfer and singlet fission couplings are more difficult to reproduce.

Comparative Study of Uracil Excited-State Photophysics in Water and Acetonitrile via RMS-CASPT2-Driven Quantum-Classical Trajectories

We present a nonadiabatic molecular dynamics study of the ultrafast processes occurring in uracil upon UV light absorption, leading to electronic excitation and subsequent nonradiative decay. Previous studies have indicated that the mechanistic details of this process are drastically different depending on whether the process takes place in the gas phase, acetonitrile, or water. However, such results have been produced using quantum chemical methods that did not incorporate both static and dynamic electron correlation. In order to assess the previously proposed mechanisms, we simulate the photodynamics of uracil in the three environments mentioned above using quantum-classical trajectories and, for solvated uracil, hybrid quantum mechanics/molecular mechanics (QM/MM) models driven by the rotated multistate complete active space second-order perturbation (RMS-CASPT2) method. To do so, we exploit the gradient recently made available in OpenMolcas and compare the results to those obtained using the complete active space self-consistent field (CASSCF) method only accounting for static electron correlation. We show that RMS-CASPT2 produces, in general, a mechanistic picture different from the one obtained at the CASSCF level but confirms the hypothesis advanced on the basis of previous ROKS and TDDFT studies thus highlighting the importance of incorporating dynamic electron correlation in the investigation of ultrafast electronic deactivation processes.

A Nonorthogonal Configuration Interaction Approach to Singlet Fission in Perylenediimide Compounds

Perylenediimide molecules constitute a family of chromophores that undergo singlet fission, a process in which an excited singlet state converts into lower energy triplets on two neighboring molecules, potentially increasing the efficiency of organic solar cells. Here, the nonorthogonal configuration interaction method is applied to study the effect of the different crystal packing of various perylenediimide derivatives on the relative energies of the singlet and triplet states, the intermolecular electronic couplings, and the relative rates for singlet fission. The analysis of the wave functions and electronic couplings reveals that charge transfer states play an important role in the singlet fission mechanism. Dimer conformations where the PDI molecules are at large displacements along the long axis and short on the short axis are posed as the most favorable for singlet fission. The role of the substituent at the imide group has been inspected concluding that, although it has no effect in the energies, for some conformations it significantly influences the electronic couplings, and therefore, replacing this substituent with hydrogen may introduce artifacts in the computational modeling of the PDI molecules.