Topological wetting states of microdroplets on closed-loop structured surfaces: Breakdown of the Gibbs equation at the microscale.

Microdroplets are a class of soft matter that has been extensively employed for chemical, biochemical, and industrial applications. However, fabricating microdroplets with largely controllable contact-area shape and apparent contact angle, a key prerequisite for their applications, is still a challenge. Here, by engineering a type of surface with homocentric closed-loop microwalls/microchannels, we can achieve facile size, shape, and contact-angle tunability of microdroplets on the textured surfaces by design. More importantly, this class of surface topologies (with universal genus value = 1) allows us to reveal that the conventional Gibbs equation (widely used for assessing the edge effect on the apparent contact angle of macrodroplets) seems no longer applicable for water microdroplets or nanodroplets (evidenced by independent molecular dynamics simulations). Notably, for the flat surface with the intrinsic contact angle ~0°, we find that the critical contact angle on the microtextured counterparts (at edge angle 90°) can be as large as >130°, rather than 90° according to the Gibbs equation. Experiments show that the breakdown of the Gibbs equation occurs for microdroplets of different types of liquids including alcohol and hydrocarbon oils. Overall, the microtextured surface design and topological wetting states not only offer opportunities for diverse applications of microdroplets such as controllable chemical reactions and low-cost circuit fabrications but also provide testbeds for advancing the fundamental surface science of wetting beyond the Gibbs equation.

A droplet-based electricity generator with high instantaneous power density.

Extensive efforts have been made to harvest energy from water in the form of raindrops1-6, river and ocean waves7,8, tides9 and others10-17. However, achieving a high density of electrical power generation is challenging. Traditional hydraulic power generation mainly uses electromagnetic generators that are heavy, bulky, and become inefficient with low water supply. An alternative, the water-droplet/solid-based triboelectric nanogenerator, has so far generated peak power densities of less than one watt per square metre, owing to the limitations imposed by interfacial effects-as seen in characterizations of the charge generation and transfer that occur at solid-liquid1-4 or liquid-liquid5,18 interfaces. Here we develop a device to harvest energy from impinging water droplets by using an architecture that comprises a polytetrafluoroethylene film on an indium tin oxide substrate plus an aluminium electrode. We show that spreading of an impinged water droplet on the device bridges the originally disconnected components into a closed-loop electrical system, transforming the conventional interfacial effect into a bulk effect, and so enhancing the instantaneous power density by several orders of magnitude over equivalent devices that are limited by interfacial effects.

Atomic imaging of the edge structure and growth of a two-dimensional hexagonal ice.

The formation and growth of water-ice layers on surfaces and of low-dimensional ice under confinement are frequent occurrences1-4. This is exemplified by the extensive reporting of two-dimensional (2D) ice on metals5-11, insulating surfaces12-16, graphite and graphene17,18 and under strong confinement14,19-22. Although structured water adlayers and 2D ice have been imaged, capturing the metastable or intermediate edge structures involved in the 2D ice growth, which could reveal the underlying growth mechanisms, is extremely challenging, owing to the fragility and short lifetime of those edge structures. Here we show that noncontact atomic-force microscopy with a CO-terminated tip (used previously to image interfacial water with minimal perturbation)12, enables real-space imaging of the edge structures of 2D bilayer hexagonal ice grown on a Au(111) surface. We find that armchair-type edges coexist with the zigzag edges usually observed in 2D hexagonal crystals, and freeze these samples during growth to identify the intermediate edge structures. Combined with simulations, these experiments enable us to reconstruct the growth processes that, in the case of the zigzag edge, involve the addition of water molecules to the existing edge and a collective bridging mechanism. Armchair edge growth, by contrast, involves local seeding and edge reconstruction and thus contrasts with conventional views regarding the growth of bilayer hexagonal ices and 2D hexagonal matter in general.

Anomalous water transport in narrow-diameter carbon nanotubes.

Carbon nanotubes (CNTs) mimicking the structure of aquaporins support fast water transport, making them strong candidates for building next-generation high-performance membranes for water treatment. The diffusion and transport behavior of water through CNTs or nanoporous graphene can be fundamentally different from those of bulk water through a macroscopic tube. To date, the nanotube-length-dependent physical transport behavior of water is still largely unexplored. Herein, on the basis of molecular dynamics simulations, we show that the flow rate of water through 0.83-nm-diameter (6,6) and 0.96-nm-diameter (7,7) CNTs exhibits anomalous transport behavior, whereby the flow rate increases markedly first and then either slowly decreases or changes slightly as the CNT length l increases. The critical range of l for the flow-rate transition is 0.37 to 0.5 nm. This anomalous water transport behavior is attributed to the l-dependent mechanical stability of the transient hydrogen-bonding chain that connects water molecules inside and outside the CNTs and bypasses the CNT orifice. The results unveil a microscopic mechanism governing water transport through subnanometer tubes, which has important implications for nanofluidic manipulation.

Strong Sliding Ferroelectricity and Interlayer Sliding Controllable Spintronic Effect in Two-Dimensional HgI2 Layers.

Exploration of two-dimensional (2D) sliding ferroelectric (FE) materials with experimentally detectable ferroelectricity and value-added novel functionalities is highly sought for the development of 2D "slidetronics". Herein, based on first-principles calculations, we identify the synthesizable van der Waals (vdW) layered crystals HgX2 (X = Br and I) as a new class of 2D sliding ferroelectrics. Both HgBr2 and HgI2 in 2D multilayered forms adopt the preferential stacking sequence, leading to room temperature stable out-of-plane (vertical) ferroelectricity that can be reversed via the sliding of adjacent monolayers. Owing to strong interlayer coupling and interfacial charge rearrangement, 2D HgI2 layers possess strong sliding ferroelectricity up to 0.16 µC/cm2, readily detectable in experiment. Moreover, robust sliding ferroelectricity and interlayer sliding controllable Rashba spin texture of FE-HgI2 layers enable potential applications as 2D spintronic devices such that the electric control of electron spin detection can be realized at the 2D regime.

Enhanced Direct Exchange Interaction and Hybridization by Single-Atom Linkers for High Curie Temperature and Superior Visible-Light Harvesting in Cr3(CN3)2.

Designing two-dimensional (2D) ferromagnetic (FM) semiconductors with elevated Curie temperature, high carrier mobility, and strong light harvesting is challenging but crucial to the development of spintronics with multifunctionalities. Herein, we show first-principles computation evidence of the 2D metal-organic framework Kagome ferromagnet Cr3(CN3)2. Monolayer Cr3(CN3)2 is predicted to be an FM semiconductor with a record-high Curie temperature of 943 K owing to the use of a single-atom linker (N), which results in strong direct d-p exchange interaction and hybridization between dyz/xz and pz of Cr and N, as well as excellent matching characteristics in energy and symmetry. The single-atom linker structural feature also leads to notable band dispersion and a relatively high carrier mobility of 420 cm2 V-1 s-1. Moreover, under the in-plane strain, 2D Cr3(CN3)2 can be tuned to possess a strong visible-light-harvesting functionality. These novel properties render monolayer Cr3(CN3)2 a distinct 2D ferromagnet with high potential for the development of multifunctional spintronics.

Tunable Surface Wettability via Terahertz Electrowave Controlled Vicinal Subnanoscale Water Layer.

Achieving timely, reversible, and long-range remote tunability over surface wettability is highly demanded across diverse fields, including nanofluidic systems, drug delivery, and heterogeneous catalysis. Herein, using molecular dynamic simulations, we show, for the first time, a theoretical design of electrowetting to achieve remotely controllable surface wettability via using a terahertz wave. The key idea driving the design is the unique terahertz collective vibration identified in the vicinal subnanoscale water layer, which is absent in bulk water, enabling efficient energy transfer from the terahertz wave to the rotational motion of the vicinal subnanoscale water layer. Consequently, a frequency-specific alternating terahertz electric field near the critical strength can significantly affect the local hydrogen-bonding network of the contact water layer on the solid surface, thereby achieving tunable surface wettability.

Evidence of Spontaneous Formation of Two-Dimensional Amorphous Clathrates on Superhydrophilic Surfaces.

Clathrate hydrates reserved in the seabed are often dispersed in the pores of coarse-grained sediments; hence, their formation typically occurs under nanoconfinement. Herein, we show the first molecular dynamics (MD) simulation evidence of the spontaneous formation of two-dimensional (2D) clathrate hydrates on crystal surfaces without conventional nanoconfinement. The kinetic process of 2D clathrate formation is illustrated via simulated single-molecule deposition. 2D amorphous patterns are observed on various superhydrophilic face-centered cubic surfaces. Notably, the formation of 2D amorphous clathrate can occur over a wide range of temperatures, even at room temperature. The strong water-surface interaction, the characteristic properties of guest-gas molecules, and the underlying surface structure dictate the formation of 2D amorphous clathrates. Semiquantitative phase diagrams of 2D clathrates are constructed where representative patterns of 2D clathrates for characteristic gas molecules on prototypical Pd(111) and Pt(111) surfaces are confirmed by independent MD simulations. A tunable pattern of 2D amorphous clathrates is demonstrated by changing the lattice strain of the underlying substrate. Moreover, ab initio MD simulations confirm the stability of 2D amorphous clathrate. The underlining physical mechanism for 2D clathrate formation on superhydrophilic surfaces is elucidated, which offers deeper insight into the crucial role of water-surface interaction.

Superdiffusive Rotation of Interfacial Water on Noble Metal Surface.

Interfacial water on a metal surface acts as an active layer through the reorientation of water, thereby facilitating the energy transfer and chemical reaction across the metal surface in various physicochemical and industrial processes. However, how this active interfacial water collectively behaves on flat noble metal substrates remains largely unknown due to the experimental limitation in capturing librational vibrational motion of interfacial water and prohibitive computational costs at the first-principles level. Herein, by implementing a machine-learning approach to train neural network potentials, we enable performing advanced molecular dynamics simulations with ab initio accuracy at a nanosecond scale to map the distinct rotational motion of water molecules on a metal surface at room temperature. The vibrational density of states of the interfacial water with two-layer profiles reveals that the rotation and vibration of water within the strong adsorption layer on the metal surface behave as if the water molecules in the bulk ice, wherein the O-H stretching frequency is well consistent with the experimental results. Unexpectedly, the water molecules within the adjacent weak adsorption layer exhibit superdiffusive rotation, contrary to the conventional diffusive rotation of bulk water, while the vibrational motion maintains the characteristic of bulk water. The mechanism underlying this abnormal superdiffusive rotation is attributed to the translation-rotation decoupling of water, in which the translation is restrained by the strong hydrogen bonding within the bilayer interfacial water, whereas the rotation is accelerated freely by the asymmetric water environment. This superdiffusive rotation dynamics may elucidate the experimentally observed large fluctuation of the potential of zero charge on Pt and thereby the conventional Helmholtz layer model revised by including the contribution of interfacial water orientation. The surprising superdiffusive rotation of vicinal water next to noble metals will shed new light on the physicochemical processes and the activity of water molecules near metal electrodes or catalysts.

Highly Efficient and Scalable p-i-n Perovskite Solar Cells Enabled by Poly-metallocene Interfaces.

Inverted p-i-n perovskite solar cells (PSCs) are easy to process but need improved interface characteristics with reduced energy loss to prevent efficiency drops when increasing the active photovoltaic area. Here, we report a series of poly ferrocenyl molecules that can modulate the perovskite surface enabling the construction of small- and large-area PSCs. We found that the perovskite-ferrocenyl interaction forms a hybrid complex with enhanced surface coordination strength and activated electronic states, leading to lower interfacial nonradiative recombination and charge transport resistance losses. The resulting PSCs achieve an enhanced efficiency of up to 26.08% for small-area devices and 24.51% for large-area devices (1.0208 cm2). Moreover, the large-area PSCs maintain >92% of the initial efficiency after 2000 h of continuous operation at the maximum power point under 1-sun illumination and 65 °C.

Low-Potential Iodide Oxidation Enables Dual-Atom CoFe─N─C Catalysts for Ultra-Stable and High-Energy-Efficiency Zn-Air Batteries.

The low energy efficiency and limited cycling life of rechargeable Zn-air batteries (ZABs) arising from the sluggish oxygen reduction/evolution reactions (ORR/OERs) severely hinder their commercial deployment. Herein, a zeolitic imidazolate framework (ZIF)-derived strategy associated with subsequent thermal fixing treatment is proposed to fabricate dual-atom CoFe─N─C nanorods (Co1 Fe1 ─N─C NRs) containing atomically dispersed bimetallic Co/Fe sites, which can promote the energy efficiency and cyclability of ZABs simultaneously by introducing the low-potential oxidation redox reactions. Compared to the mono-metallic nanorods, Co1 Fe1 ─N─C NRs exhibit remarkable ORR performance including a positive half-wave potential of 0.933 V versus reversible hydrogen electrode (RHE) in alkaline electrolyte. Surprisingly, after introducing the potassium iodide (KI) additive, the oxidation overpotential of Co1 Fe1 ─N─C NRs to reach 10 mA cm-2 can be significantly reduced by 395 mV compared to the conventional destructive OER. Theoretical calculations show that the markedly decreased overpotential of iodide oxidation can be ascribed to the synergistic effects of neighboring Co─Fe diatomic sites as the unique adsorption sites. Overall, aqueous ZABs assembled with Co1 Fe1 ─N─C NRs and KI as the air-cathode catalyst and electrolyte additive, respectively, can deliver a low charging voltage of 1.76 V and ultralong cycling stability of over 230 h with a high energy efficiency of ≈68%.

Formation of a two-dimensional helical square tube ice in hydrophobic nanoslit using the TIP5P water model.

In extreme and nanoconfinement conditions, the tetrahedral arrangement of water molecules is challenged, resulting in a rich and new phase behavior unseen in bulk phases. The unique phase behavior of water confined in hydrophobic nanoslits has been previously observed, such as the formation of a variety of two-dimensional (2D) ices below the freezing temperature. The primary identified 2D ice phase, termed square tube ice (STI), represents a unique arrangement of water molecules in 2D ice, which can be viewed as an array of 1D ice nanotubes stacked in the direction parallel to the confinement plane. In this study, we report the molecular dynamics (MD) simulations evidence of a novel 2D ice phase, namely, helical square tube ice (H-STI). H-STI is characterized by the stacking of helical ice nanotubes in the direction parallel to the confinement plane. Its structural specificity is evident in the presence of helical square ice nanotubes, a configuration unseen in both STI and single-walled ice nanotubes. A detailed analysis of the hydrogen bonding strength showed that H-STI is a 2D ice phase diverging from the Bernal-Fowler-Pauling ice rules by forming only two strong hydrogen bonds between adjacent molecules along its helical ice chain. This arrangement of strong hydrogen bonds along ice nanotube and weak bonds between the ice nanotube shows a similarity to quasi-one-dimensional van der Waals materials. Ab initio molecular dynamics simulations (over a 30 ps) were employed to further verify H-STI's stability at 1 GPa and temperature up to 200 K.

Formation of porous ice frameworks at room temperature.

Bulk crystalline ices with ultralow densities have been demonstrated to be thermodynamically metastable at negative pressures. However, the direct formation of these bulk porous ices from liquid water at negative pressures is extremely challenging. Inspired by approaches toward porous media based on host-guest chemistry, such as metal-organic frameworks and covalent organic frameworks, we herein demonstrate via molecular dynamics simulations that a class of ultralow-density porous ices with upright channels can be formed spontaneously from liquid water at 300 K with the assistance of carbon nanotube arrays. We refer to these porous ice polymorphs as water oxygen-vertex frameworks (WOFs). Notably, our simulations revealed that the liquid-WOF phase transition is first-order and occurs at room temperature. All the WOFs exhibited the unique structural feature that they can be regarded as assemblies of nanoribbons of hexagonal bilayer ice (2D ice I) at their armchair or zigzag edges. Based on density functional theory calculations, a comprehensive phase diagram of the WOFs was constructed considering both the thermodynamic and thermal stabilities of the porous ices at negative pressures. Like other types of porous media, these WOFs may be applicable to gas storage, purification, and separation. Moreover, these biocompatible porous ice networks may be exploited as medical-related carriers.

A general descriptor for single-atom catalysts decorated with axial ligand.

Decoration of an axial coordination ligand (ACL) on the active metal site is a highly effective and versatile strategy to tune activity of single-atom catalysts (SACs). However, the regulation mechanism of ACLs on SACs is still incompletely known. Herein, we investigate diversified combinations of ACL-SACs, including all 3d-5d transition metals and ten prototype ACLs. We identify that ACLs can weaken the adsorption capability of the metal atom (M) by raising the bonding energy levels of the M-O bond while enhancing dispersity of the d orbital of M. Through examination of various local configurations and intrinsic parameters of ACL-SACs, a general structure descriptor σ is constructed to quantify the structure-activity relationship of ACL-SACs which solely based on a few key intrinsic features. Importantly, we also identified the axial ligand descriptor σACL, as a part of σ, which can serve as a potential descriptor to determine the rate-limiting steps (RLS) of ACL-SACs in experiment. And we predicted several ACL-SACs, namely, CrN4-, FeN4-, CoN4-, RuN4-, RhN4-, OsN4-, IrN4- and PtN4-ACLs, that entail markedly higher activities than the benchmark catalysts of Pt and IrO2, thereby supporting that the general descriptor σ can provide a simple and cost-effective method to assess efficient electrocatalysts.

Identification of Direct Anchoring Sites for Monoatomic Dispersion of Precious Metals (Pt, Pd, Ag) on CeO2 Support.

Monoatomic dispersion of precious metals on the surface of CeO2 nanocrystals is a highly practical approach for dramatically reducing the usage of precious metals while exploiting the unique properties of single-atom catalysts. However, the specific atomic sites for anchoring precious metal atoms on the CeO2 support and underlying chemical mechanism remain partially unknown. Herein, we show that the terminal hydroxyls on the (100) surface are the most stable sites for anchoring Ag atoms on CeO2 , indicating that CeO2 nanocubes are the most efficient substrates to achieve monoatomic dispersion of Ag. Importantly, the newly identified chemical mechanism for single-metal-atom dispersion on CeO2 nanocubes appears to be generic and can thus be extended to other precious metals (Pt and Pd). In fact, our experiments also show that atomically dispersed Pt/Pd species exhibit morphology- and temperature-dependent CO selectivity in the catalytic CO2 hydrogenation reaction.

Resolving Temperature-Dependent Hydrate Nucleation Pathway: The Role of "Transition Layer".

Understanding the nucleation of natural gas hydrate (NGH) at different conditions has important implications to NGH recovery and other industrial applications, such as gas storage and separation. Herein, vast numbers of hydrate nucleation events are traced via molecular dynamics (MD) simulations at different degrees of supercooling (or driving forces). Specifically, to precisely characterize a hydrate nucleus from an aqueous system during the MD simulation, we develop an evolutionary order parameter (OP) to recognize the nucleus size and shape. Subsequently, the free energy landscapes of hydrate during nucleation are explored by using the newly developed OP. The results suggest that at 270 K (or 0.92 Tm supercooling, where Tm is the melting point), the near-rounded nucleus prevails during the nucleation, as described from the classical nucleation theory. In contrast, at relatively strong driving forces of 0.85 and 0.88 Tm, nonclassical nucleation events arise. Specifically, the pathway toward an elongated nucleus becomes as important as the pathway toward a near-rounded nucleus. To explain the distinct nucleation phenomena at different supercoolings, a notion of a "transition layer" (or liquid-blob-like layer) is proposed. Here, the transition layer is to describe the interfacial region between the nucleus and aqueous solution, and this layer entails two functionalities: (1) it tends to retain CH4 depending on the degrees of supercooling and (2) it facilitates collision among CH4, which thus promote the incorporation of CH4 into nucleus. Our simulation indicates that compared to the near-rounded nucleus, the transition layer surrounding the elongated nucleus is more evident with the higher collision rate among CH4 molecules. As such, the transition layer tends to promote the elongated nucleus pathway, while offsetting the cost of larger surface free energy associated with the elongated nucleus. At 0.92 Tm, however, the transition layer gradually disappears, and classical nucleation events dominate. Overall, the notion of "transition layer" offers deeper insight into the NGH nucleation at different degrees of supercooling and could be extended to describe other types of hydrate nucleation.

Unraveling nucleation pathway in methane clathrate formation.

Methane clathrates are widespread on the ocean floor of the Earth. A better understanding of methane clathrate formation has important implications for natural-gas exploitation, storage, and transportation. A key step toward understanding clathrate formation is hydrate nucleation, which has been suggested to involve multiple evolution pathways. Herein, a unique nucleation/growth pathway for methane clathrate formation has been identified by analyzing the trajectories of large-scale molecular dynamics (MD) simulations. In particular, ternary water-ring aggregations (TWRAs) have been identified as fundamental structures for characterizing the nucleation pathway. Based on this nucleation pathway, the critical nucleus size and nucleation timescale can be quantitatively determined. Specifically, a methane hydration layer compression/shedding process is observed to be the critical step in (and driving) the nucleation/growth pathway, which is manifested through overlapping/compression of the surrounding hydration layers of the methane molecules, followed by detachment (shedding) of the hydration layer. As such, an effective way to control methane hydrate nucleation is to alter the hydration layer compression/shedding process during the course of nucleation.

Influence of Organic-Cation Defects on Optoelectronic Properties of ASnI3 Perovskites A=HC(NH2 )2 , CH3 NH3.

Tin organic-inorganic halide perovskites (tin OIHPs) possess a desirable band gap and their power conversion efficiency (PCE) has reached 14 %. A commonly held view is that the organic cations in tin OIHPs would have little impact on the optoelectronic properties. Herein, we show that the defective organic cations with randomly dynamic characteristics can have marked effect on optoelectronic properties of the tin OIHPs. Hydrogen vacancies originated from the proton dissociation from FA [HC(NH2 )2 ] in FASnI3 can induce deep transition levels in the band gap but yield relatively small nonradiative recombination coefficients of 10-15  cm3 s-1 , whereas those from MA (CH3 NH3 ) in MASnI3 can yield much larger nonradiative recombination coefficients of 10-11  cm3 s-1 . Additional insight into the "defect tolerance" is gained by disentangling the correlations between dynamic rotation of organic cations and charge-carrier dynamics.

Photocatalytic Oxidation of NO2 on TiO2 : Evidence of a New Source of N2 O5.

N2 O5 is an important intermediate in the atmospheric nitrogen cycle. Using a flow tube reactor, N2 O5 was found to be released from the TiO2 surface during the photocatalytic oxidation of NO2 , revealing a previously unreported source of N2 O5 . The rate of N2 O5 release from TiO2 was dependent on the initial NO2 concentration, relative humidity, O2 /N2 ratio, and irradiation intensity. Experimental evidence and quantum chemical calculations showed that NO2 can react with the surface hydroxyl groups and the generated electron holes on the TiO2 , followed by combining with another NO2 molecule to form N2 O5 . The latter was physisorbed on TiO2 and had a low adsorption energy of -0.13 eV. Box model simulations indicated that the new source of N2 O5 released from TiO2 can increase the daytime N2 O5 concentration by up to 20 % in urban areas if abundant TiO2 -containing materials and high NOx concentrations were present. This joint experimental/theoretical study not only demonstrates a new chemical mechanism for N2 O5 formation but also has important implications for air quality in urban areas.

Solvation and Hydrolysis Reaction of Isocyanic Acid at the Air-Water Interface: A Computational Study.

Isocyanic acid (HNCO) is known to be inert to strong oxidants and photolysis in the atmosphere but often appears in different forms of smoke; therefore, it is linked to various smoke-related illnesses due to tobacco usage or wildfire events. To date, the major loss pathway of HNCO is believed to be through its uptake on aerosol droplets. However, the molecular mechanisms underlying such an uptake process are still incompletely understood. Herein, we use the Born-Oppenheimer molecular dynamics (BOMD) simulations to study solvation and hydrolysis reactions of HNCO on water droplets at ambient temperature. The BOMD simulations indicate that the scavenging of HNCO by water droplets is largely attributed to the preferential adsorption of HNCO at the air-water interface, rather than inside bulk water. Specifically, the H atom of HNCO interacts with the O atom of interfacial water, leading to the formation of a hydrogen bond (H-bond) of (HNCO)H···O(H2O), which prevents HNCO from evaporating. Moreover, the interfacial water can act as H-bond acceptors/donors to promote the proton transfer during the HNCO hydrolysis reaction. Compared to the gas phase, the activation barrier is lowered from 45 to 14 kcal·mol-1 on the water surface, which facilitates the formation of the key intermediate of NH2COOH. This intermediate eventually decomposes into NH3 and CO2, consistent with the previous study [ Atmos. Chem. Phys. 2016, 16, 703-714]. The new molecular insight into HNCO solvation and reaction on the water surface improves our understanding of the uptake of HNCO on aerosols.