Solvation of cationic copper clusters in molecular hydrogen.

Multiply charged superfluid helium nanodroplets are utilized to facilitate the growth of cationic copper clusters (Cun+, where n = 1-8) that are subsequently solvated with up to 50 H2 molecules. Production of both pristine and protonated cationic Cu clusters are detected mass spectrometrically. A joint effort between experiment and theory allows us to understand the nature of the interactions determining the bonding between pristine and protonated Cu+ and Cu2+ cations and molecular hydrogen. The analysis reveals that in all investigated cationic clusters, the primary solvation shell predominantly exhibits a covalent bonding character, which gradually decreases in strength, while for the subsequent shells an exclusive non-covalent behaviour is found. Interestingly, the calculated evaporation energies associated with the first solvation shell markedly surpass thermal values, positioning them within the desirable range for hydrogen storage applications. This comprehensive study not only provides insights into the solvation of pristine and protonated cationic Cu clusters but also sheds light on their unique bonding properties.

Aggregation of coronene: the effect of carboxyl and amine functional groups.

The aggregation of coronene is relevant to understand the formation of carbon nanomaterials, including graphene quantum dots (GQDs) that show exceptional photophysical properties. This article evaluates the influence of carboxyl and amine substituting groups on the aggregation of coronene by performing a global optimization study based on a new potential energy surface. The structures of clusters with substituted coronene are similar to those formed by un-substituted monomers, that is, stacked (non-stacked) motifs are favoured for small-size (large-size) clusters. Nonetheless, the presence of carboxyl and amine groups leads to an increase of the number of local minima of comparable energy. The clusters with substituted monomers have also shown to enhance the attractive component interaction, which can be attributed to weak induction and charge transfer effects and to stronger electrostatic contributions. Moreover, the calculated height of magic-number structures of the clusters in this work is compatible with the morphology of the GQDs reported in the literature.

Aggregation enhancement of coronene molecules by seeding with alkali-metal ions.

Microsolvation constitutes the first step in the formation of cluster structures of molecules that surround a solute in the bulk and it allows for a deep insight into the relationship between the structure of the solvation shells and other physical properties. We propose semiempirical potential energy functions that are able to describe the interaction between K+ or Cs+ with coronene. Such functions were calibrated through the comparison with accurate estimations of the interaction between the cation and the planar hydrocarbon, obtained by means of ab initio electronic-structure calculations. By employing the potential energy functions and an evolutionary algorithm (EA), we have investigated the structure and energetics of the clusters resulting from the microsolvation of either K+ or Cs+ with coronene molecules. The reliability of the results for smaller clusters was checked by performing geometry re-optimization exploiting a suitable DFT level of theory. This has allowed for the characterization of the first solvation shells of planar molecules of coronene around an alkali-metal ion. It has also been found that the presence of metal ion impurities considerably enhances the formation of small coronene clusters leading to much stronger binding energies for heterogeneous with respect to homogeneous aggregates. These clusters could represent relevant species involved in the early stages of soot nucleation.

The role of charge transfer in the stability and reactivity of chemical systems from experimental findings.

A variety of phenomena, of apparently different natures, are described within a unifying picture, by properly isolating the role of charge/electron transfer as an interaction component triggering chemical reactivity. This basic quantity is isolated by analyzing, with advanced theoretical methods developed by our group, experimental findings characterized with different techniques, such as double photo-ionization spectra, scattering cross sections and auto-ionization reaction probabilities. Suitable rationalization of such phenomena appears to be crucial for modeling the selectivity of basic elementary processes occurring in systems at increasing complexity of fundamental/applied interest, such as plasmas, flames, interstellar media, planetary atmospheres and biological environments.

Closer versus Long Range Interaction Effects on the Non-Arrhenius Behavior of Quasi-Resonant O2 + N2 Collisions.

We report in this paper an investigation on energy transfer processes from vibration to vibration and/or translation in thermal and subthermal regimes for the O2 + N2 system performed using quantum-classical calculations on different empirical, semiempirical, and ab initio potential energy surfaces. In particular, the paper focuses on the rationalization of the non-Arrhenius behavior (inversion of the temperature dependence) of the quasi-resonant vibration-to-vibration energy transfer transition rate coefficients at threshold. To better understand the microscopic nature of the involved processes, we pushed the calculations to the detail of the related cross sections and analyzed the impact of the medium and long-range components of the interaction on them. Furthermore, the variation with temperature of the dependence of the quasi-resonant rate coefficient on the vibrational energy gap between initial and final vibrational states and the effectiveness of quantum-classical calculations to overcome the limitations of the purely classical treatments were also investigated. These treatments, handled in an open molecular science fashion by chaining data and competencies of the various laboratories using a grid empowered molecular simulator, have allowed a rationalization of the dependence of the computed rate coefficients in terms of the distortion of the O2-N2 configuration during the diatom-diatom collisions. A way of relating such distortions to a smooth and continuous progress variable, allowing a proper evolution from both long to closer range formulation of the interaction and from its entrance to exit channel (through the strong interaction region) relaxed graphical representations, is also discussed in the paper.

Molecular Dynamics Simulations of Small Clusters and Liquid Hydrogen Sulfide at Different Thermodynamic Conditions.

A new force field for the intermolecular H2S-H2S interaction has been used to study the most relevant properties of the hydrogen sulfide system from gaseous to liquid phases by means of molecular dynamics (MD) simulations. In order to check the validity of the interaction formulation, ab initio CCSD(T)/aug-cc-pVTZ calculations, including the counterpoise correction on the H2S, (H2S)2, and (H2S)3 structures optimized at the MP2/aug-cc-pVDZ level, have been performed. The (H2S)2,3 systems have been characterized by performing NVE MD simulations at decreasing values of the temperature, while the liquid sulfide behavior has been investigated considering a NpT ensemble of 512 molecules at several thermodynamic states, defined by different pressure and temperature values. Additional calculations using an ensemble of 2197 molecules at two different temperatures have been performed to investigate the liquid/vapor interface of the system. The S-S, S-H, and H-H radial distribution functions and the coordination number, calculated at the same conditions used in X-ray and neutron diffraction experiments, and the evaluated thermodynamic and structural properties have been compared successfully with experimental data, thus confirming the reliability of the force field formulation and of the MD predictions.

Enhanced Flexibility of the O2 + N2 Interaction and Its Effect on Collisional Vibrational Energy Exchange.

Prompted by a comparison of measured and computed rate coefficients of Vibration-to-Vibration and Vibration-to-Translation energy transfer in O2 + N2 non-reactive collisions, extended semiclassical calculations of the related cross sections were performed to rationalize the role played by attractive and repulsive components of the interaction on two different potential energy surfaces. By exploiting the distributed concurrent scheme of the Grid Empowered Molecular Simulator we extended the computational work to quasiclassical techniques, investigated in this way more in detail the underlying microscopic mechanisms, singled out the interaction components facilitating the energy transfer, improved the formulation of the potential, and performed additional calculations that confirmed the effectiveness of the improvement introduced.

Low-energy structures of benzene clusters with a novel accurate potential surface.

The benzene-benzene (Bz-Bz) interaction is present in several chemical systems and it is known to be crucial in understanding the specificity of important biological phenomena. In this work, we propose a novel Bz-Bz analytical potential energy surface which is fine-tuned on accurate ab initio calculations in order to improve its reliability. Once the Bz-Bz interaction is modeled, an analytical function for the energy of the Bzn clusters may be obtained by summing up over all pair potentials. We apply an evolutionary algorithm (EA) to discover the lowest-energy structures of Bzn clusters (for n=2-25), and the results are compared with previous global optimization studies where different potential functions were employed. Besides the global minimum, the EA also gives the structures of other low-lying isomers ranked by the corresponding energy. Additional ab initio calculations are carried out for the low-lying isomers of Bz3 and Bz4 clusters, and the global minimum is confirmed as the most stable structure for both sizes. Finally, a detailed analysis of the low-energy isomers of the n = 13 and 19 magic-number clusters is performed. The two lowest-energy Bz13 isomers show S6 and C3 symmetry, respectively, which is compatible with the experimental results available in the literature. The Bz19 structures reported here are all non-symmetric, showing two central Bz molecules surrounded by 12 nearest-neighbor monomers in the case of the five lowest-energy structures.

Benzene-hydrogen bond (C6H6-HX) interactions: the influence of the X nature on their strength and anisotropy.

The intermolecular potential energy of the C6H6-SH2 and C6H6-NH3 dimers is formulated as combination of independent electrostatic and nonelectrostatic contributions. The relevant parameters of the nonelectrostatic terms, derived from molecular polarizability components, have been proved to be useful to describe in a consistent way both size repulsion and dispersion attraction forces. The representation adopted for the electrostatic contribution asymptotically reproduces the dipole quadrupole interaction. To test the validity of the proposed potential formulation, the features of the most stable configurations of the systems predicted have been compared with the available ab initio and experimental data. Moreover, the strength of the C6H6-HX interaction has been analyzed comparing the obtained results with the corresponding ones for the C6H6-H2O and C6H6-CH4 systems, investigated previously with the same methodology. Information on the relative orientation dependence of the partners, arising from the anisotropy of the intermolecular interaction, evaluated at different intermolecular distances, has been also obtained. Such information is crucial to evaluate sterodynamics effects in bimolecular collisions.

Chirality of weakly bound complexes: the potential energy surfaces for the hydrogen-peroxide-noble-gas interactions.

We consider the analytical representation of the potential energy surfaces of relevance for the intermolecular dynamics of weakly bound complexes of chiral molecules. In this paper we study the H2O2-Ng (Ng=He, Ne, Ar, Kr, and Xe) systems providing the radial and the angular dependence of the potential energy surface on the relative position of the Ng atom. We accomplish this by introducing an analytical representation which is able to fit the ab initio energies of these complexes in a wide range of geometries. Our analysis sheds light on the role that the enantiomeric forms and the symmetry of the H2O2 molecule play on the resulting barriers and equilibrium geometries. The proposed theoretical framework is useful to study the dynamics of the H2O2 molecule, or other systems involving O-O and S-S bonds, interacting by non-covalent forces with atoms or molecules and to understand how the relative orientation of the O-H bonds changes along collisional events that may lead to a hydrogen bond formation or even to selectivity in chemical reactions.

Carbon dioxide clathrate hydrates: selective role of intermolecular interactions and action of the SDS catalyst.

The ability of a single sodium dodecyl sulfate (SDS) molecule to promote the formation of CO2 clathrate hydrates in water (as it does for methane) has been investigated at the microscopic level. For this purpose, the components of the related force field were carefully formulated and assembled following the procedure previously adopted for methane. The properties of the whole system (as well as those of its components) were analyzed by carrying out extended molecular dynamics calculations. Contrary to what happens for methane, the calculations singled out the propensity of CO2 (pure) water clusters to form clathrate hydrate-like structures and the disappearance of such propensity when a single SDS molecule is added to the clusters. This feature was found to be due to the strong interaction of carbon dioxide with the additive that makes the SDS molecule lose its shape together with its ability to drive water molecules to form a suitable cage.

Microsolvation of the potassium ion with aromatic rings: comparison between hexafluorobenzene and benzene.

We employ a recently developed methodology to study structural and energetic properties of the first solvation shells of the potassium ion in nonpolar environments due to aromatic rings, which is important to understand the selectivity of several biochemical phenomena. Our evolutionary algorithm is used in the global optimization study of clusters formed of K(+) solvated with hexafluorobenzene (HFBz) molecules. The global intermolecular interaction for these clusters has been decomposed in HFBz-HFBz and in K(+)-HFBz contributions, using a potential model based on different decompositions of the molecular polarizability of hexafluorobenzene. Putative global minimum structures of microsolvation clusters up to 21 hexafluorobenzene molecules were obtained and compared with the analogous K(+)-benzene clusters reported in our previous work (J. Phys. Chem. A 2012, 116, 4947-4956). We have found that both K(+)-(Bz)n and K(+)-(HFBz)n clusters show a strong magic number around the closure of the first solvation shell. Nonetheless, all K(+)-benzene clusters have essentially the same first solvation shell geometry with four solvent molecules around the ion, whereas the corresponding one for K(+)-(HFBz)n is completed with nine HFBz species, and its structural motif varies as n increases. This is attributed to the ion-solvent interaction that has a larger magnitude for K(+)-Bz than in the case of K(+)-HFBz. In addition, the ability of having more HFBz than Bz molecules around K(+) in the first solvation shell is intimately related to the inversion in the sign of the quadrupole moment of the two solvent species, which leads to a distinct ion-solvent geometry of approach.

Pharmacology and toxicology of xylazine: quid novum?

The current opioid overdose crisis is characterized by the presence of unknown psychoactive adulterants. Xylazine is an alpha-2 receptor agonist that is not approved for human use but is commonly used in veterinary medicine due to its sedative and muscle-relaxant properties. Cases of human intoxication due to accidental or voluntary use have been reported since the 1980s. However, reports of adulteration of illicit opioids (heroin and illicit fentanyl) with xylazine have been increasing all over Western countries. In humans, xylazine causes respiratory depression, bradycardia, and hypotension-posing individuals, using xylazine-adulterated opioids. We present a narrative review of the latest intoxication cases related to xylazine, to bring awareness to readers and also to help pathologists to detect and deal with xylazine cases.

Competitive role of CH4-CH4 and CH-π interactions in C6H6-(CH4)n aggregates: the transition from dimer to cluster features.

The intermolecular methane-methane and benzene (Bz)-methane interactions formulated in this paper are suitable to investigate systems of increasing complexity. The proposed CH(4)-CH(4) and Bz-CH(4) potential energy functions are indeed applied to study some macroscopic properties of methane and important features of both small Bz-(CH(4))(n) (n > 1-10) clusters and Bz surrounded by several CH(4) molecules. Relevant parameters of the interaction, derived from molecular polarizability components, have been proved to be useful to describe in a consistent way both size repulsion and dispersion attraction forces. The proposed potential model also allows one to isolate the role of the different intermolecular energy contributions. The spatial distribution of the CH(4) molecules in the clusters is investigated by means of molecular dynamics simulations under various conditions, even when methane phase transition from liquid to gas is likely to occur. In addition, several properties, such as radial distribution functions, density values, and mean diffusion coefficients, are analyzed in detail.

Alkali-ion microsolvation with benzene molecules.

The target of this investigation is to characterize by a recently developed methodology, the main features of the first solvation shells of alkaline ions in nonpolar environments due to aromatic rings, which is of crucial relevance to understand the selectivity of several biochemical phenomena. We employ an evolutionary algorithm to obtain putative global minima of clusters formed with alkali-ions (M(+)) solvated with n benzene (Bz) molecules, i.e., M(+)-(Bz)(n). The global intermolecular interaction has been decomposed in Bz-Bz and in M(+)-Bz contributions, using a potential model based on different decompositions of the molecular polarizability of benzene. Specifically, we have studied the microsolvation of Na(+), K(+), and Cs(+) with benzene molecules. Microsolvation clusters up to n = 21 benzene molecules are involved in this work and the achieved global minimum structures are reported and discussed in detail. We observe that the number of benzene molecules allocated in the first solvation shell increases with the size of the cation, showing three molecules for Na(+) and four for both K(+) and Cs(+). The structure of this solvation shell keeps approximately unchanged as more benzene molecules are added to the cluster, which is independent of the ion. Particularly stable structures, so-called "magic numbers", arise for various nuclearities of the three alkali-ions. Strong "magic numbers" appear at n = 2, 3, and 4 for Na(+), K(+), and Cs(+), respectively. In addition, another set of weaker "magic numbers" (three per alkali-ion) are reported for larger nuclearities.

Angular and energy distribution of fragment ions in dissociative double photoionization of acetylene molecules at 39 eV.

The two-body dissociation reactions of the dication, C(2)H(2)(2+), produced by 39.0 eV double photoionization of acetylene molecules, have been studied by coupling photoelectron-photoion-photoion coincidence and ion imaging techniques. The results provide the kinetic energy and angular distributions of product ions. The analysis of the results indicates that the dissociation leading to C(2)H(+)+H(+) products occurs through a metastable dication with a lifetime of 108±22 ns, and a kinetic energy release (KER) distribution exhibiting a maximum at ∼4.3 eV with a full width at half maximum (FWHM) of about 60%. The reaction leading to CH(2)(+)+C(+) occurs in a time shorter than the typical rotational period of the acetylene molecules (of the order of 10(-12) s). The KER distribution of product ions for this reaction, exhibits a maximum at ∼4.5 eV with a FWHM of about 28%. The symmetric dissociation, leading to CH(+)+CH(+), exhibits a KER distribution with a maximum at ∼5.2 eV with a FWHM of 44%. For the first two reactions the angular distributions of ion products also indicate that the double photoionization of acetylene occurs when the neutral molecule is mainly oriented perpendicularly to the light polarization vector.

A portable intermolecular potential for molecular dynamics studies of NMA-NMA and NMA-H2O aggregates.

A recently formulated intermolecular potential has been adapted to describe the interaction of the N-methylacetamide (NMA) dimer and of the NMA-H(2)O adduct. The pure electrostatic component of the intermolecular potential functional representation is as usual expressed in terms of a set of punctual charges distributed over the molecular frames, consistently with the permanent molecular dipole values. In contrast, the remainder of the intermolecular potential is expressed in terms of Improved Lennard Jones effective pair potential functions, referred to multiple interaction centers (or sites) placed on the N-methylacetamide molecule and to a single interaction center placed on the water molecule. The characteristic of this pair potential relies on a mix of transferable and non-transferable descriptions of the parameters. The first set of parameters has a structural connotation bearing a site-site interaction nature and exploiting the molecular polarizability decomposability. The second one, depending on the particles clustering and charge distribution and transfer, bears a delocalized and ambient bulk nature. This choice has been tested against ab initio calculations and molecular dynamics simulations. The results show that the model potential is appropriate for describing the energetic of the various stable configurations of NMA-NMA and NMA-H(2)O weakly interacting aggregates, including the formation of hydrogen bonds.

Ar solvation shells in K(+)-HFBz: from cluster rearrangement to solvation dynamics.

The effect of some leading intermolecular interaction components on specific features of weakly bound clusters involving an aromatic molecule, a closed shell ion, and Ar atoms is analyzed by performing molecular dynamics simulations on potential energy surfaces properly formulated in a consistent way. In particular, our investigation focuses on the three-dimensional Ar distributions around the K(+)-hexafluorobenzene (K(+)-HFBz) dimer, in K(+)-HFBz-Ar(n) aggregates (n ≤ 15), and on the gradual evolution from cluster rearrangement to solvation dynamics when ensembles of 50, 100, 200, and 500 Ar atoms are taken into account. Results indicate that the Ar atoms compete to be placed in such a way to favor an attractive interaction with both K(+) and HFBz, occupying positions above and below the aromatic plane but close to the cation. When these positions are already occupied, the Ar atoms tend to be placed behind the cation, at larger distances from the center of mass of HFBz. Accordingly, three different groups of Ar atoms are observed when increasing n, with two of them surrounding K(+), thus, disrupting the K(+)-HFBz equilibrium geometry and favoring the dissociation of the solvated cation when the temperature increases. The selective role of the leading intermolecular interaction components directly depending on the ion size repulsion is discussed in detail by analyzing similarities and differences on the behavior of the Ar-solvated K(+)-HFBz and Cl(-)-Bz aggregates.

Dissociative double photoionization of singly deuterated benzene molecules in the 26-33 eV energy range.

This work provides new experimental and theoretical results about the formation and dissociation of benzene dication. The experiment has been carried out by using a vacuum ultraviolet radiation from a synchrotron source together with a time-of-flight spectrometer and a position sensitive ion detector. Isotopically labeled benzene molecules with a single deuterium atom have been used in order to study the symmetric dissociation of the benzene dication, not well evident in previous experiments. A threshold of 30.1 ± 0.1 eV has been observed for this dissociation reaction. Moreover, the lifetime of the dissociation of the benzene metastable dication producing CH(3)(+) and C(5)H(3)(+) has been obtained as a function of the photon energy, by the use of a Monte Carlo trajectory analysis of the coincidence distributions. The determined lifetime is independent of the photon energy and has an average value of 0.75 ± 0.22 µs. Theoretical calculations of the energy and structure of dissociation product ions have been also performed to provide crucial information about the dynamics of the charge separation reactions following the photoionization event.

Substance abuse in times of COVID-19 pandemic: a "perfect storm" for substance use disorder and mental health?

ABSTRACT: Our lives and communities have been besieged by COVID-19 for almost two years, and society and its functioning have been turned upside down and upset, through limitations and restrictions aimed at stemming the spread of the SARS-CoV-2 virus. Severe levels of anxiety, distress and uncertainty have taken a huge toll on our daily lives from the social, professional and emotional perspectives, but it is still rather unclear how serious an impact the emergency has had on a class of particularly vulnerable individuals: those with substance abuse issues. Particularly, we would like to draw attention to a uniquely troublesome development: forced isolation resulting from pandemic-related restrictions and how it has impacted drug users and their efforts to achieve recovery.