Advancement and Application of TDDFT-based Non-Adiabatic Molecular Dynamics Methods for Triplet States

非绝热分子动力学（NAMD）模拟是研究飞秒时间尺度上发生的超快光致过程的重要理论工具。以往的研究主要集中在电子跃迁过程中不涉及电子自旋的变化（单线态激发和内部转换）。涉及自旋逆转的过程（三重态激发和系统间交叉，ISC）通常在更长的时间尺度上进行。然而，越来越多的人发现了这一规则的重要例外，例如过渡金属配合物中的超快ISC过程，这与染料敏化太阳能电池或oled的发展有关。然而，在NAMD模拟中考虑三重态也会引入各种问题，其中一些问题将在本项目中解决。对于较大的系统，时间相关密度泛函理论（TDDFT）是计算NAMD模拟所需的量子化学量的流行方法，因为它提供了一个有吸引力的性价比。然而，传统的密度泛函对三重态激发能有较大的误差，这可能导致NAMD模拟结果不正确。另一方面，相对较新的局部混合泛函在描述三重态方面具有显著的优势，因此第一个目标是为这些泛函开发非绝热耦合矩阵元素，使它们可用于NAMD模拟。为了描述NAMD模拟中的ISC，必须描述潜在的自旋轨道耦合（SOC）效应。为了效率和简单起见，通常采用摄动方法。然而，对于强耦合，例如在重元素的复合体中预期的耦合，这些方法可能会失败。该项目的第二个目标是使用变分的双组分TDDFT方法，这种方法对于强耦合更可靠，可以确定NAMD模拟中的SOC矩阵元素。为了在NAMD模拟中计算控制原子核运动的力，需要激发态的核梯度。然而，当代的NAMD方法忽略了SOC项的梯度，因为它们还不能被任何量子化学程序解析计算。该项目的第三个目标是在双分量TDDFT框架中推导和实现这些梯度。这些新方法将用于探索与太阳能收集高度相关的分子中的ISC过程，但迄今为止只使用更简单的方法进行了研究，或者根本没有进行过研究。在这种情况下，还将比较不同方法的准确性。

Non-adiabatic molecular dynamics (NAMD) simulations are an important theoretical tool to study ultrafast light-induced processes occurring on a femtosecond time scale. Previous research has been predominantly focused on processes in which the electronic transitions involve no change in electron spin (singlet excitations and internal conversion). Processes involving spin reversal (triplet excitations and intersystem crossing, ISC) typically proceed on a much longer time scale. However, important exceptions to this rule are increasingly found, such as ultrafast ISC processes in transition-metal complexes, which are relevant for the development of dye-sensitized solar cells or OLEDs. However, the consideration of triplet states in NAMD simulations also introduces various problems, some of which will be addressed in this project.For larger systems, time-dependent density functional theory (TDDFT) is a popular method for calculating the quantum chemical quantities needed for NAMD simulations, as it offers an attractive cost-performance ratio. However, conventional density functionals exhibit large errors for triplet excitation energies, which can lead to incorrect results in NAMD simulations. The relatively recent class of local hybrid functionals, on the other hand, offers significant advantages in describing triplet states, so a first goal is to develop non-adiabatic coupling matrix elements for these functionals to make them available for NAMD simulations.To describe ISC in NAMD simulations, the underlying spin-orbit coupling (SOC) effect must be described. For efficiency and simplicity, this is usually done using perturbative approaches. However, for strong couplings, such as those expected in complexes with heavy elements, these approaches may fail. The second goal of this project is to instead use variational, two-component TDDFT methods that are more reliable for strong couplings to determine SOC matrix elements in NAMD simulations.For the calculation of the forces governing the motion of atomic nuclei in NAMD simulations, the nuclear gradients of the excited states are needed. However, contemporary NAMD methods neglect the gradients of the SOC terms because they are not yet analytically computable by any quantum chemical program. The third goal of this project is to derive and implement these gradients in the framework of two-component TDDFT.The new methods will be used to explore ISC processes in molecules that are highly relevant in the context of solar energy harvesting but have so far only been studied using simpler approaches or not at all. The accuracy of the different methods will also be compared in this context.