Quantum Chemical Calculations Research Articles

Indium(III)-Catalyzed Haloalkynylation Reaction of Alkynes.

Green chemistry strives for sustainability at the molecular level and is gaining increasing relevance in the development of chemical reactions. The haloalkynylation reaction is a highly atom-economical C-C coupling reaction that was previously only achieved using transition metal catalysts. It enables the introduction of an alkyne unit and a halogen atom into the target molecule. Herein, we present a haloalkynylation reaction catalyzed by indium(III) halides. The use of indium(III) bromide as a catalyst leads exclusively to the cis addition products with yields up to 86 %. In addition, iodoacetylenes can be applied for the first time for the haloalkynylation reaction of internal alkynes which is an important step forward in the development of industrially relevant and sustainable catalysts. In contrast to gold catalysis, which proceeds via a similar mechanism, the use of alkyl-substituted haloacetylenes as reagents is also possible. Based on 13C labeling experiments and quantum chemical calculations, we postulate two possible mechanisms for the indium(III)-catalyzed haloalkynylation reactions.

Cyсlopentadienone and Pyrone Derivatives as Precursors of Electron-Deficient Cycloheptatrienes: Quantum Chemical Investigation and Synthesis

AbstractUnstable tetra(methoxycarbonyl)cyclopentadienone was investigated in the synthesis of electron-deficient cycloheptatrienes via [4+2]-cycloaddition/cycloelimination reaction with cyclopropenes. The use of its stable dimer did not afford the product although similar reactions with alkynes have been reported. Quantum chemical calculation revealed that cyclopentadienone is not generated from the dimer and the reaction with alkynes proceeds via a more complicated cycloelimination/cycloaddition/cycloelimination cascade. However, the formation of cycloheptatrienes was found favorable over the formation of the dimer. Therefore, the trapping of tetra(methoxycarbonyl)cyclopentadienone upon formation was successful to give cycloheptatrienes with five ester groups. The use of methyl coumalate as a four-electron component was successful with cyclopropenes containing only one ester group to afford only two ester groups in the product.

Efficient simulation of potential energy operators on quantum hardware: a study on sodium iodide (NaI)

This study introduces a conceptually novel polynomial encoding algorithm for simulating potential energy operators encoded in diagonal unitary forms in a quantum computing machine. The current trend in quantum computational chemistry is effective experimentation to achieve high-precision quantum computational advantage. However, high computational gate complexity and fidelity loss are some of the impediments to the realization of this advantage in a real quantum hardware. In this study, we address the challenges of building a diagonal Hamiltonian operator having exponential functional form, and its implementation in the context of the time evolution problem (Hamiltonian simulation and encoding). Potential energy operators when represented in the first quantization form is an example of such types of operators. Through systematic decomposition and construction, we demonstrate the efficacy of the proposed polynomial encoding method in reducing gate complexity from O(2n)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {O}(2^n)$$\\end{document} to O∑i=1rnCr\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {O}\\left( \\sum _{i=1}^{r} {} ^nC_r \\right)$$\\end{document} (for some r≪n\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$r\\ll n$$\\end{document}). This offers a solution with lower complexity in comparison to the conventional Hadamard basis encoding approach. The effectiveness of the proposed algorithm was validated with its implementation in the IBM quantum simulator and IBM quantum hardware. This study demonstrates the proposed approach by taking the example of the potential energy operator of the sodium iodide molecule (NaI) in the first quantization form. The numerical results demonstrate the potential applicability of the proposed method in quantum chemistry problems, while the analytical bound for error analysis and computational gate complexity discussed, throw light on issues regarding its implementation.

An interactive approach to the notion of chemical substance and the case of water

AbstractFrom organic synthesis to quantum chemical calculation, chemists interact with chemical substances in a wide variety of ways. But what even is a chemical substance? My aim is to propose a notion of chemical substance that is consistent with the way in which chemical substances are individuated in chemistry, addressing gaps in previous conceptions of chemical substance. Water is employed as a case study to develop the account, not only because it is a familiar example of a chemical substance, but also because its structural peculiarities make it an ideal test case for drawing out potential issues and limitations. Examining four distinct views of chemical substance—the microstructural, thermodynamic, purification, and a functional/relational account—I conclude that each has considerable drawbacks when used as a standalone concept. However, these accounts are not rendered obsolete, but are combined into a semi-pluralist conceptual patchwork. My interactive account of chemical substance is consistent with existing substance descriptions and chemical practice.

A Tumor Environment-Activated Photosensitized Biomimetic Nanoplatform for Precise Photodynamic Immunotherapy of Colon Cancer.

Aggressive nature of colon cancer and current imprecise therapeutic scenarios simulate the development of precise and effective treatment strategies. To achieve this, a tumor environment-activated photosensitized biomimetic nanoplatform (PEG2000-SiNcTI-Ph/CpG-ZIF-8@CM) is fabricated by encapsulating metal-organic framework loaded with developed photosensitizer PEG2000-SiNcTI-Ph and immunoadjuvant CpG oligodeoxynucleotide within fusion cell membrane expressing programmed death protein 1 (PD-1) and cluster of differentiation 47 (CD47). By stumbling across, systematic evaluation, and deciphering with quantum chemical calculations, a unique attribute of tumor environment (low pH plus high concentrations of adenosine 5'-triphosphate (ATP))-activated photodynamic effect sensitized by long-wavelength photons is validated for PEG2000-SiNcTI-Ph/CpG-ZIF-8@CM, advancing the precision of cancer therapy. Moreover, PEG2000-SiNcTI-Ph/CpG-ZIF-8@CM evades immune surveillance to target CT26 colon tumors in mice mediated by CD47/signal regulatory proteins α (SIRPα) interaction and PD-1/programmed death ligand 1 (PD-L1) interaction, respectively. Tumor environment-activated photodynamic therapy realized by PEG2000-SiNcTI-Ph/CpG-ZIF-8@CM induces immunogenic cell death (ICD) to elicit anti-tumor immune response, which is empowered by enhanced dendritic cells (DC) uptake of CpG and PD-L1 blockade contributed by the nanoplatform. The photodynamic immunotherapy efficiently combats primary and distant CT26 tumors, and additionally generates immune memory to inhibit tumor recurrence and metastasis. The nanoplatform developed here provides insights for the development of precise cancer therapeutic strategies.

Electronic properties of polyethylene naphthalate as derived from photo-stimulated discharge, luminescence experiments and quantum chemical calculation

The electronic properties of thin films of poly(ethylene 2,6-naphthalate)—PEN, are investigated based on their photo-physical (optical absorption, photoluminescence) and electrical (space charge distribution, photo-stimulated discharge) behavior. Photo-stimulated currents are associated with optical absorption of the material leading to space charge dissipation as demonstrated by space charge distribution measurement. Based on this set of experimental results and quantum chemical calculation performed on PEN macromolecular system, we propose a new scheme for the electronic levels of PEN. This scheme allows understanding the mechanisms at play in photo-stimulated discharge. One of the main conclusions of our work is that photo-stimulated current measurements do not probe the energy level of traps. Detrapping of charges results from a two-step process where the photon energy is absorbed by chromophores that restitute a part of this energy to trapped charges through various mechanisms. Moreover, the new scheme allows discussing the components of the luminescence excited under different stresses, being electric field, electronic and UV irradiation, charge recombination and thermal activation.

THEORETICAL PREDICTION AND COMPREHENSIVE CHARACTERIZATION OF AN ALL-NITROGENATOMIC RING, CYCLO[18]NITROGEN (N18).

The cyclic molecule cyclo[18]carbon composed of 18 carbon atoms has been observed in condensed phase experiment in recent years and has attracted great attention. Through state-of-art quantum chemistry calculation, this study found that 18 nitrogen atoms can also form a macrocyclic system, cyclo[18]nitrogen (N18), though its lifetime is very short at room temperature and can only exist for a relatively long time at very low temperatures. We comprehensively theoretically studied properties of N18, including geometric configurations, thermal decomposition mechanism and rate, molecular dynamics behavior, energetic properties, vibrational and electronic spectra. We also discussed in depth the electronic structure of N18, including nature of the N-N bonds, lone-pairs, charge distribution characteristics, electronic delocalization, and aromaticity. This work is not only the first exploration of the macrocyclic N18 molecule, but also the first time to systematically examine a very long-chain substance fully composed of nitrogen atoms in isolated state.

Toward Real Chemical Accuracy on Current Quantum Hardware Through the Transcorrelated Method.

Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and connectivity) hamper the implementation of most quantum algorithms and call for more noise-resilient solutions. We propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach to target these major roadblocks directly. This method transfers, without any approximation, correlations from the wave function directly into the Hamiltonian, thus reducing the resources needed to achieve accurate results with noisy quantum devices. We show that the TC approach allows for shallower circuits and improves the convergence toward the complete basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, thus, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using two and four qubits, respectively. To demonstrate our approach's current and near-term potential, we perform hardware experiments, where our results confirm that the TC method paves the way toward accurate quantum chemistry calculations already on today's quantum hardware.

Adaptive Response to Solvation in Flexible Molecules: Oligo Hydrates of 4‐hydroxy‐2‐butanone

Structural changes induced by water play a pivotal role in chemistry and biology but remain challenging to predict, measure, and control at molecular level. Here we explore size‐governed gas‐phase water aggregation in the flexible molecule 4‐hydroxy‐2‐butanone, modeling the conformational adaptability of flexible substrates to host water scaffolds and the preference for sequential droplet growth. The experiment was conducted using broadband rotational spectroscopy, rationalized with quantum chemical calculations. Two different isomers were observed experimentally from the di‐ to the pentahydrates (4‐hydroxy‐2‐butanone‐(water)n = 2‐5), including the 18O isotopologues for the di‐ and trihydrates. Interestingly, to accommodate water molecules effectively, the heavy atom skeleton of 4‐hydroxy‐2‐butanone reshapes in every observed isomer and does not correspond to the stable conformer of the free monomer. All solvates initiate from the (proton donor) alcohol group but retain the carbonyl group as secondary binding point. The water scaffolds closely resemble those found in the pure water clusters, balancing between the capability of 4‐hydroxy‐2‐butanone for steering the orientation and position of the water molecules and the ability of water to modulate the monomer's conformation. The present work thus provides an accurate molecular answer on how torsionally flexible molecules dynamically adapt to water along progressing solvation.

Adaptive Response to Solvation in Flexible Molecules: Oligo Hydrates of 4-hydroxy-2-butanone.

Structural changes induced by water play a pivotal role in chemistry and biology but remain challenging to predict, measure, and control at molecular level. Here we explore size-governed gas-phase water aggregation in the flexible molecule4-hydroxy-2-butanone, modeling the conformational adaptability of flexible substrates to host water scaffolds and the preference for sequential droplet growth. The experiment was conducted using broadband rotational spectroscopy, rationalized with quantum chemical calculations. Two different isomers were observed experimentally from the di- to the pentahydrates (4-hydroxy-2-butanone-(water)n = 2-5), including the 18O isotopologues for the di- and trihydrates. Interestingly, to accommodate water molecules effectively, the heavy atom skeleton of 4-hydroxy-2-butanone reshapes in every observed isomer and does not correspond to the stable conformer of the free monomer. All solvates initiate from the (proton donor) alcohol group but retain the carbonyl group as secondary binding point. The water scaffolds closely resemble those found in the pure water clusters, balancingbetween the capability of 4-hydroxy-2-butanone for steering the orientation and position of the water molecules and the ability of water to modulate the monomer's conformation. The present work thus providesan accurate molecular answer on how torsionally flexible molecules dynamically adapt to water along progressing solvation.

Invisible photochromism and on-off switch based on thiophene-imidazole fused derivatives: Synthesis, properties and applications

Three novel diarylethene fluorescent molecules based on a thiophene-imidazole core have been successfully synthesized. Their photophysical properties, including photocyclization and ions response characteristics, as well as quantum chemical calculations were systematically investigated. All compounds exhibited significant invisible photochromism upon photocyclization, resulting in remarkable fluorescence “turn-off” properties. This phenomenon can be attributed to photoisomerization, specifically the transition from the open-ring isomer to the closed-ring isomer during this process, which was confirmed by 1H NMR data analysis. Moreover, AT-DPIo demonstrated high selectivity toward I− with pronounced changes in absorption from colorless to yellow and fluorescence from bright yellow to dark, enabling a limit of detection for I− at 9.75 × 10−7 mol L−1. These exceptional sensing properties make it suitable for rapid “naked-eye” detection of I− using test strips.

Understanding of mercury/oxygen reaction mechanism over BiOBr photocatalyst with oxygen vacancy

BiOBr with the advantage of easily generated surface oxygen vacancy, has been considered as a promising photocatalyst for Hg0 oxidation in the flue gas and the active oxygen species is crucial for the oxidation reaction. Nevertheless, the detailed mercury/oxygen reaction mechanism on the BiOBr surface is still unclear and the effect of oxygen vacancy on the reaction process was worth in-depth study. Herein, the oxygen-involved Hg0 oxidation mechanism on the BiOBr surface and the effect of oxygen vacancy were deeply understood by quantum chemistry calculation. The reactivity difference of chemisorbed and lattice oxygen species was revealed. Therein, the chemisorbed oxygen derived from the activation of gaseous O2 molecule showed the higher reactivity than lattice oxygen for Hg0 oxidation reaction. The introduction of oxygen vacancy could promote O2 activation and simultaneously reduce the energy barrier of HgO formation from 1.76 to 1.45 eV. In addition, the energy barrier of oxygen vacancy recovery (0.05 eV) was higher than that of O2 activation (0.02 eV), indicating that the initial oxygen vacancy could be retained in the reaction process. This work deeply revealed the mercury/oxygen reaction mechanism over BiOBr photocatalyst and provided the guidance for the catalyst design in the photocatalytic Hg0 oxidation.

Photo-Induced Generation of Oxygenated Quaternary Centers via EnT Enabled Singlet O2 Addition to C3-Maleimidated Quinoxaline: A Reagent-Less Approach.

Demonstrated here is an external photo-sensitizer-free (auto-sensitized) singlet oxygen-enabled solvent-dependent tertiary hydroxylation and aryl-alkyl spiro-etherification of C3-maleimidated quinoxalines. Such "reagent-less" photo-oxygenation at Csp3-H and etherification involving Csp3-H/Csp2-H are unparalleled. Possibly, the highly π-conjugated N-H tautomer allows the substrate to get excited by irradiation, and subsequently, it attains the triplet state via ISC. This excited triplet-state sensitized molecule then transfers its energy to a triplet-state oxygen (3O2) generating reactive singlet oxygen (1O2) for hydroxylation and spirocyclization depending on the solvent used. In HFIP, the generated alkoxy radical accepts a proton via HAT giving hydroxylated product. In contrast, in an aprotic PhCl it underwent a radical addition at the ortho-position of the C2 aryl to provide spiro-ether. An unprecedented orthogonal spiro-etherification was observed via the displacement of o-substitutents for ortho-substituted substrates. The order of ipso substitution follows the trend-OMe > -OEt > -F > -H > -Cl > -Br. Both these oxygenation reactions can be carried out with nearly equal ease using direct sunlight without the requirement of any elaborate reaction setup. Demonstration of large-scale synthesis and a few interesting transformations have also been realized. Furthermore, several insightful control experiments and quantum chemical computations were performed to unravel the mechanism.

N2 cleavage by silylene and formation of H2Si(μ-N)2SiH2

Fixation and functionalisation of N2 by main-group elements has remained scarce. Herein, we report a fixation and cleavage of the N ≡ N triple bond achieved in a dinitrogen (N2) matrix by the reaction of hydrogen and laser-ablated silicon atoms. The four-membered heterocycle H2Si(μ-N)2SiH2, the H2SiNN(H2) and HNSiNH complexes are characterized by infrared spectroscopy in conjunction with quantum-chemical calculations. The synergistic interaction of the two SiH2 moieties with N2 results in the formation of final product H2Si(μ-N)2SiH2, and theoretical calculations reveal the donation of electron density of Si to π* antibonding orbitals and the removal of electron density from the π bonding orbitals of N2, leading to cleave the non-polar and strong NN triple bond.

Theoretical spectroscopic electronic elucidation with polar and non-polar solvents (IEFPCM model), molecular docking and molecular dynamic studies on bendiocarb -antiallergic drug agent

The potential antiallergic compound bendiocarb was studied theoretically by polar and non-polar (six) solvents (IEFPCM model) by quantum chemical DFT calculations. Rotational potential barriers were estimated for water, acetone, THF, chloroform, benzene, and toluene are 12.62, 12.45, 12.55, 12.69, 12.74 and 12.57 kJ/mol. Analysis was done on the variations in geometries brought about by the electronegative atom’s replacement as well as the intramolecular interactions caused by the benzodioxol ring. NBO analysis as well as geometrical equilibrium have been performed using DFT employing the B3LYP technique with the 6-311G++(d, p) basis set. In vibrational analysis, CH stretching of bendiocarb was theoretically estimated within 3207–3181 cm−1 (gas), 3207–3184 cm−1 (acetone and water), 3207–3182 cm−1 (benzene and toluene), and 3207–3183 cm−1 (chloroform and THF) due to CN vibrations is challenging because of the overlapping and appearance of multiple bands in similar regions. Energy gaps for Bendiocarb in different solvents for water (5.719 eV), acetone (5.716 eV), THF (5.722 eV), chloroform (5.724 eV), benzene (5.731 eV), toluene (5.730 eV), and gas phase (5.746 eV). Lipinski’s rule of five was exploited to guide research on drug similarity and oral action. Using docking simulation against the antiallergic protein, a binding energy of −5.22 kcal/mol, 4P5M has been examined. The protein 4P5M exhibits variations between 0.12 and 0.155 nm for RMSD in molecular dynamic simulation. One can infer information about the pharmacokinetic features of the target molecule from its drug similarity, absorption, distribution, metabolism, excretion, and toxicity properties.

Experiment-free multiscale simulation of residual deformation in non-crimp fabric composites

The effectiveness of multiscale modeling, consisting of quantum chemical reaction path calculations, molecular dynamics simulations, and microscopic and macroscopic finite element analysis, developed for the process-induced residual deformation of carbon fiber-reinforced plastic was experimentally evaluated. In this study, non-crimp fabric composites with an arbitrary thermoset resin were fabricated, and process-induced deformation was measured. Process-induced residual deformations due to cure shrinkage and thermal shrinkage of cross-ply laminates under conditions similar to those in the experiment were predicted using multiscale modeling. The predicted deformation shape and maximum deformation value are in good agreement with the fabrication experiments. Multiscale modeling can consider the combined factors that affect the deformation (i.e. size, elastic moduli, coefficient of thermal expansion, and fiber volume fraction) and can provide important knowledge regarding the development of high-performance composite structures and for stable manufacturing.

K2(Bi@Pd12@Bi20)]4-: An Endohedral Inorganic Fullerene with Spherical Aromaticity.

Inorganic fullerene clusters have attracted widespread attention due to their highly symmetrical geometric structures and intrinsic electronic properties. However, cage-like clusters composed of heavy metal elements with high symmetry are rarely reported, and their synthesis is also highly challenging. In this study, we present the synthesis of a [K2(Bi@Pd12@Bi20)]4- cluster that incorporates a {Bi20} cage with pseudo-Ih symmetry, making it the largest main group metal cluster compound composed of the bismuth element to date. Magnetic characterization and theoretical calculations suggest that the spin state of the overall cluster is a quartet. Quantum chemical calculations reveal that the [Bi20]3- cluster has a similar electronic configuration to C606- and the [Bi@Pd12@Bi20]6- cluster exhibits a unique open-shell aromatic character.

Systematic Study of Different Types of Interactions in α-, β- and γ-Cyclodextrin: Quantum Chemical Investigation

In this work, comprehensive ab initio quantum chemical calculations using the DFT level of theory were performed to characterize the stabilization interactions (H-bonding and hyperconjugation effects) of two stable symmetrical conformations of α-, β-, and γ-cyclodextrins (CDs). For this purpose, we analyzed the electron density using “Atom in molecules” (AIM), “Natural Bond Orbital” (NBO), and energy decomposition method (CECA) in 3D and in Hilbert space. We also calculated the H-bond lengths and OH vibrational frequencies. In every investigated CD, the quantum chemical descriptors characterizing the strength of the interactions between the H-bonds of the primary OH (or hydroxymethyl) and secondary OH groups are examined by comparing the same quantity calculated for ethylene glycol, α-d-glucose (α-d-Glcp) and a water cluster as reference systems. By using these external standards, we can characterize more quantitatively the properties of these bonds (e.g., strength). We have demonstrated that bond critical points (BCP) of intra-unit H-bonds are absent in cyclodextrins, similar to α-d-Glcp and ethylene glycol. In contrast, the CECA analysis showed the existence of an exchange (bond-like) interaction between the interacting O…H atoms. Consequently, the exchange interaction refers to a chemical bond, namely the H-bond between two atoms, unlike BCP, which is not suitable for its detection.

Dauresorcinols A and B, two pairs of merosesquiterpenoid enantiomers with new carbon skeletons from Rhododendron dauricum

Five pairs of new merosesquiterpenoid enantiomers, named dauresorcinols A−E (1−5), were isolated from the leaves of Rhododendron dauricum. Their structures were elucidated by comprehensive spectroscopic data analysis, quantum chemical calculations, Rh2(OCOCF3)4-induced ECD, and single-crystal X-ray diffraction analysis. Dauresorcinols A (1) and B (2) possess two new merosesquiterpene skeletons bearing an unprecedented 2,6,7,10,14-pentamethyl-11-oxatetracyclo[8.8.0.02,7.012,17]octadecane and a caged 15-isohexyl-1,5,15-trimethyl-2,10-dioxatetracyclo[7.4.1.111,14.03,8]pentadecane motif, respectively. Plausible biosynthetic pathways of 1−5 are proposed involving key oxa-electrocyclization and Wagner−Meerwein rearrangement reactions. (+)/(−)-1 and 3−5 showed potent α-glucosidase inhibitory activity, 3 to 22 times stronger than acarbose, an antidiabetic drug targeting α-glucosidase. Docking results provide a basis to design and develop merosesquiterpenoids as potent α-glycosidase inhibitors.

Near-infrared two-photon absorption and excited state dynamics of a fluorescent diarylethene derivative.

Near-infrared two-photon absorption and excited state dynamics of a fluorescent diarylethene (fDAE) derivative were investigated by time-resolved absorption and fluorescence spectroscopies. Prescreening with quantum chemical calculation predicted that a derivative with methylthienyl groups (mt-fDAE) in the closed-ring isomer has a two-photon absorption cross-section larger than 1000 GM, which was experimentally verified by Z-scan measurements and excitation power dependence in transient absorption. Comparison of transient absorption spectra under one-photon and simultaneous two-photon excitation conditions revealed that the closed-ring isomer of mt-fDAE populated into higher excited states deactivates following three pathways on a timescale of ca. 200fs: (i) the cycloreversion reaction more efficient than that by the one-photon process, (ii) internal conversion into the S1 state, and (iii) relaxation into a lower state (S1' state) different from the S1 state. Time-resolved fluorescence measurements demonstrated that this S1' state is relaxed to the S1 state with the large emission probability. These findings obtained in the present work contribute to extension of the ON-OFF switching capability of fDAE to the biological window and application to super-resolution fluorescence imaging in a two-photon manner.