Density functional theory with London dispersion corrections

An electrochemical sensing potential of cobalt oxide nanoparticles towards citric acid integrated with computational approach in food and biological media

Although citric acid (CA) has antioxidant, antibacterial, and acidulating properties, chronic ingestion of CA can cause urolithiasis, hypocalcemia, and duodenal cancer, emphasizing the need for early detection. There are very few documented electrochemical-based sensing methods for CA detection due to the challenging behavior of electrode fouling caused by reactive oxidation products. In this study, a novel, non-enzymatic, and economical electrochemical sensor based on cobalt oxide nanoparticles (CoOxNPs) is successfully reported for detection CA. The CoOxNPs were synthesized through a simple thermal decomposition method and characterized by SEM, FT-IR, EDX, and XRD techniques. The proposed sensing platform was optimized by various parameters, including pH (7.0), time (15 min), and concentration of nanoparticles (100 mM) etc. In a linear range of 0.05-2500 μM, a low detection limit (LOD) of 0.13 μM was achieved. Theoretical calculations (ΔRT), confirmed hydrogen bonding and electrostatic interactions between CoOxNPs and CA. The detection method exhibited high selectivity in real media like food and biological samples, with good recovery values when compared favorably to the HPLC method. To facilitate effective on-site investigation, such a sensing platform can be assembled into a portable device.

High resolution far-infrared synchrotron spectroscopy of 2-furfural conformers: Fundamental and hot bands

In the continuity of a previous jet-cooled rovibrational study of trans and cis conformers of 2-furfural in the mid-infrared region (700-1750 cm-1) [Chawananon et al., Molecules 28 (10), 4165 (2023)], the present work investigates the far-infrared spectroscopy of 2-furfural using a long path absorption cell coupled to a high-resolution Fourier transform spectrometer and synchrotron radiation at the AILES beamline of the SOLEIL synchrotron. Guided by anharmonic calculations, vibrational energy levels and excited-state rotational constants are sufficiently predictive for a complete assignment of all fundamental and combination bands up to 700 cm-1, as well as the rovibrational analysis of 4 (1) low-frequency modes of trans-(cis-)2-furfural. A global rovibrational simulation, including far-infrared rovibrational lines and microwave and millimeter-wave rotational lines assigned in a previous study [Motiyenko et al., J. Mol. Spectrosc., 244, 9 (2007)] provides a reliable set of ground- and excited-state rotational parameters involving ring torsion, bending, and ring puckering modes of 2-furfural. In a second step, a rovibrational analysis of several hot band sequences, mainly involving the lowest frequency ring CHO torsion mode, is carried out. Reliable values of some anharmonic coefficients are obtained experimentally and could serve as a benchmark for validating advanced anharmonic calculations related to these large amplitude motions of flexible molecules.

Computational Insights into Amide Bond Formation Catalyzed by the Condensation Domain of Nonribosomal Peptide Synthetases

Nonribosomal peptide synthetases (NRPSs) are important enzymes that synthesize an array of nongenetically encoded peptides. The latter have diverse physicochemical properties and roles. NRPSs are modular enzymes in which, for example, the condensation (C-) domain catalyzes the formation of amide bonds. The NRPS tyrocidine synthetase from Brevibacillus brevis is responsible for synthesizing the cyclic-peptide antibiotic tyrocidine. The first step is formation of an amide bond between a proline and phenylalanine which is catalyzed by a C-domain. In this study, a multiscale computational approach (molecular dynamics and QM/MM) has been used to investigate substrate binding and catalytic mechanism of the C-domain of tyrocidine synthetase. Overall, the mechanism is found to proceed through three exergonic steps in which an active site Histidine, His222, acts as a base and acid. First, His222 acts as a base to facilitate nucleophilic attack of the prolyl nitrogen at the phenylalanyl's carbonyl carbon. This is also the rate-limiting step with a free energy barrier of 38.8 kJ mol-1. The second step is collapse of the resulting tetrahedral intermediate with cleavage of the S-C bond between the phenylalanyl and its Ppant arm, along with formation of the above amide bond. Meanwhile, the now protonated His222 imidazole has rotated toward the newly formed thiolate of the Ppant arm. In the final step, His222 acts as an acid, protonating the thiolate and regenerating a neutral His222. The overall mechanism is found to be exergonic with the final product complex being 46.3 kJ mol-1 lower in energy than the initial reactant complex.

Temperature-dependent Raman spectroscopy and thermodynamic behaviors of energetic NTO crystal

The vibrational and thermodynamic properties of energetic materials (EMs) are critical to understand their structure responses at finite temperature. In this work, the zero-point energy and temperature effects are incorporated into dispersion-corrected density functional theory to improve the calculated accuracy for vibrational responses and thermodynamic behaviors of 3-nitro-1,2,4-triazole-5-one (NTO). Based on temperature-dependent Raman spectroscopy, the emergence and disappearance of new peaks as well as discontinuous Raman shifts indicate the distinct changes of molecular configuration and intermolecular interactions within the temperature of 250-350 K. From Hirshfeld surface and structure analysis, the subtle changes of intermolecular hydrogen bonds (HBs) related with the shrinkage of thermal expansion coefficient, are treated as an essential step of a potential structural transformation of NTO. Moreover, the calculated heat capacity, entropy and bulk moduli could reflect the softening behavior of NTO and further enrich the thermodynamic data set of EMs. These results demonstrate the evolution of NTO molecules controlled by non-covalent interactions and provide vital insights into the thermodynamic behaviors at finite temperature.

Navigating through the coordination preferences of heavy alkaline earth metals: Laying the foundations for 223Ra- and 131/135mBa-based targeted alpha therapy and theranostics of cancer

The clinical success of [223Ra]RaCl2 (Xofigo®) for the palliative treatment of bone metastases in patients with prostate cancer has highlighted the therapeutic potential of α-particle emission. Expanding the applicability of radium-223 in Targeted Alpha Therapy of non-osseous tumors is followed up with significant interest, as it holds the potential to unveil novel treatment options in the comprehensive management of cancer. Moreover, the use of barium radionuclides, like barium-131 and -135m, is still unfamiliar in nuclear medicine applications, although they can be considered as radium-223 surrogates for imaging purposes. Enabling these applications requires the establishment of chelators able to form stable complexes with radium and barium radionuclides. Until now, only a limited number of ligands have been suggested and these molecules have been primarily inspired by existing structures known for their ability to complex large metal cations. However, a systematic inspection of chelators specifically tailored to Ra2+ and Ba2+ has yet to be conducted. This work delves into a comprehensive investigation of a series of small organic ligands, aiming to unveil the coordination preferences of both radium-223 and barium-131/135m. Electronic binding energies of both metal cations to each ligand were theoretically computed via Density Functional Theory calculations (COSMO-ZORA-PBE-D3/TZ2P), while thermodynamic stability constants were experimentally determined for Ba2+-ligand complexes by potentiometry, NMR and UV-Vis spectroscopies. The outcomes revealed malonate, 2-hydroxypyridine 1-oxide and picolinate as the most favorable building blocks to design multidentate chelators. These findings serve as foundation guidelines, propelling the development of cutting-edge radium-223- and barium-131/135m-based radiopharmaceuticals for Targeted Alpha Therapy and theranostics of cancer.

Sodium Triflate Water-in-Salt Electrolytes in Advanced Battery Applications: A First-Principles-Based Molecular Dynamics Study

Offering a compelling combination of safety and cost-effectiveness, water-in-salt (WiS) electrolytes have emerged as promising frontiers in energy storage technology. Still, there is a strong demand for research and development efforts to make these electrolytes ripe for commercialization. Here, we present a first-principles-based molecular dynamics (MD) study addressing in detail the properties of a sodium triflate WiS electrolyte for Na-ion batteries. We have developed a workflow based on a machine learning (ML) potential derived from ab initio MD simulations. As ML potentials are typically restricted to the interpolation of the data points of the training set and have hardly any predictive properties, we subsequently optimize a classical force field based on physics principles to ensure broad applicability and high performance. Performing and analyzing detailed MD simulations, we identify several very promising properties of the sodium triflate as a WiS electrolyte but also indicate some potential stability challenges associated with its use as a battery electrolyte.

Azophosphines: Synthesis, Structure and Coordination Chemistry

The conceptual replacement of nitrogen with phosphorus in common organic functional groups unlocks new properties and reactivity. The phosphorus-containing analogues of triazenes are underexplored but offer great potential as flexible and small bite-angle ligands. This manuscript explores the synthesis and characterisation of a family of air-stable azophosphine-borane complexes, and their subsequent deprotection to the free azophosphines. These compounds are structurally characterised, both experimentally and computationally, and highlight the availability of the phosphorus lone pair for coordination. This is confirmed by demonstrating that neutral azophosphines can act as ligands in Ru complexes, and can coordinate as monodentate or bidentate ligands in a controlled manner, in contrast to their nitrogen analogues.

Deciphering in-situ surface reconstruction in two-dimensional CdPS3 nanosheets for efficient biomass hydrogenation

Steering on the intrinsic active site of an electrode material is essential for efficient electrochemical biomass upgrading to valuable chemicals with high selectivity. Herein, we show that an in-situ surface reconstruction of a two-dimensional layered CdPS3 nanosheet electrocatalyst, triggered by electrolyte, facilitates efficient 5-hydroxymethylfurfural (HMF) hydrogenation to 2,5-bis(hydroxymethyl)furan (BHMF) under ambient condition. The in-situ Raman spectroscopy and comprehensive post-mortem catalyst characterizations evidence the construction of a surface-bounded CdS layer on CdPS3 to form CdPS3/CdS heterostructure. This electrocatalyst demonstrates promising catalytic activity, achieving a Faradaic efficiency for BHMF reaching 91.3 ± 2.3 % and a yield of 4.96 ± 0.16 mg/h at - 0.7 V versus reversible hydrogen electrode. Density functional theory calculations reveal that the in-situ generated CdPS3/CdS interface plays a pivotal role in optimizing the adsorption of HMF* and H* intermediate, thus facilitating the HMF hydrogenation process. Furthermore, the reconstructed CdPS3/CdS heterostructure cathode, when coupled with MnCo2O4.5 anode, enables simultaneous BHMF and formate synthesis from HMF and glycerol substrates with high efficiency.

External Electric Field Control of Exciton Motion in Porphyrin-Based Metal Organic Frameworks

Porphyrins are excellent light-harvesting complexes. Presently they are unsuitable for photovoltaic applications, as their excellent light absorbance is compensated to a large extent by their poor transport properties, where most excitons are lost by recombination. Arranging porphyrins in regular, strongly bound, lattices of surface-anchored metal-organic frameworks (PP-SURMOFs) may facilitate charge carrier dissociation, but does not significantly enhance the conductive properties. In most cases, photogenerated excitons traverse undirected, Brownian motion through a hopping process, resulting in a substantial diffusion length to reach electrodes, leading to significant exciton loss through recombination. Here, we propose to guide exciton diffusion indirectly by an external electric field. We show that electric fields, even as strong as 1 V nm-1, do not affect the HOMO-LUMO gap of the porphyrins. However, fields of 0.1 V nm-1 and even less demonstrate a notable Stark effect, with slight band gap reductions, for some PP-SURMOFs. When applied as an electric field gradient, for instance, via the substrate, it creates a unidirectional hopping pathway for the excitons. Consequently, we expect a significant reduction of exciton diffusion length leading to increased utilization of photogenerated excitons as they reach the electrodes. This strategy holds promise for integrating photoactive molecules in photovoltaic and photocatalytic applications.

Light Induced Cobalt(III) Carbene Radical Formation from Dimethyl Malonate As Carbene Precursor

Radical-type carbene transfer catalysis is an efficient method for the direct functionalization of C-H and C=C bonds. However, carbene radical complexes are currently formed via high-energy carbene precursors, such as diazo compounds or iodonium ylides. Many of these carbene precursors require additional synthetic steps, have an explosive nature, or generate halogenated waste. Consequently, the utilization of carbene radical catalysis is limited by specific carbene precursors that access the carbene radical intermediate. In this study, we generate a cobalt(III) carbene radical complex from dimethyl malonate, which is commercially available and bench-stable. EPR and NMR spectroscopy were used to identify the intermediates and showed that the cobalt(III) carbene radical complex is formed upon light irradiation. In the presence of styrene, carbene transfer occurred, forming cyclopropane as the product. With this photochemical method, we demonstrate that dimethyl malonate can be used as an alternative carbene precursor in the formation of a cobalt(III) carbene radical complex.

Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment

This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.

Theoretical Photoelectron Spectroscopy of Metal-Metal Quintuple Bonds: Relativity-Driven Reordering of Frontier Orbitals

A recent reinvestigation of the gas-phase photoelectron spectra of Group 6 metal-metal quadruple-bonded complexes with scalar-relativistic DFT calculations showed that common exchange-correlation functionals reproduce the lowest ionization potentials in a semiquantitative manner. The finding encouraged us to undertake a DFT study of metal-metal quintuple bonds in a set of bisamidinato complexes with the formula MI 2[HC(NR)2]2 (M = Cr, Mo, W; R = H, Ph, 2,6-iPr2C6H3) and idealized D 2h symmetry. Scalar-relativistic OLYP/STO-TZ2P calculations indicated significant shifts in valence orbital energies among the three metals, which translate to lower first ionization potentials, higher electron affinities, and lower HOMO-LUMO gaps for the W complexes relative to their Cr and Mo counterparts. These differences are largely attributable to substantially larger relativistic effects in the case of tungsten relative to those of its lighter congeners.

In Silico Chemical Experiments in the Age of AI: From Quantum Chemistry to Machine Learning and Back

Computational chemistry is an indispensable tool for understanding molecules and predicting chemical properties. However, traditional computational methods face significant challenges due to the difficulty of solving the Schrödinger equations and the increasing computational cost with the size of the molecular system. In response, there has been a surge of interest in leveraging artificial intelligence (AI) and machine learning (ML) techniques to in silico experiments. Integrating AI and ML into computational chemistry increases the scalability and speed of the exploration of chemical space. However, challenges remain, particularly regarding the reproducibility and transferability of ML models. This review highlights the evolution of ML in learning from, complementing, or replacing traditional computational chemistry for energy and property predictions. Starting from models trained entirely on numerical data, a journey set forth toward the ideal model incorporating or learning the physical laws of quantum mechanics. This paper also reviews existing computational methods and ML models and their intertwining, outlines a roadmap for future research, and identifies areas for improvement and innovation. Ultimately, the goal is to develop AI architectures capable of predicting accurate and transferable solutions to the Schrödinger equation, thereby revolutionizing in silico experiments within chemistry and materials science.

π-π+ Interaction Promoting the Absorption of Electroactive Chiral Selectors into the Cavity of Conductive Covalent Organic Framework for Enantioselective Sensing of Electrochemically Silent Molecules

To date, achieving enantioselective electroanalysis for electrochemically silent chiral molecules is still highly desired. Here, an ionic covalent organic framework (COF) consisting of the pyridinium cation was derived from the tripyridinium Zincke salt and 1,4-phenylenediamine in a one-pot reaction. The electrochemical measurements revealed that the ionic backbone contributed to the electron transfer with a low charge transfer resistance. Besides, the π-π+ interaction between the pyridinium cation and ferrocenyl unit can promote the absorption of electroactive chiral ferrocenyl reagents into the hole of COF, so as to afford the electrochemical signals by themselves, replacing the testing enantiomers. As a result, the electroactive complex used as an electrochemical platform was highly effective at enantiomerically recognizing amino alcohols (prolinol, valinol, leucinol, and alaninol) and amino acids (methionine, serine, and penicillamine), giving the ratios of current intensity between l- and d-enantiomers in the range of 1.46-1.72. Moreover, the density functional theory calculations determined the possible intermolecular interactions between the testing enantiomers and chiral selector: namely, hydrogen bonds and electrostatic attractions. Overall, the present work offers an effective strategy to enlarge the electrochemical scope for chiral recognition based on electroactive chiral COFs.

Inverse Hypercorroles

Ground-state and time-dependent density functional theory (TDDFT) calculations with the long-range-corrected, Coulomb-attenuating CAMY-B3LYP exchange-correlation functional and large, all-electron STO-TZ2P basis sets have been used to examine the potential "inverse hypercorrole" character of meso-p-nitrophenyl-appended dicyanidocobalt(III) corrole dianions. The effect is most dramatic for 5,15-bis(p-nitrophenyl) derivatives, where it manifests itself in intense NIR absorptions. The 10-aryl groups in these complexes play a modulatory role, as evinced by experimental UV-visible spectroscopic and electrochemical data for a series of 5,15-bis(p-nitrophenyl) dicyanidocobalt(III) corroles. TDDFT (CAMY-B3LYP) calculations ascribe these features clearly to a transition from the corrole's a2u-like HOMO (retaining the D4h irrep used for metalloporphyrins) to a nitrophenyl-based LUMO. The outward nature of this transition contrasts with the usual phenyl-to-macrocycle direction of charge transfer transitions in many hyperporphyrins and hypercorroles; thus, the complexes studied are aptly described as inverse hypercorroles.

Allyl-Allyl Coupling Promoted by Catalyst Systems with two Palladium Atoms - A Plethora of Potentially Pericyclic Processes

Recently, Huang and co-workers reported a catalytic reaction that utilizes H2 as the sole reductant for a C-C coupling of allyl groups with yields up to 96 %. Here we use computational quantum chemistry to identify several key features of this reaction that provide clarity on how it proceeds. We propose the involvement of a Pd-Pd bound dimer precatalyst, demonstrate the importance of ligand π-π interactions and counterions, and identify a new, energetically viable, mechanism involving two dimerized, outer-sphere reductive elimination transition structures that determine both the rate and selectivity. Although we rule out the previously proposed transmetalation step on energetic grounds, we show it to have an unusual aromatic transition structure in which two Pd atoms support rearranging electrons. The prevalence of potential metal-supported pericyclic reactions in this system suggests that one should consider such processes regularly, but the results of our calculations also indicate that one should do so with caution.

Enantioselective Synthesis of Phosphine Boranes via Crystallization-Induced Dynamic Resolution of Lithiated Intermediate by Understanding the Underlying Epimerization Process

Described herein is the successful crystallization-induced dynamic resolution (CIDR) of an α-lithiated phosphine borane utilizing the easily accessible and inexpensive ligand (R,R)-TMCDA. Starting from the essential P-prochiral building block dimethyl phenyl phosphine borane we were able to obtain phosphine boranes in yields up to 80 % and e.r. up to 98 : 2 by crystallization of the lithiated intermediate prior to the trapping reaction. NMR-based deuterium labeling experiments indicate that the epimerization in solution is based on the intermolecular proton transfer between nonlithiated phosphine borane and the corresponding lithiated intermediate, rendering the presence of the remaining starting compound in an optimized solvent mixture the main factor for successful enantioselective synthesis. Quantum chemical calculations using different model systems based on solid state structures confirm these experimental results. By gaining insights into the epimerization mechanism, essential principles for CIDR of lithiated phosphine boranes are elucidated that may be expanded to other important P-stereogenic compounds and simple chiral amines.

Covalent functionalization of germanene employing computational simulations

Computational simulations through density functional theory in conjunction with M06-L and HSE functional have been carried out to investigate the chemical reactivity of the germanene monolayer. It is exceptionally reactive, with an average reaction energy of -60.4 kcal mol-1 for the nineteen functional groups considered: H, F, Cl, Br, O, S, Se, Ge, OH, SH, CH3, CF3, NH, NH2, C6H5, C6H4, CCl2, CBr2, and the azomethine ylide. The results indicate that oxygen is the most reactive reagent (-96.9 kcal mol-1), followed by fluorine (-83.1 kcal mol-1). Germanene presents a rich organic chemistry, and functionalization with azomethine ylides, benzynes, and carbenes can be easily accomplished as indicated by the reaction energies computed, which lie between -45 and -65 kcal mol-1. Furthermore, germanene is significantly more reactive than graphene and hexagonal boron nitride monolayers since the reaction energy for germanene is more than 40 kcal mol-1 lower. Although, in general, germanene is slightly more reactive than black and blue phosphorene and less prone to oxidation, but its oxidation when exposed to air occurs spontaneously. The addition of functional groups works cooperatively. The reaction energies become lower as the number of functional groups increases, thus favouring the agglomeration of functional groups attached unless the steric effect alters this pattern. Finally, we analyzed the electronic properties of functionalized germanene. It is possible to fine-tune the band gap of germanene from 0.1 to 2 eV using different functional groups and coverages. For O-50% and S-50% functionalized germanene, we found that carrier recombination is the most difficult due to the considerable differences between the effective masses of holes and electrons, which is promising for optical applications.

QM/MM-Based Energy Decomposition Analysis Method for Large Systems

In this work, a QM/MM-based EDA method, called GKS-EDA(QM/MM), is proposed. As an extension of GKS-EDA, this scheme divides the total interaction energy into electrostatic, exchange-repulsion, polarization, and correlation/dispersion terms. GKS-EDA(QM/MM) can be applied to describe the interactions of large-scale systems combined with various QM/MM platforms. By using the examples of a hydrated hydronium ion complex in water solution, the barnase-barstar complex, and MMP-13-pyrimidinetrione in a metalloprotein, the capability of GKS-EDA(QM/MM) for various interactions in large systems is validated.

Theoretical Study of Cyanidin-Resveratrol Copigmentation by the Functional Density Theory

Anthocyanins are colored water-soluble plant pigments. Upon consumption, anthocyanins are quickly absorbed and can penetrate the blood-brain barrier (BBB). Research based on population studies suggests that including anthocyanin-rich sources in the diet lowers the risk of neurodegenerative diseases. The copigmentation caused by copigments is considered an effective way to stabilize anthocyanins against adverse environmental conditions. This is attributed to the covalent and noncovalent interactions between colored forms of anthocyanins (flavylium ions and quinoidal bases) and colorless or pale-yellow organic molecules (copigments). The present work carried out a theoretical study of the copigmentation process between cyanidin and resveratrol (CINRES). We used three levels of density functional theory: M06-2x/6-31g+(d,p) (d3bj); ωB97X-D/6-31+(d,p); APFD/6-31+(d,p), implemented in the Gaussian16W package. In a vacuum, the CINRES was found at a copigmentation distance of 3.54 Å between cyanidin and resveratrol. In water, a binding free energy ∆G was calculated, rendering -3.31, -1.68, and -6.91 kcal/mol, at M06-2x/6-31g+(d,p) (d3bj), ωB97X-D/6-31+(d,p), and APFD/6-31+(d,p) levels of theory, respectively. A time-dependent density functional theory (TD-DFT) was used to calculate the UV spectra of the complexes and then compared to its parent molecules, resulting in a lower energy gap at forming complexes. Excited states' properties were analyzed with the ωB97X-D functional. Finally, Shannon aromaticity indices were calculated and isosurfaces of non-covalent interactions were evaluated.

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers

Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.

Variation of Activation Volume as an Indicator of the Difference in Clusterization Phenomenon Induced by H-Bonding and F-Π Stacking Interactions in Enantiomers and a Racemate of Flurbiprofen

In this paper, X-ray diffraction (XRD), differential scanning calorimetry (DSC), broadband dielectric (BDS), and Fourier transform infrared (FTIR) spectroscopy supported by molecular dynamics (MD) simulations and quantum chemical computations were applied to investigate the structural and thermal properties, molecular dynamics, and H-bonding pattern of R-, S-, and RS-flurbiprofen (FLP). Experimental data indicated various spatial molecular arrangements in crystalline forms of examined systems, which seemed to disappear in the liquid state. Surprisingly, deeper analysis of high-pressure dielectric data revealed unexpected variation in the activation volume of pure enantiomers and a racemate. MD simulations showed that it is an effect of the clusterization phenomenon and a higher population of small associates in the former samples. Moreover, theoretical consideration exposed the particular role of unspecific F-Π interactions as a driving force underlying local molecular arrangements of molecules in the liquid and the crystal lattice of R-, S-, and RS-FLP.

Characterization of the Coordination and Solvation Dynamics of Solvated Systems─Implications for the Analysis of Molecular Interactions in Solutions and Pure H2O

The characterization of solvation shells of atoms, ions, and molecules in solution is essential to relate solvation properties to chemical phenomena such as complex formation and reactivity. Different definitions of the first-shell coordination sphere from simulation data can lead to potentially conflicting data on the structural properties and associated ligand exchange dynamics. The definition of a solvation shell is typically based on a given threshold distance determined from the respective solute-solvent pair distribution function g(r) (i.e., GC). Alternatively, a nearest neighbor (NN) assignment based on geometric properties of the coordination complex without the need for a predetermined cutoff criterion, such as the relative angular distance (RAD) or the modified Voronoi (MV) tessellation, can be applied. In this study, the effect of different NN algorithms on the coordination number and ligand exchange dynamics evaluated for a series of monatomic ions in aqueous solution, carbon dioxide in aqueous and dichloromethane solutions, and pure liquid water has been investigated. In the case of the monatomic ions, the RAD approach is superior in achieving a well separated definition of the first solvation layer. In contrast, the MV algorithm provides a better separation of the NNs from a molecular point of view, leading to better results in the case of solvated CO2. When analyzing the coordination environment in pure water, the cutoff-based GC framework was found to be the most reliable approach. By comparison of the number of ligand exchange reactions and the associated mean ligand residence times (MRTs) with the properties of the coordination number autocorrelation functions, it is shown that although the average coordination numbers are sensitive to the different definitions of the first solvation shell, highly consistent estimates for the associated MRT of the solvated system are obtained in the majority of cases.

Quantum Size-Driven Spectral Variations in Pillar[n]arene Systems: A Density Functional Theory and Wave Function Assessment

This study explores the quantum size effects on the optical properties of pillar[n]arene (n = 5, 6, 7, 8) utilizing density functional theory (DFT) and wave function analysis. The mechanisms of electron transitions in one-photon absorption (OPA) and two-photon absorption (TPA) spectra are investigated, alongside the calculation of electron circular dichroism (ECD) for these systems. Transition Density Matrix (TDM) and electron-hole pair density maps are employed to study the electron excitation characteristics, unveiling a notable size dependency. Analysis of the transition electric dipole moment (TEDM) and the transition magnetic dipole moment (TMDM) reveals the electromagnetic interaction mechanism within pillar[n]arene. Raman spectra computations further elucidate vibrational modes, while interactions with external environments are studied using electrostatic potential (ESP) analysis, and electron delocalization is assessed under an external magnetic field, providing insights into the magnetically induced current phenomena within these supramolecular structures. The thermal stability of pillar[n]arene was investigated by ab initio molecular dynamics (AIMD).

Phosphorescent Platinum(II) Complexes with a Spiro-fused Xanthene Unit: Synthesis and Photophysical Properties

Novel platinum(II) complexes, derived from the spiro[fluorene-9,9'-xanthene] (SFX) motif, were synthesized and combined with different auxiliary ligands such as acetylacetonate (acac), bis(2,4,6-trimethylphenyl)propane-1,3-dionate (mesacac) and dihydrobis(3,5-dimethylpyrazole-1-yl) borate. The final products were obtained in yields of up to 36 % and characterized by NMR, X-ray and combustion analysis. These complexes have structured green-blue emission spectra with Commission Internationale de l'Éclairage (CIExy) coordinates of (0.21;0.46). Excellent photoluminescence quantum yields (PLQYs) ranging from 87 %-91 % were found. The emission lifetimes vary from 33 μs to 43 μs. Calculations on the B3LYP/6-311++G** level of theory reveal, that the nature of the emissive state is dependent on the positional regioisomerism of the SFX motif. The 2-SFX complexes demonstrate ligand-centered (3LC) emission, while the 2'-SFX regioisomer with the mesacac ligand shows a strong 3MLCT character.

Puzzling Structure of the Key Intermediates in Gold(I)-catalyzed Cyclization Reactions of Enynes and Allenenes

We identify the dominant structures of the intermediates of gold(I)-catalyzed cyclizations of 1,5-enynes and 1,5-allenenes through computational analysis as gold(I) cyclopropylcarbenes, endocyclic vinylgold complexes and previously unreported non-classical carbocationic minima. In contrast to 1,6-enynes, the exocyclic carbocations are found to be less stable. Cyclopropylcarbene structures are consistently favoured as the most stable intermediates for all studied substitution patterns. We validate the computational methods used by using DLPNO-CCSD(T) energies as a benchmark, indicating that the B3LYP-D3 and M06-D3 functionals are most accurate for energy determination, while NPA charges are mostly insensitive to functional. The evolution of a 1,6-enyne in a single-cleavage or double-cleavage rearrangement is attributed to the barrierless evolution of a common cyclopropyl-gold(I) carbocation non-stationary geometry. Our findings provide insights into reaction pathways and substrate dependence of the cycloisomerization processes.

A Novel Abiotic Pathway for Phosphine Synthesis over Acidic Dust in Venus' Atmosphere

Recent ground-based observations of Venus have detected a single spectral feature consistent with phosphine (PH3) in the middle atmosphere, a gas which has been suggested as a biosignature on rocky planets. The presence of PH3 in the oxidized atmosphere of Venus has not yet been explained by any abiotic process. However, state-of-the-art experimental and theoretical research published in previous works demonstrated a photochemical origin of another potential biosignature-the hydride methane-from carbon dioxide over acidic mineral surfaces on Mars. The production of methane includes formation of the HC · O radical. Our density functional theory (DFT) calculations predict an energetically plausible reaction network leading to PH3, involving either HC · O or H· radicals. We suggest that, similarly to the photochemical formation of methane over acidic minerals already discussed for Mars, the origin of PH3 in Venus' atmosphere could be explained by radical chemistry starting with the reaction of ·PO with HC·O, the latter being produced by reduction of CO2 over acidic dust in upper atmospheric layers of Venus by ultraviolet radiation. HPO, H2P·O, and H3P·OH have been identified as key intermediate species in our model pathway for phosphine synthesis.

In silico investigation of falcipain-2 inhibition by hybrid benzimidazole-thiosemicarbazone antiplasmodial agents: A molecular docking, molecular dynamics simulation, and kinetics study

The emergence of artemisinin-resistant variants of Plasmodium falciparum necessitates the urgent search for novel antimalarial drugs. In this regard, an in silico study to screen antimalarial drug candidates from a series of benzimidazole-thiosemicarbazone hybrid molecules with interesting antiplasmodial properties and explore their falcipain-2 (FP2) inhibitory potentials has been undertaken herein. FP2 is a key cysteine protease that degrades hemoglobin in Plasmodium falciparum and is an important biomolecular target in the development of antimalarial drugs. Pharmacokinetic properties, ADMET profiles, MM/GBSA-based binding free energies, reaction mechanisms, and associated barrier heights have been investigated. DFT, molecular dynamics simulation, molecular docking, and ONIOM methods were used. From the results obtained, four 4N-substituted derivatives of the hybrid molecule (E)-2-(1-(5-chloro-1H-benzo[d]imidazol-2-yl)ethylidene)hydrazine-1-carbothioamide (1A) denoted 1B, 1C, 1D, and 1E are drug-like and promising inhibitors of FP2, exhibiting remarkably small inhibitory constants (5.94 × 10-14 - 2.59 × 10-04 n M) and favorable binding free energies (-30.32 to -17.17 kcal/mol). Moreover, the ONIOM results have revealed that 1B and possibly 1C and 1D may act as covalent inhibitors of FP2. The rate-determining step of the thermodynamically favorable covalent binding mechanism occurs across a surmountable barrier height of 24.18 kcal/mol in water and 28.42 kcal/mol in diethyl ether. Our findings are useful for further experimental investigations on the antimalarial activities of the hybrid molecules studied.

Insight into the interaction of host-guest structures for pyrrole-based metal compounds and C70

This study focuses on the recognition and isolation of fullerenes, which are crucial for further exploration of their physical and chemical properties. Our goal is to investigate the potential recognition of the D5h-C70 fullerene using crown-shaped metal compositions through density functional theory calculations. We assess the effectiveness of fullerene C70 recognition by studying the binding energy. Additionally, various analyses were conducted, including natural bond order charge analysis and reduced density gradient analysis, to understand the interaction mechanism between the host and guest molecules. These investigations provide valuable insights into the nature of the interaction and the stability of the host-guest system. To facilitate the release of the fullerene guest molecule, the vis-NIR spectra were simulated for the host-guest structures. This analysis offers guidance on the specific wavelengths that can be utilized to release the fullerene guest from the host-guest structures. Overall, this work proposes a new strategy for the effective recognition of various fullerene molecules and their subsequent release from host-guest systems. These findings could potentially be applied in assemblies involving fullerenes, advancing their practical applications.

High-Performance Multi-GPU Analytic RI-MP2 Energy Gradients

This article presents a novel algorithm for the calculation of analytic energy gradients from second-order Møller-Plesset perturbation theory within the Resolution-of-the-Identity approximation (RI-MP2), which is designed to achieve high performance on clusters with multiple graphical processing units (GPUs). The algorithm uses GPUs for all major steps of the calculation, including integral generation, formation of all required intermediate tensors, solution of the Z-vector equation and gradient accumulation. The implementation in the EXtreme Scale Electronic Structure System (EXESS) software package includes a tailored, highly efficient, multistream scheduling system to hide CPU-GPU data transfer latencies and allows nodes with 8 A100 GPUs to operate at over 80% of theoretical peak floating-point performance. Comparative performance analysis shows a significant reduction in computational time relative to traditional multicore CPU-based methods, with our approach achieving up to a 95-fold speedup over the single-node performance of established software such as Q-Chem and ORCA. Additionally, we demonstrate that pairing our implementation with the molecular fragmentation framework in EXESS can drastically lower the computational scaling of RI-MP2 gradient calculations from quintic to subquadratic, enabling further substantial savings in runtime while retaining high numerical accuracy in the resulting gradients.

Running Wild through Dirhodium Tetracarboxylate-Catalyzed Combined CH(C)-Functionalization/Cope Rearrangement Landscapes: Does Post-Transition-State Dynamic Mismatching Influence Product Distributions?

A special type of C-H functionalization can be achieved through C-H insertion combined with Cope rearrangement (CHCR) in the presence of dirhodium catalysts. This type of reaction was studied using density functional theory and ab initio molecular dynamics simulations, the results of which pointed to the dynamic origins of low yields observed in some experiments. These studies not only reveal intimate details of the complex reaction network underpinning CHCR reactions but also further cement the generality of the importance of nonstatistical dynamic effects in controlling Rh2L4-promoted reactions.

Unraveling Atomic Contributions to the London Dispersion Energy: Insights into Molecular Recognition and Reactivity

We present a general framework that enables quantification with atomic resolution of the overall London dispersion energy, which can be readily integrated with currently available energy decomposition schemes. This approach can be used to determine the contribution of individual atoms and functional groups to molecular recognition, conformational preferences, molecular stability, and reactivity. Its efficacy across diverse realms of molecular chemistry and biology is demonstrated with application to molecular balances in solution, asymmetric organocatalytic transformations, and a subcomplex of the F1FO ATP synthase.

Evidence of Excited-State Vibrational Mode Governing the Photorelaxation of a Charge-Transfer Complex

Modern, nonlinear, time-resolved spectroscopic techniques have opened new doors for investigating the intriguing but complex world of photoinduced ultrafast out-of-equilibrium phenomena and charge dynamics. The interaction between light and matter introduces an additional dimension, where the complex interplay between electronic and vibrational dynamics needs the most advanced theoretical-computational protocols to be fully understood on the molecular scale. In this study, we showcase the capabilities of ab initio molecular dynamics simulation integrated with a multiresolution wavelet protocol to carefully investigate the excited-state relaxation dynamics in a noncovalent complex involving tetramethylbenzene (TMB) and tetracyanoquinodimethane (TCNQ) undergoing charge transfer (CT) upon photoexcitation. Our protocol provides an accurate description that facilitates a direct comparison between transient vibrational analysis and time-resolved spectroscopic signals. This molecular level perspective enhances our understanding of photorelaxation processes confined in the adiabatic regime and offers an improved interpretation of vibrational spectra. Furthermore, it enables the quantification of anharmonic vibrational couplings between high- and low-frequency modes, specifically the TCNQ "rocking" and "bending" modes. Additionally, it identifies the primary vibrational mode that governs the adiabaticity between the ground state and the CT state. This comprehensive understanding of photorelaxation processes holds significant importance in the rational design and precise control of more efficient photovoltaic and sensor devices.

Mechanism of Nitrous Oxide Activation in C(sp2)-O Bond Formation Reactions Catalyzed by Nickel Complexes

Recent groundbreaking experimental reports demonstrated that Ni complexes bearing a bidentate- or tridentate-bipyridine-based ligand can be used to activate N2O for use as an O-transfer agent in C(sp2)-O bond formation reactions under mild experimental conditions. In this work, quantum chemical calculations are used to shed light on the mechanism through which such metal complexes catalytically activate nitrous oxide, providing new fundamental insights into the development of novel catalysts for N2O revalorization. As a case study, we consider the recent work by Cornella and co-workers (Nature, 2022, 604, 677) concerning the synthesis of phenols from aryl halides at room temperature, which requires the use of an external reducing agent. Our results suggest that the metal center remains in its Ni(II) oxidation state throughout the whole catalytic cycle, despite the presence of various redox steps in the mechanism and the Ni ability to maneuver between a number of oxidation states. This counterintuitive behavior is made possible by the ligand redox activity in the catalytic process, which involves accepting electrons from the reducing agent. Several possible pathways are systematically investigated, each associated with distinct activation modes, kinetics, and reaction outcomes. The governing factors in dictating the preferred path lie in the electronic nature of the ligand (strong vs weak field) and its geometric structure (specifically, the number of coordinating arms). These characteristics play a pivotal role in determining whether the process follows a catalytic or stoichiometric route and can be in principle modulated for the design of new metal complexes with tailored redox properties and reactivity.

Computational insight into the crystal structures of cubane and azacubanes

CONTEXT

Using quantum chemistry and atom-atom potential methods, the molecular and crystal structures of cubane 1 and all types of unsubstituted azacubanes 2-22 were calculated. Alternative possible polymorphs of cubane 1 have been proposed. The thermochemical properties of azacubanes in the gas and solid phases were assessed. Thermodynamic aspects of stability are considered, and a significant decrease in stability is revealed upon transition from cubane 1 to octaazacubane 22. It has been shown that the density and energetic properties of azacubanes depend nonlinearly on the number of nitrogen atoms in the structure and the density of octaazacubane 22 at room temperature is 1.546 g cm-3, which is significantly lower than the previously given estimate.

METHODS

In this work, DFT calculations were conducted through the software Gaussian 09 using B3LYP functional with basis set aug-cc-PVDZ and the Grimme dispersion correction D2. For crystal structure optimization, the atom-atom potential methods with PMC (packing of molecules in crystal) program were used. Charges for molecular electrostatic potential were fitted by FitMEP, and enthalpies of formation in gas phase were assessed by G3B3.

The energy landscape of the ribosome

The ribosome is a prototypical assembly that can be used to establish general principles and techniques for the study of biological molecular machines. Motivated by the fact that the dynamics of every biomolecule is governed by an underlying energy landscape, there has been great interest to understand and quantify ribosome energetics. In the present review, we will focus on theoretical and computational strategies for probing the interactions that shape the energy landscape of the ribosome, with an emphasis on more recent studies of the elongation cycle. These efforts include the application of quantum mechanical methods for describing chemical kinetics, as well as classical descriptions to characterize slower (microsecond to millisecond) large-scale (10-100 Å) rearrangements, where motion is described in terms of diffusion across an energy landscape. Together, these studies provide broad insights into the factors that control a diverse range of dynamical processes in this assembly.

The Physical Mechanism of Linear and Nonlinear Optical Properties of Nanographene-Induced Chiral Inversion

Based on density functional theory (DFT) and wave function analysis, the ultraviolet and visible spectrophotometry (UV-Vis) spectra and Raman spectra of 1-meso and 1-rac obtained by the chiral separation of chiral nanographenes are theoretically investigated. The electron excitation properties of 1-meso and 1-rac are studied by means of transition density matrix (TDM) and charge density difference (CDD) diagrams. The intermolecular interaction is discussed based on an independent gradient model based on Hirshfeld partition (IGMH). The interaction of 1-meso and 1-rac with the external environment is studied using the electrostatic potential (ESP), and the electron delocalization degree of 1-meso and 1-rac is studied based on the magnetically induced current under the external magnetic field. Through the chiral separation of 1-rac, two enantiomers, 1-(P, P) and 1-(M, M), were obtained. The electrical-magnetic interaction of the molecule is revealed by analyzing the electron circular dichroism (ECD) spectra of 1-meso, 1-(P, P) and 1-(M, M), the transition electric dipole moment (TEDM) and the transition magnetic dipole moment (TMDM). It is found that 1-(P, P) and 1-(M, M) have opposite chiral properties due to the inversion of the structure.

Probing Non-Covalent Interactions through Molecular Balances: A REG-IQA Study

The interaction energies of two series of molecular balances (1-X with X = H, Me, OMe, NMe2 and 2-Y with Y = H, CN, NO2, OMe, NMe2) designed to probe carbonyl…carbonyl interactions were analysed at the B3LYP/6-311++G(d,p)-D3 level of theory using the energy partitioning method of Interacting Quantum Atoms/Fragments (IQA/IQF). The partitioned energies are analysed by the Relative Energy Gradient (REG) method, which calculates the correlation between these energies and the total energy of a system, thereby explaining the role atoms have in the energetic behaviour of the total system. The traditional "back-of-the-envelope" open and closed conformations of molecular balances do not correspond to those of the lowest energy. Hence, more care needs to be taken when considering which geometries to use for comparison with the experiment. The REG-IQA method shows that the 1-H and 1-OMe balances behave differently to the 1-Me and 1-NMe2 balances because the latter show more prominent electrostatics between carbonyl groups and undergoes a larger dihedral rotation due to the bulkiness of the functional groups. For the 2-Y balance, REG-IQA shows the same behaviour across the series as the 1-H and 1-OMe balances. From an atomistic point of view, the formation of the closed conformer is favoured by polarisation and charge-transfer effects on the amide bond across all balances and is counterbalanced by a de-pyramidalisation of the amide nitrogen. Moreover, focusing on the oxygen of the amide carbonyl and the α-carbon of the remaining carbonyl group, electrostatics have a major role in the formation of the closed conformer, which goes against the well-known n-π* interaction orbital overlap concept. However, REG-IQF shows that exchange-correlation energies overtake electrostatics for all the 2-Y balances when working with fragments around the carbonyl groups, while they act on par with electrostatics for the 1-OMe and 1-NMe2. REG-IQF also shows that exchange-correlation energies in the 2-Y balance are correlated to the inductive electron-donating and -withdrawing trends on aromatic groups. We demonstrate that methods such as REG-IQA/IQF can help with the fine-tuning of molecular balances prior to the experiment and that the energies that govern the probed interactions are highly dependent on the atoms and functional groups involved.

Interplay between protonation and Jahn-Teller effects in a manganese vanadium cubane water oxidation catalyst

Understanding the protonation behavior of metal-oxo water oxidation catalysts is essential to improve catalyst efficiency and long-term performance, as well as to tune their properties for specific applications. In this work, we explore the basicity and protonation effects of the highly active water oxidation catalyst [(Mn4O4) (V4O13) (OAc)3]3- using density functional theory. We computed the relative free energies of protonation in a systematic fashion for all symmetry-inequivalent O atoms, where the presence of multiple oxidation states from Mn4IV to Mn4III and a rich Jahn-Teller isomerism adds a significant amount of complexity. For high oxidation states, the compound behaves like some other polyoxometalates, showing protonation preferably at the terminal and μ2-bridging O atoms of the vanadate cap. However, upon reduction, eventually, the protonation preference switches to the cubane O atoms, mostly driven by a strong increase in basicity for O atoms located along the Jahn-Teller axes. Our work further evidences that protonation can potentially lead to several chemical transformations, like disproportionation and charge transfer to vanadium, dissociation of ligands, or the opening of the cubane structure. Our simulated UV/Vis absorption spectra additionally provide valuable insights about how the protonation of the catalyst could be tracked experimentally. Overall, our analysis highlights the complexity involved in the protonation of heterometallic polyoxometalate clusters.

A physical organic strategy to predict and interpret stabilities of chemical bonds in energetic compounds for the discovery of thermal-resistant properties

CONTEXT

The in-depth understanding about the stability of chemical bonds in energetic compounds plays a central role for molecular design and safety-related evaluations. Most energetic compounds contain nitro as explosophores, and nitro cleavage is fundamental for thermal and mechanical stability. However, the quantum chemistry approach to accurately predict energy and temperature properties related to bond stability is challenging, due to the tradeoff between computational costs and deviations. Herein, the bond orders are proposed as accurate and computational-cost efficient descriptors for predicting the chemical bond stability and thermal-resistant properties. The intrinsic bond strength index (IBSI) demonstrates the best prediction for experimental homolytic bond dissociation energies (R2 > 0.996), which is on par with the results from high-precision quantum chemistry methods. The effects from bond connectivity and steric hindrance hierarchy were analyzed to reveal underlying mechanisms. Additionally, the IBSI descriptors are successfully applied to predict the thermal decomposition temperatures of 24 heat-resistant energetic compounds (R2 = 0.995), thus validating the effectiveness for the prediction and interpretation of chemical bond stability in energetic compounds via a physical organic approach.

METHODS

All DFT calculations were performed with Gaussian 09 software. To investigate the dependence of the method on functionals and basis sets, 9 DFT methods were considered (B3LYP/6-31G(d,p), B3LYP/6-311G(d,p), B3LYP/def2-TZVP, M062X/6-31G(d,p), M062X/6-311G(d,p), M062X/def2-TZVP, ωB97XD/6-31G(d,p), ωB97XD/6-311G(d,p), and ωB97XD/def2-TZVP). The bond order descriptors LBO and IBSI are obtained through the bond order analysis module in the Multiwfn software.

Decoding Van der Waals Impact on Chirality Transfer in Perovskite Structures: Density Functional Theory Insights

Chiral organic-inorganic perovskites exhibit unique physicochemical properties driven by the symmetry of monovalent organic cations. However, an atomistic understanding of how chiral cations transfer their chirality to the inorganic framework and the role played by van der Waals (vdW) interactions in this process is still incomplete. In this work, we report a theoretical investigation, based on density functional theory calculations within the Perdew-Burke-Ernzerhof (PBE) formulation for the exchange-correlation functional, into the role of the vdW interactions in the chirality transfer process. For that, we selected several vdW corrections, namely, Grimme (D2, D3, D3(BJ)), Tkatchenko-Scheffler (TS, TS+SCS, TS+HSI), density-dependent energy correction (dDsC), and many-body scattering (MBD) energy method correction. For the chiral perovskite systems, we selected a set of chiral organic-inorganic perovskites with several dimensions, namely, from zero-dimensional to three-dimensional, each having enantiomers with R and S configurations. Based on a statistical treatment of the relative errors of all lattice parameters with respect to experimental data, we found that D3, D3(BJ), TS, TS+SCS, TS+HSI, and MBD vdW are the most accurate corrections to describe the equilibrium structural properties of chiral perovskites using the PBE method. We identify chirality-induced sequential asymmetries of distorted octahedrons and propose angular descriptors to quantify them, where the orientations of these distortions depend on the R or S nature of the chiral cations. Furthermore, we demonstrate the importance of accurate vdW interactions in precisely describing these asymmetric distortions. By means of binding energies and charge-transfer analysis, we show that the impact of vdW corrections on the charge distribution leads to a subtle strengthening of hydrogen bonds between chiral cations and inorganic octahedra, resulting in an increase in the binding energy. Finally, we identified that the Rashba-Dresselhaus effect in two-dimensionality is refined by vdW interactions.

Comprehensive Study on the Adsorption and Degradation of Dichlorodiphenyltrichloroethane on Bifunctional Adsorption-Photocatalysis Material TiO2/MCM-41 Using Quantum Chemical Methods

The adsorption and degradation capacities of dichlorodiphenyltrichloroethane (DDT) on a photocatalyst composed of TiO2 supported on the mesoporous material MCM-41 (TiO2/MCM-41) were investigated using density functional theory and real-time density functional theory methods. The van der Waals interactions within the PBE functional were adjusted by using the Grimme approach. The adsorption of DDT was evaluated through analyses involving adsorption energy, Hirshfeld atomic charges, Wiberg bond orders, molecular electrostatic potential, noncovalent interaction analysis, and bond path analysis. The findings reveal that DDT undergoes physical adsorption on pristine MCM-41 or MCM-41 modified with Al or Fe due to the very small bond order (only about 0.15-0.18) as well as the change in total charge of DDT after adsorption is close to 0. However, it chemically adsorbs onto the TiO2/MCM-41 composite through the formation of Ti···Cl coordination bonds because the maximum bond order is very large, about 1.0 (it can be considered as a single bond). The adsorption process is significantly influenced by van der Waals interactions (accounting for approximately 30-40% of the interaction energy), hydrogen bonding, and halogen bonding. MCM-41 is demonstrated to concurrently function as a support for the TiO2 photocatalyst, creating a synergistic effect that enhances the photocatalytic activity of TiO2. Based on the computational results, a novel photocatalytic mechanism for the degradation of DDT on the TiO2/MCM-41 catalyst system was proposed.

Uniaxial Strain-Induced Stacking Order Change in Trilayer Graphene

The layer stacking order in two-dimensional heterostructures, like graphene, affects their physical properties and potential applications. Trilayer graphene, specifically ABC-trilayer graphene, has captured significant interest due to its potential for correlated electronic states. However, achieving a stable ABC arrangement is challenging due to its lower thermodynamic stability compared to the more stable ABA stacking. Despite recent advancements in obtaining ABC graphene through external perturbations, such as strain, the stacking transition mechanism remains insufficiently explored. In this study, we unveil a universal mechanism to achieve ABC stacking, applicable for understanding ABA to ABC stacking changes induced by any mechanical perturbations. Our approach is based on a novel strain engineering technique that induces interlayer slippage and results in the formation of stable ABC domains. We investigate the underlying interfacial mechanisms of this stacking change through computational simulations and experiments. Our findings demonstrate a highly anisotropic and significant transformation of ABA stacking to large and stable ABC domains facilitated by interlayer slippage. Through atomistic simulations and local energy analysis, we systematically demonstrate the mechanism for this stacking transition, that is dependent on specific loading orientation. Understanding such a mechanism allows this material system to be engineered by design compatible with industrial techniques on a device-by-device level. We conduct Raman studies to validate and characterize the formed ABC stacking, highlighting its distinct features compared to the ABA region. Our results contribute to a clearer understanding of the stacking change mechanism and provide a robust and controllable method for achieving stable ABC domains, facilitating their use in developing advanced optoelectronic devices.

Ab Initio Potential Energy Surface for NaCl-H2 with Correct Long-Range Behavior

We report a full dimensional ab initio potential energy surface for NaCl-H2 based on precise fitting of a large data set of CCSD(T)/aug-cc-pVTZ energies. A major goal of this fit is to describe the very long-range interaction accurately. This is done in this instance via the dipole-quadrupole interaction. The NaCl dipole and the H2 quadrupole are available through previous works over a large range of internuclear distances. We use these to obtain exact effect charges on each atom. Diffusion Monte Carlo calculations are done for the ground vibrational state using the new potential.

Long-Chain Bio-Based Nylon 514 Salt: Crystal Structure, Phase Transformation, and Polymerization

Nylon 514 is one of the new long-chain bio-based nylon materials; its raw material, 1,5-pentanediamine (PDA), is prepared by biological techniques, using biomass as the raw material. The high-performance monomer of nylon 514, 1,5-pentanediamine-tetradecanedioate (PDA-TDA) salt, was obtained through efficient crystallization methods. Here, two crystal forms of PDA-TDA, anhydrous and dihydrate, were identified and studied in this paper. From the characterization data, their crystal structures and thermal behaviors were investigated. Lattice energy was calculated to gain further insight into the relationship between thermal stability and crystal structures. The contribution of hydrogen bonds and other intermolecular interactions to the crystal structure stability have been quantified according to detailed Hirshfeld and IRI analyses. Additionally, the transformation mechanism of the anhydrate and dihydrate was established through a series of well-designed stability experiments, in which the temperature and water activity play a significant role in the structural stability of crystalline forms. Eventually, we obtained nylon 514 products with good thermal stability and low absorption using stable dihydrate powders as monomers. The properties of nylon 514 products prepared by different polymerization methods were also compared.

Biaxial compressive strain enhances calcium binding and mobility on two-dimensional Sc2C: a density functional theory investigation

In this work, we investigated calcium binding and diffusion on pristine and biaxially strained 2D Sc2C via density functional theory calculations, for potential applications in calcium-ion batteries (CIBs). We found that 2D Sc2C is metallic under PBE, HSE06, and DFT+U approximation conditions, and thus can be potentially used as an electrode material for CIBs. Results showed that pristine 2D Sc2C adsorbs calcium modestly, with relatively low binding energy on the most stable site (0.38 eV). Interestingly, this value shoots up to -1.94 eV and -3.23 eV at 5% and 10% biaxial compressive strains, respectively. Furthermore, calcium's diffusion energy barrier, which is already low (80 meV) on pristine 2D Sc2C, goes down further (to 35 meV) upon application of median biaxial compressive strain (5%). As a result of the enhanced binding of calcium on strained 2D Sc2C, the maximum stable calcium concentration was also boosted. Consequently, the calculated theoretical specific energy capacity of 2D Sc2C with biaxial compressive strain is higher compared to that of the pristine case (878.29 mA h g-1vs. 1051.84 mA h g-1). The average open circuit voltages of the two cases are high and quite close at 9.3 V (pristine) and 9.0 V (with 5% biaxial compressive strain). Our results demonstrated that biaxial compressive strain could be tapped to improve the properties of 2D MXenes, such as Sc2C, thereby enhancing the battery performance indicators of these materials, such as theoretical specific energy capacity and open circuit voltage. Such findings are of great importance in the emerging new technology of CIBs.

Structural, Electronic and Vibrational Properties of B24N24 Nanocapsules: Novel Anodes for Magnesium Batteries

We report on DFT-TDDFT studies of the structural, electronic and vibrational properties of B24N24 nanocapsules and the effect of encapsulation of homonuclear diatomic halogens (Cl2, Br2 and I2) and chalcogens (S2 and Se2) on the interaction of the B24N24 nanocapsules with the divalent magnesium cation. In particular, to foretell whether these BN nanostructures could be proper negative electrodes for magnesium-ion batteries, the structural, vibrational and electronic properties, as well as the interaction energy and the cell voltage, which is important for applications, have been computed for each system, highlighting their differences and similarities. The encapsulation of halogen and chalcogen diatomic molecules increases the cell voltage, with an effect enhanced down groups 16 and 17 of the periodic table, leading to better performing anodes and fulfilling a remarkable cell voltage of 3.61 V for the iodine-encapsulated system.

Paving the road towards automated homogeneous catalyst design

In the past decade, computational tools have become integral to catalyst design. They continue to offer significant support to experimental organic synthesis and catalysis researchers aiming for optimal reaction outcomes. More recently, data-driven approaches utilizing machine learning have garnered considerable attention for their expansive capabilities. This Perspective provides an overview of diverse initiatives in the realm of computational catalyst design and introduces our automated tools tailored for high-throughput in silico exploration of the chemical space. While valuable insights are gained through methods for high-throughput in silico exploration and analysis of chemical space, their degree of automation and modularity are key. We argue that the integration of data-driven, automated and modular workflows is key to enhancing homogeneous catalyst design on an unprecedented scale, contributing to the advancement of catalysis research.

Harnessing the electronic structure of active metals to lower the overpotential of the electrocatalytic oxygen evolution reaction

Despite substantial advancements in the field of the electrocatalytic oxygen evolution reaction (OER), the efficiency of earth-abundant electrocatalysts remains far from ideal. The difficulty stems from the complex nature of the catalytic system, which limits our fundamental understanding of the process and thus the possibility of a rational improvement of performance. Herein, we shed light on the role played by the tunable 3d configuration of the metal centers in determining the OER catalytic activity by combining electrochemical and spectroscopic measurements with an experimentally validated computational protocol. One-dimensional coordination polymers based on Fe, Co and Ni held together by an oxonato linker were selected as a case study because of their well-defined electronic and geometric structure in the active site, which can be straightforwardly correlated with their catalytic activity. Novel heterobimetallic coordination polymers were also considered, in order to shed light on the cooperativity effects of different metals. Our results demonstrate the fundamental importance of electronic structure effects such as metal spin and oxidation state evolutions along the reaction profile to modulate ligand binding energies and increase catalyst efficiency. We demonstrated that these effects could in principle be exploited to reduce the overpotential of the electrocatalytic OER below its theoretical limit, and we provide basic principles for the development of coordination polymers with a tailored electronic structure and activity.

Classical Gold Carbonyl Complexes in Tetrahedral and Trigonal-Planar Settings

A unique four-coordinate, classical gold(I)-carbonyl complex with substantial backdonation from gold has been isolated by using a B-methylated and fluorinated tris(pyridyl)borate chelator. Its lighter silver(I) and copper(I) analogs enabled a study of trends in the coinage-metal family. The B-arylated ligand version also afforded a gold-carbon monoxide complex that displays a notably low C-O stretch value, but with trigonal planar geometry at the gold. A computational analysis shows that the AuI -CO bonds of these tris(pyridyl)borate ligand-supported molecules consist of electrostatic attraction, OC→Au σ-donation, and very significant Au→CO π-back-bonding components. The latter is responsible for the observed C-O stretching frequencies, which are lower than in free CO.