Polyvalent s-block elements: A missing link challenges the periodic law of chemistry for the heavy elements.

The Periodic Law of Chemistry is one of the great discoveries in cultural history. Elements behaving chemically similar are empirically merged in groups G of a Periodic Table, each element with G valence electrons per neutral atom, and with upper limit G for the oxidation and valence numbers. Here, we report that among the usually mono- or di-valent s-block elements (G = 1 or 2), the heaviest members (87Fr, 88Ra, 119E, and 120E) with atomic numbers Z = 87, 88, 119, 120 form unusual 5- or 6-valent compounds at ambient conditions. Together with well-reported basic changes of valence at the end of the 6d-series, in the whole 7p-series, and for 5g6f-elements, it indicates that at the bottom of common Periodic Tables, the classic Periodic Law is not as straightforward as commonly expected. Specifically, we predict the feasible experimental synthesis of polyvalent [RaL-n] (n = 4, 6) compounds.

Fermionic quantum processing with programmable neutral atom arrays.

Simulating the properties of many-body fermionic systems is an outstanding computational challenge relevant to material science, quantum chemistry, and particle physics.-5.4pc]Please note that the spelling of the following author names in the manuscript differs from the spelling provided in the article metadata: D. González-Cuadra, D. Bluvstein, M. Kalinowski, R. Kaubruegger, N. Maskara, P. Naldesi, T. V. Zache, A. M. Kaufman, M. D. Lukin, H. Pichler, B. Vermersch, Jun Ye, and P. Zoller. The spelling provided in the manuscript has been retained; please confirm. Although qubit-based quantum computers can potentially tackle this problem more efficiently than classical devices, encoding nonlocal fermionic statistics introduces an overhead in the required resources, limiting their applicability on near-term architectures. In this work, we present a fermionic quantum processor, where fermionic models are locally encoded in a fermionic register and simulated in a hardware-efficient manner using fermionic gates. We consider in particular fermionic atoms in programmable tweezer arrays and develop different protocols to implement nonlocal gates, guaranteeing Fermi statistics at the hardware level. We use this gate set, together with Rydberg-mediated interaction gates, to find efficient circuit decompositions for digital and variational quantum simulation algorithms, illustrated here for molecular energy estimation. Finally, we consider a combined fermion-qubit architecture, where both the motional and internal degrees of freedom of the atoms are harnessed to efficiently implement quantum phase estimation as well as to simulate lattice gauge theory dynamics.

Informing geometric deep learning with electronic interactions to accelerate quantum chemistry.

Predicting electronic energies, densities, and related chemical properties can facilitate the discovery of novel catalysts, medicines, and battery materials. However, existing machine learning techniques are challenged by the scarcity of training data when exploring unknown chemical spaces. We overcome this barrier by systematically incorporating knowledge of molecular electronic structure into deep learning. By developing a physics-inspired equivariant neural network, we introduce a method to learn molecular representations based on the electronic interactions among atomic orbitals. Our method, OrbNet-Equi, leverages efficient tight-binding simulations and learned mappings to recover high-fidelity physical quantities. OrbNet-Equi accurately models a wide spectrum of target properties while being several orders of magnitude faster than density functional theory. Despite only using training samples collected from readily available small-molecule libraries, OrbNet-Equi outperforms traditional semiempirical and machine learning-based methods on comprehensive downstream benchmarks that encompass diverse main-group chemical processes. Our method also describes interactions in challenging charge-transfer complexes and open-shell systems. We anticipate that the strategy presented here will help to expand opportunities for studies in chemistry and materials science, where the acquisition of experimental or reference training data is costly.

Reliably assessing the electronic structure of cytochrome P450 on today's classical computers and tomorrow's quantum computers.

An accurate assessment of how quantum computers can be used for chemical simulation, especially their potential computational advantages, provides important context on how to deploy these future devices. To perform this assessment reliably, quantum resource estimates must be coupled with classical computations attempting to answer relevant chemical questions and to define the classical algorithms simulation frontier. Herein, we explore the quantum computation and classical computation resources required to assess the electronic structure of cytochrome P450 enzymes (CYPs) and thus define a classical-quantum advantage boundary. This is accomplished by analyzing the convergence of density matrix renormalization group plus n-electron valence state perturbation theory (DMRG+NEVPT2) and coupled-cluster singles doubles with noniterative triples [CCSD(T)] calculations for spin gaps in models of the CYP catalytic cycle that indicate multireference character. The quantum resources required to perform phase estimation using qubitized quantum walks are calculated for the same systems. Compilation into the surface code provides runtime estimates to compare directly to DMRG runtimes and to evaluate potential quantum advantage. Both classical and quantum resource estimates suggest that simulation of CYP models at scales large enough to balance dynamic and multiconfigurational electron correlation has the potential to be a quantum advantage problem and emphasizes the important interplay between classical computations and quantum algorithms development for chemical simulation.

Deep learning of dynamically responsive chemical Hamiltonians with semiempirical quantum mechanics.

Conventional machine-learning (ML) models in computational chemistry learn to directly predict molecular properties using quantum chemistry only for reference data. While these heuristic ML methods show quantum-level accuracy with speeds several orders of magnitude faster than traditional quantum chemistry methods, they suffer from poor extensibility and transferability; i.e., their accuracy degrades on large or new chemical systems. Incorporating quantum chemistry frameworks into the ML models directly solves this problem. Here we take the structure of semiempirical quantum mechanics (SEQM) methods to construct dynamically responsive Hamiltonians. SEQM methods use empirical parameters fitted to experimental properties to construct reduced-order Hamiltonians, facilitating much faster calculations than ab initio methods but with compromised accuracy. By replacing these static parameters with machine-learned dynamic values inferred from the local environment, we greatly improve the accuracy of the SEQM methods. Trained on molecular energies and atomic forces, these dynamically generated Hamiltonian parameters show a strong correlation with atomic hybridization and bonding. Trained with only about 60,000 small organic molecular conformers, the resulting model retains interpretability, extensibility, and transferability when testing on much larger chemical systems and predicting various molecular properties. Overall, this work demonstrates the virtues of incorporating physics-based descriptions with ML to develop models that are simultaneously accurate, transferable, and interpretable.

Exploring the Tl 2 H 2 potential energy surface: A comparative analysis with group 13 systems and experiment.

Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet Tl 2 H 2 potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug-cc-pwCVQZ-PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero-point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 M 2 H 2 (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist.

Learning organo-transition metal catalyzed reactions by graph neural networks.

Chemical reaction outcome prediction presents a fundamental challenge in synthetic chemistry. Most existing machine learning (ML) approaches focus on chemical reactions of typical elements. We developed a simple ML model focused on organo-transition metal-catalyzed reactions (OMCRs). Instead of overall reactions observed in experiments, we let the ML model learn the sequence of simplified elementary reactions. This drastically reduced the complexity of the model and helped it find common patterns from distinct reactions. We let a graph neural network learn the reactivity index of a pair of atoms. The model was able to learn a wide variety of OMCRs, and the accuracy of reaction prediction reached 97%, even though the model has extremely fewer learnable parameters than other standards. The learned reactivity indices of bonds nicely summarize the knowledge of reactions in the dataset.

Mutual influence of non-covalent interactions formed by imidazole: A systematic quantum-chemical study.

Imidazole is a five-membered heterocycle that is part of a number of biologically important molecules such as the amino acid histidine and the hormone histamine. Imidazole has a unique ability to participate in a variety of non-covalent interactions involving the NH group, the pyridine-like nitrogen atom or the π-system. For many biologically active compounds containing the imidazole moiety, its participation in formation of hydrogen bond NH⋯O/N and following proton transfer is the key step of mechanism of their action. In this work a systematic study of the mutual influence of various paired combinations of non-covalent interactions (e.g., hydrogen bonds and π-interactions) involving the imidazole moiety was performed by means of quantum chemistry (PW6B95-GD3/def2-QZVPD) for a series of model systems constructed based on analysis of available x-ray data. It is shown that for considered complexes formation of additional non-covalent interactions can only enhance the proton-donating ability of imidazole. At the same time, its proton-accepting ability can be both enhanced and weakened, depending on what additional interactions are added to a given system. The mutual influence of non-covalent interactions involving imidazole can be classified as weak geometric and strong energetic cooperativity-a small change in the length of non-covalent interaction formed by imidazole can strongly influence its strength. The latter can be used to develop methods for controlling the rate and selectivity of chemical reactions involving the imidazole fragment in larger systems. It is shown that the strong mutual influence of non-covalent interactions involving imidazole is due to the unique ability of the imidazole ring to effectively redistribute electron density in non-covalently bound systems with its participation.

A multifunctional flavoprotein monooxygenase HspB for hydroxylation and C-C cleavage of 6-hydroxy-3-succinoyl-pyridine.

Flavoprotein monooxygenases catalyze reactions, including hydroxylation and epoxidation, involved in the catabolism, detoxification, and biosynthesis of natural substrates and industrial contaminants. Among them, the 6-hydroxy-3-succinoyl-pyridine (HSP) monooxygenase (HspB) from Pseudomonas putida S16 facilitates the hydroxylation and C-C bond cleavage of the pyridine ring in nicotine. However, the mechanism for biodegradation remains elusive. Here, we refined the crystal structure of HspB and elucidated the detailed mechanism behind the oxidative hydroxylation and C-C cleavage processes. Leveraging structural information about domains for binding the cofactor flavin adenine dinucleotide (FAD) and HSP substrate, we used molecular dynamics simulations and quantum/molecular mechanics calculations to demonstrate that the transfer of an oxygen atom from the reactive FAD peroxide species (C4a-hydroperoxyflavin) to the C3 atom in the HSP substrate constitutes a rate-limiting step, with a calculated reaction barrier of about 20 kcal/mol. Subsequently, the hydrogen atom was rebounded to the FAD cofactor, forming C4a-hydroxyflavin. The residue Cys218 then catalyzed the subsequent hydrolytic process of C-C cleavage. Our findings contribute to a deeper understanding of the versatile functions of flavoproteins in the natural transformation of pyridine and HspB in nicotine degradation.IMPORTANCEPseudomonas putida S16 plays a pivotal role in degrading nicotine, a toxic pyridine derivative that poses significant environmental challenges. This study highlights a key enzyme, HspB (6-hydroxy-3-succinoyl-pyridine monooxygenase), in breaking down nicotine through the pyrrolidine pathway. Utilizing dioxygen and a flavin adenine dinucleotide cofactor, HspB hydroxylates and cleaves the substrate's side chain. Structural analysis of the refined HspB crystal structure, combined with state-of-the-art computations, reveals its distinctive mechanism. The crucial function of Cys218 was never discovered in its homologous enzymes. Our findings not only deepen our understanding of bacterial nicotine degradation but also open avenues for applications in both environmental cleanup and pharmaceutical development.

Phosphonium-substituted Diphosphaindenylide (PPI): Exploration of Biradical Character and Ligand Properties.

Starting from C6 H4 (PCl2 )2 and the TMS-substituted ylide (TMS)2 C=PR3 (TMS=trimethylsilyl, R=p-tolyl), the phosphonium-substituted diphosphaindenylide PPI was prepared in two steps. CASSCF calculations as well as the reactivity toward diphenyl acetylene suggest a notable biradical character in PPI. Reaction with [Cr(CO)3 (MeCN)3 ] affords the complex [Cr(CO)3 (η5 -PPI)] (5). This complex was employed to explore the ligand properties of PPI, which demonstrates considerable potential through the combination of strong metal-ligand interactions and the possibility of a pronounced indenyl effect.

Mono-Ortho-Beryllated Carbodiphosphoranes: Synthesis, Structure, Bonding and Reactivity.

The reaction of organoberyllium compounds with hexaphenylcarbodiphosphorane yields mono-ortho-beryllated complexes, which feature a double dative Be=C bond. The bonding situation in these compounds together with a simple carbodiphosphorane and an N-heterocyclic carbene adduct was analysed with energy decomposition analysis in combination with natural orbital for chemical valence as well as with quantum theory of atoms-in-molecules. Furthermore, the driving forces accountable for mono-ortho-beryllation were elucidated along with the reactivity of the Be=C bond.

Striking Borane Planarization in the Thermal Rearrangement (η5-C5H5)Fe(η3-B5H10)→(η5-C5H5)Fe(η5-B5H10).

In 1977 Weiss and Grimes, by means of mass spectrometry and 1H and 11B NMR spectroscopy, proposed two structures (I and II) for the ferraborane (η5-C5H5)Fe(B5H10), isoelectronic with ferrocene. In this work, by means of high-level quantum-chemical computations, we confirm the experimental structures of the two isomers with their corresponding energies, and assign the reported 1H and 11B NMR chemical shifts. A striking result from this study is the planarization (3D→2D) of the B5H10 - ligand - an unknown isolated anion, isoelectronic with aromatic cyclopentadienyl anion C5H5 - - when attached to the (η5-C5H5)Fe+ moiety, thus resulting in a more stable ferraborane isomer II.

Stability of the Protactinium(V) Mono-Oxo Cation Probed by First-Principle Calculations.

This study explores the distinctive behavior of protactinium (Z=91) within the actinide series. In contrast to neighboring elements like uranium or plutonium, protactinium in the pentavalent state diverges by not forming the typical dioxo protactinyl moiety PaO2 + in aqueous phase. Instead, it manifests as a monooxo PaO3+ cation or a Pa5+ . Employing first-principle calculations with implicit and explicit solvation, we investigate two stoichiometrically equivalent neutral complexes: PaO(OH)2 (X)(H2 O) and Pa(OH)4 (X), where X represents various monodentate and bidentate ligands. Calculating the Gibbs free energy for the reaction PaO(OH)2 (X)(H2 O)→Pa(OH)4 (X), we find that the PaO(OH)2 (X)(H2 O) complex is stabilized with Cl- , Br- , I- , NCS- , NO3 - , and SO4 2- ligands, while it is not favored with OH- , F- , and C2 O4 2- ligands. Quantum Theory of Atoms in Molecules (QTAIM) and Natural Bond Orbital (NBO) methods reveal the Pa mono-oxo bond as a triple bond, with significant contributions from the 5f and 6d shells. Covalency of the Pa mono-oxo bond increases with certain ligands, such as Cl- , Br- , I- , NCS- , and NO3 - . These findings elucidate protactinium's unique chemical attributes and provide insights into the conditions supporting the stability of relevant complexes.

Structure and Solution Behavior of Rare-Earth-Metal Complexes with Tripodal N-Donor Ligands.

Rare-earth-metal complexes (Ln=Y, La, Ce, Sm and Lu) of tripodal N-donor ligands respecting the CHON principle have been synthetized and characterized. The selectivity of the ligands through the lanthanide cations was investigated and related to their donor strength.

Investigation of Isolated IrF5 -, IrF6 - Anions and M[IrF6] (M=Na, K, Rb, Cs) Ion Pairs by Matrix-Isolation Spectroscopy and Relativistic Quantum-Chemical Calculations.

The molecular IrF5 -, IrF6 - anions and M[IrF6] (M=Na, K, Rb, Cs) ion pairs were prepared by co-deposition of laser-ablated alkali metal fluorides MF with IrF6 and isolated in solid neon or argon matrices under cryogenic conditions. The free anions were obtained as well by co-deposition of IrF6 with laser-ablated metals (Ir or Pt) as electron sources. The products were characterized in a combined analysis of matrix IR spectroscopy and electronic structure calculations using two-component quasi-relativistic DFT methods accounting for spin-orbit coupling (SOC) effects as well as multi-reference configuration-interaction (MRCI) approaches with SOC. Inclusion of SOC is crucial in the prediction of spectra and properties of IrF6 - and its alkali-metal ion pairs. The observed IR bands and the computations show that the IrF6 - anion adopts an Oh structure in a nondegenerate ground state stabilized by SOC effects, and not a distorted D4h structure in a triplet ground state as suggested by scalar-relativistic calculations. The corresponding "closed-shell" M[IrF6] ion pairs with C3v symmetry are stabilized by coordination of an alkali metal ion to three F atoms, and their structural change in the series from M=Na to Cs was proven spectroscopically. There is no evidence for the formation of IrF7, IrF7 - or M[IrF7] (M=Na, K, Rb, Cs) ion pairs in our experiments.

Oxidation-state sensitive light-induced dynamics of Ruthenium-4H-Imidazole complexes.

Oxidized molecular states are key intermediates in photo-induced redox reactions, e. g., intermolecular charge transfer between photosensitizer and catalyst in photoredox catalysis. The stability and longevity of the oxidized photosensitizer is an important factor in optimizing the respective light-driven reaction pathways. In this work the oxidized states of ruthenium(II)-4H-imidazole dyes are studied. The ruthenium complexes constitute benchmark photosensitizers in solar energy interconversion processes with exceptional chemical stability, strong visible light absorption, and favourable redox properties. To rationalize the light-induced reaction in the oxidized ruthenium(III) systems, we combine UV-vis absorption, resonance Raman, and transient absorption spectroelectrochemistry (SEC) with time-dependent density functional theory (TDDFT) calculations. Three complexes are compared, which vary with respect to their coordination environment, i. e., combining an 4H-imidazole with either 2,2'-bipyridine (bpy) or 2,2';6'2"-terpyridine (tpy) coligands, and chloride or isothiocyanate ligands. While all oxidized complexes have similar steady state absorption properties, their excited state kinetics differ significantly; the study thus opens the doorway to study the light-driven reactivity of oxidized molecular intermediates in intermolecular charge transfer cascades.

Probing the Isolobal Relation between Cp'''NiP3 and White Phosphorus by Experimental Charge Density Analysis.

An in-depth analysis of the description of bonding within Cp'''Ni-cyclo-P3 (Cp'''=1,2,4-tri-tert-butylcyclopentadienyl, [Ni]P3) employing X-ray diffraction based multipolar modeling, density functional theory (DFT) as well as an "experimental wavefunction" obtained from X-ray restrained wavefunction (XRW) fitting is presented. The results are compared to DFT calculations on white phosphorus - an isolobal analogue to [Ni]P3. A complementary bonding analysis shows insights into the reactivity of [Ni]P3. The isolobal principle is reflected in every aspect of our analysis and the employed methods seamlessly predict the differences in reactivity of [Ni]P3 and P4. Crystallographic modeling, solid-state NMR, and DFT calculations describe the dynamic behavior of the cyclo-P3 unit in the title molecule.

Exploring the Influence of the Phosphorus-Heteroatom Substitution in Nicotine on Its Electronic and Vibrational Spectroscopic Properties.

The influence of phosphorus substitution of nitrogen in heterocyclic compounds on the vibrational spectroscopy as well as frontier molecular orbitals are analyzed. Nicotine with two nitrogen atoms in its structure is taken as the sample system to be studied computationally. By replacing the nitrogen atom in one or both rings of this molecule with phosphorus, three nicotine derivatives are created. The vibrational circular dichroism and infrared spectra of these four molecules in their monomer state, as well as the assemblies up to trimers are determined. The aforementioned spectra are calculated using static quantum chemical calculations employing a cluster-weighted approach. The calculated gas phase spectra of nicotine are compared to their respective experimental spectra. It is observed that the nicotine derivatives with phosphorus in the methylpyrrolidine ring have considerably different gas phase and bulk phase vibrational circular dichroism spectra when compared to nicotine. The phosphorus substitution reduces the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital as well as altering the polarizability and reactivity of the investigated molecules.

Interactive Quantum Chemistry Enabled by Machine Learning, Graphical Processing Units, and Cloud Computing.

Modern quantum chemistry algorithms are increasingly able to accurately predict molecular properties that are useful for chemists in research and education. Despite this progress, performing such calculations is currently unattainable to the wider chemistry community, as they often require domain expertise, computer programming skills, and powerful computer hardware. In this review, we outline methods to eliminate these barriers using cutting-edge technologies. We discuss the ingredients needed to create accessible platforms that can compute quantum chemistry properties in real time, including graphical processing units-accelerated quantum chemistry in the cloud, artificial intelligence-driven natural molecule input methods, and extended reality visualization. We end by highlighting a series of exciting applications that assemble these components to create uniquely interactive platforms for computing and visualizing spectra, 3D structures, molecular orbitals, and many other chemical properties.

Vibrational Relaxation of D2O Induced by Collision with He: A Rigid Bender Close Coupling Study.

The vibrational relaxation of the first excited bending state of D2O induced by collision with He is studied at the close coupling level and using the Rigid Bender approximation. A new 4D potential energy surface is calculated and reported for this system. It is then used to determine the low-lying bound states of the D2O-He van der Waals complex and to perform scattering calculations. Collision rates are determined for pure rotational transitions as well as for rovibrational transitions within the first excited bending state. The results are compared with those obtained for the collision of D2O with other noble gases such as Ne and Ar. We also analyse the differences observed with respect to the H2O+He collisions and compare our results with experiment.