Blueshift of the CN stretching vibration of acetonitrile in solution: computational and experimental study.

Acetonitrile, a polar molecule that cannot form hydrogen bonds on its own, interacts with solvent molecules mainly through the lone pair of its nitrogen atom and the π electrons of its CN triple bond [Correction added on 17 July 2024, after first online publication: Acetole has been changed to Acetonitrile in the preceeding sentence.]. Interestingly, acetonitrile exhibits an unexpected strengthening of the triple bond's force constant in an aqueous environment, leading to an upshift (blueshift) in the corresponding stretching vibration: this effect contrasts with the usual consequence of hydrogen bonding on the vibrational frequencies of the acceptor groups, that is, frequency redshift. This investigation elucidates this phenomenon using Raman spectroscopy to examine the behavior of acetonitrile in organic solvent, water, and silver ion aqueous solutions, where an even more pronounced upshift is observed. Raman spectroscopy is particularly well suited for analyzing aqueous solutions due to the minimal scattering effect of water molecules across most of the vibrational spectrum. Computational approaches, both static and dynamical, based on Density Functional Theory and hybrid functionals, are employed here to interpret these findings, and accurately reproduce the vibrational frequencies of acetonitrile in different environments. Our calculations also allow an explanation for this unique behavior in terms of electric charge displacements. On the other hand, the study of the interaction of acetonitrile with water molecules and metal ions is relevant for the use of this molecule as a solvent in both chemical and pharmaceutical applications.

Accelerating NMR Shielding Calculations Through Machine Learning Methods: Application to Magnesium Sodium Silicate Glasses.

In this work, we have applied the Kernel Ridge Regression (KRR) method using a Least Square Support Vector Regression (LSSVR) approach for the prediction of the NMR isotropic magnetic shielding (σiso) of active nuclei (17O, 23Na, 25Mg, and 29Si) in a series of (Mg, Na) - silicate glasses. The Machine Learning (ML) algorithm has been trained by mapping the local environment of each atom described by the Smooth Overlap of Atomic Position (SOAP) descriptor with isotropic chemical shielding values computed with DFT using the Gauge-Included-Projector-Augmented-Wave (GIPAW) approach. The influence of different training datasets generated through molecular dynamics simulations at various temperatures and with different inter-atomic potentials has been tested and we demonstrate the importance of a wide exploration of the configurational space to enhance the transferability of the ML-regressor.  Finally, the trained ML-regressor has been used to simulate the 29Si MAS NMR spectra of systems containing up to 20000 atoms by averaging hundreds of configurations extracted from classical MD simulations to account for thermal vibrations. This ML approach is a powerful tool for the interpretation of NMR spectra using relatively large systems at a fraction of the computational time required by quantum mechanical calculations which are of high computational cost.

Impact of Atomic Defects on Ceria Surfaces on Chemical Mechanical Polishing of Silica Glass Surfaces.

This study investigates the impact of atomic defects, such as oxygen vacancies and Ce3+ ions, on cerium oxide (ceria) surfaces during chemical mechanical polishing (CMP) for silica glass finishing. Using density functional theory (DFT) and reactive molecular dynamics simulations, the interaction of orthosilicic molecules and silica glass with dry and wet ceria surfaces is explored. Defects alter the surface reactivity, leading to the dissociation of orthosilicic acid on oxygen vacancies, forming a strong Si-O-Ce bond. Hydroxylated surfaces exhibit easier oxygen vacancy formation and thermodynamically favored substitution of hydroxyl groups with orthosilicic acid. A new ReaxFF library for silica/ceria interfaces with defects is validated using DFT outcomes. Reactive MD simulations demonstrate that ceria surfaces with 30% Ce3+ ions on (111) planes exhibit higher polishing efficiency, attributed to increased Si-O-Ce bond formation. The simultaneous presence of oxygen vacancies and various acidic and basic sites on ceria surfaces enhances the polishing efficiency, involving acid-base reactions with silica. Defective surfaces show superior efficiency by removing silicate chains, contrasting with nondefective surfaces removing isolated orthosilicate units. This study provides insights into optimizing CMP processes for high-precision glass industry surface finishing.

Quantum-Chemistry Study of the Hydrolysis Reaction Profile in Borate Networks: A Benchmark.

This investigation involved an ab initio and Density Functional Theory (DFT) analysis of the hydrolysis mechanism and energetics in a borate network. The focus was on understanding how water molecules interact with and disrupt the borate network, an area where the experimental data are scarce and unreliable. The modeled system consisted of two boron atoms, bridging oxygen atoms, and varying numbers of water molecules. This setup allows for an exploration of hydrolysis under different environmental conditions, including the presence of OH- or H+ ions to simulate basic or acidic environments, respectively. Our investigation utilized both ab initio calculations at the MP2 and CCSD(T) levels and DFT with a range of exchange-correlation functionals. The findings indicate that the borate network is significantly more susceptible to hydrolysis in a basic environment, with respect to an acidic or to a neutral pH setting. The inclusion of explicit water molecules in the calculations can significantly affect the results, depending on the nature of the transition state. In fact, some transition states exhibited closed-ring configurations involving water and the boron-oxygen-boron network; in these cases, there were indeed more water molecules corresponding to lower energy barriers for the reaction, suggesting a crucial role of water in stabilizing the transition states. This study provides valuable insights into the hydrolysis process of borate networks, offering a detailed comparison between different computational approaches. The results demonstrate that the functionals B3LYP, PBE0, and wB97Xd closely approximated the reference MP2 and CCSD(T) calculated reaction pathways, both qualitatively in terms of the mechanism, and quantitatively in terms of the differences in the reaction barriers within the 0.1-0.2 eV interval for the most plausible reaction pathways. In addition, CAM-B3LYP also yielded acceptable results in all cases except for the most complicated pathway. These findings are useful for guiding further computational studies, including those employing machine learning approaches, and experimental investigations requiring accurate reference data for hydrolysis reactions in borate networks.

New Atomistic Insights on the Chemical Mechanical Polishing of Silica Glass with Ceria Nanoparticles.

Reactive molecular dynamics simulations have been used to simulate the chemical mechanical polishing (CMP) process of silica glass surfaces with the ceria (111) and (100) surfaces, which are predominantly found in ceria nanoparticles. Since it is known that an alteration layer is formed at the glass surface as a consequence of the chemical interactions with the slurry solutions used for polishing, we have created several glass surface models with different degrees of hydroxylation and porosity for investigating their morphology and chemistry after the interaction with acidic, neutral, and basic water solutions and the ceria surfaces. Both the chemical and mechanical effects under different pressure and temperature conditions have been studied and clarified. According to the simulation results, we have found that the silica slab with a higher degree of hydroxylation (thicker alteration layer) is more reactive, suggesting that proper chemical treatment is fundamental to augment the polishing efficiency. The reactivity between the silica and ceria (111) surfaces is higher at neutral pH since more OH groups present at the two surfaces increased the Si-O-Ce bonds formed at the interface. Usually, an outermost tetrahedral silicate unit connected to the rest of the silicate network through a single bond was removed during the polishing simulations. We observed that higher pressure and temperature accelerated the removal of more SiO4 units. However, excessively high pressure was found to be detrimental since the heterogeneous detachment of SiO4 units led to rougher surfaces and breakage of the Si-O-Si bond, even in the bulk of the glass. Despite the lower concentration of Ce ions at the surface resulting in the lower amount of Si-O-Ce formed, the (100) ceria surface was intrinsically more reactive than (111). The different atomic-scale mechanisms of silica removal at the two ceria surfaces were described and discussed.

Can DFT Calculations Provide Useful Information for SERS Applications?

Density functional theory (DFT) calculations allow us to reproduce the SERS (surface-enhanced Raman scattering) spectra of molecules adsorbed on nanostructured metal surfaces and extract the most information this spectroscopy is potentially able to provide. The latter point mainly concerns the anchoring mechanism and the bond strength between molecule and metal as well as the structural and electronic modifications of the adsorbed molecule. These findings are of fundamental importance for the application of this spectroscopic technique. This review presents and discusses some SERS-DFT studies carried out in Italy as a collaboration between the universities of Modena and Reggio-Emilia and of Florence, giving an overview of the information that we can extract with a combination of experimental SERS spectra and DFT modeling. In addition, a selection of the most recent studies and advancements on the DFT approach to SERS spectroscopy is reported with commentary.

Biasing crystallization in fused silica: An assessment of optimal metadynamics parameters.

Metadynamics (MetaD) is a useful technique to study rare events such as crystallization. It has been only recently applied to study nucleation and crystallization in glass-forming liquids such as silicates, but the optimal set of parameters to drive crystallization and obtain converged free energy surfaces is still unexplored. In this work, we systematically investigated the effects of the simulation conditions to efficiently study the thermodynamics and mechanism of crystallization in highly viscous systems. As a prototype system, we used fused silica, which easily crystallizes to ß-cristobalite through MetaD simulations, owing to its simple microstructure. We investigated the influence of the height, width, and bias factor used to define the biasing Gaussian potential, as well as the effects of the temperature and system size on the results. Among these parameters, the bias factor and temperature seem to be most effective in sampling the free energy landscape of melt to crystal transition and reaching convergence more quickly. We also demonstrate that the temperature rescaling from T > Tm is a reliable approach to recover free energy surfaces below Tm, provided that the temperature gap is below 600 K and the configurational space has been properly sampled. Finally, albeit a complete crystallization is hard to achieve with large simulation boxes, these can be reliably and effectively exploited to study the first stages of nucleation.

Development and Application of a ReaxFF Reactive Force Field for Cerium Oxide/Water Interfaces.

Ceria (CeO2) is a well-known catalytic oxide with many environmental, energy production, and industrial applications, most of them involving water as a reactant, byproduct, solvent, or simple spectator. In this work, we parameterized a Ce/O/H ReaxFF for the study of ceria and ceria/water interfaces. The parameters were fitted to an ab initio training set obtained at the DFT/PBE0 level, including the structures, cohesive energies, and elastic properties of the crystalline phases Ce, CeO2, and Ce2O3; the O-defective structures and energies of vacancy formation on CeO2 bulk and CeO2 (111) surface, as well as the absorption and reaction energies of H2 and H2O molecules on CeO2 (111). The new potential reproduced reasonably well all the fitted properties as well as the relative stabilities of the different ceria surfaces, the oxygen vacancies formation, and the energies and structures of associative and dissociative water molecules on them. Molecular dynamics simulations of the liquid water on the CeO2 (111) and CeO2 (100) surfaces were carried out to study the coverage and the mechanism of water dissociation. After equilibration, on average, 35% of surface sites of CeO2 (111) are hydroxylated whereas 15% of them are saturated with molecular water associatively adsorbed. As for the CeO2 (100) surface, we observed that water preferentially dissociates covering 90% of the available surface sites in excellent agreement with recent experimental findings.

Unraveling the internal conversion process within the Q-bands of a chlorophyll-like-system through surface-hopping molecular dynamics simulations.

The non-radiative relaxation process within the Q-bands of chlorophylls represents a crucial preliminary step during the photosynthetic mechanism. Despite several experimental and theoretical efforts performed in order to clarify the complex dynamics characterizing this stage, a complete understanding of this mechanism is still far to be reached. In this study, non-adiabatic excited-state molecular dynamic simulations have been performed to model the non-radiative process within the Q-bands for a model system of chlorophylls. This system has been considered in the gas phase and then, to have a more representative picture of the environment, with implicit and mixed implicit-explicit solvation models. In the first part of this analysis, absorption spectra have been simulated for each model in order to guide the setup for the non-adiabatic excited-state molecular dynamic simulations. Then, non-adiabatic excited-state molecular dynamic simulations have been performed on a large set of independent trajectories and the population of the Qx and Qy states has been computed as the average of all the trajectories, estimating the rate constant for the process. Finally, with the aim of investigating the possible role played by the solvent in the Qx-Qy crossing mechanism, an essential dynamic analysis has been performed on the generated data, allowing one to find the most important motions during the simulated dynamics.

Role of specific solute-solvent interactions on the photophysical properties of distyryl substituted BODIPY derivatives.

Computational spectroscopy is now a valuable tool to better understand experimental spectroscopic data. Recent advancements in computational methods allow the study of electronic excited states of medium-large molecular systems with greater accuracy, simulating spectra directly comparable with experiments. Besides the system of interest, proper modelling of the environment effects is essential to reach a sufficient level of accuracy. This is especially important for solvents with high polarity or those able to establish specific interactions with the solute. One of the most straightforward and fastest ways to deal with solvent effects is the use of implicit models, like the polarizable continuum model (PCM). However, these models cannot describe specific solute-solvent interactions. In this case, one possible solution is the adoption of a mixed model that includes the solute and few solvent molecules within a PCM cavity. With the aim of investigating the role played by specific solute-solvent interactions, the behaviour of styryl substituted BODIPY systems in methanol has been studied at a theoretical level, considering both implicit and mixed implicit-explicit models. In the first part of our analysis, vibrationally resolved electronic spectra for absorption and emission processes have been simulated, using several sets of coordinates to represent normal modes, and compared with experiments. Then, to verify if specific solute-solvent interactions play a crucial role during the ultrafast intraband relaxation processes, 2D electronic spectra were also simulated, which could provide valuable information on the complex ultrafast dynamics.

Disclosing the Interaction of Gold Nanoparticles with Aß(1-40) Monomers through Replica Exchange Molecular Dynamics Simulations.

Amyloid-ß aggregation is one of the principal causes of amyloidogenic diseases that lead to the loss of neuronal cells and to cognitive impairments. The use of gold nanoparticles treating amyloidogenic diseases is a promising approach, because the chemistry of the gold surface can be tuned in order to have a specific binding, obtaining effective tools to control the aggregation. In this paper, we show, by means of Replica Exchange Solute Tempering Molecular Simulations, how electrostatic interactions drive the absorption of Amyloid-ß monomers onto citrates-capped gold nanoparticles. Importantly, upon binding, amyloid monomers show a reduced propensity in forming ß-sheets secondary structures that are characteristics of mature amyloid fibrils.

Insights into the Effect of Curcumin and (-)-Epigallocatechin-3-Gallate on the Aggregation of Aß(1-40) Monomers by Means of Molecular Dynamics.

In this study, we compared the effects of two well-known natural compounds on the early step of the fibrillation process of amyloid-ß (1-40), responsible for the formation of plaques in the brains of patients affected by Alzheimer's disease (AD). The use of extensive replica exchange simulations up to the µs scale allowed us to characterize the inhibition activity of (-)-epigallocatechin-3-gallate (EGCG) and curcumin (CUR) on unfolded amyloid fibrils. A reduced number of ß-strands, characteristic of amyloid fibrils, and an increased distance between the amino acids that are responsible for the intra- and interprotein aggregations are observed. The central core region of the amyloid-ß (Aß(1-40)) fibril is found to have a high affinity to EGCG and CUR due to the presence of hydrophobic residues. Lastly, the free binding energy computed using the Poisson Boltzmann Surface Ares suggests that EGCG is more likely to bind to unfolded Aß(1-40) fibrils and that this molecule can be a good candidate to develop new and more effective congeners to treat AD.

Structure Model and Toxicity of the Product of Biodissolution of Chrysotile Asbestos in the Lungs.

Asbestos is a commercial term indicating six natural silicates with asbestiform crystal habit. Of these, five are double-chain silicates (amphibole) and one is a layer silicate (serpentine asbestos or chrysotile). Although all species are classified as human carcinogens, their degree of toxicity is still a matter of debate. Amphibole asbestos species are biopersistent in the human lungs and exert their chronic toxic action for decades, whereas chrysotile is not biopersistent and transforms into an amorphous silica structure prone to chemical/physical clearance when exposed to the acidic environment created by the alveolar macrophages. There is evidence in the literature of the toxicity of chrysotile, but its limited biopersistence is thought to explain the difference in toxicity with respect to amphibole asbestos. To date, no comprehensive model describing the toxic action of chrysotile in the lungs is available, as the structure and toxic action of the product formed by the biodissolution of chrysotile are unknown. This work is aimed at fulfilling this gap and explaining the toxic action in terms of structural, chemical, and physical properties. We show that chrysotile's fibrous structure induces cellular damage, mainly through physical interactions. Based on our previous work and novel findings, we propose the following toxicity model: inhaled chrysotile fibers exert their toxicity in the alveolar space by physical and biochemical action. The fibers are soon leached by the intracellular acid environment into a product with residual toxicity, and the dissolution process liberates toxic metals in the intracellular and extracellular environment.

On the simulation of vibrationally resolved electronic spectra of medium-size molecules: the case of styryl substituted BODIPYs.

BODIPY dyes are used in a variety of applications because of their peculiar spectroscopic and photo-physical properties that vary depending on the stereochemistry of the functional groups attached to the boron-dipyrromethene core structure. In this work, we have applied several computational methods, adapted for semi-rigid molecules based on the Franck-Condon principle, for the study of the optical properties of BODIPY systems and for the understanding of the influence of functional groups on their spectroscopic features. We have analyzed the electronic spectra of two styryl substituted BODIPY molecules of technological interest, properly taking into account the vibronic contribution. For comparison with recently recorded experimental data in methanol, the vibrationally resolved electronic spectra of these systems were computed using both Time-Independent (TI) and Time-Dependent (TD) formalisms. The first step toward the analysis of optical properties of the styryl modified BODIPYs was a benchmark of several density functionals, to select the most appropriate one. We have found that all benchmarked functionals provide good results in terms of band shape but some of them show strong discrepancies in terms of band position. Beyond the issue of the electronic structure calculation method, different levels of sophistication can be adopted for the calculation of vibronic transitions. In this study, the effect of mode couplings and the influence of the Herzberg-Teller terms on the theoretical spectra has been investigated. It has been found that all levels of theory considered give reproducible results for the investigated systems: band positions and shapes are similar at all levels and little improvements have been found in terms of band shape with the inclusion of Herzberg-Teller effect. Inclusion of temperature effects proved to be challenging due to the important impact of large amplitude motions. Better agreement can be achieved by adopting a suitable set of coordinates coupled with a reduced-dimensionality scheme.

SERS, XPS and DFT Study of Xanthine Adsorbed on Citrate-Stabilized Gold Nanoparticles.

We have studied the adsorption of xanthine, a nucleobase present in human tissue and fluids that is involved in important metabolic processes, on citrate-reduced gold colloidal nanoparticles by means of surface-enhanced Raman scattering (SERS), absorption, and X-ray photoelectron spectroscopy (XPS) measurements, along with density functional theory (DFT) calculations. The citrate anions stabilize the colloidal suspensions by strongly binding the gold nanoparticles. However, these anions do not impair the adsorption of xanthine on positively-charged active sites present on the metal surface. We have obtained the Fourier transform (FT)-SERS spectra of adsorbed xanthine by laser excitation in the near infrared spectral region, where interference due to fluorescence emission does not usually occur. In fact, the addition of chloride ions to the Au/xanthine colloid induces the aggregation of the gold nanoparticles, whose plasmonic band is shifted to the near infrared region where there is the exciting laser line of the FT-Raman instrument. Hence, this analytical approach is potentially suitable for spectroscopic determination of xanthine directly in body fluids, avoiding fluorescence phenomena induced by visible laser irradiation.

Multiscale Molecular Dynamics Simulation of Multiple Protein Adsorption on Gold Nanoparticles.

A multiscale molecular dynamics simulation study has been carried out in order to provide in-depth information on the adsorption of hemoglobin, myoglobin, and trypsin over citrate-capped AuNPs of 15 nm diameter. In particular, determinants for single proteins adsorption and simultaneous adsorption of the three types of proteins considered have been studied by Coarse-Grained and Meso-Scale molecular simulations, respectively. The results, discussed in the light of the controversial experimental data reported in the current experimental literature, have provided a detailed description of the (i) recognition process, (ii) number of proteins involved in the early stages of corona formation, (iii) protein competition for AuNP adsorption, (iv) interaction modalities between AuNP and protein binding sites, and (v) protein structural preservation and alteration.

Computational Insight into the Effect of Natural Compounds on the Destabilization of Preformed Amyloid-ß(1⁻40) Fibrils.

One of the principal hallmarks of Alzheimer's disease (AD) is related to the aggregation of amyloid-ß fibrils in an insoluble form in the brain, also known as amyloidosis. Therefore, a prominent therapeutic strategy against AD consists of either blocking the amyloid aggregation and/or destroying the already formed aggregates. Natural products have shown significant therapeutic potential as amyloid inhibitors from in vitro studies as well as in vivo animal tests. In this study, the interaction of five natural biophenols (curcumin, dopamine, (-)-epigallocatechin-3-gallate, quercetin, and rosmarinic acid) with amyloid-ß(1⁻40) fibrils has been studied through computational simulations. The results allowed the identification and characterization of the different binding modalities of each compounds and their consequences on fibril dynamics and aggregation. It emerges that the lateral aggregation of the fibrils is strongly influenced by the intercalation of the ligands, which modulates the double-layered structure stability.

Modeling emission features of salicylidene aniline molecular crystals: A QM/QM' approach.

A new computational protocol relying on the use of electrostatic embedding, derived from QM/QM' ONIOM calculations, to simulate the effect of the crystalline environment on the emission spectra of molecular crystals is here applied to the ß-form of salicylidene aniline (SA). The first singlet excited states (S1 ) of the SA cis-keto and trans-keto conformers, surrounded by a cluster of other molecules representing the crystalline structure, were optimized by using a QM/QM' ONIOM approach with and without electronic embedding. The model system consisting of the central salicylidene aniline molecule was treated at the DFT level by using either the B3LYP, PBE0, or the CAM-B3LYP functional, whereas the real system was treated at the HF level. The CAM-B3LYP/HF level of theory provides emission energies in good agreement with experiment with differences of -20/-32 nm (cis-keto form) and -8/-14 nm (trans-keto form), respectively, whereas notably larger differences are obtained using global hybrids. Though such differences on the optical properties arise from the density functional choice, the contribution of the electronic embedding is rather independent of the functional used. This plays in favor of a more general applicability of the present protocol to other crystalline molecular systems.

Raman and DFT study of methimazole chemisorbed on gold colloidal nanoparticles.

The adsorption of methimazole on gold colloidal nanoparticles was investigated using a combination of surface-enhanced Raman scattering and density functional theory calculations, which allowed identifying the thiolate anion as the molecular species chemically interacting with the active sites of the gold surface, modeled as zero-charge metal adatoms, only through the sulfur atom. This result can be important for the use of these ligand/metal nanohybrids in the process of drug delivery. Moreover, functionalized gold nanoparticles are able to promote the Raman enhancement in the red-light region as well as in the near-infrared, where generally no fluorescence emission occurs. This paves the way for the use of these nanosystems in a biological environment, even in vivo experiments.

Insights into structural and dynamical features of water at halloysite interfaces probed by DFT and classical molecular dynamics simulations.

Density functional theory calculations and classical molecular dynamics simulations have been used to investigate the structure and dynamics of water molecules on kaolinite surfaces and confined in the interlayer of a halloysite model of nanometric dimension. The first technique allowed us to accurately describe the structure of the tetrahedral-octahedral slab of kaolinite in vacuum and in interaction with water molecules and to assess the performance of two widely employed empirical force fields to model water/clay interfaces. Classical molecular dynamics simulations were used to study the hydrogen bond network structure and dynamics of water adsorbed on kaolinite surfaces and confined in the halloysite interlayer. The results are in nice agreement with the few experimental data available in the literature, showing a pronounced ordering and reduced mobility of water molecules at the hydrophilic octahedral surfaces of kaolinite and confined in the halloysite interlayer, with respect to water interacting with the hydrophobic tetrahedral surfaces and in the bulk. Finally, this investigation provides new atomistic insights into the structural and dynamical properties of water-clay interfaces, which are of fundamental importance for both natural processes and industrial applications.