Density Functional Molecular Dynamics Simulations Research Articles

Formation of Supernarrow Borophene Nanoribbons

AbstractBorophenes have sparked considerable interest owing to their fascinating physical characteristics and diverse polymorphism. However, borophene nanoribbons (BNRs) with widths less than 2 nm have not been achieved. Herein, we report the experimental realization of supernarrow BNRs. Combining scanning tunneling microscopy imaging with density functional theory modeling and ab initio molecular dynamics simulations, we demonstrate that, under the applied growth conditions, boron atoms can penetrate the outermost layer of Au(111) and form BNRs composed of a pair of zigzag (2,2) boron rows. The BNRs have a width self‐contained to ∼1 nm and dipoles at the edges to keep them separated. They are embedded in the outermost Au layer and shielded on top by the evacuated Au atoms, free of the need for post‐passivation. Scanning tunneling spectroscopy reveals distinct edge states, primarily attributed to the localized spin at the BNRs’ zigzag edges. This work adds a new member to the boron material family and introduces a new physical feature to borophenes.

Formation of Supernarrow Borophene Nanoribbons.

Borophenes have sparked considerable interest owing to their fascinating physical characteristics and diverse polymorphism. However, borophene nanoribbons (BNRs) with widths less than 2 nm have not been achieved. Herein, we report the experimental realization of supernarrow BNRs. Combining scanning tunneling microscopy imaging with density functional theory modeling and ab initio molecular dynamics simulations, we demonstrate that, under the applied growth conditions, boron atoms can penetrate the outermost layer of Au(111) and form BNRs composed of a pair of zigzag (2,2) boron rows. The BNRs have a width self-contained to ∼1 nm and dipoles at the edges to keep them separated. They are embedded in the outermost Au layer and shielded on top by the evacuated Au atoms, free of the need for post-passivation. Scanning tunneling spectroscopy reveals distinct edge states, primarily attributed to the localized spin at the BNRs' zigzag edges. This work adds a new member to the boron material family and introduces a new physical feature to borophenes.

Exploring the mechanism of nerve agent (Tabun and Sarin) adsorption on carbon nanocones: Computational insights

Carbon NanoCones (CNCs), distinguished by their conical shape and large surface area, are evaluated for their potential application in nerve gas [tabun (GA) and sarin (GB)] adsorption. Density Functional Theory (DFT) calculations, Molecular Dynamics (MD), and Monte Carlo (MC) simulations are employed to investigate this phenomenon comprehensively. The study findings indicate that tabun (GA) and sarin (GB) undergo spontaneous adsorption, characterized by significant negative energy values. Density Functional Theory (DFT) calculations indicate electron donation tendencies of nerve agents to CNC surfaces. Computational analysis shows higher adsorption efficiency in the inner portion (P1) of CNCs. Molecular Dynamics (MD) and Monte Carlo (MC) simulations confirm planar adsorption configurations parallel to CNC surfaces. The DFT calculations show that GA has a far greater adsorption energy on a Carbon Nanocone compared to GB (−44.96 kcal/mol for GA, whereas GB has an adsorption energy of −35.95 kcal/mol). Furthermore, the analysis reveals significant differences between Sarin and Tabun nerve agents. Sarin exhibits higher global hardness, indicating increased stability, while Tabun shows higher overall softness, suggesting greater reactivity. Despite similar electronegativities, Sarin’s greater ionization energy and electron affinity imply enhanced stability and reactivity. Differences in HOMO and LUMO energies highlight unique reactivity profiles, with Sarin being more susceptible to nucleophilic attacks and Tabun to electrophilic attack. Sarin’s larger energy gap between HOMO and LUMO orbitals signifies its superior stability under electronic disturbances.

Revealing the role of interlayer spacing in radioactive-ion sieving of functionalized graphene membranes

Functionalization of graphene enables precise control over interlayer spacing during film formation, thereby enhancing the separation efficiency of radioactive ions in graphene membranes. However, the systematic impact of interlayer spacing of graphene membranes on radioactive-ion separation remains unexplored. This study aims to elucidate how interlayer spacing in functionalized graphene membranes affects the separation of radioactive ions. Utilizing polyamidoxime (PAO) to modify graphene oxide, we controlled the interlayer spacing of graphene membranes. Experimental results indicate that tuning interlayer spacing enables control of the permeation flux of radioactive ions (UO22+ 1.01 × 10−5–8.32 × 10−5 mol/m2·h, and K+ remains stable at 3.60 × 10−4 mol/m2·h), and the K+/UO22+ separation factors up to 36.2 at an interlayer spacing of 8.8 Å. Using density functional theory and molecular dynamics simulations, we discovered that the effective separation is mainly determined via interlayer spacing and the quantity of introduced functional groups, explaining the anomalous high permeation flux of target ions at low interlayer spacing (4.3 Å). This study deepens our comprehension of interlayer spacing within nanoconfined spaces for ion separation and recovery via graphene membranes, offering valuable insights for the design and synthesis of high-performance nanomembrane materials.

Computational Insights into Chromene/pyran Derivatives: Molecular Docking, ADMET Studies, DFT Calculations, and MD Simulations as Promising Candidates for Parkinson's Disease.

Parkinson's disease (PD) is a neurodegenerative condition characterized by both motor and non-motor symptoms. Although PD is commonly associated with a decline of dopaminergic neurons in the substantia nigra, other diagnostic criteria and biomarkers also exist. In the search for novel therapeutic agents, chromene and pyran derivatives have shown potential due to their diverse pharmacological activities. This study utilizes a comprehensive computational approach to investigate the viability of chromene/pyran compounds as potential treatments for PD. The drug-likeness characteristics of these molecules were analyzed using ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) studies. Molecular docking was performed against PDB ID: 2V5Z. The best three molecules chosen were compound 7, compound 24, and compound 67 have a binding energy of -6.7, -8.6, and -10.9 kcal/mol. Molecules demonstrating positive blood-brain barrier permeability, good solubility, and favorable binding affinity were further evaluated using Density Functional Theory (DFT) calculations and Molecular Dynamics (MD) simulations to assess their electronic structure and stability. DFT calculations indicated that molecule 82 has a dipole moment of 15.70 D. RMSD and RMSF results confirmed the stability of the complexes over a 100 ns simulation, with a maximum of 3 hydrogen bonds formed.

The cage effect of electron beam irradiation damage in cryo-electron microscopy

Electron beam irradiation can cause damage to biological and organic samples, as determined via transmission electron microscopy (TEM). Cryo-electron microscopy (cryo-EM) significantly reduces such damage by quickly freezing the environmental water around organic molecules. However, there are multiple hypotheses about the mechanism of cryo-protection in cryo-EM. A lower temperature can cause less molecular dissociation in the first stage, or frozen water can have a “cage” effect by preventing the dissociated fragments from flying away. In this work, we use real-time time-dependent density functional theory molecular dynamics(rt-TDDFT-MD) simulations to study the related dynamics. We use our recently developed natural orbital branching (NOB) algorithm to describe the molecular dissociation process after the molecule is ionized. We find that despite the difference in surrounding water molecules at different temperatures, the initial dissociation process is similar. On the other hand, the dissociated fragments fly away at room temperature, while they remain in the same cage when frozen water is used. Our results provide direct support for the cage effect mechanism.

Oligo-Adenine and Cyanuric Acid Supramolecular DNA-Based Hydrogels Exhibiting Acid-Resistance and Physiological pH-Responsiveness.

Expanding the functions and applications of DNA by integrating noncanonical bases and structures into biopolymers is a continuous scientific effort. An adenine-rich strand (A-strand) is introduced as functional scaffold revealing, in the presence of the low-molecular-weight cofactor cyanuric acid (CA, pKa 6.9), supramolecular hydrogel-forming efficacies demonstrating multiple pH-responsiveness. At pH 1.2, the A-strand transforms into a parallel A-motif duplex hydrogel cross-linked by AH+-H+A units due to the protonation of adenine (pKa 3.5). At pH 5.2, and in the presence of coadded CA, a helicene-like configuration is formed between adenine and protonated CA, generating a parallel A-CA triplex cross-linked hydrogel. At pH 8.0, the hydrogel undergoes transition into a liquid state by deprotonation of CA cofactor units and disassembly of A-CA triplex into its constituent components. Density functional theory calculations and molecular dynamics simulations, supporting the structural reconfigurations of A-strand in the presence of CA, are performed. The sequential pH-stimulated hydrogel states are rheometrically characterized. The hydrogel framework is loaded with fluorescein-labeled insulin, and the pH-stimulated release of insulin from the hydrogel across the pH barriers present in the gastrointestinal tract is demonstrated. The results provide principles for future application of the hydrogel for oral insulin administration for diabetes.

Ab initio exploration of modified carbon nanotubes as potential corrosion inhibitors

In order to develop novel unexplored potential corrosion inhibitors, covalently modified single-walled carbon nanotubes (SWCNT) via benzoic (–PhCOOH) and aniline (–PhNH2) groups are being investigated as corrosion inhibitors for the first time. Utilizing a comprehensive approach, this study employed density functional theory (DFT), Monte Carlo (MC), and molecular dynamics simulations (MD) to assess the adsorption behavior of modified nanotubes as corrosion inhibitors on the Cu(111) surface within a simulated aqueous HCl corrosion medium. The results provided molecular information on the adsorption capability, geometry adsorption centers, and adsorption energies (Eads) of carbon nanotubes on the surface of Cu(111). The adsorption energy values unveiled robust interactions between SWCNT–PhCOOH and SWCNT–PhNH2 inhibitors and the Cu(111) surface, suggesting a highly effective corrosion protection mechanism. The calculated Eads values exhibited notable ranges, spanning from –260.82 to –308.18 kcal/mol for SWCNT–PhCOOH and –220.92 to –261.01 kcal/mol for SWCNT–PhNH2 with the maximum probability values, representing the most favorable adsorption scenarios, determined to be –292.96 and –229.39 kcal/mol, respectively. A key insight from Monte Carlo simulations underscored the inherent spontaneity of the adsorption process, corroborated by the consistently negative Eads values. These findings collectively underscore the substantial affinity of the inhibitors to the copper surface, contributing to a deeper comprehension of their corrosion inhibition capabilities.

The inhibitor activity of some azo compound derivatives using density functional theory and molecular dynamics simulations

Corrosion poses a significant challenge for metals and alloys, carrying considerable economic and safety consequences. Various manufacturing processes, including surface preparation techniques, further exacerbate it. While corrosion inhibitors are employed to safeguard these materials from corrosion, they can also induce oxidative damage. This study aims to determine the corrosion inhibition efficiency of the specified materials. Density Functional Theory (DFT) and molecular dynamics simulations (MDS) were used to investigate the electronic characteristics of diverse compounds. The study included key parameters including highest occupied molecular orbital energy (EHOMO), lowest unoccupied molecular orbital energy (ELUMO), energy gap (Eg) between ELUMO and EHOMO, dipole moment, global hardness, softness (σ), ionization energy (I), electron affinity (A), electronegativity (χ), back-donation energy Eb-d, global electrophilicity (ω), electron transfer, global nucleophilicity (ε), and total energy (sum of electronic and zero-point energies). The calculations were executed utilizing the DFT method, employing a 6-311G++(d,p) basis set at the B3LYP level. Our research aimed to rank materials based on their effectiveness in inhibiting corrosion by considering EHOMO, ELUMO, and Eg. The order of effectiveness was found to be as follows: S4 > S3 > S1 > Mad1 > Mad3 > Mad4 > P-azo > Bad4. Additionally, our Investigation delved into alterations in electron affinity (A) and ionization potential (I) within different solvent phases, highlighting their importance in determining the chemical reactivity and stability of the studied compounds. The reactive sites were selected using Fukui indices, and the molecular electrostatic potential (MESP) surface was characterized—a Monte Carlo simulation to further understand these chemicals' adsorption behavior on Fe (110) surfaces. Our findings revealed that these compounds exhibit corrosion inhibition potential due to their high EHOMO, A, σ, ΔN, ΔE-back-donation, and low Eg, ELUMO, I, and η. Moreover, our MESP surface analysis demonstrated their ability to donate electrons to the metal's d-orbitals and receive electrons through back-donation. The compounds investigated in this study exhibit the potential to be employed as corrosion inhibitors in various industrial sectors. As a result, through the use and determination of quantum computational chemistry parameters, the level of corrosion resistance was as follows S3 > S1 > Bad4 > Mad4 > Mad3 > S4 > Mad1 > P-azo.

Valorizing the potential of saffron petals extract for aluminium corrosion control: An integrated approach involving extraction, experimental and computational analysis

Our study focuses on enhancing the antioxidant potential of saffron petals extract (SPE) as a sustainable corrosion inhibitor for aluminium (Al) in a 3.5% NaCl solution. The phytochemical screening and antioxidant activity are identified. Various investigative methodologies, including experimental techniques such as Electrochemical Impedance Spectroscopy (EIS) and Potentiodynamic Polarization (PDP), were employed to evaluate the effectiveness of this compound in inhibiting corrosion. Computational analysis, involving Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations, were utilized to carry out a Theoretical study on the main constituents of SPE. The surface characteristics and composition of the corroded Al were examined using a combination of Scanning Electron Microscopy (SEM), Energy-Dispersive X-ray Spectroscopy (EDX) and Fourier Transform Infrared Spectroscopy (FT-IR). The Electrochemical results suggest that the effectiveness of the inhibitor was dependent on its concentration, reaching 95.11% at 500 ppm of SPE. The PDP indicates that mixed inhibition control effectively slowed down the corrosion of Al. Furthermore, the adsorption behaviour of SPE onto Al follows to the Langmuir isotherm model. SEM, EDX and FT-IR analyses further verified significant alterations in the surface morphology and roughness of Al and confirming the successful formation of a protective layer on the Al surface. Ultraviolet-Visible spectroscopy (UV-Vis) reveals evidence of interaction between Al and the SPE molecules. Computational studies substantiated these observations, demonstrating the reactivity and adsorption patterns of the evaluated major bioactive molecules of SPE on the Al surface.

Experimental and theoretical insights into Artemisia Stems aqueous extract as a sustainable and eco-friendly corrosion inhibitor for mild steel in 1M HCl environment.

The present research demonstrates an innovative investigation of environmentally friendly mild steel (M-steel) corrosion inhibition using the artemisia stems aqueous extract (ASAEx) as an inhibitor in hydrochloric acid 1M. The standard extraction technique of hydrodistillation was used for producing the aqueous solutions of ASAEx. To assess the ratios of the chemical components, phytochemical screening was used to identify the stems of this plant. We used a variety of methods and techniques in our research on corrosion inhibition, including weight loss measures, surface analysis methods like XPS and SEM/EDS, electrochemical testing like PDP and EIS, as well as computational lead compound evaluation. Maximum inhibitory efficacy was achieved with 400mg/L ASAEx in 1M HCl at 303K, i.e. 90%. The PDP investigation verified the mixed-kind inhibitor status of the ASAEx extract. To describe the surface of M-steel, fitting and synthetic data were used to identify a constant phase element (CPE). SEM surface analysis was also used to detect the ASAEx effect on the surface of M-steel. X-ray photoelectron spectroscopy (XPS) analysis shows the presence of trace molecules of ASAEx on M-steel surface characterizing the bands in Maj-ASAEx (major compound of ASAEx). Density functional theory (DFT) and molecular dynamics simulations (MDs) were used in computational chemistry to clarify the adsorption mechanism and inhibitory impact.

Modulating Thermal Conductivity via Targeted Phonon Excitation.

Thermal conductivity is a critical material property in numerous applications, such as those related to thermoelectric devices and heat dissipation. Effectively modulating thermal conductivity has become a great concern in the field of heat conduction. Here, a quantum modulation strategy is proposed to modulate the thermal conductivity/heat flux by exciting targeted phonons. It shows that the thermal conductivity of graphene can be tailored in the range of 1559 W m-1 K-1 (decreased to 49%) to 4093 W m-1 K-1 (increased to 128%), compared with the intrinsic value of 3189 W m-1 K-1. The effects are also observed for graphene nanoribbons and bulk silicon. The results are obtained through both density functional theory calculations and molecular dynamics simulations. This novel modulation strategy may pave the way for quantum heat conduction.

Alkyl chain length affecting the structure and interactions of choline-based amino acid ionic liquids: Insights from a computational study

The influences of the alkyl chain length of amino acid anions on the structure and interactions of choline-based amino acid ([Cho][AA]) ionic liquids were investigated based on density functional theory (DFT) and molecular dynamics (MD) simulations. The most stable configuration of the individual [Cho][AA] ion pair was determined, revealing negligible variations in the binding energy between choline cation ([Cho]+) and different amino acid anions ([AA]−). To further characterize the intermolecular interactions between charged entities, the atoms in molecules theory (AIM) and independent gradient model based on Hirshfeld partition (IGMH) were employed, which highlight the presence of robust hydrogen bonds between [Cho]+ and [AA] −, underscoring their significance in shaping the structure of [Cho][AA] ion pairs. The MD studies exhibited the aggregation of anions in [Cho][AA] ILs, which was attributed to side chain interactions of anions. This aggregation effect becomes more pronounced with the elongation of anionic alkyl chain. Simultaneously, radial distribution functions (RDF) analysis revealed a higher probability of hydrogen bonds occurring between [Cho]+ and [AA]− as the alkyl chain in [AA]− lengthened. Additionally, the hydrogen bond lifetimes roughly display an increasing trend as the alkyl chain of the [AA]− increases, which may be due to the reduced diffusion coefficients of ions. Hopefully, our results would contribute fundamental knowledge regarding the impact of alkyl chain length of anions on the structure and interactions in the choline-based amino acid ionic liquids.

Understanding the Interface Enhancement Mechanisms of CFRP with the Polydopamine-Polyetheramine Interphase at the Molecular Level.

In this work, polydopamine (PDA) and polyetheramine D230 were selected to construct the PDA-D230 interphase between the carbon fiber (CF) and epoxy matrix. Density functional theory (DFT) and molecular dynamics (MD) simulations were performed to explore the interface enhancement mechanisms of a carbon fiber reinforced polymer (CFRP) with the PDA-D230 interphase from the molecular level. The adsorption characteristics of a PDA molecule on the CF surface were investigated using the DFT method. The results show that stronger π-π stacking interactions are formed due to the structure and orientation preference of the PDA molecule. The interfacial structures and properties of CFRP with the PDA-D230 interphase are derived from MD simulations. The PDA-D230 interphase on the CF surface induces stronger interfacial interaction energy, leading to the better load transfer between the CF and epoxy matrix. The existence of the PDA-D230 interphase on the CF surface can decrease the mean-square displacement (MSD) value and the free volume fraction of CFRP, which restricts the movement of epoxy atoms and inhibits the translational and rotational motion of epoxy chains. Compared with the epoxy using pristine CFs as reinforcement, the interfacial shear stress (ISS) of CFRP with the PDA-D230 interphase is improved by 13.1%. Our results provide valuable insights into the interface characteristics of CFRP with the PDA-D230 interphase, which are of great significance for exploring the strengthening mechanisms for CFRPs with the PDA-D230 interphase.

Unveiling a Novel Solvatomorphism of Anti-inflammatory Flufenamic Acid: X-ray Structure, Quantum Chemical, and In Silico Studies.

This paper delves into the polymorphism of 2-[3-(trifluoromethyl)anilino]benzoic acid, commonly referred to as flufenamic acid (FA), a pharmaceutical agent employed in treating inflammatory conditions. The central focus of the study is on a newly unearthed solvatomorphic structure of FA in methanol (FAM), and a thorough comparison is conducted with the commercially available standard structure. Employing a comprehensive approach, including X-ray crystallography, Hirshfeld surface analysis, density functional theory (DFT), molecular docking, and molecular dynamics (MD) simulations, the research aims to unravel the structural and functional implications of solvatomorphism. The X-ray crystal structure analysis brings to light notable differences between the standard FA and solvatomorphic FAM, showcasing variations in intermolecular interactions and crystal packing. Key features such as hydrogen bonding, π···π stacking, and C-H···π interactions are identified as influential factors shaping the stability and conformation of the compounds. Hirshfeld surface analysis further quantifies the nature and contribution of intermolecular interactions, providing a comprehensive perspective on molecular stability. Density functional theory offers valuable electronic structure insights, highlighting disparities in frontier molecular orbitals between FA and FAM. Molecular docking studies against prostaglandin D2 11-ketoreductase explore potential drug interactions, unveiling distinct binding modes and hydrogen bonding patterns that shed light on how the solvatomorphic structure may impact drug-target interactions. In-depth molecular dynamics simulations over 100 ns investigate the stability of the protein-ligand complex, with root mean square deviation and root mean square fluctuation analyses revealing minimal deviations and affirming the stability of FAM within the active site of the target protein.

Comprehensive analysis of a thiazole-substituted corrosion inhibitor's impact on N80 carbon steel in acidic conditions: Integrating computational predictions with experimental verifications

The present investigation elucidates the effectiveness of a newly developed organic inhibitor, namely 6,6-diphenyl-2,3-dihydroimidazo[2,1-b]thiazol-5(6H)-one (PHIT) in inhibiting the corrosion of N80 carbon steel (N80-CS) in 15 wt% HCl at 303 K. An integrative approach combining theoretical insights through Density Functional Theory (DFT) and molecular dynamics (MD) simulations with empirical data derived from a suite of analytical techniques (namely, weight loss measurements, electrochemical impedance spectroscopy (EIS), potentiodynamic polarization (PDP), and scanning electron microscopy (SEM)), was employed to assess the corrosion inhibition performance. The research identified the thiazolone fragment as the predominant site of activity within the PHIT molecule. According to EIS results, a significant reduction in the effective double-layer capacitance from 174 (blank) to 36.17 μF/cm2 was registered following the addition of 5 × 10−3 mol/L PHIT to the HCl medium. Moreover, at a low concentration of 10−4 mol/L, the inhibition efficiency reached a maximum of 95 %. The immersion time had a significant effect on the inhibitor's polarization resistance, peaking at 48 h (596 Ω cm2), then experiencing a minor reduction at 72 h (457 Ω cm2), in comparison to the blank measurement at 0.5 h (231 Ω cm2). Besides, the inhibitor demonstrated a high degree of effectiveness across various temperatures with a slight decrement in efficiency from 96.5 % at 303 K to 94 % at 363 K. PDP results confirmed the tested compound as an anodic-type inhibitor. These results affirmed the potential of PHIT as a highly efficacious organic inhibitor for N80-CS in acidic conditions, encouraging further development of its derivatives.

Exploring the synthesis and application of a pyrazole derivative in corrosion protection: Theoretical modeling and experimental investigations

This research investigates the inhibitory performance of BM-02, a synthesized pyrazole derivative (R)-5-(4-methoxyphenyl)-1,3-diphenyl-4,5-dihydro-1H-pyrazole, against the corrosion of mild steel in a 1 M hydrochloric acid solution. A comprehensive array of experimental techniques, including weight loss analysis, potentiodynamic polarization, electrochemical impedance spectroscopy, and scanning electron microscopy, were employed to assess the inhibitory effectiveness. BM-02 is identified as a mixed inhibitor, showcasing a remarkable efficiency of 91.5 % at a concentration of 10−3 mol/l. The adsorption behavior of the BM-02 molecule was investigated, revealing a physical adsorption model in accordance with the Langmuir isotherm. Furthermore, ATR-IR, 1H NMR, and 13C NMR spectra were acquired, complemented by density functional theory (DFT) calculations and molecular dynamics (MD) simulations to gain deeper insights into the molecular interactions and corrosion inhibition mechanism.

Experimental and computational exploration of the biodegradable platanus acerifolia leaf extract against mild steel corrosion in hydrochloric acid

This study delved into the potential application of Platanus acerifolia leaf extract (PALE) as an environmentally benign corrosion inhibitor. Evaluation of corrosion inhibition was conducted using electrochemical methods, field emission scanning electron microscopy equipped with energy dispersive spectrometer. The results reveal that PALE offers efficient protection for Q235 steel across the temperature range of 298 K–318 K, achieving a maximum corrosion inhibition efficiency of approximately 93%. Furthermore, Density functional theory (DFT) calculations and molecular dynamics (MD) simulations were performed to elucidate the spontaneous and stable adsorption of PALE molecules onto the surface of Q235 steel. This adsorption impedes the diffusion of corrosive substances towards the steel surface, thereby endowing Q235 steel with robust corrosion resistance.

Alloy [FA,Cs]PbI3 perovskite surfaces: The role of surface cesium composition in stability and tolerance to defect formation

Halide-perovskite alloys that include cesium have achieved records of stability and efficiency in solar cells. Controlling the surface composition, defects, and electronic properties guarantees interface stability and improves performance. By using density functional theory and molecular dynamic simulations, we analyzed which surface compositions of the formamidinium (FA) and cesium (Cs) lead iodide perovskite FA1−xCsxPbI3 with 25 and 50% of Cs become more stable than pure perovskites. Structural and electronic properties and tolerance to defect formation were also evaluated. Surface energy calculations show that only the alloys with 25% Cs and FAI-enriched surfaces are more stable than pure FAPbI3 ones. The most stable alloy surface shows electronic energy levels similar to the FAPbI3 perovskite, suggesting that this alloy may also be efficient for charge transport in the cell. However, the presence of Cs on the alloy surface, although low, favors the formation of FAI vacancies, which is detrimental to the stability of the perovskite. These results suggest evaluating FA1−xCsxPbI3 alloys with small Cs compositions to mitigate the formation of defects or using a passivation scheme. This study delivers valuable information for efficiency device improvement from the perspective of interface stability.

Exploring the adsorption of five emerging pollutants on activated carbon: A theoretical approach

The identification and management of contaminants of emerging concern (CECs) in water systems is crucial for protecting public and environmental health. This paper reports a theoretical approach to studying the adsorption of five CECs: Atrazine (ATZ), Caffeine (CAF), Carbamazepine (CBZ), Sulfamethoxazole (SMX), and Ibuprofen (IBU) - onto Activated Carbon (AC). A set of computational methods, including electrostatic molecular potential maps, conceptual density functional theory, Fukui functions, thermodynamic analysis, and tight-binding molecular dynamics simulations, were employed to analyze the electronic/energetic interactions and mechanisms involved in the adsorption of CECs on AC. The theoretical methodology offered valuable predictions on reactivity sites, stability, and binding mechanisms. Results showed that adsorption primarily occurred through non-covalent interactions like π-π electron donor-acceptor interactions, van der Waals forces, and hydrophobic interactions. Thermodynamic properties suggested the adsorption process was spontaneous and exothermic. However, for the AC/SMX system, the Gibbs free energy reveals that adsorption may be unfavorably compared to the other study systems. Molecular dynamics simulations validated the kinetic stability in the following order CAF (0.13 Å)>CBZ (0.23 Å)>ATZ (0.75 Å)> IBU (1.28 Å)>SMX (1.54 Å). This exploratory theoretical study provides a deep understanding of the interactions between AC and five CECs, aiding in the rational design and optimization of AC-based treatment systems for environmental and industrial applications.