Intense Circularly Polarized Luminescence Induced by Chiral Supramolecular Assembly: The Importance of Intermolecular Electronic Coupling.

Herein we report on circularly polarized luminescence (CPL) emission originating from supramolecular chirality of organic microcrystals with a |glum| value up to 0.11. The microcrystals were prepared from highly emissive difluoroboron ß-diketonate (BF2dbk) dyes R-1 or S-1 with chiral binaphthol (BINOL) skeletons. R-1 and S-1 exhibit undetectable CPL signals in solution but manifest intense CPL emission in their chiral microcrystals. The chiral superstructures induced by BINOL skeletons were confirmed by XRD analysis. Spectral analysis and theoretical calculations indicate that intermolecular electronic coupling, mediated by the asymmetric stacking in the chiral superstructures, effectively alters excited-state electronic structures and facilitates electron transitions perpendicular to BF2bdk planes. The coupling increases cosθµ,m from 0.05 (monomer) to 0.86 (tetramer) and triggers intense optical activity of BF2bdk. The results demonstrate that optical activity of chromophores within assemblies can be regulated by both orientation and extent of intermolecular electronic couplings.

Construction of Highly Accurate Machine Learning Potential Energy Surfaces for Excited-State Dynamics Simulations Based on Low-Level Data Sets.

Machine learning is capable of effectively predicting the potential energies of molecules in the presence of high-quality data sets. Its application in the construction of ground- and excited-state potential energy surfaces is attractive to accelerate nonadiabatic molecular dynamics simulations of photochemical reactions. Because of the huge computational cost of excited-state electronic structure calculations, the construction of a high-quality data set becomes a bottleneck. In the present work, we first built two data sets. One was obtained from surface hopping dynamics simulations at the semiempirical OM2/MRCI level. Another was extracted from the dynamics trajectories at the CASSCF level, which was reported previously. The ground- and excited-state potential energy surfaces of ethylene-bridged azobenzene at the CASSCF computational level were constructed based on the former low-level data set. Although non-neural network machine learning methods can achieve good or modest performance during the training process, only neural network models provide reliable predictions on the latter external test data set. The BPNN and SchNet combined with the Δ-ML scheme and the force term in the loss functions are recommended for dynamics simulations. Then, we performed excited-state dynamics simulations of the photoisomerization of ethylene-bridged azobenzene on machine learning potential energy surfaces. Compared with the lifetimes of the first excited state (S1) estimated at different computational levels, our results on the E isomer are in good agreement with the high-level estimation. However, the overestimation of the Z isomer is unimproved. It suggests that smaller errors during the training process do not necessarily translate to more accurate predictions on high-level potential energies or better performance on nonadiabatic dynamics simulations, at least in the present case.

Quantum Computation of Conical Intersections on a Programmable Superconducting Quantum Processor.

Conical intersections (CIs) are pivotal in many photochemical processes. Traditional quantum chemistry methods, such as the state-average multiconfigurational methods, face computational hurdles in solving the electronic Schrödinger equation within the active space on classical computers. While quantum computing offers a potential solution, its feasibility in studying CIs, particularly on real quantum hardware, remains largely unexplored. Here, we present the first successful realization of a hybrid quantum-classical state-average complete active space self-consistent field method based on the variational quantum eigensolver (VQE-SA-CASSCF) on a superconducting quantum processor. This approach is applied to investigate CIs in two prototypical systems─ethylene (C2H4) and triatomic hydrogen (H3). We illustrate that VQE-SA-CASSCF, coupled with ongoing hardware and algorithmic enhancements, can lead to a correct description of CIs on existing quantum devices. These results lay the groundwork for exploring the potential of quantum computing to study CIs in more complex systems in the future.

Nonadiabatic Molecular Dynamics in Momentum Space Beyond Harmonic Approximation: Hot Electron Relaxation in Photoexcited Black Phosphorus.

We simulated hot-electron relaxation in black phosphorus using the nonadiabatic molecular dynamics (NA-MD) approach with a non-Condon effect in momentum space beyond the harmonic approximation. By comparing simulations at the Γ point in a large supercell with those using a few k-points in a smaller supercell─while maintaining the same number of electronic states within the same energy range, we demonstrate that both setups yield remarkably consistent energy relaxation times, regardless of the initial state energy. This consistency arises from the complementary effects of supercell size in real space and the number of k-points in the reciprocal space. This finding confirms that simulations at a single k-point in large size supercells are an effective approximation for NA-MD with a non-Condon effect. This approach offers significant advantages for complex photophysics, such as intervalley scattering and indirect bandgap charge recombination, and is particularly suitable for large systems without the need for a harmonic approximation.

Mechanistic Insights into Chloric Acid Production by Hydrolysis of Chlorine Trioxide at an Air-Water Interface.

Chlorine oxides play crucial roles in ozone depletion, and the final oxidation steps of chlorine oxide potentially result in the formation of chloric acid (HClO3) or perchloric acid (HClO4). Herein, the solvation and reactive uptake of three stable isomers of chlorine trioxide (Cl2O3), namely, ClOCl(O)O, ClClO3, and ClOOOCl, at the air-water interface were investigated using classical and hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) coupled with advanced free energy methods. Two distinct mechanisms were revealed for the hydrolysis of ClOCl(O)O and ClClO3: molecular and ionic mechanisms. A comparison of the computed free-energy profiles for the gaseous and air-water interfacial systems indicated that the air-water interface could markedly lower the free-energy barrier for ClO3- or HClO3 formation while stabilizing the product state. In particular, the hydrolysis of ClClO3 at the air-water interface was barrierless. In contrast, our calculations showed that the hydrolysis of ClOOOCl was very slow, indicating that ClOOOCl was inert to water at the air-water interface. This study provides theoretical evidence for the hypothesis that HClO3 is a sink for chlorine oxides and for the widespread distributions of HClO3 recently observed in the Arctic region.

Multi-fidelity machine learning for predicting bandgaps of nonlinear optical crystals.

Nonlinear optical (NLO) materials are of great importance in modern optics and industry because of their intrinsic capability of wavelength conversion. Bandgap is a key property of NLO crystals. In recent years, machine learning (ML) has become a powerful tool to predict the bandgaps of compounds before synthesis. However, the shortage of available experimental data of NLO crystals poses a significant challenge for the exploration of new NLO materials using ML. In this work, we proposed a new multi-fidelity ML approach based on the multilevel descriptors developed by us (Z.-Y. Zhang, X. Liu, L. Shen, L. Chen and W.-H. Fang, J. Phys. Chem. C, 2021, 125, 25175-25188) and the gradient boosting regression tree algorithm. The calculated and experimental bandgaps of NLO crystals were collected as the low- and high-fidelity labels, respectively. The experimental values were predicted based on chemical compositions of crystals without prior knowledge about crystal structures. The multi-fidelity ML model overcame the performance of single-fidelity predictor. Furthermore, it was observed that less accurate predictions on the low-fidelity label may result in more accurate prediction on the high-fidelity label, at least in the present case. Using the multi-fidelity ML model with the best performance in this work, the mean absolute error on the test set of experimental bandgaps was 0.293 eV, which is smaller than that using the single-fidelity model (0.355 eV). It is far from perfect but accurate enough as an effective computational tool in the first step to discover novel NLO materials.

Compression of Organic Molecules Coupled with Hydrogen Bonding Extends the Charge Carrier Lifetime in BA2SnI4.

Two-dimensional (2D) metal halide perovskites, such as BA2SnI4 (BA═CH3(CH2)3NH3), exhibit an enhanced charge carrier lifetime in experiments under strain. Experiments suggest that significant compression of the BA molecule, rather than of the inorganic lattice, contributes to this enhancement. To elucidate the underlying physical mechanism, we apply a moderate compressive strain to the entire system and subsequently introduce significant compression to the BA molecules. We then perform ab initio nonadiabatic molecular dynamics simulations of nonradiative electron-hole recombination. We observe that the overall lattice compression reduces atomic motions and decreases nonadiabatic coupling, thereby delaying electron-hole recombination. Additionally, compression of the BA molecules enhances hydrogen bonding between the BA molecules and iodine atoms, which lengthens the Sn-I bonds, distorts the [SnI6]4- octahedra, and suppresses atomic motions further, thus reducing nonadiabatic coupling. Also, the elongated Sn-I bonds and weakened antibonding interactions increase the band gap. Altogether, the compression delays the nonradiative electron-hole recombination by more than a factor of 3. Our simulations provide new and valuable physical insights into how compressive strain, accommodated primarily by the organic ligands, positively influences the optoelectronic properties of 2D layered halide perovskites, offering a promising pathway for further performance improvements.

Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation.

Herein, we have used the "on-the-fly" ring-polymer surface-hopping simulation method with the centroid approximation (RPSH-CA), in combination with the multireference OM2/MRCI electronic structure calculations to study the photoinduced dynamics of a green fluorescent protein (GFP) chromophore analogue in the gas phase, i.e., o-HBI, at 50, 100, and 300 K with 1, 5, 10, and 15 beads (3600 1 ps trajectories). The electronic structure calculations identified five new minimum-energy conical intersection (MECI) structures, which, together with the previous one, play crucial roles in the excited-state decay dynamics of o-HBI. It is also found that the excited-state intramolecular proton transfer (ESIPT) occurs in an ultrafast manner and is completed within 20 fs in all the simulation conditions because there is no barrier associated with this ESIPT process in the S1 state. However, the other excited-state dynamical results are strongly related to the number of beads. At 50 and 100 K, the nuclear quantum effects (NQEs) are very important; therefore, the excited-state dynamical results change significantly with the bead number. For example, the S1 decay time deduced from time-dependent state populations becomes longer as the bead number increases. Nevertheless, an essentially convergent trend is observed when the bead number is close to 10. In contrast, at 300 K, the NQEs become weaker and the above dynamical results converge very quickly even with 1 bead. Most importantly, the NQEs seriously affect the excited-state decay mechanism of o-HBI. At 50 and 100 K, most trajectories decay to the S0 state via perpendicular keto MECIs, whereas, at 300 K, only twisted keto MECIs are responsible for the excited-state decay. The present work not only comprehensively explores the temperature-dependent photoinduced dynamics of o-HBI, but also demonstrates the importance and necessity of NQEs in nonadiabatic dynamics simulations, especially at relatively low temperatures.

Photoinduced Dynamics of a Single-Walled Carbon Nanotube with a sp3 Defect: The Importance of Excitonic Effects.

Herein, we employed linear-response time-dependent functional theory nonadiabatic dynamic simulations to explore the photoinduced exciton dynamics of a chiral single-walled carbon nanotube CNT(6,5) covalently doped with a 4-nitrobenzyl group (CNT65-NO2). The results indicate that the introduction of a sp3 defect leads to the splitting of the degenerate VBM/VBM-1 and CBM/CBM+1 states. Both the VBM upshift and the CBM downshift are responsible for the experimentally observed redshifted E11* trapping state. The simulations reveal that the photoinduced exciton relaxation dynamics completes within 500 fs, which is consistent with the experimental work. On the other hand, we also conducted the nonadiabatic carrier (electron and hole) dynamic simulations, which completely ignore the excitonic effects. The comparison demonstrates that excitonic effects are indispensable. Deep analyses show that such effects induce several dark states, which play an important role in regulating the photoinduced dynamics of CNT65-NO2. The present work demonstrates the importance of including excitonic effects in simulating photoinduced processes of carbon nanotubes. In addition, it not only rationalizes previous experiments but also provides valuable insights that will help in the future rational design of novel covalently doped carbon nanotubes with superior photoluminescent properties.

Topological wetting states of microdroplets on closed-loop structured surfaces: Breakdown of the Gibbs equation at the microscale.

Microdroplets are a class of soft matter that has been extensively employed for chemical, biochemical, and industrial applications. However, fabricating microdroplets with largely controllable contact-area shape and apparent contact angle, a key prerequisite for their applications, is still a challenge. Here, by engineering a type of surface with homocentric closed-loop microwalls/microchannels, we can achieve facile size, shape, and contact-angle tunability of microdroplets on the textured surfaces by design. More importantly, this class of surface topologies (with universal genus value = 1) allows us to reveal that the conventional Gibbs equation (widely used for assessing the edge effect on the apparent contact angle of macrodroplets) seems no longer applicable for water microdroplets or nanodroplets (evidenced by independent molecular dynamics simulations). Notably, for the flat surface with the intrinsic contact angle ~0°, we find that the critical contact angle on the microtextured counterparts (at edge angle 90°) can be as large as >130°, rather than 90° according to the Gibbs equation. Experiments show that the breakdown of the Gibbs equation occurs for microdroplets of different types of liquids including alcohol and hydrocarbon oils. Overall, the microtextured surface design and topological wetting states not only offer opportunities for diverse applications of microdroplets such as controllable chemical reactions and low-cost circuit fabrications but also provide testbeds for advancing the fundamental surface science of wetting beyond the Gibbs equation.

Solid-State Photochemical Cascade Process Boosting Smart Ultralong Room-Temperature Phosphorescence in Bismuth Halides.

Molecular ultralong room-temperature phosphorescence (RTP), exhibiting multiple stimuli-responsive characteristics, has garnered considerable attention due to its potential applications in light-emitting devices, sensors, and information safety. This work proposes the utilization of photochemical cascade processes (PCCPs) in molecular crystals to design a stepwise smart RTP switch. By harnessing the sequential dynamics of photo-burst movement (induced by [2+2] photocycloaddition) and photochromism (induced by photogenerated radicals) in a bismuth (Bi)-based metal-organic halide (MOH), a continuous and photo-responsive ultralong RTP can be achieved. Furthermore, utilizing the same Bi-based MOH, diverse application demonstrations, such as multi-mode anti-counterfeiting and information encryption, can be easily implemented. This work thus not only serves as a proof-of-concept for the development of solid-state PCCPs that integrate photosalient effect and photochromism with light-chemical-mechanical energy conversion, but also lays the groundwork for designing new Bi-based MOHs with dynamically responsive ultralong RTP.

Significant Impact of Defect Fluctuation on Charge Dynamics in CsPbI3: A Study Combining Machine Learning with Quantum Dynamics.

In this study, we developed a machine-learned force field for CsPbI3 using a neural network potential, enabling molecular dynamics simulations (MD) with ab initio accuracy over nanoseconds. This approach, combined with ab initio MD and nonadiabatic MD, was used to study the charge trapping and recombination dynamics in both pristine and defective CsPbI3. Our simulations revealed key transitions affecting carrier lifetimes, especially in systems with iodine vacancy and interstitial iodine defects. An iodine trimer, formed when iodine replaces cesium, exhibits a high-frequency phonon mode. This mode enhances nonadiabatic coupling, accelerating charge recombination in defective systems compared to pristine ones. In the iodine vacancy system, recombination times varied significantly due to differences in NA coupling and energy gaps. The interplay between nonadiabatic coupling and pure dephasing time is crucial in determining recombination times for interstitial iodine defects. Our findings highlight the role of defect evolution in perovskites, offering insights for enhancing perovskite performance.

Unveiling Mechanistic Insights and Photocatalytic Advancements in Intramolecular Photo-(3 + 2)-Cycloaddition: A Comparative Assessment of Two Paradigmatic Single-Electron-Transfer Models.

The synthesis of 1-aminonorbornane (1-aminoNB), a potential aniline bioisostere, through photochemistry or photoredox catalysis signifies a remarkable breakthrough with implications in organic chemistry, pharmaceutical chemistry, and sustainable chemistry. However, an understanding of the underlying mechanisms involved in these reactions remains limited and ambiguous. Herein, we employ high-precision CASPT2//CASSCF calculations to elucidate the intricate mechanisms regulating the intramolecular photo-(3 + 2)-cycloaddition reactions for the synthesis of 1-aminoNB in the presence or absence of the Ir-complex-based photocatalyst. Our investigations delve into radical cascades, stereoselectivity, particularly single-electron-transfer (SET) events, etc. Furthermore, we innovatively introduce and compare two SET models integrating Marcus electron-transfer theory and transition-state theory. These models combined with kinetic data contribute to recognizing the critical control factors in diverse photocatalysis, thereby guiding the design and manipulation of photoredox catalysis as well as the improvement and modification of photocatalysts.

Dual Functionality of 6-Methylthioguanine: Synergistic Effects Enhancing the Photolability of DNA Nucleobases.

A small chemical modification of the nucleobase structure can significantly enhance the photoactivity of DNA, which may incur DNA damage, thus holding promising applications in photochemotherapy treatment of cancers or pathogens. However, single substitution confers only limited phototoxicity to DNA. Herein, we combine femtosecond and nanosecond time-resolved spectroscopy with high-level ab initio calculations to disentangle the excited-state dynamics of 6-methylthioguanine (me6-TG) under variable wavelength UVA excitation (310-330 nm). We find that double substitution of nucleobases (thionation and methylation) boosts the photoactivity by introducing more reactive channels. Intriguingly, 1nNπ*, rather than 1nSπ*, acts as the doorway state engendering the formation of the long-lived reactive triplet state in me6-TG. The 1nNπ* induces a low spin-orbit coupling of 8.3 cm-1, which increases the intersystem crossing (ISC) time (2.91 ± 0.14 ns). Despite the slowed ISC, the triplet quantum yield (ΦT) still accounts for a large fraction (0.6 ± 0.1), consistent with the potential energy surface that favors excited-state bifurcation to 1nNπ*min (3.36 ± 0.15 ps) rather than 1ππ*min (5.05 ± 0.26 ps), such that the subsequent ISC to triplet via 1nNπ*min constitutes the main relaxation pathway in me6-TG. Although this ΦT is inferior to its single-substituted predecessor 6-thioguanine (6-TG, 0.8 ± 0.2), the effect of thionation in synergy with methylation opens a unique C-S bond cleavage pathway through crossing to a repulsive 1πσ* state, generating thiyl radicals as highly reactive intermediates that may invoke biological damage. This photodissociation channel is extremely difficult for conventional nucleobases. These findings demonstrate the synergistic effects of double functionality substitution in modulating excited-state dynamics and enhancing the photolabile character of DNA nucleobases, providing inspirations for the rational design of advanced photodynamic and photochemotherapy approaches.

Mechanistic insight into the competition between interfacial and bulk reactions in microdroplets through N2O5 ammonolysis and hydrolysis.

Reactive uptake of dinitrogen pentaoxide (N2O5) into aqueous aerosols is a major loss channel for NOx in the troposphere; however, a quantitative understanding of the uptake mechanism is lacking. Herein, a computational chemistry strategy is developed employing high-level quantum chemical methods; the method offers detailed molecular insight into the hydrolysis and ammonolysis mechanisms of N2O5 in microdroplets. Specifically, our calculations estimate the bulk and interfacial hydrolysis rates to be (2.3 ± 1.6) × 10-3 and (6.3 ± 4.2) × 10-7 ns-1, respectively, and ammonolysis competes with hydrolysis at NH3 concentrations above 1.9 × 10-4 mol L-1. The slow interfacial hydrolysis rate suggests that interfacial processes have negligible effect on the hydrolysis of N2O5 in liquid water. In contrast, N2O5 ammonolysis in liquid water is dominated by interfacial processes due to the high interfacial ammonolysis rate. Our findings and strategy are applicable to high-chemical complexity microdroplets.

Integrating Full-Color 2D Optical Waveguide and Heterojunction Engineering in Halide Microsheets for Multichannel Photonic Logical Gates.

Ensuring information security has emerged as a paramount concern in contemporary human society. Substantial advancements in this regard can be achieved by leveraging photonic signals as the primary information carriers, utilizing photonic logical gates capable of wavelength tunability across various time and spatial domains. However, the challenge remains in the rational design of materials possessing space-time-color multiple-resolution capabilities. In this work, a facile approach is proposed for crafting metal-organic halides (MOHs) that offer space-time-color resolution. These MOHs integrate time-resolved room temperature phosphorescence and color-resolved excitation wavelength dependencies with both space-resolved ex situ optical waveguides and in situ heterojunctions. Capitalizing on these multifaceted properties, MOHs-based two-dimensional (2D) optical waveguides and heterojunctions exhibit the ability to tune full-color emissions across the spectra from blue to red, operating within different spatial and temporal scales. Therefore, this work introduces an effective methodology for engineering space-time-color resolved MOH microstructures, holding significant promise for the development of high-density photonic logical devices.

Physics-Constrained Hardware-Efficient Ansatz on Quantum Computers That Is Universal, Systematically Improvable, and Size-Consistent.

Variational wave function ansätze are at the heart of solving quantum many-body problems in physics and chemistry. Previous designs of hardware-efficient ansatz (HEA) on quantum computers are largely based on heuristics and lack rigorous theoretical foundations. In this work, we introduce a physics-constrained approach for designing HEA with rigorous theoretical guarantees by imposing a few fundamental constraints. Specifically, we require that the target HEA to be universal, systematically improvable, and size-consistent, which is an important concept in quantum many-body theories for scalability but has been overlooked in previous designs of HEA. We extend the notion of size-consistency to HEA and present a concrete realization of HEA that satisfies all these fundamental constraints while only requiring linear qubit connectivity. The developed physics-constrained HEA is superior to other heuristically designed HEA in terms of both accuracy and scalability, as demonstrated numerically for the Heisenberg model and some typical molecules. In particular, we find that restoring size-consistency can significantly reduce the number of layers needed to reach a certain accuracy. In contrast, the failure of other HEA to satisfy these constraints severely limits their scalability to larger systems with more than 10 qubits. Our work highlights the importance of incorporating physical constraints into the design of HEA for efficiently solving many-body problems on quantum computers.

A Complete Unveiling of the Mechanism and Chirality in Photoisomerization of Arylazopyrazole 3pzH: Combined Electronic Structure Calculations and AIMS Dynamic Simulations.

The arylazopyrazole 3pzH as a novel photoswitch exhibits quantitative switching and high thermal stability. In this work, combined electronic structure calculations and ab initio multiple spawning (AIMS) dynamic simulations were performed to systemically investigate the cis ↔ trans photoisomerization mechanism and the chiral preference after photoexcitation of 3pzH to the first excited singlet state (S1). Unlike most of the azoheteroarene photoswitches reported previously, many twisted and T-shaped cis isomers were found to be stable for 3pzH in the S0 state, owing to the moderate interaction between the hydrogen atom and π electrons of the aromatic ring. Two twisted cis isomers with different chirality ((M)-Z1 and (P)-Z1), the most stable T-shaped cis isomer ((T)-Z2), and the most stable planar trans isomer (E2) were selected as the initial structures to carry out the AIMS nonadiabatic dynamic simulations. Following excitation to the S1 state, all of the cis isomers decayed to conical intersection (CI) regions via the same bicycle pedal mechanism, while the evolution of the trans isomers to their CI regions was achieved via rotation around the N═N bond. More importantly, chiral preferences were found for the twisted cis isomers in the S1 state through the AIMS dynamic simulations due to the steric effect and static electronic repulsion. Notably, chirality was also observed in S1 isomerization starting from the planar E2 isomer because of the dynamic effect. After the nonadiabatic transition to the S0 state, the bicycle pedal mechanism was found to play a crucial role in cis ↔ trans photoisomerization. The simulated photoisomerization productivities were generally consistent with past experimental observations. Our calculations not only uncover the underlying reason for the excellent photoswitching properties of 3pzH but also enrich the knowledge of photoisomerization for azoheteroarene photoswitches, which will surely benefit their rational design.

Distributed Multi-GPU Ab Initio Density Matrix Renormalization Group Algorithm with Applications to the P-Cluster of Nitrogenase.

The presence of many degenerate d/f orbitals makes polynuclear transition-metal compounds, such as iron-sulfur clusters in nitrogenase, challenging for state-of-the-art quantum chemistry methods. To address this challenge, we present the first distributed multi-graphics processing unit (GPU) ab initio density matrix renormalization group (DMRG) algorithm suitable for modern high-performance computing (HPC) infrastructures. The central idea is to parallelize the most computationally intensive part─the multiplication of O(K2) operators with a trial wave function, where K is the number of spatial orbitals, by combining operator parallelism for distributing the workload with a batched algorithm for performing contractions on GPU. With this new implementation, we are able to reach an unprecedentedly large bond dimension D = 14,000 on 48 GPUs (NVIDIA A100 80 GB SXM) for an active space model (114 electrons in 73 active orbitals) of the P-cluster, which is nearly 3 times larger than the bond dimensions reported in previous DMRG calculations for the same system using only central processing units (CPUs).

Photocatalytic Reduction of CO2 to HCOOH and CO by a Phosphine-Bipyridine-Phosphine Ir(III) Catalyst: Photophysics, Nonadiabatic Effects, Mechanism, and Selectivity.

Photocatalytic CO2 reduction is one of the best solutions to solve the global energy crisis and to realize carbon neutralization. The tetradentate phosphine-bipyridine (bpy)-phosphine (PNNP)-type Ir(III) photocatalyst, Mes-IrPCY2, was reported with a high HCOOH selectivity but the photocatalytic mechanism remains elusive. Herein, we employ electronic structure methods in combination with radiative, nonradiative, and electron transfer rate calculations, to explore the entire photocatalytic cycle to either HCOOH or CO, based on which a new mechanistic scenario is proposed. The catalytic reduction reaction starts from the generation of the precursor metal-to-ligand charge transfer (3 MLCT) state. Subsequently, the divergence happens from the 3 MLCT state, the single electron transfer (SET) and deprotonation process lead to the formation of one-electron-reduced species and Ir(I) species, which initiate the reduction reaction to HCOOH and CO, respectively. Interestingly, the efficient occurrence of proton or electron transfer reduces barriers of critical steps. In addition, nonadiabatic transitions play a nonnegligible role in the cycle. We suggest a lower free-energy barrier in the reaction-limiting step and the very efficient SET in 3 MLCT are cooperatively responsible for a high HCOOH selectivity. The gained mechanistic insights could help chemists to understand, regulate, and design photocatalytic CO2 reduction reaction of similar function-integrated molecular photocatalyst.