The SHARK integral generation and digestion system

A perspective on the future of quantum chemical software: the example of the ORCA program package

The field of computational chemistry has made an impressive impact on contemporary chemical research. In order to carry out computational studies on actual systems, sophisticated software is required in form of large-scale quantum chemical program packages. Given the enormous diversity and complexity of the methods that need to be implementation in such packages, it is evident that these software pieces are very large (millions of code lines) and extremely complex. Most of the packages in widespread use by the computational chemistry community have had a development history of decades. Given the rapid progress in the hardware and a lack of resources (time, workforce, money), it is not possible to keep redesigning these program packages from scratch in order to keep up with the ever more quickly shifting hardware landscape. In this perspective, some aspects of the multitude of challenges that the developer community faces are discussed. While the task at hand - to ensure that quantum chemical program packages can keep evolving and make best use of the available hardware - is daunting, there are also new evolving opportunities. The problems and potential cures are discussed with the example of the ORCA package that has been developed in our research group.

Insights into the adsorption mechanism of chlorpyrifos on activated carbon derived from prickly pear seeds waste: An experimental and DFT modeling study

The removal of chlorpyrifos (CPF) from water was achieved using activated carbon (AC) derived from prickly pear seeds (PPS) wastes, developed through chemical activation with phosphoric acid. Several physico-chemical characterization methods were employed. The determination of surface functions using the Boehm assay indicated that the processed AC predominantly possesses acidic functions. The results obtained from the Boehm assay were corroborated by the pH value of the point of zero charge (pHpzc), which was equal to 2.5. Specific area calculation by the BET (Brunauer Emmett Teller) method revealed a large specific area (SBET) of 1077.66 m2 g⁻1. Adsorption experiments of CPF on AC demonstrated that the pseudo-second order (PSO) model and the Freundlich model were the most suitable for kinetic and isothermal modeling, respectively. The maximum CPF adsorption capacity of the PPS AC was found to be approximately 35 mg g⁻1. A theoretical study employing the density functional theory (DFT) was conducted using the B3LYP/6-311G (d, p) method. The most reliable adsorption energy (Eads) and Gibbs free energy (ΔGads) values between CPF and the functional groups on the AC surface were calculated. Results indicated a strong interaction between the lactone group of AC and CPF (ΔGads = -7.15 kcal mol⁻1, ΔEads = -21.55 kcal mol⁻1) and the hydroxyl group (ΔGads = -6.61 kcal mol⁻1, ΔEads = -20.66 kcal mol⁻1). This study demonstrates that activated carbon possesses significant adsorption power, making it highly effective for depolluting water contaminated by pesticides. The application of the theoretical DFT method enhances the understanding of the adsorption phenomenon of CPF on AC.

1.3 Å Crystal Structure of E. coli Peptidyl-Prolyl Isomerase B with Uniform Substitution of Valine by (2S,3S)-4-Fluorovaline Reveals Structure Conservation and Multiple Staggered Rotamers of CH2F Groups

(2S,3S)-4-Fluorovaline (FVal) is an analogue of valine, where a single CH3 group is substituted by a CH2F group. In the absence of valine, E. coli valyl-tRNA synthetase uses FVal as a substitute, enabling the production of proteins uniformly labeled with FVal. Here, we describe the production and analysis of E. coli peptidyl-prolyl isomerase B where all 16 valine residues have been replaced by FVal synthesized with a 13C-labeled CH2F group. Although the melting temperature is lower by about 11 °C relative to the wild-type protein, the three-dimensional protein structure is almost completely conserved, as shown by X-ray crystallography. The CH2F groups invariably populate staggered rotamers. Most CH2F groups populate two different rotamers. The increased space requirement of fluorine versus hydrogen does not prohibit rotamers that position fluorine next to a backbone carbonyl carbon. 19F NMR spectra show a signal dispersion over 25 ppm. The most high-field shifted 19F resonances correlate with large 3JHF coupling constants, confirming the impact of the γ-gauche effect on the signal dispersion. The present work is the second experimental verification of the effect and extends its validity to fluorovaline. The abundance of valine in proteins and structural conservation with FVal renders this valine analogue attractive for probing proteins by 19F NMR spectroscopy.

Assessment of the similarity-transformed equation of motion (STEOM) for open-shell organic and transition metal molecules

This study benchmarks the newly re-implemented single-reference excited-state methods, IP-EOM-CCSD, EA-EOM-CCSD, and STEOM-CCSD, in ORCA6.0, with a focus on open-shell systems. We compare STEOM against EOM-CCSD, CC3, and CCSDT across a range of systems, including small organic radicals, hydrated transition metal (TM) ions, and TM diatomic systems with both closed and open-shell configurations. For organic radicals, STEOM and EOM-CCSD show comparable performance, aligning closely with CC3 and CCSDT results. In the case of hydrated TM ions, IP-EOM closely matches DLPNO-CCSD results, while deviations from DLPNO-CCSD(T) are consistent. For open-shell TM systems, IP-EOM exhibits a blueshift relative to both the DLPNO-CCSD methods, while EA-EOM-CCSD shows better agreement. When comparing STEOM and CC3 to CCSDT, STEOM shows slightly larger deviations in closed-shell systems but shows excellent agreement in open-shell systems. Computational efficiency is also assessed, revealing a significant speedup in ORCA 6.0 compared to ORCA 5.0, with optimizations improving computation times. This study provides valuable insights into the performance and efficiency of STEOM in various chemical environments, highlighting its strengths and limitations.

Rapidly Gelled Lipoic Acid-Based Supramolecular Hydrogel for 3D Printing of Adhesive Bandage

Developing a strongly adhesive, easily removable, and robust bandage is valuable in trauma emergencies. Poly(lipoic acid) (PLA)-based adhesives with good mechanical properties have been well-developed through a thermal ring-opening polymerization (ROP) method that is easiness. However, the additive manufacturing of PLA-based adhesives remains a challenge. Herein, α-lipoic acid (LA) and trometamol (Tris) are found to rapidly form a supramolecular hydrogel at room temperature with injectability and 3D printing potential. Meanwhile, the synthesized LA-grafted hyaluronic acid and cellulose nanocrystals are involved not only to optimize the extrusion of 3D printing but also to effectively promote fidelity and prevent the inverse closed-loop depolymerization of PLA in water. The hydrogel bandage exhibits strong adhesion to skin while it can be removed with no residue by water flushing, showing protection to neo-tissue during dressing replacement. The in vivo application of the hydrogel bandage significantly promoted wound healing by closing the wound, forming a physical barrier, and providing an anti-inflammatory effect, showing great potential in future clinical applications.

Tuning Electron Transfer Coupling and Exchange Interaction in Bis-triarylamine Radical Cations and Dications by Bridge Electron Density

The influence of the electron density of a bridge connecting two redox centers on both the intervalence hole transfer and the magnetic superexchange was investigated in a series of bridged bis-triarylamine mono- and dications. In this series, the bridge was 2,7-fluorenyl, where the bridge electron density was modified by substituents at the 9-position. For the mixed-valence monocations, the observation of both an intervalence charge transfer (IVCT) band and an absorption band associated with an electron transfer from the bridging fluorene to the triarylamine radical cation centers allowed determination of the electron transfer couplings in the framework of the three-state generalized Mulliken-Hush theory. Comparison of the derived couplings with those obtained from a classical two-state approach demonstrates an enhancement of the electronic coupling which increases with decreasing bridge state energy. For the dicationic diradical counterparts, the singlet-triplet gap (exchange interaction) was determined both experimentally and by quantum chemical methods. Hereby, an increase of antiferromagnetic coupling with a lowering of the bridge state energy by electron donating substituents was observed. Analysis of the involved molecular orbitals suggests that the ferromagnetic coupling is inversely proportional to the square of the bridge energy, which is also supported by the experimental findings. This influence of the bridge state energy on both types of interactions, electron transfer and magnetic exchange, provides a design guideline for fine-tuning the properties of electronically coupled organic redox dyads by variation of the bridge electron density.

Automatic Differentiation for Explicitly Correlated MP2

Automatic differentiation (AD) offers a route to achieve arbitrary-order derivatives of challenging wave function methods without the use of analytic gradients or response theory. Currently, AD has been predominantly used in methods where first- and/or second-order derivatives are available, but it has not been applied to methods lacking available derivatives. The most robust approximation of explicitly correlated MP2, MP2-F12/3C(FIX)+CABS, is one such method. By comparing the results of MP2-F12 computed with AD versus finite-differences, it is shown that (a) optimized geometries match to about 10-3 Å for bond lengths and a 10-6 degree for angles, and (b) dipole moments match to about 10-6 D. Hessians were observed to have poorer agreement with numerical results (10-5), which is attributed to deficiencies in AD implementations currently. However, it is notable that vibrational frequencies match within 10-2 cm-1. The use of AD also allowed the prediction of MP2-F12/3C(FIX)+CABS IR intensities for the first time.

Elucidation of the exchange interaction in photoexcited three-spin systems - a second-order perturbational approach

Photogenerated three-spin systems show great potential for applications in the field of molecular spintronics. In these systems, the exchange interaction in the electronically excited state dictates their magnetic properties. To design such molecules for specific applications, it is thus important to understand how the sign and magnitude of the exchange interaction can be controlled. For this purpose, we developed a perturbational approach, based on previous work by the groups of de Loth and Malrieu, that allows for the direct calculation of the exchange interaction and its individual contributions up to the second order and implemented it within the ORCA program package. Within this manuscript, we present the derivation of the individual second-order contributions, provide an overview of the implementation of the code and illustrate its performance. We show that, using this perturbational approach in combination with state-averaged orbitals from minimal active space calculations, accurate values for the exchange interaction can be computed for organic nitroxides. Further, we demonstrate that the weight of the ionic determinants in the orbital optimisation of the CASSCF procedure is crucial for the computation of accurate exchange couplings. In the case of photoexcited chromophore-radical systems, we find that the dynamic spin polarisation effect constitutes the most important contribution to the exchange interaction, whereby the sign of this contribution determines the sign of the exchange interaction.

Solar Light CO2 Photoreduction Enhancement by Mononuclear Rhenium(I) Complexes: Characterization and Mechanistic Insights

The catalytic efficacy of a novel mononuclear rhenium(I) complex in CO2 reduction is remarkable, with a turnover number (TONCO) of 1517 in 3 h, significantly outperforming previous Re(I) catalysts. This complex, synthesized via a substitution reaction on an aromatic ring to form a bromo-bipyridine derivative, L1 = 2-bromo-6-(1H-pyrazol-1-yl)pyridine, and further reacting with [Re(CO)5Cl], results in the facial-tricarbonyl complex [ReL1(CO)3Cl] (1). The light green solid was obtained with an 80% yield and thoroughly characterized using cyclic voltammetry, nuclear magnetic resonance (NMR), Fourier transform infrared (FTIR) spectroscopy, and ultraviolet-visible (UV-vis) spectroscopy. Cyclic voltammetry under CO2 atmosphere revealed three distinct redox processes, suggesting the formation of new electroactive compounds. The studies on photoreduction highlighted the ability of the catalyst to reduce CO2, while NMR, FTIR, and electrospray ionization (ESI) mass spectrometry provided insights into the mechanism, revealing the formation of solvent-coordinated complexes and new species under varying conditions. Additionally, computational studies (DFT) were undertaken to better understand the electronic structure and reactivity patterns of 1, focusing on the role of the ligand, the spectroscopic features, and the redox behavior. This comprehensive approach provides insights into the intricate dynamics of CO2 photoreduction, showcasing the potential of Re(I) complexes in catalysis.

Nature of Charge Transfer Effects in Complexes of Dopamine Derivatives Adsorbed on Graphene-Type Nanostructures

Continuing the investigation started for dopamine (DA) and dopamine-o-quinone (DoQ) (see, the light absorption and charge transfer properties of the dopamine zwitterion (called dopamine-semiquinone or DsQ) adsorbed on the graphene nanoparticle surface is investigated using the ground state and linear-response time-dependent density functional theories, considering the ωB97X-D3BJ/def2-TZVPP level of theory. In terms of the strength of molecular adsorption on the surface, the DsQ form has 50% higher binding energy than that found in our previous work for the DA or DoQ cases (-20.24 kcal/mol vs. -30.41 kcal/mol). The results obtained for electronically excited states and UV-Vis absorption spectra show that the photochemical behavior of DsQ is more similar to DA than that observed for DoQ. Of the three systems analyzed, the DsQ-based complex shows the most active charge transfer (CT) phenomenon, both in terms of the number of CT-like states and the amount of charge transferred. Of the first thirty electronically excited states computed for the DsQ case, eleven are purely of the CT type, and nine are mixed CT and localized (or Frenkel) excitations. By varying the adsorption distance between the molecule and the surface vertically, the amount of charge transfer obtained for DA decreases significantly as the distance increases: for DoQ it remains stable, for DsQ there are states for which little change is observed, and for others, there is a significant change. Furthermore, the mechanistic compilation of the electron orbital diagrams of the individual components cannot describe in detail the nature of the excitations inside the complex.

Synthesis, biological evaluation, theoretical calculations, QSAR and molecular docking studies of novel arylaminonaphthols as potent antioxidants and BChE inhibitors

A completely green protocol was developed for the synthesis of a series of arylaminonaphthol derivatives in the presence of N-ethylethanolamine (NEEA) as a catalyst under ultrasonic irradiation and solventless conditions. The major assets of this methodology were the use of non-toxic organic medium, available catalyst, mild reaction condition, and good to excellent yield of desired products. All of the synthesized products were screened for their in vitro antioxidant activity using DPPH, ABTS, and Ferric-phenanthroline assays and it was found that most of them are potent antioxidant agents. Also, their butyrylcholinesterase inhibitory activity has been investigated in vitro. All tested compounds exhibited potential inhibitory activity toward BuChE when compared to standard reference drug galantamine, however, compounds 4r, 4u, 4 g and 4x gave higher butyrylcholinesterase inhibitory with IC50 values of 14.78 ± 0.65 µM, 16.18 ± 0.50 µM, 20.00 ± 0.50 µM, and 20.28 ± 0.08 µM respectively. On the other hand, we employed density functional theory (DFT), calculations to analyze molecular geometry and global reactivity descriptors, and MESP analysis to predict electrophilic and nucleophilic attacks. A quantitative structure-activity relationship (QSAR) investigation was conducted on the antioxidant and butyrylcholinesterase properties of 25 arylaminonaphthol derivatives, resulting in robust and satisfactory models. To evaluate their anti-Alzheimer's activity, compounds 4 g, 4q, 4r, 4u, and 4x underwent docking simulations at the active site of the acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), revealing why these compounds displayed superior activity, consistent with the biological findings.

dxtb-An efficient and fully differentiable framework for extended tight-binding

Automatic differentiation (AD) emerged as an integral part of machine learning, accelerating model development by enabling gradient-based optimization without explicit analytical derivatives. Recently, the benefits of AD and computing arbitrary-order derivatives with respect to any variable were also recognized in the field of quantum chemistry. In this work, we present dxtb-an open-source, fully differentiable framework for semiempirical extended tight-binding (xTB) methods. Developed entirely in Python and leveraging PyTorch for array operations, dxtb facilitates extensibility and rapid prototyping while maintaining computational efficiency. Through comprehensive code vectorization and optimization, we essentially reach the speed of compiled xTB programs for high-throughput calculations of small molecules. The excellent performance also scales to large systems, and batch operability yields additional benefits for execution on parallel hardware. In particular, energy evaluations are on par with existing programs, whereas the speed of automatically differentiated nuclear derivatives is only 2 to 5 times slower compared to their analytical counterparts. We showcase the utility of AD in dxtb by calculating various molecular and spectroscopic properties, highlighting its capacity to enhance and simplify such evaluations. Furthermore, the framework streamlines optimization tasks and offers seamless integration of semiempirical quantum chemistry in machine learning, paving the way for physics-inspired end-to-end differentiable models. Ultimately, dxtb aims to further advance the capabilities of semiempirical methods, providing an extensible foundation for future developments and hybrid machine learning applications. The framework is accessible at https://github.com/grimme-lab/dxtb.

Trapping an Elusive Fe(IV)-Superoxo Intermediate Inside a Self-Assembled Nanocage in Water at Room Temperature

Molecular cavities that mimic natural metalloenzymes have shown the potential to trap elusive reaction intermediates. Here, we demonstrate the formation of a rare yet stable Fe(IV)-superoxo intermediate at room temperature subsequent to dioxygen binding at the Fe(III) site of a (Et4N)2[FeIII(Cl)(bTAML)] complex confined inside the hydrophobic interior of a water-soluble Pd6L412+ nanocage. Using a combination of electron paramagnetic resonance, Mössbauer, Raman/IR vibrational, X-ray absorption, and emission spectroscopies, we demonstrate that the cage-encapsulated complex has a Fe(IV) oxidation state characterized by a stable S = 1/2 spin state and a short Fe-O bond distance of ∼1.70 Å. We find that the O2 reaction in confinement is reversible, while the formed Fe(IV)-superoxo complex readily reacts when presented with substrates having weak C-H bonds, highlighting the lability of the O-O bond. We envision that such optimally trapped high-valent superoxos can show new classes of reactivities catalyzing both oxygen atom transfer and C-H bond activation reactions.

The sugar cube: Network control and emergence in stereoediting reactions

Stereochemical editing strategies have recently enabled the transformation of readily accessible substrates into rare and valuable products. Typically, site selectivity is achieved by minimizing kinetic complexity by using protecting groups to suppress reactivity at undesired sites (substrate control) or by using catalysts with tailored shapes to drive reactivity at the desired site (catalyst control). We propose "network control," a contrasting paradigm that exploits hidden interactions between rate constants to greatly amplify modest intrinsic biases and enable precise multisite editing. When network control is applied to the photochemical isomerization of hexoses, six of the eight possible diastereomers can be selectively obtained. The amplification effect can be viewed as a mesoscale phenomenon between the limiting regimes of kinetic control in simple chemical systems and metabolic regulation in complex biological systems.

3-center and 4-center 2-particle Gaussian AO integrals on modern accelerated processors

We report an implementation of the McMurchie-Davidson (MD) algorithm for 3-center and 4-center 2-particle integrals over Gaussian atomic orbitals (AOs) with low and high angular momenta l and varying degrees of contraction for graphical processing units (GPUs). This work builds upon our recent implementation of a matrix form of the MD algorithm that is efficient for GPU evaluation of 4-center 2-particle integrals over Gaussian AOs of high angular momenta (l ≥ 4) [A. Asadchev and E. F. Valeev, J. Phys. Chem. A 127, 10889-10895 (2023)]. The use of unconventional data layouts and three variants of the MD algorithm allow for the evaluation of integrals with double precision and sustained performance between 25% and 70% of the theoretical hardware peak. Performance assessment includes integrals over AOs with l ≤ 6 (a higher l is supported). Preliminary implementation of the Hartree-Fock exchange operator is presented and assessed for computations with up to a quadruple-zeta basis and more than 20 000 AOs. The corresponding C++ code is part of the experimental open-source LibintX library available at https://github.com/ValeevGroup/libintx.

Restricted Open-Shell Hartree-Fock Method for a General Configuration State Function Featuring Arbitrarily Complex Spin-Couplings

In this work, we present a general spin restricted open-shell Hartree-Fock (ROHF) implementation that is able to generate self-consistent field (SCF) wave functions for an arbitrary configuration state function (CSF). These CSFs can contain an arbitrary number of unpaired electrons in arbitrary spin-couplings. The resulting method is named CSF-ROHF. We demonstrate that starting from the ROHF energy expression, for example, the one given by Edwards and Zerner, it is possible to obtain the values of the ROHF vector-coupling coefficients by setting up an open-shell for each group of consecutive parallel-coupled spins dictated by the unique spin-coupling pattern of any given CSF. To achieve this important and nontrivial goal, we employ the machinery of the iterative configuration expansion configuration interaction (ICE-CI) method, which is able to tackle general CI problems on the basis of spin-adapted CSFs. This development allows for the efficient generation of SCF spin-eigenfunctions for systems with complex spin-coupling patterns, such as polymetallic chains and metal clusters, while maintaining SCF scaling with system size (quadratic or less, depending on the specific algorithm and approximations chosen).

Code generation in ORCA: progress, efficiency and tight integration

An improved version of ORCA's automated generator environment (ORCA-AGE II) is presented. The algorithmic improvements and the move to C++ as the programming language lead to a performance gain of up to two orders of magnitude compared to the previously developed PYTHON toolchain. Additionally, the restructured modular design allows for far more complex code engines to be implemented readily. Importantly, we have realised an extremely tight integration with the ORCA host program. This allows for a workflow in which only the wavefunction Ansatz is part of the source code repository while all actual high-level code is generated automatically, inserted at the appropriate place in the host program before it is compiled and linked together with the hand written code parts. This construction ensures longevity and uniform code quality. Furthermore the new developments allow ORCA-AGE II to generate parallelised production-level code for highly complex theories, such as fully internally contracted multireference coupled-cluster theory (fic-MRCC) with an enormous number of contributing tensor contractions. We also discuss the automated implementation of nuclear gradients for arbitrary theories. All these improvements enable the implementation of theories that are too complex for the human mind and also reduce development times by orders of magnitude. We hope that this work enables researchers to concentrate on the intellectual content of the theories they develop rather than be concerned with technical details of the implementation.

Intramolecular Charge Transfer and Stimuli-Responsive Emission in Cholesterol-Appended Phenothiazine-Cyanostyryl-Based Donor-Acceptor Systems

Organic fluorescent molecules have received considerable attention owing to their various optoelectronic applications. Herein, we report the design and synthesis of two cholesterol-functionalized cyanostyrene-phenothiazine-based D-π-A systems that are emissive in both the solution and solid states. The newly synthesized cholesterol-appended phenothiazine-cyanostyrene diads PTCS-1 and PTCS-2 vary in the N-alkylation of phenothiazine, respectively, with─octyl and─hexyl chains. Both molecules are highly fluorescent and show reasonably good quantum yields in nonpolar solvents because of twisted intramolecular charge transfer (TICT). The molecules exhibit aggregation-induced emission in the solid state. Due to the presence of flexible alkyl chains in the phenothiazine and cholesterol moieties, PTCS-1 and PTCS-2 show mechanochromic luminescence switching in response to external shear stress and emission recovery under methanol vapor. Powder X-ray diffraction studies prove that the emission switching on the applied stimuli in both PTCS-1 and PTCS-2 is attributed to the reversible transformation between the crystalline and amorphous states. Time-dependent density functional theory (TD-DFT) studies are carried out to gain insight into the ICT interactions. TD-DFT analysis at the TD-M06-2X/def2-TZVP level further revealed that in both molecules, the lowest unoccupied molecular orbital (LUMO) + 2, LUMO, highest occupied molecular orbital (HOMO), and HOMO - 1 orbitals are responsible for the charge transfer interactions. These ICT interactions are identified as π-π* type interactions.

Metal free synthesis of thermally stable blue fluorescent m-terphenyls by ring transformation of 2H-pyran-2-ones: chemical synthesis, spectroscopy and computational studies

A simple yet convenient nucleophile-induced synthetic route for the construction of thermally stable fluorescent active functionalized stilbenes has been delineated. The nucleophile-induced base encouraged synthetic protocol was performed under mild conditions without harming the environment and products were achieved in good yields. The synthesized stilbenes showed amazing emission properties and good thermal stability. Synthesized products displayed interesting positive solvatochromism in different solvents based on variation in polarity. Further, we present a comprehensive analysis of the eight molecules, leveraging a combination of Density Functional Tight Binding (DFTB), Density Functional Theory (DFT) calculations, and Molecular Dynamics (MD) simulations. This integrated approach allowed for an in-depth exploration of the electronic structures, reactivity profiles, and dynamic behaviors of these complex molecular systems. Our findings reveal significant insights into the physicochemical properties of the synthesized molecules, contributing to a deeper understanding of their potential applications in various fields.

Optimal Oriented External Electric Fields to Trigger a Barrierless Oxaphosphetane Ring Opening Step of the Wittig Reaction

The Wittig reaction is one of the most important processes in organic chemistry for the asymmetric synthesis of olefinic compounds. In view of the increasingly acknowledged potentiality of the electric fields in promoting reactions, here we will consider the effect of the oriented external electric field (OEEF) on the second step of Wittig reaction (i. e. the ring opening oxaphosphetane) in a model system for non-stabilized ylides. In particular, we have determined the optimal direction and strength of the electric field that should be applied to annihilate the reaction barrier of the ring opening through the polarizable molecular electric dipole (PMED) model that we have recently developed. We conclude that the application of the optimal external electric field for the oxaphosphetane ring opening favours a Bestmann-like mechanism.

The synthesis and characterization of an iron(VII) nitrido complex

Complexes of iron in high oxidation states are captivating research subjects due to their pivotal role as active intermediates in numerous catalytic processes. Structural and spectroscopic studies of well-defined model complexes often provide evidence of these intermediates. In addition to the fundamental molecular and electronic structure insights gained by these complexes, their reactivity also affects our understanding of catalytic reaction mechanisms for small molecule and bond-activation chemistry. Here, we report the synthesis, structural and spectroscopic characterization of a stable, octahedral Fe(VI) nitrido complex and an authenticated, unique Fe(VII) species, prepared by one-electron oxidation. The super-oxidized Fe(VII) nitride rearranges to an Fe(V) imide through an intramolecular amination mechanism and ligand exchange, which is characterized spectroscopically and computationally. This enables combined reactivity and stability studies on a single molecular system of a rare high-valent complex redox pair. Quantum chemical calculations complement the spectroscopic parameters and provide evidence for a diamagnetic (S = 0) d 2 Fe(VI) and a genuine S = 1/2, d 1 Fe(VII) configuration of these super-oxidized nitrido complexes.

Extraction-less green spectrofluorimetric method for determination of mercury using an Astra Phloxine fluorophore: Comprehensive experimental and theoretical studies

A new sensitive and selective extraction-less method for the spectrofluorimetric determination of Hg2+ has been developed using the polymethine dye Astra Phloxine (AP) with requirements for green analytical chemistry. The principle of this method for analytical use is the formation of an ion associate (IA) in the presence of mercury ions and AP, which causes a change in the electronic structure of the dye cation and its subsequent fluorescence quenching. The IA created during the process was sufficiently stable in aqueous solutions and did not require the use of surfactants or organic solvents, which are typically used for similar analytical systems. The system indicated high selectivity towards Hg2+ in the presence of at least 100-fold higher concentrations of Co2+, Fe2+, Mg2+, Ni2+, Zn2+, Ca2+, Ba2+, K+, Na+. The calibration curve was linear in the range from 0.003 to 0.1 mg L-1. Precision and accuracy were expressed as intra-day and inter-day analysis with RSD values of 2.9 - 4.8 % and recovery of 97.1 - 110.0 %. The proposed method was used for the determination of Hg2+ on the model and real water samples. Theoretical calculations were used to predict the IA structure and to explain the experimentally observed fluorescence spectra. A few theoretical methods were tested for predicting IA emission energies, and the CIS(D) method was found to be the most accurate.

Towards routine organic structure determination using Raman microscopy

Raman microscopy can reveal a compound-specific vibrational "fingerprint" from micrograms of material with no sample preparation. We expect this increasingly available instrumentation to routinely assist synthetic chemists in structure determination; however, interpreting the information-dense spectra can be challenging for unreported compounds. Appropriate theoretical calculations using the highly efficient r2SCAN-3c method can accurately predict peak positions but are less precise in matching peak heights. To limit incorrect biases while comparing experimental and theoretical spectra, we introduce a user-friendly software that gives a match score to assist with structure determination. The capabilities and limitations of this approach are demonstrated for several proof-of-concept examples including the characterization of intermediates in the total synthesis of deoxyaspidodispermine.

High-Performance Evaluation of High Angular Momentum 4-Center Gaussian Integrals on Modern Accelerated Processors

We present a high-performance evaluation method for 4-center 2-particle integrals over Gaussian atomic orbitals with high angular momenta (l ≥ 4) and arbitrary contraction degrees on graphical processing units (GPUs) and other accelerators. The implementation uses the matrix form of McMurchie-Davidson recurrences. Evaluation of the four-center integrals over four l = 6 (i) Gaussian AOs in double precision (FP64) on an NVIDIA V100 GPU outperforms the reference implementation of the Obara-Saika recurrences (Libint) running on a single Intel Xeon core by more than a factor of 1000, easily exceeding the 73:1 ratio of the respective hardware peak FLOP rates while reaching almost 50% of the V100 peak. The approach can be extended to support AOs with even higher angular momenta; for lower angular momenta (l ≤ 3), additional improvements will be reported elsewhere. The implementation is part of an open-source LibintX library freely available at github.com:ValeevGroup/LibintX.

A Cooperative Cobalt-Driven System for One-Carbon Extension in the Synthesis of (Z)-Silyl Enol Ethers from Aldehydes: Unlocking Regio- and Stereoselectivity

The research presented herein explores a cobalt-based catalytic system, distinctively featuring a cooperative boron-centric element within its intricate ligand architecture. This system is strategically engineered to enable the integration of a singular carbon atom into aldehydes, a process culminating in the production of (Z)-silyl enol ethers. Beyond offering an efficient one-pot synthesis route, this method adeptly overcomes challenges inherent to conventional techniques, such as the need for large amounts of additives, restrictive functional group tolerance, and extreme reaction temperatures. Initial mechanistic studies suggest the potential role of a cobalt-carbene complex as a catalytically significant species and underscore the importance of the borane segment. Collectively, these observations highlight the potential of this system in advancing complex bond activation pursuits.

Modulating Aluminum Solvation with Ionic Liquids for Improved Aqueous-Based Aluminum-Ion Batteries

Aqueous-based Al-ion batteries are attractive alternatives to Li-ion batteries due to their safety, high volumetric energy density, abundance, and recyclability. Although aluminum-ion batteries are attractive, there are major challenges to overcome, which include understanding the nature of the passive layer of aluminum oxide on the aluminum anode, the narrow electrochemical window of aqueous electrolytes, and lack of suitable cathodes. Here, we report using experiments in conjunction with DFT simulations to clarify the role of ionic liquids (ILs) in altering the Al solvation dynamics, which in turn affects the aluminum electrochemistry and aqueous-based battery performance significantly. DFT calculations showed that the addition of 1-ethyl-3-methylimidazolium trifluoromethylsulfonate (EMIMTfO) changes the aluminum solvation structure in the aqueous (Al(TfO)3) electrolyte to lower coordinated solvation shells, thereby influencing and improving Al deposition/stripping on the Zn/Al alloy anode. Furthermore, the addition of an IL reduces the strain in manganese oxide during intercalation/deintercalation, thereby improving the Zn/Al-MnOx battery performance. By optimizing the electrolyte composition, a battery potential of >1.7 V was achieved for the Zn/Al-MnOx system.

Probing the magnetic and magneto-optical properties of a radical-bridged Tb4 single-molecule magnet

Reaction of the 1,2,4,5-tetrazine (tz˙-) radical and {Cp*2Tb}+ has yielded a tetranuclear radical-bridged TbIII single-molecule magnet. The strong Ln-radical coupling and the electronic differences of the TbIII ions in [(Cp*2Tb)4(tz˙-)4]·3C6H6 (1) are probed via magnetic, magneto-optical and computational studies.

An algorithm to find the optimal oriented external electrostatic field for annihilating a reaction barrier in a polarizable molecular system

The use of oriented external electric fields (OEEFs) to promote and control chemical reactivity has motivated many theoretical and computational studies in the last decade to model the action of OEEFs on a molecular system and its effects on chemical processes. Given a reaction, a central goal in this research area is to predict the optimal OEEF (oOEEF) required to annihilate the reaction energy barrier with the smallest possible field strength. Here, we present a model rooted in catastrophe and optimum control theories that allows us to find the oOEEF for a given reaction valley in the potential energy surface (PES). In this model, the effective (or perturbed) PES of a polarizable molecular system is constructed by adding to the original, non-perturbed, PES a term accounting for the interaction of the OEEF with the intrinsic electric dipole and polarizability of the molecular system, so called the polarizable molecular electric dipole (PMED) model. We demonstrate that the oOEEF can be established by locating a point in the original PES with unique topological properties: the optimal barrier breakdown or bond-breaking point (oBBP). The essential feature of the oBBP structure is the fact that this point maintains its topological properties for all the applied OEEFs, also for the unperturbed PES, thus becoming much more relevant than the commonly used minima and transition state structures. The PMED model proposed here has been implemented in an open access package and is shown to successfully predict the oOEEF for two processes: an isomerization reaction of a cumulene derivative and the Huisgen cycloaddition reaction.

Accurate Evaluation of Dispersion Energies at Coupled Cluster Level to Understand the Substituent Effects in Am(III) and Eu(III) Complexes

The effect of cyclic and aromatic substituents on the complexation behavior of phosphine oxide ligands with Am(III) and Eu(III) was investigated at density functional theory (DFT) and domain-based local pair natural orbital coupled-cluster (DLPNO-CC) levels. Combining DFT with accurate coupled cluster methods, we have evaluated the dispersion energy contributions to the complexation energies for trivalent Am and Eu complexes for the first time. Irrespective of the nature of substituents on the P atom, the electronic structure of the P═O group remains identical in all of the ligands. The study reveals the importance of dispersion interactions during complexation and is estimated to be more significant for Am(III) than for Eu(III) complexes.

Memory-Efficient Recursive Evaluation of 3-Center Gaussian Integrals

To improve the efficiency of Gaussian integral evaluation on modern accelerated architectures, FLOP-efficient Obara-Saika-based recursive evaluation schemes are optimized for the memory footprint. For the 3-center 2-particle integrals that are key for the evaluation of Coulomb and other 2-particle interactions in the density-fitting approximation, the use of multiquantal recurrences (in which multiple quanta are created or transferred at once) is shown to produce significant memory savings. Other innovations include leveraging register memory for reduced memory footprint and direct compile-time generation of optimized kernels (instead of custom code generation) with compile-time features of modern C++/CUDA. Performance of conventional and CUDA-based implementations of the proposed schemes is illustrated for both the individual batches of integrals involving up to Gaussians with low and high angular momenta (up to L = 6) and contraction degrees, as well as for the density-fitting-based evaluation of the Coulomb potential. The computer implementation is available in the open-source LibintX library.

A catastrophe theory-based model for optimal control of chemical reactions by means of oriented electric fields

The effect of oriented external electric fields (OEEF) on chemical reactivity has been studied theoretically and computationally in the last decades. A central goal in this research area is to predict the orientation and the smallest amplitude electric field that renders a barrierless chemical process with the smallest possible strength. Recently, a model to find the optimal electric field has been proposed and described (Bofill JM et al., J. Chem. Theory Comput. 18:935, 2022). We here proof that this model is based on catastrophe and optimum control theories. Based on both theories a technical treatment of the model is given and applied to a two-dimensional generic example that provides insight into its nature and capability. Finally, the model is applied to determine the optimal OEEF for the trans-to-cis isomerization of a [3]cumulene derivative.

Photocatalytic Activity of the V2O5 Catalyst toward Selected Pharmaceuticals and Their Mixture: Influence of the Molecular Structure on the Efficiency of the Process

Due to the inability of conventional wastewater treatment procedures to remove organic pharmaceutical pollutants, active pharmaceutical components remain in wastewater and even reach tap water. In terms of pharmaceutical pollutants, the scientific community focuses on β-blockers due to their extensive (over)usage and moderately high solubility. In this study, the photocatalytic activity of V2O5 was investigated through the degradation of nadolol (NAD), pindolol (PIN), metoprolol (MET), and their mixture under ultraviolet (UV) irradiation in water. For the preparation of V2O5, facile hydrothermal synthesis was used. The structural, morphological, and surface properties and purity of synthesized V2O5 powder were investigated by scanning electron microscopy (SEM), X-ray, and Raman spectroscopy. SEM micrographs showed hexagonal-shaped platelets with well-defined morphology of materials with diameters in the range of 10−65 µm and thickness of around a few microns. X-ray diffraction identified only one crystalline phase in the sample. The Raman scattering measurements taken on the catalyst confirmed the result of XRPD. Degradation kinetics were monitored by ultra-fast liquid chromatography with diode array detection. The results showed that in individual solutions, photocatalytic degradation of MET and NAD was relatively insignificant (<10%). However, in the PIN case, the degradation was significant (64%). In the mixture, the photodegradation efficiency of MET and NAD slightly increased (15% and 13%). Conversely, it reduced the PIN to the still satisfactory value of 40%. Computational analysis based on molecular and periodic density functional theory calculations was used to complement our experimental findings. Calculations of the average local ionization energy indicate that the PIN is the most reactive of all three considered molecules in terms of removing an electron from it.