The mixed-reference spin−flip time-dependent density functional theory (MRSF-TD-DFT) method eliminates the erroneous spin contamination of the SF-TD-DFT methodology, while retaining the conceptual and practical simplicity of the latter. The availability of the analytic gradient of the energy of the MRSF-TD-DFT response states enables automatic geometry optimization of the targeted states. Here, we apply the new method to optimize the geometry of several S 1 /S 0 conical intersections occurring in typical organic molecules. We demonstrate that MRSF-TD-DFT is capable of producing the correct double-cone topology of the intersections and describing the geometry of the lowest-energy conical intersections and their relative energies with accuracy matching that of the best multireference wavefunction ab initio methods. In this regard, MRSF-TD-DFT differs from many popular singlereference methods, such as, e.g., the linear response TD-DFT method, which fail to produce the correct topology of the intersections. As the new methodology completely eliminates the ambiguity with the identification of the response states as proper singlets or triplets, which is plaguing the SF-TD-DFT calculations, it can be used for automatic geometry optimization and molecular dynamic simulations not requiring constant human intervention.

Two-dimensional (2D) ferromagnets with high spin-polarization ratio and high Curie temperature are crucial for developing next-generation spintronic nanodevices. Using first-principles calculations, we predict two polymorphic modifications ( t-VN and h-VN) of 2D VN lattices that have robust intrinsic ferromagnetic properties and high Curie temperatures. Whereas t-VN has 99.9% of spin polarization at the Fermi level, h-VN possesses a half-metallic type of conductivity and keeps it after contact with semiconducting MoS, which can be used as the substrate for h-VN synthesis and valley polarized contacts. Magnetocrystalline anisotropy energy of 2D VN polymorphs is found to be at least an order larger than those of Fe and Ni bulks. The phonon spectra and ab initio molecular dynamic simulation prove that 2D VN lattices have a high thermodynamic stability. These advantages demonstrate that the VN monolayers should be promising candidates for low-dimensional spintronic devices.

Excited-state aromatization dynamics in the photochemical ring opening of dihydroazulene (DHA) is investigated by nonadiabatic molecular dynamics simulations in connection with the mixed-reference spin-flip (MRSF)-TDDFT method. It is found that, in the main reaction channel, the ring opening occurs in the excited state in a sequence of steps with increasing aromaticity. The first stage lasting ca. 200 fs produces an 8π semiaromatic S 1 minimum (S 1, min ) through an ultrafast damped bond length alternation (BLA) movement synchronized with a partial planarization of the cycloheptatriene ring. An additional ca. 200 fs are required to gain the vibrational energy needed to overcome a ring-opening transition state characterized by an enhanced Baird aromaticity. Unlike other BLA motions of ππ* state, it was shown that their damping is a characteristic feature of aromatic bond-equalization process. In addition, some minor channels of the reaction have also been discovered, where noticeably higher barriers of the S 1 non/antiaromatic transition structures must be surmounted. These anti-Baird channels led to reformation of DHA or other closed-ring products. The observed competition between the Baird and anti-Baird channels suggests that the quantum yield of photochemical products can be controllable by tipping their balance. Hence, here we suggest including the concept of anti-Baird, which would expand the applicability of Baird rule to much broader situations.

With the help of quantum mechanical methods, the formation of H2SO4 by the oxidation of H2SO3 with H2O2 was studied theoretically. Both stepwise and concerted mechanisms were calculated. It was found that the direct oxidation of H2SO3 by H2O2 alone requires prohibitive activation energies of >38.6 kcal/mol. However, the addition of one water molecule exhibits a strong catalytic effect that dramatically reduces the overall reaction barrier to 6.2 (2.3 with PCM) kcal/mol. The deprotonated HSO3 – species also reduces the overall reaction barrier to 5.6 (−5.8 with PCM) kcal/mol. Both of these proceed via concerted pathways. On the other hand, the stepwise mechanisms generally produce intermediates with a hydroperoxy group (−O–O–H), which is a result of a nucleophilic attack by the oxygens of H2O2. While studying the catalytic effect of water, a previously unknown hydroperoxy intermediate (HOO)­S­(OH)3, where sulfur is coordinated with three OH groups, was found. This work also reveals a rearrangement step of another hydroperoxy intermediate (HOO)­SO2 – to HSO4 – that was found in earlier experimental studies. For all of the mechanisms calculated, the final H2SO4 is formed with a significant exothermicity of >60 kcal/mol. In general, even without sunlight, it was found that the formation of sulfuric acid by hydrogen peroxide can occur in a heterogeneous moisturized environment.