Search

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.

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.

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.

Fully ab initio treatment of complex solid systems needs computational software which is able to efficiently take advantage of the growing power of high performance computing (HPC) architectures. Recent improvements in CRYSTAL, a periodic ab initio code that uses a Gaussian basis set, allows treatment of very large unit cells for crystalline systems on HPC architectures with high parallel efficiency in terms of running time and memory requirements. The latter is a crucial point, due to the trend toward architectures relying on a very high number of cores with associated relatively low memory availability. An exhaustive performance analysis shows that density functional calculations, based on a hybrid functional, of low-symmetry systems containing up to 100,000 atomic orbitals and 8000 atoms are feasible on the most advanced HPC architectures available to European researchers today, using thousands of processors., (Copyright © 2012 Wiley Periodicals, Inc.)