The 50 most-cited articles in The Journal of Physical Chemistry A, 2022

This Feature Article starts highlighting some recent experimental and theoretical advances in the field of IR and Raman spectroscopy, giving a taste of the breadth and dynamics of this striving field. The local mode theory is then reviewed, showing how local vibrational modes are derived from fundamental normal modes. New features are introduced that add to current theoretical efforts: (i) a unique measure of bond strength based on local mode force constants ranging from bonding in single molecules in different environments to bonding in periodic systems and crystals and (ii) a new way to interpret vibrational spectra by pinpointing and probing interactions between particular bond stretching contributions to the normal modes. All of this represents a means to work around the very nature of normal modes, namely that the vibrational motions in polyatomic molecules are delocalized. Three current focus points of the local mode analysis are reported, demonstrating how the local mode analysis extracts important information hidden in vibrational spectroscopy data supporting current…

Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent…

High Resolution Image Download MS PowerPoint Slide Recent experimental and theoretical studies have shown several new organic molecules that violate Hund’s rule and have the first singlet excited state lower in energy than the first triplet excited state. While many correlated single reference wave function methods have successfully predicted excited-state energetics of these low-lying states, conventional linear-response time-dependent density functional theory (TDDFT) fails to predict the correct excited-state energy ordering. In this article, we have explored the performance of combined DFT and wave function methods like doubles-corrected TDDFT and multiconfiguration pair-density functional theory for the calculation of inverted singlet–triplet gaps. We have also tested the performance of the excited-state DFT (eDFT) method for this problem. Our results have shown that it is possible to obtain inverted singlet–triplet gaps both by using doubles-corrected TDDFT with a proper choice of double-hybrid functionals or by using eDFT.

Machine learning has become more and more popular in computational chemistry, as well as in the important field of spectroscopy. In this concise review, we walk the reader through a short summary of machine learning algorithms and a comprehensive discussion on the connection between machine learning methods and vibrational spectroscopy, particularly for the case of infrared and Raman spectroscopy. We also briefly discuss state-of-the-art molecular representations which serve as meaningful inputs for machine learning to predict vibrational spectra. In addition, this review provides an overview of the transferability and best practices of machine learning in the prediction of vibrational spectra as well as possible future research directions.

Quantitative estimates of reaction barriers are essential for developing kinetic mechanisms and predicting reaction outcomes. However, the lack of experimental data and the steep scaling of accurate quantum calculations often hinder the ability to obtain reliable kinetic values. Here, we train a directed message passing neural network on nearly 24,000 diverse gas-phase reactions calculated at CCSD(T)-F12a/cc-pVDZ-F12//ωB97X-D3/def2-TZVP. Our model uses 75% fewer parameters than previous studies, an improved reaction representation, and proper data splits to accurately estimate performance on unseen reactions. Using information from only the reactant and product, our model quickly predicts barrier heights with a testing MAE of 2.6 kcal mol –1 relative to the coupled-cluster data, making it more accurate than a good density functional theory calculation. Furthermore, our results show that future modeling efforts to estimate reaction properties would significantly benefit from fine-tuning calibration using a transfer learning technique. We anticipate this model will accelerate and…

The power of structural information for informing biological mechanisms is clear for stable folded macromolecules, but similar structure-function insight is more difficult to obtain for highly dynamic systems such as intrinsically disordered proteins (IDPs) which must be described as structural ensembles. Here, we present IDPConformerGenerator, a flexible, modular open-source software platform for generating large and diverse ensembles of disordered protein states that builds conformers that obey geometric, steric, and other physical restraints on the input sequence. IDPConformerGenerator samples backbone phi (φ), psi (ψ), and omega (ω) torsion angles of relevant sequence fragments from loops and secondary structure elements extracted from folded protein structures in the RCSB Protein Data Bank and builds side chains from robust Monte Carlo algorithms using expanded rotamer libraries. IDPConformerGenerator has many user-defined options enabling variable fractional sampling of secondary structures, supports Bayesian models for assessing the agreement of IDP ensembles for consistency…

We present the first benchmark set focusing on the conformational energies of highly flexible, long n-alkane chains, termed ACONFL. Unbranched alkanes are ubiquitous building blocks in nature, so the goal is to be able to calculate their properties most accurately to improve the modeling of, e.g, complex (biological) systems. Very accurate DLPNO-CCSD(T1)/CBS reference values are provided, which allow for a statistical meaningful evaluation of even the best available density functional methods. The performance of established and modern (dispersion corrected) density functionals is comprehensively assessed. The recently introduced r²SCAN-V functional shows excellent performance, similar to efficient composite DFT methods like B97-3c and r²SCAN-3c, which provide an even better cost-accuracy ratio, while almost reaching the accuracy of much more computationally demanding hybrid or double hybrid functionals with large QZ AO basis sets. In addition, we investigated the performance of common wavefunction methods, where MP2/CBS surprisingly performs worse compared to simple D4 dispersion…

We have conducted an extensive search for nitrogen-, oxygen-, and sulfur-bearing heterocycles toward Taurus Molecular Cloud 1 (TMC-1) using the deep, broadband centimeter-wavelength spectral line survey of the region from the GOTHAM large project on the Green Bank Telescope. Despite their ubiquity in terrestrial chemistry, and the confirmed presence of a number of cyclic and polycyclic hydrocarbon species in the source, we find no evidence for the presence of any heterocyclic species. Here, we report the derived upper limits on the column densities of these molecules obtained by Markov Chain Monte Carlo (MCMC) analysis and compare this approach to traditional single-line upper limit measurements. We further hypothesize why these molecules are absent in our data, how they might form in interstellar space, and the nature of observations that would be needed to secure their detection.

Singlet–triplet gap (Δ E S – T ) and spin–orbit coupling (SOC) primarily govern intersystem crossing (ISC)-mediated photo- and electro-luminescence processes. Structural-twist in organic molecules is known to improve ISC efficiency. However, how and to what extent a twist affects the Δ E S – T and SOC are not yet fully understood. In this work, the impact of molecular-twist on these energetics governing ISC is unveiled in a series of highly fluorescent prototype perylenediimides (PDIs) in dichloromethane implementing reliable quantum-chemical calculations. While S 1 → T 1 ISC remains suppressed with increasing twist, a relatively larger decrease in Δ E S – T together with a modest increase in SOC results in enhanced S 1 → T 2 ISC. Significantly modulated ISC rates are predicted in a few experimentally relevant −CN- and −Br-substituted PDIs, where twist of varied extent arises naturally depending on substituent’s chemical nature, numbers, and positions. This study uncovers the critical role of molecular-twist in tailoring ISC and thereby helps designing functional organic…

The characteristics of the CH stretching and out-of-plane bending modes in polycyclic aromatic hydrocarbon molecules are investigated using anharmonic density functional theory (DFT) coupled to a vibrational second-order perturbation treatment taking resonance effects into account. The results are used to calculate the infrared emission spectrum of vibrationally excited species in the collision-less environment of interstellar space. This model follows the energy cascade as the molecules relax after the absorption of a UV photon in order to calculate the detailed profiles of the infrared bands. The results are validated against elegant laboratory spectra of polycyclic aromatic hydrocarbon absorption and emission spectra obtained in molecular beams. The factors which influence the peak position, spectral detail, and relative strength of the CH stretching and out-of-plane bending modes are investigated, and detailed profiles for these modes are derived. These are compared to observations of astronomical objects in space, and the implications for our understanding of the…

In this paper, a multiscale physics-informed neural network (MPINN) approach is proposed based on the regular physics-informed neural network (PINN) for solving stiff chemical kinetic problems with governing equations of stiff ordinary differential equations (ODEs). In MPINNs, chemical species with different time scales are grouped and trained by multiple corresponding neural networks with the same structure. The adaptive weight based on a key performance indicator is assigned to each loss term when calculating the summation of loss residues. With this structure, MPINNs provide a framework to solve challenging stiff chemical kinetic problems without any stiffness-removal artifacts before training. In addition, by introducing a small number of ground truth data (GTD) points (less than 10% of the number required for residual loss calculation) and adding data loss terms into loss functions, MPINNs show superior ability to represent stiff ODE solutions at any desired time. The accuracy of MPINNs is tested with classical chemical kinetic problems, by comparing with the regular PINN and…

A novel method of spectroscopy in ion traps termed leak-out spectroscopy (LOS) is presented. Here, mass-selected, cold ions are excited by an infrared laser. In a subsequent collision with a neutral buffer gas particle, their internal energy is then transferred to kinetic energy. As a result, these ions leak out from the ion trap and are detected. The LOS scheme is generally applicable, very sensitive, and close to background-free when operated at low temperature. The potential of this method is demonstrated and characterized here for the first time by recording the rotationally resolved spectrum of the C–H stretching vibration ν 1 of linear C 3 H + . Besides performing high-resolution spectroscopy, this method opens up the way for analyzing the composition of trap content, for example, determining isomer ratios, by selectively expelling isomers or other isobaric ions from the trap. Likewise, LOS can be used to prepare clean samples of structural and nuclear spin isomers.

Dark chamber experiments were conducted to study the SOA formed from the oxidation of α-pinene and Δ-carene under different peroxy radical (RO 2 ) fate regimes: RO 2 + NO 3, RO 2 + RO 2, and RO 2 + HO 2 . SOA mass yields from α-pinene oxidation were <1 to ∼25% and strongly dependent on available OA mass up to ∼100 μg m –3 . The strong yield dependence of α-pinene oxidation is driven by absorptive partitioning to OA and not by available surface area for condensation. Yields from Δ-carene + NO 3 were consistently higher, ranging from ∼10–50% with some dependence on OA for <25 μg m –3 . Explicit kinetic modeling including vapor wall losses was conducted to enable comparisons across VOC precursors and RO 2 fate regimes and to determine atmospherically relevant yields. Furthermore, SOA yields were similar for each monoterpene across the nominal RO 2 + NO 3, RO 2 + RO 2, or RO 2 + HO 2 regimes; thus, the volatility basis sets (VBS) constructed were independent of the chemical regime. Elemental O/C ratios of ∼0.4–0.6 and nitrate/organic mass ratios of ∼0.15 were observed in the particle…

Resistively heated silicon carbide microreactors are widely applied as continuous sources to selectively prepare elusive and reactive intermediates with astrochemical, catalytic, or combustion relevance to measure their photoelectron spectrum. These reactors also provide deep mechanistic insights into uni- and bimolecular chemistry. However, the sampling conditions and effects have not been fully characterized. We use cation velocity map imaging to measure the velocity distribution of the molecular beam signal and to quantify the scattered, rethermalized background sample. Although translational cooling is efficient in the adiabatic expansion from the reactor, the breakdown diagrams of methane and chlorobenzene confirm that the molecular beam component exhibits a rovibrational temperature comparable with that of the reactor. Thus, rovibrational cooling is practically absent in the expansion from the microreactor. The high rovibrational temperature also affects the threshold photoelectron spectrum of both benzene and the allyl radical in the molecular beam, but to different degrees.…

Chemical transformations in aerosol impact the lifetime of particle-phase species, the fate of atmospheric pollutants, and both climate and health-relevant aerosol properties. Timescales for multiphase reactions of ozone in atmospheric aqueous phases are governed by coupled kinetic processes between the gas-phase, the particle interface and its bulk, which respond dynamically to reactive consumption of O3. However, models of atmospheric aerosol reactivity often do not account for the coupled nature of multiphase processes. To examine these dynamics, we use new and prior experimental observations of aqueous droplet reaction kinetics, including three systems with a range of surface affinities and ozonolysis rate coefficients (trans-aconitic acid (C6H6O6), maleic acid (C4H4O4) and sodium nitrite (NaNO2)). Using literature rate coefficients and thermodynamic properties, we constrain a simple two-compartment stochastic kinetic model, which resolves the interface from the particle bulk and represents O3 partitioning, diffusion and reaction as a coupled kinetic system. Our kinetic model…

An extension of the exact two-component theory with atomic mean-field integrals (the X2CAMF scheme) to the treatment of the Breit term together with efficient implementation using an atomic Dirac-Coulomb-Breit Hartree-Fock program is reported. The accuracy of the X2CAMF scheme for treating the contributions from the Breit term to the molecular properties is demonstrated using benchmark calculations of equilibrium bond lengths, harmonic frequencies, and dipole moments for molecules containing elements across the periodic table. Calculations of the properties for molecules containing period four elements aiming at high accuracy as well as for Th- and U-containing molecules are also presented and compared with experimental results to demonstrate the usefulness of the X2CAMF scheme in combination with accurate treatments of electron correlation by the coupled-cluster (CC) methods. The combination of CC methods and the X2CAMF scheme shows potential to extend the accuracy of CC calculations to heavy elements, e.g., to computational heavy-element thermochemistry.

SCAN-3c implementation with STOs is a promising alternative to the original implementation with GTOs and can be generally used for a broad field of quantum chemical problems.

values of modified nucleobases give insight into their structural and energetic impacts within nucleic acids.

electronic structure calculations, minimum-energy path refinements, and transition state searches. Using representative examples from homogeneous catalysis and interstellar chemistry, we highlight how these schemes increasingly act as "virtual reaction vessels" for interrogating mechanistic questions. Finally, we highlight where challenges remain, including issues of chemical accuracy and calculation speeds, as well as the inherent challenge of dealing with the vast size of accessible chemical reaction space.

To better understand the thermochemical kinetics and mechanism of a specific chemical reaction, an accurate estimation of barrier heights (forward and reverse) and reaction energies is vital. Because of the large size of reactants and transition state structures involved in real-life mechanistic studies (e.g., enzymatically catalyzed reactions), density functional theory remains the workhorse for such calculations. In this paper, we have assessed the performance of 91 density functionals for modeling the reaction energies and barrier heights on a large and chemically diverse data set (BH9) composed of 449 organic chemistry reactions. We have shown that range-separated hybrid functionals perform better than the global hybrids for BH9 barrier heights and reaction energies. Except for the PBE-based range-separated nonempirical double hybrids, range separation of the exchange term helps improve the performance for barrier heights and reaction energies. The 16-parameter Berkeley double hybrid, ωB97M(2), performs remarkably well for both properties. However, our minimally empirical…

High Resolution Image Download MS PowerPoint Slide Thermally activated delayed fluorescence (TADF) materials are commonly used in various apparatus, including organic light-emitting device-based displays, as they remarkably improve the internal quantum efficiencies. Although there is a wide range of donor–acceptor-based compounds possessing TADF properties, in this computational study, we investigated TADF and some non-TADF chromophores, containing benzophenone or its structural derivatives as the acceptor core, together with various donor moieties. Following the computational modeling of the emitters, several excited state properties, such as the absorption spectra, singlet–triplet energy gaps (Δ E ST ), natural transition orbitals, and the topological Φ S indices, have been computed. Along with the donor–acceptor torsion angles and spin-orbit coupling values, these descriptors have been utilized to investigate potential TADF efficiency. Our study has shown that on the one hand, our photophysical/structural descriptors and computational methodologies predict the experimental…

High Resolution Image Download MS PowerPoint Slide The excitation energy transfer (ET) pathway and mechanism from an organic antenna to a lanthanide ion has been the subject of discussion for many decades. In the case of europium (Eu 3+ ), it has been suggested that the transfer originates from the ligand singlet state or a triplet state. Taking the lanthanide complex Eu(TTA) 3 (H 2 O) 2 as an example, we have investigated the spectra and luminescence kinetics, mainly at room temperature and 77 K, to acquire the necessary experimental data. We put forward an experimental and theoretical approach to measure the energy transfer rates from the antenna to different Eu 3+ levels using the Dexter formulation. We find that transfer from the ligand singlet state to Eu 3+ may account for the ET pathway, by combined electric dipole–electric dipole (ED–ED) and ED-electric quadrupole (EQ) mechanisms. The contributions from the triplet state by these mechanisms are very small. An independent systems rate equation approach can effectively model the experimental kinetics results. The model…

Assessing and improving the performance of organic light-emitting diode (OLED) materials require quantitative prediction of rate coefficients for the intersystem crossing (ISC) and reverse ISC (RISC) processes, which are determined not only by the energy gap and the direct spin–orbit coupling (SOC) between the first singlet and triplet excited-states at a thermal equilibrium position of the initial electronic state but also by the non-Condon effects such as the Herzberg–Teller-like vibronic coupling (HTVC) and the spin-vibronic coupling (SVC). Here we apply the time-dependent correlation function approaches to quantitatively calculate the vibrationally resolved absorption and fluorescence spectra and ISC/RISC rates of a newly synthesized multiple-resonance-type (MR-type) thermally activated delayed fluorescence (TADF) emitter, 7-phenylquinolino[3,2,1- de ]acridine-5,9-dione (7-PhQAD), with the inclusion of the Franck–Condon (FC), HTVC, and Duschinsky rotation (DR) effects. The SVC effect on the rates has also been approximately evaluated. We find that the experimentally measured ISC…

The enantiopurification of racemic mixtures of chiral molecules is important for a range of applications. Recent work has shown that chiral group-directed photoisomerization is a promising approach to enantioenrich racemic mixtures of BINOL, but increased control of the diasteriomeric excess ( de ) is necessary for its broad utility. Here we develop a cavity quantum electrodynamics (QED) generalization of time-dependent density functional theory and demonstrate computationally that strong light–matter coupling can alter the de of the chiral group-directed photoisomerization of BINOL. The relative orientation of the cavity mode polarization and the molecules in the cavity dictates the nature of the cavity interactions, which either enhance the de of the ( R )-BINOL diasteriomer (from 17% to ≈40%) or invert the favorability to the ( S )-BINOL derivative (to ≈34% de ). The latter outcome is particularly remarkable because it indicates that the preference in diasteriomer can be influenced via orientational control, without changing the chirality of the directing group. We demonstrate…

Using perturbation theory within the framework of conceptual density functional theory, we derive a lower bound for the lattice energy of the ionic solids. The main element of the lower bound is the Fukui potential in the nuclei of the molecule corresponding to the unit formula of the solid. Thus, we propose a model to calculate the lattice energy in terms of the Fukui potential. Our method, which is extremely simple, performs well as other methods using the crystal structure information of alkali halide solids. The method proposed here correlates surprisingly well with the experimental data on the lattice energy of a diverse series of solids having even a non-negligible covalent characteristic. Finally, the validity of the maximum hardness principle (MHP) is assessed, showing that in this case, the MHP is limited.

The broader availability of cost-effective methodologies like second-order vibrational perturbational theory (VPT2), also in general-purpose quantum chemical programs, has made the inclusion of anharmonic effects in vibrational calculations easier, paving the way to more accurate simulations. Combined with modern computing hardware, VPT2 can be used on relatively complex molecular systems with dozen of atoms. However, the problem of resonances and their corrections remains a critical pitfall of perturbative methods. Recent works have highlighted the sensitivity of band intensities to even subtle resonance effects, underlying the importance of a correct treatment to predict accurate spectral bandshapes. This aspect is even more critical with chiroptical spectroscopies whose signal is weak. This has motivated the present work in exploring robust methods and criteria to identify resonances not only in energy calculations but also on the transition moments. To study their performance, three molecules of representative sizes ranging from ten to several dozens of atoms were chosen. The…

Recent ambient atmospheric measurements have detected highly oxygenated organic molecules (HOMs) at many sites and are a consequence of autoxidation processes occurring at ambient temperatures. Monoterpenes in particular have a propensity to autoxidize although they exhibit a wide range of HOM yields, which may be due to a variety of reasons including reactions with different oxidants like OH and O 3, differing hydrogen (H) atom transfer or peroxy radical cyclization rates, numbers of available reaction pathways, and/or energy loss processes for activated HO-monoterpene or O 3 -monoterpene adducts. In this work, the autoxidation mechanisms of (+)-α-pinene, (+)-β-pinene, and (+)-limonene following initial OH oxidation and three successive O 2 additions are examined using density functional theory (DFT) to understand what accounts for the disparity. Rates of different potential autoxidation pathways initiated by OH addition or abstraction reactions are quantified using transition-state theory (TST) and master equation approaches using the lowest-energy conformers. OH abstraction…

In this work, we report the benchmark binding energies of the seven complexes within the L7 data set, six host–guest complexes from the S12L data set, a C 60 dimer, the DNA–ellipticine intercalation complex, and the largest system of the study, the HIV–indinavir system, which contained 343 atoms or 139 heavy atoms. The high-quality values reported were obtained via a focal point method that relies on the canonical form of second-order Møller–Plesset theory and the domain-based local pair natural orbital scheme for the coupled cluster with single double and perturbative triple excitations [DLPNO-CCSD(T)] extrapolated to the complete basis set (CBS) limit. The results in this work not only corroborate but also improve upon some previous benchmark values for large noncovalent complexes albeit at a relatively steep cost. Although local CCSD(T) and the largely successful fixed-node diffusion Monte Carlo (FN-DMC) have been shown to generally agree for small- to medium-size systems, a discrepancy in their reported binding energy values arises for large complexes, where the magnitude of the…

Mercury, a highly toxic metal, is emitted to the atmosphere mostly as gaseous Hg(0). Atmospheric Hg(0) enters ecosystems largely through via uptake by vegetation, while Hg(II) largely enters ecosystems in oceans and via rainfall. Consequently, the redox chemistry of atmospheric mercury strongly influences its fate and its global biogeochemical cycling. Here we report on the oxidation and reduction of Hg(I) (BrHg and HOHg radicals) in reactions with ozone, and how the electronic structure of these Hg(I) species affects the kinetics of these reactions. The oxidation reactions lead to XHgO• + O2 (X=Br and OH), while the reduction reaction produces Hg(0), HOX, and O2. According to our calculations with CCSD(T), NEVPT2, and CAM-B3LYP-D3BJ, the kinetics of both oxidation reactions are very similar. These two oxidation reactions are much faster than their reduction counterparts, and this effect is remarkably stronger for the oxidation of HOHg(I) by ozone. Modeling of field data supports the idea that OH and/or O3 (rather than Br) dominates Hg(II) production in the continental boundary…

A highly accurate CO 2 ab initio dipole moment surface (DMS), Ames-2021, is reported along with 12 C 16 O 2 infrared (IR) intensity comparisons approaching a 1–4‰ level of agreement and uncertainty. The Ames-2021 DMS was accurately fitted from CCSD(T) finite-field dipoles computed with the aug-cc-pVXZ (X = T, Q, 5) basis for C atom and the d-aug-cc-pVXZ (X = T, Q, 5) basis for O atoms, and extrapolated to the one particle basis set limit. Fitting σ rms is 3.8 × 10 –7 au for 4443 geometries below 15 000 cm –1 . The corresponding IR intensity, S Ames-2021, are computed using the Ames-2 potential energy surface (PES), which is the best PES available for CO 2 . Compared to high accuracy IR studies for 2001 i –00001 and 3001 i –00001 bands, S Ames-2021 matches NIST experiment-based intensities [ S NIST-HIT16 or S HIT20 ] to −1.0 ± 1.3‰, or matches DLR experiment-based intensities [ S DLR-HIT16/UCL/Ames ] to 1.9 ± 3.7‰. This indicates the systematic deviations and uncertainties have been significantly reduced in S Ames-2021 . The S UCL2015 (or S HITRAN2016 ) have larger deviations (vs S…

CrSi n – ( n = 14–18) cluster anions have been investigated by a combination of photoelectron spectroscopy (PES) and first-principles calculations. The lowest-lying structures of the clusters have been determined by a global minimum search based on the genetic algorithm, combined with density functional theory (DFT) calculations. The simulated PES spectra of the lowest-energy isomers are in agreement with the experimental results, which gives strong evidence that the correct structures have been found. While sizes n = 14 and n = 15 prefer cage-like structures based on multi-center bonding within the cage, the larger sizes adopt structures based on fullerene-type cages around the Cr atom, with the additional atoms attached to the cage surface. A Hirshfeld analysis shows that the Cr atoms act as electron donors in all clusters, thus enhancing the electron count in the cage. It also reveals that the magnetic moment of 1μ B shown by all clusters is mainly contributed by the Cr atom. One interesting exception is size 17, where the Cr atom contributes a small moment antiparallel to that…

Carbon-doped aluminum cluster anions, Al n C – ( n = 6–15), were generated by laser vaporization and investigated by mass-selected anion photoelectron spectroscopy. The geometric structures of Al n C – ( n = 6–15) anions were determined by the comparison of theoretical calculations with the experimental results. It is found that the most stable structure of Al 6 C – is a carbon endohedral triangular prism. The Al 7 C – anion is a magic cluster with high stability. The structures of Al 7-9 C – can be viewed as the additional aluminum atoms attached around the triangular prism Al 6 C – . Two isomers of Al 10 C – have been detected in the experiments. The most stable one has a planar tetracoordinate carbon structure. The second one derives from Al 9 C – with the carbon atom located in a pentagonal bipyramid. The Al 11 C – anion has a bilayer structure composed of one planar tetracoordinate carbon and one aluminum-centered hexagon, in which the major interactions between two layers are multicenter bonds. The structures of Al 12-14 C – can be viewed as evolving from Al 11 C – by adding…

High Resolution Image Download MS PowerPoint Slide The ability to tune the optoelectronic properties of quantum dots (QDs) makes them ideally suited for the use as fluorescence sensing probes. The vast structural diversity in terms of the composition and size of QDs can make designing a QD for a specific sensing application a challenging process. Quantum chemical calculations have the potential to aid this process through the characterization of the properties of QDs, leading to their in silico design. This is explored in the context of QDs for the fluorescence sensing of dopamine based upon density functional theory and time-dependent density functional theory (TDDFT) calculations. The excited states of hydrogenated carbon, silicon, and germanium QDs are characterized through TDDFT calculations. Analysis of the molecular orbital diagrams for the isolated molecules and calculations of the excited states of the dopamine-functionalized quantum dots establish the possibility of a photoinduced electron-transfer process by determining the relative energies of the electronic states formed…

Entangled two-photon absorption (ETPA) has recently become a topic of lively debate, mainly due to the apparent inconsistencies in the experimentally reported ETPA cross sections of organic molecules obtained by a number of groups. In this work, we provide a thorough experimental study of ETPA in the organic molecules Rhodamine B (RhB) and zinc tetraphenylporphirin (ZnTPP). Our contribution is 3-fold: first, we reproduce previous results from other groups; second, we on the one hand determine the effects of different temporal correlations─introduced as a controllable temporal delay between the signal and idler photons to be absorbed─on the strength of the ETPA signal, and on the other hand, we introduce two concurrent and equivalent detection systems with and without the sample in place as a useful experimental check; third, we introduce, and apply to our data, a novel method to quantify the ETPA rate based on taking into account the full photon-pair behavior rather than focusing on singles or coincidence counts independently. Through this experimental setup we find that,…

Accurate and efficient simulation of liquids, such as water and salt solutions, using high-level wave function theories is still a formidable task for computational chemists owing to the high computational costs. In this study, we develop a deep machine learning potential based on fragment-based second-order Møller–Plesset perturbation theory (DP-MP2) for water through neural networks. We show that the DP-MP2 potential predicts the structural, dynamical, and thermodynamic properties of liquid water in better agreement with the experimental data than previous studies based on density functional theory (DFT). The nuclear quantum effects (NQEs) on the properties of liquid water are also examined, which are noticeable in affecting the structural and dynamical properties of liquid water under ambient conditions. This work provides a general framework for quantitative predictions of the properties of condensed-phase systems with the accuracy of high-level wave function theory while achieving significant computational savings compared to ab initio simulations.

This work presents a detailed investigation of the energy-transfer and dissociation mechanisms in N 2 (X 1 Σ g + ) + O( 3 P) and NO(X 2 Π) + N( 4 S) systems using rovibrational-specific quasiclassical trajectory (QCT) and master equation analyses. The complete set of state-to-state kinetic data, obtained via QCT, allows for an in-depth investigation of the Zel’dovich mechanism leading to the formation of NO molecules at microscopic and macroscopic scales. The master equation analysis demonstrates that the low-lying vibrational states of N 2 and NO have dominant contributions to the NO formation and the corresponding extinction of N 2 through the exchange process. For the considered temperature range, it is found that nearly 50% of the dissociation processes for N 2 and NO molecules occur in the quasi-steady-state (QSS) regime, while for the Zel’dovich reaction, the distribution of the reactants does not reach the QSS conditions. Furthermore, using the QSS approximation to model the Zel’dovich mechanism leads to overestimating NO production by more than a factor of 4 in the…

New approaches for the sensitive and accurate quantification of aerosol optical properties are needed to improve the current understanding of the unique physical chemistry of airborne particles and to explore their roles in fields as diverse as chemical manufacturing, healthcare, and atmospheric science. We have pioneered the use of cavity ring-down spectroscopy (CRDS), with concurrent angularly resolved elastic light scattering measurements, to interrogate the optical properties of single aerosol particles levitated in optical and electrodynamic traps. This approach enables the robust quantification of optical properties such as extinction cross sections for individual particles of known size. Our measurements can now distinguish the scattering and absorption contributions to the overall light extinction, from which the real and imaginary components of the complex refractive indices can be retrieved and linked to chemical composition. In this Feature Article, we show that this innovative measurement platform enables accurate and precise optical measurements for spherical and…

High Resolution Image Download MS PowerPoint Slide Cyclo[ n ]carbons ( n = 5, 7, 9, ..., 29) composed from an odd number of carbon atoms are studied computationally at density functional theory (DFT) and ab initio complete active space self-consistent field (CASSCF) levels of theory to get insight into their electronic structure and aromaticity. DFT calculations predict a strongly delocalized carbene structure of the cyclo[ n ]carbons and an aromatic character for all of them. In contrast, calculations at the CASSCF level yield geometrically bent and electronically localized carbene structures leading to an alternating double aromaticity of the odd-number cyclo[ n ]carbons. CASSCF calculations yield a singlet electronic ground state for the studied cyclo[ n ]carbons except for C 25, whereas at the DFT level the energy difference between the lowest singlet and triplet states depends on the employed functional. The BHandHLYP functional predicts a triplet ground state of the larger odd-number cyclo[ n ]carbons starting from n = 13. Current-density calculations at the BHandHLYP level…

Materials exhibiting higher mobility than conventional organic semiconducting materials, such as fullerenes and fused thiophenes, are in high demand for applications in printed electronics. To discover new molecules that might show improved charge mobility, the adaptive design of experiments (DoE) to design molecules with low reorganization energy was performed by combining density functional theory (DFT) methods and machine learning techniques. DFT-calculated values of 165 molecules were used as an initial training dataset for a Gaussian process regression (GPR) model, and five rounds of molecular designs applying the GPR model and validation via DFT calculations were executed. As a result, new molecules whose reorganization energy is smaller than the lowest value in the initial training dataset were successfully discovered.

Extraction of uranium from seawater is one of the important ways to solve the shortage of terrestrial uranium resources. Thereinto, the competition between uranyl and vanadium cations is a significant challenge in the commonly used amidoxime-based adsorbents for extracting uranium from seawater. An in-depth understanding of the extraction behaviors of modified amidoxime groups with uranyl and vanadium ions is one of the effective means to design and develop efficient adsorbents for selective uranium sequestration. In this work, we have designed and systematically investigated the alkyl and amino functionalized amidoxime, ( Z )-2-amino- N′ -hydroxy- N,N -dimethylbenzimidamide (L 1 ), and its phenyl and methoxy derivatives (( Z )-3-amino- N′ -hydroxy- N,N -dimethyl-2-naphthimidamide (L 2 ) and ( Z )-2-amino- N ′-hydroxy-4-methoxy- N,N -dimethylbenzimidamide (L 3 )) by quantum chemistry calculations. In the uranyl complexes, the amidoxime groups prefer to act as η 2 -coordinated ligands as the amidoximes increase, and there exist substantial hydrogen bond interactions, which are…

High Resolution Image Download MS PowerPoint Slide The recombination (“dimerization”) of peroxyl radicals (RO 2 •) is one of the pathways suggested in the literature for the formation of peroxides (ROOR′, often referred to as dimers or accretion products in the literature) in the atmosphere. It is generally accepted that these dimers play a major role in the first steps of the formation of submicron aerosol particles. However, the precise reaction pathways and energetics of RO 2 • + R′O 2 • reactions are still unknown. In this work, we have studied the formation of tetroxide intermediates (RO 4 R′): their formation from two peroxyl radicals and their decomposition to triplet molecular oxygen ( 3 O 2 ) and a triplet pair of alkoxyl radicals (RO•). We demonstrate this mechanism for several atmospherically relevant primary and secondary peroxyl radicals. The potential energy surface corresponds to an overall singlet state. The subsequent reaction channels of the alkoxyl radicals include, but are not limited to, their dimerization into ROOR′. Our work considers the multiconfigurational…

We explored the hypothesis that on the nanoscale level, acids and bases might exhibit different behavior than in bulk solution. Our study system consisted of sulfuric acid, formic acid, ammonia, and water. We calculated highly accurate Domain-based Local pair-Natural Orbital- Coupled-Cluster/Complete Basis Set (DLPNO-CCSD(T)/CBS) energies on DFT geometries and used the resulting Gibbs free energies for cluster formation to compute the overall equilibrium constants for every possible cluster. The equilibrium constants combined with the initial monomer concentrations were used to predict the formation of clusters at the top and the bottom of the troposphere. Our results show that formic acid is as effective as ammonia at forming clusters with sulfuric acid and water. The structure of formic acid is uniquely suited to form hydrogen bonds with sulfuric acid. Additionally, it can partner with water to form bridges from one side of sulfuric acid to the other, hence demonstrating that hydrogen bonding topology is more important than acid/base strength in these atmospheric prenucleation…

Kinetic compensation is a strong, positive correlation between the Arrhenius activation energy E and the frequency factor A for a reaction between the same reactants under similar experimental conditions or similar reactants under the same conditions, even though these parameters are supposed to be independent. The kinetic compensation effect (KCE) is demonstrated by a linear relationship between ln[ A ] and E / R in the eponymous Constable plot and has been the subject of more than 50 000 publications over the past 100 years, with no consensus opinion about the cause of this effect. In this paper, it is suggested that the linear relationship between ln[ A ] and E is the result of a real or spurious path dependence of the reaction history between the initial state of the pure reactant(s) and the final state of the pure product(s) having standard enthalpy and entropy differences, Δ H ° and Δ S °, respectively. The single-step rate law approximation of a reversible reaction leads to T 0 = H °/Δ S ° as the dynamic thermal (thermodynamic) equilibrium temperature and 1/ T 0 = (ln[ A̅ / k…

Chemical bonds and noncovalent interactions are extraordinarily important concepts in chemistry and beyond. Using density-based quantities to describe them has a long history in the literature, yet none can satisfactorily describe the entire spectrum of interactions from strong chemical bonds to weak van der Waals forces. In this work, employing Pauli energy as the theoretical foundation, we fill in that knowledge gap. Our results show that the newly established density-based index can describe single and multiple covalent bonds, ionic bonds, metallic bonds, and different kinds of noncovalent interactions, all with unique and readily identifiable signature shapes. Two new descriptors, BNI (bonding and noncovalent interaction) index and USI (ultra-strong interaction) index, have been introduced in this work. Together with NCI (noncovalent interaction) and SCI (strong covalent interaction) indexes already available in the literature, a density-based description of both chemical bonds and noncovalent interactions is accomplished.

Machine learning predictions of molecular thermochemistry, such as formation enthalpy, have been limited for large and complicated species because of the lack of available training data. Such predictions would be important in the prediction of reaction thermodynamics and the construction of kinetic models. Herein, we introduce a graph-based deep learning approach that can separately learn the enthalpy contribution of each atom in its local environment with the effect of the overall molecular structure taken into account. Because this approach follows the additivity scheme of increment theory, it can be generalized to larger and more complicated species not present in the training data. By training the model on molecules with up to 11 heavy atoms, it can predict the formation enthalpy of testing molecules with up to 42 heavy atoms with a mean absolute error of 2 kcal/mol, which is less than half of the error of the conventional increment theory. We expect that this approach will also enable rapid prediction of other extensive properties of large molecules that are difficult to derive…

A model is developed to describe trace gas uptake and reaction with applications to aerosols and microdroplets. Gas uptake by the liquid is formulated as a coupled equilibria that links gas, surface and bulk regions of the droplet or solution. Previously, this framework was used in explicit stochastic reaction-diffusion simulations to predict the reactive uptake kinetics of ozone with droplets containing aqueous aconitic acid, maleic acid and sodium nitrite. Using prior data and simulation results, a new equation for the uptake coefficient is derived, which accounts for both surface and bulk reactions. Lambert W functions are used to obtain closed form solutions to the integrated rate laws for the multiphase kinetics; similar to previous expressions that describe Michaelis–Menten enzyme kinetics. Together these equations couple interface and bulk processes over a wide range of conditions and don’t require many of the limiting assumptions needed to apply resistor model formulations to explain trace gas uptake and reaction.

Cells are heterogeneous on every length and time scale; cytosol contains thousands of proteins, lipids, nucleic acids, and small molecules, and molecular interactions within this crowded environment determine the structure, dynamics, and stability of biomolecules. For decades, the effects of crowding at the atomistic scale have been overlooked in favor of more tractable models largely based on thermodynamics. Crowding can affect the conformations and stability of biomolecules by modulating water structure and dynamics within the cell, and these effects are nonlocal and environment dependent. Thus, characterizing water's hydrogen-bond (H-bond) networks is a critical step toward a complete microscopic crowding model. This perspective provides an overview of molecular crowding and describes recent time-resolved spectroscopy approaches investigating H-bond networks and dynamics in crowded or otherwise complex aqueous environments. Ultrafast spectroscopy combined with atomistic simulations has emerged as a powerful combination for studying H-bond structure and dynamics in heterogeneous…

Controlling the distribution of electrons in materials is the holy grail of chemistry and material science. Practical attempts at this feat are common but are often reliant on simplistic arguments based on electronegativity. One challenge is knowing when such arguments work, and which other factors may play a role. Ultimately, electrons move to equalize chemical potentials. In this work, we outline a theory in which chemical potentials of atoms and molecules are expressed in terms of reinterpretations of common chemical concepts and some physical quantities: electronegativity, chemical hardness, and the sensitivity of electronic repulsion and core levels with respect to changes in the electron density. At the zero-temperature limit, an expression of the Fermi level emerges that helps to connect several of these quantities to a plethora of material properties, theories and phenomena predominantly explored in condensed matter physics. Our theory runs counter to Sanderson's postulate of electronegativity equalization and allows a perspective in which electronegativities of bonded atoms…

A recent study has examined the accuracy of NIST heats of formation for a set of C, H, and O-containing species with a proposed low-cost quantum chemistry approach. In the present study, we have used high-level methods such as W1X-2 to obtain these data more rigorously, which we have then used to assess the NIST and the previously computed values. We find that many of these NIST data that are as suggested to be unreliable by the previous study are indeed inconsistent with our high-level reference values. However, we also find substantial deviations for the previously computed values from our benchmark. Thus, we have assessed the performance of alternative low-cost methods. In our assessment, we have additionally examined C, H, N, and O-containing species for which heats of formation are available from the NIST database. We find the ωB97M-V/ma-def2-TZVP, DSD-PBEP86/ma-def2-TZVP, and CCSD(T)-F12b/aug′-cc-pVDZ methods to be adequate for obtaining heats of formation with the atomization approach, once their atomic energies are optimized with our benchmark. Notably, the low-cost ωB97M-V…

The structures of tautomers and rotameric forms of curcumin, the bioactive compound present in spice plant turmeric, have been investigated using ion mobility mass spectrometry (IMMS) in conjunction with high-performance liquid chromatography (HPLC) and UV–visible spectroscopy. Two tautomeric forms of this β-diketone compound, keto–enol and diketo, have been chromatographically separated, and the electronic absorption spectra for these two tautomeric forms in methanol solution have been recorded separately for the first time. The molecular identity of the HPLC-separated solution fractions is established unambiguously by recording the mass and fragmentation spectra simultaneously. The ion mobility spectrum for the deprotonated curcumin anion, [Cur-H] −, corresponding to the diketo tautomer, displays only one peak (P), and the collision cross-section (CCS) value is measured to be 185.9 Å 2 . However, the ion mobility spectrum corresponding to the HPLC-separated keto–enol tautomer shows three distinctly separated peaks, P, Q, and R, with CCS values of 185.9, 194.8, and 203.4 Å 2,…