Energetics of water expulsion from intervening space between two particles during aggregation.

Solvent expulsion away from an intervening region between two approaching particles plays important roles in particle aggregation yet remains poorly understood. In this work, we use metadynamics molecular simulations to study the free energy landscape of removing water molecules from gibbsite and pyrophyllite slit pores representing the confined spaces between two approaching particles. For gibbsite, removing water from the intervening region is both entropically and enthalpically unfavorable. The closer the particles approach each other, the harder it is to expel water molecules. For pyrophyllite, water expulsion is spontaneous, which is different from the gibbsite system. A smaller pore makes the water removal more favorable. When water is being drained from the intervening region, single chains of water molecules are observed in gibbsite pore, while in pyrophyllite pore water cluster is usually observed. Water-gibbsite hydrogen bonds help stabilize water chains, while water forms clusters in pyrophyllite pore to maximize the number of hydrogen bonds among themselves. This work provides the first assessment into the energetics and structure of water being drained from the intervening region between two approaching particles during oriented attachment and aggregation.

Confinement-induced clustering of H2 and CO2 gas molecules in hydrated nanopores.

Gas molecule clustering within nanopores holds significance in the fields of nanofluidics, biology, gas adsorption/desorption, and geological gas storage. However, the intricate roles of nanoconfinement and surface chemistry that govern the formation of gas clusters remain inadequately explored. In this study, through free energy calculation in molecular simulations, we systematically compared the tendencies of H2 and CO2 molecules to aggregate within hydrated hydrophobic pyrophyllite and hydrophilic gibbsite nanopores. The results indicate that nanoconfinement enhances gas dimer formation in the nanopores, irrespective of surface chemistry. However, surface hydrophilicity prohibits the formation of gas clusters larger than dimers, while large gas clusters form easily in hydrophobic nanopores. Despite H2 and CO2 both being non-polar, the larger quadrupole moment of CO2 leads to a stronger preference for dimer/cluster formation compared to H2. Our results also indicate that gases prefer to enter the nanopores as individual molecules, but exit the nanopores as dimers/clusters. This investigation provides a mechanistic understanding of gas cluster formation within nanopores, which is relevant to various applications, including geological gas storage.

Capturing CO2 in Quadrupolar Binding Pockets: Broadband Microwave Spectroscopy of Pyrimidine-(CO2)n, n = 1,2.

Pyrimidine has two in-plane CH(δ+)/N̈(δ-)/CH(δ+) binding sites that are complementary to the (δ-/2δ+/δ-) quadrupole moment of CO2. We recorded broadband microwave spectra over the 7.5-17.5 GHz range for pyrimidine-(CO2)n with n = 1 and 2 formed in a supersonic expansion. Based on fits of the rotational transitions, including nuclear hyperfine splitting due to the two 14N nuclei, we have assigned 313 hyperfine components across 105 rotational transitions for the n = 1 complex and 208 hyperfine components across 105 rotational transitions for the n = 2 complex. The pyrimidine-CO2 complex is planar, with CO2 occupying one of the quadrupolar binding sites, forming a structure in which the CO2 is stabilized in the plane by interactions with the C-H hydrogens adjacent to the nitrogen atom. This structure is closely analogous to that of the pyridine-CO2 complex studied previously by (Doran, J. L. J. Mol. Struct. 2012, 1019, 191-195). The fit to the n = 2 cluster gives rotational constants consistent with a planar cluster of C2v symmetry in which the second CO2 molecule binds in the second quadrupolar binding pocket on the opposite side of the ring. The calculated total binding energy in pyrimidine-CO2 is -13.7 kJ mol-1, including corrections for basis set superposition error and zero-point energy, at the CCSD(T)/ 6-311++G(3df,2p) level, while that in pyrimidine-(CO2)2 is almost exactly double that size, indicating little interaction between the two CO2 molecules in the two binding sites. The enthalpy, entropy, and free energy of binding are also calculated at 300 K within the harmonic oscillator/rigid-rotor model. This model is shown to lack quantitative accuracy when it is applied to the formation of weakly bound complexes.

Theoretical and Experimental Advances in High-Pressure Behaviors of Nanoparticles.

Using compressive mechanical forces, such as pressure, to induce crystallographic phase transitions and mesostructural changes while modulating material properties in nanoparticles (NPs) is a unique way to discover new phase behaviors, create novel nanostructures, and study emerging properties that are difficult to achieve under conventional conditions. In recent decades, NPs of a plethora of chemical compositions, sizes, shapes, surface ligands, and self-assembled mesostructures have been studied under pressure by in-situ scattering and/or spectroscopy techniques. As a result, the fundamental knowledge of pressure-structure-property relationships has been significantly improved, leading to a better understanding of the design guidelines for nanomaterial synthesis. In the present review, we discuss experimental progress in NP high-pressure research conducted primarily over roughly the past four years on semiconductor NPs, metal and metal oxide NPs, and perovskite NPs. We focus on the pressure-induced behaviors of NPs at both the atomic- and mesoscales, inorganic NP property changes upon compression, and the structural and property transitions of perovskite NPs under pressure. We further discuss in depth progress on molecular modeling, including simulations of ligand behavior, phase-change chalcogenides, layered transition metal dichalcogenides, boron nitride, and inorganic and hybrid organic-inorganic perovskites NPs. These models now provide both mechanistic explanations of experimental observations and predictive guidelines for future experimental design. We conclude with a summary and our insights on future directions for exploration of nanomaterial phase transition, coupling, growth, and nanoelectronic and photonic properties.

Control of the Structural Charge Distribution and Hydration State upon Intercalation of CO2 into Expansive Clay Interlayers.

Numerous experimental investigations indicated that expansive clays such as montmorillonite can intercalate CO2 preferentially into their interlayers and therefore potentially act as a material for CO2 separation, capture, and storage. However, an understanding of the energy-structure relationship during the intercalation of CO2 into clay interlayers remains elusive. Here, we use metadynamics molecular dynamics simulations to elucidate the energy landscape associated with CO2 intercalation. Our free energy calculations indicate that CO2 favorably partitions into nanoconfined water in clay interlayers from a gas phase, leading to an increase in the CO2/H2O ratio in clay interlayers as compared to that in bulk water. CO2 molecules prefer to be located at the centers of charge-neutral hydrophobic siloxane rings, whereas interlayer spaces close to structural charges tend to avoid CO2 intercalation. The structural charge distribution significantly affects the amount of CO2 intercalated in the interlayers. These results provide a mechanistic understanding of CO2 intercalation in clays for CO2 separation, capture, and storage.

Hydrophobic Nanoconfinement Enhances CO2 Conversion to H2CO3.

Understanding the formation of H2CO3 in water from CO2 is important in environmental and industrial processes. Although numerous investigations have studied this reaction, the conversion of CO2 to H2CO3 in nanopores, and how it differs from that in bulk water, has not been understood. We use ReaxFF metadynamics molecular simulations to demonstrate striking differences in the free energy of CO2 conversion to H2CO3 in bulk and nanoconfined aqueous environments. We find that nanoconfinement not only reduces the energy barrier but also reverses the reaction from endothermic in bulk water to exothermic in nanoconfined water. Also, charged intermediates are observed more often under nanoconfinement than in bulk water. Stronger solvation and more favorable proton transfer with increasing nanoconfinement enhance the thermodynamics and kinetics of the reaction. Our results provide a detailed mechanistic understanding of an important step in the carbonation process, which depends intricately on confinement, surface chemistry, and CO2 concentration.

Atomistic Mismatch Defines Energy-Structure Relationships during Oriented Attachment of Nanoparticles.

Oriented attachment is an important crystal growth pathway in nature and has been extensively exploited to develop hierarchically structured crystalline materials. Atomistic mismatch in the crystal structure of two particles in the solvent-separated state creates forces that drive particle motions enabling solvent expulsion and coalescence, but the relative magnitudes of the energy barriers for approach, rotation, and translation are not well-known. Here we use classical molecular simulations to calculate the potential of mean force for these three different motions for basal surface encounters of gibbsite nanoplatelets separated by one water layer. In all cases, the highest energy barrier is associated with removing this last water layer to enable jump to contact, even when coaligned. Mutual rotation is more probable than sliding motion, which are both much more probable than jump to contact. This work provides the first comparison on an equal footing of the energy-structure relationships for multiple alignment paths between solvent-separated particles in bulk aqueous solution.

Carbon dioxide-enhanced metal release from kerogen.

Heavy metals released from kerogen to produced water during oil/gas extraction have caused major enviromental concerns. To curtail water usage and production in an operation and to use the same process for carbon sequestration, supercritical CO2 (scCO2) has been suggested as a fracking fluid or an oil/gas recovery agent. It has been shown previously that injection of scCO2 into a reservoir may cause several chemical and physical changes to the reservoir properties including pore surface wettability, gas sorption capacity, and transport properties. Using molecular dynamics simulations, we here demonstrate that injection of scCO2 might lead to desorption of physically adsorbed metals from kerogen structures. This process on one hand may impact the quality of produced water. On the other hand, it may enhance metal recovery if this process is used for in-situ extraction of critical metals from shale or other organic carbon-rich formations such as coal.

Control of Structural Hydrophobicity and Cation Solvation on Interlayer Water Transport during Clay Dehydration.

Swelling clay hydration/dehydration is important to many environmental and industrial processes. Experimental studies usually probe equilibrium hydration states in an averaged manner and thus cannot capture the fast water transport and structural change in interlayers during hydration/dehydration. Using molecular simulations and thermogravimetric analyses, we observe a two-stage dehydration process. The first stage is controlled by evaporation at the edges: water molecules near hydrophobic sites and the first few water molecules of the hydration shell of cations move fast to particle edges for evaporation. The second stage is controlled by slow desorption of the last 1-2 water molecules from the cations and slow transport through the interlayers. The two-stage dehydration is strongly coupled with interlayer collapse and the coordination number changes of cations, all of which depend on layer charge distribution. This mechanistic interpretation of clay dehydration can be key to the coupled chemomechanical behavior in natural/engineered barriers.

Molecular Origin of Wettability Alteration of Subsurface Porous Media upon Gas Pressure Variations.

Upon extraction/injection of a large quantity of gas from/into a subsurface system in shale gas production or carbon sequestration, the gas pressure varies remarkably, which may significantly change the wettability of porous media involved. Mechanistic understanding of such changes is critical for designing and optimizing a related subsurface engineering process. Using molecular dynamics simulations, we have calculated the contact angle of a water droplet on various solid surfaces (kerogen, pyrophyllite, calcite, gibbsite, and montmorillonite) as a function of CO2 or CH4 gas pressure up to 200 atm at a temperature of 300 K. The calculation reveals a complex behavior of surface wettability alteration by gas pressure variation depending on surface chemistry and structure, and molecular interactions of fluid molecules with surfaces. As the CO2 gas pressure increases, a partially hydrophilic kerogen surface becomes highly hydrophobic, while a calcite surface becomes more hydrophilic. Considering kerogen and calcite being the major components of a shale formation, we postulate that the wettability alteration of a solid surface induced by a gas pressure change may play an important role in fluid flows in shale gas production and geological carbon sequestration.

Molecular-level understanding of gibbsite particle aggregation in water.

Using molecular dynamics simulations, we investigate the molecular scale origin of crystal face selectivity when one gibbsite particle attaches to another in water. A comparison of the free energy per unit surface area of particle-particle attachment indicates that particle attachment through edge surfaces, where the edge surfaces are either (1 0 0) or (1 1 0) crystal faces, is more energetically favorable compared to attachment between two basal surfaces (i.e., (0 0 1) crystal faces) or between the basal surface of one particle and the edge surface of another. This result suggests that gibbsite crystals with low basal/edge surface area ratio will preferentially attach through edge surfaces, potentially helping the crystals grow laterally. However, for larger gibbsite particles (high basal/edge surface area ratio) the total free energy, not normalized by surface area, of particle attachment through the basal surfaces is lower (more negative) than attachment through the edge surfaces, indicating that larger gibbsite particles will preferentially aggregate through basal surface attachments. The short-range electrostatic interactions including the interparticle hydrogen bonds from surface -OH groups drive particle attachment, and the dominant contribution to the free energy minimum is enthalpic rather than entropic. However, the enthalpy of basal-edge attachment is significantly offset by the entropy leading to a higher free energy (less negative) compared to that of basal-basal attachment. Study of the free energy for a few imperfect attachments of two particles indicates a higher free energy (i.e., less negative, less stable), compared to a perfect attachment.

Revealing Transition States during the Hydration of Clay Minerals.

A molecular-scale understanding of the transition between hydration states in clay minerals remains a challenging problem because of the very fast stepwise swelling process observed from X-ray diffraction (XRD) experiments. XRD profile modeling assumes the coexistence of multiple hydration states in a clay sample to fit the experimental XRD pattern obtained under humid conditions. While XRD profile modeling provides a macroscopic understanding of the heterogeneous hydration structure of clay minerals, a microscopic model of the transition between hydration states is still missing. Here, for the first time, we use molecular dynamics simulation to investigate the transition states between a dry interlayer, one-layer hydrate, and two-layer hydrate. We find that the hydrogen bonds that form across the interlayer at the clay particle edge make an important contribution to the energy barrier to interlayer hydration, especially for initial hydration.

Enhancement of oil flow in shale nanopores by manipulating friction and viscosity.

Understanding the viscosity and friction of a fluid under nanoconfinement is the key to nanofluidics research. Existing work on nanochannel flow enhancement has been focused on simple systems with only one to two fluids considered such as water flow in carbon nanotubes, and large slip lengths have been found to be the main factor for the massive flow enhancement. In this study, we use molecular dynamics simulations to study the fluid flow of a ternary mixture of octane-carbon dioxide-water confined within two muscovite and kerogen surfaces. The results indicate that, in a muscovite slit, supercritical CO2 (scCO2) and H2O both enhance the flow of octane due to (i) a decrease in the friction of octane with the muscovite wall because of the formation of thin layers of H2O and scCO2 near the surfaces; and (ii) a reduction in the viscosity of octane in nanoconfinement. Water reduces octane viscosity by weakening the interaction of octane with the muscovite surface, while scCO2 reduces octane viscosity by weakening both octane-octane and octane-surface interactions. In a kerogen slit, water does not play any significant role in changing the friction or viscosity of octane. In contrast, scCO2 reduces both the friction and the viscosity of octane, and the enhancement of octane flow is mainly caused by the reduction of viscosity. Our results highlight the importance of multicomponent interactions in nanoscale fluid transport. The results presented here also have a direct implication in enhanced oil recovery in unconventional reservoirs.

Structural Properties of Aqueous Solutions at the (100) and (101) Goethite Surfaces by Molecular Dynamics Simulation.

Trace metal concentrations in soils and sediments are often controlled by adsorption to iron oxides such as goethite in both natural and contaminated systems. Because of goethite's importance as an adsorbent, its interaction with aqueous solutions has been studied extensively. Nonetheless, despite the use of numerous analytical and computational tools, the properties of goethite-aqueous solution interfaces are not fully understood. In this research, we investigate the interaction of water and aqueous NaCl, MgCl2, and BaCl2 solutions ranging in concentration from 0.1 to 4 M, with two goethite surfaces, (100) and (101), using classical molecular dynamics simulation. In the past, the (100) surface has been studied the most because of its simplicity; however, goethite crystals in the environment exhibit other prominent surfaces like the (101) surface which may exhibit very different adsorption properties than the (100) surface. The (100) surface has three surface sites; one is an under-coordinated Fe3+ which interacts strongly with water affecting the interfacial water structure, another site remains deprotonated and forms a hydrogen bond to the only hydroxylated surface site. The (101) surface is terminated with five hydroxyl groups that do not interact as strongly with water as the under-coordinated Fe3+ ion and that form a more corrugated surface structure. As a result, the (101) goethite-solution interface exhibits less water structure, weaker electric double layer oscillations, and more inner-sphere ion adsorption, especially for Cl- and Ba2+ ions. The fundamental differences in interfacial properties for these surfaces suggest that the adsorption properties of one goethite surface cannot be averaged to represent goethite interfaces present in soils and sediments.

Supercritical CO2-induced atomistic lubrication for water flow in a rough hydrophilic nanochannel.

A fluid flow in a nanochannel highly depends on the wettability of the channel surface to the fluid. The permeability of the nanochannel is usually very low, largely due to the adhesion of fluid at the solid interfaces. Using molecular dynamics (MD) simulations, we demonstrate that the flow of water in a nanochannel with rough hydrophilic surfaces can be significantly enhanced by the presence of a thin layer of supercritical carbon dioxide (scCO2) at the water-solid interfaces. The thin scCO2 layer acts like an atomistic lubricant that transforms a hydrophilic interface into a super-hydrophobic one and triggers a transition from a stick- to- a slip boundary condition for a nanoscale flow. This work provides an atomistic insight into multicomponent interactions in nanochannels and illustrates that such interactions can be manipulated, if needed, to increase the throughput and energy efficiency of nanofluidic systems.

Concerted Metal Cation Desorption and Proton Transfer on Deprotonated Silica Surfaces.

The adsorption equilibrium constants of monovalent and divalent cations to material surfaces in aqueous media are central to many technological, natural, and geochemical processes. Cation adsorption-desorption is often proposed to occur in concert with proton transfer on hydroxyl-covered mineral surfaces, but to date this cooperative effect has been inferred indirectly. This work applies density functional theory-based molecular dynamics simulations of explicit liquid water/mineral interfaces to calculate metal ion desorption free energies. Monodentate adsorption of Na+, Mg2+, and Cu2+ on partially deprotonated silica surfaces are considered. Na+ is predicted to be unbound, while Cu2+ exhibits binding free energies to surface SiO- groups that are larger than those of Mg2+. The predicted trends agree with competitive adsorption measurements on fumed silica surfaces. As desorption proceeds, Cu2+ dissociates one of the H2O molecules in its first solvation shell, turning into Cu2+(OH-)(H2O)3, while Mg remains Mg2+(H2O)6. The protonation state of the SiO- group at the initial binding site does not vary monotonically with cation desorption.

Enhanced Ion Adsorption on Mineral Nanoparticles.

Classical molecular dynamics simulation was used to study the adsorption of Na+, Ca2+, Ba2+, and Cl- ions on gibbsite edge (1 0 0), basal (0 0 1), and nanoparticle (NP) surfaces. The gibbsite NP consists of both basal and edge surfaces. Simulation results indicate that Na+ and Cl- ions adsorb on both (1 0 0) and (0 0 1) surfaces as inner-sphere species (i.e., no water molecules between an ion and the surface). Outer-sphere Cl- ions (i.e., one water molecule between an ion and the surface) were also found on these surfaces. On the (1 0 0) edge, Ca2+ ions adsorb as inner-sphere and outer-sphere complexes, whereas on the (0 0 1) surface, outer-sphere Ca2+ ions are the dominant species. Ba2+ ions were found as inner-sphere and outer-sphere complexes on both surfaces. Calculated ion surface coverages indicate that, for all ions, surface coverages are always higher on the basal surface compared to those on the edge surface. More importantly, surface coverages for cations on the gibbsite NP are always higher than those calculated for the (1 0 0) and (0 0 1) surfaces. This enhanced ion adsorption behavior for the NP is due to the significant number of inner-sphere cations found at NP corners. Outer-sphere cations do not contribute to the enhanced surface coverage. In addition, there is no ion adsorption enhancement observed for the Cl- ion. Our work provides a molecular-scale understanding of the relative significance of ion adsorption onto gibbsite basal versus edge surfaces and demonstrates the corner effect on ion adsorption on NPs.

Latanoprost for open-angle glaucoma (UKGTS): a randomised, multicentre, placebo-controlled trial.

BACKGROUND: Treatments for open-angle glaucoma aim to prevent vision loss through lowering of intraocular pressure, but to our knowledge no placebo-controlled trials have assessed visual function preservation, and the observation periods of previous (unmasked) trials have typically been at least 5 years. We assessed vision preservation in patients given latanoprost compared with those given placebo. METHODS: In this randomised, triple-masked, placebo-controlled trial, we enrolled patients with newly diagnosed open-angle glaucoma at ten UK centres (tertiary referral centres, teaching hospitals, and district general hospitals). Eligible patients were randomly allocated (1:1) with a website-generated randomisation schedule, stratified by centre and with a permuted block design, to receive either latanoprost 0·005% (intervention group) or placebo (control group) eye drops. Drops were administered from identical bottles, once a day, to both eyes. The primary outcome was time to visual field deterioration within 24 months. Analyses were done in all individuals with follow-up data. The Data and Safety Monitoring Committee (DSMC) recommended stopping the trial on Jan 6, 2011 (last patient visit July, 2011), after an interim analysis, and suggested a change in primary outcome from the difference in proportions of patients with incident progression between groups to time to visual field deterioration within 24 months. This trial is registered, number ISRCTN96423140. FINDINGS: We enrolled 516 individuals between Dec 1, 2006, and March 16, 2010. Baseline mean intraocular pressure was 19·6 mm Hg (SD 4·6) in 258 patients in the latanoprost group and 20·1 mm Hg (4·8) in 258 controls. At 24 months, mean reduction in intraocular pressure was 3·8 mm Hg (4·0) in 231 patients assessed in the latanoprost group and 0·9 mm Hg (3·8) in 230 patients assessed in the placebo group. Visual field preservation was significantly longer in the latanoprost group than in the placebo group: adjusted hazard ratio (HR) 0·44 (95% CI 0·28-0·69; p=0·0003). We noted 18 serious adverse events, none attributable to the study drug. INTERPRETATION: This is the first randomised placebo-controlled trial to show preservation of the visual field with an intraocular-pressure-lowering drug in patients with open-angle glaucoma. The study design enabled significant differences in vision to be assessed in a relatively short observation period. FUNDING: Pfizer, UK National Institute for Health Research Biomedical Research Centre.

Capacitance enhancement via electrode patterning.

The necessity of increasing the energy density in electric double layer capacitors to meet current demand is fueling fundamental and applied research alike. We report here molecular dynamics simulation results for aqueous electrolytes near model electrodes. Particular focus is on the effect of electrode patterning on the structure of interfacial electrolytes, and on the potential drop between the solid electrodes and the bulk electrolytes. The latter is estimated by numerically integrating the Poisson equation using the charge densities due to water and ions accumulated near the interface as input. We considered uniform and patterned electrodes, both positively and negatively charged. The uniformly charged electrodes are modeled as graphite. The patterned ones are obtained by removing carbon atoms from the top-most graphene layer, yielding nanoscopic squares and stripes patterns. For simplicity, the patterned electrodes are effectively simulated as insulators (the charge remains localized on the top-most layer of carbon atoms). Our simulations show that the patterns alter the structure of water and the accumulation of ions at the liquid-solid interfaces. Using aqueous NaCl solutions, we found that while the capacitance calculated for three positively charged electrodes did not change much, that calculated for the negatively charged electrodes significantly increased upon patterning. We find that both water structure and orientation, as well as ion accumulation affect the capacitance. As electrode patterning affects differently water structure and ion accumulation, it might be possible to observe ion-specific effects. These results could be useful for advancing our understanding of electric double layer capacitors, capacitive desalination processes, as well as of fundamental interfacial electrolytes properties.

Simulation insights for graphene-based water desalination membranes.

Molecular dynamics simulations were employed to study the transport of water and ions through pores created on the basal plane of one graphene sheet (GS). Graphene pore diameters ranged from 7.5 to 14.5 Å. Different pore functionalities were considered, obtained by tethering various functional groups to the terminal carbon atoms. The ease of ion and water translocation across the pores was monitored by calculating the potential of mean force along the direction perpendicular to the GS pore. The results indicate that effective ion exclusion can be achieved only using nonfunctionalized (pristine) pores of diameter ~7.5 Å, whereas the ions can easily penetrate pristine pores of diameters ~10.5 and 14.5 Å. Carboxyl functional groups can enhance ion exclusion for all pores considered, but the effect becomes less pronounced as both the ion concentration and the pore diameter increase. When compared to a carbon nanotube of similar pore diameter, our results suggest that GS pores functionalized with COO(-) groups are more effective in excluding Cl(-) ions from passing through the membrane. Our results suggest that narrow graphene pores functionalized with hydroxyl groups remain effective at excluding Cl(-) ions even at moderate solution ionic strength. The results presented could be useful for the design of water desalination membranes.