Density Functional Theory Method Research Articles

Electrochemical Oxidation of UV Filters: A First-Principles Molecular Dynamics Study.

A theoretical model is proposed to study the oxidation mechanisms of the organic UV filters BP3 and BP4 during electrochemical water treatment utilizing Car-Parrinello molecular dynamics. Factors such as the amount of solvent to be included and how to design the system with the least possible intervention are discussed. The proposed model consist of the optimization of the geometries by density functional theory methods, the equilibration of the structure immersed in a water box, the inclusion of the reactive species, and the analysis of the reaction energies of each reaction pathway. The ab-initio molecular dynamics simulations lead to several products, and some trends can be identified, in accordance with the well-known reactivity rules of organic chemistry. The products proposed in this work are intermediates in longer oxidative pathways.

DFT Study in an Asymmetric Organocatalytic Homologation: Mechanism, Origin of Stereoselectivity

AbstractIn the present study, the mechanism, origin of stereoselectivity, and role of catalyst of the organocatalytic homologation reaction between phenyl boronate and trifluoromethyl diazomethane have been theoretically investigated using the density functional theory (DFT) method. Based on the calculations, the in situ generated ethanol plays significant role in triggering the reaction. Moreover, the [1,2]‐boron migration process is the stereoselectivity‐determining step, with the S‐configured isomer being generated predominantly. In addition, NCI analysis is performed to disclose the origin of stereoselectivity.

Force-Assisted Orbital Crossing in Mechanochemical Oxirane Ring Opening.

Polymer mechanochemistry induces chemical reactivity by applying a directed force, which can lead to unexpected reaction mechanisms. Strained cyclic molecules are often used in force-sensitive motifs because of the strong force coupling of ring-opening reactions. In this computational study, the force dependence of the ring-opening reactions of oxirane will be investigated. Density functional theory and multireference methods were used to investigate the electronic character of both symmetry-allowed and symmetry-forbidden reactions. In the latter case, an orbital crossing occurs during the reaction course, forcing the Woodward-Hoffmann-forbidden reaction to proceed via a diradical pathway. The performance of broken-symmetry density functional theory is evaluated and compares well to high-accuracy CASPT2, MRCI, and ic-MRCC computations. Due to the high ring strain, the barrier heights of both ring-opening reactions are steeply reduced by the application of an external force. Furthermore, the use of unsaturated linkers was shown to yield a significant reduction of the barrier heights, explaining previous experimental findings. Finally, we show through analysis of the PES topology how the external force transforms characteristic points such as saddle points and bifurcations, providing insights into force-dependent mechanism changes.

Transition Metals Doped into g-C3N4 via N,O Coordination as Efficient Electrocatalysts for the Carbon Dioxide Reduction Reaction.

The electrochemical carbon dioxide reduction reaction (CO2RR) is a potential and efficient method that can directly convert CO2 into high-value-added chemicals under mild conditions. Owing to the exceptionally high activation barriers of CO2, catalysts play a pivotal role in CO2RR. In this study, the transition metal (TM = Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, and Zn) is doped into g-C3N4 with a unique N,O-coordination environment, namely, TM-N1O2/g-C3N4. Herein, the catalytic performance and reaction mechanism for the CO2RR on TM-N1O2/g-C3N4 are systematically investigated by density functional theory methods. Especially, through the calculation of ΔG*H and ΔG*COOH/ΔG*OCHO, the catalysts with preference for the CO2RR over the hydrogen evolution reaction (HER) are selected for further study. Furthermore, Gibbs free energy computation results of each elementary step for the CO2RR on these catalysts indicate that Ti-N1O2/g-C3N4 has significant catalytic activity and selectivity for reducing CO2 to methanol (CH3OH) with the limiting potential (UL) of -0.55 V. Finally, through frontier molecular orbital theory and charge transfer analyses, the introduction of the O atoms illustrates that it is instrumental in regulating the electron distribution of the catalytic active site, thereby improving the catalytic performance. This work provides insight into the design of single-atom catalysts with unique coordination structures for the CO2RR.

Synthesis of New 1,4-Naphthoquinone Fluorosulfate Derivatives and the Study of Their Biological and Electrochemical Properties

This study presents the synthesis of new fluorosulfate derivatives of 1,4-naphthoquinone by the SuFEx reaction. Anticancer properties of obtained compounds were studied on PC-3 (prostate adenocarcinoma), SKOV-3 (ovarian cancer), MCF-7 (breast cancer), and Jurkat cell lines. All the studied compounds showed higher cytotoxic effects than Cisplatin. The DFT method was applied to determine the electronic structure characteristics of 1,4-naphthoquinone derivatives associated with cytotoxicity. A method of determination of 2,3-dichloro-1,4-naphthoquinone (NQ), 3-chloro-2-((4-hydroxyphenylamino)-1,4-naphthoquinone (NQ1), and 4-((3-chloro-1,4-naphthoquinon-2-yl)amino)phenyl fluorosulfate (NQS) in a pharmaceutical substance using an impregnated graphite electrode (IMGE) was developed. The morphology of the IMGE surface was studied using scanning electron microscopy (SEM). The electrochemical behavior of NQ, NQ1, and NQS was studied by cyclic voltammetry (CV) in 0.1 M NaClO4 (96% ethanol solution) at pH 4.0 in a potential range from −1 to +1.2 V. Electrochemical redox mechanisms for the investigated compounds were proposed based on the determining main features of the electrochemical processes. Calibration curves were obtained by linear scan voltammetry in the first derivative mode (LSVFD) with the detection limit (LOD) 7.2 × 10−6 mol·L−1 for NQ, 8 × 10−7 mol·L−1 for NQ1, and 8.6 × 10−8 mol·L−1 for NQS, respectively.

Enhanced Optoelectronic Properties of Few-Layer Graphene for Transparent Conducting Electrode: Density Functional Theory Perspective

The challenge in utilizing graphene for transparent conducting electrodes (TCE) is to increase the electronic conductivity without significantly decreasing its transmittance. Forming a few-layer structure is an effective way to modify the electronic properties of graphene. In this paper, we used the density functional theory (DFT) method to study the effect of the number of few-layer graphene (FLG) layers on the electronic conductivity and transmittance. We find that increasing the number of layers leads to an increase in electronic conductivity. On the other hand, a decrease in transmittance was observed when the number of layers was increased. Based on the performance analysis, the optoelectronic properties of FLG with more than three layers were better than those of conventional TCE materials, such as indium tin oxide (ITO) and aluminium-doped zinc oxide (AZO). The calculated electronic conductivity and transmittance of a 4-layer FLG were 2.23×1019 (1/Ω.m.s) and 95.69%, respectively. This finding can be used as guidance for experimental research in developing graphene-based TCE materials. Keywords: density functional theory, electronic conductivity, few-layer graphene, transmittance

Potassium and Boron Co-Doping of g-C3N4 Tuned CO2 Reduction Mechanism for Enhanced Photocatalytic Performance: A First-Principles Investigation

Graphitic phase carbon nitride (g-C3N4, abbreviated as CN) can be used as a photocatalyst to reduce the concentration of atmospheric carbon dioxide. However, there is still potential for improvement in the small band gap and carrier migration properties of intrinsic materials. K-B co-doped CN (KBCN) was investigated as a promising photocatalyst for carbon dioxide reduction via the Density Functional Theory (DFT) method. The electronic and optical properties of CN and KBCN indicate that doping K and B can improve the catalytic performance of CN by promoting charge migration and separation. In terms of the Gibbs free energy change, the CO2 reduction reaction catalysed by KBCN results in CH3OH, and its optimal pathway is CO2 → *CO2 → *COOH → CO → *OCH → HCHO → *OCH3 → CH3OH. Compared with CN, the doping elements K and B shift the rate-determining step from CO2 → *CO2 to *CO2 → *COOH. The K and B elements co-doping tuned the charge distribution between the catalyst and the adsorbate and reduced the Gibbs free energy of the rate-determining step from 1.571 to 0.861 eV, suggesting that the CO2 reduction activity of KBCN is superior to that of CN. Our work provides useful insights for the design of metallic–nonmetallic co-doped CN for photocatalytic CO2 reduction (CO2PR) reactions.

Homobimetallic Ruthenium(II) Complexes Catalysed Selective Transfer Hydrogenation of Aldehydes in Water.

Herein we report chemoselective transfer hydrogenation (TH) of aldehydes in aqueous medium using a series of homobimetallic Ru(II) catalysts. Two homobimetallic complexes (Ru1 and Ru3) and one monometallic complex (Ru2) have been employed in the catalytic reduction of aldehydes. Bimetallic complex [(p-cymene)2(RuCl)2L3] (Ru3) is obtained from the reaction of Schiff base ligand 2,2'-((1E,1'E)-((3,3',5,5'-tetraisopropyl-[1,1'-biphenyl]-4,4'diyl)bis(azaneylylidene))bis(methaneylylidene))bis(4-bromophenol) (H2L3) and characterized by various spectroscopic and analytical techniques. The use of formic acid/formate buffer as the hydride source and a catalyst loading of 0.01 mol% of Ru1 or Ru3 resulted in the conversion of various aldehydes to the corresponding alcohols in good to excellent yield. This method is very efficient for selective reduction of aldehydes in the presence of other reducible functional groups. A loading of 0.0001 mol% of Ru1 catalyst is sufficient to achieve a turnover frequency (TOF) of 5.5× 105 h-1. Furthermore, the catalyst can been recycled and reused for six consecutives cycles without sacrificing the efficiency. A comparison of results obtained between bimetallic and monometallic complexes offers valuable insights into the distinct reactivity patterns of the bimetallic complexes, presumably originating from a cooperative effect. We have explored the mechanistic pathway using DFT methods.

From germolane to germylenes: a theoretical DFT study of thermal decomposition pathways and reactivity

This study addresses the often-overlooked chemistry of germanium compared to the extensively researched carbon and silicon. Using advanced DFT methods, we investigated the thermal decomposition of germolane (germacyclopentane). The suggested mechanisms include a 1,2-H shift and 1,1-H2 elimination to form a pentacyclic germylene (1λ2-germolane). The other pathway involves a stepwise [3 + 2] cycloreversion to form a diradical followed by ethene and germirane. Under M062X/def2-TZVP level of theory, the activation barriers in terms of Gibbs energy (ΔG ‡ 298) for the 1,2-H shift and 1,1-H2 elimination pathways were 240.9 and 236.3 kJ/mol, respectively. The reaction energies (ΔG° 298) for the initiative steps in the 1,2-H shift and 1,1-H2 elimination mechanisms were 102.9 and 96.2 kJ/mol, indicating a thermodynamic and kinetic competition between the two routes. Temperature dependence analysis from 300 to 1200 K reveals that the 1,1-H2 elimination dominates at higher temperatures and is expected to become spontaneous above 1000 K.

Unveiling the Novel Mechanistic Insights and Role of Base in Zn‐Catalyzed Csp–Csp2 Cross‐Coupling Reaction

ABSTRACTA detailed mechanistic investigation of the Zn (II)‐catalyzed Csp–Csp2 (Sonogashira‐type) cross‐coupling reaction is reported herein, using the Density Functional Theory method. The present study unveiled an unconventional non‐redox mechanism for Zn‐catalyzed cross‐coupling reaction, where the oxidation state of Zn remains intact throughout the catalytic cycle. Our study further revealed the significant role of the base in controlling the feasibility of cross‐coupling reactions that are catalyzed by electron‐deficient metal centers. Our study indicates that K3PO4 acts as an ancillary ligand (Lewis base) for the electron‐deficient Zn (II) catalytic center rather than as a proton abstractor for the nucleophilic coupling partner (phenylacetylene) in this reaction. The active catalyst was identified to be a four‐coordinate bis‐DMEDA Zn (II) complex. The mechanism proceeds via the initial activation of the nucleophilic coupling partner (phenylacetylene), followed by the electrophilic coupling partner (organic halide) activation liberating the cross‐coupled product by a concerted nucleophilic substitution pathway. The turn‐over limiting step was identified to be the activation of the electrophilic coupling partner. The activation barrier obtained for the reaction, 31.0 kcal/mol concords well with experimental temperature requirements (125°C). The coordination by base is found to stabilize the rate‐determining intermediates and transition states involved in the reaction. The mechanistic insights gained from this study could aid in the rational design and development of sustainable cross‐coupling reactions using zinc as the catalyst.

A theoretical exploration of the influence of oxygen element on photo-induced behaviours for flavonoid derivatives with Me2N substitution

The exceptional properties of flavonoid derivatives are a great inspiration for our research. In the fields of medicine and chemistry, flavonoids present the superior ability to enhance non-specific immune function and humoral immune response. Herein, our primary focus lies in the theoretical exploration of the photo-induced characteristics exhibited by the novel Me2N-substituted flavonoid (MEF) fluorophore. Considering the potential photo-induced properties associated with chalcogen atomic electronegativity in photoexcitation, using density functional theory (DFT) and time-dependent DFT (TDDFT) methods, we primarily pay attention to probing into photo-induced hydrogen bonding effects and concomitant charge reorganisation processes. Given changes in chemical geometries and the infrared (IR) vibrational shifts of three MEF derivatives (MEF-O, MEF-S and MEF-Se), we confirm that intramolecular hydrogen bonding should be strengthened by facilitating excited state intramolecular proton transfer (ESIPT) reactions. Particularly, the redistribution of charge further shows the ESIPT behaviour. By constructing potential energy curves (PECs) and identifying transitional reaction state (TS) structures, we ultimately validate the electronegativity-dependent ESIPT mechanisms of chalcogen elements in MEF derivatives. We sincerely hope that our work can accelerate future development and applications of flavonoid derivatives.

Molecular designing and structural tuning of derivatives of 4,7-divinyl-1H-benzo[b]silole for optoelectronic properties using DFT and TD-DFT methods

Molecular designing and structural tuning of derivatives of 4,7-divinyl-1H-benzo[b]silole for optoelectronic properties using DFT and TD-DFT methods

Machine Learning Model for the Prediction of Hubbard U Parameters and Its Application to Fe-O Systems.

Without incurring additional computational cost, the Hubbard model can prevalently address the electron self-interaction problems of the local or semilocal exchange-correlation functions within density functional theory. However, determining the value of the Hubbard parameter, U, promptly, efficiently, and accurately has been a long-standing challenge. Here, we develop a method for predicting the Hubbard U of iron oxides by establishing a potential relationship through machine learning fitting of structural fingerprints and the U evaluated by the linear response-constrained density functional theory method. This method performs well in calculating the properties of wüstite, hematite, and magnetite, aligning with experimental measurements or more costly hybrid functional results. Using this method, we redefine the convex hulls of the Fe-O system at 0, 50, and 100 GPa; the obtained results are in good agreement with experimental observations. We also provide insights into the debates surrounding the high-pressure phases of Fe2O3 and Fe3O4.

Copper‐Catalyzed Simultaneous Dehydrogenation and Enantioselective Allylic Oxidation of Alkanes Using Chiral Heterogeneous Ligands to Produce Allylic Esters and Theoretical Investigation of the Mechanism

ABSTRACTA new asymmetric method for preparing chiral allylic esters from various unactivated alkanes through simultaneous dehydrogenation and enantioselective allylic oxidation in the presence of chiral heterogeneous oxazoline–based ligands is reported for the first time. For this purpose, chiral amino oxazoline ligands, synthesized from chiral amino alcohols and cyanogen bromide, were immobilized on modified MCM‐41 mesoporous silica. These chiral heterogeneous ligands were then characterized using Fourier‐transform infrared spectroscopy, thermogravimetric analysis, X‐ray diffraction, scanning electron microscopy, energy‐dispersive X‐ray spectroscopy, and Brunauer–Emmett–Teller and Barrett–Joyner–Halenda techniques. Finally, copper complexes of these ligands were employed in the synthesis of chiral allylic esters, achieving the best result with a 90% yield and 77% enantiomeric excess. Evaluation of the recyclability of the chiral heterogeneous catalysts revealed that they can be reused three times without a noticeable reduction in the results. In addition, the mechanism of formation of the chiral allylic esters was investigated using a density functional theory method.

Enabling Visible Light Sensitization of YbIII, NdIII and ErIII in Dimeric LnIII/GaIII Metallacrowns through Functionalization with RuII Complexes for NIR‐II Multiplex Imaging

AbstractMultiplex imaging in the second near‐infrared window (NIR‐II, 1000–1700 nm) provides exciting opportunities for more precise understanding of biological processes and more accurate diagnosis of diseases by enabling real‐time acquisition of images with improved contrast and spatial resolution in deeper tissues. Today, the number of imaging agents suitable for this modality remains very scarce. In this work, we have synthesized and fully characterized, including theoretical calculations, a series of dimeric LnIII/GaIII metallacrowns bearing RuII polypyridyl complexes, LnRu‐3 (Ln=YIII, YbIII, NdIII, ErIII). Relaxed structures of YRu‐3 in the ground and the excited electronic states have been calculated using dispersion‐corrected density functional theory methods. Detailed photophysical studies of LnRu‐3 have demonstrated that characteristic emission signals of YbIII, NdIII and ErIII in the NIR‐II range can be sensitized upon excitation in the visible range through RuII‐centered metal‐to‐ligand charge transfer (MLCT) states. We have also showed that these NIR‐II signals are unambiguously detected in an imaging experiment using capillaries and biological tissue‐mimicking phantoms. This work opens unprecedented perspectives for NIR‐II multiplex imaging using LnIII‐based molecular compounds.

Synthesis, Characterization and Theoretical Calculations of a Novel Organic Corrosion Inhibitor for Mild Steel

Abstract Q235 is a carbon structural steel sheet that is widely employed in the construction and engineering sectors due to its cost-effectiveness and durability. Here, A novel organic corrosion inhibitor for Q235 sheet based on furan-thiadiazole (5-amino-2-furyl-1, 3, 4-thiadiazole, AFTZ) has been designed, successfully synthesized and characterized in this work. Results of experimental data and theoretical calculations indicated that AFTZ has an excellent corrosion inhibition property for Q235 sheet in 1.0 mol/L H2SO4. The highest efficiency of corrosion inhibitor tested by weight loss experiment can reach 91.17% and electrochemical experiment was up to 95.85% when AFTZ increased to 2 mg/mL in H2SO4 solution. The morphological characteristics of Q235 sheets with or without corrosion inhibitor showed that the target inhibitor has good corrosion inhibition performance. Quantum computation using the DFT method of guassian09 software and FT-IR date attribute the excellent corrosion inhibition to the film-forming nature of AFTZ.

DFT molecular modeling investigation and optical properties characterization of CR-39 films

Abstract In this study, the structural and optical characteristics of CR-39 films were investigated both theoretically and experimentally. A number of parameters for the CR-39 structure were theoretically computed by the density functional theory (DFT) method at B3LYP/6–311 G (p, d) level of theory. FTIR and UV/Vis spectroscopy were used to determine the structure compositions and the optical parameters of the CR39 film, respectively. The computed and experimental findings show good concordance. It was found that the CR-39 compound’s total energy, dipole moment, and energy difference between the LUMO and HOMO were, respectively, −994.575 a.u., 2.772 Debye, and 7.045 eV. The MEP showed a progressive change in color, with blue signifying a low electron density and red denoting a high electron density linked to nucleophilic reactivity and electrophilic reactivity. Further, the positive charges on all hydrogen atoms, which range from 0.16506 to 0.22032 a.u., imply that they are acceptors. The positive H34 is strongly produced when electrons from the negatively charged C17 are taken up. Because they are donor atoms, part of the carbon atoms in the structure is negative, while the other atoms are positive. The highest electronegative atoms H23 and H24 were substituted, resulting in the extremely negative carbon atom C7 (−0.32436). According to the experimentally determined absorption coefficient and energy gap values, CR-39 films may find application in optoelectronic devices.

Implementation of time-dependent Hartree–Fock in real space

Abstract Time-dependent Hartree–Fock (TDHF) is one of the fundamental post-Hartree–Fock (HF) methods to describe excited states. In its Tamm-Dancoff form, equivalent to Configuration Interaction Singles, it is still widely used and particularly applicable to big molecules where more accurate methods may be unfeasibly expensive. However, it is rarely implemented in real space, mostly because of the expensive nature of the exact-exchange potential in real space. Compared to widely used Gaussian-type orbitals (GTO) basis sets, real space often offers easier implementation of equations and more systematic convergence of Rydberg states, as well as favorable scaling, effective domain parallelization, flexible boundary conditions, and ability to treat model systems. We implemented TDHF in the Octopus real-space code as a step toward linear-response hybrid time-dependent density-functional theory (TDDFT), other post-HF methods, and ensemble density-functional theory methods involving exact exchange. Calculation of HF’s non-local exact exchange is very expensive in real space. We overcome this limitation with Octopus’ implementation of Adaptively Compressed Exchange, and find the appropriate mixing scheme and starting point to complete the ground-state calculation in a practical amount of time, and thus enable TDHF. We compared our results to those from GTOs on a set of small molecules and confirmed close agreement of results, though with larger deviations than in the case of semi-local TDDFT. We find that convergence of TDHF demands a finer real-space grid than semi-local TDDFT. We also present the subtleties in benchmarking a real-space calculation against GTOs, relating to Rydberg and vacuum states.

Slope of the Delocalization Function Is Proportional to Analytical Hardness.

Conceptual Density Functional Theory (CDFT) has been extended beyond its traditional role in elucidating chemical reactivity to the development of density functional theory methods, e.g., the investigation of the delocalization error. This delocalization error causes the dependence of the energy on the number of electrons (N) to deviate from its exact piecewise linear behavior, an error which is the basis of many well-known limitations of commonly used density-functional approximations (DFAs). Following our previous work on the analytical hardness η± for pure functionals, we extend its application to hybrid and range-separated functionals. A comparison is made between the analytical hardness and the slope of the delocalization function introduced by Hait and Head-Gordon. Our results show that there is a linear relationship between its slope and the analytical hardness. An approximate scheme is presented to construct the energy vs N curve without fractional occupation number calculations. The extension to densities is discussed.

Exploring nonlinear correlations among transition metal nanocluster properties using deep learning: a comparative analysis with LOO-CV method and cosine similarity

A novel approach is introduced for the rapid and accurate correlation analysis of nonlinear properties in Transition Metal (TM) clusters utilizing the Deep Leave-One-Out Cross-Validation technique. This investigation demonstrates that the Deep Neural Network (DNN)-based approach offers a more efficient predictive method for various properties of fourth-row TM nanoclusters compared to conventional Density Functional Theory methods, which are computationally intensive and time-consuming. The feature space, also known as descriptors, is established based on a broad spectrum of electronic and physical characteristics. Leveraging the similarities among these clusters, the DNN-based model is employed to explore the correlations among TM cluster properties. The proposed method, in conjunction with cosine similarity, achieves remarkable accuracy up to 10-9 for predicting total energy, lowest vibrational mode, binding energy, and HOMO-LUMO energy gap of TM2, TM3, and TM4nanoclusters. By analyzing correlation errors, the most closely coupled TM clusters are identified. Notably, Mn and Ni clusters exhibit the highest and lowest levels of energy coupling with other TMs, respectively. Generally, energy prediction for TM2, TM3, and TM4clusters exhibit similar trends, while an alternating behavior is observed for vibrational modes and binding energies. Furthermore, Ti, V, and Co demonstrate the highest binding energy correlations with TM2, TM3, and TM4sets, respectively. Regarding energy gap predictions, Ni exhibits the strongest correlation in the smallest TM2clusters, while Cr shows the highest dependence in TM3and TM4sets. Lastly, Zn displays the largest error in HOMO-LUMO energy gap across all sets, indicating distinctive independent energy gap characteristics.