Exchange-correlation Contributions Research Articles

Pentafluoroorthotellurate Uncovered: Theoretical Perspectives on an Extremely Electronegative Group.

The pentafluoroorthotellurate group (-OTeF5, teflate) exhibits high electron-withdrawing properties. Indeed, it is often used as a bulky substitute for fluoride due to its high chemical stability and larger size, which reduces its tendency to act as a bridging ligand. These characteristics make it a valuable ligand in synthetic chemistry, facilitating the preparation of molecular structures analogous to polymeric fluoride-based compounds. In this study, we explore the electronic structure of the teflate group by using advanced Quantum Chemical Topology (QCT) methods to better understand its bonding nature and compare its group electronegativity with that of the halogens. For that, we examine XOTeF5 systems (X = F, Cl, Br, I) and decompose X-OTeF5 interactions into classical (ionic) and exchange-correlation (covalent) contributions by using interacting quantum atoms (IQA) energy decomposition scheme. We also conduct a detailed analysis of electron distribution by utilizing the statistical framework of electron distribution functions (EDFs) and examine the electron localization function (ELF), electron density, and reduced density gradient scalar functions, as well as delocalization indices and QTAIM charges. The results show that the electron-withdrawing properties of the teflate group are comparable to those of fluorine, albeit slightly lower. Moreover, its internal bonding is primarily ionic. Additionally, we compare -OTeF5 with other O-donor groups, demonstrating that the electron-withdrawing properties within OEF5 (E = S, Se, Te) systems are nearly identical, and these groups show a higher group electronegativity than OCF3, OC(CF3)3, and OC6F5.

Exploring the Reactivity of Donor-Acceptor Systems through a Combined Conceptual and Constrained DFT Approach.

In the context of the conceptual density functional theory (cDFT) and based on the computational efficiency of the constrained DFT (CDFT), we demonstrate that chemical reactivity can be governed by the difference between the local interacting chemical potentials of the reactants (referred as Edual), in agreement with Sanderson's equalization principle. In a proof-of-concept study, we investigated illustrative examples involving typical non-covalent donor-acceptor systems and reactive systems are provided. For the selected systems, our approach reveals significant mimicking between Edual and the DFT-computed intermolecular interaction energy profiles. We further evaluate the influence of the Coulomb and exchange-correlation contributions in Edual. These latter results suggest that numerous potential energy surfaces of clusters can be explored using a Sanderson-like model only based on classical interactions between molecular orbitals domains. To conclude, this study achieved a deeper understanding of the principles of cDFT and assessed, in a wider context, its efficiency in predicting the chemical reactivity.

Giant Spatial Redistribution of Electrons in a Wide Quantum Well Induced by Quantizing Magnetic Field

In samples of field-effect transistors based on GaAs/AlGaAs heterostructures with an electron system in a single 50-nm-wide GaAs quantum well, a transition stimulated by a quantizing magnetic field has been detected from a bilayer state of the system in zero magnetic field to a single-layer state when only the lowest Landau level is filled. In contrast to the results for the 60-nm-wide quantum well obtained in [S. I. Dorozhkin, A. A. Kapustin, I. V. Fedorov, V. Umansky, and J. H. Smet, Phys. Rev. V 102, 235307 (2020)], the single-layer state is observed not only in incompressible quantum Hall effect states of the electron system at filling factors of 1 and 2, but also in compressible states between these filling factors. The spatial location of the single-layer system in the quantum well has been established; it appears to be independent of the electron distribution over the layers in a low magnetic field. A possible qualitative explanation for this observation has been proposed. The detected transition is supposedly due to the negative compressibility of two-dimensional electron systems caused by exchange-correlation contributions to the electron−electron interaction.

Ab initio study on structure, IR spectrum, and chemical bonds of AnOSF2 (An=Th, Pa, U, Np, and Pu) molecules

The structures, vibrational spectra, stability, and chemical bonds for a series of AnOSF2 (An=Th-Pu) molecules were studied. Four IR bands were fingerprinted for An–O, An-S, symmetric and antisymmetric F–An–F stretchings, respectively. The Gibbs free energy of the reaction obtained at DLPNOCCSD (T) level suggests an exothermic reaction. Interaction quantum atom (IQA) analyses indicate that the quantum exchange-correlation contributions are great, especially for U-S, Np-S, and Pu-S interactions. The delocalization characters in the U-S and U-O interactions are greatly higher than those of U-F.

Interpreting the different emissive properties of cyclic triimidazole-based CuI and AgI coordination polymers: a QTAIM and IQA study

A QTAIM and IQA investigation on model compounds of two isostructural AgI and CuI coordination polymers (CPs) based on cyclic triimidazole (L), i.e. the [MIL] n 1D double-stranded stair chain and the [MClL] n 3D network (M = Cu, Ag), has allowed light to be shed on the different emissive behaviour associated with the two metal ions. According to a previously reported investigation [Malpicci et al. (2021). Inorg. Chem. Front. 8, 1312–1323], AgI CPs showed both fluorescence and multiple ligand-centred room-temperature phosphorescences, whereas CuI CPs displayed non-thermally equilibrated halogen and metal-to-ligand charge transfer and two ligand-centred phosphorescences, the latter observed only by their selective activation. Analysis of both local and integral QTAIM descriptors, including delocalization indices and source function, of the Ag—N and Cu—N bonds reveals a higher covalent and local character for the latter, explaining the greater metal–ligand electronic communication observed for the Cu compounds. Moreover, IQA investigation shows that the Cu—N bond is characterized by higher interaction energy, due to both higher electrostatic and exchange-correlation contributions. Analysis on the M—X (M = Ag, Cu; X = I, Cl) bonds, also present in these structures, highlights a much higher covalent and local character with respect to the M—N bonds.

Interatomic exchange‐correlation interaction energy from a measure of quantum theory of atoms in molecules topological bonding: A diatomic case

The potential relations between the measure of topological interatomic bonding-integrals of electron density with respect to internuclear axis over the corresponding quantum theory of atoms in molecules (QTAIM)-defined interatomic surface (IAS)-and interatomic exchange-correlation contributions from the interacting quantum atoms approach are discussed. The quantum chemical computations of 38 equilibrium diatomic systems at different levels of theory (HF, MP2, MP4SDQ, and CCSD) are invoked to support abstract considerations. Parameters of excellent correlations between IAS integrals and interatomic exchange-correlation energy are found by the optimization. The performance of these trends depends on the accuracy of the electronic correlation treatment. The resulting trends are a unique feature of equilibrium states, whereas more complicated dependencies are explored for several systems at non-equilibrium conditions. The relations of established trends with other IAS-based estimations of strength of bonding interactions between topological atoms and issues explored for multiatomic systems are briefly discussed.

Feynman Force Analysis of Chemical Processes in Terms of Topological Atomic Contributions

The nature of the chemical bond is analyzed in terms of the atomic contributions to the Feynman forces using the Quantum Theory of Atoms in Molecules and the Interacting Quantum Atoms method. This approach provides a means for quantifying the relationship between the atomic electronic reorganization and the evolution of functional group interactions with the forces exerted on the nuclear framework during a chemical transformation. Using this decomposition scheme, the forces driving a chemical process are locally assigned to atoms or functional group contributions. The interatomic component of the forces can be ascribed as bonding forces; their exchange-correlation and electrostatic contributions reveal the nature of the interactions affecting the forces on the nuclei. This method is used to analyze the chemical interactions involved in the formation of ground and excited state diatomic molecules, the prototropism of formamide, the Diels-Alder cycloaddition of 1,3-butadiene with ethylene, and the Jahn-Teller effect of hydrated transition metal complexes.

Anionic behavior and single-molecule crystal in fullerene confinements: A contribution from DFT energy decomposition and cooperativity analyses.

The recently proposed systems of various anions (A) confined inside C60 , A- @ C60 , which in turn behave as large and stable anions, (A @ C60 )- , can find potential applications in various fields. On the other hand, it has earlier been shown that from the dihalogens (X2 ) encapsulated C60 , X2 @ C60 , only F2 @ C60 can be introduced as a system in which the cage acts as a cation C60 + and interacts with an endohedral anion, F2 - , forming the F2 - @ C60 + as a single-molecule crystal compound. In this work, two density functional theory energy decomposition analysis (EDA) schemes, where in one of them the noninteracting kinetic, electrostatic, and exchange-correlation energies come into play while another scheme, called as EDA-SBL, includes the steric, electrostatic, and quantum effects as essential ingredients (S. Liu, J. Chem. Phys. 2007, 126, 244103), are utilized to find out what energetic components govern the unique characteristics of the (A @ C60 )- and X2 @ C60 confinements. It is shown that the noninteracting kinetic energy and steric energies have important contributions to the total interaction energies for the considered systems. However, there are other confinements for which the electrostatic and exchange-correlation contributions play also imperative roles. Furthermore, we find reasonable correlations between interaction energies and their components as well as the energetic components themselves, leading to an alternative EDA scheme including the noninteracting kinetic, steric, and electrostatic energies for investigations on other endohedral fullerenes. Extending our analyses to large size confinements, Cl- @ Cn with n up to 90 as illustrative examples, the quantitative cooperativity concept is also explored, where the positive and negative cooperativity profiles unveil a specific size of the anionic confinements to form the most stable large anion.

Reactivity of α,ω-Dihydrofluoropolyethers toward OH Predicted by Multiconformer Transition State Theory and the Interacting Quantum Atoms Approach.

We report rate constants for the tropospheric reaction between the OH radical and α,ω-dihydrofluoropolyethers, which represent a specific class of the hydrofluoropolyethers family with the formula HF2C(OCF2CF2)p(OCF2)qOCF2H. Four cases were considered: p0q2, p0q3, p1q0, and p1q1 (pxqy denoting p = x and q = y) with the calculations performed by a cost-effective protocol developed for bimolecular hydrogen-abstraction reactions. This protocol is based on multiconformer transition state theory and relies on computationally accessible M08-HX/apcseg-2//M08-HX/pcseg-1 calculations. Within the protocol's approximations, the results show that (1) the calculated rate constants are within a factor of five of the experimental results (p1q0 and p1q1) and (2) the chain length and composition have a negligible effect on the rate constants, which is consistent with the experimental work. The interacting quantum atoms energy decomposition scheme is used to analyze the observed trends and extract chemical information related to the imaginary frequencies and barrier heights that are key parameters controlling the reactivity of the reaction. The intramolecular exchange-correlation contributions in the reactants and transition states were found to be the dominating factor.

A step in the direction of resolving the paradox of Perdew-Zunger self-interaction correction.

Self-interaction (SI) error, which results when exchange-correlation contributions to the total energy are approximated, limits the reliability of many density functional approximations. The Perdew-Zunger SI correction (PZSIC), when applied in conjunction with the local spin density approximation (LSDA), improves the description of many properties, but overall, this improvement is limited. Here, we propose a modification to PZSIC that uses an iso-orbital indicator to identify regions where local SICs should be applied. Using this local-scaling SIC (LSIC) approach with LSDA, we analyze predictions for a wide range of properties including, for atoms, total energies, ionization potentials, and electron affinities and, for molecules, atomization energies, dissociation energy curves, reaction energies, and reaction barrier heights. LSIC preserves the results of PZSIC-LSDA for properties where it is successful and provides dramatic improvements for many of the other properties studied. Atomization energies calculated using LSIC are better than those of the Perdew, Burke, and Ernzerhof generalized gradient approximation (GGA) and close to those obtained with the strongly constrained and appropriately normed meta-GGA. LSIC also restores the uniform gas limit for the exchange energy that is lost in PZSIC-LSDA. Further performance improvements may be obtained by an appropriate combination or modification of the local scaling factor and the particular density functional approximation.

Exchange-correlation contributions to thermodynamic properties of the homogeneous electron gas from a cumulant Green's function approach

Thermodynamic properties of the interacting homogeneous electron gas are calculated using a finite-temperature cumulant Green's function approach over a broad range of densities and temperatures up to the warm dense matter regime $T\ensuremath{\sim}{T}_{F}$, where ${T}_{F}$ is the Fermi degeneracy temperature. These properties can be separated into independent particle and exchange-correlation contributions, and our focus here is on the latter. Our approach is based on the Galitskii-Migdal-Koltun and electron number sum rules from the finite temperature many-body Green's function formalism, together with an extension of the cumulant Green's function to finite temperature. Previously this approach yielded exchange-correlation energies and potentials in good agreement with quantum Monte-Carlo calculations. Here the method is extended for various thermodynamic quantities including the chemical potential, total energy, Helmholtz free-energy, electronic equation of state, specific heat, and isothermal compressibility, which optionally include spin dependence. We find that the exchange-correlation contributions are weakly varying at low temperature but exhibit significant temperature dependence in the WDM regime, as well as a crossover from exchange- to correlation-dominated behavior. In contrast to the $T={0}^{+}$ limit, we also find that renormalization effects are largely but not completely suppressed at finite temperature. Comparisons with other approaches at various levels of approximation are also discussed.

THz plasma wave instability in field effect transistor with electron diffusion current density

The plasma wave instability in rectangle field effect transistors (FETs) is studied with electron diffusion current density by quantum hydrodynamic model in this paper. General dispersion relation including effects of electrical thermal motion, external friction associated with electron scattering effect, electron exchange-correlation contributions and quantum effects were obtained for rectangle FETs. The electron diffusion current density term is considered for further analysis in this paper. It is found that the quantum effects, the electron diffusion current density and electrical thermal motion enhance the radiation power and frequencies. But the electron exchange-correlation effects and the electron scattering effects reduce the radiation power and frequencies. Results showed that a transistor has advantages for the realization of practical terahertz sources.

When hydroquinone meets methoxy radical: Hydrogen abstraction reaction from the viewpoint of interacting quantum atoms

Interacting Quantum Atoms methodology is used for a detailed analysis of hydrogen abstraction reaction from hydroquinone by methoxy radical. Two pathways are analyzed, which differ in the orientation of the reactants at the corresponding transition states. Although the discrepancy between the two barriers amounts to only 2 kJ/mol, which implies that the two pathways are of comparable probability, the extent of intra-atomic and inter-atomic energy changes differs considerably. We thus demonstrated that Interacting Quantum Atoms procedure can be applied to unravel distinct energy transfer routes in seemingly similar mechanisms. Identification of energy components with the greatest contribution to the variation of the overall energy (intra-atomic and inter-atomic terms that involve hydroquinone's oxygen and the carbon atom covalently bound to it, the transferring hydrogen and methoxy radical's oxygen), is performed using the Relative energy gradient method. Additionally, the Interacting Quantum Fragments approach shed light on the nature of dominant interactions among selected fragments: both Coulomb and exchange-correlation contributions are of comparable importance when considering interactions of the transferring hydrogen atom with all other atoms, whereas the exchange-correlation term dominates interaction between methoxy radical's methyl group and hydroquinone's aromatic ring. This study represents one of the first applications of Interacting Quantum Fragments approach on first order saddle points. © 2018 Wiley Periodicals, Inc.

Ab Initio Optimized Effective Potentials for Real Molecules in Optical Cavities: Photon Contributions to the Molecular Ground State

We introduce a simple scheme to efficiently compute photon exchange-correlation contributions due to the coupling to transversal photons as formulated in the newly developed quantum-electrodynamical density-functional theory (QEDFT).1−5 Our construction employs the optimized-effective potential (OEP) approach by means of the Sternheimer equation to avoid the explicit calculation of unoccupied states. We demonstrate the efficiency of the scheme by applying it to an exactly solvable GaAs quantum ring model system, a single azulene molecule, and chains of sodium dimers, all located in optical cavities and described in full real space. While the first example is a two-dimensional system and allows to benchmark the employed approximations, the latter two examples demonstrate that the correlated electron-photon interaction appreciably distorts the ground-state electronic structure of a real molecule. By using this scheme, we not only construct typical electronic observables, such as the electronic ground-state density, but also illustrate how photon observables, such as the photon number, and mixed electron-photon observables, for example, electron–photon correlation functions, become accessible in a density-functional theory (DFT) framework. This work constitutes the first three-dimensional ab initio calculation within the new QEDFT formalism and thus opens up a new computational route for the ab initio study of correlated electron–photon systems in quantum cavities.

Anion−π Catalysis on Fullerenes

Anion-π interactions on fullerenes are about as poorly explored as the use of fullerenes in catalysis. However, strong exchange-correlation contributions and the localized π holes on their surface promise unique selectivities. To elaborate on this promise, tertiary amines are attached nearby. Dependent on their positioning, the resulting stabilization of anionic transition states on fullerenes is shown to accelerate disfavored enolate addition and exo Diels-Alder reactions enantioselectively. The found selectivities are consistent with computational simulations, particularly concerning the discrimination of differently planarized and charge-delocalized enolate tautomers by anion-π interactions. Enolate-π interactions on fullerenes are much shorter than standard π-π interactions and anion-π interactions on planar surfaces, and alternative cation-π interactions are not observed. These findings open new perspectives with regard to anion-π interactions in general and the use of carbon allotropes in catalysis.

Strain engineering of H/transition metal systems

The effects of both compressive and tensile surface strain on the hydrogenated low-index faces of Fe, Co, Ni, Cu, Ru, Rh, Pd, Ag, Os, Ir, Pt and Au has been investigated using density functional theory (DFT). Changes in the preferred H binding site have been observed on the Pd, Ir and Pt surfaces when the surface lattice constant is strained by up to 2%, and on Fe, Rh, Ag and Os surfaces for larger strains of up to 5%. A complete discussion of the variance of the hydrogen binding energy, charge, density of states and local geometry with strain is presented. The exchange-correlation, electrostatic and kinetic contributions to the binding energy are delineated and their respective contributions are discussed. The mechanism which determines the preferred binding site for each system under strain is shown to be complex.

Cooperative and anticooperative effects in resonance assisted hydrogen bonds in merged structures of malondialdehyde.

We analyzed non-additive effects in resonance assisted hydrogen bonds (RAHBs) in different β-enolones, which are archetypal compounds of these types of interactions. For this purpose, we used (i) potential energy curves to compute the formation energy, ΔE, of the RAHBs of interest in different circumstances along with (ii) tools offered by quantum chemical topology, namely, the Quantum Theory of Atoms In Molecules (QTAIM) and the Interacting Quantum Atoms (IQA) electronic energy partition. We established the effect that a given H-bond exerts over ΔE associated with another RAHB, determining in this way the cooperativity or the anticooperativity of these interactions. The mesomeric structures and the QTAIM delocalisation indices are consistent with the determined cooperative or anticooperative character of two given RAHBs. The HB cooperativity and anticooperativity studied herein are directly reflected in the IQA interaction energy E, but they are modulated by the surrounding hydrocarbon chain. The IQA decomposition of ΔEcoop, a measure of the cooperativity between a pair of interacting RAHBs, indicates that the analyzed H-bond cooperative/anticooperative effects are associated with greater/smaller (i) strengthening of the pseudo-bicyclic structure of the compounds of interest and (ii) electron localisations with its corresponding changes in the intra and intermolecular exchange-correlation contributions to ΔE. Overall, we expect that this investigation will provide valuable insights into the interplay among hydrogen bonded atoms and the π system in RAHBs contributing in this way to the understanding of the general features of H-bonds.

Equation of state of warm dense deuterium and its isotopes from density-functional theory molecular dynamics.

Of the two approaches of density-functional theory molecular dynamics, quantum molecular dynamics is limited at high temperature by computational cost whereas orbital-free molecular dynamics, based on an approximation of the kinetic electronic free energy, can be implemented in this domain. In the case of deuterium, it is shown how orbital-free molecular dynamics can be regarded as the limit of quantum molecular dynamics at high temperature for the calculation of the equation of state. To this end, accurate quantum molecular dynamics calculations are performed up to 20 eV at mass densities as low as 0.5g/cm^{3} and up to 10 eV at mass densities as low as 0.2g/cm^{3}. As a result, the limitation in temperature so far attributed to quantum molecular dynamics is overcome and an approach combining quantum and orbital-free molecular dynamics is used to construct an equation of state of deuterium. The thermodynamic domain addressed is that of the fluid phase above 1 eV and 0.2g/cm^{3}. Both pressure and internal energy are calculated as functions of temperature and mass density, and various exchange-correlation contributions are compared. The generalized gradient approximation of the exchange-correlation functional, corrected to approximately include the influence of temperature, is retained and the results obtained are compared to other approaches and to experimental shock data; in parts of the thermodynamic domain addressed, these results significantly differ from those obtained in other first-principles investigations which themselves disagree. The equations of state of hydrogen and tritium above 1 eV and above, respectively, 0.1g/cm^{3} and 0.3g/cm^{3}, can be simply obtained by mass density scaling from the results found for deuterium. This ab initio approach allows one to consistently cover a very large domain of temperature on the domain of mass density outlined above.

Exploring the Interaction Natures in Plutonyl (VI) Complexes with Topological Analyses of Electron Density.

The interaction natures between Pu and different ligands in several plutonyl (VI) complexes are investigated by performing topological analyses of electron density. The geometrical structures in both gaseous and aqueous phases are obtained with B3LYP functional, and are generally in agreement with available theoretical and experimental results when combined with all-electron segmented all-electron relativistic contracted (SARC) basis set. The Pu– bond orders show significant linear dependence on bond length and the charge of oxygen atoms in plutonyl moiety. The closed-shell interactions were identified for Pu-Ligand bonds in most complexes with quantum theory of atoms in molecules (QTAIM) analyses. Meanwhile, we found that some Pu–Ligand bonds, like Pu–OH−, show weak covalent. The interactive nature of Pu–ligand bonds were revealed based on the interaction quantum atom (IQA) energy decomposition approach, and our results indicate that all Pu–Ligand interactions is dominated by the electrostatic attraction interaction as expected. Meanwhile it is also important to note that the quantum mechanical exchange-correlation contributions can not be ignored. By means of the non-covalent interaction (NCI) approach it has been found that some weak and repulsion interactions existed in plutonyl(VI) complexes, which can not be distinguished by QTAIM, can be successfully identified.

Homogeneity properties of the embedding potential in frozen-density embedding theory

ABSTRACTIn numerical simulations, based on frozen-density embedding theory, the independent variables describing the total system are the embedded wave function (ΨA) and the density representing the environment. Due to inhomogeneity of the non-electrostatic component of the total energy: ), the expectation value of the embedding potential is not equal to the corresponding component of the total energy. The differences are evaluated using local and semi-local approximations for the functional EnadxcT[ρA, ρB] in two model systems representing embedded species weakly interacting with the environment. It is found that ΔnadxcT is typically one order of magnitude smaller than EnadxcT[ρA, ρB] and decreases with the overlap between and . The kinetic- and exchange-correlation contributions to ΔnadxcT cancel partially reducing its magnitude to mHartrees. Compared to local approximation for EnadxcT[ρA, ρB], the inhomogeneity is more pronounced in semi-local functionals.