Chemical interaction between the tips and molecules is one of the main contributing mechanisms to tip-enhanced Raman spectroscopy (TERS). In this work, we calculate the TERS spectra of the biphenylene (BP) dimer at 13 nonequivalent tip sites by means of density functional theory and explore the influence of the TERS tip on vibrational mode characters and Raman intensity. The Raman intensity of the vibrational mode involving the antisymmetric stretching of tetra-rings is found to be specifically enhanced. We attribute this specific enhancement to the electronic sensitive atom vibrational character of the mode and infer that the vibrational strength of atoms can be tuned by the TERS tip. The results provide an intuitive interpretation on the effects of tip-induced electronic redistributions on specific vibrational modes in TERS and indicate the possibility to further improve the TERS resolution.

Recently, the development of machine learning (ML) potentials has made it possible to perform large-scale and long-time molecular simulations with the accuracy of quantum mechanical (QM) models. However, for different levels of QM methods, such as density functional theory (DFT) at the meta-GGA level and/or with exact exchange, quantum Monte Carlo, etc., generating a sufficient amount of data for training an ML potential has remained computationally challenging due to their high cost. In this work, we demonstrate that this issue can be largely alleviated with Deep Kohn–Sham (DeePKS), an ML-based DFT model. DeePKS employs a computationally efficient neural network-based functional model to construct a correction term added upon a cheap DFT model. Upon training, DeePKS offers closely matched energies and forces compared with high-level QM method, but the number of training data required is orders of magnitude less than that required for training a reliable ML potential. As such, DeePKS can serve as a bridge between expensive QM models and ML potentials: one can generate a decent amount of high-accuracy QM data to train a DeePKS model and then use the DeePKS model to label a much larger amount of configurations to train an ML potential. This scheme for periodic systems is implemented in a DFT package ABACUS, which is open source and ready for use in various applications.

Proton transfer modulation in an organic diradical is apparently the most conspicuously attractive phenomenon. In this work, we have computationally designed the trans and cis forms of photochromic azobenzene- (AB-) bridged diradicals by considering AB as coupler and two nitroxide (NO) as spin sources and a −OH attaching at the ortho site as modulator. Our object is that through intramolecular proton transfer to protonate the azo-unit (−N═N−) the magnetic coupling characteristics of the designed diradicals can be modulated in their photocontrolled trans and cis forms. The calculated results indicate that PT can significantly regulate the magnetic spin coupling constants, J = −701.3 cm^–1 ↔ −286.2 cm^–1 for the trans form and −544.1 cm^–1 ↔ −328.1 cm^–1 for the cis form. In particular, we discover that these designed magnetic molecules can undergo magnetic conversion between antiferromagnetic and ferromagnetic coupling through PT, besides there is considerable increase in the magnitude of their magnetic coupling constants J, (e.g., −59.97 to 172.4 cm^–1) for the trans-mode at the m/m linking site. Moreover, we discover that the nitroxide radicals at different linking positions have a significant impact and remarkably alter the magnetic spin coupling characteristics of AB-based diradicals. Besides, various radical groups are used as spin sources which corroborated our assumptions and tended to the same conclusion. This work offers a novel understanding of the spin interaction mechanism and a viable approach for the rational design of new AB-based magnets which are beneficial for further applications in the future.

Study of sulfur (S) centered hydrogen bonding (SCHB) interactions in the literature is mostly limited to the molecular systems where S acts as a hydrogen-bond acceptor. It has been found that this unconventional SCHB is similar in strength to any conventional hydrogen bonding interaction involving electronegative atoms. However, SCHB involving S as a hydrogen-bond donor is not explored much in the literature. Herein, we have studied the nature and strength of an unconventional S–H···O hydrogen bond in a 1:1 complex of 2-fluorothiophenol (2-FTP) and H₂O using gas-phase electronic and IR spectroscopy in combination with quantum chemistry calculations. Both of the two conformers of 2-FTP···H₂O observed in the experiment are found to be stabilized primarily by S–H···O hydrogen bonding interaction. O–H···S hydrogen-bonded conformers of the complex, which are higher in energy, are not observed in the experiment. There is a nice agreement between the theoretical and experimental IR spectra of the two observed conformers. The observed IR red-shift of 25–30 cm^–1 in the S–H stretching frequency of both the conformers of the complex with respect to that of the 2-FTP monomer bespeaks that the S–H···O hydrogen bond present in 2-FTP···H₂O is weak in nature. The present work demonstrates that the S–H···O hydrogen bond can have preference over the O–H···S hydrogen bond depending on the pK_a values or proton affinities of the hydrogen bonding partners in a complex.

Heats of formation and gas phase acidities for the simple acids and their deprotonated anions (A^– = F^–, Cl^–, Br^–, I^–, OH^–, SH^–, SeH^–, TeH^–, OCl^–, OBr^–, and OI^–) were calculated using the Feller–Peterson–Dixon (FPD) method with large basis sets including Douglass–Kroll scalar relativistic corrections. Hydration of the neutral and anionic species was predicted using the supermolecule-continuum approach, resulting in absolute hydration free energies that, when combined with calculated gas phase acidities, produce aqueous acidities and pK_a values for these simple acids that are, in general, in excellent agreement with experimental literature values. Absolute hydration free energy values converged quickly in terms of the experimental values for neutral species, requiring only four explicit H₂O molecules. HI is anomalous in that it fully dissociates ionically in a water tetramer and was treated without explicit water molecules. The hydration energies of anionic species converged more slowly and were modeled with up to 16 explicit H₂O molecules. Calculated values for ΔH_f and ΔG_gas agree with experimental values within ca. 1.2 kcal/mol, and ΔG_aq and ΔΔG_hyd agree with experimental values within ca. 2 kcal/mol in most cases.

This study examined two pK_a calculation approaches (direct and proton exchange schemes) that employ high-level quantum chemical methods and implicit solvent models to predict aqueous Brønsted acidities of a large set of sulfonamides. For gas-phase deprotonation energies, the DSD-PBEP86-D3(BJ) double-hybrid functional provided the best agreement with the LNO-CCSD(T)/CBS benchmark with a mean absolute deviation less than 2 kJ mol^–1 when the aug-cc-pVTZ or larger basis sets are used. For a large test set of 54 primary and secondary sulfonamides, the use of the DSD-PBEP86-D3(BJ)/aug-cc-pVTZ level of theory in conjunction with SM12 solvation free energies predict their pK_a values with a mean accuracy of 0.9 units. In comparison, the SMD and ADF-COSMO-RS models have slightly higher mean errors of 1.4 and 1.1 pK_a units provided that the proton exchange scheme was employed to cancel the systematic errors in these models. The performance of these protocols was less ideal when applied to sulfonic acids, sulfamates, and N-substituted sulfonamides, indicating that the degree of error cancellation is sensitive to the chemical environment around the −NH₂ head group. The validated protocols were then used to estimate the pK_a values of arylsulfonamide carbonic anhydrase inhibitors, which are used to correct their experimentally measured binding free energies to account for deprotonation of the sulfonamide group upon binding to the enzyme. These corrected values did not have a significant impact on the correlation with MMGBSA binding free energies obtained from classical MD simulations where the ligand is usually considered in the deprotonated form.

In this paper, potential energy curves of Λ–S and Ω states of SBr+ are reported for the first time, and the spectrum data of some low excited bound states are obtained. The differences in the spectrum properties of main-group molecules and SBr+ were compared and analyzed, providing a sufficient theoretical basis for the subsequent study of main-group molecules. The avoided crossing that occurs in the Ω state is analyzed, and finally it is concluded that this phenomenon mainly occurs in the energy region between 20,000 and 40,000 cm–1 that is relative to the minimum energy value. Potential transitions in the Ω state capable of achieving laser cooling of SBr+ are explored. The Franck–Condon factor, radiation lifetime, and Einstein coefficient between X3Σ0+– and b1Σ0++ are calculated. From the calculation results, we concluded that direct laser cooling of SBr+ is not feasible. What we have studied in this paper provides a theoretical basis for subsequent computational exploration of the spectrum properties of SBr+.

In mixed-phase or ice clouds, ice can be formed through heterogeneous nucleation. A major type of ice-nucleating particle (INP) in the atmosphere are mineral dust particles. For mixed-phase clouds, the pH of water droplets can vary widely and influence ice nucleation by altering the surface of some INPs, including mineral dust. Kaolinite is a commonly occurring clay mineral, and laboratory experiments, as well as molecular dynamics (MD) simulations, have demonstrated its ice-nucleating efficiency at neutral pH. We examine the influence of pH on the ice-nucleating efficiency of kaolinite, in the immersion freezing mode, through both droplet freezing experiments and MD simulations. Droplet freezing experiments using KGa-1b kaolinite samples are reported under both acidic (HNO₃ solutions) and basic (NaOH solutions) conditions, covering the measured pH range 0.18–13.26. These experiments show that the ice-nucleating efficiency of kaolinite is not significantly influenced by the presence of acid but is reduced in extremely basic conditions. We report MD simulations aimed at gaining a microscopic understanding of the pH dependence of ice nucleation by kaolinite. The Al(001), Si(001), and three edge surfaces of kaolinite are considered, but ice nucleation was observed only for the Al(001) surface. The hydroxy groups exposed on the Al(001) surface can be deprotonated in a basic solution or dual-protonated in an acidic solution, which can influence ice nucleation efficiency. The protonation state of the Al(001) surface for a particular pH can be estimated using previously measured pK_a values. We find that the monoprotonated Al(001) surface expected to be stable at near-neutral pH is the most effective ice-nucleating surface. In MD simulations, the ice nucleation efficiency persists for dual-protonation but decreases significantly with increasing deprotonation, qualitatively consistent with the experimental observations. Taken together, our experimental and MD results for a wide range of pH values support the suggestion that the Al(001) surface may be important for ice nucleation by kaolinite. Additionally, the deprotonation of hydroxy groups on INP surfaces can have a significant effect on their ice-nucleating ability.

The H^– + C₂H₂ → H₂ + C₂H^– reaction is important in understanding the production mechanisms of anionic molecules in interstellar environments. Herein, the rate coefficients for the H^– + C₂H₂ → H₂ + C₂H^– reaction were calculated using ring-polymer molecular dynamics (RPMD), classical molecular dynamics (MD), and quasi-classical trajectory (QCT) approaches on a newly developed ab initio potential energy surface (PES) in full dimensions. PES was constructed by fitting a large number of ab initio energy points and their gradients using the permutationally invariant polynomial basis set method. There was no barrier in the reaction coordinates, which was a collinear-dominated reaction, and the reaction proceeded exothermically. It is found that the fitted PES provides the appropriate thermal rate coefficients based on all RPMD, classical MD, and QCT simulations at higher temperatures. The evaluation of the rate coefficients at lower temperatures should be conducted carefully because the fitting of the PES associated with the long-range interaction should be further improved. The spatial distribution of the nucleus allows a more effective attraction between the reactants.

Rotationally excited dimerization of aromatic moieties, a mechanism proposed recently to explain the initial steps of soot particle inception in combustion and pyrolysis of hydrocarbons, produces a molecular structure, termed E-bridge, combining the two aromatics via five-membered aromatic rings sharing a common bond. The present study investigates a hydrogen-mediated addition of acetylene to the fused five-membered ring part of the E-bridge forming a seven-membered ring. The carried out quantum-mechanical and rate theoretical calculations indicate the plausibility of such capping reactions, and kinetic Monte Carlo simulations demonstrate their frequent occurrence. The capping frequency, however, is limited by “splitting” the fused five-membered bridge due to five-membered ring migration. A similar migration of edge seven-membered rings is shown to be also rapid but short, as their encounter with five-membered rings converts them both into six-membered rings.

A new technique is reported to determine absolute photodissociation quantum yields, ϕ_diss, in a molecular beam. The technique relies on a molecule having two available product channels, where a species in channel A can be converted photolytically to a species in channel B. The relative decrease in the species from channel A and the relative increase in species from B provide a direct measure of the relative product yield of each channel, with no external calibration required. In the event that only channels A and B exist, or at least dominate, then the sum rule ϕ_A + ϕ_B = 1 can be used to convert relative quantum yields into absolute yields. The technique is demonstrated using the well-understood and characterized photochemistry of HCHO. Formaldehyde photolysis at wavelengths near 310 nm produces either H + HCO (channel A) or H₂ + CO (channel B). HCO can then be photolyzed with high efficiency into H + CO. The product state distributions for HCO from channel A, CO from channel B, and CO from the secondary HCO photolysis event are all well-known; this is not a requirement but is utilized here to demonstrate the veracity of the technique. The zero-pressure quantum yields of HCO from HCHO photolysis via the 2³4¹ and 2¹5¹ states of HCHO are determined to be 0.66 and 0.74, respectively, which are in excellent agreement with the established quantum yields at atmospheric pressure and support the conclusion that HCHO quantum yields at these photolysis energies are not pressure dependent.

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 impact of resonances, as well as the accuracy achievable once they are properly treated, is illustrated by the changes in spectral bandshapes, including chiroptical spectroscopies.

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 that the observed effect stems from cavity-induced changes to the Kohn–Sham orbitals of the ground state.

LModeAGen, a new protocol for the automatic determination of a nonredundant, complete set of local vibrational modes is reported, which is based on chemical graph concepts. Whereas local mode properties can be calculated for a selection of parameters targeting specific local modes of interest, a complete set of nonredundant local mode parameters is requested for the adiabatic connection scheme (ACS), relating each local vibrational mode with a normal mode counterpart, and for the decomposition of normal modes (CNM) in terms of local mode contributions, a unique way to analyze vibrational spectra. So far, nonredundant parameter sets have been generated manually following chemical intuition or from a set of redundant parameters in a trial-and-error fashion, which has hampered the study of larger systems with hundreds of parameters. LModeAGen was successfully applied for a test set of 11 systems, ranging from small molecules to the large QM (>100 atoms) subsystem of carbomonoxy-neuroglobin protein, described with a hybrid QM/MM method. The ωB97X-D/aug-cc-pVDZ, M06L/def2-TZVP, and QM/MM ωB97X-D/6-31G(d,p)/AMBER model chemistries were adopted for the description of the molecules in the test set. Our new protocol is an important step forward for a routine ACS and CNM analysis of the vibrational spectra of complex and large systems with hundreds of atoms, providing new access to important encoded electronic structure information.