Removing the Deadwood from DFT/MRCI Wave Functions: The p-DFT/MRCI Method

Minimum Energy Conical Intersection Optimization Using DFT/MRCI(2)

The combined density functional theory and multireference configuration interaction (DFT/MRCI) method is a semiempirical electronic structure approach that is both computationally efficient and has predictive accuracy for the calculation of electronic excited states and for the simulation of electronic spectroscopies. However, given that the reference space is generated via a selected-CI procedure, a challenge arises in the construction of smooth potential energy surfaces. To address this issue, we treat the local discontinuities that arise as noise within the Gaussian progress regression framework and learn the surfaces by explicitly incorporating and optimizing a white-noise kernel. The characteristic polynomial coefficient surfaces of the potential matrix, which are smooth functions of nuclear coordinates even at conical intersections, are learned in place of the adiabatic energies and are used to optimize the DFT/MRCI(2) minimum energy conical intersection geometries for representative intersection motifs in the molecules ethylene, butadiene, and fulvene. One consequence of explicitly treating the noise in the surfaces is that the energy difference cannot be made arbitrarily small at points of nominal intersection. Despite the limitations, however, we find the structures as well as the branching spaces to compare well with ab initio MRCI and conclude that this approach is a viable method to learn a smooth representation of DFT/MRCI(2) surfaces.

A DFT/MRCI Hamiltonian parameterized using only ab initio data. II. Core-excited states

A newly parameterized combined density functional theory and multi-reference configuration interaction (DFT/MRCI) Hamiltonian, termed core-valence separation (CVS)-QE12, is defined for the computation of K-shell core-excitation and core-ionization energies. This CVS counterpart to the recently reported QE8 Hamiltonian [Costain et al., J. Chem. Phys, 160, 224106 (2024)] is parameterized by fitting to benchmark quality ab initio data. The definition of the CVS-QE12 and QE8 Hamiltonians differ from previous CVS-DFT/MRCI parameterizations in three primary ways: (i) the replacement of the BHLYP exchange-correlation functional with QTP17 to yield a balanced description of both core and valence excitation energies, (ii) the adoption of a new, three-parameter damping function, and (iii) the introduction of separate scaling of the core-valence and valence-valence Coulombic interactions. Crucially, the parameters of the CVS-QE12 Hamiltonian are obtained via fitting exclusively to highly accurate ab initio vertical core-excitation and ionization energies computed at the CVS-EOM-CCSDT level of theory. The CVS-QE12 Hamiltonian is validated against further benchmark computations and is found to furnish K-edge core vertical excitation and ionization energies exhibiting absolute errors ≤0.5 eV at low computational cost.

A new parameterization of the DFT/CIS method with applications to core-level spectroscopy

Time-dependent density functional theory (TD-DFT) within a restricted excitation space is an efficient means to compute core-level excitation energies using only a small subset of the occupied orbitals. However, core-to-valence excitation energies are significantly underestimated when standard exchange-correlation functionals are used, which is partly traceable to systemic issues with TD-DFT's description of Rydberg and charge-transfer excited states. To mitigate this, we have implemented an empirically modified combination of configuration interaction with single substitutions (CIS) based on Kohn-Sham orbitals, which is known as "DFT/CIS." This semi-empirical approach is well-suited for simulating x-ray near-edge spectra, as it contains sufficient exact exchange to model charge-transfer excitations yet retains DFT's low-cost description of dynamical electron correlation. Empirical corrections to the matrix elements enable semi-quantitative simulation of near-edge x-ray spectra without the need for significant a posteriori shifts; this should be useful in complex molecules and materials with multiple overlapping x-ray edges. Parameter optimization for use with a specific range-separated hybrid functional makes this a black-box method intended for both core and valence spectroscopy. Results herein demonstrate that realistic K-edge absorption and emission spectra can be obtained for second- and third-row elements and 3d transition metals, with promising results for L-edge spectra as well. DFT/CIS calculations require absolute shifts that are considerably smaller than what is typical in TD-DFT.

Calculation of quasi-diabatic states within the DFT/MRCI(2) framework: The QD-DFT/MRCI(2) method

We describe a procedure for the calculation of quasi-diabatic states within the recently introduced DFT/MRCI(2) framework [S. P. Neville and M. S. Schuurman, J. Chem. Phys. 157, 164103 (2022)]. Based on an effective Hamiltonian formalism, the proposed procedure, which we term QD-DFT/MRCI(2), has the advantageous characteristics of being simultaneously highly efficient and effectively black box in nature while directly yielding both quasi-diabatic potentials and wave functions of high quality. The accuracy and efficiency of the QD-DFT/MRCI(2) formalism are demonstrated via the simulation of the vibronic absorption spectra of furan and chlorophyll a.

A DFT/MRCI Hamiltonian parameterized using only ab initio data: I. valence excited states

A new combined density functional theory and multi-reference configuration interaction (DFT/MRCI) Hamiltonian parameterized solely using the benchmark ab initio vertical excitation energies obtained from the QUEST databases is presented. This new formulation differs from all previous versions of the method in that the choice of the underlying exchange-correlation (XC) functional employed to construct the one-particle (orbital) basis is considered, and a new XC functional, QTP17, is chosen for its ability to generate a balanced description of core and valence vertical excitation energies. The ability of the new DFT/MRCI Hamiltonian, termed QE8, to furnish accurate excitation energies is confirmed using benchmark quantum chemistry computations, and a mean absolute error of 0.16 eV is determined for the wide range of electronic excitations included in the validation dataset. In particular, the QE8 Hamiltonian dramatically improves the performance of DFT/MRCI for doubly excited states. The performance of fast approximate DFT/MRCI methods, p-DFT/MRCI and DFT/MRCI(2), is also evaluated using the QE8 Hamiltonian, and they are found to yield excitation energies in quantitative agreement with the parent DFT/MRCI method, with the two methods exhibiting a mean difference of 0.01 eV with respect to DFT/MRCI over the entire benchmark set.

A laboratory frame density matrix for ultrafast quantum molecular dynamics

In most cases, the ultrafast dynamics of resonantly excited molecules are considered and almost always computed in the molecular frame, while experiments are carried out in the laboratory frame. Here, we provide a formalism in terms of a lab frame density matrix, which connects quantum dynamics in the molecular frame to those in the laboratory frame, providing a transparent link between computation and measurement. The formalism reveals that in any such experiment, the molecular frame dynamics vary for molecules in different orientations and that certain coherences, which are potentially experimentally accessible, are rejected by the orientation-averaged reduced vibronic density matrix. Instead, molecular angular distribution moments are introduced as a more accurate representation of experimentally accessible information. Furthermore, the formalism provides a clear definition of a molecular frame quantum tomography and specifies the requirements to perform such a measurement enabling the experimental imaging of molecular frame vibronic dynamics. Successful completion of such a measurement fully characterizes the molecular frame quantum dynamics for a molecule at any orientation in the laboratory frame.

Carbon Nanodots from an In Silico Perspective

Carbon nanodots (CNDs) are the latest and most shining rising stars among photoluminescent (PL) nanomaterials. These carbon-based surface-passivated nanostructures compete with other related PL materials, including traditional semiconductor quantum dots and organic dyes, with a long list of benefits and emerging applications. Advantages of CNDs include tunable inherent optical properties and high photostability, rich possibilities for surface functionalization and doping, dispersibility, low toxicity, and viable synthesis (top-down and bottom-up) from organic materials. CNDs can be applied to biomedicine including imaging and sensing, drug-delivery, photodynamic therapy, photocatalysis but also to energy harvesting in solar cells and as LEDs. More applications are reported continuously, making this already a research field of its own. Understanding of the properties of CNDs requires one to go to the levels of electrons, atoms, molecules, and nanostructures at different scales using modern molecular modeling and to correlate it tightly with experiments. This review highlights different in silico techniques and studies, from quantum chemistry to the mesoscale, with particular reference to carbon nanodots, carbonaceous nanoparticles whose structural and photophysical properties are not fully elucidated. The role of experimental investigation is also presented. Hereby, we hope to encourage the reader to investigate CNDs and to apply virtual chemistry to obtain further insights needed to customize these amazing systems for novel prospective applications.