Molecular dynamics simulations are used to study nonlinear dielectric responses of a confined aqueous film in a planar nanopore under perpendicular electric fields at varied voltages between confining graphene sheets. Dielectric saturation reminiscent of the bulk phase behavior is prevalent at very strong fields, whereas we observe a nonmonotonic permittivity dependence on electric field at intermediate strengths where field-alignment and spontaneous polarization of interfacial water are of comparable magnitudes. The coupling between the two effects results in distinct dielectric responses at opposite confinement walls. The normal component of both the differential dielectric constant and dielectric difference constant tensors averaged over the region closer to the wall under incoming electric field (field pointing from the liquid to the solid phase) initially increases with the strength of the imposed field. The differential permittivity peaks at a field strength previously shown to offset the surface-induced orientation bias of hydration molecules at this wall. Further strengthening of the field results in conventional saturation behavior. At the opposite wall (subject to outgoing field) and in the central region of water slab, the nonlinear dielectric response resembles bulklike saturation. The conditions at the permittivity extremum coincide with the window of accelerated reorientation rates of interfacial water molecules under incoming field uncovered in earlier molecular dynamics analyses.
more »
« less
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 lead 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, by using a non-polar 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.
more »
« less
- Award ID(s):
- 1800120
- NSF-PAR ID:
- 10170653
- Date Published:
- Journal Name:
- ACS Nano
- ISSN:
- 1936-0851
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Rigid wet cooling media is a key component of direct and indirect evaporative cooling systems. Evaporation is the process of a substance in a liquid state changing to a gaseous state. When water evaporates only water molecules get evaporated and the other chemicals in the water are left behind on the surface as residue. Many studies have been conducted on how the change in air flow velocity, media depth, porosity and water distribution affect performance of the cooling system. The operational efficiency of the cooling media varies over its life cycle and depends primarily on temperature and speed of inlet air, water distribution system, type of pad and dimension of the pad.Although evaporative cooling when implemented with air-side economization enables efficiency gains, a trade-off between the system maintenance and its operational efficiency exists. In this study, the primary objective is to determine how calcium scale affects the overall performance of the cooling pad and the water system. Areas of the pad that are not wetted effectively allow air to pass through without being cooled and the edges between wetted and dry surface establish sites for scale formation. An Accelerated Degradation Testing (ADT) by rapid wetting and drying on the media pads at elevated levels of calcium is designed and conducted on the cellulose wet cooling media pad. This research focuses on monitoring the degradation that occurs over its usage and establish a key maintenance parameter for water used in media pad.As a novel study, preliminary tests were mandatory because there were no established standards for media pad degradation testing. Sump water conductivity is identified as the key maintenance parameter for monitoring sump replenishing and draining cycles which will result in reduced water usage. The average water conductivity in the sump during wetting cycles increases monotonically when ADT was performed on a new media pad. An empirical relationship between sump water conductivity and number of wetting cycles is proposed. This information will be very helpful for the manufacturers to guide their customers for maintenance of the media pad and sump water drain cycles.more » « less
-
more » « less
This data set for the manuscript entitled "Design of Peptides that Fold and Self-Assemble on Graphite" includes all files needed to run and analyze the simulations described in the this manuscript in the molecular dynamics software NAMD, as well as the output of the simulations. The files are organized into directories corresponding to the figures of the main text and supporting information. They include molecular model structure files (NAMD psf or Amber prmtop format), force field parameter files (in CHARMM format), initial atomic coordinates (pdb format), NAMD configuration files, Colvars configuration files, NAMD log files, and NAMD output including restart files (in binary NAMD format) and trajectories in dcd format (downsampled to 10 ns per frame). Analysis is controlled by shell scripts (Bash-compatible) that call VMD Tcl scripts or python scripts. These scripts and their output are also included.
Version: 2.0
Changes versus version 1.0 are the addition of the free energy of folding, adsorption, and pairing calculations (Sim_Figure-7) and shifting of the figure numbers to accommodate this addition.
Conventions Used in These Files
===============================Structure Files
----------------
- graph_*.psf or sol_*.psf (original NAMD (XPLOR?) format psf file including atom details (type, charge, mass), as well as definitions of bonds, angles, dihedrals, and impropers for each dipeptide.)- graph_*.pdb or sol_*.pdb (initial coordinates before equilibration)
- repart_*.psf (same as the above psf files, but the masses of non-water hydrogen atoms have been repartitioned by VMD script repartitionMass.tcl)
- freeTop_*.pdb (same as the above pdb files, but the carbons of the lower graphene layer have been placed at a single z value and marked for restraints in NAMD)
- amber_*.prmtop (combined topology and parameter files for Amber force field simulations)
- repart_amber_*.prmtop (same as the above prmtop files, but the masses of non-water hydrogen atoms have been repartitioned by ParmEd)Force Field Parameters
----------------------
CHARMM format parameter files:
- par_all36m_prot.prm (CHARMM36m FF for proteins)
- par_all36_cgenff_no_nbfix.prm (CGenFF v4.4 for graphene) The NBFIX parameters are commented out since they are only needed for aromatic halogens and we use only the CG2R61 type for graphene.
- toppar_water_ions_prot_cgenff.str (CHARMM water and ions with NBFIX parameters needed for protein and CGenFF included and others commented out)Template NAMD Configuration Files
---------------------------------
These contain the most commonly used simulation parameters. They are called by the other NAMD configuration files (which are in the namd/ subdirectory):
- template_min.namd (minimization)
- template_eq.namd (NPT equilibration with lower graphene fixed)
- template_abf.namd (for adaptive biasing force)Minimization
-------------
- namd/min_*.0.namdEquilibration
-------------
- namd/eq_*.0.namdAdaptive biasing force calculations
-----------------------------------
- namd/eabfZRest7_graph_chp1404.0.namd
- namd/eabfZRest7_graph_chp1404.1.namd (continuation of eabfZRest7_graph_chp1404.0.namd)Log Files
---------
For each NAMD configuration file given in the last two sections, there is a log file with the same prefix, which gives the text output of NAMD. For instance, the output of namd/eabfZRest7_graph_chp1404.0.namd is eabfZRest7_graph_chp1404.0.log.Simulation Output
-----------------
The simulation output files (which match the names of the NAMD configuration files) are in the output/ directory. Files with the extensions .coor, .vel, and .xsc are coordinates in NAMD binary format, velocities in NAMD binary format, and extended system information (including cell size) in text format. Files with the extension .dcd give the trajectory of the atomic coorinates over time (and also include system cell information). Due to storage limitations, large DCD files have been omitted or replaced with new DCD files having the prefix stride50_ including only every 50 frames. The time between frames in these files is 50 * 50000 steps/frame * 4 fs/step = 10 ns. The system cell trajectory is also included for the NPT runs are output/eq_*.xst.Scripts
-------
Files with the .sh extension can be found throughout. These usually provide the highest level control for submission of simulations and analysis. Look to these as a guide to what is happening. If there are scripts with step1_*.sh and step2_*.sh, they are intended to be run in order, with step1_*.sh first.
CONTENTS
========The directory contents are as follows. The directories Sim_Figure-1 and Sim_Figure-8 include README.txt files that describe the files and naming conventions used throughout this data set.
Sim_Figure-1: Simulations of N-acetylated C-amidated amino acids (Ac-X-NHMe) at the graphite–water interface.
Sim_Figure-2: Simulations of different peptide designs (including acyclic, disulfide cyclized, and N-to-C cyclized) at the graphite–water interface.
Sim_Figure-3: MM-GBSA calculations of different peptide sequences for a folded conformation and 5 misfolded/unfolded conformations.
Sim_Figure-4: Simulation of four peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) at the graphite–water interface at 370 K.
Sim_Figure-5: Simulation of four peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) at the graphite–water interface at 295 K.
Sim_Figure-5_replica: Temperature replica exchange molecular dynamics simulations for the peptide cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) with 20 replicas for temperatures from 295 to 454 K.
Sim_Figure-6: Simulation of the peptide molecule cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) in free solution (no graphite).
Sim_Figure-7: Free energy calculations for folding, adsorption, and pairing for the peptide CHP1404 (sequence: cyc(GTGSGTG-GPGG-GCGTGTG-SGPG)). For folding, we calculate the PMF as function of RMSD by replica-exchange umbrella sampling (in the subdirectory Folding_CHP1404_Graphene/). We make the same calculation in solution, which required 3 seperate replica-exchange umbrella sampling calculations (in the subdirectory Folding_CHP1404_Solution/). Both PMF of RMSD calculations for the scrambled peptide are in Folding_scram1404/. For adsorption, calculation of the PMF for the orientational restraints and the calculation of the PMF along z (the distance between the graphene sheet and the center of mass of the peptide) are in Adsorption_CHP1404/ and Adsorption_scram1404/. The actual calculation of the free energy is done by a shell script ("doRestraintEnergyError.sh") in the 1_free_energy/ subsubdirectory. Processing of the PMFs must be done first in the 0_pmf/ subsubdirectory. Finally, files for free energy calculations of pair formation for CHP1404 are found in the Pair/ subdirectory.
Sim_Figure-8: Simulation of four peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) where the peptides are far above the graphene–water interface in the initial configuration.
Sim_Figure-9: Two replicates of a simulation of nine peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) at the graphite–water interface at 370 K.
Sim_Figure-9_scrambled: Two replicates of a simulation of nine peptide molecules with the control sequence cyc(GGTPTTGGGGGGSGGPSGTGGC) at the graphite–water interface at 370 K.
Sim_Figure-10: Adaptive biasing for calculation of the free energy of the folded peptide as a function of the angle between its long axis and the zigzag directions of the underlying graphene sheet.
This material is based upon work supported by the US National Science Foundation under grant no. DMR-1945589. A majority of the computing for this project was performed on the Beocat Research Cluster at Kansas State University, which is funded in part by NSF grants CHE-1726332, CNS-1006860, EPS-1006860, and EPS-0919443. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1548562, through allocation BIO200030. -
more » « less
Thermodynamically‐driven exfoliation and self‐assembly of pristine graphene sheets is shown to provide thermally and electrically functional polymer composites. The spreading of graphene sheets at a high energy liquid/liquid interface is driven by lowering the overall energy of the system, and provides for the formation of water‐in‐oil emulsions stabilized by overlapping graphene sheets. Polymerization of the oil phase, followed by removal of the dispersed water phase, produces inexpensive and porous composite foams. Contact between the graphene‐stabilized water droplets provides a pathway for electrical and thermal transport through the composite. Unlike other graphene foams, the graphite used to synthesize these composites is natural flake material, with no oxidation, reduction, sonication, high temperature thermal treatment, addition of surfactants, or high shear mixing required. The result is an inexpensive, low‐density material that exhibits Joule heating and displays increasing electrical conductivity with decreasing thermal conductivity.
-
more » « less
Peptide self-assembly, wherein molecule A associates with other A molecules to form fibrillar β-sheet structures, is common in nature and widely used to fabricate synthetic biomaterials. Selective coassembly of peptide pairs A and B with complementary partial charges is gaining interest due to its potential for expanding the form and function of biomaterials that can be realized. It has been hypothesized that charge-complementary peptides organize into alternating ABAB-type arrangements within assembled β-sheets, but no direct molecular-level evidence exists to support this interpretation. We report a computational and experimental approach to characterize molecular-level organization of the established peptide pair, CATCH. Discontinuous molecular dynamics simulations predict that CATCH(+) and CATCH(−) peptides coassemble but do not self-assemble. Two-layer β-sheet amyloid structures predominate, but off-pathway β-barrel oligomers are also predicted. At low concentration, transmission electron microscopy and dynamic light scattering identified nonfibrillar ∼20-nm oligomers, while at high concentrations elongated fibers predominated. Thioflavin T fluorimetry estimates rapid and near-stoichiometric coassembly of CATCH(+) and CATCH(−) at concentrations ≥100 μM. Natural abundance13C NMR and isotope-edited Fourier transform infrared spectroscopy indicate that CATCH(+) and CATCH(−) coassemble into two-component nanofibers instead of self-sorting. However,13C–13C dipolar recoupling solid-state NMR measurements also identify nonnegligible AA and BB interactions among a majority of AB pairs. Collectively, these results demonstrate that strictly alternating arrangements of β-strands predominate in coassembled CATCH structures, but deviations from perfect alternation occur. Off-pathway β-barrel oligomers are also suggested to occur in coassembled β-strand peptide systems.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1021/acsnano.9b09321
