Sulfites detection by surface-enhanced Raman spectroscopy: A feasibility study.

The exhaustive control required for the correct wine production needs of many chemical analysis throughout the process. The most extended investigations for wine production control are focused on the quantification of total and free SO2. Most methods described in the literature have an adequate detection limit, but they usually lack reproducibility and require a previous sample treatment for the extraction of the SO2 from the wine-matrix. In this context, Surface-Enhanced Raman Spectroscopy (SERS) can be a promising technique for free SO2 determination without the need for any sample pre-processing. This work describes a proof of concept of a new methodology based on SERS and supported by Density Functional Theory (DFT) calculations to identify the active vibrational modes of the key molecules that contribute to the concentration of free SO2 in solution. Theoretical predictions and experimental outcomes are brought together to chemometrics to get a simple and real-time free SO2 monitoring. This general procedure could pave the way towards an implementation of a portable SERS detection module for in-field measurements.

Extracting Quantitative Information at Quantum Mechanical Level from Noncovalent Interaction Index Analyses.

The noncovalent interaction (NCI) index is nowadays a well-known strategy to detect NCIs in molecular systems. Even though it initially provided only qualitative descriptions, the technique has been recently extended to also extract quantitative information. To accomplish this task, integrals of powers of the electron distribution were considered, with the requirement that the overall electron density can be clearly decomposed as sum of distinct fragment contributions to enable the definition of the (noncovalent) integration region. So far, this was done by only exploiting approximate promolecular electron densities, which are given by the sum of spherically averaged atomic electron distributions and thus represent too crude approximations. Therefore, to obtain more quantum mechanically (QM) rigorous results from NCI index analyses, in this work, we propose to use electron densities obtained through the transfer of extremely localized molecular orbitals (ELMOs) or through the recently developed QM/ELMO embedding technique. Although still approximate, the electron distributions resulting from the abovementioned methods are fully QM and, above all, are again partitionable into subunit contributions, which makes them completely suitable for the NCI integral approach. Therefore, we benchmarked the integrals resulting from NCI index analyses (both those based on the promolecular densities and those based on ELMO electron distributions) against interaction energies computed at a high quantum chemical level (in particular, at the coupled cluster level). The performed test calculations have indicated that the NCI integrals based on ELMO electron densities outperform the promolecular ones. Furthermore, it was observed that the novel quantitative NCI-(QM/)ELMO approach can be also profitably exploited both to characterize and evaluate the strength of specific interactions between ligand subunits and protein residues in protein-ligand complexes and to follow the evolution of NCIs along trajectories of molecular dynamics simulations. Although further methodological improvements are still possible, the new quantitative ELMO-based technique could be already exploited in situations in which fast and reliable assessments of NCIs are crucial, such as in computational high-throughput screenings for drug discovery.

Addressing molecular optomechanical effects in nanocavity-enhanced Raman scattering beyond the single plasmonic mode.

The description of surface-enhanced Raman scattering (SERS) as a molecular optomechanical process has provided new insights into the vibrational dynamics and nonlinearities of this inelastic scattering process. In earlier studies, molecular vibrations have typically been assumed to couple with a single plasmonic mode of a metallic nanostructure, ignoring the complexity of the plasmonic response in many configurations of practical interest such as in metallic nanojunctions. By describing the plasmonic fields as a continuum, we demonstrate here the importance of considering the full plasmonic response to properly address the molecule-cavity optomechanical interaction. We apply the continuum-field model to calculate the Raman signal from a single molecule in a plasmonic nanocavity formed by a nanoparticle-on-a-mirror configuration, and compare the results of optomechanical parameters, vibrational populations, and Stokes and anti-Stokes signals of the continuum-field model with those obtained from the single-mode model. Our results reveal that high-order non-radiative plasmonic modes significantly modify the optomechanical behavior under strong laser illumination. Moreover, Raman linewidths, lineshifts, vibrational populations, and parametric instabilities are found to be sensitive to the energy of the molecular vibrational modes. The implications of adopting the continuum-field model to describe the plasmonic cavity response in molecular optomechanics are relevant in many other nanoantenna and nanocavity configurations commonly used to enhance SERS.

Steric clash in real space: biphenyl revisited.

A textbook case of twisted structure due to hydrogen-hydrogen steric clash, the biphenyl molecule, has been studied in real space from a new perspective. Long-term discrepancies regarding the origin of the steric repulsion are now reconciled under the NCI (Non Covalent Interaction) method, which reflects in 3D the balance between attractive and repulsive interactions taking place in the region between the phenyl rings. The NCI method confirms that the steric repulsion does not merely come from the H-H interaction itself, but from the many-atom interactions arising from the Cortho-H region, therefore providing rigorous physical grounds for the steric clash. This method allows a continuous scan of all the subtle changes on the electron density on going from the planar to the perpendicular biphenyl structure. The NCI results agree with other topological approaches (IQA, ELF) and are in line with previous findings in the literature regarding controversial H-H interactions in steric clash situations: H-H interactions are attractive, but repulsion appears between (Cortho-H)(Cortho-H), raising the intraatomic energy of the ortho H. ELF is also used to support these conclusions. Indeed, deformations are observed in compressed basins that allow to visualize the intraatomic effect of steric repulsion. These results can be easily extrapolated to systems with similar topological features in which steric clash is claimed to be the reason for instability.

NCIPLOT4: Fast, Robust, and Quantitative Analysis of Noncovalent Interactions.

The NonCovalent Interaction index (NCI) enables identification of attractive and repulsive noncovalent interactions from promolecular densities in a fast manner. However, the approach remained up to now qualitative, only providing visual information. We present a new version of NCIPLOT, NCIPLOT4, which allows quantifying the properties of the NCI regions (volume, charge) in small and big systems in a fast manner. Examples are provided of how this new twist enables characterization and retrieval of local information in supramolecular chemistry and biosystems at the static and dynamic levels.

Exchange Rules for Diradical π-Conjugated Hydrocarbons.

A variety of planar π-conjugated hydrocarbons such as heptauthrene, Clar's goblet and, recently synthesized, triangulene have two electrons occupying two degenerate molecular orbitals. The resulting spin of the interacting ground state is often correctly anticipated as S = 1, extending the application of Hund's rules to these systems, but this is not correct in some instances. Here we provide a set of rules to correctly predict the existence of zero mode states as well as the spin multiplicity of both the ground state and the low-lying excited states, together with their open- or closed-shell nature. This is accomplished using a combination of analytical arguments and configuration interaction calculations with a Hubbard model, both backed by quantum chemistry methods with a larger Gaussian basis set. Our results go beyond the well established Lieb's theorem and Ovchinnikov's rule, as we address the multiplicity and the open-/closed-shell nature of both ground and excited states.

Nine questions on energy decomposition analysis.

The paper collects the answers of the authors to the following questions: Is the lack of precision in the definition of many chemical concepts one of the reasons for the coexistence of many partition schemes? Does the adoption of a given partition scheme imply a set of more precise definitions of the underlying chemical concepts? How can one use the results of a partition scheme to improve the clarity of definitions of concepts? Are partition schemes subject to scientific Darwinism? If so, what is the influence of a community's sociological pressure in the "natural selection" process? To what extent does/can/should investigated systems influence the choice of a particular partition scheme? Do we need more focused chemical validation of Energy Decomposition Analysis (EDA) methodology and descriptors/terms in general? Is there any interest in developing common benchmarks and test sets for cross-validation of methods? Is it possible to contemplate a unified partition scheme (let us call it the "standard model" of partitioning), that is proper for all applications in chemistry, in the foreseeable future or even in principle? In the end, science is about experiments and the real world. Can one, therefore, use any experiment or experimental data be used to favor one partition scheme over another? © 2019 Wiley Periodicals, Inc.

Bending Carbon Nanoforms for Supramolecular Recognition: A Topological Study on Hemifullerene-Based Aggregates.

Buckybowls have risen as appealing fullerene fragment derivatives. Their intrinsic curvature has been exploited in the generation of host-guest supramolecular assemblies, not only through concave-convex complementarity but also through less-known concave-concave staggered arrangements. Whereas the stabilization of bowl-in-bowl dispositions has been ascribed to efficient π-π forces together with favorable dipole-dipole interactions, a detailed analysis on the forces guiding the formation of the staggered arrangements is missing so far. Herein, we present a thorough theoretical characterization of bowl-in-bowl vs staggered hemifullerene-based homodimers and heterodimers with the electron-donor truxTTF molecule, as test cases, under the density functional theory and by means of chemical bonding techniques. Our results clearly reveal strong and localized noncovalent signatures, together with an enhanced orbital interaction, associated with CH-π and sulfur-mediated interactions governing the staggered formation. Bending the fullerene fragment is demonstrated to favor the stabilization in both homo- and heterodimers, in good accord with the depletion in the π-electron density calculated upon increasing the buckybowl curvature. The optimal buckybowl curvature for the highest interaction energy is, however, dependent on the type of supramolecular assembly (bowl-in-bowl vs staggered) and the concave region to which hemifullerene approaches truxTTF. Interestingly, two regimes are found as a function of buckybowl curvature for hemifullerene homodimers: bowl-in-bowl dispositions are calculated more stable at low curvatures whereas staggered dimers prevail for highly curved buckybowls. Our results highlight the potential of discrete CH-π and sulfur-mediated interactions to generate unconventional staggered supramolecular arrangements toward the development of a new and unexplored host-guest chemistry.

How strong are the metallocene-metallocene interactions? Cases of ferrocene, ruthenocene, and osmocene.

An exhaustive exploration of the potential energy surfaces of ferrocene, ruthenocene and osmocene dimers has been performed. Our computations involving dispersion show that only four different isomers are present in each metallocene dimer. The collective action of small interaction energies of dispersive nature leads to a dissociation energy of 7.5 kcal mol(-1) for the ferrocene dimer. Dispersion has strong effects on the geometrical parameters, reducing the M···M distances by almost 1 Å. Our results also reveal that inclusion of entropic factors modifies the relative stability of the complexes. The nature of bonding is examined using the energy decomposition analysis and the non-covalent interaction index. Both analyses indicate that dispersion is the major contributing factor in stabilizing a metallocene dimer.

Characterizing Molecular Interactions in Chemical Systems.

Interactions between atoms have a major influence on the chemical properties of molecular systems. While covalent interactions impose the structural integrity of molecules, noncovalent interactions govern more subtle phenomena such as protein folding, bonding or self assembly. The understanding of these types of interactions is necessary for the interpretation of many biological processes and chemical design tasks. While traditionally the electron density is analyzed to interpret the quantum chemistry of a molecular system, noncovalent interactions are characterized by low electron densities and only slight variations of them--challenging their extraction and characterization. Recently, the signed electron density and the reduced gradient, two scalar fields derived from the electron density, have drawn much attention in quantum chemistry since they enable a qualitative visualization of these interactions even in complex molecular systems and experimental measurements. In this work, we present the first combinatorial algorithm for the automated extraction and characterization of covalent and noncovalent interactions in molecular systems. The proposed algorithm is based on a joint topological analysis of the signed electron density and the reduced gradient. Combining the connectivity information of the critical points of these two scalar fields enables to visualize, enumerate, classify and investigate molecular interactions in a robust manner. Experiments on a variety of molecular systems, from simple dimers to proteins or DNA, demonstrate the ability of our technique to robustly extract these interactions and to reveal their structural relations to the atoms and bonds forming the molecules. For simple systems, our analysis corroborates the observations made by the chemists while it provides new visual and quantitative insights on chemical interactions for larger molecular systems.