Relative Energetics Research Articles

Creating Benchmarks for Lithium Clusters and Using Them for Testing and Validation.

Metal clusters often have a variety of possible structures, and they are calculated by a wide range of methods; however, fully converged benchmarks on the energy differences of structures and spin states that could be used to test or validate these methods are rare or nonexistent. Small lithium clusters are good candidates for such benchmarks to test different methods against well-converged relative energetics for qualitatively different structures because they have a small number of electrons. The present study provides fully converged benchmarks for Li4 and Li5 clusters and uses them to test a diverse group of approximation methods. To create a dataset of well-converged single-point energies for Li4 and Li5, stationary structures were optimized by Kohn-Sham density functional theory (KS-DFT) and then single-point energy calculations at these structures were carried out by two quite different beyond-CCSD(T) methods. To test other methods single-point energy calculations at these structures were carried out by KS-DFT, Mo̷ller-Plesset (MP) theory, coupled cluster (CC) theory, five composite methods (Gaussian-4, the Wuhan-Minnesota (WM) composite method, and the W2X, W3X, and W3X-L composite methods of Radom and co-workers), multiconfiguration pair-density functional theory (MC-PDFT), complete active space second-order perturbation theory (CASPT2), and n-electron valence state second-order perturbation theory (NEVPT2). Our results show that rhomboid and trigonal bipyramid (TBP) geometries are the most stable structures for Li4 and Li5, respectively. Using the W3X-L method to obtain our best estimates, the mean unsigned deviations were calculated for other methods for several structures and spin states of both Li4 and Li5. Binding energies and M diagnostics were calculated for all structures. The data in this paper are valuable for assessing the reliability of current electronic structure theories and also developing new density functionals and machine learned models.

Neural network potential for Zr-H

The introduction of Hydrogen (H) into Zirconium (Zr) influences many mechanical properties, especially due to low H solubility and easy formation of Zirconium hydride phases. Understanding the various effects of H requires studies with atomistic resolution but at scales that incorporate defects such as cracks, interfaces, and dislocations. Such studies thus demand accurate interatomic potentials. Here, a neural network potential (NNP) for the Zr-H system is developed within the Behler-Parrinello framework. The Zr-H NNP retains the accuracy of a recent NNP for hcp Zr and exhibits excellent agreement with first-principles density functional theory (DFT) for (i) H interstitials and their diffusion in hcp Zr, (ii) formation energies, elastic constants, and surface energies of relevant Zr hydrides, and (iii) energetics of a common Zr/Zr-H interface. The Zr-H NNP shows physical behavior for many different crack orientations in the most-stable ε-hydride, and structures and reasonable relative energetics for the 〈a〉 screw dislocation in pure Zr. This Zr-H NNP should thus be very powerful for future study of many phenomena driving H degradation in Zr that require atomistic detail at scales far above those accessible by first-principles.

Tailoring Interface Energies via Phosphonic Acids to Grow and Stabilize Cubic FAPbI3 Deposited by Thermal Evaporation.

Coevaporation of formamidinium lead iodide (FAPbI3) is a promising route for the fabrication of highly efficient and scalable optoelectronic devices, such as perovskite solar cells. However, it poses experimental challenges in achieving stoichiometric FAPbI3 films with a cubic structure (α-FAPbI3). In this work, we show that undesired hexagonal phases of both PbI2 and FAPbI3 form during thermal evaporation, including the well-known 2H-FAPbI3, which are detrimental for optoelectronic performance. We demonstrate the growth of α-FAPbI3 at room temperature via thermal evaporation by depositing phosphonic acids (PAc) on substrates and subsequently coevaporating PbI2 and formamidinium iodide. We use density-functional theory to develop a theoretical model to understand the relative growth energetics of the α and 2H phases of FAPbI3 for different molecular interactions. Experiments and theory show that the presence of PAc molecules stabilizes the formation of α-FAPbI3 in thin films when excess molecules are available to migrate during growth. This migration of molecules facilitates the continued presence of adsorbed organic precursors at the free surface throughout the evaporation, which lowers the growth energy of the α-FAPbI3 phase. Our theoretical analyses of PAc molecule-molecule interactions show that ligands can form hydrogen bonding to reduce the migration rate of the molecules through the deposited film, limiting the effects on the crystal structure stabilization. Our results also show that the phase stabilization with molecules that migrate is long-lasting and resistant to moist air. These findings enable reliable formation and processing of α-FAPbI3 films via vapor deposition.

Approaching the Complete Basis Set Limit for Spin-State Energetics of Mononuclear First-Row Transition Metal Complexes.

Convergence to the complete basis set (CBS) limit is analyzed for the problem of spin-state energetics in mononuclear first-row transition metal (TM) complexes by taking under scrutiny a benchmark set of 18 energy differences between spin states for 13 chemically diverse TM complexes. The performance of conventional CCSD(T) and explicitly correlated CCSD(T)-F12a/b calculations in approaching the CCSD(T)/CBS limits is systematically studied. An economic computational protocol is developed based on the CCSD-F12a approximation and (here proposed) modified scaling of the perturbative triples term (T#). This computational protocol recovers the relative spin-state energetics of the benchmark set in excellent agreement with the reference CCSD(T)/CBS limits (mean absolute deviation of 0.4, mean signed deviation of 0.2, and maximum deviation of 0.8 kcal/mol) and enables performing canonical CCSD(T) calculations for mononuclear TM complexes sized up to ca. 50 atoms, which is illustrated by application to heme-related metalloporphyrins. Furthermore, a good transferability of the basis set incompleteness error (BSIE) is demonstrated for spin-state energetics computed using CCSD(T) and other wave function methods (MP2, CASPT2, CASPT2/CC, NEVPT2, and MRCI + Q), which justifies efficient focal-point approximations and simplifies the construction of multimethod benchmark studies.

CuAAC-Based Synthesis, Copper-Catalyzed Aldehyde-Forming Hydrolytic Fission and Antiproliferative Evaluation of Novel Ferrocenoylamino-Substituted Triazole-Tethered Quinine-Chalcone Hybrids.

A series of novel triazole-tethered ferrocenoylamino-substituted cinchona-chalcone hybrids along with two representative benzoylamino-substituted reference compounds were prepared by three methods of CuAAC chemistry. In line with the limited success or complete failure of attempted conversions with low catalyst loadings, by means of DFT modeling studies, we demonstrated that a substantial part of the Cu(I) ions can be chelated and thus trapped in the aroylamino-substituted cinchona fragment and all of the accessible coordinating sites of the chalcone residues. Accordingly, increased amounts of catalysts were used to achieve acceptable yields; however, the cycloadditions with para-azidochalcones were accompanied by partial or complete aldehyde-forming hydrolytic fission of the enone C=C bond in a substituent-, solvent- and copper load-dependent manner. The experienced hydrolytic stability of the hybrids obtained by cycloadditions with ortho-azidochalcones was interpreted in terms of relative energetics, DFT reactivity indices and MO analysis of simplified models of two isomer copper-enone complexes. The novel hybrids were evaluated on HeLa, MDA-MB-231 and A2780 cell lines and showed substantial activity at low-to-submicromolar concentrations. An organometallic model carrying 3,4,5-trimethoxyphenyl residue in the enone part with a para-disubstituted benzene ring in the central skeletal region was identified as the most potent antiproliferative lead, characterized by submicromolar IC50 values measured on the three investigated cells. The biological assays also disclosed that this ferrocenoylamino-containing lead compound displays a ca. two- to five-fold more substantial antiproliferative effect than its benzoylamino-substituted counterpart.

Kiosk 8R-TA-11 - Relative Vortex Strength and Energetics Associate with Pulmonary Artery Pressure and Cardiac Function in Pulmonary Arterial Hypertension

Kiosk 8R-TA-11 - Relative Vortex Strength and Energetics Associate with Pulmonary Artery Pressure and Cardiac Function in Pulmonary Arterial Hypertension

Analytical and Parameter-Free Hückel Theory Made Possible for Symmetric Hx Clusters.

It is widely accepted that energetics of chemical bond breaking and formation can be described with simple mathematical forms only at the expense of extensive parameterization. In this work, the discovery of a simple tight-binding-type mathematical framework that can accurately predict the relative energetics of regular Hx polygons (2 ≤ x ≤ 15) in the ground states with their respective spin multiplicities using no parameters has been reported. The framework recasts Hückel theory in a density functional theory form by making use of Anderson and Adams-Gilbert theories of localized orbitals. For the systems examined, the method exhibits mean absolute errors of ∼0.02 Å (edge lengths) and ∼0.15 eV/atom (energy minima) relative to correlated-electron quantum chemistry calculations. Its accuracy is found to be comparable to the generalized gradient approximation and superior to standard parameterized tight binding and reactive potentials applied to Hx structures. Generalization of the theoretical framework to systems of many-electron atoms is presented, along with the comparison of the method to existing semiempirical tight binding and bond order potential approaches.

Characterizing the Semiconductor Nanocrystal Surface through Chemical Reactivity.

ConspectusMany desirable and undesirable properties of semiconductor nanocrystals (NCs) can be traced to the NC surface due to the large surface-to-volume ratio. Therefore, precise control of the NC surface is imperative to achieve NCs with the desired qualities. Ligand-specific reactivity and surface heterogeneity make it difficult to accurately control and tune the NC surface. Without a molecular-level appreciation of the NC surface chemistry, modulating the NC surface is impossible and the risk of introducing deleterious surface defects is imminent. To gain a more comprehensive understanding of the surface reactivity, we have utilized a variety of spectroscopic techniques and analytical methods in concert.This Account describes our use of robust characterization techniques and ligand exchange reactions in effort to establish a molecular-level understanding of NC surface reactivity. The utility of NCs in target applications such as catalysis and charge transfer hangs on precise tunability of NC ligands. Modulating the NC surface requires the necessary tools to monitor chemical reactions. One commonly utilized analytical method to achieve targeted surface compositions is 1H nuclear magnetic resonance (NMR) spectroscopy. Here we describe our use of 1H NMR spectroscopy to monitor chemical reactions at CdSe and PbS NC surfaces to identify ligand specific reactivity. However, seemingly straightforward ligand exchange reactions can vary widely depending on the NC materials and anchoring group. Some non-native X-type ligands will irreversibly displace native ligands. Other ligands exist in equilibrium with native ligands. Depending on the application, it is important to understand the nature of exchange reactions. This level of understanding can be obtained by extracting exchange ratios, exchange equilibrium, and reaction mechanism information from 1H NMR spectroscopy to establish precise NC reactivity.Reactivity that occurs through multiple, parallel ligand exchange mechanisms can involve both the liberation of metal-based Z-type ligands in addition to reactivity of X-type ligands. In these reactions, 1H NMR spectroscopy fails to discern between an X-type oleate or a Z-type Pb(oleate)2 because only the alkene resonance of the organic constituent is probed by this method. Multiple, parallel reaction pathways occur when thiol ligands are introduced to oleate-capped PbS NCs. This necessitated the use of synergistic characterization methods including 1H NMR spectroscopy, Fourier-transform infrared (FTIR) spectroscopy, and inductively coupled plasma mass spectrometry (ICP-MS) to characterize both surface-bound and liberated ligands.Similar analytical methods have been employed to probe the NC topology, which is an important, but often overlooked, component to NC reactivity given the facet-specific reactivity of PbS NCs. Through the tandem use of NMR spectroscopy and ICP-MS, we have monitored the liberation of Pb(oleate)2 as an L-type ligand is titrated to the NC to determine the quantity and equilibrium of Z-type ligands. By studying a variety of NC sizes, we correlated the number of liberated ligands with the size-dependent topology of PbS NCs.Lastly, we incorporate redox-active chemical probes into our toolbox to study NC surface defects. We describe how the site-specific reactivity and relative energetics of redox-active surface-based defects are elucidated using redox probes and show that this reactivity is highly dependent on the surface composition. This Account is designed to encourage readers to consider the necessary characterization techniques needed establish a molecular-level understanding of NC surfaces in their own work.

Influence of Dispersion Interactions on the Polymorphic Stability of Crystalline Oxides.

The accurate determination of relative phase stabilities using DFT methods is a significant challenge when some of these can vary by only a few kJ/mol. Here, we demonstrate that for a selection of oxides (TiO2, MnO2, and ZnO) the inclusion of dispersion interactions, accomplished using the DFT-D3 correction scheme, allows for the correct ordering and an improved calculation of the energy differences between polymorphic phases. The energetic correction provided is of the same order of magnitude as the energy difference between phases. D3-corrected hybrid functionals systematically yield results closest to experiment. We propose that the inclusion of dispersion interactions makes a significant contribution to the relative energetics of polymorphic phases, especially those with different densities, and should therefore be included for calculations of relative energies using DFT methods.

Anisotropic Phonon Bands in H-Bonded Molecular Crystals: The Instructive Case of α-Quinacridone

Phonons play a crucial role in the thermodynamic and transport properties of solid materials. Nevertheless, rather little is known about phonons in organic semiconductors. Thus, we employ highly reliable quantum mechanical calculations for studying the phonons in the α-polymorph of quinacridone. This material is particularly interesting, as it has highly anisotropic properties with distinctly different bonding types (H-bonding, π-stacking, and dispersion interactions) in different spatial directions. By calculating the overlaps of modes in molecular quinacridone and the α-polymorph, we associate Γ-point phonons with molecular vibrations to get a first impression of the impact of the crystalline environment. The situation becomes considerably more complex when analyzing phonons in the entire 1st Brillouin zone, where, due to the low symmetry of α-quinacridone, a multitude of avoided band crossings occur. At these, the character of the phonon modes typically switches, as can be inferred from mode participation ratios and mode longitudinalities. Notably, avoided crossings are observed not only as a function of the length but also as a function of the direction of the phonon wave vector. Analyzing these avoided crossings reveals how it is possible that the highest frequency acoustic band is always the one with the largest longitudinality, although longitudinal phonons in different crystalline directions are characterized by fundamentally different molecular displacements. The multiple avoided crossings also give rise to a particularly complex angular dependence of the group velocities, but combining the insights from the various studied quantities still allows drawing general conclusions, e.g., on the relative energetics of longitudinal vs transverse deformations (i.e., compressions and expansions vs slips of neighboring molecules). They also reveal how phonon transport in α-quinacridone is impacted by the reinforcing H-bonds and by π-stacking interactions (resulting from a complex superposition of van der Waals, charge penetration, and exchange repulsion).

Development and Benchmarking of Open Force Field 2.0.0: The Sage Small Molecule Force Field.

We introduce the Open Force Field (OpenFF) 2.0.0 small molecule force field for drug-like molecules, code-named Sage, which builds upon our previous iteration, Parsley. OpenFF force fields are based on direct chemical perception, which generalizes easily to highly diverse sets of chemistries based on substructure queries. Like the previous OpenFF iterations, the Sage generation of OpenFF force fields was validated in protein-ligand simulations to be compatible with AMBER biopolymer force fields. In this work, we detail the methodology used to develop this force field, as well as the innovations and improvements introduced since the release of Parsley 1.0.0. One particularly significant feature of Sage is a set of improved Lennard-Jones (LJ) parameters retrained against condensed phase mixture data, the first refit of LJ parameters in the OpenFF small molecule force field line. Sage also includes valence parameters refit to a larger database of quantum chemical calculations than previous versions, as well as improvements in how this fitting is performed. Force field benchmarks show improvements in general metrics of performance against quantum chemistry reference data such as root-mean-square deviations (RMSD) of optimized conformer geometries, torsion fingerprint deviations (TFD), and improved relative conformer energetics (ΔΔE). We present a variety of benchmarks for these metrics against our previous force fields as well as in some cases other small molecule force fields. Sage also demonstrates improved performance in estimating physical properties, including comparison against experimental data from various thermodynamic databases for small molecule properties such as ΔHmix, ρ(x), ΔGsolv, and ΔGtrans. Additionally, we benchmarked against protein-ligand binding free energies (ΔGbind), where Sage yields results statistically similar to previous force fields. All the data is made publicly available along with complete details on how to reproduce the training results at https://github.com/openforcefield/openff-sage.

Energetic vs. entropic stabilization between a Remdesivir analogue and cognate ATP upon binding and insertion into the active site of SARS-CoV-2 RNA dependent RNA polymerase.

SARS-CoV-2 RNA dependent RNA polymerase (RdRp) serves as a highly promising antiviral drug target such as for a Remdesivir nucleotide analogue (RDV-TP or RTP). In this work, we mainly used alchemical all-atom simulations to characterize relative binding free energetics between the nucleotide analogue RTP and natural cognate substrate ATP upon initial binding and pre-catalytic insertion into the active site of SARS-CoV-2 RdRp. Natural non-cognate substrate dATP and mismatched GTP were also examined for computation control. We first identified significant differences in dynamical responses between nucleotide initial binding and subsequent insertion configurations to the open and closed active sites of the RdRp, respectively, though the RdRp protein conformational changes between the active site's open and closed states are subtle. Our alchemical simulations indicated that upon initial binding (active site open), RTP and ATP show similar binding free energies to the active sites while in the insertion state (active site closed), ATP is more stabilized (∼-2.4 kcal mol-1) than RTP in free energetics. Additional analyses show, however, that RTP is more stabilized in binding energetics than ATP, in both the insertion and initial binding states, with RTP more stabilized due to the electrostatic energy in the insertion state and due to vdW energy in the initial binding state. Hence, it appears that natural cognate ATP still excels at association stability with the RdRp active site due to that ATP maintains sufficient flexibilities e.g., in base pairing with the template, which exemplifies an entropic contribution to the cognate substrate stabilization. These findings highlight the importance of substrate flexibilities in addition to energetic stabilization in antiviral nucleotide analogue design.

A data-driven interpretation of the stability of organic molecular crystals.

Due to the subtle balance of intermolecular interactions that govern structure-property relations, predicting the stability of crystal structures formed from molecular building blocks is a highly non-trivial scientific problem. A particularly active and fruitful approach involves classifying the different combinations of interacting chemical moieties, as understanding the relative energetics of different interactions enables the design of molecular crystals and fine-tuning of their stabilities. While this is usually performed based on the empirical observation of the most commonly encountered motifs in known crystal structures, we propose to apply a combination of supervised and unsupervised machine-learning techniques to automate the construction of an extensive library of molecular building blocks. We introduce a structural descriptor tailored to the prediction of the binding (lattice) energy and apply it to a curated dataset of organic crystals, exploiting its atom-centered nature to obtain a data-driven assessment of the contribution of different chemical groups to the lattice energy of the crystal. We then interpret this library using a low-dimensional representation of the structure-energy landscape and discuss selected examples of the insights into crystal engineering that can be extracted from this analysis, providing a complete database to guide the design of molecular materials.

Computational NMR investigation of mixed-metal (Al,Sc)-MIL-53 and its phase transitions.

Compositionally complex metal-organic frameworks (MOFs) have properties that depend on local structure that is often difficult to characterise. In this paper a density functional theory (DFT) computational study of mixed-metal (Al,Sc)-MIL-53, a flexible MOF with several different forms, was used to calculate the relative energetics of these forms and to predict NMR parameters that can be used to evaluate whether solid-state NMR spectroscopy can be used to differentiate, identify and characterise the forms adopted by mixed-metal MOFs of different composition. The NMR parameters can also be correlated with structural features in the different forms, giving fundamental insight into the nature and origin of the interactions that affect nuclear spins. Given the complexity of advanced NMR experiments required, and the potential need for expensive and difficult isotopic enrichment, the computational work is invaluable in predicting which experiments and approaches are likely to give the most information on the disorder, local structure and pore forms of these mixed-metal MOFs.

Unambiguous spectral characterization on triplet energy transfer from quantum dots mediated by hole transfer competing with other carrier dynamics.

Triplet generation by quantum dots (QDs)-sensitized molecules emerges great potential in many applications. However, the mechanism of triplet energy transfer (TET) is still fuzzy especially due to the complicated energy level alignment of QDs and molecules or trap states in QDs. Here, CdSe QDs and 5-tetracene carboxylic acid (TCA) molecules are selected as the triplet donor and acceptor, respectively, to form a TET system. By tuning the band gap of CdSe, the CdSe-TCA complex is exactly designed to present a Type-II like alignment of relative energetics. Coupling the transient absorption and time-resolved fluorescence spectra, all carrier dynamics is distinctly elucidated. Quantitative analysis demonstrates that hole transfer persisting for ∼ 2 ps outcompetes all other carrier dynamics such as electron trapping (∼100 ps level), charge recombination (∼ 5 ns) and the so-called "back transfer charge recombination" (∼50 ns), and thus leads to a hole-transfer-mediated TET process. The low TET yield (∼34.0%) ascribed to electron behavior can be further improved if electron trapping and charge recombination are efficiently suppressed. The observation on distinguishable carrier dynamics attributed to legitimate design of energy level alignment facilitates a better understanding of the TET mechanism from QDs to molecules as well as further development of photoelectronic devices based on such TET systems.

DFT Investigation of the catalytic conversion of acetic acid to acetone on the zeolite H-ZSM5

Biomaterials processing has become increasingly important in the chemical industry. Identification and optimization of processes to facilitate biomass conversion is therefore recognized as being of great importance. Carboxylic acids derived from biomaterials are important building blocks that can be used in a wide variety of industrial applications. As such, methods to process them in an efficient and cost-effective manner are highly desirable. In this study we report the use of theoretical methods to explore the catalytic conversion of acetic acid to acetone on the zeolite H-ZSM5. We have employed a 46T DFT cluster model to explore mechanistic proposals reported in the literature. We investigate the relative energetics associated with the formation of the proposed intermediates, including acyl-zeolite complexes, enols, acylium cations, ketenes, anhydrides, and beta-keto acids, that could potentially connect acetic acid to the desired product, acetone. This assessment would allow us to identify the most probable mechanism connecting the reactant to products. We predict a low energy pathway starting with the generation of a surface acyl, followed by an anhydride, with the rate determining step involving methyl group migration. The reaction is predicted to be bi-molecular and involves C-C bond formation, in line with proposals based on isotopic labelling experiments.

Probing fragmentation mechanisms of deprotonated isomaltotriose: Charge-remote or charge-directed?

Unimolecular dissociation mechanisms of polysaccharide ions are typically separated into three groups: cross-ring cleavages (forming A/X ions), and those resulting from cleaving different sides of the glycosidic bond (forming B/Y ions) and (C/Z ions). In this study, the relative energetics of the first steps of the cross-ring cleavage and both glycosidic bond cleavage channels for isomaltotriose, [glc(α1-6)glc(α1-6)glc], as well as a minor water loss channel were explored using density functional theory (DFT) calculations at the B3LYP/6-31+g(d) level of theory. RRKM approximations of the competing rate constants indicate that charge-remote reactions have energy barriers and activation entropies that deny them competition with charge-directed reactions in this example.

Sulfite-Catalyzed Nucleophilic Substitution Reactions with Thiamin and Analogous Pyrimidine Donors Proceed via an SNAE Mechanism.

When treated with SO32-, thiamin undergoes a substitution reaction to release a thiazole leaving group and the corresponding sulfonate. Although this reaction could proceed via a simple SN2-like mechanism, a multistep addition-elimination (SNAE) mechanism involving the addition of SO32- to C6' of the 4-aminopyrimidine of thiamin has also been proposed. Although this reaction has potential utility in the synthesis of substituted pyrimidines and provides a direct analogue to reactions catalyzed by thiaminases, a detailed mechanistic picture of the SO32--catalyzed cleavage of thiamin has remained elusive. Here, DFT calculations have been used to probe the relative energetics and the factors that shape the potential energy surfaces that define the possible mechanisms of substitution. These calculations provide clear support for the SNAE mechanism over an SN2-like process and illustrate that the unique ability of SO32- to activate thiamin toward nucleophilic displacement is due to the combined nucleophilicity and relatively poor leaving group ability of SO32-. Both of these factors favor the forward partitioning of the sulfite adduct toward the cleavage products whereas adducts formed with other nucleophiles overwhelmingly revert to reactants. Calculations performed with a range of substrates with various electrophilicities and nucleofugalities consistently suggest that the SNAE pathway is significantly lower in energy than the direct substitution, illustrating that this SO32--catalyzed multistep process is likely to be broadly applicable both in solution and in catalysis by thiaminases.

Understanding the Reactivity of Combination Reactions of Intramolecular Geminal Group 13 Element/Phosphorus and Gallium/Group 15 Element Frustrated Lewis Pairs with CS2.

The reactions of CS2 captured by intramolecular geminal G13/P-based (G13 = group 13 elements) and Ga/G15-based (G15 = group 15 elements) frustrated Lewis pairs have been theoretically examined by using density functional theory (DFT) computations. With regard to the nine FLP-related compounds, our DFT calculated results reveal that only Al/P-Rea and Ga/P-Rea can kinetically and thermodynamically precede the energetically feasible combination reactions with CS2 to form the five-membered heterocyclic adducts. Our activation strain model analyses on the nine aforementioned model molecules indicate that the atomic radius of the Lewis acceptor (G13) and the Lewis donor (G15) plays a role in controlling their barrier heights to obtain good orbital overlaps among G13/P-Rea, Ga/G15-Rea, and CS2. Our theoretical observations based on the energy decomposition analysis-natural orbitals for chemical valence (EDA-NOCV) approach strongly indicate that the donor-acceptor bonding (i.e., singlet-singlet bonding) rather than the electron-sharing bonding (i.e., triplet-triplet bonding) plays a central role in determining the bonding conditions of the transition states, G13/P-TS and Ga/G15-TS. In addition, the theoretical evidence obtained by the frontier molecular orbital theory and EDA-NOCV analyses reveals that the best description for the bonding natures of the combination reactions of intramolecular geminal G13/P-Rea and Ga/G15-Rea with CS2 is the lone pair(G15) → p-π*(C) interaction rather than the p-π*(G13) ← p-π(S) interaction. Moreover, our present DFT computations concerning the calculated structures and corresponding relative energetics of the stationary points connected with the aforementioned sophisticated approaches are in accordance with the Hammond postulate.

Insights into the photoinduced degradation of terbuthylazine from aqueous solution: The synergic effects generated from hydrogen-bond interactions

In this work, we examined the stability of terbuthylazine (TBA) as well as the mechanisms related to its photoinduced degradation considering two conditions: i) when exposed (solely) to UV radiation and ii) when exposed to UV radiation in presence of H2O2. Ground state properties (such as bond lengths, hydrogen-bond interactions, and relative energetics) as well as excited state parameters (vertical excitation energies and generalized oscillator strengths, GOS) were determined through the use of standard and well documented density functional theory based electronic structure methods and, subsequently, explored for accomplishing the purpose of the work; results were obtained considering both the gas-phase and a (water) solvent environment. In terms of ground state, the presence of the additional hydrogen-bond interaction when the H2O2 molecule is nearing the N atom located in the opposite side to the Cl atom in TBA suggested the system to be more stable (for instance, 9.93 kJ/mol at the CAM-B3LYP/aug-cc-pVTZ level of theory in water) than its counterpart with the H2O2 molecule neighboring the N and Cl atoms. The excitation energies determined considering solvation were found to be in fairly good agreement when compared to the corresponding obtained in the gas-phase. Interestingly and more important, a different situation is observed in terms of GOS when solvation is considered. In this case, two excited states (versus no state in the gas-phase) were found to have considerably large GOS and, thus, are expected to be accessible at the UV region. Therefore, the interactions between TBA, H2O2, and water molecules (from the solvent environment) play an important role in regards to the excited states to be accessed. By coexisting with the UV photolysis of the H2O2 molecules in aqueous solution, the synergic effects from the TBA-H2O2-H2O interactions may add to the enhancement in the yield of the photoinduced degradation of TBA.