Mapping the emergence of molecular vibrations mediating bond formation.

Fundamental studies of chemical reactions often consider the molecular dynamics along a reaction coordinate using a calculated or suggested potential energy surface1-5. But fully mapping such dynamics experimentally, by following all nuclear motions in a time-resolved manner-that is, the motions of wavepackets-is challenging and has not yet been realized even for the simple stereotypical bimolecular reaction6-8: A-B + C â†’ A + B-C. Here we track the trajectories of these vibrational wavepackets during photoinduced bond formation of the gold trimer complex [Au(CN)2-]3 in an aqueous monomer solution, using femtosecond X-ray liquidography9-12 with X-ray free-electron lasers13,14. In the complex, which forms when three monomers A, B and C cluster together through non-covalent interactions15,16, the distance between A and B is shorter than that between B and C. Tracking the wavepacket in three-dimensional nuclear coordinates reveals that within the first 60 femtoseconds after photoexcitation, a covalent bond forms between A and B to give A-B + C. The second covalent bond, between B and C, subsequently forms within 360 femtoseconds to give a linear and covalently bonded trimer complex A-B-C. The trimer exhibits harmonic vibrations that we map and unambiguously assign to specific normal modes using only the experimental data. In principle, more intense X-rays could visualize the motion not only of highly scattering atoms such as gold but also of lighter atoms such as carbon and nitrogen, which will open the door to the direct tracking of the atomic motions involved in many chemical reactions.

Phosphorus-Ligand Redox Cooperative Catalysis: Unraveling Four-Electron Dioxygen Reduction Pathways and Reactive Intermediates.

The reduction of dioxygen to water is crucial in biology and energy technologies, but it is challenging due to the inertness of triplet oxygen and complex mechanisms. Nature leverages high-spin transition metal complexes for this, whereas main-group compounds with their singlet state and limited redox capabilities exhibit subdued reactivity. We present a novel phosphorus complex capable of four-electron dioxygen reduction, facilitated by unique phosphorus-ligand redox cooperativity. Spectroscopic and computational investigations attribute this cooperative reactivity to the unique electronic structure arising from the geometry of the phosphorus complex bestowed by the ligand. Mechanistic study via spectroscopic and kinetic experiments revealed the involvement of elusive phosphorus intermediates resembling those in metalloenzymes. Our result highlights the multielectron reactivity of phosphorus compound emerging from a carefully designed ligand platform with redox cooperativity. We anticipate that the work described expands the strategies in developing main-group catalytic reactions, especially in small molecule fixations demanding multielectron redox processes.

Quantum chemical calculations of electron affinities of alkaline earth metal atoms (Ca, Sr, Ba, and Ra).

We performed high-level ab initio quantum chemical calculations, incorporating higher-order excitations, spin-orbit coupling (SOC), and the Gaunt interaction, to calculate the electron affinities (EAs) of alkaline earth (AE) metal atoms (Ca, Sr, Ba, and Ra), which are notably small. The coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method is insufficient to accurately calculate the EAs of AE metal atoms. Higher-order excitations proved crucial, with the coupled-cluster singles, doubles, and triples with perturbative quadruples [CCSDT(2)Q] method effectively capturing dynamic electron correlation effects. The contributions of SOC (ΔESOs) to the EAs calculated using the multireference configuration interaction method with the Davidson correction, including SOC, positively enhance the EAs; however, these contributions are overestimated. The Dirac-Hartree-Fock (DHF)-CCSD(T) method addresses this overestimation and provides reasonable values for ΔESO (ΔESO-D). Employing additional sets of diffuse and core-valence correlation basis sets is critical for accurately calculating the EAs of AE metal atoms. The contributions of the Gaunt interaction (ΔEGaunt) to the EAs of AE metal atoms are negligible. Notably, the CCSDT(2)Q with the complete basis set limit + ΔESO-D + ΔEGaunt produced EA values for Ca, Sr, and Ba that closely aligned with experimental data and achieved accuracy exceeding the chemical accuracy. Based on our findings, the accurately proposed EA for Ra is 9.88 kJ/mol.

Photochemical Reaction Mechanism of Intramolecular H Transfer Reaction of H2 C3 O+ ⋠Radical Cation.

The photochemical reaction mechanism underlying the intramolecular H-transfer of the H2 C3 O+ ⋅ radical cation to the H2 CCCO+ ⋅ methylene ketene cation was elucidated using time-dependent density functional theory and high-level ab initio methods. Once the D1 state of H2 C3 O+ ⋅ is populated, the reaction proceeds to form an intermediate (IM) in the D1 state (IM4D1 ). The molecular structure of the conical intersection (CI) was optimized using a multiconfigurational ab initio method. The CI is readily accessible because it lies slightly above the IM4D1 in energy. In addition, the gradient difference vector of the CI is almost parallel to the intramolecular H-transfer reaction coordinate. Once the vibration mode of IM4D1 which is parallel to the reaction coordinate is populated, the degeneracy of the CI is readily lifted and H2 CCCO+ ⋅ was formed via a relaxation pathway in the D0 state. Our calculated results clearly describe the photochemical intramolecular H transfer reaction reported in a recent study.

Ultrafast photoisomerization mechanism of azaborine revealed by nonadiabatic molecular dynamics simulations.

1,2-Dihydro-1,2-azaborine is an isoelectronic analog of benzene with a B-N substitution, and its unique photoisomerization behavior, which is distinct from that of benzene, has drawn significant attention. To understand the detailed mechanism of azaborine photochemistry considering the dynamical effect and gain a comprehensive understanding of photochemical reactions, we investigated the photoisomerization dynamics of azaborine using nonadiabatic molecular dynamics simulations with Tully's surface hopping algorithm. Herein, the structural and energetic analyses of the trajectories revealed three different paths: direct relaxation (path 1), relaxation via a prefulvene-like intermediate (path 2), and formation of the Dewar isomer as a photoproduct (path 3). Our results confirmed that the photoisomerization of azaborine follows the energetically favored pathway predicted by the previous minimum energy path (MEP) calculations, exclusively forming the Dewar isomer, which is consistent with the experimental observations. Additionally, despite the low quantum yield found in our simulations, the high-level excitation energy calculations support the complete conversion observed in the experiments.

Do Planar Tetracoordinate Fluorine Atoms Exist? Revisiting a Theoretical Prediction.

We re-examined the existence of planar tetracoordinate F (ptF) atoms, which was proposed recently by using high-level ab initio methods such as coupled-cluster singles and doubles with perturbative triples (CCSD(T)) with large basis sets. Our calculations indicate that the planar structures of FIn4+ (D4h), FTl4+ (D4h), FGaIn3+ (C2V), FIn2Tl2+ (D2h), FIn3Tl+ (C2V), and FInTl3+ (C2V) are not the minimum energy states; by contrast, they are transition states. Density functional theory calculations overestimate the size of the cavity formed by the four peripheral atoms, leading to erroneous conclusions regarding the existence of ptF atoms. Our analysis suggests that the preference for non-planar structures in the six cations studied is not due to the pseudo Jahn-Teller effect. Additionally, spin-orbit coupling does not alter the main conclusion that the ptF atom does not exist. If sufficiently large cavity formation by group 13 elements to accommodate the central F- ion is guaranteed, then the existence of ptF atoms is plausible.

High level ab initio and density functional study of TeF6+ and TeCl6+: Attainability of +7 oxidation state for tellurium.

Discovery of a new oxidation state for an element expands its chemistry. A high oxidation state, such as +7, is rare for sp-block elements except for halogens. In this study, we determined that Te can attain a +7 oxidation state through the existence of a distorted octahedron (DOH) structure of TeCl6+ based on coupled cluster singles and doubles with perturbative triples calculations. We propose a new type of isomerization that resembles pseudorotation. The octahedron structure of TeF6+ bearing one elongated axial bond isomerizes to a DOH via an associated pseudorotation.

Increasing Chemical Diversity of B2 N2 Anthracene Derivatives by Introducing Continuous Multiple Boron-Nitrogen Units.

Increasing the chemical diversity of organic semiconductors is essential to develop efficient electronic devices. In particular, the replacement of carbon-carbon (C-C) bonds with isoelectronic boron-nitrogen (B-N) bonds allows precise modulation of the electronic properties of semiconductors without significant structural changes. Although some researchers have reported the preparation of B2 N2 anthracene derivatives with two B-N bonds, no compounds with continuous multiple BN units have been prepared yet. Herein, we report the synthesis and characterization of a B2 N2 anthracene derivative with a BNBN unit formed by converting the BOBN unit at the zigzag edge. Compared to the all-carbon analogue 2-phenylanthracene, BNBN anthracene exhibits significant variations in the C-C bond length and a larger highest occupied molecular orbital-lowest unoccupied molecular orbital energy gap. The experimentally determined bond lengths and electronic properties of BNBN anthracene are confirmed through theoretical calculations. The BOBN anthracene organic light-emitting diode, used as a blue host, exhibits a low driving voltage. The findings of this study may facilitate the development of larger acenes with multiple BN units and potential applications in organic electronics.

Optical Kerr Effect of Liquid Acetonitrile Probed by Femtosecond Time-Resolved X-ray Liquidography.

Optical Kerr effect (OKE) spectroscopy is a method that measures the time-dependent change of the birefringence induced by an optical laser pulse using another optical laser pulse and has been used often to study the ultrafast dynamics of molecular liquids. Here we demonstrate an alternative method, femtosecond time-resolved X-ray liquidography (fs-TRXL), where the microscopic structural motions related to the OKE response can be monitored using a different type of probe, i.e., X-ray solution scattering. By applying fs-TRXL to acetonitrile and a dye solution in acetonitrile, we demonstrate that different types of molecular motions around photoaligned molecules can be resolved selectively, even without any theoretical modeling, based on the anisotropy of two-dimensional scattering patterns and extra structural information contained in the q-space scattering data. Specifically, the dynamics of reorientational (libration and orientational diffusion) and translational (interaction-induced motion) motions are captured separately by anisotropic and isotropic scattering signals, respectively. Furthermore, the two different types of reorientational motions are distinguished from each other by their own characteristic scattering patterns and time scales. The measured time-resolved scattering signals are in excellent agreement with the simulated scattering signals based on a molecular dynamics simulation for plausible molecular configurations, providing the detailed structural description of the OKE response in liquid acetonitrile.

Direct observation of bond formation in solution with femtosecond X-ray scattering.

The making and breaking of atomic bonds are essential processes in chemical reactions. Although the ultrafast dynamics of bond breaking have been studied intensively using time-resolved techniques, it is very difficult to study the structural dynamics of bond making, mainly because of its bimolecular nature. It is especially difficult to initiate and follow diffusion-limited bond formation in solution with ultrahigh time resolution. Here we use femtosecond time-resolved X-ray solution scattering to visualize the formation of a gold trimer complex, [Au(CN)2(-)]3 in real time without the limitation imposed by slow diffusion. This photoexcited gold trimer, which has weakly bound gold atoms in the ground state, undergoes a sequence of structural changes, and our experiments probe the dynamics of individual reaction steps, including covalent bond formation, the bent-to-linear transition, bond contraction and tetramer formation with a time resolution of ∼500 femtoseconds. We also determined the three-dimensional structures of reaction intermediates with sub-ångström spatial resolution. This work demonstrates that it is possible to track in detail and in real time the structural changes that occur during a chemical reaction in solution using X-ray free-electron lasers and advanced analysis of time-resolved solution scattering data.

Theoretical Investigation of Excited-State Intramolecular Proton Transfer and Photoisomerization of 2-(Iminomethyl)phenol.

Time-dependent density functional theory and high-level ab initio calculations were performed to investigate the excited-state intramolecular proton transfer (ESIPT) and subsequent isomerization of 2-(iminomethyl)phenol (IMP). According to the results of the correlated theoretical methods, ESIPT is a barrierless process; subsequently, the isomerization (rotation of the torsion angle) of IMP also readily occurs. Fictitious intermediates are found due to an insufficient theoretical level. The molecular structure of the conical intersection (CI) during the isomerization process is optimized, and its branching plane is characterized. Both the gradient difference vector and the derivative coupling vector are significantly correlated to the C═O and H2N-C antiparallel stretching coordinates, and the dynamic electron correlation effect is crucial to optimize the molecular structure of the real CI of IMP. The relaxation pathway from the CI in the S0 state was examined; the dominant pathway proceeds to the trans-keto form of IMP. However, if the C═O and H2N-C antiparallel stretching mode is sufficiently populated, then the reaction proceeds to the cis-keto form of IMP and can eventually recover to the cis-enol form.

Electronic and Structural Dynamics of Dicyanoaurate Trimer in Excited State.

A trimer of dicyanoaurate has been studied as a model system of the covalent chemical bond formation. Here, we report the dynamics of dicyanoaurate trimer in water upon photoexcitation by femtosecond time-resolved luminescence (TL) and luminescence spectra at cyrogenic temperature. Temperature was varied as a means to control the medium flexibility as well as the population of isomers. A unique parallelism between the luminescence spectrum vs. time and vs. temperature was observed, which enables unambiguous luminescence band assignments and facilitates investigation of the dynamics. Upon photoexcitation to S1, intersystem crossing proceeds in an ultrafast manner within 20 fs due to the large spin-orbit coupling followed by a structural change from a loose bent to a tight linear form in 1.5 ps. Higher oligomerization occurs above the melting temperature. TL reveals a strong coherent excitation of the symmetric Au-Au stretching vibration at 74 cm-1 through the non-Condon effect.

Fate of transient isomer of CH2I2: Mechanism and origin of ionic photoproducts formation unveiled by time-resolved x-ray liquidography.

Diiodomethane, CH2I2, in a polar solvent undergoes a unique photoinduced reaction whereby I2 - and I3 - are produced from its photodissociation, unlike for other iodine-containing haloalkanes. While previous studies proposed that homolysis, heterolysis, or solvolysis of iso-CH2I-I, which is a major intermediate of the photodissociation, can account for the formation of I2 - and I3 -, there has been no consensus on its mechanism and no clue for the reason why those negative ionic species are not observed in the photodissociation of other iodine-containing chemicals in the same polar solvent, for example, CHI3, C2H4I2, C2F4I2, I3 -, and I2. Here, using time-resolved X-ray liquidography, we revisit the photodissociation mechanism of CH2I2 in methanol and determine the structures of all transient species and photoproducts involved in its photodissociation and reveal that I2 - and I3 - are formed via heterolysis of iso-CH2I-I in the photodissociation of CH2I2 in methanol. In addition, we demonstrate that the high polarity of iso-CH2I-I is responsible for the unique photochemistry of CH2I2.

Mechanistic Investigation of Thermal and Photoreactions between Boron and Silane.

Density functional theory and high-level ab initio calculations were performed to elucidate the detailed reaction mechanism from B and SiH4 to a structure with two bridging H atoms (Si(µ-H2)BH2, silicon tetrahydroborate). On the basis of the calculated results, this reaction mechanism includes both thermal and photochemical reactions. Especially, thermal conversion of silylene dihydroborate (H2B═SiH2) to Si(µ-H2)BH2 is not feasible because two high energetic barriers must be overcome. In contrast, the reverse reaction is feasible because it is effectively only necessary to overcome a single barrier. The characteristics of the excited states of H2B═SiH2 and Si(µ-H2)BH2 have been identified. Two successive conical intersections (CIs) are involved in the photochemical reaction. The BSiH4 bending coordinate is almost parallel to the reaction coordinate near the regions from the second CI to Si(µ-H2)BH2. The activated BSiH4 bending mode lift the degeneracy of the second CI, thereby the reaction readily proceeds to Si(µ-H2)BH2. All calculated results in this work reasonably well describe the recent experimental observations.

Benchmark Study of Density Functional Theory for Neutral Gold Clusters, Aun (n = 2-8).

Neutral gold clusters, Aun (n = 2-8), were optimized using coupled cluster singles and doubles with perturbative triples (CCSD(T)) with a triple-ζ-level basis set to develop reliable reference values for their structural and energy parameters in order to assess the performance of density functionals. The performance of 44 density functional theory (DFT) methods for calculating molecular structures and relative energies is assessed with respect to CCSD(T). In addition, their performance when calculating vertical ionization potentials (vIPs) of Aun (n = 2-8) is also assessed by comparison with experimental values. The revTPSS functional shows good performance for calculating both the structural and energy properties of Aun (n = 2-8), whereas B3P86 shows a remarkable performance in calculating the vIPs. The quadruple-ζ-level valence basis set is necessary for obtaining accurate energy values in CCSD(T) calculations.

Spin-Orbit Effect on the Molecular Properties of TeXn (X = F, Cl, Br, and I; n = 1, 2, and 4): A Density Functional Theory and Ab Initio Study.

Density functional theory (DFT) and ab initio calculations, including spin-orbit coupling (SOC), were performed to investigate the spin-orbit (SO) effect on the molecular properties of tellurium halides, TeXn (X = F, Cl, Br, and I; n = 1, 2, and 4). SOC elongates the Te-X bond and slightly reduces the vibrational frequencies. Consideration of SOC leads to better agreement with experimental values. Møller-Plesset second-order perturbation theory (MP2) seriously underestimates the Te-X bond lengths. In contrast, B3LYP significantly overestimates them. SO-PBE0 and multireference configuration interactions with the Davidson correction (MRCI+Q), which include SOC via a state-interaction approach, give the Te-I bond length of TeI2 that matches the experimental value. On the basis of the calculated thermochemical energy and optimized molecular structure, TeI4 is unlikely to be stable. The use of PBE0 including SOC is strongly recommended for predicting the molecular properties of Te-containing compounds.

Ab Initio Investigation of the Ground States of F2P(S)N, F2PNS, and F2PSN.

A recent spectroscopic experiment identified difluorothiophosphoryl nitrene (F2P(S)N) and found that it showed rich photochemistry. However, a discrepancy between the experimental results and the quantum chemical calculations was reported. Thus, high-level ab initio calculations using the coupled cluster singles and doubles with perturbative triples and second-order multiconfigurational perturbation theory were performed to elucidate this inconsistency. The discrepancy arose due to the failure to consider the triplet state of difluoro(thionitroso)phosphine (F2PNS). In this work, we identify that the global minimum of the system is the triplet state of F2PNS, which allows us to explain the inconsistency between the experimental and theoretical results. All calculated results give consistent results with the recent experimental results.

Performance of Density Functional Theory and Relativistic Effective Core Potential for Ru-Based Organometallic Complexes.

Herein a performance assessment of density functionals used for calculating the structural and energetic parameters of bi- and trimetallic Ru-containing organometallic complexes has been performed. The performance of four popular relativistic effective core potentials (RECPs) has also been assessed. On the basis of the calculated results, the MN12-SX (range-separated hybrid functional) demonstrates good performance for calculating the molecular structures, while MN12-L (local functional) performs well for calculating the energetics, including that of the Ru-Ru bond breaking process. The choice of appropriate density functional is a crucial factor for calculating the energetics. The LANL08 demonstrates the lowest performance of the RECPs for calculating the molecular structures, especially the Ru-Ru bond length.

Theoretical Investigation of the Reaction Mechanism of the Photoisomerization of 1,2-Dihydro-1,2-azaborine.

The photoisomerization of 1,2-dihydro-1,2-azaborine was investigated by high-level multireference ab initio and density functional theory calculations. The intermediates (IMs) and transition states (TSs) on the S(0) and S(1) states were optimized using the state-averaged complete active space self-consistent field method. The multireference configuration interaction method with the Davidson correction was used to obtain accurate energetics. Moreover, the conical intersections (CIs), which play a crucial role in photoisomerization, were also optimized. On the basis of the calculation results, the most favorable proposed reaction pathway is as follows: reactant→Franck-Condon region→TS(1) →CI→IM(0) →TS(0P) →product. The product was not directly formed through the CI, and the IM(0) existed on the S(0) state. These results show that the isomerization of 1,2-dihydro-1,2-azaborine involves both photoreactions and thermal reactions. The calculated results clarify recent experimental observations.