Nonorthogonal Multireference Wave Function Description of Triplet-Triplet Energy Transfer Couplings.

In this study, the use of self-consistent field quasi-diabats is investigated for calculation of triplet energy transfer diabatic coupling elements. It is proposed that self-consistent field quasi-diabats are particularly useful for studying energy transfer (EnT) processes because orbital relaxation in response to changes in electron configuration is implicitly built into the model. The conceptual model that is developed allows for the simultaneous evaluation of direct and charge-transfer mechanisms to establish the importance of the different possible EnT mechanisms. The method's performance is evaluated using two model systems: the ethylene dimer and ethylene with the methaniminium cation. While states that mediate the charge-transfer mechanism were found to be higher in energy than the states involved in the direct mechanism, the coupling elements that control the kinetics were found to be significantly larger in the charge-transfer mechanism. Subsequently, we discuss the advantage of the approach in the context of practical difficulties with the use of established approaches.

Protonation state control of electric field induced molecular switching mechanisms.

Scanning tunneling microscopy tip-induced deprotonation has been demonstrated experimentally and can be used as an additional control mechanism in electric-field induced molecular switching. The goal of the current work is to establish whether (de)protonation can be used to inhibit or enhance the electric field controlled thermal and photoisomerization processes. Dihydroxyazobenzene is used as a model system, where protonation/deprotonation of the free hydroxyl moiety changes the azo bond order, and so modifies the rate of electric field induced isomerization. Through the combined action of deprotonation and applied field, it was found that the cis-to-trans thermal isomerization barrier could be completely removed, changing the isomerization half-life from the order of several months. In addition, due to the presence of multiple isomerization mechanisms, electric fields could modify the isomerization kinetics by increasing the number of energetically viable isomerization pathways, rather than reducing the activation barrier of the lowest energy pathway. Excited state calculations indicated that the protonation state and electric field could be used together to control the presence of electronic degeneracies along the rotation pathway between S0/S1, and along all three pathways between S1/S2. This work provides insight into the mechanisms that enable the use of protonation state, light, and electric fields in concert to control molecular switches.

Nonorthogonal Active Space Decomposition of Wave Functions with Multiple Correlation Mechanisms.

The nonorthogonal active space decomposition (NO-ASD) methodology is proposed for describing systems containing multiple correlation mechanisms. NO-ASD partitions the wave function by a correlation mechanism, such that the interactions between different correlation mechanisms are treated with an effective Hamiltonian approach, while interactions between correlated orbitals in the same correlation mechanism are treated explicitly. As a result, the determinant expansion scales polynomially with the number of correlation mechanisms rather than exponentially, which significantly reduces the factorial scaling associated with the size of the correlated orbital space. Despite the nonorthogonal framework of NO-ASD, the approach can take advantage of computational efficient matrix element evaluation when performing nonorthogonal coupling of orthogonal determinant expansions. In this work, we introduce and examine the NO-ASD approach in comparison to complete active space methods to establish how the NO-ASD approach reduces the problem dimensionality and the extent to which it affects the amount of correlation energy recovered. Calculations are performed on ozone, nickel-acetylene, and isomers of µ-oxo dicopper ammonia.

Role of Exact Exchange in Difference Projected Double-Hybrid Density Functional Theory for Treatment of Local, Charge Transfer, and Rydberg Excitations.

Difference approaches to the study of excited states have undergone a renaissance in recent years, with the development of a plethora of algorithms for locating self-consistent field approximations to excited states. Density functional theory is likely to offer the best balance of cost and accuracy for difference approaches, and yet there has been little investigation of how the parametrization of density functional approximations affects performance. In this work, we aim to explore the role of the global Hartree-Fock exchange parameter in tuning accuracy of different excitation types within the framework of the recently introduced difference projected double-hybrid density functional theory approach and contrast the performance with conventional time-dependent double-hybrid density functional theory. Difference projected double-hybrid density functional theory was demonstrated to give vertical excitation energies with average error and standard deviation with respect to multireference perturbation theory comparable to more expensive linear-response coupled cluster approaches ( J. Chem. Phys.2020, 153, 074103). However, despite benchmarking of local excitations, there has been no investigation of the methods performance for charge transfer or Rydberg excitations. In this work we report a new benchmark of charge transfer, Rydberg, and local excited state vertical excitation energies and examine how the exact Hartree-Fock exchange affects the benchmark performance to provide a deeper understanding of how projection and nonlocal correlation balance differing sources of error in the ground and excited states.

Oriented External Electric Field Tuning of Unsubstituted Azoheteroarene Thermal Isomerization Half-Lives.

Azoheteroarenes are relatively new photoswitchable compounds, where one of the phenyl rings of an azobenzene molecule is replaced by a heteroaromatic five-membered ring. Recent findings on methylated azoheteroarenes show that these photoswitches have potential in various optically addressable applications. The thermal stability of molecular switches is one of the primary factors considered in the design process. For molecular memory or energy storage devices, long thermal relaxation times are required. However, inducing a short thermal isomerization lifetime is required to release stored energy or as an alternative to photoswitching to avoid overlapping absorption spectra that reduce switching fidelity. In this study, we investigate how oriented external electric fields can be used to tune the thermal isomerization properties of three unsubstituted heteroaryl azo compounds-azoimidazole, azopyrazole, and azopyrrole. We show that favorable electric field orientations can increase the thermal half-life of studied molecules by as much as 60 times or reduce it from tens of days to seconds, compared to their half-life values in the field-free environment. A deeper understanding of the relationship between structure and kinetic properties provides insight as to how molecular switches can be designed for their electric field response in switching applications.

Difference projection-after-variation double-hybrid density functional theory applied to the calculation of vertical excitation energies.

The use of projection-after-variation double-hybrid density functional theory is proposed and examined as a difference method for the calculation of excited states. The strengths and weaknesses of the proposed method are discussed with particular reference to connections with linear response coupled-cluster theory. Vertical excitation energies are computed for the 28 molecule benchmark of Schreiber and co-workers in order to compare how the model performs with linear response coupled-cluster theories and multireference perturbation theory. The findings of this study show that the proposed method can achieve standard deviations in the error of computed vertical excitation energies compared to complete active space second-order perturbation theory of similar size to linear response coupled-cluster theories.

Global Elucidation of Self-Consistent Field Solution Space Using Basin Hopping.

Reliable global elucidation of (subsets of) self-consistent field solutions is required for continued development and application of computational approaches that utilize these solutions as reference wavefunctions. We report the derivation and implementation of a stochastic approach to perform global elucidation of self-consistent field solutions by exploiting the connection between global optimization and global elucidation problems. We discuss the design of the algorithm through combining basin-hopping search algorithms with a Lie algebraic approach to linearize self-consistent field solution space, while also allowing preservation of desired spin-symmetry properties of the wavefunction. The performance of the algorithm is demonstrated on minimal basis C2v H4 due to its use as a model system for global self-consistent field solution exploration algorithms. Subsequently, we show that the model is capable of successfully identifying low-lying self-consistent solutions of benzene and NO2 with polarized double-zeta and triple-zeta basis sets and examine the properties of these solutions.

Effect of Oriented External Electric Fields on the Photo and Thermal Isomerization of Azobenzene.

Azobenzene is a prototype molecule with potential applications in molecular switches, solar thermal batteries, sensors, photoresponsive membranes, molecular electronics, data storage, and nonlinear optics. Photo and thermal isomerization pathways exhibit different charge-transfer character and dipole moments, implying that the use of electric fields can be used to modulate the reactivity of azobenzene. This article examines the differential effect of orientated electric fields on the rotation and inversion thermal and photoisomerization pathways of azobenzene to explore the feasibility of using electric fields in the design of azobenzene-based molecular devices. Our findings demonstrate that the application of orientated electric fields modifies the accessibility of the S0/S1 seam of electronic degeneracy, as well as changes the energetically favored relaxation pathway in the branching space to yield different photoproducts. In addition, we observed strong-field dipole-inversion effects that cause a topographical change in the response of the potential energy surface to the applied field and can result in geometric minima that do not exist under field-free conditions. On the S0 surface, transition barriers can be modified on the order of ±10 kcal mol-1, enabling control of thermal isomerization rates.