Accurate QM/MM Molecular Dynamics for Periodic Systems in GPU4PySCF with Applications to Enzyme Catalysis

We present an implementation of the quantum mechanics/molecular mechanics (QM/MM) method for periodic systems using GPU accelerated QM methods, a distributed multipole formulation of the electrostatics, and a pseudobond treatment of the QM/MM boundary. We demonstrate that our method has well-controlled errors, stable self-consistent QM convergence, and energy-conserving dynamics. We further describe an application to the catalytic kinetics of chorismate mutase. Using an accurate hybrid functional reparametrized to coupled cluster energetics, our QM/MM simulations highlight the sensitivity in the calculated rate to the choice of quantum method, quantum region selection, and local protein conformation. Our work is provided through the open-source PySCF package using acceleration from the GPU4PySCF module.

Analytical derivative approaches for vibro-polaritonic structures and properties. I. Formalism and implementation

Vibro-polaritons are hybrid light-matter states that arise from the strong coupling between the molecular vibrational transitions and the photons in an optical cavity. Developing theoretical and computational methods to describe and predict the unique properties of vibro-polaritons is of great significance for guiding the design of new materials and experiments. Here, we present the ab initio cavity Born-Oppenheimer density functional theory (CBO-DFT) and formulate the analytic energy gradient and Hessian as well as the nuclear and photonic derivatives of dipole and polarizability within the framework of CBO-DFT to efficiently calculate the harmonic vibrational frequencies, infrared absorption, and Raman scattering spectra of vibro-polaritons as well as to explore the critical points on the cavity potential energy surface. The implementation of analytic derivatives into the electronic structure package is validated by a comparison with the finite-difference method and with other reported computational results. By adopting appropriate exchange-correlation functionals, CBO-DFT can better describe the structure and properties of molecules in the cavity than CBO-Hartree-Fock method. It is expected that CBO-DFT is a useful tool for studying the polaritonic structures and properties.

Data-Driven Improvement of Local Hybrid Functionals: Neural-Network-Based Local Mixing Functions and Power-Series Correlation Functionals

Local hybrid functionals (LHs) use a real-space position-dependent admixture of exact exchange (EXX), governed by a local mixing function (LMF). The systematic construction of LMFs has been hampered over the years by a lack of exact physical constraints on their valence behavior. Here, we exploit a data-driven approach and train a new type of "n-LMF" as a relatively shallow neural network. The input features are of meta-GGA character, while the W4-17 atomization-energy and BH76 reaction-barrier test sets have been used for training. Simply replacing the widely used "t-LMF" of the LH20t functional by the n-LMF provides the LH24n-B95 functional. Augmented by DFT-D4 dispersion corrections, LH24n-B95-D4 remarkably improves the WTMAD-2 value for the large GMTKN55 test suite of general main-group thermochemistry, kinetics, and noncovalent interactions (NCIs) from 4.55 to 3.49 kcal/mol. As we found the limited flexibility of the B95c correlation functional to disfavor much further improvement on NCIs, we proceeded to replace it by an optimized B97c-type power-series expansion. This gives the LH24n functional. LH24n-D4 gives a WTMAD-2 value of 3.10 kcal/mol, the so far lowest value of a rung 4 functional in self-consistent calculations. The new functionals perform moderately well for organometallic transition-metal energetics while leaving room for further data-driven improvements in that area. Compared to complete neural-network functionals like DM21, the present more tailored approach to train just the LMF in a flexible but well-defined human-designed LH functional retains the possibility of graphical LMF analyses to gain deeper understanding. We find that both the present n-LMF and the recent x-LMF suppress the so-called gauge problem of local hybrids without adding a calibration function as required for other LMFs. LMF plots show that this can be traced back to large LMF values in the small-density region between the interacting atoms in NCIs for n- and x-LMFs and low values for the t-LMF. We also find that the trained n-LMF has relatively large values in covalent bonds without deteriorating binding energies. The current approach enables fast and efficient routine self-consistent calculations using n-LMFs in Turbomole.

Diels-Alder Reactivity of Triisopropylsilyl Ethynyl Substituted Acenes

We investigated the Diels-Alder reaction of 6,13-bis(triisopropylsilylethynyl)pentacene (1) with small dienophiles such as (bridged) dihydronaphthalenes/cyclohexenes that yielded adducts at the central ring, the other dienophiles predominantly or exclusively attacked the unsubstituted off-center ring. The difference in regioselectivity was investigated by DFT calculations. Apart from dispersion interactions, it is due to the steric demand of the dienophiles, which need to fit in between the silylethynyl substituents to react at the central ring. Epoxynaphthalene adducts of 1 as well as its anthracene and tetracene congeners were deoxygenated, easily furnishing triarenobarrelenes with TIPS-ethynyl substituents at the bridgeheads, attractive building blocks for porous solids and higher acene-based trimers.

Theoretical Study of Extensive Hydrogen Abstraction Reactions for 2-Hydroxyethylhydrazine (HEH)

Energetic ionic liquids have a high potential to replace the traditional monopropellant hydrazine as a high-energy green propellant and can be widely used in aerospace technology. A high-energy ionic liquid─HEHN has also gained extensive attention from researchers. To explore the reaction mechanism of HEHN and establish a chemical kinetic model for high-energy ionic liquid propellants, 28 hydrogen abstraction reactions of HEH, which is the main decomposition product of HEHN, were investigated in this study. Seven abstractors were involved, including •H, •OH, NO2, HO2•, •CH3, CH3O•, and CH3O2•. In the case of ab initio calculations, the M06-2X/6311++G(d,p) approach was utilized for geometry optimization, determination of vibrational frequencies, and dihedral scans. The CCSD/cc-pVXZ (X = T, Q) level of theory was used to calculate the single-point energies (SPEs). The rate coefficients of all 28 reactions and the thermochemical parameters of all involved species were determined. The results indicate that the rate of hydrogen abstraction at the -NH site is faster than that at other sites at relatively low temperatures. For all four abstraction sites, HEH + •H, •OH, and CH3O• have higher reaction rates than HEH + CH3O2• and HO2•. In particular, NO2 systems at the -NH and -NH2 sites even begin to show higher reactivity than the •H, •OH, and CH3O• systems when the temperature is above ∼1100 K.

Photochemical pathways in astronomical ices: A computational study of singlet oxygen reactions with hydrocarbons

Complex organic molecules are widespread in different areas of the interstellar medium, including cold areas, such as molecular clouds, where chemical reactions occur in ice. Among the observed molecules are oxygen-bearing organic molecules, which are of high interest given their significant role in astrobiology. Despite the observed rich chemistry, the underlying molecular mechanisms responsible for molecular formation in such cold dilute areas are still not fully understood. In this paper, we study the unique chemistry taking place in astronomically relevant ices, where UV radiation is a central driving force for chemical reactions. Photofragmentation of ice components gives rise to highly reactive species, such as the O(1D) atom. These species provide a pathway for chemical complexity even in cold areas. Using quantum chemistry calculations, we demonstrate that O(1D) reacts barrierlessly with hydrocarbons. Moreover, photoprocessing of the reaction products (and other components of the ice), followed by radical recombination, is found to be an essential part of the overall mechanism. In ice containing O(1D) and hydrocarbons, the formation of formaldehyde in methane ice, acetaldehyde in ethane ice, and carbon monoxide in acetylene ice, and the consumption of alcohol in all systems, was predicted in agreement with experimental results.

Quick-and-Easy Validation of Protein-Ligand Binding Models Using Fragment-Based Semiempirical Quantum Chemistry

Electronic structure calculations in enzymes converge very slowly with respect to the size of the model region that is described using quantum mechanics (QM), requiring hundreds of atoms to obtain converged results and exhibiting substantial sensitivity (at least in smaller models) to which amino acids are included in the QM region. As such, there is considerable interest in developing automated procedures to construct a QM model region based on well-defined criteria. However, testing such procedures is burdensome due to the cost of large-scale electronic structure calculations. Here, we show that semiempirical methods can be used as alternatives to density functional theory (DFT) to assess convergence in sequences of models generated by various automated protocols. The cost of these convergence tests is reduced even further by means of a many-body expansion. We use this approach to examine convergence (with respect to model size) of protein-ligand binding energies. Fragment-based semiempirical calculations afford well-converged interaction energies in a tiny fraction of the cost required for DFT calculations. Two-body interactions between the ligand and single-residue amino acid fragments afford a low-cost way to construct a "QM-informed" enzyme model of reduced size, furnishing an automatable active-site model-building procedure. This provides a streamlined, user-friendly approach for constructing ligand binding-site models that requires neither a priori information nor manual adjustments. Extension to model-building for thermochemical calculations should be straightforward.

Convergent Protocols for Computing Protein-Ligand Interaction Energies Using Fragment-Based Quantum Chemistry

Fragment-based quantum chemistry methods offer a means to sidestep the steep nonlinear scaling of electronic structure calculations so that large molecular systems can be investigated using high-level methods. Here, we use fragmentation to compute protein-ligand interaction energies in systems with several thousand atoms, using a new software platform for managing fragment-based calculations that implements a screened many-body expansion. Convergence tests using a minimal-basis semiempirical method (HF-3c) indicate that two-body calculations, with single-residue fragments and simple hydrogen caps, are sufficient to reproduce interaction energies obtained using conventional supramolecular electronic structure calculations, to within 1 kcal/mol at about 1% of the computational cost. We also demonstrate that the HF-3c results are illustrative of trends obtained with density functional theory in basis sets up to augmented quadruple-ζ quality. Strategic deployment of fragmentation facilitates the use of converged biomolecular model systems alongside high-quality electronic structure methods and basis sets, bringing ab initio quantum chemistry to systems of hitherto unimaginable size. This will be useful for generation of high-quality training data for machine learning applications.

3,3-Difluorooxetane-A Versatile Functional Group for Bioisosteric Replacements in Drug Discovery

Functional group (FG) is one of the cornerstone concepts in organic chemistry and related areas. The wide spread of bioisosterism ideas in medicinal chemistry and beyond caused a striking rise in demand for novel FGs with a defined impact on the developed compound properties. In this work, the evaluation of the 3,3-difluorooxetane unit (3,3-diFox) as a functional group for bioisosteric replacements is disclosed. A comprehensive experimental study (including multigram building block synthesis, quantification of steric and electronic properties, measurements of pKa, LogP, chemical stability, and biological evaluation of the 3,3-diFox-derived bioisostere of a drug candidate) revealed a prominent behavior of the 3,3-diFox fragment as a versatile substituent for early drug discovery programs.

Canonical coupled cluster binding benchmark for nanoscale noncovalent complexes at the hundred-atom scale

In this study, we introduce two datasets for nanoscale noncovalent binding, featuring complexes at the hundred-atom scale, benchmarked using coupled cluster with single, double, and perturbative triple [CCSD(T)] excitations extrapolated to the complete basis set (CBS) limit. The first dataset, L14, comprises 14 complexes with canonical CCSD(T)/CBS benchmarks, extending the applicability of CCSD(T)/CBS binding benchmarks to systems as large as 113 atoms. The second dataset, vL11, consists of 11 even larger complexes, evaluated using the local CCSD(T)/CBS method with stringent thresholds, covering systems up to 174 atoms. We compare binding energies obtained from local CCSD(T) and fixed-node diffusion Monte Carlo (FN-DMC), which have previously shown discrepancies exceeding the chemical accuracy threshold of 1 kcal/mol in large complexes, with the new canonical CCSD(T)/CBS results. While local CCSD(T)/CBS agrees with canonical CCSD(T)/CBS within binding uncertainties, FN-DMC consistently underestimates binding energies in π-π complexes by over 1 kcal/mol. Potential sources of error in canonical CCSD(T)/CBS are discussed, and we argue that the observed discrepancies are unlikely to originate from CCSD(T) itself. Instead, the fixed-node approximation in FN-DMC warrants further investigation to elucidate these binding discrepancies. Using these datasets as reference, we evaluate the performance of various electronic structure methods, semi-empirical approaches, and machine learning potentials for nanoscale complexes. Based on computational accuracy and stability across system sizes, we recommend MP2+aiD(CCD), PBE0+D4, and ωB97X-3c as reliable methods for investigating noncovalent interactions in nanoscale complexes, maintaining their promising performance observed in smaller systems.

Analysis of Interaction Features of Cyclo[13]carbon with Small Molecules and Formation Mechanism of Its Dimer

The newly discovered cyclo[13]carbon, the first artificially synthesized odd-numbered carbon ring, is an intriguing carbon isomer that provides a valuable subject for studying low-symmetry carbon materials. In this work, we employed first-principles calculations to explore the geometric structure and electronic properties of cyclo[13]carbon through various techniques such as vibrational mode analysis, bond order analysis, spin density analysis, electron localization analysis, electrostatic potential and van der Waals potential analysis, visualization of weak interactions, and energy decomposition analysis. We investigated the interaction characteristics of cyclo[13]carbon with small molecules and examined its dimer formation mechanism and dynamics features using ab initio molecular dynamics. Our study reveals the unique physicochemical properties of this novel carbon ring system. The antiaromaticity of the low-symmetry cyclo[13]carbon sets it apart from previously synthesized even-numbered carbon rings, with van der Waals interactions playing a crucial role in its binding with small molecules and in the formation of C13 dimers. This research provides theoretical insights that complement experimental observations and theoretical studies, aiding further investigation into the diverse properties of fresh carbon material isomers and promoting the synthesis and application of novel molecular materials in molecular electronics and nanotechnology.

How accurate can Kohn-Sham density functional be for both main-group and transition metal reactions

Achieving chemical accuracy in describing reactions involving both main-group elements and transition metals poses a substantial challenge for density functional approximations (DFAs), primarily due to the significantly different behaviors for electrons moving in the s,p-orbitals or in the d,f-orbitals. MOR41, a representative dataset of transition metal chemistry, has highlighted the PWPB95-D3(BJ) functional, a B2PLYP-type doubly hybrid (bDH) approximation equipped with an empirical dispersion correction, as the leading functional thus far (Dohm et al., J Chem Theory Comput 2018;14: 2596-2608). However, this functional is not among the top bDH methods for main-group chemistry (Goerigk et al., Phys Chem Chem Phys. 2017;19: 32184). Conversely, bDH methods such as DSD-BLYP-D3, proficient in main-group chemistry, often falter for transition metal chemistry. Herein, taking advantage of the home-made Rust-based Electronic-Structure Toolkits, we examine a suite of XYG3-type doubly hybrid (xDH) methods. We confirm that the trade-off in descriptive accuracy between main-group and transition metal systems persists within the realm of perturbation theory (PT2)-based xDH methods. Notably, however, our study ushers in a pivotal advance with the recently proposed renormalized xDH method, R-xDH7-SCC15. This method not only distinguishes itself among the elite methods for main-group chemistry, but also achieves an unprecedented accuracy for the MOR41 dataset, outperforming all other reported DFAs. The efficacy of R-xDH7-SCC15 stems from the successful integration of a renormalized PT2 correlation model (rPT2) and a machine-learning strong-correlation correction (SCC15), marking a significant step forward in the realm of computational chemistry.

Spin parameter optimization for spin-polarized extended tight-binding methods

We present an optimization strategy for atom-specific spin-polarization constants within the spin-polarized GFN2-xTB framework, aiming to enhance the accuracy of molecular simulations. We compare a sequential and global optimization of spin parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. Sensitivity analysis using Sobol indices guides the identification of the most influential parameters for a given reference dataset, allowing for a nuanced understanding of their impact on diverse molecular properties. In the case of the W4-11 dataset, substantial error reduction was achieved, demonstrating the potential of the optimization. Transferability of the optimized spin-polarization constants over different properties, however, is limited, as we demonstrate by applying the optimized parameters on a set of singlet-triplet gaps in carbenes. Further studies on ionization potentials and electron affinities highlight some inherent limitations of current extended tight-binding methods that can not be resolved by simple parameter optimization. We conclude that the significantly improved accuracy strongly encourages the present re-optimization of the spin-polarization constants, whereas the limited transferability motivates a property-specific optimization strategy.

Advanced Charge Extended Hückel (CEH) Model and a Consistent Adaptive Minimal Basis Set for the Elements Z = 1-103

The Charge Extended Hückel (CEH) model, initially introduced for adaptive atomic orbital (AO) basis set construction (J. Chem. Phys. 2023, 159, 164108), has been significantly revised to enhance accuracy and robustness, particularly in challenging electronic situations. This revision includes an extension toward f-elements, covering actinoids with their f-electrons in the valence space. We present a novel noniterative approximation for the electrostatic contribution to the effective Fock matrix, which substantially improves performance in polar or charged systems. Additionally, the training data set for elements Z = 1-103 has been expanded to encompass even more chemically diverse reference molecules as well as dipole moments and shell populations in addition to atomic charges. It includes a greater variety of "mindless" molecules (MLMs) as well as more complex electronic structures through open-shell and highly charged species. The revised method achieves mean absolute errors for atomic charges q of approximately 0.02 e- for randomly selected (mostly organic) molecules and 0.09 e- for MLMs, outperforming both classical charge models and established tight-binding methods. Furthermore, the revised CEH model has been validated through density functional theory calculations with the updated adaptive q-vSZP AO basis set on common thermochemical databases. Consistent with the extension of the CEH model, q-vSZP has also been variationally optimized and tested for elements Z = 58-71 and 87-103. The original versions of both CEH and q-vSZP are now considered deprecated.

The diradicaloid electronic structure of dialumenes: a benchmark study at the Full-CI limit

Multiply-bonded main group compounds of groups 13-15 are attracting significant interest not only because they provide fundamental insight into the nature of metal-metal bonding, but also for their potential in small molecule bond activation and catalysis. This includes dialumenes, neutral Al(I) compounds that contain AlAl double bonds, which display high reactivity owing to their intrinsic diradicaloid character. The electronic structure of the simplest dialumene, Al2H2, is here analyzed up to a practical Full-CI limit using DMRG and selected CI methods for the bond dissociation energy (BDE), geometry and properties of the electron density (difference density, ELF). Acquiring Full-CI reference values for the simplest dialumene (but possessing the highest diradical character) allows for a rigorous benchmarking of simpler correlated wavefunction theory (WFT) methods and density functional methods in treating the electronic structure of such systems. Single-reference coupled cluster theory using a RHF reference is found to reliably converge to the Full-CI limit and CCSD(T) is fully capable of capturing the diradical character, while multi-reference methods offer no clear advantages. Density functional methods struggle to fully describe the electronic structure complexity although non-hybrid functionals such as TPSS come close. Solving the inverse Kohn-Sham problem for a Full-CI-quality density revealed minimal density-driven errors in the TPSS-calculated BDE unlike high-percentage hybrids such as M06-2X. No density functional, however, predicts accurate relative energies. Fractional-occupation density plots at the TPSS level correlate well with WFT-based diradical character metrics, a useful result for determining diradical character in larger systems.

Performance of quantum chemistry methods for a benchmark set of spin-state energetics derived from experimental data of 17 transition metal complexes (SSE17)

Accurate prediction of spin-state energetics for transition metal (TM) complexes is a compelling problem in applied quantum chemistry, with enormous implications for modeling catalytic reaction mechanisms and computational discovery of materials. Computed spin-state energetics are strongly method-dependent and credible reference data are scarce, making it difficult to conduct conclusive computational studies of open-shell TM systems. Here, we present a novel benchmark set of first-row TM spin-state energetics, which is derived from experimental data of 17 complexes containing FeII, FeIII, CoII, CoIII, MnII, and NiII with chemically diverse ligands. The estimates of adiabatic or vertical spin-state splittings, which are obtained from spin crossover enthalpies or energies of spin-forbidden absorption bands, suitably back-corrected for the vibrational and environmental effects, are employed as reference values for benchmarking density functional theory (DFT) and wave function methods. The results demonstrate a high accuracy of the coupled-cluster CCSD(T) method, which features the mean absolute error (MAE) of 1.5 kcal mol-1 and maximum error of -3.5 kcal mol-1, and outperforms all the tested multireference methods: CASPT2, MRCI+Q, CASPT2/CC and CASPT2+δMRCI. Switching from Hartree-Fock to Kohn-Sham orbitals is not found to consistently improve the CCSD(T) accuracy. The best performing DFT methods are double-hybrids (PWPB95-D3(BJ), B2PLYP-D3(BJ)) with the MAEs below 3 kcal mol-1 and maximum errors within 6 kcal mol-1, whereas the DFT methods so far recommended for spin states (e.g., B3LYP*-D3(BJ) and TPSSh-D3(BJ)) are found to perform much worse with the MAEs of 5-7 kcal mol-1 and maximum errors beyond 10 kcal mol-1. This work is the first such extensive benchmark study of quantum chemistry methods for TM spin-state energetics making use of experimental reference data. The results are relevant for the proper choice of methods to characterize TM systems in computational catalysis and (bio)inorganic chemistry, and may also stimulate new developments in quantum-chemical or machine learning approaches.

Data Quality in the Fitting of Approximate Models: A Computational Chemistry Perspective

Empirical parametrization underpins many scientific methodologies including certain quantum-chemistry protocols [e.g., density functional theory (DFT), machine-learning (ML) models]. In some cases, the fitting requires a large amount of data, necessitating the use of data obtained using low-cost, and thus low-quality, means. Here we examine the effect of using low-quality data on the resulting method in the context of DFT methods. We use multiple G2/97 data sets of different qualities to fit the DFT-type methods. Encouragingly, this fitting can tolerate a relatively large proportion of low-quality fitting data, which may be attributed to the physical foundations of the DFT models and the use of a modest number of parameters. Further examination using "ML-quality" data shows that adding a large amount of low-quality data to a small number of high-quality ones may not offer tangible benefits. On the other hand, when the high-quality data is limited in scope, diversification by a modest amount of low-quality data improves the performance. Quantitatively, for parametrizing DFT (and perhaps also quantum-chemistry ML models), caution should be taken when more than 50% of the fitting set contains questionable data, and that the average error of the full set is more than 20 kJ mol-1. One may also follow the recently proposed transferability principles to ensure diversity in the fitting set.

Assessment of the applicability of DFT methods to [Cp*Rh]-catalyzed hydrogen evolution processes

The present computational study provides a benchmark of density functional theory (DFT) methods in describing hydrogen evolution processes catalyzed by [Cp*Rh]-containing organometallic complexes. A test set was composed of 26 elementary reactions featuring chemical transformations and bonding situations essential for the field, including the emerging concept of non-innocent Cp* behavior. Reference values were obtained from a highly accurate 3/4 complete basis set and 6/7 complete PNO space extrapolated DLPNO-CCSD(T) energies. The performance of lower-level extrapolation procedures was also assessed. We considered 84 density functionals (DF) (including 13 generalized gradient approximations (GGA), nine meta-GGAs, 33 hybrids, and 29 double-hybrids) and three composite methods (HF-3c, PBEh-3c, and r2SCAN-3c), combined with different types of dispersion corrections (D3(0), D3BJ, D4, and VV10). The most accurate approach is the PBE0-DH-D3BJ (MAD of 1.36 kcal mol-1) followed by TPSS0-D3BJ (MAD of 1.60 kcal mol-1). Low-cost r2SCAN-3c composite provides a less accurate but much faster alternative (MAD of 2.39 kcal mol-1). The widely used Minnesota-family M06-L, M06, and M06-2X DFs should be avoided (MADs of 3.70, 3.94, and 4.01 kcal mol-1, respectively).

Delocalization error poisons the density-functional many-body expansion

The many-body expansion is a fragment-based approach to large-scale quantum chemistry that partitions a single monolithic calculation into manageable subsystems. This technique is increasingly being used as a basis for fitting classical force fields to electronic structure data, especially for water and aqueous ions, and for machine learning. Here, we show that the many-body expansion based on semilocal density functional theory affords wild oscillations and runaway error accumulation for ion-water interactions, typified by F-(H2O) N with N ≳ 15. We attribute these oscillations to self-interaction error in the density-functional approximation. The effect is minor or negligible in small water clusters, explaining why it has not been noticed previously, but grows to catastrophic proportion in clusters that are only moderately larger. This behavior can be counteracted with hybrid functionals but only if the fraction of exact exchange is ≳50%, whereas modern meta-generalized gradient approximations including ωB97X-V, SCAN, and SCAN0 are insufficient to eliminate divergent behavior. Other mitigation strategies including counterpoise correction, density correction (i.e., exchange-correlation functionals evaluated atop Hartree-Fock densities), and dielectric continuum boundary conditions do little to curtail the problematic oscillations. In contrast, energy-based screening to cull unimportant subsystems can successfully forestall divergent behavior. These results suggest that extreme caution is warranted when the many-body expansion is combined with density functional theory.

Direct Covalent Functionalization of H-Terminated 2D Germanane with Thiolated Molecules: Passivation and Tuning of Optoelectronic Properties

Covalent molecular functionalization allows the physicochemical properties of 2D materials to be precisely tuned and modulated on-demand. Nonetheless, research on the molecular functionalization of 2D monoelemental graphene-like materials─known as Xenes─remains scarce, being mainly restricted to a specific type of solid-state chemical reaction based on the topotactic transformation of bulkier Zintl phases. Herein, a robust and general chemical approach is reported for the direct functionalization of commercially available H-terminated 2D germanene (2D-GeH) with thiolated molecules (R-SH) via Ge-S bond formation. While the material characterization data provide direct experimental evidence of the Ge-S chemical bonding, density functional theory (DFT) calculations also predict its existence. Remarkably, the anchored thiolated molecules also favor the passivation of the 2D Xene against air oxidation, enlarging its benefits for real implementation. As a proof-of-principle, a redox-responsive molecular moiety such as 6-(ferrocenyl)hexanethiol (Fc6-SH) has been exploited to induce changes in the optoelectronic properties of the resulting 2D-GeFc6 heterostructure by simply modulating the external bias potential, making it possible to optically and electrically read out a molecular switch on 2D Xene via implanting molecular responsiveness. Remarkably, the ON/OFF ratio has been shown to be dependent on the distance between the redox-responsive Fc moiety and the 2D Xene surface through the alkyl chain length. Overall, the reported a-la-carte molecular engineering approach provides the basis toward the rapid development of stable 2D-GeR derivatives exhibiting molecule-programmable properties.

Molecular Decoration and Unconventional Double Bond Migration in Irumamycin Biosynthesis

Irumamycin (Iru) is a complex polyketide with pronounced antifungal activity produced by a type I polyketide (PKS) synthase. Iru features a unique hemiketal ring and an epoxide group, making its biosynthesis and the structural diversity of related compounds particularly intriguing. In this study, we performed a detailed analysis of the iru biosynthetic gene cluster (BGC) to uncover the mechanisms underlying Iru formation. We examined the iru PKS, including the domain architecture of individual modules and the overall spatial structure of the PKS, and uncovered discrepancies in substrate specificity and iterative chain elongation. Two potential pathways for the formation of the hemiketal ring, involving either an olefin shift or electrocyclization, were proposed and assessed using 18O-labeling experiments and reaction activation energy calculations. Based on our findings, the hemiketal ring is likely formed by PKS-assisted double bond migration and TE domain-mediated cyclization. Furthermore, putative tailoring enzymes mediating epoxide formation specific to Iru were identified. The revealed Iru biosynthetic machinery provides insight into the complex enzymatic processes involved in Iru production, including macrocycle sculpting and decoration. These mechanistic details open new avenues for a targeted architecture of novel macrolide analogs through synthetic biology and biosynthetic engineering.

The Best DFT Functional Is the Ensemble of Functionals

The development of better density functional theory (DFT) methods is one of the most active research areas, given the importance of DFT for ubiquitous molecular and materials simulations. However, this research primarily focuses on improving a specific exchange-correlation Kohn-Sham density functional. Here, a robust procedure is proposed for constructing transferable ensembles of density functionals that perform superior to any constituent individual density functional. It is shown that such ensembles built only with the density functionals predating the GMTKN55 benchmark of 2017 can reach a record-low weighted error of 1.62 kcal mol-1 on this benchmark compared to 3.08 kcal mol-1 of the best constituent density functional. The DENS24 density functional ensembles are also introduced as practical DFT methods with consistently accurate performance for various simulations at affordable cost. DENS24 ensembles are open-source and can be used for simulations online. Additionally, it is shown that the ensembles can be integrated into the SCF procedure by creating mixed DENS24 functionals, which have the same accuracy but are faster than ensembles of independent functionals.

An Acceptor-Substituted N-Heterocyclic ortho-Quinodimethane: Pushing the Boundaries of Polarization in Donor-Acceptor-Substituted Polyenes

We report the synthesis, isolation, and characterization of a stable donor-acceptor substituted ortho-quinodimethane (oQDM). This system with an imidazolidine scaffold as the donor can also be referred to as acceptor-substituted ortho-N-heterocyclic quinodimethane (oNHQ). We have examined the extent of polarization of the conjugated π-system using single-crystal X-ray diffraction, NMR and UV/vis spectroscopy, cyclic voltammetry, and DFT computations. The bond lengths in the phenyl linker do not exhibit the alternation typical of oQDMs. In addition, the 13C and 15N NMR shifts suggest significant charge separation, an interpretation supported by the diatropic ring current determined by NICSZZ(r) computations, which is characteristic of aromatic compounds. DFT calculations show that polarization is an electronic effect that is amplified by steric influences. More strikingly, the oxidation and reduction potentials of the push-pull substituted oQDM are virtually identical to those of authenticated anionic and cationic derivatives. The results therefore indicate that an aromatic zwitterionic structure represents the electronic structure more accurately than a neutral quinoidal Lewis structure, which indicates that the acceptor-substituted oNHQ is a rare example of an organic zwitterion in which the centers of charge are in conjugation. The ambiphilic reactivity of the acceptor-substituted oNHQ, which is evidenced by the dehydrogenation of ammonia borane and the addition of phenylacetylene via heterolytic C-H bond cleavage, further supports its notation as an organic zwitterion and is reminiscent of frustrated Lewis pairs (FLPs). Thus, the acceptor-substituted oNHQ can be considered to be an intramolecular carbogenic FLP in terms of its reactivity.

Thermodynamic Assessment of Sacrificial Oxidant Potential, H2O/O2 Potential, and Rate-Overpotential Relationship to Examine Catalytic Water Oxidation in Nonaqueous Solvents

The water oxidation reaction (WOR), which is pivotal to storing energy in chemical bonds, requires a catalyst to overcome its inherent kinetic barrier. In bulk solutions, sacrificial oxidants (SOs) can regenerate the catalysts to ensure that the homogeneous WOR can be operated with long-term consistent performance. To implement this strategy for organic WOR systems, we modified four common SOs with tetra-n-butylammonium ([NBu4]+)─[NBu4]2[Ce(NO3)6], [NBu4][IO4], [NBu4][HSO5], and [NBu4]2[S2O8]─and examined their chemical stability and electrochemical behaviors in various organic solvents. We also derived the organic-solvent-associated redox potential of H2O/O2 in organic media (EH2O/O2(org)) using open-circuit potential measurements of the H+/H2 redox couple and the related thermochemical cycle. Our findings indicate that the EH2O/O2(org) varies with solvent identity and can be adjusted by changing the [H2O], [acid], and [base] levels; thus, the SO should be carefully selected for WOR, because the innate redox potentials of SOs are not always higher than EH2O/O2(org) under the studied conditions. Lastly, we obtained catalyst-performance-related insights via a rate-overpotential free-energy relationship by calculating the overpotentials of previously studied WOR systems in organic media.

Calculated ionization energies, orbital eigenvalues (HOMO), and related QSAR descriptors of organic molecules: a set of 61 experimental values enables elimination of systematic errors and provides realistic error estimates

Ionization energies (IEs) of organic compounds come in different forms-adiabatic, vertical, as electrode potentials, or as orbital eigenvalues via Koopmans' theorem. They have been linked to the reactivity towards electrophiles and have been used to quantitatively describe electron transfer processes. The de novo prediction of IEs is only meaningful when an estimate of the prediction's uncertainty is included. Pulsed-field ionization (PFI) experiments have quantified adiabatic IEs with unprecedented precision. In this work, a new set of PFI-derived IEs is compiled from the literature as a benchmark for prediction methods. This set includes many common functional groups, a size range from diatomics to two aromatic rings, and IEs between 7 and 14 eV. The first-principles CCSD(T)/CBS protocol presently used reproduces these values within 0.05 eV. For adiabatic IEs and vertical IEs/orbital eigenvalues predicted using approximate density functional theory (DFT), linear regression models are proposed, so that IEs calculated using different methods can be directly compared on a physical scale. This elimination of systematic errors improves the error statistics and allows the performance of predicted IEs to be evaluated if used in quantitative structure-property or -activity relationships, as the latter implicitly correct a descriptor's bias. Owing to the structural scope of the test set, the minimum and maximum deviations from experiment should correspond to those expected for common organic molecules. Deviations from reference values found for orbital eigenvalues but also for IEs calculated explicitly with HF or semi-empirical MO methods were as large as 0.5 eV to 2.0 eV. Such large errors could also propagate into quantitative structure-property models, as shown in illustrative examples of oxidation rate constants in solution.

On the Interplay Between Force, Temperature, and Electric Fields in the Rupture Process of Mechanophores

The use of oriented external electric fields (OEEFs) shows promise as an alternative approach to chemical catalysis. The ability to target a specific bond by aligning it with a bond-weakening electric field may be beneficial in mechanochemical reactions, which use mechanical force to selectively rupture bonds. Previous computational studies have focused on a static description of molecules in OEEFs, neglecting to test the influence of thermal oscillations on molecular stability. Here, we performed ab initio molecular dynamics (AIMD) simulations based on density functional theory (DFT) to investigate the behaviour of a model mechanophore under the simultaneous influence of thermal and electric field effects. We show that the change in bond length caused by a strong electric field is largely independent of the temperature, both without and with mechanical stretching forces applied to the molecule. The amplitude of thermal oscillations increases with increasing field strength and temperature, but at low temperatures, the application of mechanical force leads to an additional increase in amplitude. Our research shows that methods for applying mechanical force and OEEFs can be safely combined and included in an AIMD simulation at both low and high temperatures, allowing researchers to computationally investigate mechanochemical reactions in realistic application scenarios.

Near equivalence of polarizability and bond order flux metrics for describing covalent bond rearrangements

Identification of the breaking point for the chemical bond is essential for our understanding of chemical reactivity. The current consensus is that a point of maximal electron delocalization along the bonding axis separates the different bonding regimes of reactants and products. This maximum transition point has been investigated previously through the total position spread and the bond-parallel components of the static polarizability tensor for describing covalent bond breaking. In this paper, we report that the first-order change of the Wiberg and Mayer bond index with respect to the reaction coordinate, the bond flux, is similarly maximized and is nearly equivalent with the bond breaking points determined by the bond-parallel polarizability. We investigate the similarites and differences between the two bonding metrics for breaking the nitrogen triple bond, twisting around the ethene double bond, and a set of prototypical reactions in the hydrogen combustion reaction network. The Wiberg-Mayer bond flux provides a simpler approach to calculating the point of bond dissociation and formation and can yield greater chemical insight through bond specific information for certain reactions where multiple bond changes are operative.

Triantaspirols A-C and Paraphaeolactone Cs from Paraphaeosphaeria sp. KT4192: Sensitivity of CP3 in Distinguishing Close NMR Signals

Hybridized spirobisnaphthalene derivatives, triantaspirols A-C (1-3) and paraphaeolactones C1 and C2 (4 and 5), were identified from the culture broth of the fungus Paraphaeosphaeria sp. KT4192. The NMR spectra of 2 and 3, as well as 4 and 5, closely resembled each other, indicating that these were pairs of diastereomers. Although this NMR spectral resemblance made it challenging to distinguish their relative configurations, detailed analysis of the electronic circular dichroism (ECD) spectra and NOE correlations allowed us to deduce them. The CP3 metric with the DFT-based NMR chemical shifts was found to distinguish configurations of diastereomers in a highly sensitive and accurate manner that DP4 could not account for because of the very close chemical shift differences in the experimental NMR spectra. The reliability of this method was assessed using 23 published examples which could not be distinguished by DP4 protocol.

Sequence-Selective Recognition of the d(GGCGCC)2 DNA Palindrome by Oligopeptide Derivatives of Mitoxantrone. Enabling for Simultaneous Targeting of the Two Guanine Bases Upstream from the Central Intercalation Site in Both Grooves and along Both Strands

The d(GGCGCC)2 palindrome is encountered in several oncogenic and retroviral sequences. In order to target it, we previously designed several oligopeptide derivatives of the mitoxantrone and ametantrone anticancer intercalators. These have two arms with a cationic side-chain in the major groove, each destined to bind along each strand O6/N7 of the two successive guanine bases (G1-G2/G1'-G2') upstream from the central anthraquinone intercalation site. We retained from a previous study (El Hage et al., 2022) a tris-intercalating molecule with two outer 9-aminoacridine (9-AA) intercalators, denoted as III. We sought enhancements in both affinity and selectivity by simultaneously targeting the minor groove of the extracyclic -NH2 groups of these bases and G4-G4' of the intercalation site. We considered derivatives of distamycin, having each pyrrole ring replaced by an imidazole to act as an in-register electron acceptor from the -NH2 group of a target guanine. We substituted the C6 and C7 carbons of anthraquinone, or the C8 and C9 ones of anthracycline, by an (imidazole-amide)3 chain. Four different derivatives of III were designed with different connectors to the anthraquinone/anthracycline and 9-AA. Polarizable molecular dynamics simulations of their complexes with a double-stranded DNA 18-mer with a central d(C GGGC GCCC G)2 palindrome sequence showed in-register minor groove binding to -NH2 of G1-G2/G1'-G2' to coexist with major groove recognition of O6/N7. Up to 12 H-bonds could be stabilized in the minor groove coexisting with four bidentate interactions of the alkyl diammonium moieties in the major groove. Since there is no mutual interference, the binding enthalpies, ΔH, contributed by each groove could add up and enable significant enhancements of the affinity constants. As was the case for their Lys precursor, these derivatives are amenable to chemical syntheses and in vitro and in vivo tests, for which the present results provide an incentive. The construction of derivatives III-A-III-D is modular. For in vitro experiments, this should enable unraveling the most important structural elements to further optimize both ΔH and TΔS and sequence selectivity and how this could translate to in vivo tests.

Fast and Robust Modeling of Lanthanide and Actinide Complexes, Biomolecules, and Molecular Crystals with the Extended GFN-FF Model

Lanthanides (Ln) and actinides (An) have recently become important tools in biomedical and materials science. However, the development of computational methods able to describe such elements in various environments has not kept up with the pace of the field. Addressing this challenge, this work introduces and showcases an extension of the GFN-FF to An alongside a reparameterization for Ln. This development fills a gap for fast computational methods that are out-of-the-box applicable to large f-element-containing systems with thousands of atoms. We discuss the reparameterization of the charge model and the covalent topology setup and showcase the model through various applications: Molecular dynamics simulations, optimization of Ln-containing biomolecules, and optimizations of several periodic structures. With the presented improvements, GFN-FF is a powerful method that routinely delivers robust and accurate geometries for large Ln/An systems with thousands of atoms.

A Critical Evaluation of the Hybrid KS DFT Functionals Based on the KS Exchange-Correlation Potentials

We have developed a critical methodology for the evaluation of the quality of hybrid exchange-correlation (XC) density functional approximations (DFAs) based on very fundamental quantities, i.e., Kohn-Sham (KS) XC potentials, self-consistent electron densities, first ionization potentials (IPs), and total energies. Since the XC potentials, the primary objects in the current study, are not directly accessible for the hybrids, we calculate them by inverting the KS electron densities. Utilizing this methodology, we tested 155 hybrid DFAs available in the LIBXC library using FCI and CCSD(T) methods as a reference. We have found that a group of functionals produces very decent XC potentials, mainly those with a large mixture of Hartree-Fock exchange. Moreover, the value of IP strongly depends on the XC potential quality. On the other hand, we show that the XC energy is dominated by functional-driven error, which in some cases leads to substantial errors in electronic densities. The study shows new directions for constructing more accurate XC functionals within the KS-DFT framework.

Cu(II) Coordination Polymers for the Selective Oxidation of Biomass-Derived Veratryl Alcohol in Green Solvents: A Sustainable Catalytic Approach

Four air-stable one-dimensional copper(II) coordination polymers (CP1-CP4) with azide linkers were synthesized using tridentate NNS and NNN ligands. Single-crystal X-ray diffraction (XRD) analysis confirmed the molecular structures of CP1, CP3, and CP4. In the presence of TEMPO, all four coordination polymers demonstrated effective catalytic activity for the selective aerobic oxidation of veratryl alcohol, a biomass model compound, under base-free conditions. CP4 exhibited the best catalytic efficiency. Oxidations were conducted at ambient temperature (40 °C) utilizing air as a sustainable oxidant. Selective oxidation of veratryl alcohol to veratraldehyde was also conducted in the presence of a catalytic amount of base (5 mol %), and enhanced reactivity was observed. The green solvents, acetone, and water, were used to maximize sustainability. The optimized reaction conditions were applied to broaden the substrate scope of various lignin model alcohols and substituted benzylic alcohols with wide electronic variability. CP4 exhibited high recyclability, consistently providing quantitative yields even after ten consecutive runs. The catalytic protocol demonstrated sustainability and environmental compatibility, as evidenced by a low E-factor (4.29) and a high Eco-scale score (90). Based on experimental evidence and theoretical calculations, a plausible catalytic cycle was proposed. Finally, the sustainability credentials of the different optimized reaction protocols were evaluated using the CHEM21 green metrics toolkit.

Exploring non-covalent interactions in excited states: beyond aromatic excimer models

Time-dependent density functional theory (TD-DFT) offers a relatively accurate and inexpensive approach for excited state calculations. However, conventional TD-DFT may suffer from the same poor description of non-covalent interactions (NCIs) which is known from ground-state DFT. In this work we present a comprehensive benchmark study of TD-DFT for excited-state NCIs. This is achieved by calculating dissociation curves for excited complexes ('exciplexes'), whose binding strength depends on excited-state NCIs including electrostatics, Pauli repulsion, charge-transfer, and London dispersion. Reference dissociation curves are calculated with the reasonably accurate wave function method SCS-CC2/CBS(3,4) which is used to benchmark a range of TD-DFT methods. Additionally, we test the effect of ground-state dispersion corrections, DFT-D3(BJ) and VV10, for exciplex binding. Overall, we find that TD-DFT methods generally under-bind exciplexes which can be explained by the missing dispersion forces. Underbinding errors reduce going up the rungs of Jacob's ladder. Further, the D3(BJ) dispersion correction is essential for good accuracy in most cases. Likewise, the VV10-type non-local kernel yields relatively low errors and has comparable performance in either its fully self-consistent implementation or as a post-SCF additive correction, but its impact is solely on ground-state energies and not on excitation energies. From our analysis, the most robust TD-DFT methods for exciplexes with localised excitations in their equilibrium and non-equilibrium geometries are the double hybrids B2GP-PLYP-D3(BJ) and B2PLYP-D3(BJ). Their range-separated versions ωB2(GP-)PLYP-D3(BJ) or the spin-opposite scaled, range-separated double hybrid SOS-ωB88PP86 can be recommended when charge transfer plays a role in the excitations. We also identify the need for a state-specific dispersion correction as the next step for improved TD-DFT performance.

Graphene in Water is Hardly Ever Neutral

Graphene in water is electrically charged in most conditions. The level of charge can be large enough to stabilize single (or few) layer graphene colloidal dispersions in water, without the need of using any other additive. In this work, potentiometric titration, isothermal titration calorimetry, electrokinetic measurements, Density Functional Theory calculations, Raman Spectroscopy, and direct force measurements using Atomic Force Microscopy to investigate this charge and explore its origin are combined. The body of data collected suggests that this charge is a consequence of the interaction between water ions (hydroxide and hydronium) and graphene, and can be conveniently tuned (in magnitude and sign) by changing the pH of water.

Enhanced ClNO2 Formation at the Interface of Sea-Salt Aerosol

The reactive uptake of N2O5 on sea-spray aerosol plays a key role in regulating the NOx concentration in the troposphere. Despite numerous field and laboratory studies, a microscopic understanding of its heterogeneous reactivity remains unclear. Here, we use molecular simulation and theory to elucidate the chlorination of N2O5 to form ClNO2, the primary reactive channel within sea-spray aerosol. We find that the formation of ClNO2 is markedly enhanced at the air-water interface due to the stabilization of the charge-delocalized transition state, as evident from the formulation of bimolecular rate theory in heterogeneous environments. We explore the consequences of the enhanced interfacial reactivity in the uptake of N2O5 using numerical solutions of molecular reaction-diffusion equations as well as their analytical approximations. Our results suggest that the current interpretation of aerosol branching ratios needs to be revisited.

Toward Reliable Conformational Energies of Amino Acids and Dipeptides─The DipCONFS Benchmark and DipCONL Datasets

Simulating peptides and proteins is becoming increasingly important, leading to a growing need for efficient computational methods. These are typically semiempirical quantum mechanical (SQM) methods, force fields (FFs), or machine-learned interatomic potentials (MLIPs), all of which require a large amount of accurate data for robust training and evaluation. To assess potential reference methods and complement the available data, we introduce two sets, DipCONFL and DipCONFS, which cover large parts of the conformational space of 17 amino acids and their 289 possible dipeptides in aqueous solution. The conformers were selected from the exhaustive PeptideCS dataset by Andris et al. [ J. Phys. Chem. B 2022, 126, 5949-5958]. The structures, originally generated with GFN2-xTB, were reoptimized using the accurate r2SCAN-3c density functional theory (DFT) composite method including the implicit CPCM water solvation model. The DipCONFS benchmark set contains 918 conformers and is one of the largest sets with highly accurate coupled cluster conformational energies so far. It is employed to evaluate various DFT and wave function theory (WFT) methods, especially regarding whether they are accurate enough to be used as reliable reference methods for larger datasets intended for training and testing more approximated SQM, FF, and MLIP methods. The results reveal that the originally provided BP86-D3(BJ)/DGauss-DZVP conformational energies are not sufficiently accurate. Among the DFT methods tested as an alternative reference level, the revDSD-PBEP86-D4 double hybrid performs best with a mean absolute error (MAD) of 0.2 kcal mol-1 compared with the PNO-LCCSD(T)-F12b reference. The very efficient r2SCAN-3c composite method also shows excellent results, with an MAD of 0.3 kcal mol-1, similar to the best-tested hybrid ωB97M-D4. With these findings, we compiled the large DipCONFL set, which includes over 29,000 realistic conformers in solution with reasonably accurate r2SCAN-3c reference conformational energies, gradients, and further properties potentially relevant for training MLIP methods. This set, also in comparison to DipCONFS, is used to assess the performance of various SQM, FF, and MLIP methods robustly and can complement training sets for those.

van der Waals Radii of Free and Bonded Atoms from Hydrogen (Z = 1) to Oganesson (Z = 118)

Reliable numerical values of van der Waals (vdW) radii are required for constructing empirical force fields, vdW-inclusive density functional, and quantum-chemical methods, as well as for implicit solvent models. However, multiple definitions exist for vdW radii, involving either equilibrium or the closest contact distances between free or bonded atoms within molecules or crystals. For the paradigmatic case of the hydrogen atom, its reported vdW radius fluctuates between 2.15 and 3.70 Bohr depending on the definition, leading to a high uncertainty in calculations and different conceptual interpretations of noncovalent interactions. In this work, we systematically review different definitions and methodologies to establish the free and bonded vdW radii for hydrogen, based on equilibrium vdW distances in noncovalently bonded molecules, enveloping electron density cutoffs, noncovalent positron bonds in hydrogen anion dimer, vacuum virtual photon cloud caused by the hydrogen atom, and atomic dipole polarizability. By doing so, we show that the vdW radius of the free hydrogen atom is 3.16 ± 0.06 Bohr. By employing the most general and elegant definition of atomic vdW radius as a function of the atomic polarizability, we tabulate consistent values of vdW radii for all atoms in the periodic table up to Z = 118.

Organic Chimeras based on Selenosugars, Steroids, and Fullerenes as Potential Inhibitors of the β-amyloid Peptide Aggregation

The aggregation of β-amyloid peptide (Aβ) is associated with neurodegenerative diseases such as Alzheimer's disease (AD). Several therapies aimed at reducing the aggregation of this peptide have emerged as potential strategies for the treatment of AD. This paper describes the design and preparation of new hybrid molecules based on steroids, selenosugars, and [60]fullerene as potential inhibitors of Aβ oligomerization. These moieties were selected based on their antioxidant properties and possible areas of interaction with the Aβ. Cyclopropanations between C60 and malonates bearing different steroid and selenosugar moieties using the Bingel-Hirsch protocol have enabled the synthesis of functionalized molecular hybrids. The obtained derivatives were characterized by physical and spectroscopic techniques. Theoretical calculations for all the selenium compounds were performed using the density functional theory DFT/B3LYP-D3(BJ)/6-311G(2d,p) predicting the most stable conformations of the synthesized derivatives. Relevant geometrical parameters were investigated to relate the stereochemical behavior and the spectroscopic data obtained. The affinity of the compounds for Aβ-peptide was estimated by molecular docking simulation, which predicted an increase in affinity and interactions for Aβ for the hybrids containing the C60 core. In addition, parameters such as lipophilicity, polar surface area, and dipole moment were calculated to predict their potential interaction with membrane cells.

Toward Efficient and Unified Treatment of Static and Dynamic Correlations in Generalized Kohn-Sham Density Functional Theory

Accurate description of the static correlation poses a persistent challenge in electronic structure theory, particularly when it has to be concurrently considered with the dynamic correlation. We develop here a method in the generalized Kohn-Sham density functional theory (DFT) framework, named R-xDH7-SCC15, which achieves an unprecedented accuracy in capturing the static correlation, while maintaining a good description of the dynamic correlation on par with the state-of-the-art DFT and wave function theory methods, all grounded in the same single-reference black-box methodology. Central to R-xDH7-SCC15 is a general-purpose static correlation correction (SCC) model applied to the renormalized XYG3-type doubly hybrid method (R-xDH7). The SCC model development involves a hybrid machine learning strategy that integrates symbolic regression with nonlinear parameter optimization, aiming to achieve a balance between generalization capability, numerical accuracy, and interpretability. Extensive benchmark studies confirm the robustness and broad applicability of R-xDH7-SCC15 across a diverse array of main-group chemical scenarios. Notably, it displays exceptional aptitude in accurately characterizing intricate reaction kinetics and dynamic processes in regions distant from equilibrium, where the influence of static correlation is most profound. Its capability to consistently and efficiently predict the whole energy profiles, activation barriers, and reaction pathways within a user-friendly "black-box" framework represents an important advance in the field of electronic structure theory.

Dynamic and solid-state behaviour of bromoisotrianglimine

Solid-state materials formed from discrete imine macrocycles have potential in industrial separations, but dynamic behaviour during both synthesis and crystallisation makes them challenging to exploit. Here, we explore opportunities for structural control by investigating the dynamic nature of a C-5 brominated isotrianglimine in solution and under crystallisation conditions. In solution, the equilibrium between the [3 + 3] and the less reported [2 + 2] macrocycle was investigated, and both macrocycles were fully characterised. Solvent templating during crystallisation was used to form new packing motifs for the [3 + 3] macrocycle and a previously unreported [4 + 4] macrocycle. Finally, chiral self-sorting was used to demonstrate how crystallisation conditions can not only influence packing arrangements but also shift the macrocycle equilibrium to yield new structures. This work thus exemplifies three strategies for exploiting dynamic behaviour to form isotrianglimine materials, and highlights the importance of understanding the dynamic behaviour of a system when designing and crystallising functional materials formed using dynamic covalent chemistry.

Getting the intermolecular forces correct: introducing the ASTA strategy for a water model

Having a force field for water providing good bulk properties is paramount for modern studies of most biological systems. Some of the most common three-site force fields are TIP3, SPC/ε or OPC3, providing a decent range of bulk properties. That does not mean though, that they have realistic inter-atomic forces. These force fields have been parameterized with a top-down approach, meaning, by fitting the force field parameters to the experimental bulk properties. This approach has been the governing strategy also for many variants of four- and more-site models. We test a bottom-up approach, in which the force field is parameterized by optimizing the non-bonded inter-atomic forces. Our philosophy is that correct inter-atomic forces lead to correct geometrical and dynamical properties. The first system we try to optimize with the accurately system tailored atomic (ASTA) approach is water, but we aim to eventually probe other systems in the future as well. We applied our ASTA strategy to find a good set of parameters providing accurate bulk properties for the simple three-site force field forms, and also for AMOEBA, a more detailed and polarizable force field. Even though our bottom-up approach did not provide satisfactory results for the simple three-site force fields (with fixed charges), for the case of the AMOEBA force field it led to a modification of the original strategy, giving very good intra- and inter-molecular forces, as compared to accurate quantum chemically calculated reference forces. At the same time, important bulk properties, in this study restricted to the density and diffusion, were accurately reproduced with respect to the experimental values.

Combined Multireference-Multiscale Approach to the Description of Photosynthetic Reaction Centers

A first-principles description of the primary photochemical processes that drive photosynthesis and sustain life on our planet remains one of the grand challenges of modern science. Recent research established that explicit incorporation of protein electrostatics in excited-state calculations of photosynthetic pigments, achieved for example with quantum-mechanics/molecular-mechanics (QM/MM) approaches, is essential for a meaningful description of the properties and function of pigment-protein complexes. Although time-dependent density functional theory has been used productively so far in QM/MM approaches for the study of such systems, this methodology has limitations. Here we pursue for the first time a QM/MM description of the reaction center in the principal enzyme of oxygenic photosynthesis, Photosystem II, using multireference wave function theory for the high-level QM region. We identify best practices and establish guidelines regarding the rational choice of active space and appropriate state-averaging for the efficient and reliable use of complete active space self-consistent field (CASSCF) and the N-electron valence state perturbation theory (NEVPT2) in the prediction of low-lying excited states of chlorophyll and pheophytin pigments. Given that the Gouterman orbitals are inadequate as a minimal active space, we define specific minimal and extended active spaces for the NEVPT2 description of electronic states that fall within the Q and B bands. Subsequently, we apply our multireference-QM/MM protocol to the description of all pigments in the reaction center of Photosystem II. The calculations reproduce the electrochromic shifts induced by the protein matrix and the ordering of site energies consistent with the identity of the primary donor (ChlD1) and the experimentally known asymmetric and directional electron transfer. The optimized protocol sets the stage for future multireference treatments of multiple pigments, and hence for multireference studies of charge separation, while it is transferable to the study of any photoactive embedded tetrapyrrole system.

Machine learning for accuracy in density functional approximations

Machine learning techniques have found their way into computational chemistry as indispensable tools to accelerate atomistic simulations and materials design. In addition, machine learning approaches hold the potential to boost the predictive power of computationally efficient electronic structure methods, such as density functional theory, to chemical accuracy and to correct for fundamental errors in density functional approaches. Here, recent progress in applying machine learning to improve the accuracy of density functional and related approximations is reviewed. Promises and challenges in devising machine learning models transferable between different chemistries and materials classes are discussed with the help of examples applying promising models to systems far outside their training sets.

A nucleophilic beryllyl complex via metathesis at [Be-Be]2

Owing to its high toxicity, the chemistry of element number four, beryllium, is poorly understood. However, as the lightest elements provide the basis for fundamental models of chemical bonding, there is a need for greater insight into the properties of beryllium. In this context, the chemistry of the homo-elemental Be-Be bond is of fundamental interest. Here the ligand metathesis chemistry of diberyllocene (1; CpBeBeCp)-a stable complex with a Be-Be bond-has been investigated. These studies yield two complexes with Be-Be bonds: Cp*BeBeCp (2) and [K{(HCDippN)2BO}2]BeBeCp (3; Dipp = 2,6-diisopropylphenyl). Quantum chemical calculations indicate that the Be-Be bond in 3 is polarized to such an extent that the complex could be formulated as a mixed-oxidation state Be0/BeII complex. Correspondingly, it is demonstrated that 3 can transfer the 'beryllyl' anion, [BeCp]-, to an organic substrate, by analogy with the reactivity of sp2-sp3 diboranes. Indeed, this work reveals striking similarities between the homo-elemental bonding linkages of beryllium and boron, despite the respective metallic and non-metallic natures of these elements.

Evaluating Noncovalent Interactions in Halogenated Molecules with Double-Hybrid Functionals and a Dedicated Small Basis Set

We present here an extension of our recently developed PBE-QIDH/DH-SVPD basis set to halogen atoms, with the aim of obtaining, for weakly interacting halogenated molecules, interaction energies close to those provided by a large basis set (def2-TZVPP) coupled to empirical dispersion potential. The core of our approach is the split-valence basis set, DH-SVPD, that has been developed for F, Cl, Br, and I atoms using a self-consistent formula, containing only energy terms computed for dimers and the corresponding monomers at the same level of theory. The basis set developed considering four systems, one for each halogen atoms, has been then tested on the X40, X4 × 10 benchmarks as well as on other two, less standard, data sets. Finally, a large system (380 atoms) has been also considered as a "crash" test. Our results show that the simple and nonempirical PBE-QIDH/DH-SVPD approach is able to provide accurate results for interaction energies of all the considered systems and can thus be considered as a cheaper alternative to DH functionals paired with empirical dispersion corrections and a large basis set of triple-ζ quality.

A Vacuum Ultraviolet Photoionization Mass Spectrometry and Density Functional Calculation Study of Formic Acid-Water Clusters

The interaction between formic acid (FA) and water (W) holds significant importance in various chemical processes. Our study combines vacuum-ultraviolet photoionization mass spectrometry with density functional calculations to investigate formic acid water clusters generated in supersonic molecular beams. The mass spectra obtained reveal the formation of protonated clusters as the major product. Enhanced intensities are observed in the mass spectra for a number of clusters holding the following composition, FA1W5H+, FA2W4H+, FA3W3H+, FA4W2H+, FA5W1H+ and FA6W2H+ compared to their neighbors with one less or one more water component. Our calculations shed light on these potentially stable structures, highlighting cyclic arrangements with molecules enclosed within the ring as the most stable structures, and demonstrate a decrease in the stability upon the addition of a water molecule. Comparing experimental appearance energies with calculated ionization energies suggests that fragmentation can occur from clusters of various sizes.

Influence of Electronic Correlations on Electron-Phonon Interactions of Molecular Systems with the GW and Coupled Cluster Methods

Electron-phonon interactions are of great importance to a variety of physical phenomena, and their accurate description is an important goal for first-principles calculations. Isolated examples of materials and molecular systems have emerged where electron-phonon coupling is enhanced over density functional theory (DFT) when using the Green's-function-based ab initio GW method, which provides a more accurate description of electronic correlations. It is, however, unclear how general this enhancement is and how employing high-end quantum chemistry methods, which further improve the description of electronic correlations, might further alter electron-phonon interactions over GW or DFT. Here, we address these questions by computing the renormalization of the highest occupied molecular orbital energies of Thiel's set of organic molecules by harmonic vibrations using DFT, GW, and equation-of-motion coupled-cluster calculations. We find that, depending on the amount of exact exchange included in the DFT starting point, GW can increase the magnitude of the electron-phonon coupling across Thiel's set of molecules by an average factor of 1.1-1.8 compared to the underlying DFT, while equation-of-motion coupled-cluster leads to an increase of 1.4-2. The electron-phonon coupling predicted with the ab initio GW method is generally in much closer agreement to coupled cluster values compared to DFT, establishing GW as a promising route for accurately computing electron-phonon phenomena in molecules and beyond at a much lower computational cost than higher-end quantum chemistry techniques.

Quantum chemical calculations of nitrosamine activation and deactivation pathways for carcinogenicity risk assessment

N-nitrosamines and nitrosamine drug substance related impurities (NDSRIs) became a critical topic for the development and safety of small molecule medicines following the withdrawal of various pharmaceutical products from the market. To assess the mutagenic and carcinogenic potential of different N-nitrosamines lacking robust carcinogenicity data, several approaches are in use including the published carcinogenic potency categorization approach (CPCA), the Enhanced Ames Test (EAT), in vivo mutagenicity studies as well as read-across to analogue molecules with robust carcinogenicity data. We employ quantum chemical calculations as a pivotal tool providing insights into the likelihood of reactive ion formation and subsequent DNA alkylation for a selection of molecules including e.g., carcinogenic N-nitrosopiperazine (NPZ), N-nitrosopiperidine (NPIP), together with N-nitrosodimethylamine (NDMA) as well as non-carcinogenic N-nitrosomethyl-tert-butylamine (NTBA) and bis (butan-2-yl) (nitros)amine (BBNA). In addition, a series of nitroso-methylaminopyridines is compared side-by-side. We draw comparisons between calculated reaction profiles for structures representing motifs common to NDSRIs and those of confirmed carcinogenic and non-carcinogenic molecules with in vivo data from cancer bioassays. Furthermore, our approach enables insights into reactivity and relative stability of intermediate species that can be formed upon activation of several nitrosamines. Most notably, we reveal consistent differences between the free energy profiles of carcinogenic and non-carcinogenic molecules. For the former, the intermediate diazonium ions mostly react, kinetically controlled, to the more stable DNA adducts and less to the water adducts via transition-states of similar heights. Non-carcinogenic molecules yield stable carbocations as intermediates that, thermodynamically controlled, more likely form the statistically preferred water adducts. In conclusion, our data confirm that quantum chemical calculations can contribute to a weight of evidence approach for the risk assessment of nitrosamines.

Pushing the Limits of Organometallic Redox Chemistry with an Isolable Mn(-I) Dianion

We report an incredibly reducing and redox-active Mn-I dianion, [Mn(CO)3(Ph2B(tBuNHC)2)]2- (NHC = N-heterocyclic carbene), furnished via 2e- reduction of the parent 16e- MnI complex with Na0 or K0. Cyclic voltammograms show a Mn0/-I redox couple at -3.13 V vs Fc+/0 in tetrahydrofuran (THF), -3.06 V in 1,2-dimethoxyethane, and -2.85 V in acetonitrile. The diamagnetic Mn-I dianion is stable in solution and solid-state at room temperature, tolerating a wide range of countercations ([M(2.2.2)crypt]+, [M(18-crown-6)]+, [nBu4N]+; M = Na, K). Countercation identity does not significantly alter 13C NMR spectral signatures with [nBu4N]+ and Na+, suggesting minimal ion pairing in solution. IR spectroscopy reveals a significant decrease in CO stretching frequencies from MnI to Mn-I (ca. 240 cm-1), consistent with a drastic increase in electron density at Mn. State-of-the-art DFT calculations are in excellent agreement with the observed IR spectral data. Moreover, the Mn-I dianion behaves as a chemical reductant, smoothly releasing 1e- or 2e- to regenerate the oxidized Mn0 or MnI species in solution. The reducing potential of [Mn(CO)3(Ph2B(tBuNHC)2)]2- surpasses the naphthalenide anion in THF (-3.09 V) and represents one of the strongest isolable chemical redox agents.