Perspectives on Basis Sets Beautiful: Seasonal Plantings of Diffuse Basis Functions

The dioxasilepanyl group as a versatile organometallic unit: studies on stability, reactivity, and utility

Organic synthesis is performed based on precise choices of functional groups and reactions employed. In a multistep synthesis, an ideal functional group should be compatible with various reaction conditions and unaltered until it is subjected to a selective conversion. The current study was set out to search for a silicon functionality that meets these criteria. Here we have established a new silicon-based synthetic methodology centred on a bulky 7-membered dialkoxysilyl group (2,4,4,7,7-pentamethyl-1,3,2-dioxasilepan-2-yl) that uniquely has both stability and on-demand reactivity. The exceptional stability of this functional group was corroborated by both experimental and computational studies which demonstrated that key factors for its stability were a 7-membered structure and steric hindrance. In turn, the dioxasilepanyl group was found to become reactive and to be easily transformed in the presence of appropriate activators. Combined with the development of easy and robust methods to introduce the dioxasilepanyl group onto aryl rings, these findings have allowed a shorter and more efficient synthesis of a bioactive molecule, thus demonstrating the potential utility of the easily accessible dioxasilepanyl group in organic synthesis.

Alkylammonium Cation Affinities of Nitrogenated Organobases: The Roles of Hydrogen Bonding and Proton Transfer

Alkylammonium cation affinities of 64 nitrogen-containing organobases, as well as the respective proton transfer processes from the alkylammonium cations to the base, have been computed in the gas phase by using DFT methods. The guanidine bases show the highest proton transfer values (191.9-233 kJ mol^-1 ) whereas the cis-2,2'-biimidazole presents the largest affinity towards the alkylammonium cations (>200 kJ mol^-1 ) values. The resulting data have been compared with the experimentally reported proton affinities of the studied nitrogen-containing organobases revealing that the propensity of an organobase for the proton transfer process increases linearly with its proton affinity. This work can provide a tool for designing senors for bioactive compounds containing amino groups that are protonated at physiological pH.

General Perturb-Then-Diagonalize Model for the Vibrational Frequencies and Intensities of Molecules Belonging to Abelian and Non-Abelian Symmetry Groups

In this paper, we show that the standard second-order vibrational perturbation theory (VPT2) for Abelian groups can be used also for non-Abelian groups without employing specific equations for two- or threefold degenerate vibrations but rather handling in the proper way all the degeneracy issues and deriving the peculiar spectroscopic signatures of non-Abelian groups (e.g., l -doubling) by a posteriori transformations of the eigenfunctions. Comparison with the results of previous conventional implementations shows a perfect agreement for the vibrational energies of linear and symmetric tops, thus paving the route to the transparent extension of the equations already available for asymmetric tops to the energies of spherical tops and the infrared and Raman intensities of molecules belonging to non-Abelian symmetry groups. The whole procedure has been implemented in our general engine for vibro-rotational computations beyond the rigid rotor/harmonic oscillator model and has been validated on a number of test cases.

Toward an Accurate Ab Initio Description of Low-Lying Singlet Excited States of Polyenes

The low-lying excited states of carotenoids play a crucial role in many important biophysical processes such as photosynthesis. Most of these excited states are strongly correlated, which makes them both challenging for a qualitative ab initio description and an engaging model system for trying out emerging multireference methods. Among these methods, driven similarity renormalization group (DSRG) and its perturbative version (DSRG-MRPT2) are especially attractive in terms of both accuracy and moderate numerical complexity. In this paper, we applied density matrix renormalization group (DMRG) followed by DSRG-MRPT2 for the calculation of vertical and adiabatic excitation energies into the 2A_g^-, 1B_u^-, and 1B_u⁺ electronic states of polyenes containing from 8 to 13 conjugating double bonds acting as a model for natural carotenoids. It was shown that the DSRG flow parameter should be adjusted to ensure both the energy convergence with respect to it and the agreement with the experimental data. With the increased flow parameter, the proposed combination of methods provides a reasonable agreement with the experiment. The deviations of the adiabatic excitation energies are less than 1000 cm^-1 for the 2A_g^- and less than 3000 cm^-1 for the excited states of the B_u symmetry, which in terms of accuracy significantly outperforms the N-electron valence state perturbation theory. At the same time, DSRG-MRPT2 is shown to be robust with respect to variation of quality of the DMRG reference wave function such as the orbital optimization or the number of electronic states in the averaging.

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low-energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. Intensity borrowing by these spurious states can affect intensities of the valence excitations, altering electronic bandshapes. At higher excitation energies, it is difficult to distinguish spurious charge-transfer states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces "natural charge-transfer orbitals" that provide a means to isolate orbitals that are most likely to participate in a CTTS excitation. Projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation and be used as an automated subspace selection algorithm for projection-based embedding of a high-level description of the CTTS state in a lower-level description of its environment. We apply this method to an ab initio molecular dynamics trajectory of I^-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with the experiment, and only one-third of the water molecules in the I^-(H₂O)₉₆ simulation cell need to be described with LR-TDDFT to obtain excitation energies that are converged to <0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments.

SERS Investigation on Oligopeptides Used as Biomimetic Coatings for Medical Devices

The surface-enhanced Raman scattering (SERS) spectra of three amphiphilic oligopeptides derived from EAK16 (AEAEAKAK)2 were examined to study systematic amino acid substitution effects on the corresponding interaction with Ag colloidal nanoparticles. Such self-assembling molecular systems, known as "molecular Lego", are of particular interest for their uses in tissue engineering and as biomimetic coatings for medical devices because they can form insoluble macroscopic membranes under physiological conditions. Spectra were collected for both native and gamma-irradiated samples. Quantum mechanical data on two of the examined oligopeptides were also obtained to clarify the assignment of the prominent significative bands observed in the spectra. In general, the peptide-nanoparticles interaction occurs through the COO- groups, with the amide bond and the aliphatic chain close to the colloid surface. After gamma irradiation, mimicking a free oxidative radical attack, the SERS spectra of the biomaterials show that COO- groups still provide the main peptide-nanoparticle interactions. However, the spatial arrangement of the peptides is different, exhibiting a systematic decrease in the distance between aliphatic chains and colloid nanoparticles.

Gaseous cyclodextrin-closo-dodecaborate complexes χCD·B12X122- (χ = α, β, and γ; X = F, Cl, Br, and I): electronic structures and intramolecular interactions

A fundamental understanding of cyclodextrin-closo-dodecaborate inclusion complexes is of great interest in supramolecular chemistry. Herein, we report a systematic investigation on the electronic structures and intramolecular interactions of perhalogenated closo-dodecaborate dianions B12X122- (X = F, Cl, Br and I) binding to α-, β-, and γ-cyclodextrins (CDs) in the gas phase using combined negative ion photoelectron spectroscopy (NIPES) and density functional theory (DFT) calculations. The vertical detachment energy (VDE) of each complex and electronic stabilization of each dianion due to the CD binding (ΔVDE, relative to the corresponding isolated B12X122-) are determined from the experiments along α-, β- and γ-CD in the form of VDE (ΔVDE): 4.00 (2.10), 4.33 (2.43), and 4.30 (2.40) eV in X = F; 4.09 (1.14), 4.64 (1.69), and 4.69 (1.74) eV in X = Cl; 4.11 (0.91), 4.58 (1.38), and 4.70 (1.50) eV in X = Br; and 3.54 (0.74), 3.88 (1.08), and 4.05 (1.25) eV in X = I, respectively. All complexes have significantly higher VDEs than the corresponding isolated dodecaborate dianions with ΔVDE spanning from 0.74 eV at (α, I) to 2.43 eV at (β, F), sensitive to both host CD size and guest substituent X. DFT-optimized complex structures indicate that all B12X122- prefer binding to the wide openings of CDs with the insertion depth and binding motif strongly dependent on the CD size and halogen X. Dodecaborate anions with heavy halogens, i.e., X = Cl, Br, and I, are found outside of α-CD, while B12F122- is completely wrapped by γ-CD. Partial embedment of B12X122- into CDs is observed for the other complexes via multipronged B-XH-O/C interlocking patterns. The simulated spectra based on the density of states agree well with those of the experiments and the calculated VDEs well reproduce the experimental trends. Molecular orbital analyses suggest that the spectral features at low binding energies originated from electrons detached from the dodecaborate dianion, while those at higher binding energies are derived from electron detachment from CDs. Energy decomposition analyses reveal that the electrostatic interaction plays a dominating role in contributing to the host-guest interactions for the X = F series partially due to the formation of a O/C-HX-B hydrogen bonding network, and the dispersion forces gradually become important with the increase of halogen size.

Optimized damping parameters for empirical dispersion corrections to symmetry-adapted perturbation theory

Symmetry-adapted perturbation theory (SAPT) has become an invaluable tool for studying the fundamental nature of non-covalent interactions by directly computing the electrostatics, exchange (steric) repulsion, induction (polarization), and London dispersion contributions to the interaction energy using quantum mechanics. Further application of SAPT is primarily limited by its computational expense, where even its most affordable variant (SAPT0) scales as the fifth power of system size [O(N⁵)] due to the dispersion terms. The algorithmic scaling of SAPT0 is reduced from O(N⁵)→O(N⁴) by replacing these terms with the empirical D3 dispersion correction of Grimme and co-workers, forming a method that may be termed SAPT0-D3. Here, we optimize the damping parameters for the -D3 terms in SAPT0-D3 using a much larger training set than has previously been considered, namely, 8299 interaction energies computed at the complete-basis-set limit of coupled cluster through perturbative triples [CCSD(T)/CBS]. Perhaps surprisingly, with only three fitted parameters, SAPT0-D3 improves on the accuracy of SAPT0, reducing mean absolute errors from 0.61 to 0.49 kcal mol^-1 over the full set of complexes. Additionally, SAPT0-D3 exhibits a nearly 2.5× speedup over conventional SAPT0 for systems with ∼300 atoms and is applied here to systems with up to 459 atoms. Finally, we have also implemented a functional group partitioning of the approach (F-SAPT0-D3) and applied it to determine important contacts in the binding of salbutamol to G-protein coupled β₁-adrenergic receptor in both active and inactive forms. SAPT0-D3 capabilities have been added to the open-source Psi4 software.

Mechanism-Based Insights into Removing the Mutagenicity of Aromatic Amines by Small Structural Alterations

Aromatic and heteroaromatic amines (ArNH₂) are activated by cytochrome P450 monooxygenases, primarily CYP1A2, into reactive N-arylhydroxylamines that can lead to covalent adducts with DNA nucleobases. Hereby, we give hands-on mechanism-based guidelines to design mutagenicity-free ArNH₂. The mechanism of N-hydroxylation of ArNH₂ by CYP1A2 is investigated by density functional theory (DFT) calculations. Two putative pathways are considered, the radicaloid route that goes via the classical ferryl-oxo oxidant and an alternative anionic pathway through Fenton-like oxidation by ferriheme-bound H₂O₂. Results suggest that bioactivation of ArNH₂ follows the anionic pathway. We demonstrate that H-bonding and/or geometric fit of ArNH₂ to CYP1A2 as well as feasibility of both proton abstraction by the ferriheme-peroxo base and heterolytic cleavage of arylhydroxylamines render molecules mutagenic. Mutagenicity of ArNH₂ can be removed by structural alterations that disrupt geometric and/or electrostatic fit to CYP1A2, decrease the acidity of the NH₂ group, destabilize arylnitrenium ions, or disrupt their pre-covalent transition states with guanine.

Competition of Intra- and Intermolecular Forces in Anthraquinone and Its Selected Derivatives

Intra- and intermolecular forces competition was investigated in the 9,10-anthraquinone (1) and its derivatives both in vacuo and in the crystalline phase. The 1,8-dihydroxy-9,10-anthraquinone (2) and 1,8-dinitro-4,5-dihydroxy-anthraquinone (3) contain Resonance-Assisted Hydrogen Bonds (RAHBs). The intramolecular hydrogen bonds properties were studied in the electronic ground and excited states employing Møller-Plesset second-order perturbation theory (MP2), Density Functional Theory (DFT) method in its classical formulation as well as its time-dependent extension (TD-DFT). The proton potential functions were obtained via scanning the OH distance and the dihedral angle related to the OH group rotation. The topological analysis was carried out on the basis of theories of Atoms in Molecules (AIM-molecular topology, properties of critical points, AIM charges) and Electron Localization Function (ELF-2D maps showing bonding patterns, calculation of electron populations in the hydrogen bonds). The Symmetry-Adapted Perturbation Theory (SAPT) was applied for the energy decomposition in the dimers. Finally, Car-Parrinello molecular dynamics (CPMD) simulations were performed to shed light onto bridge protons dynamics upon environmental influence. The vibrational features of the OH stretching were revealed using Fourier transformation of the autocorrelation function of atomic velocity. It was found that the presence of OH and NO2 substituents influenced the geometric and electronic structure of the anthraquinone moiety. The AIM and ELF analyses showed that the quantitative differences between hydrogen bonds properties could be neglected. The bridged protons are localized on the donor side in the electronic ground state, but the Excited-State Intramolecular Proton Transfer (ESIPT) was noticed as a result of the TD-DFT calculations. The hierarchy of interactions determined by SAPT method indicated that weak hydrogen bonds play modifying role in the organization of these crystal structures, but primary ordering factor is dispersion. The CPMD crystalline phase results indicated bridged proton-sharing in the compound 2.

Interplay of stereo-electronic, vibronic and environmental effects in tuning the chiroptical properties of an Ir(III) cyclometalated N-heterocyclic carbene

Chiroptical spectra are among the most suitable techniques for investigating the ground and excited electronic states of chiral systems, but their interpretation is not straightforward and strongly benefits from quantum chemical simulations, provided that the employed computational model is sufficiently accurate and deals properly with stereo-electronic, vibrational averaging and environmental effects. Since the synergy among all these effects is only rarely accounted for, especially for large and flexible organometallic systems, the main aim of this contribution is to illustrate the latest developments of computational approaches rooted into the density functional theory for describing stereo-electronic effects and complemented by effective techniques to deal with vibrational modulation effects and solvatochromic shifts. In this connection, chiral iridium complexes offer an especially suitable case study in view of their bright phosphorescence, which is particularly significant for building effective light emitting diodes (OLEDs) and biomarkers and can be finely tuned by the nature of the metal ligands. For instance, a recently synthesized family of cycloiridiated complexes, KC and KD, bearing a pentahelicenic N-heterocyclic carbene (KB), has shown an enhanced long-lasting, bright phosphorescence. Deeper insights into the still unclear nature and origin of the enhancement could be gained by the interpretation of the chiroptical spectra, which is quite challenging in view of the presence of two sources of chirality, the chiral center on Ir and the chiral axis related to the helicene ligand, in addition to the relativistic effects related to the presence of the Ir center. At the same time, the large dimensions of KC and KD hamper the use of the most sophisticated (but prohibitively expensive) computational models, so that more approximate approaches must be validated on a suitable model compound. To this end, after optimizing the computational scheme on a model system devoid of the helicene moiety (KA), we have performed a comprehensive investigation of the KC and KD spectra, whose interpretation is further aided by novel graphical tools. The discussion and analysis of the results will not be focused on the theoretical background, but, rather, on practical details (specific functional, basis set, vibronic model, solvent regime) with the aim of providing general guidelines for the use of last-generation computational spectroscopy tools also by non-specialists.

4-Fluoro-Threonine: From Diastereoselective Synthesis to pH-Dependent Conformational Equilibrium in Aqueous Solution

4-Fluoro-threonine, the only fluoro amino acid of natural origin discovered so far, is an interesting target for both synthetic and theoretical investigations. In this work, we lay the foundation for spectroscopic characterization of 4-fluoro-threonine. First, we report a diastereoselective synthetic route, which is suitable to produce synthetic material for experimental characterization. The addition of the commercially available ethyl isocyanoacetate to benzyloxyacetaldehyde led to the corresponding benzyloxy-oxazoline, which was hydrolyzed and transformed into ethyl (4S*,5S*)-5-hydroxymethyl-2-oxo-4-oxazolidinecarboxylate in a few steps. Fluorination with diethylamino sulfur trifluoride (DAST) afforded ethyl (4S*,5S*)-5-fluoromethyl-2-oxo-4-oxazolidinecarboxylate, which was deprotected to give the desired diastereomerically pure 4-fluoro-threonine, in 8-10% overall yield. With the synthetic material in our hands, acid-base titrations have been carried out to determine acid dissociation constants and the isoelectric point, which is the testing ground for the theoretical analysis. We have used machine learning coupled with quantum chemistry at the state-of-the-art to analyze the conformational space of 4-fluoro-threonine, with the aim of gaining insights from the comparison of computational and experimental results. Indeed, we have demonstrated that our approach, which couples a last-generation double-hybrid density functional including empirical dispersion contributions with a model combining explicit first-shell molecules and a polarizable continuum for describing solvent effects, provides results and trends in remarkable agreement with experiments. Finally, the conformational analysis applied to fluoro amino acids represents an interesting study for the effect of fluorine on the stability and population of conformers.

A computational study of the reaction mechanism involved in the fast cleavage of an unconstrained amide bond assisted by an amine intramolecular nucleophilic attack

In the present work, the fast amide bond cleavage of [3-((1R,5S,7s)-3-azabicyclo[3.3.1]nonane-7-carbonyl)-3-azabicyclo[3.3.1]nonane-7-carboxylic acid (bi-ATDO)], through an intramolecular nucleophilic attack of an amine group is evaluated. First, six possible peptide bond cleavage mechanisms, two of them including a water molecule, are described at the ωB97XD/6-311 + G(d,p)//MP2/6-311 + G(d,p) level of theory. The reaction consisting of an intramolecular nitrogen nucleophilic attack followed by a proton transfer and the amide bond cleavage is determined as the most favorable mechanism. The activation free energy computed for the latter is 20.5 kcal mol^-1 , which agrees with the reported experimental result of 24.8 kcal mol^-1 . Inclusion of a water molecule to assist the first step of the reaction results in an activation free energy increase of about 17 kcal mol^-1 . All the steps in the most favorable mechanism are studied more in detail employing intrinsic reaction coordinate as well as the reaction force and reaction electronic flux analysis.

Multiconfiguration Pair-Density Functional Theory

Kohn-Sham density functional theory with the available exchange-correlation functionals is less accurate for strongly correlated systems, which require a multiconfigurational description as a zero-order function, than for weakly correlated systems, and available functionals of the spin densities do not accurately predict energies for many strongly correlated systems when one uses multiconfigurational wave functions with spin symmetry. Furthermore, adding a correlation functional to a multiconfigurational reference energy can lead to double counting of electron correlation. Multiconfiguration pair-density functional theory (MC-PDFT) overcomes both obstacles, the second by calculating the quantum mechanical part of the electronic energy entirely by a functional, and the first by using a functional of the total density and the on-top pair density rather than the spin densities. This allows one to calculate the energy of strongly correlated systems efficiently with a pair-density functional and a suitable multiconfigurational reference function. This article reviews MC-PDFT and related background information.

Mechanistic insights into paracetamol transformation in UV/NH2Cl process: Experimental and theoretical study

The UV/monochloramine (NH₂Cl) process is an advanced oxidation process that can effectively remove emerging contaminants (ECs). However, the degradation mechanisms of reactive radicals with ECs are not clear. In this work, we combined theoretical calculations with experimental studies to investigate the kinetics and mechanism of radical-mediated degradation of paracetamol (AAP) in UV/NH₂Cl process. The degradation of AAP in UV/NH₂Cl process accords with the pseudo first-order kinetics. Impact factors including NH₂Cl dose, pH, natural organic matter, HCO₃^-, and NO₃^- were evaluated. The reaction mechanisms of AAP with hydroxyl radical (HO^·), reactive chlorine species (RCS), and reactive nitrogen species (RNS) were discussed in detail. Specifically, HO^· attacked AAP mainly through hydrogen atom transfer (HAT) and radical adduct formation (RAF), while Cl₂^·- play a certain role through single electron transfer (SET). ^·NH₂ and Cl^· destructed AAP mainly through HAT. Based on the mechanism analysis, the second-order rate constants of AAP reacts with HO^·, Cl^·, ^·NH₂, ClO^·, Cl₂^·- and ^·NO₂ were calculated through transition state theory as 2.66×10⁹ M^-1 s^-1, 2.61×10⁹ M^-1 s^-1, 1.02×10⁷ M^-1 s^-1, 7.74×10⁶ M^-1 s^-1, 1.32×10⁶ M^-1 s^-1, 1.48×10³ M^-1 s^-1 respectively. The second-order rate constants were then used to distinguish the contribution of radicals to the degradation of AAP. Thirteen transformation products were identified by high-resolution mass spectrometry. Combined active sites with potential energy surface, the detailed reaction pathways were proposed. Overall, this study provides deep insights into the mechanism of radical-mediated degradation of AAP.

Tuning the first hyperpolarizability of hexaphyrins with different connections of mislinked pyrrole units: a theoretical study

In the satisfactory design and synthesis of high-performance nonlinear optical (NLO) materials, for meeting the rapidly expanding demands of optoelectronic devices, a deeper understanding of the relationship between the structures and NLO properties has become a key issue. Herein, five novel mislinked expanded hexaphyrins with different connections of pyrrole units are selected to study the relationship between the structures and NLO properties. These five mislinked expanded hexaphyrins are neo-fused, neo-confused hexaphyrins, singly, doubly, and triply N-confused hexaphyrins. From theoretical results, the order of the static first hyperpolarizability (β₀) values is found to be: neo-fused hexaphyrin (β₀ = 4163 a.u.) < neo-confused hexaphyrin (β₀ = 5494 a.u.) < singly N-confused hexaphyrin (β₀ = 6510 a.u.) < doubly N-confused hexaphyrin (β₀ = 15 130 a.u.) < triply N-confused hexaphyrin (β₀ = 26 095 a.u.). Furthermore, β₀ values of the doubly and triply N-confused hexaphyrins are improved 2.1 and 3.7 times over that of their usual parent hexaphyrin (β₀ = 7120 a.u.), respectively. It is worth noting that increasing mislinked connection numbers and changing mislinked connection ways of the pyrrole units in these mislinked expanded hexaphyrins plays a crucial role in the tune of their second-order NLO responses.

Predicting OH stretching fundamental wavenumbers of alcohols for conformational assignment: different correction patterns for density functional and wave-function-based methods

A model is presented for the prediction of OH stretching fundamental wavenumbers of alcohol conformers in the gas phase by application of a small set of empirical anharmonicity corrections to calculations in the harmonic approximation. In contrast to the popular application of a uniform scaling factor, the local chemical structure of the alcohol is taken into account to greatly improve accuracy. Interestingly, different correction patterns emerge for results of hybrid density functional (B3LYP-D3 and PBE0-D3) and wave-function-based methods (SCS-LMP2, LCCSD(T*)-F12a and CCSD(T)-F12a 1D). This raises questions about electronic structure deficiencies in these methods and differences in anharmonicity between alcohols. After its initial construction on the basis of literature assignments the model is tested with Raman jet spectroscopy of propargyl alcohol, cyclohexanol, borneol, isopinocampheol and 2-methylbutan-2-ol. For propargyl alcohol a spectral splitting attributed to tunneling is resolved. PBE0-D3 is identified as a well performing and broadly affordable electronic structure method for this model. A mean absolute error of 1.3 cm-1 and a maximum absolute error of 3 cm-1 result for 46 conformers of 24 alcohols in a 60 cm-1 range, when a single parameter is adjusted separately for each alcohol substitution class (methanol, primary, secondary, tertiary).

From luminescence quenching to high-efficiency phosphorescence: a theoretical study on the monomeric and dimeric forms of platinum(II) complexes with both 2-pyridylimidazol-2-ylidene and bipyrazolate chelates

To develop solid-state light-emitting materials with high luminescence efficiency, determining the potential photophysics and luminescence mechanisms of the aggregation state remains a challenge and a priority. Here, we apply density functional theory to study the photophysical properties of a series of square planar Pt(ii) complexes in both monomeric and dimeric forms. We reveal that four monomeric Pt(ii) complexes are dominated by triplet ligand-to-ligand charge-transfer, and the lack of the triplet metal-to-ligand charge-transfer feature results in weak spin-orbit coupling (SOC), which leads to limited radiative rates; moreover, calculated nonradiative transition rates are one or two orders of magnitude higher than those radiative rates because a large amount of reorganization energy caused by the vibration of the bipyrazolate (bipz) ligand cannot be readily suppressed in the monomeric form. Therefore, four monomers exhibit photoluminescence quenching in CH2Cl2 solution in both theoretical calculations and experiments. However, in the solid state, the intense luminescence phenomenon indicates obviously distinct properties between the monomer and aggregation. We carried out a dimer model to interpret that the interaction of PtPt induces a metal-metal-to-ligand charge-transfer excimeric state, which leads more metal components to participate in the charge transfer and enhance the SOC effect. At the same time, the ligand vibration can be significantly reduced by the shortened distance, and there is a strong π-π packing interaction in the dimer; thus, an excellent quantum yield can be achieved in aggregation. In addition, we disclose that introducing bulky substituents bearing electron-donating groups at R' and R'' positions have little effect on the properties of the monomers; however, there is a benefit of restricting the internal reorganization energy through the intermolecular interaction when packing in the solid state. Therefore, substitutions can be tuned to improve the properties of monomers (such as emission energy and reorganization energy). We hope that our work will shine some light on Pt(ii) emitters in the fabrication of efficient OLEDs.

Efficient Density-Fitted Explicitly Correlated Dispersion and Exchange Dispersion Energies

The leading-order dispersion and exchange-dispersion terms in symmetry-adapted perturbation theory (SAPT), E_disp⁽²⁰⁾ and E_exch-disp⁽²⁰⁾, suffer from slow convergence to the complete basis set limit. To alleviate this problem, explicitly correlated variants of these corrections, E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12, have been proposed recently. However, the original formalism (M., Kodrycka , J. Chem. Theory Comput. 2019, 15, 5965-5986), while highly successful in terms of improving convergence, was not competitive to conventional orbital-based SAPT in terms of computational efficiency due to the need to manipulate several kinds of two-electron integrals. In this work, we eliminate this need by decomposing all types of two-electron integrals using robust density fitting. We demonstrate that the error of the density fitting approximation is negligible when standard auxiliary bases such as aug-cc-pVXZ/MP2FIT are employed. The new implementation allowed us to study all complexes in the A24 database in basis sets up to aug-cc-pV5Z, and the E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12 values exhibit vastly improved basis set convergence over their conventional counterparts. The well-converged E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12 numbers can be substituted for conventional E_disp⁽²⁰⁾ and E_exch-disp⁽²⁰⁾ ones in a calculation of the total SAPT interaction energy at any level (SAPT0, SAPT2+3, ...). We show that the addition of F12 terms does not improve the accuracy of low-level SAPT treatments. However, when the theory errors are minimized in high-level SAPT approaches such as SAPT2+3(CCD)δMP2, the reduction of basis set incompleteness errors thanks to the F12 treatment substantially improves the accuracy of small-basis calculations.

Non-Covalent Interactions Atlas Benchmark Data Sets 3: Repulsive Contacts

The new R739×5 data set from the Non-Covalent Interactions Atlas series (www.nciatlas.org) focuses on repulsive contacts in molecular complexes, covering organic molecules, sulfur, phosphorus, halogens, and noble gases. Information on the repulsive parts of the potential energy surface is crucial for the development of robust empirically parametrized computational methods. We use the new data set of highly accurate CCSD(T)/CBS interaction energies to test selected density functional theory (DFT) and semiempirical quantum-mechanical methods. The double-hybrid functionals were the best performing, with the revDSD-PBEP86-D3 being the most accurate DFT method, followed by the range-separated ωB97X functionals. Out of semiempirical methods, GFN2-xTB yielded the best results. On the example of the PM6 method, we analyze the source of error and its relation to the difficulties in the description of conformational energies, and we also devise an immediately applicable correction that fixes the most serious uncorrected issues previously encountered in practical calculations.

The Devil's Triangle of Kohn-Sham density functional theory and excited states

Exchange-correlation (XC) functionals from Density Functional Theory (DFT) developed under the rigorous arguments of Correlated Orbital Theory (COT) address the Devil's Triangle of prominent errors in Kohn-Sham (KS) DFT. At the foundation of this triangle lie the incorrect one-particle spectrum, the lack of integer discontinuity, and the self-interaction error. At the top level, these failures manifest themselves in incorrect charge transfer and Rydberg excitation energies, along with poor activation barriers. Accordingly, the Quantum Theory Project (QTP) XC functionals have been created to address several of the long-term issues encountered in KS theory and its Time Dependent DFT (TDDFT) variant for electronic excitations. Recognizing that COT starts with a correct one-particle spectrum, a condition imposed on the QTP functionals by means of minimum parameterization, the question that arises is how does this affect the electronically excited states? Among up to 28 XC functionals considered, the QTP family provides one of the smallest mean absolute deviations for charge-transfer excitations while also showing excellent results for Rydberg states. However, there is some room for improvement in the case of excitation energies to valence states, which are systematically underestimated by all functionals investigated. An alternative path for better treatment of excitation energies, mainly for valence states, is offered by using orbital energies from QTP functionals, especially by CAM-QTP-02 and LC-QTP. In this case, the deviations from the reference data can be reduced approximately by half.

Electrostatics, Charge Transfer, and the Nature of the Halide-Water Hydrogen Bond

Binary halide-water complexes X^-(H₂O) are examined by means of symmetry-adapted perturbation theory, using charge-constrained promolecular reference densities to extract a meaningful charge-transfer component from the induction energy. As is known, the X^-(H₂O) potential energy surface (for X = F, Cl, Br, or I) is characterized by symmetric left and right hydrogen bonds separated by a C_2v-symmetric saddle point, with a tunneling barrier height that is <2 kcal/mol except in the case of F^-(H₂O). Our analysis demonstrates that the charge-transfer energy is correspondingly small (<2 kcal/mol except for X = F), considerably smaller than the electrostatic interaction energy. Nevertheless, charge transfer plays a crucial role determining the conformational preferences of X^-(H₂O) and provides a driving force for the formation of quasi-linear X··· H-O hydrogen bonds. Charge-transfer energies correlate well with measured O-H vibrational redshifts for the halide-water complexes and also for OH^-(H₂O) and NO₂^-(H₂O), providing some indication of a general mechanism.

Energetics and kinetics of various cyano radical hydrogen abstractions

The cyano radical (CN) is an abundant, open-shell molecule found in a variety of environments, including the atmosphere, the interstellar medium and combustion processes. In these environments, it often reacts with small, closed-shell molecules via hydrogen abstraction. Both carbon and nitrogen atoms of the cyano radical are reactive sites, however the carbon is more reactive with reaction barrier heights generally between 2-15 kcal mol^-1 lower than those of the analogous nitrogen. The CN + HX → HCN/HNC + X, with X = H, CH₃, NH₂, OH, F, SiH₃, PH₂, SH, Cl, C₂H, CN reactions have been studied at a high-level of theory, including CCSD(T)-F12a. Furthermore, kinetics were obtained over the 100-1000 K temperature range, showing excellent agreement with those rate constants that have been determined experimentally.

Interplay between non-covalent interactions in 1D supramolecular polymers based on 1,4-bis(iodoethynyl)benzene

A study of the solution-phase, solid-state structures of halogen-bonded co-crystals of 1,4-bis(iodoethynyl)benzene (p-BIB) with three salts, namely, decyltrimethylammonium bromide (DTMABr), tetrapropylammonium bromide (TPABr), and tetrabutylammonium bromide (TBABr), has been carried out, along with theoretical calculations. Isothermal titration calorimetry (ITC) showed that the binding constant of bromide with p-BIB in THF is not strongly dependent on the cation, and that the entropic term clearly dominates the enthalpic one in the free energy of binding. In the three crystal structures, the bromide anion acts as a doubly connected node for halogen bonding interactions, which results in linear or angular open chains. The intrachain angles (IBr^-I) of the 1D supramolecular polymers based on p-BIB depend on the geometry and size of the cation and vary from 180° for DMTABr to 75° for TBABr. Non-covalent interaction (NCI) analysis of selected motifs and optimized crystals demonstrates that the balance between halogen bonds, hydrogen bonds, and van der Waals interactions, especially type-I halogenhalogen contacts, determines the crystal structures.

Isomerization and Fragmentation Reactions on the [C2SH4] Potential Energy Surface: The Metastable Thione S-Methylide Isomer

Thione S-methylide, parent species of the thiocarbonyl ylide family, is a 1,3-dipolar species on the [C₂SH₄] potential energy surface, not so much studied as its isomers, thiirane, vinyl thiol, and thioacetaldehyde. The conrotatory ring-closure reaction toward thiirane was studied in the 90s, but no complete analysis of the potential energy surface is available. In this paper, we report a computational study of the reaction scheme linking all species. We employed several computational methods (density functional theory, CCSD(T) composite schemes, and CASSCF/CASPT2 multireference procedures) to find the best description of thione S-methylide, its isomers, and transition states. The barrier from thiirane to thione S-methylide amounts to 52.2 kcal mol^-1 (against 17.6 kcal mol^-1 for the direct one), explaining why thiocarbonyl ylides cannot be prepared from thiiranes. Conversion of thiirane to vinyl thiol implies a large barrier, supporting why the reaction has been observed only at high temperatures. Fragmentations of thiirane to S(³P) or S(¹D) and ethylene as well as decomposition to hydrogen sulfide plus acetylene were also explored. Triplet and singlet open-shell species were identified as intermediates in the fragmentations, with energies lower than the transition state between thiirane and vinyl thiol, explaining the preference of the latter at low temperatures.

Quantitative characterisation of the ring normal modes. Pyridine as a study case

In the present work, the vibrational normal modes (NM) of pyridine were revisited. Quantum Chemical calculations were performed to help understand the true nature of some ring related vibrational normal modes (RNM) and how they may be correlated with the electronic structure on the ring. The 27 vibrational normal modes were decomposed into the molecular internal coordinates, and the interest was focused on 7 of them, involving the in-plane ring motion. The electronic structure was analysed through frontier Molecular Orbitals (MO), maps of Molecular Electrostatic Potential surfaces (MEPs) and Natural Bond Orbital (NBO) analysis in a dynamic manner, wherein, each vibration was scanned. The present investigation is aimed to provide the Reader with a quantitative characterisation of the RNMs of pyridine.

London Dispersion Interactions Rather than Steric Hindrance Determine the Enantioselectivity of the Corey-Bakshi-Shibata Reduction

The well-known Corey-Bakshi-Shibata (CBS) reduction is a powerful method for the asymmetric synthesis of alcohols from prochiral ketones, often featuring high yields and excellent selectivities. While steric repulsion has been regarded as the key director of the observed high enantioselectivity for many years, we show that London dispersion (LD) interactions are at least as important for enantiodiscrimination. We exemplify this through a combination of detailed computational and experimental studies for a series of modified CBS catalysts equipped with dispersion energy donors (DEDs) in the catalysts and the substrates. Our results demonstrate that attractive LD interactions between the catalyst and the substrate, rather than steric repulsion, determine the selectivity. As a key outcome of our study, we were able to improve the catalyst design for some challenging CBS reductions.

Predicting Absorption and Emission Maxima of Polycyclic Aromatic Azaborines: Reliable Transition Energies and Character

Polycyclic aromatic azaborines have potential applications as luminophores, novel fluorescent materials, organic light-emitting diodes, and fluorescent sensors. Additionally, their relative structural simplicity should allow the use of computational techniques to design and screen novel compounds in a rapid manner. Herein, the absorption and emission maxima of twelve polycyclic aromatic BN-1,2-azaborine analogues containing the N-BOH moiety were examined to determine a methodology for reliably predicting both the energy and character (local excitation [LE] vs charge transfer [CT]) of the absorption and emission maxima for these compounds. The necessity of implicit solvation models was also investigated. The cam-QTP(01) functional with a small, double-ζ quality basis set provides reliable data compared to EOM-CCSD/cc-pVDZ single-point computations. Of note, commonly used functionals for these applications (B3LYP and ωB97xD) struggle to provide reliable results for both the energy and LE character of the transitions relative to EOM-CCSD computations.

Extended gate-type organic transistor functionalized by molecularly imprinted polymer for taurine detection

Molecularly imprinted polymers (MIPs) are a fascinating technology for the sensitive and selective detection of target molecules. However, in most situations, the need for complicated and expensive analytical devices for reading the responses of MIPs greatly limits their applications. For exploring low-cost and easy-to-use applications of MIPs, herein we have developed a MIP-modified extended-gate type organic field-effect transistor (MIP-OFET). Taurine was selected as a demonstrative analyte due to its biological roles and utility as a nutrient. We explored the rational design of the novel MIP with the aid of density functional theory and wave function calculations and characterized the electrochemically synthesized MIP using differential pulse voltammetry and electrochemical impedance spectroscopy. The mechanism of taurine detection by the MIP-OFET can be explained by the changes in the surface potential of the MIP-functionalized extended-gate electrode accompanied with the capture of taurine. The detection limit of taurine in complete aqueous media was estimated to be 0.33 μM, which was lower or comparable to those calculated by high-performance liquid chromatography. Furthermore, taurine in a commercial drink without any extraction was also successfully detected using the fabricated MIP-OFET. This study would broaden the scope of the applications of MIP-OFETS as chemical sensors for on-site detection of various daily nutrients.

A new anthraquinone derivative as a near UV and visible light photoinitiator for free-radical, thiol–ene and cationic polymerizations

We report a new photoinitiator Q4 which has excellent initiated properties, and can initiate diverse photopolymerization modes. This novel photoinitiator may find use in many specific curing applications due to its unique performance.

Lactic Acid Spectroscopy: Intra- and Intermolecular Interactions

Lactic acid, a relevant molecule in biology and the environment, is an α-hydroxy acid with a high propensity to form hydrogen bonds, both internally and to other hydrogen-bond-accepting molecules. This work includes the novel recording of infrared spectra of gas-phase lactic acid using Fourier transform infrared spectroscopy, and the vibrational absorption features of lactic acid are assigned with the aid of computationally simulated vibrational spectra with anharmonic corrections. Theoretical chemistry methods are used to relate intramolecular hydrogen-bond strengths to the relative stability of lactic acid conformers. The formation of hydrogen-bonded lactic acid dimers and 1:1 water complexes is investigated by simulated vibrational spectra and calculated thermodynamic parameters for the lactic acid monomer and dimer and its water complex in the gas phase. The results of this study are discussed in the context of environmental chemistry with an emphasis on indoor environments.

Comment on "Exploring nature and predicting strength of hydrogen bonds: A correlation analysis between atoms-in-molecules descriptors, binding energies, and energy components of symmetry-adapted perturbation theory"

We evaluate the correlation between binding energy (BE) and electron density ρ(r) at the bond critical point for 28 neutral hydrogen bonds, recently reported by Emamian and co-workers (J. Comput. Chem., 2019, 40, 2868). As an efficient tool, we use local stretching force constant k HB a derived from the local vibrational mode theory of Konkoli and Cremer. We compare the physical nature of BE versus k HB a , and provide an important explanation for cases with significant deviation in the BE- k HB a relation as well as in the BE-ρ(r) correlation. We also show that care has to be taken when different hydrogen bond strength measures are compared. The BE is a cumulative hydrogen bond strength measure while k HB a is a local measure of hydrogen bond strength covering different aspects of bonding. A simplified and unified description of hydrogen bonding is not always possible and needs an in-depth understanding of the systems involved.

Electron-Passing Neural Networks for Atomic Charge Prediction in Systems with Arbitrary Molecular Charge

Atomic charges are critical quantities in molecular mechanics and molecular dynamics, but obtaining these quantities requires heuristic choices based on atom typing or relatively expensive quantum mechanical computations to generate a density to be partitioned. Most machine learning efforts in this domain ignore total molecular charges, relying on overfitting and arbitrary rescaling in order to match the total system charge. Here, we introduce the electron-passing neural network (EPNN), a fast, accurate neural network atomic charge partitioning model that conserves total molecular charge by construction. EPNNs predict atomic charges very similar to those obtained by partitioning quantum mechanical densities but at such a small fraction of the cost that they can be easily computed for large biomolecules. Charges from this method may be used directly for molecular mechanics, as features for cheminformatics, or as input to any neural network potential.

Compressed-State Multistate Pair-Density Functional Theory

Multiconfiguration pair-density functional theory (MC-PDFT) is a multireference method that can be used to calculate excited states. However, MC-PDFT potential energy surfaces have the wrong topology at conical intersections because the last step of MC-PDFT is not a diagonalization of a model-space Hamiltonian matrix, as done in, for example, multistate second-order perturbation theory (MS-CASPT2). We have previously proposed methods that solve this problem by diagonalizing a model-space effective Hamiltonian matrix, where the diagonal elements are MC-PDFT energies for intermediate states, and the off-diagonal elements are evaluated by wave function theory. One previous method is called variational multistate PDFT (VMS-PDFT), whose intermediate states maximize the trace of the effective Hamiltonian, namely, the sum of the MC-PDFT energies of the model-space states; the VMS-PDFT is very robust but is more computationally expensive than another method, extended multistate PDFT (XMS-PDFT), in which the transformation to intermediate states is accomplished without needing any density functional evaluations. However, although VMS-PDFT was accurate in all cases tested, XMS-PDFT was accurate in only some of them. In the present paper, we propose a new method, called compressed-state multistate PDFT (CMS-PDFT), that is as efficient as XMS-PDFT and as accurate as VMS-PDFT. The new method maximizes the trace of the classical Coulomb energy of the intermediate states such that the electron densities of the intermediate states are compressed. We show that CMS-PDFT performs robustly even where XMS-PDFT fails.

Mechanistic Study of the Production of NOx Gases from the Reaction of Copper with Nitric Acid

Copper dissolution in nitric acid is a historic reaction playing a central role in many industrial processes, particularly for metal recovery from the electronics to nuclear industries. The mechanism through which this process occurs is debated. In order to better understand this process, quantum chemical calculations were performed to elucidate the key steps in the mechanism of copper dissolution in nitric acid. We combine both Kohn-Sham density functional theory and ab initio molecular dynamics simulations to understand the mechanism of the formation of the key products: NO₂, HNO₂, and NO. Our calculations suggest that the mechanisms of formation of NO₂, HNO₂, and NO are interconnected.

Multi-state pair-density functional theory

Multi-configuration pair-density functional theory (MC-PDFT) has previously been applied successfully to carry out ground-state and excited-state calculations. However, because they include no interaction between electronic states, MC-PDFT calculations in which each state's PDFT energy is calculated separately can give an unphysical double crossing of potential energy surfaces (PESs) in a region near a conical intersection. We have recently proposed state-interaction pair-density functional theory (SI-PDFT) to treat nearly degenerate states by creating a set of intermediate states with state interaction; although this method is successful, it is inconvenient because two SCF calculations and two sets of orbitals are required and because it puts the ground state on an unequal footing with the excited states. Here we propose two new methods, called extended-multi-state-PDFT (XMS-PDFT) and variational-multi-state-PDFT (VMS-PDFT), that generate the intermediate states in a balanced way with a single set of orbitals. The former uses the intermediate states proposed by Granovsky for extended multi-configuration quasi-degenerate perturbation theory (XMC-QDPT); the latter obtains the intermediate states by maximizing the sum of the MC-PDFT energies for the intermediate states. We also propose a Fourier series expansion to make the variational optimizations of the VMS-PDFT method convenient, and we implement this method (FMS-PDFT) both for conventional configuration-interaction solvers and for density-matrix-renormalization-group solvers. The new methods are tested for eight systems, exhibiting avoided crossings among two to six states. The FMS-PDFT method is successful for all cases for which it has been tested (all cases in this paper except O3 for which it was not tested), and XMS-PDFT is successful for all eight cases except the mixed-valence case. Since both XMS-PDFT and VMS-PDFT are less expensive than XMS-CASPT2, they will allow well-correlated calculations on much larger systems for which perturbation theory is unaffordable.

Understanding benzyl alcohol aggregation by chiral modification: the pairing step

A combination of linear infrared and Raman spectroscopy in supersonic slit jet expansions is used to clarify the conformational preferences in the dimer of the transiently chiral benzyl alcohol (phenylmethanol) under vacuum isolation. By experimentally exploring close analogies with the permanently chiral 1-phenylethanol, which is conformationally locked in the jet through intramolecular chirality induction, and by computational analysis of their conformational energy landscapes, several conclusions are drawn. The lowest energy dimer is confirmed to be cooperatively OHOHπ-bonded and shown to be homochiral. Its heterochiral counterpart is slightly higher in energy and can be spectrally assigned as a minor constituent. A metastable heterochiral OHπ/OHπ structure with weakly coupled hydrogen bonds is efficiently trapped behind a Ci symmetry-enhanced barrier and can be assigned by IR/Raman mutual exclusion. Its homochiral counterpart is kinetically less stable but might be addressed by rotational spectroscopy. Ratings of standard density functionals with a standard basis set in terms of reproducing these experimental chirality synchronization benchmarks are presented.