Binding energies of ethanol and ethylamine on interstellar water ices: synergy between theory and experiments.

Experimental and computational chemistry are two disciplines used to conduct research in astrochemistry, providing essential reference data for both astronomical observations and modeling. These approaches not only mutually support each other, but also serve as complementary tools to overcome their respective limitations. Leveraging on such synergy, we characterized the binding energies (BEs) of ethanol (CH3CH2OH) and ethylamine (CH3CH2NH2), two interstellar complex organic molecules (iCOMs), on crystalline and amorphous water ices through density functional theory (DFT) calculations and temperature-programmed desorption (TPD) experiments. Experimentally, CH3CH2OH and CH3CH2NH2 behave similarly, in which desorption temperatures are higher on the water ices than on a bare gold surface. Computed cohesive energies of pure ethanol and ethylamine bulk structures allow describing of the BEs of the pure species deposited on the gold surface, as extracted from the TPD curve analyses. The BEs of submonolayer coverages of CH3CH2OH and CH3CH2NH2 on the water ices cannot be directly extracted from TPD due to their co-desorption with water, but they are computed through DFT calculations, and found to be greater than the cohesive energy of water. The behaviour of CH3CH2OH and CH3CH2NH2 is different when depositing adsorbate multilayers on the amorphous ice, in that, according to their computed cohesive energies, ethylamine layers present weaker interactions compared to ethanol and water. Finally, from the computed BEs of ethanol, ethylamine and water, we can infer that the snow-lines of these three species in protoplanetary disks will be situated at different distances from the central star. It appears that a fraction of ethanol and ethylamine is already frozen on the grains in the water snow-lines, causing their incorporation in water-rich planetesimals.

Adsorption of HCN on cosmic silicates: a periodic quantum mechanical study.

Hydrogen cyanide (HCN) represents a small but widely distributed fraction of the interstellar molecules, and it has been observed in all the environments characterizing the formation of a new planetary system. HCN can polymerize to form biomolecules, including adenine (H5C5N5), and it has drawn attention as a possible precursor of several building blocks of life due to the presence of its polymerization products in meteorites, comets and other asteroidal bodies. To elucidate the potential catalytic role that cosmic silicates have played in these processes, we have investigated, at DFT-PBE level inclusive of a posteriori dispersion correction, the energetic and spectroscopic features of the adsorption of HCN molecules on the most relevant crystalline surfaces of the mineral forsterite (Mg2SiO4), a common silicate constituent of the interstellar core grains and planetary rocky bodies. The results reveal that HCN adsorbs both in molecular and dissociative ways, within a wide range of adsorption energies (-29.4 to -466.4 kJ mol-1). Thermodynamic and kinetic results show that dissociative adsorption is dominant already at low temperatures, a fact particularly relevant at the protoplanetary conditions (i.e., the latest stages in the star system formation process). The simulated spectroscopic features of the studied adducts show a wide range of different degrees of perturbation of C-H and CN bonds. This finding agrees with previous experimental works, and our results confirm that a complex chemistry is observed when this astrochemically-relevant molecule interacts with Mg2SiO4, which may be associated with a considerable potential reactivity towards the formation of relevant prebiotic compounds.

Nearly free surface silanols are the critical molecular moieties that initiate the toxicity of silica particles.

Inhalation of silica particles can induce inflammatory lung reactions that lead to silicosis and/or lung cancer when the particles are biopersistent. This toxic activity of silica dusts is extremely variable depending on their source and preparation methods. The exact molecular moiety that explains and predicts this variable toxicity of silica remains elusive. Here, we have identified a unique subfamily of silanols as the major determinant of silica particle toxicity. This population of "nearly free silanols" (NFS) appears on the surface of quartz particles upon fracture and can be modulated by thermal treatments. Density functional theory calculations indicates that NFS locate at an intersilanol distance of 4.00 to 6.00 Å and form weak mutual interactions. Thus, NFS could act as an energetically favorable moiety at the surface of silica for establishing interactions with cell membrane components to initiate toxicity. With ad hoc prepared model quartz particles enriched or depleted in NFS, we demonstrate that NFS drive toxicity, including membranolysis, in vitro proinflammatory activity, and lung inflammation. The toxic activity of NFS is confirmed with pyrogenic and vitreous amorphous silica particles, and industrial quartz samples with noncontrolled surfaces. Our results identify the missing key molecular moieties of the silica surface that initiate interactions with cell membranes, leading to pathological outcomes. NFS may explain other important interfacial processes involving silica particles.

Physico-Chemical Approaches to Investigate Surface Hydroxyls as Determinants of Molecular Initiating Events in Oxide Particle Toxicity.

The study of molecular recognition patterns is crucial for understanding the interactions between inorganic (nano)particles and biomolecules. In this review we focus on hydroxyls (OH) exposed at the surface of oxide particles (OxPs) which can play a key role in molecular initiating events leading to OxPs toxicity. We discuss here the main analytical methods available to characterize surface OH from a quantitative and qualitative point of view, covering thermogravimetry, titration, ζ potential measurements, and spectroscopic approaches (NMR, XPS). The importance of modelling techniques (MD, DFT) for an atomistic description of the interactions between membranes/proteins and OxPs surfaces is also discussed. From this background, we distilled a new approach methodology (NAM) based on the combination of IR spectroscopy and bioanalytical assays to investigate the molecular interactions of OxPs with biomolecules and membranes. This NAM has been already successfully applied to SiO2 particles to identify the OH patterns responsible for the OxPs' toxicity and can be conceivably extended to other surface-hydroxylated oxides.

On the stability of peptide secondary structures on the TiO2 (101) anatase surface: a computational insight.

The biological activity of proteins is partly due to their secondary structures and conformational states. Peptide chains are rather flexible so that finding ways inducing protein folding in a well-defined state is of great importance. Among the different constraint techniques, the interaction of proteins with inorganic surfaces is a fruitful strategy to stabilize selected folded states. Surface-induced peptide folding can have potential applications in different biomedicine areas, but it can also be of fundamental interest in prebiotic chemistry since the biological activity of a peptide can turn-on when folded in a given state. In this work, periodic quantum mechanical simulations (including implicit solvation effects) at the PBE-D2* level have been carried out to study the adsorption and the stability of the secondary structures (α-helix and ß-sheet) of polypeptides with different chemical composition (i.e., polyglycine, polyalanine, polyglutamic acid, polylysine, and polyarginine) on the TiO2 (101) anatase surface. The computational cost is reduced by applying periodic boundary conditions to both the surface and the peptides, thus obtaining full periodic polypeptide/TiO2 surface systems. At variance with polyglycine, the interaction of the other polypeptides with the surface takes place with the lateral chain functionalities, leaving the secondary structures almost undistorted. Results indicate that the preferred conformation upon adsorption is the α-helix over the ß-sheet, with the exception of the polyglutamic acid. According to the calculated adsorption energies, the affinity trend of the polypeptides with the (101) anatase surface is: polyarginine ≈ polylysine > polyglutamic acid > polyglycine ≈ polyalanine, both when adsorbed in gas phase and in presence of the implicit water solvent, which is very similar to the trend for the single amino acids. A set of implications related to the areas of surface-induced peptide folding, biomedicine and prebiotic chemistry are finally discussed.

From gaseous HCN to nucleobases at the cosmic silicate dust surface: an experimental insight into the onset of prebiotic chemistry in space.

HCN in the gas form is considered as a primary nitrogen source for the synthesis of prebiotic molecules in extraterrestrial environments. Nevertheless, the research mainly focused on the reactivity of HCN and its derivatives in aqueous systems, often using external high-energy supply in the form of cosmic rays or high energy photons. Very few studies have been devoted to the chemistry of HCN in the gas phase or at the gas/solid interphase, although they represent the more common scenarios in the outer space. In this paper we report about the reactivity of highly pure HCN in the 150-300 K range at the surface of amorphous and crystalline Mg2SiO4 (forsterite olivine), i.e. of solids among the constituents of the core of cosmic dust particles, comets, and meteorites. Amorphous silica and MgO were also studied as model representatives of Mg2SiO4 structural building blocks. IR spectroscopic results and the HR-MS analysis of the reaction products revealed Mg2+O2- acid/base pairs at the surface of Mg2SiO4 and MgO to be key in promoting the formation of HCN oligomers along with imidazole and purine compounds, already under very mild temperature and HCN pressure conditions, i.e. in the absence of external energetic triggers. Products include adenine nucleobase, a result which supports the hypothesis that prebiotic molecular building blocks can be easily formed through surface catalytic processes in the absence of high-energy supply.

Tracing the Primordial Chemical Life of Glycine: A Review from Quantum Chemical Simulations.

Glycine (Gly), NH2CH2COOH, is the simplest amino acid. Although it has not been directly detected in the interstellar gas-phase medium, it has been identified in comets and meteorites, and its synthesis in these environments has been simulated in terrestrial laboratory experiments. Likewise, condensation of Gly to form peptides in scenarios resembling those present in a primordial Earth has been demonstrated experimentally. Thus, Gly is a paradigmatic system for biomolecular building blocks to investigate how they can be synthesized in astrophysical environments, transported and delivered by fragments of asteroids (meteorites, once they land on Earth) and comets (interplanetary dust particles that land on Earth) to the primitive Earth, and there react to form biopolymers as a step towards the emergence of life. Quantum chemical investigations addressing these Gly-related events have been performed, providing fundamental atomic-scale information and quantitative energetic data. However, they are spread in the literature and difficult to harmonize in a consistent way due to different computational chemistry methodologies and model systems. This review aims to collect the work done so far to characterize, at a quantum mechanical level, the chemical life of Gly, i.e., from its synthesis in the interstellar medium up to its polymerization on Earth.

First-Principles Modeling of Protein/Surface Interactions. Polyglycine Secondary Structure Adsorption on the TiO2 (101) Anatase Surface Adopting a Full Periodic Approach.

Computational modeling of protein/surface systems is challenging since the conformational variations of the protein and its interactions with the surface need to be considered at once. Adoption of first-principles methods to this purpose is overwhelming and computationally extremely expensive so that, in many cases, dramatically simplified systems (e.g., small peptides or amino acids) are used at the expenses of modeling nonrealistic systems. In this work, we propose a cost-effective strategy for the modeling of peptide/surface interactions at a full quantum mechanical level, taking the adsorption of polyglycine on the TiO2 (101) anatase surface as a test case. Our approach is based on applying the periodic boundary conditions for both the surface model and the polyglycine peptide, giving rise to full periodic polyglycine/TiO2 surface systems. By proceeding this way, the considered complexes are modeled with a drastically reduced number of atoms compared with the finite-analogous systems, modeling the polypeptide structures at the same time in a realistic way. Within our modeling approach, full periodic density functional theory calculations (including implicit solvation effects) and ab initio molecular dynamics (AIMD) simulations at the PBE-D2* theory level have been carried out to investigate the adsorption and relative stability of the different polyglycine structures (i.e., extended primary, ß-sheet, and α-helix) on the TiO2 surface. It has been found that, upon adsorption, secondary structures become partially denatured because the peptide C═O groups form Ti-O═C dative bonds. AIMD simulations have been fundamental to identify these phenomena because thermal and entropic effects are of paramount importance. Irrespective of the simulated environments (gas phase and implicit solvent), adsorption of the α-helix is more favorable than that of the ß-sheet because in the former, more Ti-O═C bonds are formed and the adsorbed secondary structure results less distorted with respect to the isolated state. Under the implicit water solvent, additionally, adsorbed ß-sheet structures weaken with respect to their isolated states as the H-bonds between the strands are longer due to solvation effects. Accordingly, the results indicate that the preferred conformation upon adsorption is the α-helix over the ß-sheet.

Infrared harmonic features of collagen models at B3LYP-D3: From amide bands to the THz region.

In this paper, we have studied the vibrational spectral features for the collagen triple helix using a dispersion corrected hybrid density functional theory (DFT-D) approach. The protein is simulated by an infinite extended polymer both in the gas phase and in a water micro-solvated environment. We have adopted proline-rich collagen models in line with the high content of proline in natural collagens. Our scaled harmonic vibrational spectra are in very good agreement with the experiments and allow for the peak assignment of the collagen amide I and III bands, supporting or questioning the experimental interpretation by means of vibrational normal modes analysis. Furthermore, we demonstrated that IR spectroscopy in the THz region can detect the small variations inherent to the triple helix helicity (10/3 over 7/2), thus elucidating the packing state of the collagen. So far, identifying the collagen helicity is only possible by means of crystal x-ray diffraction.

A quantum mechanical study of dehydration vs. decarbonylation of formamide catalysed by amorphous silica surfaces.

Formamide is abundant in the interstellar medium and was also present during the formation of the Solar system through the accretion process of interstellar dust. Under the physicochemical conditions of primordial Earth, formamide could have undergone decomposition, either via dehydration (HCN + H2O) or via decarbonylation (CO + NH3). The first reactive channel provides HCN, which is an essential molecular building block for the formation of RNA/DNA bases, crucial for the emergence of life on Earth. In this work, we studied, at the CCSD(T)/cc-pVTZ level, the two competitive routes of formamide decomposition, i.e. dehydration and decarbonylation, either in liquid formamide (by using the polarization continuum model technique) or at the interface between liquid formamide and amorphous silica. Amorphous silica was adopted as a convenient model of the crystalline silica phases ubiquitously present in the primordial (and actual) Earth's crust, and also due to its relevance in catalysis, adsorption and chromatography. Results show that: (i) silica surface sites catalyse both decomposition channels by reducing the activation barriers by about 100 kJ mol-1 with respect to the reactions in homogeneous medium, and (ii) the dehydration channel, giving rise to HCN, is strongly favoured from a kinetic standpoint over decarbonylation, the latter being, instead, slightly favoured from a thermodynamic point of view.

Monitoring the Reactivity of Formamide on Amorphous SiO2 by In-Situ UV-Raman Spectroscopy and DFT Modeling.

Formamide has been recognized in the literature as a key species in the formation of the complex molecules of life, such as nucleobases. Furthermore, several studies reported the impact of mineral phases as catalysts for its decomposition/polymerization processes, increasing the conversion and also favoring the formation of specific products. Despite the progresses in the field, in situ studies on these mineral-catalyzed processes are missing. In this work, we present an in situ UV-Raman characterization of the chemical evolution of formamide over amorphous SiO2 samples, selected as a prototype of silicate minerals. The experiments were carried out after reaction of formamide at 160 °C on amorphous SiO2 (Aerosil OX50) either pristine or pre-calcined at 450 °C, to remove a large fraction of surface silanol groups. Our measurements, interpreted on the basis of density functional B3LYP-D3 calculations, allow to assign the spectra bands in terms of specific complex organic molecules, namely, diaminomaleonitrile (DAMN), 5-aminoimidazole (AI), and purine, showing the role of the mineral surface on the formation of relevant prebiotic molecules.

Carbon monoxide adsorption at forsterite surfaces as models of interstellar dust grains: An unexpected bathochromic (red) shift of the CO stretching frequency.

Carbon monoxide (CO) is one of the most abundant species in the interstellar medium (ISM). In the colder regions of the ISM, it can directly adsorb onto exposed Mg cations of forsterite (Fo, Mg2SiO4), one of the main constituents of the dust grains. Its energetic of adsorption can strongly influence the chemico-physical evolution of cold interstellar clouds; thus, a detailed description of this process is desirable. We recently simulated the CO adsorption on crystalline Fo surfaces by computer ab initio methods and, surprisingly, reported cases where the CO stretching frequency underwent a bathochromic (red) shift (i.e., it is lowered with respect to the CO gas phase frequency), usually not experimentally observed for CO adsorbed onto oxides with non-d cations, like the present case. Here, we elucidate in deep when and under which conditions this case may happen and concluded that this red shift may be related to peculiar surface sites occurring at the morphologically complex Fo surfaces. The reasons for the red shift are linked to both the quadrupolar nature of the CO molecule and the role of dispersion interactions with surfaces of complex morphology. The present work, albeit speculative, suggests that, at variance with CO adsorption on simple oxides like MgO, the CO spectrum may exhibit features at lower frequencies than the reference gas frequency when CO is adsorbed on complex oxides, even in the absence of transition metal ions.

When the Surface Matters: Prebiotic Peptide-Bond Formation on the TiO2 (101) Anatase Surface through Periodic DFT-D2 Simulations.

The mechanism of the peptide-bond formation between two glycine (Gly) molecules has been investigated by means of PBE-D2* and PBE0-D2* periodic simulations on the TiO2 (101) anatase surface. This is a process of great relevance both in fundamental prebiotic chemistry, as the reaction univocally belongs to one of the different organizational events that ultimately led to the emergence of life on Earth, as well as from an industrial perspective, since formation of amides is a key reaction for pharmaceutical companies. The efficiency of the surface catalytic sites is demonstrated by comparing the reactions in the gas phase and on the surface. At variance with the uncatalyzed gas-phase reaction, which involves a concerted nucleophilic attack and dehydration step, on the surface these two steps occur along a stepwise mechanism. The presence of surface Lewis and Brönsted sites exerts some catalytic effect by lowering the free energy barrier for the peptide-bond formation by about 6 kcal mol-1 compared to the gas-phase reaction. Moreover, the co-presence of molecules acting as proton-transfer assistants (i.e., H2 O and Gly) provide a more significant kinetic energy barrier decrease. The reaction on the surface is also favorable from a thermodynamic standpoint, involving very large and negative reaction energies. This is due to the fact that the anatase surface also acts as a dehydration agent during the condensation reaction, since the outermost coordinatively unsaturated Ti atoms strongly anchor the released water molecules. Our theoretical results provide a comprehensive atomistic interpretation of the experimental results of Martra et al. (Angew. Chem. Int. Ed. 2014, 53, 4671), in which polyglycine formation was obtained by successive feedings of Gly vapor on TiO2 surfaces in dry conditions and are, therefore, relevant in a prebiotic context envisaging dry and wet cycles occurring, at mineral surfaces, in a small pool.

Working at the membrane interface: Ligand-induced changes in dynamic conformation and oligomeric structure in human aromatase.

Aromatase catalyzes the biosynthesis of estrogens from androgens. Owing to the physiological importance of this conversion of lipophilic substrates, the interaction with the lipid bilayer for this cytochrome P450 is crucial for its dynamics that must allow an easy access to substrates and inhibitors. Here, the aromatase-anastrozole interaction is studied by combining computational methods to identify possible access/egress routes with the protein inserted in the membrane and experimental tools aimed at the investigation of the effect of the inhibitor on the protein conformation. By means of molecular dynamics simulations of the protein inserted in the membrane, two channels, not detected in the starting crystal structure, are found after a 20-nSec simulation. Trypsin digestion on the recombinant protein shows that the enzyme is strongly protected by the presence of the substrate and even more by the inhibitor. DSC experiments show an increase in the melting temperature of the protein in complex with the substrate (49.3 °C) and the inhibitor (58.7 °C) compared to the ligand-free enzyme (45.9 °C), consistent with a decrease of flexibility of the protein. The inhibitor anastrozole enters the active site of the protein through a channel different from that used from the substrate and promotes a conformational change that stiffens the protein conformation and decreases the protein-protein interaction between different aromatase molecules.

Models for biomedical interfaces: a computational study of quinone-functionalized amorphous silica surface features.

A density functional theory (PBE functional) investigation is carried out, in which a model of an amorphous silica surface is functionalized by ortho-benzoquinone. Surface functionalization with catechol and quinone-based compounds is relevant in biomedical fields, from prosthetic implants to dentistry, to develop multifunctional coatings with antimicrobial properties. The present study provides atomistic information on the specific interactions between the functionalizing agent and the silanol groups at the silica surface. The distinct configurations of the functional groups, the hydrogen bond pattern, the role of dispersion forces and the simulated IR spectra provide detailed insight into the features of this model surface coating. Ab initio molecular dynamics gives further insights into the mobility of the functionalizing groups. As a final step, we studied the condensation reaction with allylamine, via Schiff base formation, to ground subsequent simulations on condensation with model peptides of antimicrobial activity.

Modeling hydroxylated nanosilica: Testing the performance of ReaxFF and FFSiOH force fields.

We analyze the performance of the FFSiOH force field and two parameterisations of the ReaxFF force field for modeling hydroxylated nanoscale silica (SiO2). Such nanosystems are fundamental in numerous aspects of geochemistry and astrochemistry and also play a key role during the hydrothermal synthesis of technologically important nanoporous silicas (e.g., catalysts, absorbents, and coatings). We consider four aspects: structure, relative energies, vibrational spectra, and hydroxylation energies, and compare the results with those from density functional calculations employing a newly defined dataset (HND: Hydroxylated Nanosilica Dataset). The HND consists of three sets of (SiO2)16(H2O)N nanoparticles (NPs), each with a different degree of hydroxylation and each containing between 23 and 26 distinct isomers and conformers. We also make all HND reference data openly available. We further consider hydroxylated silica NPs of composition (SiO2)M(H2O)N with M = 4, 8, 16, and 24 and infinite surface slabs of amorphous silica, both with variable hydroxylation. For energetics, both ReaxFF and FFSiOH perform well for NPs with an intermediate degree of hydroxylation. For increased hydroxylation, the performance of FFSiOH begins to significantly decline. Conversely, for the lower degree of hydroxylation both parameterisations of ReaxFF do not perform well. For vibrational frequencies, FFSiOH performs particularly well and significantly better than ReaxFF. This feature also opens the door to inexpensively calculating Gibbs free energies of the hydroxylated nanosilica systems in order to efficiently correct density functional theory calculated electronic energies. We also show how some small changes to FFSiOH could improve its performance for higher degrees of hydroxylation.

Computational Study of Acidic and Basic Functionalized Crystalline Silica Surfaces as a Model for Biomaterial Interfaces.

In silico modeling of acidic (CH2COOH) or basic (CH2NH2) functionalized silica surfaces has been carried out by means of a density functional approach based on a gradient-corrected functional to provide insight into the characterization of experimentally functionalized surfaces via a plasma method. Hydroxylated surfaces of crystalline cristobalite (sporting 4.8 OH/nm(2)) mimic an amorphous silica interface as unsubstituted material. To functionalize the silica surface we transformed the surface Si-OH groups into Si-CH2COOH and Si-CH2NH2 moieties to represent acidic/basic chemical character for the substitution. Structures, energetics, electronic, and vibrational properties were computed and compared as a function of the increasing loading of the functional groups (from 1 to 4 per surface unit cell). Classical molecular dynamics simulations of selected cases have been performed through a Reax-FF reactive force field to assess the mobility of the surface added chains. Both DFT and force field calculations identify the CH2NH2 moderate surface loading (1 group per unit cell) as the most stable functionalization, at variance with the case of the CH2COOH group, where higher loadings are preferred (2 groups per unit cell). The vibrational fingerprints of the surface functionalities, which are the ν(C═O) stretching and δ(NH2) bending modes for acidic/basic cases, have been characterized as a function of substitution percentage in order to guide the assignment of the experimental data. The final results highlighted the different behavior of the two types of functionalization. On the one hand, the frequency associated with the ν(C═O) mode shifts to lower wavenumbers as a function of the H-bond strength between the surface functionalities (both COOH and SiOH groups), and on the other hand, the δ(NH2) frequency shift seems to be caused by a subtle balance between the H-bond donor and acceptor abilities of the NH2 moiety. Both sets of data are in general agreement with experimental measurements on the corresponding silica-functionalized materials and provide finer details for a deeper interpretation of experimental spectra.

Interstellar H adsorption and H2 formation on the crystalline (010) forsterite surface: a B3LYP-D2* periodic study.

The physisorption/chemisorption of atomic hydrogen on a slab model of the Mg2SiO4 forsterite (010) surface mimicking the interstellar dust particle surface has been modeled using a quantum mechanical approach based on periodic B3LYP-D2* density functional calculations (DFT) combined with flexible polarized Gaussian type basis sets, which allows a balanced description of the hydrogen/surface interactions for both minima and activated complexes. Physisorption of hydrogen is barrierless, very weak and occurs either close to surface oxygen atoms or on Mg surface ions. The contribution of dispersion interactions accounts for almost half of the adsorption energy. Both the hydrogen adsorption energy and barrier to hydrogen jump between equivalent surface sites are overestimated compared to experimental results meant to simulate the interstellar conditions in the laboratory. The hydrogen atom exclusively chemisorbs at the oxygen site of the forsterite (010) surface, forming a SiOH surface group and its spin density being entirely transferred to the neighboring Mg ion. Barrier for chemisorption allows rapid attachment of H at the surface at 100 K, but prevents the same process from occurring at 10 K. From this H-chemisorbed state, the second hydrogen chemisorption mainly occurs on the neighboring Mg ion, thus forming a Mg-H surface group, giving rise to a surface species stabilized by favorable electrostatic interactions between the OHH-Mg pair. The formation of molecular hydrogen at the (010) forsterite surface adopting a Langmuir-Hinshelwood mechanism takes place either starting from two physisorbed H atoms with an almost negligible kinetic barrier through a spin-spin coupling driven reaction or from two chemisorbed H atoms with a barrier surmountable even at T higher than 10 K. We also suggest that a nanosized model of the interstellar dust built from a replica of the forsterite unit cell is able to adsorb half the energy released by the H2 formation by increasing its temperature by about 50 K which could then radiate in about 0.02 s.

B3LYP periodic study of the physicochemical properties of the nonpolar (010) Mg-pure and fe-containing olivine surfaces.

B3LYP periodic simulations have been carried out to study some physicochemical properties of the bulk structures and the corresponding nonpolar (010) surfaces of Mg-pure and Fe-containing olivine systems; i.e., Mg2SiO4 (Fo) and Mg1.5Fe0.5SiO4 (Fo75). A detailed structural analysis of the (010) Fo and Fo75 surface models shows the presence of coordinatively unsaturated metal cations (Mg(2+) and Fe(2+), respectively) with shorter metal-O distances compared to the bulk ones. Energetic analysis devoted to the Fe(2+) electronic spin configuration and to the ion position in the surfaces reveals that Fe(2+) in its quintet state and placed at the outermost positions of the slab constitutes the most stable Fe-containing surface, which is related to the higher stability of high spin states when Fe(2+) is coordinatively unsaturated. Comparison of the simulated IR and the corresponding reflectance spectra indicates that Fe(2+) substitution induces an overall bathochromic shift of the spectra due to the larger mass of Fe compared to Mg cation. In contrast, the IR spectra of the surfaces are shifted to upper values and exhibit more bands compared to the corresponding bulk systems due to the shorter metal-O distances given in the coordinatively unsaturated metals and to symmetry reduction which brings nonequivalent motions between the outermost and the internal modes, respectively.

Silica-based materials as drug adsorbents: first principle investigation on the role of water microsolvation on Ibuprofen adsorption.

Silica-based materials find applications as excipients and, particularly for those of mesoporous nature, as drug delivery agents for pharmaceutical formulations. Their performance can be crucially affected by water moisture, as it can modify the behavior of these formulations, by limiting their shelf life. Here we describe the role of water microsolvation on the features of ibuprofen adsorbed on a model of amorphous silica surface by means of density functional theory (DFT) simulations. Starting from the results of the simulation of ibuprofen in interaction with a dry hydrophobic amorphous silica surface, a limited number of water molecules has been added to study the configurational landscape of the microsolvated system. Structural and energetics properties, as well as the role of dispersive forces, have been investigated. Our simulations have revealed that the silica surface exhibits a higher affinity for water than for ibuprofen, even if several structures coexist at room temperature, with an active competition of ibuprofen and water for the exposed surface silanols. Dispersive interactions play a key role in this system, as pure DFT fails to correctly describe its potential energy surface. Indeed, van der Waals forces are the leading contribution to adsorption, independently of whether the drug is hydrogen-bonded directly to the surface or via water molecules.