Complete Active Space Self-consistent Field Research Articles

MBE-CASSCF Approach for the Accurate Treatment of Large Active Spaces.

We present a novel implementation of the complete active space self-consistent field (CASSCF) method that makes use of the many-body expanded full configuration interaction (MBE-FCI) method to incrementally approximate electronic structures within large active spaces. On the basis of a hybrid first-order algorithm employing both Super-CI and quasi-Newton strategies for the optimization of molecular orbitals, we demonstrate both computational efficacy and high accuracy of the resulting MBE-CASSCF method. We assess the performance of our implementation on a set of established numerical tests before applying MBE-CASSCF in the investigation of the triplet-quintet spin gap of iron(II) porphyrin with active spaces as large as 50 electrons in 50 orbitals.

Chelate Complexes of 3d Transition Metal Ions─A Challenge for Electronic-Structure Methods?

Different electronic-structure methods were assessed for their ability to predict two important properties of the industrially relevant chelating agent nitrilotriacetic acid (NTA): its selectivity with respect to six different first-row transition metal ions and the spin-state energetics of its complex with Fe(III). The investigated methods encompassed density functional theory (DFT), the random phase approximation (RPA), coupled cluster (CC) theory, and the auxiliary-field quantum Monte Carlo (AFQMC) method, as well as the complete active space self-consistent field (CASSCF) method and the respective on-top methods: second-order N-electron valence state perturbation theory (NEVPT2) and multiconfiguration pair-density functional theory (MC-PDFT). Different strategies for selecting active spaces were explored, and the density matrix renormalization group (DMRG) approach was used to solve the largest active spaces. Despite somewhat ambiguous multi-reference diagnostics, most methods gave relatively good agreement with experimental data for the chemical reactions connected to the selectivity, which only involved transition-metal complexes in their high-spin state. CC methods yielded the highest accuracy followed by range-separated DFT and AFQMC. We discussed in detail that even higher accuracies can be obtained with NEVPT2, under the prerequisite that consistent active spaces along the entire chemical reaction can be selected, which was not the case for reactions involving Fe(III). A bigger challenge for electronic-structure methods was the prediction of the spin-state energetics, which additionally involved lower spin states that exhibited larger multi-reference diagnostics. Conceptually different, typically accurate methods ranging from CC theory via DMRG-NEVPT2 in combination with large active spaces to AFQMC agreed well that the high-spin state is energetically significantly favored over the other spin states. This was in contrast to most DFT functionals and RPA which yielded a smaller stabilization and some common DFT functionals and MC-PDFT even predicting the low-spin state to be energetically most favorable.

Electronic structure and optical spectral analysis of the MnO42- anion with consideration of site and Jahn-Teller distortion.

The ground state of 3d1 MnO42- was studied by density functional theory (DFT) and complete active space self-consistent field (CASSCF) methods in terms of a variety of molecular point group structures to ascertain the site and Jahn-Teller (JT) distortion effect. Modeling results from UB3LYP/6-31+G(d) calculations with natural bond orbital analysis show the four Mn-O bonds are coordinate covalent. The one-electron matrix elements from CASSCF(AILFT) (abinitio ligand field theory) with a second order perturbation treatment were used to calculate the parameters of the angular overlap model. These allowed the calculation of JT stabilization energies and first and second order JT coupling constants for MnO42- with C2v and D2d symmetry. Absorption spectra and excited state transition energies were calculated assuming the lattice distorted geometry was C2v with time-dependent DFT (TDDFT) using the unrestricted coulomb-attenuating UCAM-B3LYP density functional and with the spin-adapted spin-flip DFT using the UBH&HLYP density functional, both with the AUG-CC-PVTZ basis set. The mean absolute deviation of the calculations from experimental excited states for the ligand field and ligand-to-metal charge transfer (LMCT) bands was better than 0.1eV. Several new assignments for LMCT excited states were made on the basis of the TDDFT excitation results.

A MASH simulation of the photoexcited dynamics of cyclobutanone.

In response to a community prediction challenge, we simulate the nonadiabatic dynamics of cyclobutanone using the mapping approach to surface hopping (MASH). We consider the first 500fs of relaxation following photoexcitation to the S2 state and predict the corresponding time-resolved electron-diffraction signal that will be measured by the planned experiment. 397 abinitio trajectories were obtained on the fly with state-averaged complete active space self-consistent field using a (12,11) active space. To obtain an estimate of the potential systematic error, 198 of the trajectories were calculated using an aug-cc-pVDZ basis set and 199 with a 6-31+G* basis set. MASH is a recently proposed independent trajectory method for simulating nonadiabatic dynamics, originally derived for two-state problems. As there are three relevant electronic states in this system, we used a newly developed multi-state generalization of MASH for the simulation: the uncoupled spheres multi-state MASH method (unSMASH). This study, therefore, serves both as an investigation of the photodissociation dynamics of cyclobutanone, and also as a demonstration of the applicability of unSMASH to abinitio simulations. In line with previous experimental studies, we observe that the simulated dynamics is dominated by three sets of dissociation products, C3H6 + CO, C2H4 + C2H2O, and C2H4 + CH2 + CO, and we interpret our predicted electron-diffraction signal in terms of the key features of the associated dissociation pathways.

Self-Consistent Field Approach for the Variational Quantum Eigensolver: Orbital Optimization Goes Adaptive.

We present a self-consistent field (SCF) approach within the adaptive derivative-assembled problem-tailored ansatz variational quantum eigensolver (ADAPT-VQE) framework for efficient quantum simulations of chemical systems on near-term quantum computers. To this end, our ADAPT-VQE-SCF approach combines the idea of generating an ansatz with a small number of parameters, resulting in shallow-depth quantum circuits with a direct minimization of an energy expression that is correct to second order with respect to changes in the molecular orbital basis. Our numerical analysis, including calculations for the transition-metal complex ferrocene [Fe (C5H5)2], indicates that convergence in the self-consistent orbital optimization loop can be reached without a considerable increase in the number of two-qubit gates in the quantum circuit by comparison to a VQE optimization in the initial molecular orbital basis. Moreover, the orbital optimization can be carried out simultaneously within each iteration of the ADAPT-VQE cycle. ADAPT-VQE-SCF thus allows us to implement a routine analogous to the complete active space SCF, a cornerstone of state-of-the-art computational chemistry, in a hardware-efficient manner on near-term quantum computers. Hence, ADAPT-VQE-SCF paves the way toward a paradigm shift for quantitative quantum-chemistry simulations on quantum computers by requiring fewer qubits and opening up for the use of large and flexible atomic orbital basis sets in contrast to earlier methods that are predominantly based on the idea of full active spaces with minimal basis sets.

A CASSCF/MRCI trajectory surface hopping simulation of the photochemical dynamics and the gas phase ultrafast electron diffraction patterns of cyclobutanone.

We present the simulation of the photochemical dynamics of cyclobutanone induced by the excitation of the 3s Rydberg state. For this purpose, we apply the complete active space self-consistent field method together with the spin-orbit multireference configuration interaction singles treatment, combined with the trajectory surface hopping for the inclusion of nonadiabatic effects. The simulations were performed in the spin-adiabatic representation, including nine electronic states derived from three singlet and two triplet spin-diabatic states. Our simulations reproduce the two previously observed primary dissociation channels: the C2 pathway yielding C2H4 + CH2CO and the C3 pathway producing c-C3H6 + CO. In addition, two secondary products, CH2 + CO from the C2 pathway and C3H6 from the C3 pathway, both of them previously reported, are also observed in our simulation. We determine the ratio of the C3:C2 products to be about 2.8. Our findings show that most of the trajectories reach their electronic ground state within 200fs, with dissociation events finished after 300fs. We also identify the minimum energy conical intersections that are responsible for the relaxation and provide an analysis of the photochemical reaction mechanism based on multidimensional scaling. Furthermore, we demonstrate a minimal impact of triplet states on the photodissociation mechanism within the observed timescale. In order to provide a direct link to experiments, we simulate the gas phase ultrafast electron diffraction patterns and connect their features to the underlying structural dynamics.

Can Iron-Porphyrins Behave as Single-Molecule Magnets?

The study of magnetic properties, especially the magnetic anisotropy of iron-porphyrin complexes employing multiconfigurational methods, is quite challenging due to many strongly correlated electrons in nearly degenerate orbitals. However, a prerequisite for observing the magnetic anisotropy and slow magnetization relaxation, the zero-field splitting parameter, D, was experimentally observed decades ago for halide-based axially ligated penta-coordinate Fe(III)-porphyrins. In these complexes, the signs of D were reported mostly as positive; in a few cases, inconclusive signs of the D parameter were also mentioned. However, no ab initio calculations have been reported to shed light on this. Deciphering the electronic structure of these penta-coordinated complexes employing the complete active space self-consistent field method and N-electron valence second-order perturbation theory, we confirm the positive D values. However, a negative D value is highly desired to observe the single-molecule magnet properties without an external magnetic field, which we observed in the Fe(II)-porphyrin complexes with axial imidazole ligands instead of halide ligands. The detailed analysis of the multireference wave functions unravels the role of axial ligands in determining the sign and magnitude of the D parameters.

Zero-Field Splitting Tensor of the Triplet Excited States of Aromatic Molecules: A Valence Full-π Complete Active Space Self-Consistent Field Study.

A method to predict the D tensor in the molecular frame with multiconfigurational wave functions in large active space was proposed, and the spin properties of the lowest triplets of aromatic molecules were examined with full-π active space; such calculations were challenging because the size of active space grows exponentially with the number of π electrons. In this method, the exponential growth of complexity is resolved by the density matrix renormalization group (DMRG) algorithm. From the D tensor, we can directly determine the direction of the magnetic axes and the ZFS parameters, D- and E-values, of the phenomenological spin Hamiltonian with their signs, which are not usually obtained in ESR experiments. The method using the DMRG-CASSCF wave function can give correct results even when the sign of D- and E-values is sensitive to the accuracy of the prediction of the D tensor and existing methods fail to predict the correct magnetic axes.

Photochemical Mechanisms of Hydroxyquinoline Benzimidazole: Insights from Electronic Structure Calculations and Nonadiabatic Dynamics Simulations.

Excited-state intramolecular double proton transfer (ESIDPT) has received much attention because of its widespread existence in the life reactions of living organisms, and materials with this property are significant for their special luminescent properties. In this work, the complete active space self-consistent field (CASSCF) and OM2/multireference configuration interaction (OM2/MRCI) methods have been employed to study the static electronic structure calculations of the photochemistry and the possibility of ESIDPT process of hydroxyquinoline benzimidazole (HQB) molecule, along with the nonadiabatic dynamics simulations. The computational results show that the HQB molecule is relaxed to the S1-ENOL minimum after being excited to the Franck-Condon point in the S1 state. Subsequently, during the nonadiabatic deactivation process, the OH···N proton transfer and the twisting of benzimidazole occur before arriving at the single proton transfer conical intersection S1S0-KETO. Finally, the system can either return to the initial ground-state structure S0-ENOL or to the single proton transfer ground-state structure S0-KETO, both of which have almost the same probability. The dynamics simulations also show that no double proton transfer occurs. The excited-state lifetime of HQB is fitted to 1.1 ps, and only 64% of the dynamic trajectories return to the ground state within the 2.0 ps simulation time. We hope the detailed reaction mechanism of the HQB molecule will provide new insights into similar systems.

Deciphering the Internal Conversion Processes Involved in the Photochemical Ring-Opening of 1,3-Cyclohexadiene by Symmetric sp2-Carbon Substitutions.

Chemical substituents hold the potential to markedly influence the photochemical behavior in molecular systems and assist in gaining a comprehensive understanding of nonadiabatic phenomena. In this study, we have conducted a comparative analysis of the influence of chemical substituents on the photochemical ring-opening of 1,3-cyclohexadiene (CHD), considering four systems: CHD, 2,3-dimethylcyclohexadiene (CHD-Me2-1), 1,4-dimethylcyclohexadiene (CHD-Me2-2), and 1,2,3,4-tetramethylcyclohexadiene (CHD-Me4), using electronic structure theory calculations and nonadiabatic molecular dynamics simulations. Employing extended multistate complete active space second-order perturbation (XMS-CASPT2) theory, we optimized reactants, S1 states, conical intersections (CIns), and products, revealing structural and energetic variations consistent with prior research. Nonadiabatic molecular dynamics simulation was used to gain insights into photochemical dynamics at state-averaged complete active space self-consistent field (SA-CASSCF) theory. CHD-Me4 exhibited reduced carbon-carbon single bond rupture rates, responsible for ring-opening, due to substituent proximity. Further, CHD-Me2-2 and CHD-Me4 displayed prolonged excited-state relaxation times, highlighting notable substituents' impact. Analysis of kinetic energy profiles of specific carbon atoms also revealed restrained atomic displacements, particularly in CHD-Me2-2 and CHD-Me4. These findings advance our understanding of how substituents modulate photochemical reactions in cyclohexadiene derivatives, guiding new molecular design and future research.

Metal-Free Activation of Molecular Oxygen by Quaternary Ammonium-Based Ionic Liquid: A Detail Mechanistic Study.

Most oxidation processes in common organic synthesis and chemical biology require transition metal catalysts or metalloenzymes. Herein, we report a detailed mechanistic study of a metal-free oxygen (O2) activation protocol on benzylamine/alcohols using simple quaternary alkylammonium-based ionic liquids to produce products such as amide, aldehyde, imine, and in some cases, even aromatized products. NMR and various control experiments established the product formation and reaction mechanism, which involved the conversion of molecular oxygen into a hydroperoxyl radical via a proton-coupled electron transfer process. Detection of hydrogen peroxide in the reaction medium using colorimetric analysis supported the proposed mechanism of oxygen activation. Furthermore, first-principles calculations using density functional theory (DFT) revealed that reaction coordinates and transition state spin densities have a unique spin conversion of triplet oxygen leading to formation of singlet products via a minimum energy crossing point. In addition to DFT, domain-based local pair natural orbital coupled cluster, (DLPNO-CCSD(T)), and complete active space self-consistent field, CASSCF(20,14) methods complemented the above findings. Partial density of states analysis showed stabilization of π* orbital of oxygen in the presence of ionic liquid, making it susceptible to hydrogen abstraction in a mild, metal-free condition. Inductively coupled plasma atomic emission spectroscopic (ICP-AES) analysis of reactant and ionic liquids clearly showed the absence of any significant transition metal contamination. The current results described the origin of O2 activation within the context of molecular orbital (MO) theory and opened up a new avenue for the use of ionic liquids as inexpensive, multifunctional and high-performance alternative to metal-based catalysts for O2 activation.

Isomerization dynamics of a novel cis/trans‐only merocyanine

AbstractMerocyanines (MC) usually adopt ring opened zwitterionic structures that are interconvertible with their ring‐closed spiropyran photoisomers. By methylating the phenolate oxygen, and thereby blocking the ring‐closure reaction, a cis/trans‐only MC photoswitch was obtained, yielding a perfect candidate for a detailed examination of the cis/trans isomerization mechanism for this class of compounds. This photoswitch displays outstanding properties including excellent photoreaction quantum yields and photoswitching turnovers. Due to the central polymethine bridge of MC, in principle eight cis (C)/trans/(T) isomers are possible. Density Functional Theory (DFT) calculations revealed the CCT and TTT‐isomers of the studied compound as most stable cis and trans ground state isomers, respectively. UV/vis transient absorption studies combined with conical intersection computations with the complete active space self‐consistent field (CASSCF) method show that both trans/cis‐ and cis/trans‐photoisomerizations are initiated by a rotation of the central doubled bond fragment. A hot ground state species is then formed, which undergoes a second isomerization. Thus, the cis/trans reaction proceeds via a CCT‐CTT‐TTT sequence and the reverse reaction via TTT‐TCT‐CCT.

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.

Spectroscopic and Theoretical Investigation of High-Spin Square-Planar and Trigonal Fe(II) Complexes Supported by Fluorinated Alkoxides.

The electronic structures and spectroscopic behavior of three high-spin FeII complexes of fluorinated alkoxides were studied: square-planar {K(DME)2}2[Fe(pinF)2] (S) and quasi square-planar {K(C222)}2[Fe(pinF)2] (S') and trigonal-planar {K(18C6)}[Fe(OC4F9)3] (T) where pinF = perfluoropinacolate and OC4F9 = tris-perfluoro-t-butoxide. The zero-field splitting (ZFS) and hyperfine structure parameters of the S = 2 ground states were determined using field-dependent 57Fe Mössbauer and high-field and -frequency electron paramagnetic resonance (HFEPR) spectroscopies. The spin Hamiltonian parameters were analyzed with crystal field theory and corroborated by density functional theory (DFT) and ab initio complete active space self-consistent field (CASSCF) calculations. Whereas the ZFS tensor of S has a small rhombicity, E/D = 0.082, and a positive D = 15.17 cm-1, T exhibits a negative D = -9.16 cm-1 and a large rhombicity, E/D = 0.246. Computational investigation of the structural factors suggests that the ground-state electronic configuration and geometry of T's Fe site are determined by the interaction of [Fe(OC4F9)3]- with {K(18C6)}+. In contrast, two distinct countercations of S/S' have a negligible influence on their [Fe(pinF)2]2- moieties. Instead, the distortions in S' are likely induced by the chelate ring conformation change from δλ, observed for S, to the δδ conformation, determined for S'.

Dinitrogen Activation Mediated by the (P2PPh)Fe Complex: Electronic Structure, Dimerization Mechanism, and Magnetic Coupling.

Herein, we report the estimation of the extent of dinitrogen activation by different charged and structural forms of (P2PPh)Fe biomimetic catalysts, which, in the presence of light, exhibit significant yield in the N2-to-NH3 conversion. Complete active space self-consistent field (CASSCF) calculations have been used to determine the electronic structure of different reduced forms of the mononuclear complexes: the neutral (P2PPh)Fe(N2)2 adduct and the anionic [(P2PPh)Fe(N2)]- and [(P2PPh)Fe(N2)]2- complexes. These calculations also revealed that the extent of reduction of a dinitrogen molecule reaches up to one electron (N21-) due to the back-bonding from the Fe center, in agreement with the changes observed in the vibration frequency of the N-N bond in these complexes. In addition, the energy profile of the dimerization of the mononuclear (P2PPh)Fe(N2)2 complex to the dinuclear mono-N2-bridged [(P2PPh)Fe]2(μ-N2) complex has been determined by means of density functional theory (DFT) calculations. A three-step mechanism has been proposed for the dimerization, favored by both kinetics and thermodynamics criteria. Finally, the magnetic coupling constant in the diiron (μ-N2) complex is estimated from CASSCF/NEVPT2 calculations. Such a dinuclear complex presents a strong antiferromagnetic coupling resulting from the interaction between two S = 1 d6 Fe2+ ions, bridged by a highly activated dinitrogen molecule (N22-) with two electrons on the π* orbitals.

A Local Diabatisation Method for Two-State Adiabatic Conical Intersections.

A methodology to locally characterize conical intersections (CIs) between two adiabatic electronic states for which no nonadiabatic coupling (NAC) vectors are available is presented. Based on the Hessian and gradient at the CI, the branching space coordinates are identified. The potential energy surface around the CI in the branching space is expressed in the diabatic representation, from which the NAC vectors can be calculated in a wave-function-free, energy-based approach. To demonstrate the universality of the developed methodology, the minimum-energy CI (MECI) between the first (S1) and second (S2) singlet excited states of formamide is investigated at the state-averaged complete active space self-consistent field (SA-CASSCF) and extended multistate complete active space second-order perturbation theory (XMS-CASPT2) levels of theory. In addition, the asymmetrical MECI between the ground state (S0) and S1 of cyclopropanone is evaluated using SA-CASSCF, as well as (ME)CIs between the S1 and S2 states of benzene using SA-CASSCF and time-dependent density functional theory (TDDFT). Finally, a CI between the S1 and S2 excited states of thiophene was analyzed using TDDFT.

Study of electronic structure and optical transition properties of low-lying excited states of AuB molecules based on configuration interaction method

High-level configuration interaction method including the spin-orbit coupling is used to investigate the low-lying excited electronic states of AuB that is not reported experimentally. The electronic structure in our work is preformed through the three steps stated below. First of all, Hartree-Fock method is performed to compute the singlet-configuration wavefunction as the initial guess. Next, we generate a multi-reference wavefunction by using the state-averaged complete active space self-consistent field (SACASSCF). Finally, the wavefunctions from CASSCF are utilized as reference, the exact energy point values are calculated by the explicitly correlated dynamic multi-reference configuration interaction method (MRCI). The Davidson correction (+Q) is put forward to solve the size-consistence problem caused by the MRCI method. To ensure the accuracy, the spin-orbit effect and correlation for inner shell electrons and valence shell electrons are considered in our calculation. The potential energy curves of 12 Λ-S electronic states are obtained. According to the explicit potential energy curves, we calculate the spectroscopic constants through solving radial Schrödinger equation numerically. We analyze the influence of electronic state configuration on the dipole moment by using the variation of dipole moment with nuclear distance. The spin-orbit matrix elements for parts of low-lying exciting states are computed, and the relation between spin-orbit coupling and predissociation is discussed. The predissociation is analyzed by using the obtained spin-orbit matrix elements of the 4 Λ-S states which spilt into 12 Ω states. It indicates that due to the absence of the intersections between the curves of spin-orbit matrix elements related with the 4 low-lying Λ-S states, the predissociation for these low-lying exciting states will not occur. Finally, the properties of optical transition between the ground Ω state <inline-formula><tex-math id="M3">\begin{document}$ {\rm A}^{1}{{{\Pi}}}_{1} $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="5-20231347_M3.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="5-20231347_M3.png"/></alternatives></inline-formula> and first excited Ω state <inline-formula><tex-math id="M4">\begin{document}$ {{\mathrm{X}}}^{1}{{{\Sigma }}}_{{0}^{+}} $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="5-20231347_M4.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="5-20231347_M4.png"/></alternatives></inline-formula> are discussed in laser-cooling filed by analyzing the Franck-Condon factors and radiative lifetime. And the transition dipole moment is also calculated. But our results reveal that the AuB is not an ideal candidate for laser-cooling. In conclusion, this work is helpful in deepening the understanding of AuB, especially the structures of electronic states, interaction between excited states, and optical transition properties. All the data presented in this paper are openly available at <ext-link ext-link-type="uri" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="https://www.doi.org/10.57760/sciencedb.j00213.00009">https://www.doi.org/10.57760/sciencedb.j00213.00009</ext-link>.

Spinless formulation of linearized adiabatic connection approximation and its comparison with second order N-electron valence state perturbation theory

The adiabatic connection (AC) approximation, along with its linearized variant AC0, was introduced as a method of obtaining dynamic correlation energy. When using a complete active space self-consistent field (CASSCF)...

Changing aromatic properties through stacking: the face-to-face dimer of Ni(II) bis(pentafluorophenyl)norcorrole.

Nuclear magnetic resonance (NMR) shielding constants have been calculated for Ni(II) bis(pentafluorophenyl)norcorrole and its face-to-face stacked dimer at the Hartree-Fock (HF), second-order Møller-Plesset perturbation theory (MP2), complete-active-space self-consistent-field (CASSCF) levels as well as at density functional theory (DFT) levels using several functionals. The calculated 1H NMR shielding constants agree rather well with the experimental ones. The shielding constants of N and Ni calculated at DFT, HF, and MP2 levels differ from those obtained in the CASSCF calculations due to near-degeneracy effects at the Ni atom. The calculated magnetically induced current densities show that the monomer is antiaromatic, sustaining a strong global paratropic ring current, and the dimer is aromatic, sustaining a strong diatropic ring current. Qualitatively the same current density is obtained at the employed levels of theory. The most accurate ring-current strengths are probably obtained at the MP2 level. The aromatic dimer has a short intermolecular distance of less than 3 Å. The intermolecular interaction changes the nature of the frontier orbitals leading to a formal double bond between the norcorrole macrocycles.

Performance of point charge embedding schemes for excited states in molecular organic crystals.

Modeling excited state processes in molecular crystals is relevant for several applications. A popular approach for studying excited state molecular crystals is to use cluster models embedded in point charges. In this paper, we compare the performance of several embedding models in predicting excited states and S1-S0 optical gaps for a set of crystals from the X23 molecular crystal database. The performance of atomic charges based on ground or excited states was examined for cluster models, Ewald embedding, and self-consistent approaches. We investigated the impact of various factors, such as the level of theory, basis sets, embedding models, and the level of localization of the excitation. We consider different levels of theory, including time-dependent density functional theory and Tamm-Dancoff approximation (TDA) (DFT functionals: ωB97X-D and PBE0), CC2, complete active space self-consistent field, and CASPT2. We also explore the impact of selection of the QM region, charge leakage, and level of theory for the description of different kinds of excited states. We implemented three schemes based on distance thresholds to overcome overpolarization and charge leakage in molecular crystals. Our findings are compared against experimental data, G0W0-BSE, periodic TDA, and optimally tuned screened range-separated functionals.