Wettability of pristine and alkyl-functionalized graphane

Effectiveness of inhibitors to prevent asphaltene aggregation: Insights from atomistic and molecular simulations

The technological landscape for industrial processes handling asphaltene is evolving at a rapid pace due to the increase in the extraction of heavy crude oil. The main underlying challenges in this regard are the flow assurance, the recovery of the spent solvent, and the sophisticated extractor setup required to develop the process to an industrial scale. The number of studies focused on the handling of the asphaltene at the atomic and molecular scales is growing enormously in order to identify new sustainable solvents for the effective extraction of asphaltene from heavy crude oil or oil-bearing sands. This Perspective focuses on the importance of density functional theory and molecular dynamics simulations to explore the broader range of asphaltene inhibitors, e.g., nanoparticles, ionic liquids, and deep eutectic solvents, to prevent asphaltene precipitation. We provide a concise overview of the major accomplishments, analyze the aspects that require attention, and highlight the path-breaking studies having a significant impact on the process of chemical enhanced oil recovery from heavy crude oil reservoirs primarily based on atomistic and molecular simulations.

Reversible Surface Energy Storage in Molecular-Scale Porous Materials

Forcible wetting of hydrophobic pores represents a viable method for energy storage in the form of interfacial energy. The energy used to fill the pores can be recovered as pressure-volume work upon decompression. For efficient recovery, the expulsion pressure should not be significantly lower than the pressure required for infiltration. Hysteresis of the wetting/drying cycle associated with the kinetic barrier to liquid expulsion results in energy dissipation and reduced storage efficiency. In the present work, we use open ensemble (Grand Canonical) Monte Carlo simulations to study the improvement of energy recovery with decreasing diameters of planar pores. Near-complete reversibility is achieved at pore widths barely accommodating a monolayer of the liquid, thus minimizing the area of the liquid/gas interface during the cavitation process. At the same time, these conditions lead to a steep increase in the infiltration pressure required to overcome steric wall/water repulsion in a tight confinement and a considerable reduction in the translational entropy of confined molecules. In principle, similar effects can be expected when increasing the size of the liquid particles without altering the absorbent porosity. While the latter approach is easier to follow in laboratory work, we discuss the advantages of reducing the pore diameter, which reduces the cycling hysteresis while simultaneously improving the stored-energy density in the material.

Wettability of graphene oxide functionalized with N-alkylamines: a molecular dynamics study

The wettability of graphene oxide functionalized with N-alkylamines was studied by molecular dynamics simulations. Six different N-alkylamines and two functionalization degrees were reviewed. The nucleophilic ring-opening reaction mechanism between the N-alkylamines and epoxy functional groups of graphene oxide was considered to generate the atomistic models. Water contact angles increased with both the alkyl chain length and substitution degree. The Wenzel model was used to access the effect of both the surface roughness and alkyl chain length on wettability. The results indicated that functionalization introduces an important increase of surface roughness but its effect on wettability is countered by the alkyl chain length.

Water Adsorption Control by Surface Nanostructures on Graphene-Related Materials by Grand Canonical Monte Carlo Simulations

The interfaces of carbon materials play an important role in various technological and scientific research fields. Graphene is the fundamental unit of sp² carbon allotropes, and the evaluation of the interfacial properties of graphene-related materials is thus essential to clarify the molecular mechanisms occurring at the interfaces. Ideally, graphene is exclusively composed of sp² carbon atoms; however, some parts of graphene normally contain sp³ carbon atoms with oxygen functional groups, vacancy, and grain boundary defects, and these structural characteristics need to be considered to reveal the interfacial properties. Herein, we demonstrate the interfacial properties of graphene-related materials by analyzing the water adsorption properties of graphene, hydrogenated graphene (graphane), and partially oxidized graphene (named as graphoxide) using grand canonical Monte Carlo simulations. The hydrophobicity evaluated from the simulated water adsorption isotherms followed the order: graphane > graphene > graphoxide with 1% oxygen atomic ratio > graphoxide with 3% oxygen atomic ratio > graphoxide with 5% oxygen atomic ratio. The potential calculations between a single water molecule and graphoxides revealed that the presence of oxygen functional groups enhanced the hydrophilicity of graphoxide. This study also disclosed some differences between the hydrophobic interfaces of graphene and graphane, which have been rarely evaluated. Surprisingly, the hydrophobicity of graphane was higher than that of graphene despite the similar potential well depths between a water molecule and graphene/graphane. This was caused by the restriction of water orientation on graphane; water was preferentially adsorbed on the honeycomb center or concave sites in the initial adsorption, and it was hard to interact with neighboring water molecules. The different structures revealed for the graphene-related materials with nanoscale roughness played important roles in controlling the water vapor adsorption mechanism.

Evaluation of the work of adhesion at the interface between a surface-modified metal oxide and an organic solvent using molecular dynamics simulations

Advancing the practical applications of surface-modified nanoparticles requires that their dispersion in solvents can be controlled. The degree of dispersion depends on the affinity between surface-modified nanoparticles and solvents, which can be quantified using the work of adhesion at the interface. Herein, the affinity between a surface-modified inorganic solid and an organic solvent was evaluated by calculating the work of adhesion at the interface. The phantom-wall method, which is a thermodynamic route for evaluating the work of adhesion at an interface using molecular dynamics simulations, was applied to the decanoic acid-modified Al₂O₃/hexane interface. Molecular dynamics simulations were performed for flat interface systems to focus on the interactions between substances that affect the affinity on the surface. As a result, the surface coverage of decanoic acid was found to affect the work of adhesion, with a maximum value of 45.66 ± 0.75 mJ/m² at a surface coverage of 75%. An analysis of the mass density profiles of Al₂O₃, decanoic acid, and hexane in the vicinity of the interface showed that the increase in the work of adhesion with the surface coverage was due to the penetration of hexane molecules into the decanoic acid layer on the Al₂O₃ surface. At a surface coverage of 75%, some hexane molecules were trapped in the layer of oriented decanoic acid molecules. These results suggested that the interfacial affinity can be enhanced by controlling the surface modification so that the solvent can penetrate the layer of the modifier.

Solvent-Solvent Correlations across Graphene: The Effect of Image Charges

Wetting experiments show pure graphene to be weakly hydrophilic, but its contact angle (CA) also reflects the character of the supporting material. Measurements and molecular dynamics simulations on suspended and supported graphene often reveal a CA reduction due to the presence of the supporting substrate. A similar reduction is consistently observed when graphene is wetted from both sides. The effect has been attributed to transparency to molecular interactions across the graphene sheet; however, the possibility of substrate-induced graphene polarization has also been considered. Computer simulations of CA on graphene have so far been determined by ignoring the material's conducting properties. We improve the graphene model by incorporating its conductivity according to the constant applied potential molecular dynamics. Using this method, we compare the wettabilities of suspended graphene and graphene supported by water by measuring the CA of cylindrical water drops on the sheets. The inclusion of graphene conductivity and concomitant polarization effects leads to a lower CA on suspended graphene, but the CA reduction is significantly bigger when the sheets are also wetted from the opposite side. The stronger adhesion is accompanied by a profound change in the correlations among water molecules across the sheet. While partial charges on water molecules interacting across an insulator sheet attract charges of the opposite sign, apparent attraction among like charges is manifested across the conducting graphene. The change is associated with graphene polarization, as the image charges inside the conductor attract equally signed partial charges of water molecules on both sides of the sheet. Additionally, using a nonpolar liquid (diiodomethane), we affirm a detectable wetting translucency when liquid-liquid forces are dominated by dispersive interactions. Our findings are important for predictive modeling toward a variety of applications including sensors, fuel cell membranes, water filtration, and graphene-based electrode materials in high-performance supercapacitors.

Asphaltene Aggregation in Aqueous Solution Using Different Water Models: A Classical Molecular Dynamics Study

The aggregation behavior of asphaltene in aqueous solution is systematically investigated based on a classical molecular dynamics study. In this work, a novel approach is adopted in order to investigate the structural and dynamical properties of the asphaltene nanoaggregates using different water models. The end-to-end distance of the asphaltene molecule is probed in order to understand the aggregation behavior in aqueous solution. The accuracy of different water models, that is, simple point charge, TIP4P-D, and TIP5P, is thoroughly investigated. In order to probe the dynamical properties of the asphaltene nanoaggregates, the transport coefficients, namely, diffusion coefficient and shear viscosity, are computed. The obtained results highlight the importance of using the appropriate water model in order to accurately study the aggregation behavior of asphaltene in aqueous solution.

Molecular polarizability in open ensemble simulations of aqueous nanoconfinements under electric field

Molecular polarization at aqueous interfaces involves fast degrees of freedom that are often averaged-out in atomistic-modeling approaches. The resulting effective interactions depend on a specific environment, making explicit account of molecular polarizability particularly important in solutions with pronounced anisotropic perturbations, including solid/liquid interfaces and external fields. Our work concerns polarizability effects in nanoscale confinements under electric field, open to an unperturbed bulk environment. We model aqueous molecules and ions in hydrophobic pores using the Gaussian-charge-on-spring BK3-AH representation. This involves nontrivial methodology developments in expanded ensemble Monte Carlo simulations for open systems with long-ranged multibody interactions and necessitates further improvements for efficient modeling of polarizable ions. Structural differences between fixed-charge and polarizable models were captured in molecular dynamics simulations for a set of closed systems. Our open ensemble results with the BK3 model in neat-aqueous systems capture the ∼10% reduction of molecular dipoles within the surface layer near the hydrophobic pore walls in analogy to reported quantum mechanical calculations at water/vapor interfaces. The polarizability affects the interfacial dielectric behavior and weakens the electric-field dependence of water absorption at pragmatically relevant porosities. We observe moderate changes in thermodynamic properties and atom and charged-site spatial distributions; the Gaussian distribution of mobile charges on water and ions in the polarizable model shifts the density amplitudes and blurs the charge-layering effects associated with increased ion absorption. The use of polarizable force field indicates an enhanced response of interfacial ion distributions to applied electric field, a feature potentially important for in silico modeling of electric double layer capacitors.

Dynamical insights into the mechanism of a droplet detachment from a fiber

Quantifying the detachment behavior of a droplet from a fiber is important in many applications such as fog harvesting, oil-water separation, or water management in fuel cells. When the droplets are forcibly removed from hydrophilic fibers, the ease of detachment strongly depends on droplet volume and the rate of the process controlled by the applied force. Experiments, conducted on a ferrofluid under magnetic force, as well as continuum level calculations from fluid mechanics have so far been unable to resolve the time-dependent dynamics of droplet detachment and, most importantly, to assess the role of the applied force as the key determinant of the volume of the droplet residue remaining on the fiber after detachment. In the present work, we study the mechanism of water droplet detachment and retention of residual water on smooth hydrophilic fibers using nonequilibrium molecular dynamics simulations. We investigate how the applied force affects the breakup of a droplet and how the minimal detaching force per unit mass decreases with droplet size. We extract scaling relations that allow extrapolation of our findings to larger length scales that are not directly accessible by molecular models. We find that the volume of the residue on a fiber varies nonmonotonically with the detaching force, reaching the maximal size at an intermediate force and associated detachment time. The strength of this force decreases with the size of the drop, while the maximal residue increases with the droplet volume, V, sub-linearly, in proportion to the V2/3.

Fluorographene and Graphane as an Excellent Platform for Enzyme Biocatalysis

The application of enzymes is a crucial issue for current biotechnological application in pharmaceutical, as well as food and cosmetic industry. Effective platforms for enzyme immobilization are necessary for their industrial use in various biosynthesis procedures. Such platforms must provide high yield of immobilization and retain high activity at various conditions for their large-scale applications. Graphene derivatives such as hydrogenated graphene (graphane) and fluorographene can be applied for enzyme immobilization with close to 100 % yield that can result to activities of the composites significantly exceeding activity of free enzymes. The hydrophobic properties of graphene stoichiometric derivatives allowed for excellent non-covalent bonding of enzymes and their use in various organic solvents. The immobilized enzymes retain their high activities even at elevated temperatures. These findings show excellent application potential of enzyme biocatalysts immobilized on graphene stoichiometric derivatives.

Extent of Surface Force Additivity on Chemically Heterogeneous Substrates at Varied Orientations

Surface interactions between chemically mixed surfaces, as well as those among dissolved biomolecules, comprise distinct contributions from polar and hydrophobic moieties. These contributions are often context dependent. Approximate compliance to the Cassie additivity equation for the wetting free energies on mixed surfaces in water is, however, indicative of similarly additive forces between individual surface elements, suggesting a quadratic interpolation model for total force from the forces between pure surfaces. We use molecular dynamics/umbrella sampling simulations of parallel and nonparallel mixed surfaces with demonstrable Cassie-like behavior to verify how well the total surface force between the heterogeneous, molecularly rough surfaces can be approximated as a combination of forces among the homogeneous ones. When accounting for dissimilar distances of approach between functional groups of different types, our results for graphene surfaces with mixed methyl and nitrile coating show such a superposition to provide a reasonable first order approximation of interactions between the platelets. Deviations from additivity are more prominent in parallel-plate configurations, at high content of hydrophobic groups, and small separations. The inclusion of water polarizability does not visibly alter the observed behavior regardless of platelet orientations. The outcome of this study determines the necessary molecular conditions for observing force additivity that emphasize the context dependence of hydrophobic interaction in the presence of polar groups. This notion provides guidelines for the syntheses of new, chemically heterogeneous materials with tailored function-oriented properties in aqueous media.

Interaction between a water drop and holey graphene: retarded imbibition and generation of novel water–graphene wetting states

Water nanodrop imbibition in holey graphene is studied unraveling novel fiber-like wetting state that enhances water–accessible graphene surface area.

Water nanodroplet on a graphene surface-a new old system

The major subject of our study is the accuracy of contact angle calculations. Reporting new simulation data for graphene-water systems, we show that the majority of previously reported data should be treated with caution, since the proper contact angle can be recorded only after a sufficiently long simulation time. It has been proven that-if one wants to gain accuracy greater than 0.1°-long calculations (exceeding 50 ns) are required. Finally, we also show, using both a Groningen Machine for Chemical Simulations (GROMACS) package and our new molecular dynamics (MD) code, that the changes in the contact angle, caused by graphene bottom layer rotation, are within the range of calculation error. We also propose a novel definition of the bottom of the droplet as the height where the density is half the density of liquid water. This new definition is applied in the method of the contact angle calculation from the MD simulation data.

Universal Repulsive Contribution to the Solvent-Induced Interaction Between Sizable, Curved Hydrophobes

In addition to the direct attraction, sizable hydrophobes in water experience an attractive force mediated by interfacial water. Using simple geometric arguments, we identify the conditions at which the water-induced interaction between curved hydrocarbon surfaces becomes repulsive. The repulsive contribution arises from the thermodynamic penalty due to the emergence of the liquid/vapor boundary created as water gets expelled between curved hydrophobes. By augmenting the mean field approach with atomistic simulations of pristine and alkyl-coated graphitic nanoparticles in three distinct geometries, spherical, cylindrical and planar, immersed in water, we show the macroscopic thermodynamics remarkably works down to the molecular scale. The new insights improve the prediction and control of wetting and dispersion properties for a broad class of nonpolar nanoparticles.

Wetting behaviour and contact angles anisotropy of nematic nanodroplets on flat surfaces

We have studied the wetting behaviour of liquid crystal nanodroplets deposited on a planar surface, modelling the mesogens with Gay-Berne ellipsoids and the support surface with a slab of Lennard-Jones (LJ) spherical particles whose mesogen-surface affinity can be tuned. A crystalline and an amorphous planar surface, both showing planar anchoring, have been investigated: the first is the (001) facet of a LJ fcc crystal, the second is obtained from a disordered LJ glass. In both cases we find that the deposited nanodroplet is, in general, elongated and that the contact angle changes around its contour. Simulations for the crystalline substrate show that the angle of contact turns reversibly from anisotropic to isotropic when crossing the clearing transition. As far as we know this is a novel, not yet explored effect for thermotropic liquid crystals, that we hope will stimulate experimental investigations.