van der Waals forces in density functional theory: a review of the vdW-DF method

Elucidation of the suppression of photoinduced segregation in 2D mixed halide, A2PbI2Br2: Critical role of A2PbBr4 photostability

2D lead halide perovskites are used to improve device operational stability due to increased environmental stability and reduced ion migration compared to 3D perovskites. However, the relationship between the 2D perovskite stability under illumination and spacer cation is still not well understood. Thus, we examine photoinduced halide segregation (PHS) in different 2D mixed-halide perovskites and show that PHS is suppressed in the materials which have photostable bromide halide phase. As the spacer cations provide the barrier to ion migration in 2D perovskites, PHS would be facilitated by the loss of spacer cations through their interactions with various mobile oxidized halide species, resulting in organic ammonium deprotonation and spacer cation vacancy formation. The existence of a photostable A2PbBr4 phase, which does not exhibit spacer cation loss, results in the suppression of PHS in 2D mixed halide perovskites due to reduced spacer vacancy formation and consequently reduced halide ion migration under illumination.

Structural Evolution Leading to the Thermosalient Phase Transition of Oxitropium Bromide

This study investigates the thermosalient effect in oxitropium bromide, with a focus on the role of anisotropic thermal expansion, elastic properties, and sound propagation in driving this phenomenon. Variable-temperature X-ray powder diffraction (VTXRPD) revealed significant anisotropic thermal expansion, including negative thermal expansion (NTE) along the c-axis in the low-temperature Form A. Density functional theory (DFT) calculations were used to analyze elastic properties of oxitropium bromide and confirmed that it does not exhibit negative compressibility, emphasizing thermal anisotropy as the primary factor in the phase transition. Studies of elastic constants and sound propagation demonstrated a preferred pathway for energy transfer along the z-direction, enabling rapid strain release during the phase transition. These findings confirmed that the thermosalient effect arises from cooperative molecular motion, resulting in an abrupt and energetic transformation driven by the interplay of structural anisotropy and elastic properties.

Exploring atom-pairwise and many-body dispersion corrections for the BEEF-vdW functional

The Bayesian error estimation functional (BEEF-vdW) is widely used in surface science and catalysis, because it provides a balanced description of molecular, surface, and solid state systems, along with reliable error estimates. However, the nonlocal van-der-Waals density functional (vdW-DF2) employed in BEEF-vdW can be computationally costly and displays relatively low accuracy for molecular systems. Therefore, this work explores whether atom-pairwise and many-body dispersion treatments represent viable alternatives to using the vdW-DF2 functional with BEEF-vdW. To this end, we investigate the performance of commonly used atom-pairwise corrections [i.e., the Tkatchenko-Scheffler (TS) and the exchange-hole dipole moment (XDM) approaches] and many-body dispersion (MBD) treatments for molecular, surface, and solid-state systems. The results indicate that atom-pairwise methods such as TS and particularly XDM provide a good balance of cost and accuracy across all systems.

Many-Body Effects at Heterogeneous Interfaces from First-Principles: Progress, Challenges, and Opportunities

Heterogeneous interfaces are pivotal in numerous nanoscale devices and applications. First-principles approaches based on quantum mechanics and atomistic structures provide critical insights into structure-property relationships, enabling the informed design of materials and devices. Accurate first-principles methods must reliably capture many-body effects, i.e., electron-electron interactions, which significantly influence system properties compared to the predictions from models using free or noninteracting electrons. In this Perspective, we survey a few computational tools in this context and attempt to be forward looking by discussing the current challenges and emerging research opportunities. We examine diverse manifestations of many-body effects across several domains: geometries and reaction barriers (total-energy properties), orbital energy levels and band alignments (one-particle properties), and optical excited states (two-particle properties).

Non-Covalent Self-Assembly Behaviors Based on Racemic Binaphthol Scaffolds

Understanding the fundamental mechanisms involved in the construction and organization of multi-scale structures is crucial for the design and manufacture of complex functional systems with long-range molecular arrangements. In this paper, a series of compounds have been synthesized using racemic binaphthols as the skeleton and a Suzuki coupling reaction for derivatization at the 6,6' positions, which resulted in various structures bearing different functional groups. Control over the self-assembly of these racemic binaphthol derivatives was successfully achieved by adjusting the types and positions of the substituents in the parent binaphthol compound, which revealed the key factors influencing the types of the non-covalent interactions and the self-assembly process. For example, the single-crystal structures of the resulting compounds indicated that assembly structures such as single helix and double helix based on non-traditional hydrogen bond motifs could be obtained, and fascinating non-covalent self-assembly structures such as molecular ladders and catenane discovered.

Machine Learning Recognition of Artificial DNA Sequence with Quantum Tunneling Nanogap Junction

Artificially synthesized DNA holds significant promise in addressing fundamental biochemical questions and driving advancements in biotechnology, genetics, and DNA digital data storage. Rapid and precise electric identification of these artificial DNA strands is crucial for their effective application. Herein, we present a comprehensive investigation into the electric recognition of eight artificial synthesized DNA (xDNA and yDNA) nucleobases using quantum tunneling transport and machine learning (ML) techniques. By embedding these nucleobases within a solid-state nanogap junction, we calculated their fingerprint transmission and current readouts and also analyzed the influence of electronic coupling and molecular orbital delocalization on these properties. The trained ML model achieved a predictive basecalling accuracy of up to 100% for xDNA nucleobases and 99.80% for yDNA transmission readout data sets. ML explainability study revealed that normalized descriptors have a greater impact on nucleobase prediction than the original transmission function, proving more effective in disentangling overlapping artificial DNA nucleobase signals. Quaternary classification results highlighted higher recognition accuracy for xDNA nucleobases than for yDNA nucleobases. Furthermore, precise calling of complementary, purine, and pyrimidine base pair combinations was demonstrated with high sensitivity and an F1 score. Our findings reveal the feasibility of highly sensitive and precise electrical recognition of artificial DNA nucleobases, which can transform genetic research and spur advancements in genetic data storage, synthetic biology, and diagnostics.

A General Nonbonded Force Field Based on Accurate Quantum Mechanics Calculations for Elements H-La and Hf-Rn

Noncovalent interactions (NCI) play a central role in numerous physical, chemical, and biological phenomena. An accurate description of NCI is the key to success for any theoretical study in such areas. Although quantum mechanics (QM) methods such as dispersion-corrected density functional theory are sufficiently accurate, their applications are practical only for <300 atoms and <100 ps of simulation time. Thus, empirical force fields (FF) have generally been the only choice for systems with thousands to millions of atoms and for nanoseconds and longer. We want to develop a FF that can be applied to applications of thousands to millions of atoms with an accuracy comparable to QM methods. As the first step, we develop here a new general nonbonded potential (GNB) based on a novel functional form with four adjustable parameters for each element. We report here parameters for elements H-La, Hf-Rn (excluding lanthanides and actinides) by fitting the interaction energy of molecular complexes to QM calculations using the accurate Head-Gordon ωB97M-V density functional. We performed extensive testing of GNB for organic molecules, organometallic molecules, and metal organic-framework (MOF) systems. The mean absolute errors of GNB are 0.37 kcal/mol for the dispersion and mixed groups of the S66 × 8 benchmark set, 0.35 kcal/mol for CO2 adsorption on MOF materials, and 4.53 kcal/mol for the XTMC43 benchmark. GNB outperforms existing FF and in many cases has accuracy comparable to that of QM methods such as PBE-D3. GNB can potentially replace the nonbonded part of existing FFs in a wide range of applications.

Dynamical Disorder in the Mesophase Ferroelectric HdabcoClO4: A Machine-Learned Force Field Study

Hybrid molecular ferroelectrics with orientationally disordered mesophases offer significant promise as lead-free alternatives to traditional inorganic ferroelectrics owing to properties such as room temperature ferroelectricity, low-energy synthesis, malleability, and potential for multiaxial polarization. The ferroelectric molecular salt HdabcoClO4 is of particular interest due to its ultrafast ferroelectric room-temperature switching. However, so far, there is limited understanding of the nature of dynamical disorder arising in these compounds. Here, we employ the neural network NeuralIL to train a machine-learned force field (MLFF) with training data generated using density functional theory. The resulting MLFF-MD simulations exhibit phase transitions and thermal expansion in line with earlier reported experimental results, for both a low-temperature phase transition coinciding with the orientational disorder of ClO4 - and the onset of rotation of both Hdabco+ and ClO4 - in a high-temperature phase transition. We also find proton transfer even in the low-temperature phase, which increases with temperature and leads to associated proton disorder as well as the onset of disorder in the direction of the hydrogen-bonded chains.

Stabilization of small organic molecules on V2C and V2CO2 MXenes: first-principles insights into the performance of van der Waals functionals and the effect of oxygen vacancies

The adsorption of small organic molecules on pristine V2C MXene and its derivatives is investigated by first-principles density functional theory calculations. By employing state-of-the-art van der Waals (vdW) density functionals, the binding affinity of studied molecules, i.e., CH4, CO2, and H2O on MXene adsorbents is well described by more recent vdW functionals, i.e., SCAN-rvv10. Although both CH4 and CO2 are nonpolar molecules, on pristine and oxygen-vacancy surfaces, they show a different range of adsorption energies, in which CH4 is more inert and has weaker binding than CO2. CO2 stays intact in its molecular forms for most of the tested functionals, except for the case of the vdW-DF functional, where CO2 exhibits a dissociation regardless of its initial adsorption geometry. For full surface terminations, the adsorption affinity of all involved species is comparable within the same range, varying from -0.10 to -0.20 eV, attributed to either weak dispersion interactions or hydrogen bonds. The binding of H2O is much more pronounced compared to CO2 and CH4 in the presence of oxygen vacancies with the highest adsorption energy of -1.33 eV, vs. -0.67, and -0.20 eV obtained for H2O, CO2, and CH4 respectively. H2O can dissociate with a small activation energy barrier of 0.40 eV, much smaller than its molecular adsorption energy, to further saturate itself on the surface. At high oxygen-vacancy concentrations, stronger bindings of adsorbates are found due to a preferred attachment of adsorbates to induced undercoordinated metal sites. The findings propose a potential scheme for greenhouse gas separation based on the surface modification of novel two-dimensional structures.

The Electronic Properties and Adsorption Performance of LDH/Graphene, and LDH/g-C3N4 for the Removal of Pharmaceutical Contaminants: A Molecular Dynamics Simulation

Water shortages and pharmaceutical pollution are two interconnected crises that pose severe threats to global health, environmental sustainability, and economic stability. Pharmaceutical pollution is widespread and has reached potentially toxic levels in over 258 rivers in 104 countries. So far, more interest has been paid towards efficient water treatment processes in recent years. In this study, we explore the efficacy of layered double hydroxide (LDH) nanocomposites with graphene and graphitic carbon nitride (g-C3N4) as promising adsorbents of pharmaceutical contaminants. The LDH nanocomposite has been designed and simulated for the first time, consisting of two layers of sodium hydroxide with a layer of graphene and g-C3N4. We investigated the adsorption performance of LDH, specifically LDH/graphene and LDH/g-C3N4, for the removal of pharmaceutical contaminants including acetaminophen (AC), caffeine (CAF), and sulfamethoxazole (SMZ). Through comprehensive molecular dynamics simulations using the reactive forcefield (ReaxFF) software, we investigated the adsorption mechanisms, kinetics, and adsorption capacity of pharmaceutical contaminants onto these nanocomposite surfaces. Our findings showed that the combination of LDH/graphene had a higher adsorption capacity for the removal of pharmaceutical contaminants than LDH/g-C3N4. At 70 Picoseconds (Ps), 124, 129, and 142 molecules of each of the pharmaceutical contaminants AC, CAF and SMZ, respectively, had been adsorbed by LDH/graphene, with a higher exothermic energy equating to -1111, -1015, and -1150 × 103 kJ/mol, respectively. On the other hand, for LDH/g-C3N4 at 70 Ps, 108, 110, and 120 molecules of AC, CAF and SMZ, respectively, had been adsorbed, with exothermic energy equating to -978, -948, and -1173 × 103 kJ/mol, respectively. Finally, we calculated the electronic properties, including the band gap and density of state of the nanocomposite materials, to check their effect on the adsorption process. In addition, the results showed that the adsorption kinetics followed a pseudo-first-order model, while the adsorption isotherms for AC, CAF and SMZ adhered to the Langmuir model.

Advancements and Challenges in the Integration of Indium Arsenide and Van der Waals Heterostructures

The strategic integration of low-dimensional InAs-based materials and emerging van der Waals systems is advancing in various scientific fields, including electronics, optics, and magnetics. With their unique properties, these InAs-based van der Waals materials and devices promise further miniaturization of semiconductor devices in line with Moore's Law. However, progress in this area lags behind other 2D materials like graphene and boron nitride. Challenges include synthesizing pure crystalline phase InAs nanostructures and single-atomic-layer 2D InAs films, both vital for advanced van der Waals heterostructures. Also, diverse surface state effects on InAs-based van der Waals devices complicate their performance evaluation. This review discusses the experimental advances in the van der Waals epitaxy of InAs-based materials and the working principles of InAs-based van der Waals devices. Theoretical achievements in understanding and guiding the design of InAs-based van der Waals systems are highlighted. Focusing on advancing novel selective area growth and remote epitaxy, exploring multi-functional applications, and incorporating deep learning into first-principles calculations are proposed. These initiatives aim to overcome existing bottlenecks and accelerate transformative advancements in integrating InAs and van der Waals heterostructures.

Layer-number and strain effects on the structural and electronic properties of PtSe2material

Bandgap engineering of low-dimensional materials forms a robust basis for advancements in optoelectronic technologies. Platinum diselenide (PtSe2) material exhibits a transition from semi-metal to semiconductor (SM-SC) when going from bulk to monolayer. In this work, density functional theory (DFT) with various van der Waals (vdW) corrections has been tested to study the effect of the layer-number on the structural and electronic properties of the PtSe2material. The considered vdW corrections gave different results regarding the number of layers at which the SM-SC transition occurs. This variation is due to the different interlayer distances found for each correction, revealing the sensitivity of the bandgap to this distance in addition to the layer number. In fact, the bandgap increases with the increasing of the interlayer distance, due to the energy shift of conduction and valence bands dominated by Se-pzorbitals. According to the comparison with the available experimental data, the vdW corrections vdW-DF and rVV10 gave the most accurate results. Moreover, the control of the interlayer distance via vertical compressive strain led to the bandgap tuning of semiconductor PtSe2BL. Indeed, a semi-metal character of PtSe2BL can be obtained under 17% vertical strain. Our work shows a deep understanding of the correlation between the structural and electronic properties, and thus a possibility to tune the bandgap by strain means.

Doxorubicin-Intercalated Li-Al-Based LDHs as Potential Drug Delivery Nanovehicle with pH-Responsive Therapeutic Cargo for Tumor Treatment

Clinical oncology is currently experiencing a technology bottleneck due to the expeditious evolution of therapy defiance in tumors. Although drugs used in chemotherapy work for a sort of cell death with potential clinical application, the effectiveness of chemotherapy-inducing drugs is subject to several endogenous conditions when used alone, necessitating the urgent need for controlled mechanisms. A tumor-targeted drug delivery therapy using Li-Al (M+/M3+)-based layered double hydroxide (LDHs) family has been proposed with the general chemical formula [M+1-x M3+x (OH)]2x+[(Am-)2x/m. n(H2O)]2x-, which is fully biodegradable and works in connection with the therapeutic interaction between LDH nanocarriers and anticancerous doxorubicin (DOX). Compositional variation of Li and Al in LDHs has been used as a nanoplatform, which provides a functional balance between circulation lifetime, drug loading capacity, encapsulation efficiency, and tumor-specific uptake to act as self-regulatory therapeutic cargo to be released intracellularly. First-principle analyses based on DFT have been employed to investigate the interaction of bonding and electronic structure of LDH with DOX and assess its capability and potential for a superior drug carrier. Following the internalization into cancer cells, nanoformulations are carried to the nucleus via lysosomes, and the mechanistic pathways have been revealed. Additionally, in vitro along with in vivo therapeutic assessments on melanoma-bearing mice show a dimensional effect of nanoformulation for better biocompatibility and excellent synergetic anticancer activity. Further, the severe toxic consequences associated with traditional chemotherapy have been eradicated by using injectable hydrogel placed just beneath the tumor site, and regulated release of the drug has been confirmed through protein expression applying various markers. However, Li-Al-based LDH nanocarriers open up new design options for multifunctional nanomedicine, which has intriguing potential for use in cancer treatment through sustained drug delivery.

Slip Opacity and Fast Osmotic Transport of Hydrophobes at Aqueous Interfaces with Two-Dimensional Materials

We investigate the interfacial transport of water and hydrophobic solutes on van der Waals bilayers and heterostructures formed by stacking graphene, hBN, and MoS2 using extensive ab initio molecular dynamics simulations. We compute water slippage and the diffusio-osmotic transport coefficient of hydrophobic particles at the interface by combining hydrodynamics and the theory of the hydrophobic effect. We find that slippage is dominated by the layer that is in direct contact with water and only marginally altered by the second layer, leading to a so-called "slip opacity". The screening of the lateral forces, where the liquid does not feel the forces coming from the second nearest layer, is one of the factors leading to the "slip opacity" in our systems. The diffusio-osmotic transport of small hydrophobes (with a radius below 2.5 Å) is also affected by the slip opacity, being dramatically enhanced by slippage. Furthermore, the direction of diffusio-osmotic flow is controlled by the solute size, with the flow in the opposite direction of the concentration gradient for smaller hydrophobes, and vice versa for larger ones. We connect our findings to the wetting properties of two-dimensional materials, and we propose that slippage and wetting can be controlled separately: whereas the slippage is mostly determined by the layer in closer proximity to water, wetting can be finely tuned by stacking different two-dimensional materials. Our study advances the computational design of two-dimensional materials and van der Waals heterostructures, enabling precise control over wetting and slippage properties for applications in coatings and water purification membranes.

Putting error bars on density functional theory

Predicting the error in density functional theory (DFT) calculations due to the choice of exchange-correlation (XC) functional is crucial to the success of DFT, but currently, there are limited options to estimate this a priori. This is particularly important for high-throughput screening of new materials. In this work, the structure and elastic properties of binary and ternary oxides are computed using four XC functionals: LDA, PBE-GGA, PBEsol, and vdW-DF with C09 exchange. To analyze the systemic errors inherent to each XC functional, we employed materials informatics methods to predict the expected errors. The predicted errors were also used to better the DFT-predicted lattice parameters. Our results emphasize the link between the computed errors and the electron density and hybridization errors of a functional. In essence, these results provide "error bars" for choosing a functional for the creation of high-accuracy, high-throughput datasets as well as avenues for the development of XC functionals with enhanced performance, thereby enabling the accelerated discovery and design of new materials.

Theoretical and Experimental Studies of Molecular Interactions between Engineered Graphene and Phosphate Ions for Graphene-Based Phosphate Sensing

Fundamental understanding of the interactions of nanoscale materials with molecules of interest is essential for the development of electronic devices, such as sensors. In particular, structures and molecular interaction properties of engineered graphenes are still largely unexplored, despite these materials' great potential to be used as molecular sensors. As an example of end user application, the detection of phosphorus in the form of phosphate in a soil environment is important for soil fertility and plant growth. However, due to the lack of an affordable technology, it is currently hard to measure the amount of phosphate directly in the soil; therefore, suitable sensor technologies need to be developed for phosphate sensors. In this work, pristine graphene and several modified graphene materials (oxygenated graphene, graphene with vacancies, and curved graphene) were studied as candidates for phosphate sensor materials using density functional theory (DFT) calculations. Our calculations showed that both pristine graphene and functionalized graphene were able to adsorb phosphate species strongly. In addition, these graphene nanomaterials showed selectivity of adsorption of phosphate with respect to nitrate, with stronger adsorption energies for phosphate. Furthermore, our calculations showed significant changes in electrical conductivities of pristine graphene and functionalized graphenes after phosphate species adsorption, in particular, on graphene with oxygen (hydroxyl and epoxide) functional groups. Experimental measurements of electrical resistivity of graphene before and after adsorption of dihydrogen phosphate showed an increase in resistivity upon adsorption of phosphate, consistent with the theoretical predictions. Our results recommend graphene and functionalized graphene-based nanomaterials as good candidates for the development of phosphate sensors.

Theoretical investigation of 2D/2D van der Waals SbPO4/BiOClxBr1-x heterojunctions for photocatalytic water splitting

Bismuth halogenoxide (BiOX)-based heterojunctions have garnered considerable attention recently due to their potential to enhance photocatalytic performance. However, the predominant focus on II-type heterojunctions has posed challenges in achieving the requisite band edge positions for efficient water splitting. In this investigation, stable van der Waals SbPO4/BiOClxBr1-x heterojunctions were constructed theoretically by using density-functional theory (DFT). Our findings demonstrate that SbPO4 can modulate the formation of Z-scheme heterojunctions with BiOClxBr1-x. The structural properties of BiOX were preserved, while reaching excellent photocatalytic capabilities with high redox capacities. Further investigation unveiled that the band edge positions of the heterojunctions fully satisfy the oxidation-reduction potential of water. Moreover, these heterojunctions exhibit notable absorption efficiency in the visible range, with absorption increasing as x decreases. Our research provides valuable theoretical insights for the experimental synthesis of high-performance BiOX-based photocatalysts for water splitting, leveraging the unique properties of SbPO4. These insights contribute to the advancement of clean energy technology.

Atomic-precision control of plasmon-induced single-molecule switching in a metal-semiconductor nanojunction

Atomic-scale control of photochemistry facilitates extreme miniaturisation of optoelectronic devices. Localised surface plasmons, which provide strong confinement and enhancement of electromagnetic fields at the nanoscale, secure a route to achieve sub-nanoscale reaction control. Such local plasmon-induced photochemistry has been realised only in metallic structures so far. Here we demonstrate controlled plasmon-induced single-molecule switching of peryleneanhydride on a silicon surface. Using a plasmon-resonant tip in low-temperature scanning tunnelling microscopy, we can selectively induce the dissociation of the O-Si bonds between the molecule and surface, resulting in reversible switching between two configurations within the nanojunction. The switching rate can be controlled by changing the tip height with 0.1-Å precision. Furthermore, the plasmon-induced reactivity can be modified by chemical substitution within the molecule, suggesting the importance of atomic-level design for plasmon-driven optoelectronic devices. Thus, metal-single-molecule-semiconductor junctions may serve as a prominent controllable platform beyond conventional nano-optoelectronics.

London Dispersion Effects on the Structure and Properties of Nonlinear Optical BiB3O6 Crystal

α-BiB3O6 (BiBO) is an important nonlinear optical (NLO) material with high efficiency for applications in harmonic generations and quantum technology. Owing to its low symmetry and cooperative Bi3+ lone pair arrangement, it has also exceptional large piezoelectric and electro-optic coefficients and strong anisotropies on other material characteristics. Previous theoretical calculations on its physical (mainly optical) properties often gave confusing results. It is found here that London dispersion (LD) tends to stabilize structures with closer pack entities like lone pair heavy ion Bi3+ with large polarizabilities, which is ignored in most previous density functional theory (DFT) calculations. Present study shows that without considering the LD effect, the structure of α-BiB3O6 (BiBO) was predicted with an over-estimated (by over 10 %) unique b-axis while underestimates a and overestimates c in a less amount. Consequently it is not possible to use the calculated structure to obtain meaningful properties of this important material. By applying a modified post-DFT LD correction based on linear combination of atomic orbitals (LCAO) and B3LYP functional, the experimental structure is well reproduced with the theoretical optimized one. Many important material property tensors of BiBO crystal are calculated in unprecedented precisions, including: dielectric constants (static and in THz range), elastic and elasto-optic constants, piezoelectric constants, refractive indices, NLO and electro-optic (EO) coefficients. Among them, theoretical calculation of the refractive indices in the THz range by diagonalizing the clamped-ion dielectric constants was firstly achieved at least for BiBO crystal. The calculation also confirms that BiBO has an exceptional large piezoelectric constant d22=40 pC/N and largest free EO coefficientsγ 12 T ${{\gamma }_{12}^{T}}$ ,γ 22 T ${{\gamma }_{22}^{T}}$ ,γ 32 T ${{\gamma }_{32}^{T}}$ on the order of 10 pm/V among borate crystals. The calculation also reveals that the large free EO coefficients are largely originated from the piezoelectric induced photo-elastic effect and for practical high speed applications only the clamped-ion EO coefficients take effect. The clamped ion EO coefficient ofγ 53 S ${{\gamma }_{53}^{S}}$ =-4.17 pm/V,γ b 1 S ${{\gamma }_{{\rm b}1}^{S}}$ =-2.61 pm/V are obtained for the first time and may be consulted if one seeks to design BiBO crystal as a high-speed EO modulator. Furthermore, full tensor matrix of the elasto-optic constants was obtained on the first time. Together with the calculated elastic constants, it can help to design acoustic optic modulating devices with preferable figure of merits 10 times that of traditional quartz crystal.

Synthesis and Polymorph Manipulation of FeSe2 Monolayers

Polymorph engineering involves the manipulation of material properties through controlled structural modification and is a candidate technique for creating unique two-dimensional transition metal dichalcogenide (TMDC) nanodevices. Despite its promise, polymorph engineering of magnetic TMDC monolayers has not yet been demonstrated. Here we grow FeSe2 monolayers via molecular beam epitaxy and find that they have great promise for magnetic polymorph engineering. Using scanning tunneling microscopy (STM) and spectroscopy (STS), we find that FeSe2 monolayers predominantly display a 1T' structural polymorph at 5 K. Application of voltage pulses from an STM tip causes a local, reversible transition from the 1T' phase to the 1T phase. Density functional theory calculations suggest that this single-layer structural phase transition is accompanied by a magnetic transition from an antiferromagnetic to a ferromagnetic configuration. These results open new possibilities for creating functional magnetic devices with TMDC monolayers via polymorph engineering.

Ultralow Detection of Mancozeb Using Two-Dimensional Cobalt Telluride (CoTe2)

Pesticides are crucial in modern agriculture because they reduce pests and boost yield, but they also represent major risks to human health and the environment; therefore, it is important to monitor their presence in food. Reliable and precise detection techniques are possible ways to address this issue. In this work, we utilize atomically thin (two-dimensional) cobalt telluride (CoTe2) with a high surface area and charge as a template material to detect mancozeb using spectroscopic and electrochemical techniques. When mancozeb (MNZ) molecules interact with 2D CoTe2, spectroscopic analyses reveal distinctive spectral shifts that clarify the underlying chemical interactions and binding mechanisms. Furthermore, CoTe2's electroactive sites and their manipulation for improved sensitivity and selectivity toward certain MNZ molecules are investigated by electrochemical studies. The CoTe2/GCE electrode exhibits enhanced electrochemical activity toward the electrooxidation of MNZ. The developed sensing electrode shows a linear range from 0.184 mM to 18.48 μM and a limit of detection of about 0.18 μM. In addition, we employ density functional theory (DFT) first-principles calculations to validate the experimental findings and comprehend the mechanism behind the interaction between CoTe2 and MNZ molecules. The study highlights the effectiveness of 2D CoTe2 as a dual-mode sensing platform for qualitative and quantitative assessment of MNZ pollutants, demonstrated by the integration of electrochemistry and spectroscopy and the critical role that 2D CoTe2-based sensors can play in accurate and efficient pesticide detection, which is required for agricultural safety protocols and environmental monitoring.

Control of proton transport and hydrogenation in double-gated graphene

The basal plane of graphene can function as a selective barrier that is permeable to protons1,2 but impermeable to all ions3,4 and gases5,6, stimulating its use in applications such as membranes1,2,7,8, catalysis9,10 and isotope separation11,12. Protons can chemically adsorb on graphene and hydrogenate it13,14, inducing a conductor-insulator transition that has been explored intensively in graphene electronic devices13-17. However, both processes face energy barriers1,12,18 and various strategies have been proposed to accelerate proton transport, for example by introducing vacancies4,7,8, incorporating catalytic metals1,19 or chemically functionalizing the lattice18,20. But these techniques can compromise other properties, such as ion selectivity21,22 or mechanical stability23. Here we show that independent control of the electric field, E, at around 1 V nm-1, and charge-carrier density, n, at around 1 × 1014 cm-2, in double-gated graphene allows the decoupling of proton transport from lattice hydrogenation and can thereby accelerate proton transport such that it approaches the limiting electrolyte current for our devices. Proton transport and hydrogenation can be driven selectively with precision and robustness, enabling proton-based logic and memory graphene devices that have on-off ratios spanning orders of magnitude. Our results show that field effects can accelerate and decouple electrochemical processes in double-gated 2D crystals and demonstrate the possibility of mapping such processes as a function of E and n, which is a new technique for the study of 2D electrode-electrolyte interfaces.

Molecular Mechanistic Insights into Dipeptidyl Peptidase-IV Inhibitory Peptides to Decipher the Structural Basis of Activity

Dipeptidyl peptidase-IV (DPP-IV) inhibiting peptides have attracted increased attention because of their possible beneficial effects on glycemic homeostasis. However, the structural basis underpinning their activities has not been well understood. This study combined computational and in vitro investigations to explore the structural basis of DPP-IV inhibitory peptides. We first superimposed the Xaa-Pro-type peptide-like structures from several crystal structures of DPP-IV ligand-protein complexes to analyze the recognition interactions of DPP-IV to peptides. Thereafter, a small set of Xaa-Pro-type peptides was designed to explore the effect of key interactions on inhibitory activity. The intramolecular interaction of Xaa-Pro-type peptides at the first and third positions from the N-terminus was pivotal to their inhibitory activities. Residue interactions between DPP-IV and residues of the peptides at the fourth and fifth positions of the N-terminus contributed significantly to the inhibitory effect of Xaa-Pro-type tetrapeptides and pentapeptides. Based on the interaction descriptors, quantitative structure-activity relationship (QSAR) studies with the DPP-IV inhibitory peptides resulted in valid models with high R2 values (0.90 for tripeptides; 0.91 for tetrapeptides and pentapeptides) and Q2 values (0.33 for tripeptides; 0.68 for tetrapeptides and pentapeptides). Taken together, the structural information on DPP-IV and peptides in this study facilitated the development of novel DPP-IV inhibitory peptides.

Meta-Substituted Asymmetric Azobenzenes: Insights into Structure-Property Relationship

This article presents a comprehensive investigation into the functionalization of methoxyphenylazobenzene using electron-directing groups located at the meta position relative to the azo group. Spectroscopic analysis of meta-functionalized azobenzenes reveals that the incorporation of electron-withdrawing units significantly influences the absorption spectra of both E and Z isomers, while electron-donating functionalities lead to more subtle changes. The thermal relaxation process from Z to E result in almost twice as prolonged for electron-withdrawing functionalized azobenzenes compared to their electron-rich counterparts. Computational analysis contributes a theoretical understanding of the electronic structure and properties of meta-substituted azobenzenes. This combined approach, integrating experimental and computational techniques, yields significant insights into the structure-property relationship of meta-substituted asymmetrical phenolazobenzenes.

Predicting two-dimensional semiconductors using conductivity effective mass

In this paper we investigate the relationship between the conductivity effective mass and exfoliation energy of materials to assess whether automatic sampling of the electron band structure can predict the presence of and ease of separating chemically bonded layers. We assess 22 976 materials from the Materials Project database, screen for only those that are thermodynamically stable and identify the 1000 materials with the highest standard deviation for p-type and the 1000 materials with the highest standard deviation for n-type internal conductivity effective mass tensors. We calculate the exfoliation energy of these 2000 materials and report on the correlation between effective mass and exfoliation energy. A relationship is found which is used to identify a previously unconsidered two-dimensional material and could streamline the modelling of other two-dimensional materials in the future.

Morphology and selectivity of hydrated alkali metal ions as depth of discharge in the 1T-MoS2 electrode with aqueous electrolytes

Aqueous ion batteries have great commercial potential in green power and energy storage due to their green nature, safety and high ionic conductivities. Different from organic electrolytes, alkali ions (Li+, Na+, and K+) inevitably bring water molecules into the electrodes during the charging/discharging process due to the hydration of ions with water molecules. The selectivity of alkali ions and the mechanism of how water molecules are involved in the ion extraction/insertion process in the electrodes have not been clarified. In this study, we focus on the characteristics of the intra-layer distribution of different hydrated ions (Li+, Na+, and K+) and the quantitative analysis of the selectivity of hydrated cations in aqueous batteries. We found that the concentration of hydrated ions greatly affects their distribution within the 1T-MoS2 layers, and the presence of hydrogen bonding and O-O repulsive forces between water molecules causes the hydrated ions to gradually form chains from the dispersed state under the effect of hydrogen bonding and ionic bonding, then further form strips, and ultimately be densely dispersed within the whole layer. In addition, the chemical potential difference of hydrated ions is the key to the competitive reaction, and we quantitatively analyze the selectivity relationship between hydrated cations throughout the charging and discharging process; hydrated sodium ions will have better performance than lithium and potassium ions in aqueous batteries.

Elucidating the Non-Covalent Interactions that Trigger Interdigitation in Lead-Halide Layered Hybrid Perovskites

Two-dimensional (2D) hybrid organic-inorganic perovskites constitute a versatile class of materials applied to a variety of optoelectronic devices. These materials are composed of alternating layers of inorganic lead halide octahedra and organic ammonium cations. Most perovskite research studies so far have focused on organic sublattices based on phenethylammonium and alkylammonium cations, which are packed by van der Waals cohesive forces. Here, we report a more complex organic sublattice containing benzotriazole-based ammonium cations packed through interdigitated π-π stacking and hydrogen bonding. Single crystals and thin films of four perovskite derivatives are studied in depth with optical spectroscopy and X-ray diffraction, supported by density-functional theory calculations. We quantify the lattice stabilization of interdigitation, dipole-dipole interactions, and inter- as well as intramolecular hydrogen bonding. Furthermore, we investigate the driving force behind interdigitation by defining a steric occupancy factor σ and tuning the composition of the organic and inorganic sublattice. We relate the phenomenon of interdigitation to the available lattice space and to weakened hydrogen bonding to the inorganic octahedra. Finally, we find that the stabilizing interactions in the organic sublattice slightly improve the thermal stability of the perovskite. This work sheds light on the design rules and structure-property relationships of 2D layered hybrid perovskites.

Temperature-Dependent Frictional Behavior of MoS2 in Humid Environments: Insights from Water Molecule Adsorption and DFT Analyses

This study investigates the temperature-dependent frictional behavior of MoS2 in humid environments within the context of a long-standing debate over increased friction due to oxidation processes or molecular adsorption. By combining sliding friction experiments and density functional theory (DFT)-based first-principles simulations, it aims to clarify the role of water molecule adsorption in influencing frictional properties of MoS2, addressing the challenge of identifying interfacial bonding behavior responsible for friction in such conditions. Sliding experiments revealed that magnetron-sputtered MoS2 exhibits a reduction in the coefficient of friction (COF) with an increase in temperature from 25 to 100 °C under 20 and 40% relative humidity. This change in the COF obeys the Arrhenius law, presenting an energy barrier of 0.165 eV, indicative of the temperature-dependent nature of these frictional changes and suggests a consistent frictional mechanism. DFT simulations showed that H2O molecules are adsorbed at MoS2 vacancy defects with adsorption energies ranging from -0.56 to -0.17 eV, which align with the experimentally determined energy barrier. Adsorptive interactions, particularly the formation of stable H···S interfacial hydrogen bonds at defect sites, increase the interlayer adhesion and impede layer shearing. TEM analysis confirms that although MoS2 layers align parallel to the sliding direction in humid conditions, the COF remains at 0.12, as opposed to approximately 0.02 in dry air. This demonstrates that parallel layer alignment does not notably decrease the COF, underscoring humidity's significant role in MoS2's COF values, a result also supported by the Arrhenius analysis. The reversibility of the physisorption process, demonstrated by the recovery of the COF in high-temperature humid environments, suggests its dynamic nature. This study yields fundamental insights into MoS2 interfaces for environments with variable humidity and temperature, crucial for demanding tribological applications.

Formic acid stability in different solvents by DFT calculations

CONTEXT

New technologies have been developed toward the use of green energies. The production of formic acid (FA) from carbon dioxide (CO[Formula: see text]) hydrogenation with H[Formula: see text] is a sustainable process for H[Formula: see text] storage. However, the FA adduct stabilization is thermodynamically dependent on the type of solvent and thermodynamic conditions. The results suggest a wide range of dielectric permittivity values between the dimethyl sulfoxide (DMSO) and water solvents to stabilize the FA in the absence of base. The thermodynamics analysis and the infrared and charge density difference results show that the formation of the FA complex with H[Formula: see text]O is temperature dependent and has a major influence on aqueous solvents compared to the FA adduct with amine, in good agreement with the experiment. In these conditions, the stability thermodynamic of the FA molecule may be favorable at non-organic solvents and dielectric permittivity values closer to water. Therefore, a mixture of aqueous solvents with possible ionic composition could be used to increase the thermodynamic stability of H[Formula: see text] storage in CO[Formula: see text] conversion processes.

METHODS

Using the Quantum ESPRESSO package, density functional theory (DFT) calculations were performed with periodic boundary conditions, and the electronic wave functions were expanded in plane waves. For the exchange-correlation functional, we use the vdW-DF functional with the inclusion of van der Waals (vdW) forces. Electron-ion interactions are treated by the projector augmented wave (PAW) method with pseudopotentials available in the PSlibrary repository. The wave functions and the electronic densities were expanded employing accurate cut-off energies of 6.80[Formula: see text]10[Formula: see text] and 5.44[Formula: see text]10[Formula: see text] eV, respectively. The electronic density was computed from the wave functions calculated at the [Formula: see text]-point in the first Brillouin-zone. Each structural optimization was minimized according to the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm, with force and energy convergence criteria of 25 meV[Formula: see text]Å[Formula: see text] and 1.36 meV, respectively. The electrostatic solvation effects were performed by the [Formula: see text] package with the Self-Consistent Continuum Solvation (SCCS) approach.

Thermodynamics and electronic structure of adsorbed and intercalated plumbene in graphene/hexagonal SiC heterostructures

Graphene-covered hexagonal SiC substrates have been frequently discussed to be appropriate starting points for epitaxial overlayers of Xenes, such as plumbene, or even their deposition as intercalates between graphene and SiC. Here, we investigate, within density functional theory, the plumbene deposition for various layer orderings and substrate terminations. By means of total energy studies we demonstrate the favorization of the intercalation versus the epitaxy for both C-terminated and Si-terminated 4H-SiC substrates. These results are explained in terms of chemical bonding and by means of layer-resolved projected band structures. Our results are compared with available experimental findings.

Metal-Organic Framework-Based Composites for the Adsorption Removal of Per- and Polyfluoroalkyl Substances from Water

The increasing health risks posed by per- and polyfluoroalkyl substances (PFASs) in the environment highlight the importance of implementing effective removal techniques. Conventional wastewater treatment processes are inadequate for removing persistent organic pollutants. Recent studies have increasingly demonstrated that metal-organic frameworks (MOFs) are capable of removing PFASs from water through adsorption techniques. However, there is still constructive discussion on the potential of MOFs in adsorbing and removing PFASs for large-scale engineering applications. This review systematically investigates the use of MOFs as adsorbents for the removal of PFAS in water treatment. This primarily involved a comprehensive analysis of existing literature to understand the adsorption mechanisms of MOFs and to identify factors that enhance their efficiency in removing PFASs. We also explore the critical aspects of regeneration and stability of MOFs, assessing their reusability and long-term performance, which are essential for large-scale water treatment applications. Finally, our study highlights the challenges of removing PFASs using MOFs. Especially, the efficient removal of short-chain PFASs with hydrophilicity is a major challenge, while medium- to long-chain PFASs are frequently susceptible to being captured from water by MOFs through multiple synergistic effects. The ion-exchange force may be the key to solving this difficulty, but its susceptibility to ion interference in water needs to be addressed in practical applications. We hope that this review can provide valuable insights into the effective removal and adsorption mechanisms of PFASs as well as advance the sustainable utilization of MOFs in the field of water treatment, thereby presenting a novel perspective.

Quantum biochemical analysis of the binding interactions between a potential inhibitory drug and the Ebola viral glycoprotein

Ebola virus disease (EVD) causes outbreaks and epidemics in West Africa that persist until today. The envelope glycoprotein of Ebola virus (GP) consists of two subunits, GP1 and GP2, and plays a key role in anchoring or fusing the virus to the host cell in its active form on the virion surface. Toremifene (TOR) is a ligand that mainly acts as an estrogen receptor antagonist; however, a recent study showed a strong and efficient interaction with GP. In this context, we aimed to evaluate the energetic affinity features involved in the interaction between GP and toremifene by computer simulation techniques using the Molecular Fractionation Method with Conjugate Caps (MFCC) scheme and quantum-mechanical (QM) calculations, as well as missense mutations to assess protein stability. We identified ASP522, GLU100, TYR517, THR519, LEU186, LEU515 as the most attractive residues in the EBOV glycoprotein structure that form the binding pocket. We divided toremifene into three regions and evaluated that region i was more important than region iii and region ii for the formation of the TOR-GP1/GP2 complex, which might control the molecular remodeling process of TOR. The mutations that caused more destabilization were ARG134, LEU515, TYR517 and ARG559, while those that caused stabilization were GLU523 and ASP522. TYR517 is a critical residue for the binding of TOR, and is highly conserved among EBOV species. Our results may help to elucidate the mechanism of drug action on the GP protein of the Ebola virus and subsequently develop new pharmacological approaches against EVD.Communicated by Ramaswamy H. Sarma.

Differences in interaction of graphene/graphene oxide with bacterial and mammalian cell membranes

Graphene, a single layer, hexagonally packed two-dimensional carbon sheet is an attractive candidate for diverse applications including antibacterial potential and drug delivery. One of the knowledge gaps in biomedical application of graphene is the interaction of these materials with the cells. To address this, we investigated the interaction between graphene materials (graphene and graphene oxide) and plasma membranes of cells (bacterial and mammalian cells). The interactions of four of the most abundant phospholipids in bacteria and mammalian plasma membranes with graphene materials were studied using density functional theory (DFT) at the atomic level. The calculations showed that the mammalian phospholipids have stronger bonding to each other compared to bacterial phospholipids. When the graphene/graphene oxide sheet is approaching the phospholipid pairs, the bacterial pairs exhibit less repulsive interactions, thereby a more stable system with the sheets was found. We also assembled bacterial and mammalian phospholipids into liposomes. We further observed that the bacterial liposomes and cells let the graphene flakes penetrate the membrane. The differential scanning calorimetry measurements of liposomes revealed that the bacterial liposomes have the lowest heat capacity; this strengthens the theoretical predictions of weaker interaction between the bacterial phospholipids compared to the mammalian phospholipids. We further demonstrated that graphene oxide could be internalized into the mammalian liposomes without disrupting the membrane integrity. The results suggest that the weak bonding among bacteria phospholipids and less repulsive force when graphene materials approach, result in graphene materials interacting differently with the bacteria compared to mammalian cells.

Development of a machine learning finite-range nonlocal density functional

Kohn-Sham density functional theory has been the most popular method in electronic structure calculations. To fulfill the increasing accuracy requirements, new approximate functionals are needed to address key issues in existing approximations. It is well known that nonlocal components are crucial. Current nonlocal functionals mostly require orbital dependence such as in Hartree-Fock exchange and many-body perturbation correlation energy, which, however, leads to higher computational costs. Deviating from this pathway, we describe functional nonlocality in a new approach. By partitioning the total density to atom-centered local densities, a many-body expansion is proposed. This many-body expansion can be truncated at one-body contributions, if a base functional is used and an energy correction is approximated. The contribution from each atom-centered local density is a single finite-range nonlocal functional that is universal for all atoms. We then use machine learning to develop this universal atom-centered functional. Parameters in this functional are determined by fitting to data that are produced by high-level theories. Extensive tests on several different test sets, which include reaction energies, reaction barrier heights, and non-covalent interaction energies, show that the new functional, with only the density as the basic variable, can produce results comparable to the best-performing double-hybrid functionals, (for example, for the thermochemistry test set selected from the GMTKN55 database, BLYP based machine learning functional gives a weighted total mean absolute deviations of 3.33 kcal/mol, while DSD-BLYP-D3(BJ) gives 3.28 kcal/mol) with a lower computational cost. This opens a new pathway to nonlocal functional development and applications.

Second quantization of many-body dispersion interactions for chemical and biological systems

The many-body dispersion (MBD) framework is a successful approach for modeling the long-range electronic correlation energy and optical response of systems with thousands of atoms. Inspired by field theory, here we develop a second-quantized MBD formalism (SQ-MBD) that recasts a system of atomic quantum Drude oscillators in a Fock-space representation. SQ-MBD provides: (i) tools for projecting observables (interaction energy, transition multipoles, polarizability tensors) on coarse-grained representations of the atomistic system ranging from single atoms to large structural motifs, (ii) a quantum-information framework to analyze correlations and (non)separability among fragments in a given molecular complex, and (iii) a path toward the applicability of the MBD framework to molecular complexes with even larger number of atoms. The SQ-MBD approach offers conceptual insights into quantum fluctuations in molecular systems and enables direct coupling of collective plasmon-like MBD degrees of freedom with arbitrary environments, providing a tractable computational framework to treat dispersion interactions and polarization response in intricate systems.

Direct Observation of Twin van der Waals Molecular Chains

The van der Waals (vdW) assemblies are the most common structures of materials. However, direct mapping of intermolecular electron clouds of a vdW assembly has never been obtained, even though the intramolecular electron clouds were visualized by atomic-resolution techniques. In this report, we unprecedentedly mapped the intermolecular electron cloud of the assemblies of ethanol molecules via ethyl groups with high-resolution atomic force microscopy and scanning tunneling microscopy at 5 K, leading to the first visualization of vdW molecular chains, in which ethanol molecules assemble into twin vdW molecular chains in a reverse parallel configuration on the Ag(111) plane. Furthermore, spontaneous order-disorder transitions in the chain were dynamically observed, suggesting its unusual properties different from those of 2D vdW materials. These findings provide an "eye" to see the atomic world of vdW materials.

libMBD: A general-purpose package for scalable quantum many-body dispersion calculations

Many-body dispersion (MBD) is a powerful framework to treat van der Waals (vdW) dispersion interactions in density-functional theory and related atomistic modeling methods. Several independent implementations of MBD with varying degree of functionality exist across a number of electronic structure codes, which both limits the current users of those codes and complicates dissemination of new variants of MBD. Here, we develop and document libMBD, a library implementation of MBD that is functionally complete, efficient, easy to integrate with any electronic structure code, and already integrated in FHI-aims, DFTB+, VASP, Q-Chem, CASTEP, and Quantum ESPRESSO. libMBD is written in modern Fortran with bindings to C and Python, uses MPI/ScaLAPACK for parallelization, and implements MBD for both finite and periodic systems, with analytical gradients with respect to all input parameters. The computational cost has asymptotic cubic scaling with system size, and evaluation of gradients only changes the prefactor of the scaling law, with libMBD exhibiting strong scaling up to 256 processor cores. Other MBD properties beyond energy and gradients can be calculated with libMBD, such as the charge-density polarization, first-order Coulomb correction, the dielectric function, or the order-by-order expansion of the energy in the dipole interaction. Calculations on supramolecular complexes with MBD-corrected electronic structure methods and a meta-review of previous applications of MBD demonstrate the broad applicability of the libMBD package to treat vdW interactions.

A theoretical study of novel orthorhombic group-IVB nitride halide monolayers for photocatalytic overall water splitting

Hydrogen energy is very important as a new clean energy source to combat the growing environmental problems. In this regard, novel photocatalyst materials for water splitting have a wide range of applications. Using first principles calculations, we theoretically studied three orthorhombic group-IVB nitride halide monolayers, Hf2N2Br2, Janus HfZrN2Br2 and Janus Hf2N2ClBr. The energy, dynamic and thermal stabilities are demonstrated for all three monolayers. Using the HSE hybrid functional, the calculations reveal that they are direct band gap semiconductors with suitable band edge positions, good optical absorptions, and anisotropic carrier mobilities, which makes them promising for water splitting applications. Importantly, the photogenerated carriers provide enough driving force to trigger the hydrogen evolution reaction (HER) and the oxygen evolution reaction (OER) within wide pH ranges, and then overall water splitting can be achieved spontaneously. We conclude that orthorhombic group-IVB nitride halide monolayers have potential applications in photocatalytic nanodevices.

Machine Learning Density Functionals from the Random-Phase Approximation

Kohn-Sham density functional theory (DFT) is the standard method for first-principles calculations in computational chemistry and materials science. More accurate theories such as the random-phase approximation (RPA) are limited in application due to their large computational cost. Here, we use machine learning to map the RPA to a pure Kohn-Sham density functional. The machine learned RPA model (ML-RPA) is a nonlocal extension of the standard gradient approximation. The density descriptors used as ingredients for the enhancement factor are nonlocal counterparts of the local density and its gradient. Rather than fitting only RPA exchange-correlation energies, we also include derivative information in the form of RPA optimized effective potentials. We train a single ML-RPA functional for diamond, its surfaces, and liquid water. The accuracy of ML-RPA for the formation energies of 28 diamond surfaces reaches that of state-of-the-art van der Waals functionals. For liquid water, however, ML-RPA cannot yet improve upon the standard gradient approximation. Overall, our work demonstrates how machine learning can extend the applicability of the RPA to larger system sizes, time scales, and chemical spaces.

Separating Xylene Isomers with a Calcium Metal-Organic Framework

The purification of p-xylene (pX) from its xylene isomers represents a challenging but important industrial process. Herein, we report the efficient separation of pX from its ortho- and meta- isomers by a microporous calcium-based metal-organic framework material (HIAM-203) with a flexible skeleton. At 30 °C, all three isomers are accommodated but the adsorption kinetics of o-xylene (oX) and m-xylene (mX) are substantially slower than that of pX, and at an elevated temperature of 120 °C, oX and mX are fully excluded while pX can be adsorbed. Multicomponent column breakthrough measurements and vapor-phase/liquid-phase adsorption experiments have demonstrated the capability of HIAM-203 for efficient separation of xylene isomers. Ab initio calculations have provided useful information for understanding the adsorption mechanism.

Luminescent Metal-Organic Framework for the Selective Detection of Aldehydes

The detection of toxic, hazardous chemical species is an important task because they pose serious risks to either the environment or human health. Luminescent metal-organic frameworks (LMOFs) as alternative sensors offer rapid and sensitive detection of chemical species. Interactions between chemical species and LMOFs result in changes in the photoluminescence (PL) profile of the LMOFs which can be readily detected using a simple fluorometer. Herein, we report the use of a robust, Zn-based LMOF, [Zn5(μ3-OH)2(adtb)2(H2O)5·5 DMA] (Zn-adtb, LMOF-341), for the selective detection of benzaldehyde. Upon exposure to benzaldehyde, Zn-adtb experiences significant luminescent quenching, as characterized through PL experiments. Photoluminescent titration experiments reveal that LMOF-341 has a detection limit of 64 ppm and a Ksv value of 179 M-1 for benzaldehyde. Furthermore, we study the guest-host interactions that occur between LMOF-341 and benzaldehyde through in situ Fourier transform infrared and computational modeling employing density functional theory. The results show that benzaldehyde interacts more strongly with LMOF-341 compared to formaldehyde and propionaldehyde. Our combined studies also reveal that the mechanism of luminescence quenching originates from an electron-transfer process.

Ultrahigh-printing-speed photoresists for additive manufacturing

Printing technology for precise additive manufacturing at the nanoscale currently relies on two-photon lithography. Although this methodology can overcome the Rayleigh limit to achieve nanoscale structures, it still operates at too slow of a speed for large-scale practical applications. Here we show an extremely sensitive zirconium oxide hybrid-(2,4-bis(trichloromethyl)-6-(4-methoxystyryl)-1,3,5-triazine) (ZrO2-BTMST) photoresist system that can achieve a printing speed of 7.77 m s-1, which is between three and five orders of magnitude faster than conventional polymer-based photoresists. We build a polygon laser scanner-based two-photon lithography machine with a linear stepping speed approaching 10 m s-1. Using the ZrO2-BTMST photoresist, we fabricate a square raster with an area of 1 cm2 in ~33 min. Furthermore, the extremely small chemical components of the ZrO2-BTMST photoresist enable high-precision patterning, leading to a line width as small as 38 nm. Calculations assisted by characterizations reveal that the unusual sensitivity arises from an efficient light-induced polarity change of the ZrO2 hybrid. We envisage that the exceptional sensitivity of our organic-inorganic hybrid photoresist may lead to a viable large-scale additive manufacturing nanofabrication technology.

Template-free scalable growth of vertically-aligned MoS2 nanowire array meta-structural films towards robust superlubricity

Two-dimensional (2D) molybdenum disulfide exhibits a variety of intriguing behaviors depending on its orientation layers. Therefore, developing a template-free atomic layer orientation controllable growth approach is of great importance. Here, we demonstrate scalable, template-free, well-ordered vertically-oriented MoS2 nanowire arrays (VO-MoS2 NWAs) embedded in an Ag-MoS2 matrix, directly grown on various substrates (Si, Al, and stainless steel) via one-step sputtering. In the meta-structured film, vertically-standing few-layered MoS2 NWAs of almost micron length (∼720 nm) throughout the entire film bulk. While near the surface, MoS2 lamellae are oriented in parallel, which are beneficial for caging the bonds dangling from the basal planes. Owing to the unique T-type topological characteristics, chemically inert Ag@MoS2 nano-scrolls (NSCs) and nano-crystalline Ag (nc-Ag) nanoparticles (NPs) are in situ formed under the sliding shear force. Thus, incommensurate contact between (002) basal planes and nc-Ag NPs is observed. As a result, robust superlubricity (friction coefficient μ = 0.0039) under humid ambient conditions is reached. This study offers an unprecedented strategy for controlling the basal plane orientation of 2D transition metal dichalcogenides (TMDCs) via substrate independence, using a one-step solution-free easily scalable process without the need for a template, which promotes the potential applications of 2D TMDCs in solid superlubricity.

Non-Perturbative Determination of Isotope-induced Anomalous Vibrational Physics

In general, vibrational physics has been well described by quantum perturbation theory to provide footprint characteristics for common crystals. However, despite weak phonon anharmonicity, the recently discovered cubic crystals have shown anomalous vibrational dynamics with elusive fundamental origin. Here, we developed a non-perturbative ab initio approach, in together with spectroscopy and high-pressure experiments, to successfully determine the exact dynamic evolutions of the vibrational physics for the first time. We found that the local fluctuation and coupling isotopes significantly dictate the vibrational spectra, through the Brillouin zone folding that has been previously ignored in literature. By decomposing vibrational spectra into individual isotope eigenvectors, we observed both positive and negative contributions to Raman intensity from constitutional atoms (10B, 11B, 75As or 31P). Importantly, our non-perturbative theory predicts that a novel vibrational resonance appears at high hydrostatic pressure due to broken translational symmetry, which was indeed verified by experimental measurement under a pressure up to 31.5 GPa. Our study develops fundamental understandings for the anomalous lattice physics under the failure of quantum perturbation theory and provides a new approach in exploring novel transport phenomena for materials of extreme properties.

Physically-based simulation of elastic-plastic fusion of 3D bioprinted spheroids

Spheroids are microtissues containing cells organized in a spherical shape whose diameter is usually less than a millimetre. Depending on the properties of the environment they are placed in, some nearby spheroids spontaneously fuse and generate a tissue. Given their potential to mimic features typical of body parts and their ability to assemble by fusing in permissive hydrogels, they have been used as building blocks to 3D bioprint human tissue parts. Parameters controlling the shape and size of a bioprinted tissue using fusing spheroid cultures include cell composition, hydrogel properties, and their relative initial position. Hence, simulating, anticipating, and then controlling the spheroid fusion process is essential to control the shape and size of the bioprinted tissue. This study presents the first physically-based framework to simulate the fusion process of bioprinted spheroids. The simulation is based on elastic-plastic solid and fluid continuum mechanics models. Both models use the 'smoothed particle hydrodynamics' method, which is based on discretizing the continuous medium into a finite number of particles and solving the differential equations related to the physical properties (e.g. Navier-Stokes equation) using a smoothing kernel function. To further investigate the effects of such parameters on spheroid shape and geometry, we performed sensitivity and morphological analysis to validate our simulations within-vitrospheroids. Through ourin-silicosimulations by changing the aforementioned parameters, we show that the proposed models appropriately simulate the range of the elastic-plastic behaviours ofin-vitrofusing spheroids to generate tissues of desired shapes and sizes. Altogether, this study presented a physically-based simulation that can provide a framework for monitoring and controlling the geometrical shape of spheroids, directly impacting future research using spheroids for tissue bioprinting.

A Density Functional Theory (DFT) Perspective on Optical Absorption of Modified Graphene Interacting with the Main Amino Acids of Spider Silk

We investigated the possible adsorption of each of the main building blocks of spider silk: alanine, glycine, leucine, and proline. This knowledge could help develop new biocompatible materials and favors the creation of new biosensors. We used ab initio density functional theory methods to study the variations in the optical absorption, reflectivity, and band structure of a modified graphene surface interacting with these four molecules. Four modification cases were considered: graphene with vacancies at 5.55% and fluorine, nitrogen, or oxygen doping, also at 5.55%. We found that, among the cases considered, graphene with vacancies is the best candidate to develop optical biosensors to detect C=O amide and differentiate glycine and leucine from alanine and proline in the visible spectrum region. Finally, from the projected density of states, the main changes occur at deep energies. Thus, all modified graphene's electronic energy band structure undergoes only tiny changes when interacting with amino acids.

Flat Γ Moiré Bands in Twisted Bilayer WSe_{2}

The recent observation of correlated phases in transition metal dichalcogenide moiré systems at integer and fractional filling promises new insight into metal-insulator transitions and the unusual states of matter that can emerge near such transitions. Here, we combine real- and momentum-space mapping techniques to study moiré superlattice effects in 57.4° twisted WSe_{2} (tWSe_{2}). Our data reveal a split-off flat band that derives from the monolayer Γ states. Using advanced data analysis, we directly quantify the moiré potential from our data. We further demonstrate that the global valence band maximum in tWSe_{2} is close in energy to this flat band but derives from the monolayer K states which show weaker superlattice effects. These results constrain theoretical models and open the perspective that Γ-valley flat bands might be involved in the correlated physics of twisted WSe_{2}.

Influence of Dispersion Interactions on the Polymorphic Stability of Crystalline Oxides

The accurate determination of relative phase stabilities using DFT methods is a significant challenge when some of these can vary by only a few kJ/mol. Here, we demonstrate that for a selection of oxides (TiO₂, MnO₂, and ZnO) the inclusion of dispersion interactions, accomplished using the DFT-D3 correction scheme, allows for the correct ordering and an improved calculation of the energy differences between polymorphic phases. The energetic correction provided is of the same order of magnitude as the energy difference between phases. D3-corrected hybrid functionals systematically yield results closest to experiment. We propose that the inclusion of dispersion interactions makes a significant contribution to the relative energetics of polymorphic phases, especially those with different densities, and should therefore be included for calculations of relative energies using DFT methods.

Oxide- and Silicate-Water Interfaces and Their Roles in Technology and the Environment

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by "anomalies" or "non-idealities" such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide- and silicate-water interfaces. This critical review discusses how science progresses from understanding ideal solid-water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. We anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.