Advances in Two-Dimensional Ion-Selective Membranes: Bridging Nanoscale Insights to Industrial-Scale Salinity Gradient Energy Harvesting.

Salinity gradient energy, often referred to as the Gibbs free energy difference between saltwater and freshwater, is recognized as "blue energy" due to its inherent cleanliness, renewability, and continuous availability. Reverse electrodialysis (RED), relying on ion-selective membranes, stands as one of the most prevalent and promising methods for harnessing salinity gradient energy to generate electricity. Nevertheless, conventional RED membranes face challenges such as insufficient ion selectivity and transport rates and the difficulty of achieving the minimum commercial energy density threshold of 5 W/m2. In contrast, two-dimensional nanostructured materials, featuring nanoscale channels and abundant functional groups, offer a breakthrough by facilitating rapid ion transport and heightened selectivity. This comprehensive review delves into the mechanisms of osmotic power generation within a single nanopore and nanochannel, exploring optimal nanopore dimensions and nanochannel lengths. We subsequently examine the current landscape of power generation using two-dimensional nanostructured materials in laboratory-scale settings across various test areas. Furthermore, we address the notable decline in power density observed as test areas expand and propose essential criteria for the industrialization of two-dimensional ion-selective membranes. The review concludes with a forward-looking perspective, outlining future research directions, including scalable membrane fabrication, enhanced environmental adaptability, and integration into multiple industries. This review aims to bridge the gap between previous laboratory-scale investigations of two-dimensional ion-selective membranes in salinity gradient energy conversion and their potential large-scale industrial applications.

Interlink between Abnormal Water Imbibition in Hydrophilic and Rapid Flow in Hydrophobic Nanochannels.

Nanoscale extension and refinement of the Lucas-Washburn model is presented with a detailed analysis of recent experimental data and extensive molecular dynamics simulations to investigate rapid water flow and water imbibition within nanocapillaries. Through a comparative analysis of capillary rise in hydrophilic nanochannels, an unexpected reversal of the anticipated trend, with an abnormal peak, of imbibition length below the size of 3 nm was discovered in hydrophilic nanochannels, surprisingly sharing the same physical origin as the well-known peak observed in flow rate within hydrophobic nanochannels. The extended imbibition model is applicable across diverse spatiotemporal scales and validated against simulation results and existing experimental data for both hydrophilic and hydrophobic nanochannels.

Transcriptomic and physiological analysis of atractylodes chinensis in response to drought stress reveals the putative genes related to sesquiterpenoid biosynthesis.

BACKGROUND: Atractylodes chinensis (DC) Koidz., a dicotyledonous and hypogeal germination species, is an important medicinal plant because its rhizome is enriched in sesquiterpenes. The development and production of A. chinensis are negatively affected by drought stress, especially at the seedling stage. Understanding the molecular mechanism of A. chinensis drought stress response plays an important role in ensuring medicinal plant production and quality. In this study, A. chinensis seedlings were subjected to drought stress treatment for 0 (control), 3 (D3), and 9 days (D9). For the control, the sample was watered every two days and collected on the second morning after watering. The integration of physiological and transcriptomic analyses was carried out to investigate the effects of drought stress on A. chinensis seedlings and to reveal the molecular mechanism of its drought stress response. RESULTS: The malondialdehyde, proline, soluble sugar, and crude protein contents and antioxidative enzyme (superoxide dismutase, peroxidase, and catalase) activity were significantly increased under drought stress compared with the control. Transcriptomic analysis indicated a total of 215,665 unigenes with an average length of 759.09 bp and an N50 of 1140 bp. A total of 29,449 differentially expressed genes (DEGs) were detected between the control and D3, and 14,538 DEGs were detected between the control and D9. Under drought stress, terpenoid backbone biosynthesis had the highest number of unigenes in the metabolism of terpenoids and polyketides. To identify candidate genes involved in the sesquiterpenoid and triterpenoid biosynthetic pathways, we observed 22 unigene-encoding enzymes in the terpenoid backbone biosynthetic pathway and 15 unigene-encoding enzymes in the sesquiterpenoid and triterpenoid biosynthetic pathways under drought stress. CONCLUSION: Our study provides transcriptome profiles and candidate genes involved in sesquiterpenoid and triterpenoid biosynthesis in A. chinensis in response to drought stress. Our results improve our understanding of how drought stress might affect sesquiterpenoid and triterpenoid biosynthetic pathways in A. chinensis.

Unexpected Behavior in Thermal Conductivity of Confined Monolayer Water.

Monolayer water can be formed under extreme confinement and will present distinctive thermodynamic properties compared with bulk water. In this work, we perform molecular dynamics simulations to study the thermal conductivity of monolayer water confined in graphene channels, finding an unexpected way of thermal conductivity of monolayer water dependent on its number density, which has a close correlation with the structure of water. The monolayer water is in an amorphous state, and its thermal conductivity increases linearly with the area density when the water density is low at first. Then, the thermal conductivity increases as the number density of water rises, which is attributed to the formation of a crystal structure and the reduction of crystal defects as the number of water molecules increases. After reaching the zenith, the thermal conductivity decreases rapidly owing to the formation of a wrinkle structure of monolayer water with excessive water molecules, which weakens the phonon dispersion. Moreover, we further investigate the remarkable effects of the channel height on both the structure and thermal conductivity of monolayer water. In summary, this study demonstrates the close connection between the thermal conductivity of monolayer water and its structure, contributing to not only expanding the understanding of the thermodynamic property of nanoconfined water but also benefiting the engineering applications for nanofluidics.

CsTRM5 regulates fruit shape via mediating cell division direction and cell expansion in cucumber.

Fruit shape and size are important appearance and yield traits in cucumber, but the underlying genes and their regulatory mechanisms remain poorly understood. Here we identified a mutant with spherical fruits from an Ethyl Methane Sulfonate (EMS)-mutagenized library, named the qiu mutant. Compared with the cylindrical fruit shape in 32X (wild type), the fruit shape in qiu was round due to reduced fruit length and increased fruit diameter. MutMap analysis narrowed the candidate gene in the 6.47 MB range on Chr2, harboring the FS2.1 locus reported previously. A single-nucleotide polymorphism (SNP) (11359603) causing a truncated protein of CsaV3_2G013800, the homolog of tomato fruit shape gene SlTRM5, may underlie the fruit shape variation in the qiu mutant. Knockout of CsTRM5 by the CRISPR-Cas9 system confirmed that CsaV3_2G013800/CsTRM5 was the causal gene responsible for qiu. Sectioning analysis showed that the spherical fruit in qiu resulted mainly from increased and reduced cell division along the transverse and longitudinal directions, respectively. Meanwhile, the repressed cell expansion contributed to the decreased fruit length in qiu. Transcriptome profiling showed that the expression levels of cell-wall-related genes and abscisic acid (ABA) pathway genes were significantly upregulated in qiu. Hormone measurements indicated that ABA content was greatly increased in the qiu mutant. Exogenous ABA application reduced fruit elongation by inhibiting cell expansion in cucumber. Taken together, these data suggest that CsTRM5 regulates fruit shape by affecting cell division direction and cell expansion, and that ABA participates in the CsTRM5-mediated cell expansion during fruit elongation in cucumber.

A hydrogen bond-modulated soft nanoscale water channel for ion transport through liquid-liquid interfaces.

Ion transport through interfaces is of ubiquitous importance in many fields such as electrochemistry, emulsion stabilization, phase transfer catalysis, liquid-liquid extraction and enhanced oil recovery. However, the knowledge of interfacial structures that significantly affect ion transport through liquid-liquid interfaces is still lacking due to the difficulty of observing nanoscale interfaces. We studied here the evolution of interfacial structures during ion transport through the decane-water interface under different ionic concentrations and external forces using molecular dynamics simulations. The roles of hydrogen bonds in ion transport through interfaces are revealed. We identified a soft nanoscale channel during ion transport through liquid-liquid interfaces and the decane phase under specific external force. The stability of the water channel and the ion transport velocity both increase with ionic concentration due to the layered ordering structures of the water near the channel surface. We observed that the stability and connectivity of the water channel in the decane phase are remarkably improved both by the high increase of the number of hydrogen bonds in the water channel with increasing ionic concentration, and by the conformational change in water molecules near the water channel surface. Our discovery of a soft nanoscale water channel by molecular simulations implies that there is a potential stable passage for ion transport through liquid-liquid interfaces.

Comparative transcriptome analysis reveals sesquiterpenoid biosynthesis among 1-, 2- and 3-year old Atractylodes chinensis.

BACKGROUND: Atractylodes chinensis (DC.) Koidz is a well-known medicinal plant containing the major bioactive compound, atractylodin, a sesquiterpenoid. High-performance liquid chromatography (HPLC) analysis demonstrated that atractylodin was most abundant in 3-year old A. chinensis rhizome, compared with those from 1- and 2-year old rhizomes, however, the molecular mechanisms underlying accumulation of atractylodin in rhizomes are poorly understood. RESULTS: In this study, we characterized the transcriptomes from rhizomes of 1-, 2- and 3-year old (Y1, Y2 and Y3, respectively) A. chinensis, to identify differentially expressed genes (DEGs). We identified 240, 169 and 131 unigenes encoding the enzyme genes in the mevalonate (MVA), methylerythritol phosphate (MEP), sesquiterpenoid and triterpenoid biosynthetic pathways, respectively. To confirm the reliability of the RNA sequencing analysis, eleven key gene encoding factors involved in the sesquiterpenoid and triterpenoid biosynthetic pathway, as well as in pigment, amino acid, hormone and transcription factor functions, were selected for quantitative real time PCR (qRT-PCR) analysis. The results demonstrated similar expression patterns to those determined by RNA sequencing, with a Pearson's correlation coefficient of 0.9 between qRT-PCR and RNA-seq data. Differential gene expression analysis of rhizomes from different ages revealed 52 genes related to sesquiterpenoid and triterpenoid biosynthesis. Among these, seven DEGs were identified in Y1 vs Y2, Y1 vs Y3 and Y2 vs Y3, of which five encoded four key enzymes, squalene/phytoene synthase (SS), squalene-hopene cyclase (SHC), squalene epoxidase (SE) and dammarenediol II synthase (DS). These four enzymes directly related to squalene biosynthesis and subsequent catalytic action. To validate the result of these seven DEGs, qRT-PCR was performed and indicated most of them displayed lower relative expression in 3-year old rhizome, similar to transcriptomic analysis. CONCLUSION: The enzymes SS, SHC, SE and DS down-regulated expression in 3-year old rhizome. This data corresponded to the higher content of sesquiterpenoid in 3-year old rhizome, and confirmed by qRT-PCR. The results of comparative transcriptome analysis and identified key enzyme genes laid a solid foundation for investigation of production sesquiterpenoid in A. chinensis.

Extending the Classical Continuum Theory to Describe Water Flow through Two-Dimensional Nanopores.

Water flow through two-dimensional nanopores has attracted significant attention owing to the promising water purification technology based on atomically thick membranes. However, the theoretical description of water flow in nanopores based on the classical continuum theory is very challenging owing to the pronounced entrance/exit effects. Here, we extend the classical Hagen-Poiseuille equation for describing the relationship between flow rate and pressure loss in laminar tube flow to two-dimensional nanopores. A totally theoretical model is established by appropriately considering the velocity slip on pore surfaces both in the friction pressure loss and entrance/exit pressure loss. Based on molecular dynamics simulations of water flow through graphene nanopores, it is shown that the model can not only well predict the overall flow rate but also give a good estimation of the velocity profiles. As the pore radius and length increase, the model can reduce to the equations applicable to the fluid flow in infinitely/finitely long nanotubes, thin orifices, and macroscale tubes, showing an accurate prediction of the existing experimental and simulation data of the water flow through nanotubes and nanopores in the literature. Namely, the presented model is a unified model that can uniformly describe the fluid flow from nanoscales to macroscales by modifying the classical continuum theory.

Theoretical description of molecular permeation via surface diffusion through graphene nanopores.

We establish a theoretical model to describe the surface molecular permeation through two-dimensional graphene nanopores based on the surface diffusion equation and Fick's law. The model is established by considering molecular adsorption and desorption from the surface adsorption layer and the molecular diffusion and concentration gradient on the graphene surface. By comparing with the surface flux obtained from molecular dynamics simulations, it is shown that the model can predict well the overall permeation flux especially for strongly adsorbed molecules (i.e. CO2 and H2S) on graphene surfaces. Although good agreement between the theoretical and simulated density distribution is hard to achieve owing to the large uncertainty in the calculation of surface diffusion coefficients based on the Einstein equation, the model itself is very competent to describe the surface molecular permeation both from the aspects of the overall permeation flux and detailed density distribution. This model is believed to supplement the theoretical description of molecular permeation through graphene nanopores and provide a good reference for the description of mass transport through two-dimensional porous materials.

Wall friction should be decoupled from fluid viscosity for the prediction of nanoscale flow.

The accurate determination of fluid viscosity based on the microscopic information of molecules is very crucial for the prediction of nanoscale flow. Despite the challenge of this problem, researchers have done a lot of meaningful work and developed several distinctive methods. However, one of the common approaches to calculate the fluid viscosity is using the Green-Kubo formula by considering all the fluid molecules in nanospace, inevitably causing the involvement of the frictional interaction between fluid and the wall into the fluid viscosity. This practice is certainly not appropriate because viscosity is essentially related only to the interactions among fluid molecules. Here, we clarify that the wall friction should be decoupled from fluid viscosity by distinguishing the frictional region and the viscous region for the accurate prediction of nanoscale flow. By comparing the fluid viscosities calculated from the Green-Kubo formula in the whole region and viscous region and the viscosity obtained from the velocity profile through the Hagen-Poiseuille equation, it is found that only the calculated viscosity in the viscous region agrees well with the viscosity from the velocity profile. To demonstrate the applicability of this clarification, the Lennard-Jones fluid and water confined between Lennard-Jones, graphene, and silica walls, even with different fluid-wall interactions, are extensively tested. This work clearly defines the viscosity of fluids at nanoscales from the inherent nature of physics, aiming at the accurate prediction of nanoscale flow from the classical continuum hydrodynamic theory.

Hierarchical thermal transport in nanoconfined water.

The structure of nanoconfined fluids is particularly non-uniform owing to the wall interaction, resulting in the distinctive characteristic of thermal transport compared to bulk fluids. We present the molecular simulations on the thermal transport of water confined in nanochannels with a major investigation of its spatial distribution under the effects of wall interaction. The results show that the thermal conductivity of nanoconfined water is inhomogeneous and its layered distribution is very similar to the density profile. The layered thermal conductivity is the coupling result of inhomogeneous density and energy distributions that are generally diametrical, and their contributions to the thermal conductivity compensate with each other. However, the accumulative effect of water molecules is really dominating, resulting in a high thermal conductivity in the high-density layers with the low-energy molecules, and vice versa. Moreover, it is found that the adsorptive and repulsive interactions from solid walls have different roles in the hierarchical thermal transport in nanoconfined water. The adsorptive interaction is only responsible for the layered distribution of thermal conductivity, while the repulsive interaction is responsible for the overall thermal conductivity; accordingly, the thermal conductivity is independent of the strength of water-solid interactions. The identified hierarchical thermal transport in nanoconfined water and its underlying mechanisms have a great significance for the understanding of nanoscale thermal transport and even the mass and energy transport of nanoconfined fluids.

Asymmetric Two-Layer Porous Membrane for Gas Separation.

We present that the porous two-layer membranes of graphene and hexagonal boron nitride (h-BN) are promising for gas mixture separation. For the two-layer membranes, the mechanisms of the gas separation are (i) the different adsorption properties of gases on two membranes inducing a permeation flux difference from one side to the other and (ii) the asymmetric potential energy curves (potential energy of a gas molecule vs distance between the pore center and a gas molecule) of a two-layer membrane leading to a potential energy difference, which can affect gas permeation through the pore. As a concrete example, we explore the gas separation of CO2 and CH4 by the two-layer membrane using molecular dynamics simulations. Finally, on the basis of the distinctive permeation rates in the two directions, a gas separation system with two back-to-back arrayed graphene/h-BN membranes with big pores is designed to realize gas separation.

Nanoconfined Fluids: What Can We Expect from Them?

Nanoconfined fluids (NCFs), which are confined in nanospaces, exhibit distinctive nanoscale effects, including surface effects, small-size effects, quantum effects, and others. The continuous medium hypothesis in fluid mechanics is not valid in this context because of the comparable characteristic length of spaces and molecular mean free path, and accordingly, the classical continuum theories developed for the bulk fluids usually cannot describe the mass and energy transport of NCFs. In this Perspective, we summarize the nanoscale effects on the thermodynamics, mass transport, flow dynamics, heat transfer, phase change, and energy transport of NCFs and highlight the related representative works. The applications of NCFs in the fields of membrane separation, oil and gas production, energy harvesting and storage, and biological engineering are especially indicated. Currently, the theoretical description framework of NCFs is still missing, and it is expected that this framework can be established by adopting the classical continuum theories with the consideration of nanoscale effects.

Evidence for water ridges at oil-water interfaces: implications for ion transport.

Understanding ion transport across interfaces is of fundamental importance in many processes such as liquid-liquid extraction, phase transfer catalysis, enhanced oil recovery and emulsion stabilisation. However, the factors that control ion transport across interfaces are poorly known due to a lack of knowledge of structural changes at interfaces. We studied here the effects of ionic concentration and external force on the transport of ions across the decane-water interface using classical molecular dynamics simulations. The results show that the evolution of interfacial structures during ion transfer across the interface is controlled by hydrogen bonding and ionic interactions at the interface. We also identified a new mode of ion transfer across the interface at low ionic concentrations, involving a 'water ridge', rather that the classical 'water finger'. In the water ridge mode, hydrogen bonds are not broken due to low ion levels, and the water ridge induces gradual interface deformation. Whereas, at high ionic concentrations, hydrogen bonds are broken by the strong ion electrostatic repulsion, thus inducing the formation of a water finger. We also found that the variation of the Gibbs free energy during ion transfer is directly relevant to the ionic concentration. The water ridge at low ionic concentrations, which displaces more water molecules towards the decane phase, induces less free energy variation than the water finger at high ionic concentrations.

Effects of Molecular Chain Length on the Contact Line Movement in Water/n-Alkane/Solid Systems.

The movement of the contact line in liquid-liquid-solid systems is a major phenomenon in natural and industrial processes. In particular, n-alkanes are widely occurring in the oil, soil pollution, and chemical industries, yet there is little knowledge on the effects of molecular chain length on the contact line movement. Here, we studied the effects of molecular chain length on the contact line movement in water/n-alkane/solid systems with different surface wettabilities. We used n-heptane (C7), n-decane (C10), and n-hexadecane (C16) as alkanes and α-quartz as the solid surface. We calculated the time-variation contact line moving velocity and also analyzed the jump frequency and the mean distance of the molecular displacement occurring within the contact line zone by molecular-kinetic theory. Molecular dynamics simulation results show that the contact line velocity decreases with increasing the chain length, originally caused by the decreasing the jump frequency and mean distance. These variations with the molecular chain length are related to the more torsions and deformations of the molecules with a longer chain length. In addition, the moving mechanism of the contact line on the same solid surface does not change at different molecular chain lengths, implying that the moving mechanism mainly depends on the three-phase wettability.

Selective Molecular Sieving through a Large Graphene Nanopore with Surface Charges.

The precise control of the pore sizes at an atomic level has proved to be the biggest challenge of all for nanoporous graphene membranes for gas separation. Here, we propose a simple method to realize the selective molecular sieving through originally nonselective graphene nanopores by adding charges on the graphene surfaces. Molecular dynamic simulations show that the CO2/N2 selectivity of the graphene nanopore with a diameter of 0.52 nm increases up to 22.78 for a surface charge density of only -5.934 e/nm2. The selectivity improvement is related to the distinctive adsorption intensities of CO2 and N2 molecules on the charge-loaded graphene surfaces. This work points toward a promising road to tune the selectivity of graphene nanopores and therefore promotes the realization of porous graphene membranes and other two-dimensional porous membranes by accepting the pores with a wide size distribution and reducing the requirements in the control of pore sizes.

A Functional Allele of CsFUL1 Regulates Fruit Length through Repressing CsSUP and Inhibiting Auxin Transport in Cucumber.

Fruit length is a prominent agricultural trait during cucumber (Cucumis sativus) domestication and diversifying selection; however, the regulatory mechanisms of fruit elongation remain elusive. We identified two alleles of the FRUITFULL (FUL)-like MADS-box gene CsFUL1 with 3393C/A Single Nucleotide Polymorphism variation among 150 cucumber lines. Whereas CsFUL1A was specifically enriched in the long-fruited East Asian type cucumbers (China and Japan), the CsFUL1C allele was randomly distributed in cucumber populations, including wild and semiwild cucumbers. CsFUL1A knockdown led to further fruit elongation in cucumber, whereas elevated expression of CsFUL1A resulted in significantly shorter fruits. No effect on fruit elongation was detected when CsFUL1C expression was modulated, suggesting that CsFUL1A is a gain-of-function allele in long-fruited cucumber that acts as a repressor during diversifying selection of East Asian cucumbers. Furthermore, CsFUL1A binds to the CArG-box in the promoter region of SUPERMAN, a regulator of cell division and expansion, to repress its expression. Additionally, CsFUL1A inhibits the expression of auxin transporters PIN-FORMED1 (PIN1) and PIN7, resulting in decreases in auxin accumulation in fruits. Together, our work identifies an agriculturally important allele and suggests a strategy for manipulating fruit length in cucumber breeding that involves modulation of CsFUL1A expression.

Moving mechanisms of the three-phase contact line in a water-decane-silica system.

The movement of the three-phase contact line with chain molecules in the liquid phase displays more complex mechanisms compared to those in the usual liquid-liquid-solid systems and even to the gas-liquid-solid systems controlled by the traditional single-molecule adsorption-desorption mechanisms. By introducing decane molecules with chain structures, we demonstrate from molecular dynamics insights that the moving mechanism of the contact line in a water-decane-silica system is totally different from traditional mechanisms. Three different wettability-related moving mechanisms including "Roll up", "Piston" and "Shear" are revealed corresponding to the hydrophilic, intermediate and hydrophobic three-phase wettability, respectively. In the "Roll up" mechanism, the decane molecules are rolled up by the competitively adsorbed water molecules and then move forward under the driving force; when the "Piston" mechanism happens, the decane molecules are pushed by the piston-like water phase owing to the comparable adsorption interactions of the two liquids on the solid surface; in the "Shear" mechanism, the contact line is hard to drive due to the stronger decane-silica interactions but the decane molecules far away from the solid surface will move forward. Besides, the time-averaged velocity of the moving contact line is greatly related to the moving mechanisms. For the "Roll up" mechanism, the contact line velocity increases first and then reaches a steady value; for the "Piston" mechanism, the contact line velocity has a maximum value at the start-up stage and then decreases to a stable value; for the "Shear" mechanism, the contact line velocity fluctuates around zero due to the thermal fluctuation of the molecules. Additionally, the mean distance from Molecular Kinetics Theory increases with decreasing hydrophilicity and the displacement frequency in "Roll up" mechanism is 2 orders of magnitude higher than that in the "Piston" mechanism, further demonstrating the different moving mechanisms from a quantitative point of view.

Nanoparticle-induced ion-sensitive reduction in decane-water interfacial tension.

The synergistic effect of ions and nanoparticles on the interfacial tension is of great significance for extensive applications in interface-related industrial processes. However, its mechanisms are still unclear owing to a lack of understanding on the interaction between nanoparticles/ions at the interface. Here, we employ the molecular dynamics method to explore the synergistic effect of ions and nanoparticles on reducing the decane-water interfacial tension and reveal the dominant role of the three-phase contact angle and the interaction between nanoparticles. The results show that the reduction of interfacial tension is sensitive to cation species and temperature. The stronger hydration of cations induces an increased three-phase contact angle, weakening the interaction between nanoparticles and water molecules at the interface. Hence, the virial term of interfacial tension decreases. Meanwhile, the potential of mean force between nanoparticles at the interface indicates that the order of interaction strength between nanoparticles for different cations is Ca2+ > Mg2+ > Na+. The strong interaction between nanoparticles restricts the motion of nanoparticles and water molecules at the interface, inducing a reduced kinetic energy term of interfacial tension. Therefore, the interfacial tension decreases after adding the nanoparticles. Besides, as temperature rises, the difference in the adsorption ability of nanoparticles on water molecules causes a falling interfacial tension with a characteristic stage.

Oil Contact Angles in a Water-Decane-Silicon Dioxide System: Effects of Surface Charge.

Oil wettability in the water-oil-rock systems is very sensitive to the evolution of surface charges on the rock surfaces induced by the adsorption of ions and other chemical agents in water flooding. Through a set of large-scale molecular dynamics simulations, we reveal the effects of surface charge on the oil contact angles in an ideal water-decane-silicon dioxide system. The results show that the contact angles of oil nano-droplets have a great dependence on the surface charges. As the surface charge density exceeds a critical value of 0.992 e/nm2, the contact angle reaches up to 78.8° and the water-wet state is very apparent. The variation of contact angles can be confirmed from the number density distributions of oil molecules. With increasing the surface charge density, the adsorption of oil molecules weakens and the contact areas between nano-droplets and silicon dioxide surface are reduced. In addition, the number density distributions, RDF distributions, and molecular orientations indicate that the oil molecules are adsorbed on the silicon dioxide surface layer-by-layer with an orientation parallel to the surface. However, the layered structure of oil molecules near the silicon dioxide surface becomes more and more obscure at higher surface charge densities.