DFT study of the CO oxidation on the Au(321) surface

The CO oxidation on the Au(321) surface was investigated using spin polarized density functional theory based calculations within the GGA-PW91 exchange-correlation functional. This was done by studying separately the adsorption of isolated CO or CO2 and also the coadsorption of CO + O or CO + O-2 on the Au(321) surface. A periodic supercell approach was used to model the gold surface. The kinetic profile of the oxidation reaction was determined with the climbing image-nudged elastic band method and also with the dimer approach. It was found that CO adsorbs on the clean surface preferably at the kinks, and the same preference exists if atomic or molecular oxygen is coadsorbed on the Au(321) surface. CO2 is weakly adsorbed on Au(321) and appears at large distance from the metal surface. Importantly, the formation of carbonate species or of four atoms compounds, OCOO, adsorbed on the Au(321) surface is thermodynamically favorable from CO and O-2. The reaction of CO oxidation by atomic oxygen occurs almost without any energy cost on a reconstructed surface, whereas moderate barriers of similar to 0.6 eV were computed for the direct reaction with molecular oxygen occurring at the surface steps. These results suggest that the predissociation of the molecular oxygen on the Au(321) surface for the CO oxidation is energetically less favorable than the direct reaction with molecular oxygen. Finally, the products of the oxidation reaction are much more stable than the four atoms compound.publishe

New Mechanistic Insights on Carbon Nanotubes’ Nanotoxicity Using Isolated Submitochondrial Particles, Molecular Docking, and Nano-QSTR Approaches

Single-walled carbon nanotubes can induce mitochondrial F0F1-ATPase nanotoxicity through inhibition. To completely characterize the mechanistic effect triggering the toxicity, we have developed a new approach based on the combination of experimental and computational study, since the use of only one or few techniques may not fully describe the phenomena. To this end, the in vitro inhibition responses in submitochondrial particles (SMP) was combined with docking, elastic network models, fractal surface analysis, and Nano-QSTR models. In vitro studies suggest that inhibition responses in SMP of F0F1-ATPase enzyme were strongly dependent on the concentration assay (from 3 to 5 µg/mL) for both pristine and COOH single-walled carbon nanotubes types (SWCNT). Besides, both SWCNTs show an interaction inhibition pattern mimicking the oligomycin A (the specific mitochondria F0F1-ATPase inhibitor blocking the c-ring F0 subunit). Performed docking studies denote the best crystallography binding pose obtained for the docking complexes based on the free energy of binding (FEB) fit well with the in vitro evidence from the thermodynamics point of view, following an affinity order such as: FEB (oligomycin A/F0-ATPase complex) = −9.8 kcal/mol > FEB (SWCNT-COOH/F0-ATPase complex) = −6.8 kcal/mol ~ FEB (SWCNT-pristine complex) = −5.9 kcal/mol, with predominance of van der Waals hydrophobic nano-interactions with key F0-ATPase binding site residues (Phe 55 and Phe 64). Elastic network models and fractal surface analysis were performed to study conformational perturbations induced by SWCNT. Our results suggest that interaction may be triggering abnormal allosteric responses and signals propagation in the inter-residue network, which could affect the substrate recognition ligand geometrical specificity of the F0F1-ATPase enzyme in order (SWCNT-pristine > SWCNT-COOH). In addition, Nano-QSTR models have been developed to predict toxicity induced by both SWCNTs, using results of in vitro and docking studies. Results show that this method may be used for the fast prediction of the nanotoxicity induced by SWCNT, avoiding time- and money-consuming techniques. Overall, the obtained results may open new avenues toward to the better understanding and prediction of new nanotoxicity mechanisms, rational drug design-based nanotechnology, and potential biomedical application in precision nanomedicineThis research was funded by FCT/MCTES through national funds (Michael González-Durruthy, Riccardo Concu, and M. Natália D.S. Cordeiro), grant UID/QUI/50006/2020, as well as by Xunta de Galicia (Juan M. Ruso), grant ED41E2018/08S

On the thickness of the double layer in ionic liquids

In this study, we examined the thickness of the electrical double layer (EDL) in ionic liquids using density functional theory (DFT) calculations and molecular dynamics (MD) simulations. We focused on the BF4- anion adsorption from 1-ethyl-3-methylimidazolium tetrafluoroborate (EMImBF4) ionic liquid on the Au(111) surface. At both DFT and MD levels, we evaluated the capacitance-potential dependence for the Helmholtz model of the interface. Using MD simulations, we also explored a more realistic, multilayer EDL model accounting for the ion layering. Concurrent analysis of the DFT and MD results provides a ground for thinking whether the electrical double layer in ionic liquids is one- or multi-ionic-layer thick

Dermic diffusion and stratum corneum: a state of the art review of mathematical models

Transdermal biotechnologies are an ever increasing field of interest, due to the medical and pharmaceutical applications that they underlie. There are several mathematical models at use that permit a more inclusive vision of pure experimental data and even allow practical extrapolation for new dermal diffusion methodologies. However, they grasp a complex variety of theories and assumptions that allocate their use for specific situations. Models based on Fick's First Law found better use in contexts where scaled particle theory Models would be extensive in time-span but the reciprocal is also true, as context of transdermal diffusion of particular active compounds changes. This article reviews extensively the various theoretical methodologies for studying dermic diffusion in the rate limiting dermic barrier, the stratum corneum, and systematizes its characteristics, their proper context of application, advantages and limitations, as well as future perspectives

Coupling of Cyclic Voltammetry and Electrochemical Impedance Spectroscopy for Probing the Thermodynamics of Facilitated Ion Transfer Reactions Exhibiting Chemical Kinetic Hindrances

Mathematical models under conditions of cyclic staircase voltammetry and electrochemical impedance spectroscopy (EIS), which consider the kinetic effects due to the complexation reaction by the facilitated transfer of metal ions at polarized interfaces, are presented. Criteria for qualitative recognition of these kinetic effects from the features of simulated cyclic voltammograms are given. In case of the existence of these effects, only the EIS can bring access to the thermodynamics and kinetics of the complexation chemical reaction. Analytical equations for estimating the thermodynamic parameters by such systems under EIS conditions are evaluated. The theoretical results are compared with the experimental results of the facilitated Cu2+ transfer at the polarized water-1,2-dichlorethane interface, assisted by two phenanthroline-containing macrocycles. In the experimental case where kinetic effects due to the complexation step exist, we show how elegantly EIS can be used as a tool for estimation of the complexation constant of Cu2+ and 5-oxo-2,8-dithia [9],(2,9)-1,10-phenanthrolinophane (PhenOS2)

Quinoxaline, its derivatives and applications: a state of the art review

Quinoxaline derivatives are an important class of heterocycle compounds, where N replaces some carbon atoms in the ring of naphthalene. Its molecular formula is C8H6N2, formed by the fusion of two aromatic rings, benzene and pyrazine. It is rare in natural state, but their synthesis is easy to perform. In this review the State of the Art will be presented, which includes a summary of the progress made over the past years in the knowledge of the structure and mechanism of the quinoxaline and quinoxaline derivatives, associated medical and biomedical value as well as industrial value. Modifying quinoxaline structure it is possible to obtain a wide variety of biomedical applications, namely antimicrobial activities and chronic and metabolic diseases treatment

Targeting Beta-Blocker Drug–Drug Interactions with Fibrinogen Blood Plasma Protein: A Computational and Experimental Study

In this work, one of the most prevalent polypharmacology drug–drug interaction events that occurs between two widely used beta-blocker drugs—i.e., acebutolol and propranolol—with the most abundant blood plasma fibrinogen protein was evaluated. Towards that end, molecular docking and Density Functional Theory (DFT) calculations were used as complementary tools. A fibrinogen crystallographic validation for the three best ranked binding-sites shows 100% of conformationally favored residues with total absence of restricted flexibility. From those three sites, results on both the binding-site druggability and ligand transport analysis-based free energy trajectories pointed out the most preferred biophysical environment site for drug–drug interactions. Furthermore, the total affinity for the stabilization of the drug–drug complexes was mostly influenced by steric energy contributions, based mainly on multiple hydrophobic contacts with critical residues (THR22: P and SER50: Q) in such best-ranked site. Additionally, the DFT calculations revealed that the beta-blocker drug–drug complexes have a spontaneous thermodynamic stabilization following the same affinity order obtained in the docking simulations, without covalent-bond formation between both interacting beta-blockers in the best-ranked site. Lastly, experimental ultrasound density and velocity measurements were performed and allowed us to validate and corroborate the computational obtained resultsThis research was funded by FCT/MCTES through national funds (Michael González-Durruthy, Riccardo Concu, and M. Natália D.S. Cordeiro), grant UID/QUI/50006/2020, as well as by Xunta de Galicia (Juan M. Ruso), grant ED41E2018/08S

Effect of the exchange-correlation potential on the transferability of Bronsted-Evans-Polanyi relationships in heterogeneous catalysis

As more and more accurate density functional methods emerge, the transferability of Bronsted-Evans-Polanyi (BEP) relationships obtained with previous models is an open question. In this work, BEP relationships derived from different density functional theory based calculations are analyzed to answer this question. In particular, BEP relationships linking the activation energy of O-H bond breaking reactions taking place on metallic surfaces with the adsorption energy of the reaction products are chosen as a case study. These relationships are obtained with the widely used Perdew-Wang (PW91) generalized gradient approximation (GGA) exchange-correlation functional and with the more accurate meta-GGA Tao-Perdew Staroverov-Scuseria (TPSS) one. We provide compelling evidence that BEP relationships derived from PW91 and TPSS functionals are essentially coincidental. This finding validates previously published BEP relationships and indicates that the reaction activation energy barrier can be obtained by the determination of the energy reaction descriptor value at the less computationally demanding GGA level; an important aspect to consider in future studies aimed at the computational design of catalysts :with improved characteristics

Molecular simulations of the synthesis of periodic mesoporous silica phases at high surfactant concentrations

Molecular dynamics simulations of a coarse-grained model are used to study the formation mechanism of periodic mesoporous silica over a wide range of cationic surfactant concentrations. This follows up on an earlier study of systems with low surfactant concentrations. We started by studying the phase diagram of the surfactant-water system and found that our model shows good qualitative agreement with experiments with respect to the surfactant concentrations where various phases appear. We then considered the impact of silicate species upon the morphologies formed. We have found that even in concentrated surfactant systems --in the concentration range where pure surfactant solutions yield a liquid crystal phase -- the liquid-crystal templating mechanism is not viable because the preformed liquid crystal collapses as silica monomers are added into the solution. Upon the addition of silica dimers, a new phase-separated hexagonal array is formed. The preformed liquid crystals were found to be unstable in the presence of monomeric silicates. In addition, the silica dimer is found to be essential for mesoscale ordering at both low and high surfactant concentrations. Our results support the view that a cooperative interaction of anionic silica oligomers and cationic surfactants determines the mesostructure formation in the M41S family of materials.publishe

Multiscale model for the templated synthesis of mesoporous silica: the essential role of silica oligomers

A detailed theoretical understanding of the synthesis mechanism of periodic mesoporous silica has not yet been achieved. We present results of a multiscale simulation strategy that, for the first time, describes the molecular-level processes behind the formation of silica/surfactant mesophases in the synthesis of templated MCM-41 materials. The parameters of a new coarse-grained explicit-solvent model for the synthesis solution are calibrated with reference to a detailed atomistic model, which itself is based on quantum mechanical calculations. This approach allows us to reach the necessary time and length scales to explicitly simulate the spontaneous formation of mesophase structures while maintaining a level of realism that allows for direct comparison with experimental systems. Our model shows that silica oligomers are a necessary component in the formation of hexagonal liquid crystals from low-concentration surfactant solutions. Because they are multiply charged, silica oligomers are able to bridge adjacent micelles, thus allowing them to overcome their mutual repulsion and form aggregates. This leads the system to phase separate into a dilute solution and a silica/surfactant-rich mesophase, which leads to MCM-41 formation. Before extensive silica condensation takes place, the mesophase structure can be controlled by manipulation of the synthesis conditions. Our modeling results are in close agreement with experimental observations and strongly support a cooperative mechanism for synthesis of this class of materials. This work paves the way for tailored design of nanoporous materials using computational models