Implementation of Nonadiabatic Molecular Dynamics for Intersystem Crossing Based on a Time-Dependent Density-Functional Tight-Binding Method.

Intersystem crossing (ISC) and internal conversion (IC) are types of nonadiabatic transitions that play important roles in a wide range of fields, including photochemistry, photophysics, and photobiology. The nonadiabatic molecular dynamics (NA-MD) method is a powerful tool for computational simulations of dynamic phenomena involving nonadiabatic transitions. In this study, we implemented the NA-MD method, which treats ISC and IC on an equal footing, where the electronic structure is treated at the level of the time-dependent (TD) density-functional tight-binding (DFTB) method, a low-cost semiempirical analog of TD density functional theory (DFT). In particular, the spin-orbit coupling calculation algorithm was implemented in the TD-DFTB framework, and the results showed trends similar to those obtained using TD-DFT. In addition, the NA-MD method successfully reproduced ultrafast ISC of 2-nitronaphthalene.

Excitation configuration analysis for divide-and-conquer excited-state calculation method using dynamical polarizability.

The authors previously developed a divide-and-conquer (DC)-based non-local excited-state calculation method for large systems using dynamical polarizability [Nakai and Yoshikawa, J. Chem. Phys. 146, 124123 (2017)]. This method evaluates the excitation energies and oscillator strengths using information on the dynamical polarizability poles. This article proposes a novel analysis of the previously developed method to obtain further configuration information on excited states, including excitation and de-excitation coefficients of each excitation configuration. Numerical applications to simple molecules, such as ethylene, hydrogen molecule, ammonia, and pyridazine, confirmed that the proposed analysis could accurately reproduce the excitation and de-excitation coefficients. The combination with the DC scheme enables both the local and non-local excited states of large systems with an excited nature to be treated.

Quantum mechanical assessment on the optical properties of capsanthin conformers.

As optical properties, the ultraviolet-visible (UV-Vis) absorption spectra of capsanthin-based red natural dye are a decisive parameter for their usage in various applications. Thus, accurately predicting the maximum UV-Vis wavelength ( λ max ) values is critical in designing dye-conjugated material. Extensive metadynamics simulations were carried out to generate capsanthin conformers at various levels of the extended tight-binding method. Benchmarking the time-dependent density-functional theory (TD-DFT) methods help understand the results of a particular functional and allows a comparison between results obtained with different functional. The long-range correction (LC) scheme in LC-TD-DFT-D4/ωB97X/def2-SVP has been found to reproduce the experimental λ max , and exhibited the effect of conformational changes to the calculated wavelengths. On the other hand, an inexpensive yet efficient LC-TD-DFTB method reproduced the experimental λ max insensitive to conformational changes.

Divide-and-Conquer Linear-Scaling Quantum Chemical Computations.

Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.

Species-selective nanoreactor molecular dynamics simulations based on linear-scaling tight-binding quantum chemical calculations.

Here, extensions to quantum chemical nanoreactor molecular dynamics simulations for discovering complex reactive events are presented. The species-selective algorithm, where the nanoreactor effectively works for the selected desired reactants, was introduced to the original scheme. Moreover, for efficient simulations of large model systems with the modified approach, the divide-and-conquer linear-scaling density functional tight-binding method was exploited. Two illustrative applications of the polymerization of propylene and cyclopropane mixtures and the aggregation of sodium chloride from aqueous solutions indicate that species-selective quantum chemical nanoreactor molecular dynamics is a promising method to accelerate the sampling of multicomponent chemical processes proceeding under relatively mild conditions.

Enabling large-scale quantum path integral molecular dynamics simulations through the integration of Dcdftbmd and i-PI codes.

A large-scale quantum chemical calculation program, Dcdftbmd, was integrated with a Python-based advanced atomistic simulation program, i-PI. The implementation of a client-server model enabled hierarchical parallelization with respect to replicas and force evaluations. The established framework demonstrated that quantum path integral molecular dynamics simulations can be executed with high efficiency for systems consisting of a few tens of replicas and containing thousands of atoms. The application of the framework to bulk water systems, with and without an excess proton, demonstrated that nuclear quantum effects are significant for intra- and inter-molecular structural properties, including oxygen-hydrogen bond distance and radial distribution function around the hydrated excess proton.

Neutral-to-ionic photoinduced phase transition of tetrathiafulvalene-p-chloranil by electronic and vibrational excitation: A real-time nuclear-electronic dynamics simulation study.

Tetrathiafulvalene-p-chloranil exhibits photoinduced phase transition (PIPT) between neutral (N) and ionic (I) phases, in which the constituent molecules are approximately charge-neutral and ionic, respectively. In addition to visible-light irradiation, which can induce both N → I and I → N PIPTs, infrared irradiation has been reported to induce the N → I PIPT. These results suggest that N → I and I → N PIPTs can be driven by electronic excitation, and the I → N PIPT can also be driven by vibrational excitation. However, the feasibility of the N → I PIPT using vibrational excitation remains an open question. In this study, we address this issue by simulating the PIPT processes using a nonadiabatic molecular dynamics approach combined with real-time electron dynamics at the level of a semiempirical quantum chemical model, density-functional tight binding. The results show the importance of vibronic interactions in the PIPT processes, thereby suggesting the possibility of N → I PIPT by vibrational excitations with infrared irradiation.

Unveiling controlling factors of the S0/S1 minimum-energy conical intersection (3): Frozen orbital analysis based on the spin-flip theory.

Conical intersections (CIs), which indicate the crossing of two or more adiabatic electronic states, are crucial in the mechanisms of photophysical, photochemical, and photobiological processes. Although various geometries and energy levels have been reported using quantum chemical calculations, the systematic interpretation of the minimum energy CI (MECI) geometries is unclear. A previous study [Nakai et al., J. Phys. Chem. A 122, 8905 (2018)] performed frozen orbital analysis (FZOA) based on time-dependent density functional theory (TDDFT) at the MECI formed between the ground and first electronic excited states (S0/S1 MECI), thereby inductively clarifying two controlling factors. However, one of the factors that the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) energy gap became close to the HOMO-LUMO Coulomb integral was not valid in the case of spin-flip TDDFT (SF-TDDFT), which is frequently used as a means of the geometry optimization of MECI [Inamori et al., J. Chem. Phys. 152, 144108 (2020)]. This study revisited the controlling factors using FZOA for the SF-TDDFT method. Based on spin-adopted configurations within a minimum active space, the S0-S1 excitation energy is approximately represented by the HOMO and LUMO energy gap ΔεHL, a contribution from Coulomb integrals JHL″ and that from the HOMO-LUMO exchange integral KHL″. Furthermore, numerical applications of the revised formula at the SF-TDDFT method confirmed the control factors of S0/S1 MECI.

Direct Near Infrared Light-Activatable Phthalocyanine Catalysts.

The high penetration of near-infrared (NIR) light makes it effective for use in selective reactions under light-shielded conditions, such as in sealed reactors and deep tissues. Herein, we report the development of phthalocyanine catalysts directly activated by NIR light to transform small organic molecules. The desired photocatalytic properties were achieved in the phthalocyanines by introducing the appropriate peripheral substituents and central metal. These phthalocyanine photocatalysts promote cross-dehydrogenative-coupling (CDC) under irradiation with 810 nm NIR light. The choice of solvent is important, and a mixture of a reaction-accelerating (pyridine) and -decelerating (methanol) solvents was particularly effective. Moreover, we demonstrate photoreactions under visible-light-shielded conditions through the transmission of NIR light. A combined experimental and computational mechanistic analysis revealed that this NIR reaction does not involve a photoredox-type mechanism with electron transfer, but instead a singlet-oxygen-mediated mechanism with energy transfer.

Hydroxide Ion Mechanism for Long-Range Proton Pumping in the Third Proton Transfer of Bacteriorhodopsin.

In bacteriorhodopsin, representative light-driven proton pump, five proton transfers yield vectorial active proton translocation, resulting in a proton gradient in microbes. Third proton transfer occurs from Asp96 to the Schiff base on the photocycle, which is expected to be a long-range proton transfer via the Grotthuss mechanism through internal water molecules. Here, large-scale quantum molecular dynamics simulations are performed for the third proton transfer, where all the atoms (∼50000 atoms) are treated quantum-mechanically. The simulations demonstrate that two reaction paths exist along the water wire, namely, via hydronium and via hydroxide ions. The free energy analysis confirms that the path via hydroxide ions is considerably favorable and consistent with the observed lifetime of the transient water wire. Therefore, the proposed hydroxide ion mechanism, as in the first proton transfer, is responsible for the third long-range proton transfer.

Experimental and Theoretical Evidence for Relativistic Catalytic Activity in C-H Activation of N-Phenylbenzamide Using a Cationic Iridium Complex.

This study elucidates that relativistic effect plays a key role in catalytic C-H activation using a cationic Ir complex. Experiments show that the cationic Ir(I)-diphosphine catalyst can be used for the deuterium substitution of N-phenylbenzamide, whereas a cationic Rh(I)-diphosphine catalyst is scarcely effective. Density functional theory calculations, including the relativistic effect, demonstrate a large difference in the reaction energy diagrams for the C-H activation of N-phenylbenzamide between the cationic Ir and Rh catalysts. In particular, the relatively low reaction barrier and considerably stabilized product obtained for the Ir catalysts are rationalized by strong Ir-C and Ir-H interactions, which originate from the relativistic self-consistent d-orbital expansion of Ir.

Multiple protonation states in ligand-free SARS-CoV-2 main protease revealed by large-scale quantum molecular dynamics simulations.

The main protease (Mpro) in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) catalyzes the cleavage of polyproteins for viral replication. Here, large-scale quantum molecular dynamics and metadynamics simulations for ligand-free Mpro were performed, where all the atoms were treated quantum-mechanically, focusing on elucidation of the controversial active-site protonation state. The simulations clarified that the interconverting multiple protonation states exist in unliganded Mpro, and the catalytically relevant ion-pair state is more stable than the neutral state, which is consistent with neutron crystallography. The results highlight the importance of the ion-pair state for repurposing or discovering antiviral drugs that target Mpro.

Tetraaryldiborane(4) Can Emit Dual Fluorescence Responding to the Structural Change around the B-B Bond.

We report the successful synthesis of tetramesityldiborane(4) (Mes4 B2 ) through the reductive coupling of a dimesitylborinium ion. Owing to the steric protection conferred by the mesityl groups, Mes4 B2 shows exceptional chemical stability and remains intact in water. Single-crystal X-ray analysis revealed that Mes4 B2 has an orthogonal geometry, where the B-B center is completely hidden by the mesityl groups. Remarkably, Mes4 B2 emits dual fluorescence at 460 and 620 nm, both in solution and in the solid state. Theoretical calculations showed that Mes4 B2 in the excited S1 state adopts a twisted or planar geometry, which is responsible for the shorter- or longer-wavelength fluorescence, respectively. The intensity ratio of the dual fluorescence is sensitive to the viscosity of the medium, which suggests that Mes4 B2 has potential as a ratiometric viscosity sensor.

Picture-change correction in relativistic density functional theory.

Relativistic quantum chemical calculations are performed based on one of two physical pictures, namely the Dirac picture and the Schrödinger picture. With regard to the latter, the so-called picture-change effect (PCE) and picture-change correction (PCC) have been studied. The PCE, which is the change in the expectation value associated with the transformation, is not commonly a minor effect. The electron density, which is given by the expectation value of the density operator, is a fundamental variable in relativistic density functional theory (RDFT). Thus, performing the PCC in RDFT calculations is essential not only in terms of numerical agreement with the Dirac picture, but also from the viewpoint of fundamental theory. This paper explains theories and numerical studies of PCE and PCC in RDFT after overviewing those in properties, which involves the authors' works on the development of RDFT in the Schrödinger picture and relativistic exchange-correlation functionals based on picture-change-corrected variables.

Theoretical prediction by DFT and experimental observation of heterocation-doping effects on hydrogen adsorption and migration over the CeO2(111) surface.

Hydrogen (H) atom adsorption and migration over the CeO2-based materials surface are of great importance because of its wide applications to catalytic reactions and electrochemical devices. Therefore, comprehensive knowledge for controlling the H atom adsorption and migration over CeO2-based materials is crucially important. For controlling H atom adsorption and migration, we investigated irreducible divalent, trivalent, and quadrivalent heterocation-doping effects on H atom adsorption and migration over the CeO2(111) surface using density functional theory (DFT) calculations. Results revealed that the electron-deficient lattice oxygen (Olat) and the flexible CeO2 matrix played key roles in strong adsorption of H atoms. Heterocations with smaller valence and smaller ionic radius induced the electron-deficient Olat. In addition, smaller cation doping enhanced the CeO2 matrix flexibility. Moreover, we confirmed the influence of H atom adsorption controlled by doping on surface proton migration (i.e. surface protonics) and catalytic reaction involving surface protonics (NH3 synthesis in an electric field). Results confirmed clear correlation between H atom adsorption energy and surface protonics.

Assessing locally range-separated hybrid functionals from a gradient expansion of the exchange energy density.

Locally range-separated hybrid (LRSH) functionals feature a real-space-dependent range separation function (RSF) instead of a system-independent range-separation parameter, which thus enables a more flexible admixture of exact exchange than conventional range-separated hybrid functionals. In particular, the development of suitable RSF models and exploring the capabilities of the LRSH approach, in general, are tasks that require further investigations and will be addressed in this work. We propose a non-empirical scheme based on a detailed scaling analysis with respect to a uniform coordinate scaling and on a short-range expansion of the range-separated exchange energy density to derive new RSF models from a gradient expansion of the exchange energy density. After optimizing a small set of empirical parameters introduced to enhance their flexibility, the resulting second- and fourth-order RSFs are evaluated with respect to atomic exchange energies, atomization energies, and transition barrier heights.

An Air- and Water-Stable B4 N4 -Heteropentalene Serving as a Host Material for a Phosphorescent OLED.

Replacement of the carbon-carbon bonds of antiaromatic compounds with polar boron-nitrogen bonds often provides isoelectronic BN compounds with excellent thermodynamic stability and interesting photophysical properties. By this element-substitution strategy, we synthesized a new B4 N4 -heteropentalene derivative, 1, which is fully substituted with mesityl groups. Owing to kinetic protection by the sterically bulky substituents, 1 is remarkably stable toward air and even water. Single-crystal X-ray analysis of 1 revealed the bonding characteristics of the B4 N4 -heteropentalene structure. In a glassy matrix, 1 emitted short-wavelength phosphorescence with an onset at 350 nm, indicating that the triplet energy is substantially high. DFT calculations reasonably explained the ground- and excited-state electronic structures of 1 as well as its emission properties. Motivated by the high-energy triplet state of 1, we used it as a host material to fabricate a phosphorescent organic light-emitting diode with an external quantum efficiency of 15 %.

An Element-Substituted Cyclobutadiene Exhibiting High-Energy Blue Phosphorescence.

1,3,2,4-Diazadiboretidine, an isoelectronic heteroanalogue of cyclobutadiene, is an interesting chemical species in terms of comparison with the carbon system, whereas its properties have never been investigated experimentally. According to Baird's rule, Hückel antiaromatic cyclobutadiene acquires aromaticity in the lowest triplet state. Here we report experimental and theoretical studies on the ground- and excited-state antiaromaticity/aromaticity as well as the photophysical properties of an isolable 1,3,2,4-diazadiboretidine derivative. The crystal structure of the diazadiboretidine derivative revealed that the B2 N2 ring adopts a planar rhombic geometry in the ground state. Yet, theoretical calculations showed that the B2 N2 ring turns to a square geometry with a nonaromatic character in the lowest triplet state. Notably, the diazadiboretidine derivative has the lowest singlet and triplet states lying at close energy levels and displays blue phosphorescence.

Hierarchical parallelization of divide-and-conquer density functional tight-binding molecular dynamics and metadynamics simulations.

Massively parallel divide-and-conquer density functional tight-binding (DC-DFTB) molecular dynamics and metadynamics simulations are efficient approaches for describing various chemical reactions and dynamic processes of large complex systems via quantum mechanics. In this study, DC-DFTB simulations were combined with multi-replica techniques. Specifically, multiple walkers metadynamics, replica exchange molecular dynamics, and parallel tempering metadynamics methods were implemented hierarchically into the in-house Dcdftbmd program. Test simulations in an aqueous phase of the internal rotation of formamide and conformational changes of dialanine showed that the newly developed extensions increase the sampling efficiency and the exploration capabilities in DC-DFTB configuration space.

Spin-flip approach within time-dependent density functional tight-binding method: Theory and applications.

A spin-flip time-dependent density functional tight-binding (SF-TDDFTB) method is developed that describes target states as spin-flipping excitation from a high-spin reference state obtained by the spin-restricted open shell treatment. Furthermore, the SF-TDDFTB formulation is extended to long-range correction (LC), denoted as SF-TDLCDFTB. The LC technique corrects the overdelocalization of electron density in systems such as charge-transfer systems, which is typically found in conventional DFTB calculations as well as density functional theory calculations using pure functionals. The numerical assessment of the SF-TDDFTB method shows smooth potential curves for the bond dissociation of hydrogen fluoride and the double-bond rotation of ethylene and the double-cone shape of H3 as the simplest degenerate systems. In addition, numerical assessments of SF-TDDFTB and SF-TDLCDFTB for 39 S0 /S1 minimum energy conical intersection (MECI) structures are performed. The SF-TDDFTB and SF-TDLCDFTB methods drastically reduce the computational cost with accuracy for MECI structures compared with SF-TDDFT.