In silico studies of the interactions between propofol and fentanyl using Gaussian accelerated molecular dynamics

Fentanyl is a potent opioid analgesic, which for decades has been used routinely in surgical and therapeutic applications. In addition to its analgesic properties, fentanyl also possesses anesthetic properties, which are not well understood. Fentanyl is used in the general anesthesia process to induce and maintain anesthesia in combination with the general anesthetic propofol, which fentanyl is known to potentiate. As the atomic-level mechanism behind the potentiation of propofol is unclear, we have used classical molecular dynamics simulations to study the interactions of these drugs with the Gloeobacter violaceus ion channel (GLIC). This ion channel has been identified as a target for many anesthetic drugs. We identified multiple binding sites using flooding style and Gaussian accelerated molecular dynamics (GaMD) simulations, showing fentanyl acting as a stabiliser that holds propofol within binding sites. Our extensive GaMD simulations were also able to show the pathway by which propofol blocks the channel pore, which has previously been suggested as a mechanism for ion channel modulation. General anesthesia is a multi-drug process and this study provides the first insight into the interactions between two different drugs in the anesthesia process in a relevant biological environment.Communicated by Ramaswamy H. Sarma.

Adsorbate-Induced Segregation of Cobalt from PtCo Nanoparticles: Modeling Au Doping and Core AuCo Alloying for the Improvement of Fuel Cell Cathode Catalysts

Platinum, when used as a cathode material for the oxygen reduction reaction, suffers from high overpotential and possible dissolution, in addition to the scarcity of the metal and resulting cost. Although the introduction of cobalt has been reported to improve reaction kinetics and decrease the precious metal loading, surface segregation or complete leakage of Co atoms causes degradation of the membrane electrode assembly, and either of these scenarios of structural rearrangement eventually decreases catalytic power. Ternary PtCo alloys with noble metals could possibly maintain activity with a higher dissolution potential. First-principles-based theoretical methods are utilized to identify the critical factors affecting segregation in Pt-Co binary and Pt-Co-Au ternary nanoparticles in the presence of oxidizing species. With a decreasing share of Pt, surface segregation of Co atoms was already found to become thermodynamically viable in the PtCo systems at low oxygen concentrations, which is assigned to high charge transfer between species. While the introduction of gold as a dopant caused structural changes that favor segregation of Co, creation of CoAu alloy core is calculated to significantly suppress Co leakage through modification of the electronic properties. The theoretical framework of geometrically different ternary systems provides a new route for the rational design of oxygen reduction catalysts.

Interaction of hydrogen with actinide dioxide (011) surfaces

The corrosion and oxidation of actinide metals, leading to the formation of metal-oxide surface layers with the catalytic evolution of hydrogen, impacts the management of nuclear materials. Here, the interaction of hydrogen with actinide dioxide (AnO₂, An = U, Np, or Pu) (011) surfaces by Hubbard corrected density functional theory (PBEsol+U) has been studied, including spin-orbit interactions and non-collinear 3k anti-ferromagnetic behavior. The actinide dioxides crystalize in the fluorite-type structure, and although the (111) surface dominates the crystal morphology, the (011) surface energetics may lead to more significant interaction with hydrogen. The dissociative adsorption of hydrogen on the UO₂ (0.44 eV), NpO₂ (-0.47 eV), and PuO₂ (-1.71 eV) (011) surfaces has been calculated. It is found that hydrogen dissociates on the PuO₂ (011) surface; however, UO₂ (011) and NpO₂ (011) surfaces are relatively inert. Recombination of hydrogen ions is likely to occur on the UO₂ (011) and NpO₂ (011) surfaces, whereas hydroxide formation is shown to occur on the PuO₂ (011) surface, which distorts the surface structure.

Atomistic Molecular Dynamics Simulations of Propofol and Fentanyl in Phosphatidylcholine Lipid Bilayers

Atomistic molecular dynamics (MD) and steered MD simulations in combination with umbrella sampling methodology were utilized to study the general anesthetic propofol and the opioid analgesic fentanyl and their interaction with lipid bilayers, which is not yet fully understood. These molecules were inserted into two different fully hydrated phospholipid bilayers, namely, dioleoylphosphatidylcholine (DOPC) and dipalmitoylphosphatidylcholine (DPPC), to investigate the effects that these drugs have on the bilayer. We determined the role of the lipid chain length and saturation on the behavior of the two drugs. Pure, fully hydrated DOPC and DPPC bilayers were also simulated, and the results were in excellent agreement with the experimental values. Various structural and mechanical properties of each system, such as the area per lipid, area compressibility modulus, order parameter, lateral lipid diffusion, hydrogen bonds, and radial distribution functions, have been calculated to assess how the drug molecules affect the different bilayers. From the calculated results, we show that fentanyl and propofol generally follow similar trends in each bilayer but adopt different favorable positions close to the headgroup/chain interface at the carbonyl groups. Propofol was shown to selectively form hydrogen bonds at the carbonyl carbon in each bilayer, whereas fentanyl interacts with water molecules at the headgroup interface. From the calculated free-energy profiles, we determined that both molecules show a preference for the low-density, low-order acyl chain region of the bilayers and both significantly preferred the DOPC bilayer with propofol and fentanyl having energy minima at -6.66 and -43.07 kcal mol^-1, respectively. This study suggests that different chain lengths and levels of saturation directly affect the properties of these two important molecules, which are seen to work together to control anesthesia in surgical applications.

Towards a morphology of cobalt nanoparticles: size and strain effects

Cobalt nanoparticles with diameters of 8 nm have recently shown promising performance for biomedical applications. However, it is still unclear how the shape of cobalt clusters changes with size when reaching the nanoparticle range. In the present work, density functional theory calculations have been employed to compare the stabilities of two non-crystalline (icosahedron and decahedron) shapes, and three crystalline motifs (hcp, fcc, and bcc) for magic numbered cobalt clusters with up to 1500 atoms, based on the changes in the cohesive energies, coordination numbers, and nearest-neighbour distances arising from varying geometries. Obtained trends were extrapolated to a 10⁴ size range, and an icosahedral shape was predicted for clusters up to 5500 atoms. Larger sized clusters adopt hcp stacking, in correspondence with the bulk phase. To explain the crystalline/non-crystalline crossovers, the contributions of the elastic strain density and twin boundary from the specimen surfaces to the cohesive energy of different motifs were evaluated. These results are expected to aid the design and synthesis of cobalt nanoparticles for applications ranging from catalysis to biomedical treatments.

Ethylene carbonate adsorption on the major surfaces of lithium manganese oxide Li1-xMn2O4 spinel (0.000 < x < 0.375): a DFT+U-D3 study

Understanding the surface reactivity of the commercial cathode material LiMn₂O₄ towards the electrolyte is important to improve the cycling performance of secondary lithium-ion batteries and to prevent manganese dissolution. In this work, we have employed spin-polarized density functional theory calculations with on-site Coulomb interactions and long-range dispersion corrections [DFT+U-D3-(BJ)] to investigate the adsorption of the electrolyte component ethylene carbonate (EC) onto the (001), (011) and (111) surfaces of the fully lithiated and partially delithiated Li_1-xMn₂O₄ spinel (0.000 < x < 0.375). The surface interactions were investigated by evaluating the adsorption energies of the EC molecule and the surface free energies. Furthermore, we analyzed the impact of EC adsorption on the Wulff crystal morphologies, the molecular vibrational frequencies and the adsorbate/surface charge transfers. The adsorption energies indicate that the EC molecule strongly adsorbs on the (111) facet, which is attributed to a bidentate binding configuration. We found that EC adsorption enhances the stability of the (111) facet, as shown by the Wulff crystal morphologies. Although a negligible charge transfer was calculated between the spinel surfaces and the EC molecule, a large charge rearrangement takes place within the surfactant upon adsorption. The wavenumbers of the C[double bond, length as m-dash]O stretching mode for the interacting EC molecule are red-shifted with respect to the isolated adsorbate, suggesting that this bond becomes weaker. The surface free energies show that both the fully lithiated and partially delithiated forms of the LiMn₂O₄ surfaces are stabilized by the EC molecule.

Electronic structure, ion diffusion and cation doping in the Na4VO(PO4)2 compound as a cathode material for Na-ion batteries

Sodium-ion batteries are considered one of the most promising alternatives to lithium-ion batteries owing to the low cost and wide abundance of sodium. Phosphate compounds are promising materials for sodium-ion batteries because of their high structural stability, energy densities and capacities. Vanadium phosphates have shown high energy densities, but their sodium-ion diffusion and cation doping properties are not fully rationalized. In this work, we combine density functional theory calculations and molecular dynamics simulations to study the electronic structure, ion diffusion and cation doping properties of the Na₄VO(PO₄)₂ compound. The calculated Na-ion activation energy of this compound is 0.49 eV, which is typical for Na-based cathode materials, and the simulations predict a Na-ion diffusion coefficient of 5.1 × 10^-11 cm² s^-1. The cell voltage trends show a voltage of 3.3 V vs. Na/Na⁺. Partial substitution of vanadium atoms by other metals (Al³⁺, Co²⁺, Fe³⁺, Mn⁴⁺, Ni²⁺ or Ti⁴⁺) increases the cell voltage up to 1.1 V vs. Na/Na⁺. These new insights will help us to understand the ion transport and electrochemical behaviour of potential phosphate cathode materials for sodium-ion batteries.

Binding modes of carboxylic acids on cobalt nanoparticles

Owing to their high saturation magnetisation, cobalt nanoparticles hold significant potential for the hyperthermia treatment of tumours. Covalent binding of carboxylic acids to the nanoparticles can induce biocompatibility, whilst also preventing the formation of surface oxides which reduce the magnetic properties of cobalt. Understanding the origin of the acid-metal interaction is key, yet probably the most experimentally challenging step, for the rational design of such entities. In this density functional theory study, we use static calculations to establish that a 57-atom Co cluster is the smallest model able to reproduce the adsorption behaviour of carboxylic acids, and ab initio metadynamics to obtain the structure and the free energy landscape for its interaction with valeric acid. Our simulations show that a bridging bidentate binding mode has a stronger affinity compared to monodentate binding, with energetically high transition barriers between the two. A chelate interaction mode of two carboxyl oxygen atoms can be formed as an intermediate. These results clarify the organic-inorganic interactions in the cobalt-acid system, providing a basis for the rational design of biocompatible metallic nanoparticles.

Nanostructured zeolite with brain-coral morphology and tailored acidity: a self-organized hierarchical porous material with MFI topology

Nano-zeolite with brain-coral morphology formed by self-organization of ultra-small nanospheres, exhibits micro/meso porosity with high surface area, distributed acid sites, and reduced diffusion resistance making it a promising solid acid catalyst.

Combined density functional theory and molecular dynamics study of Sm0.75A0.25Co1−xMnxO2.88 (A = Ca, Sr; x = 0.125, 0.25) cathode material for next generation solid oxide fuel cell

Computational study of novel next-generation SOFC cathode Sm_0.75(Ca,Sr)_0.25Mn_xCo_1−xO_2.88 showing fast electronic and ionic conduction in bulk.

Interaction of H2O with the Platinum Pt (001), (011), and (111) Surfaces: A Density Functional Theory Study with Long-Range Dispersion Corrections

Platinum is a noble metal that is widely used for the electrocatalytic production of hydrogen, but the surface reactivity of platinum toward water is not yet fully understood, even though the effect of water adsorption on the surface free energy of Pt is important in the interpretation of the morphology and catalytic properties of this metal. In this study, we have carried out density functional theory calculations with long-range dispersion corrections [DFT-D3-(BJ)] to investigate the interaction of H₂O with the Pt (001), (011), and (111) surfaces. During the adsorption of a single H₂O molecule on various Pt surfaces, it was found that the lowest adsorption energy (E _ads) was obtained for the dissociative adsorption of H₂O on the (001) surface, followed by the (011) and (111) surfaces. When the surface coverage was increased up to a monolayer, we noted an increase in E _ads/H₂O with increasing coverage for the (001) surface, while for the (011) and (111) surfaces, E _ads/H₂O decreased. Considering experimental conditions, we observed that the highest coverage was obtained on the (011) surface, followed by the (111) and (001) surfaces. However, with an increase in temperature, the surface coverage decreased on all the surfaces. Total desorption occurred at temperatures higher than 400 K for the (011) and (111) surfaces, but above 850 K for the (001) surface. From the morphology analysis of the Pt nanoparticle, we noted that, when the temperature increased, only the electrocatalytically active (111) surface remained.

Modulation of the Gloeobacter violaceus Ion Channel by Fentanyl: A Molecular Dynamics Study

Fentanyl is an opioid analgesic, which is routinely used in general surgery to suppress the sensation of pain and as the analgesic component in the induction and maintenance of anesthesia. Fentanyl is also used as the main component to induce anesthesia and as a potentiator to the general anesthetic propofol. The mechanism by which fentanyl induces its anesthetic action is still unclear, and we have therefore employed fully atomistic molecular dynamics simulations to probe this process by simulating the interactions of fentanyl with the Gloeobacter violaceus ligand-gated ion channel (GLIC). In this paper, we identify multiple extracellular fentanyl binding sites, which are different from the transmembrane general anesthetic binding sites observed for propofol and other general anesthetics. Our simulations identify a novel fentanyl binding site within the GLIC that results in conformational changes that inhibit conduction through the channel.

Dependence of electron transfer dynamics on the number of graphene layers in π-stacked 2D materials: insights from ab initio nonadiabatic molecular dynamics

Recent time-resolved transient absorption studies demonstrated that the rate of photoinduced interfacial charge transfer (CT) from Zn-phthalocyanine (ZnPc) to single-layer graphene (SLG) is faster than to double-layer graphene (DLG), in contrast to the expectation from Fermi's golden rule. We present the first time-domain non-adiabatic molecular dynamics (NA-MD) study of the electron injection process from photoexcited ZnPc molecules into SLG and DLG substrates. Our calculations suggest that CT occurs faster in the ZnPc/SLG system than in the ZnPc/DLG system, with 580 fs and 810 fs being the fastest components of the observed CT timescales, respectively. The computed timescales are in close agreement with those reported in the experiment. The computed CT timescales are determined largely by the magnitudes of the non-adiabatic couplings (NAC), which we find to be 4 meV and 2 meV, for the ZnPc/SLG and ZnPc/DLG systems, respectively. The transitions are driven mainly by the ZnPc out-of-plane bending mode at 1100 cm-1 and an overtone of fundamental modes in graphene at 2450 cm-1. We find that dephasing occurs on the timescale of 20 fs and is similar in both systems, so decoherence does not notably change the qualitative trends in the CT timescales. We highlight the importance of proper energy level alignment for capturing the qualitative trends in the CT dynamics observed in experiment. In addition, we illustrate several methodological points that are important for accurately modeling nonadiabatic dynamics in the ZnPc/FLG systems, such as the choice of surface hopping methodology, the use of phase corrections, NAC scaling, and the inclusion of Hubbard terms in the density functional and molecular dynamics calculations.

The Origin of High Activity of Amorphous MoS2 in the Hydrogen Evolution Reaction

Molybdenum disulfide (MoS₂ ) and related transition metal chalcogenides can replace expensive precious metal catalysts such as Pt for the hydrogen evolution reaction (HER). The relations between the nanoscale properties and HER activity of well-controlled 2H and Li-promoted 1T phases of MoS₂ , as well as an amorphous MoS₂ phase, have been investigated and a detailed comparison is made on Mo-S and Mo-Mo bond analysis under operando HER conditions, which reveals a similar bond structure in 1T and amorphous MoS₂ phases as a key feature in explaining their increased HER activity. Whereas the distinct bond structure in 1T phase MoS₂ is caused by Li⁺ intercalation and disappears under harsh HER conditions, amorphous MoS₂ maintains its intrinsic short Mo-Mo bond feature and, with that, its high HER activity. Quantum-chemical calculations indicate similar electronic structures of small MoS₂ clusters serving as models for amorphous MoS₂ and the 1T phase MoS₂ , showing similar Gibbs free energies for hydrogen adsorption (ΔG_H* ) and metallic character.

Glucosepane is associated with changes to structural and physical properties of collagen fibrils

Collagen glycation, and in particular the formation of advanced glycation end-product (AGE) crosslinks, plays a central role in the ageing process and in many of the long-term complications of diabetes. Glucosepane, the most abundant and relevant AGE crosslink, has been suggested to increase the stiffness of tissue and reduce its solubility, although no evidence is available concerning the mechanisms. We have used a combination of computational and experimental techniques to study a collagen-rich tissue with a relatively simple organisation to further our understanding of the impact of glucosepane on the structural and physical properties of collagen fibrils. Our work shows that glucosepane levels increase dramatically in aged tendon tissue and are associated with the reduced density of collagen packing and increased porosity to water molecules. Our studies provide the basis to understand many of the tissue dysfunctions associated with ageing and diabetes across a range of different tissues types.

•

Levels of the advanced glycation end-product glucosepane increase in human tendon with increasing chronological age.

•

Glucosepane results in a less tightly held helical structure in the collagen molecule and increased porosity to water.

•

Water content is higher in Achilles and anterior tibialis tendon tissue from older individuals compared to young people.

•

The denaturation temperature of collagen increases in the older age group suggesting a more highly cross-linked structure.

•

The enthalpy of collagen denaturation decreases in older donors suggesting molecules are less confined within the fibril.

Tuning ZnO Sensors Reactivity toward Volatile Organic Compounds via Ag Doping and Nanoparticle Functionalization

Nanomaterials for highly selective and sensitive sensors toward specific gas molecules of volatile organic compounds (VOCs) are most important in developing new-generation of detector devices, for example, for biomarkers of diseases as well as for continuous air quality monitoring. Here, we present an innovative preparation approach for engineering sensors, which allow for full control of the dopant concentrations and the nanoparticles functionalization of columnar material surfaces. The main outcome of this powerful design concept lies in fine-tuning the reactivity of the sensor surfaces toward the VOCs of interest. First, nanocolumnar and well-distributed Ag-doped zinc oxide (ZnO:Ag) thin films are synthesized from chemical solution, and, at a second stage, noble nanoparticles of the required size are deposited using a gas aggregation source, ensuring that no percolating paths are formed between them. Typical samples that were investigated are Ag-doped and Ag nanoparticle-functionalized ZnO:Ag nanocolumnar films. The highest responses to VOCs, in particular to (CH₃)₂CHOH, were obtained at a low operating temperature (250 °C) for the samples synergistically enhanced with dopants and nanoparticles simultaneously. In addition, the response times, particularly the recovery times, are greatly reduced for the fully modified nanocolumnar thin films for a wide range of operating temperatures. The adsorption of propanol, acetone, methane, and hydrogen at various surface sites of the Ag-doped Ag₈/ZnO(0001) surface has been examined with the density functional theory (DFT) calculations to understand the preference for organic compounds and to confirm experimental results. The response of the synergistically enhanced sensors to gas molecules containing certain functional groups is in excellent agreement with density functional theory calculations performed in this work too. This new fabrication strategy can underpin the next generation of advanced materials for gas sensing applications and prevent VOC levels that are hazardous to human health and can cause environmental damages.

Thermal Properties and Segregation Behavior of Pt Nanowires Modified with Au, Ag, and Pd Atoms: A Classical Molecular Dynamics Study

Platinum nanowires (NWs) have been reported to be catalytically active toward the oxygen reduction reaction (ORR). The edge modification of Pt NWs with metals M (M = Au, Ag, or Pd) may have a positive impact on the overall ORR activity by facilitating diffusion of adsorbed oxygen, O_ads, and hydroxyl groups, OH_ads, between the {001} and {111} terraces. In the present study, we have employed classical molecular dynamics simulations to investigate the segregation behavior of Au, Ag, and Pd decorating the edges of Pt NWs. We observe that, under vacuum conditions, Pd prefers to diffuse toward the core rather than stay on the NW surface. Ag and Au atoms are mobile at temperatures as low as 900 K; they remain on the surface but do not appear to be preferentially more stable at edge sites. To effect segregation of Au and Ag atoms toward the edge, we propose annealing in the presence of different reactive gas environments. Overall, our study suggests potential experimental steps required for the synthesis of Pt nanowires and nanoparticles with improved O_ads and OH_ads interfacet diffusion rates and consequently an improved ORR activity.

Fe(ii) and Fe(iii) dithiocarbamate complexes as single source precursors to nanoscale iron sulfides: a combined synthetic and in situ XAS approach

Nanoparticulate iron sulfides have many potential applications and are also proposed to be prebiotic catalysts for the reduction of CO₂ to biologically important molecules, thus the development of reliable routes to specific phases with controlled sizes and morphologies is important. Here we focus on the use of iron dithiocarbamate complexes as single source precursors (SSPs) to generate greigite and pyrrhotite nanoparticles. Since these minerals contain both iron(iii) and iron(ii) centres, SSPs in both oxidation states, [Fe(S₂CNR₂)₃] and cis-[Fe(CO)₂(S₂CNR₂)₂] respectively, have been utilised. Use of this Fe(ii) precursor is novel and it readily loses both carbonyls in a single step (as shown by TGA measurements) providing an in situ source of the extremely air-sensitive Fe(ii) dithiocarbamate complexes [Fe(S₂CNR₂)₂]. Decomposition of [Fe(S₂CNR₂)₃] alone in oleylamine affords primarily pyrrhotite, although by careful control of reaction conditions (ca. 230 °C, 40-50 nM SSP) a window exists in which pure greigite nanoparticles can be isolated. With cis-[Fe(CO)₂(S₂CNR₂)₂] we were unable to produce pure greigite, with pyrrhotite formation dominating, a similar situation being found with mixtures of Fe(ii) and Fe(iii) precursors. In situ X-ray absorption spectroscopy (XAS) studies showed that heating [Fe(S₂CNⁱBu₂)₃] in oleylamine resulted in amine coordination and, at ca. 60 °C, reduction of Fe(iii) to Fe(ii) with (proposed) elimination of thiuram disulfide (S₂CNR₂)₂. We thus carried out a series of decomposition studies with added thiuram disulfide (R = ⁱBu) and found that addition of 1-2 equivalents led to the formation of pure greigite nanoparticles between 230 and 280 °C with low SSP concentrations. Average particle size does not vary significantly with increasing concentration, thus providing a convenient route to ca. 40 nm greigite nanoparticles. In situ XAS studies have been carried out and allow a decomposition pathway for [Fe(S₂CNⁱBu₂)₃] in oleylamine to be established; reduction of Fe(iii) to Fe(ii) reduction triggers substitution of the secondary amide backbone by oleylamine (RNH₂) resulting in the in situ formation of a primary dithiocarbamate derivative [Fe(RNH₂)₂(S₂CNHR)₂]. This in turn extrudes RNCS to afford molecular precursors of the observed FeS nanomaterials. The precise role of thiuram disulfide in the decomposition process is unknown, but it likely plays a part in controlling the Fe(iii)-Fe(ii) equilibrium and may also act as a source of sulfur allowing control over the Fe : S ratio in the mineral products.

Unraveling the Role of Lithium in Enhancing the Hydrogen Evolution Activity of MoS2: Intercalation versus Adsorption

Molybdenum disulfide (MoS₂) is a highly promising catalyst for the hydrogen evolution reaction (HER) to realize large-scale artificial photosynthesis. The metallic 1T'-MoS₂ phase, which is stabilized via the adsorption or intercalation of small molecules or cations such as Li, shows exceptionally high HER activity, comparable to that of noble metals, but the effect of cation adsorption on HER performance has not yet been resolved. Here we investigate in detail the effect of Li adsorption and intercalation on the proton reduction properties of MoS₂. By combining spectroscopy methods (infrared of adsorbed NO, ⁷Li solid-state nuclear magnetic resonance, and X-ray photoemission and absorption) with catalytic activity measurements and theoretical modeling, we infer that the enhanced HER performance of Li _x MoS₂ is predominantly due to the catalytic promotion of edge sites by Li.

First-principles DFT insights into the structural, elastic, and optoelectronic properties of α and β-ZnP2: implications for photovoltaic applications

Binary II-V semiconductors are highly optically active materials, possess high intrinsic mechanical and chemical durability, and have electronic properties ideal for optoelectronic applications. Among them, zinc diphosphide (ZnP₂) is a promising earth-abundant absorber material for solar energy conversion. We have investigated the structural, mechanical, and optoelectronic properties of both the tetragonal (α) and monoclinic (β) phases of ZnP₂ using standard, Hubbard-corrected and screened hybrid density functional theory methods. Through the analysis of bond character, band gap nature, and absorption spectra, we show that there exist two polymorphs of the β phase (denoted as β ₁ and β ₂) with distinct differences in the photovoltaic potential. While β ₁ exhibits the characteristics of metallic compounds, β ₂ is a semiconductor with predicted thin-film photovoltaic absorbing efficiency of almost 10%. The α phase is anticipated to be an indirect gap material with a calculated efficiency limited to only 1%. We have also analysed and gained insights into the electron localization function, projected density of states and projected crystal orbital Hamilton populations for the analogue bonds between the α and β-ZnP₂. In light of these calculations, a number of previous discrepancies have been solved and a solid ground for future employment of zinc diphosphides in photovoltaics has been established.

Computational study of the mixed B-site perovskite SmBxCo1-xO3-d (B = Mn, Fe, Ni, Cu) for next generation solid oxide fuel cell cathodes

SmCoO3 is a promising perovskite material for the next generation of intermediate temperature solid oxide fuel cells (SOFC), but its potential application is directly linked to, and dependent on, the presence of dopant ions. Doping on the Co-site is suggested to improve the catalytic and electronic properties of this cathode material. Fe, Mn, Ni, and Cu have been proposed as possible dopants and experimental studies have investigated and confirmed the potential of these materials. Here we present a systematic DFT+U study focused on the changes in electronic, magnetic, and physical properties with B-site doping of SmCoO3 to allow cathode optimization. It is shown that doping generally leads to distortion in the system, thereby inducing different electron occupations of the Co d-orbitals, altering the electronic and magnetic structure. From these calculations, the 0 K electronic conductivity (σe) was obtained, with SmMnxCo1-xO3 having the highest σe, and SmFexCo1-xO3 the lowest σe, in agreement with experiment. We have also investigated the impact of dopant species and concentration on the oxygen vacancy formation energy (Ef), which is related to the ionic conductivity (σO). We found that the Ef values are lowered only when SmCoO3 is doped with Cu or Ni. Finally, thermal expansion coefficients were calculated, with Mn-doping showing the largest decrease at low x and at x = 0.75. Combining these results, it is clear that Mn-doping in the range x = 0.125-0.25 would imbue SmCoO3 with the most favorable properties for IT-SOFC cathode applications.

Advances in Sustainable Catalysis: A Computational Perspective

The enormous challenge of moving our societies to a more sustainable future offers several exciting opportunities for computational chemists. The first principles approach to "catalysis by design" will enable new and much greener chemical routes to produce vital fuels and fine chemicals. This prospective outlines a wide variety of case studies to underscore how the use of theoretical techniques, from QM/MM to unrestricted DFT and periodic boundary conditions, can be applied to biocatalysis and to both homogeneous and heterogenous catalysts of all sizes and morphologies to provide invaluable insights into the reaction mechanisms they catalyze.

Hydrogen adsorption on transition metal carbides: a DFT study

Transition metal carbides are a class of materials widely known for both their interesting physical properties and catalytic activity. In this work, we have used plane-wave DFT methods to study the interaction with increasing amounts of molecular hydrogen on the low-index surfaces of four major carbides - TiC, VC, ZrC and NbC. Adsorption is found to be generally exothermic and occurs predominantly on the surface carbon atoms. We identify trends over the carbides and their surfaces for the energetics of the adsorption as a function of their electronic and geometrical characteristics. An ab initio thermodynamics formalism is used to study the properties of the slabs as the hydrogen coverage is increased.