skip to main content

This content will become publicly available on March 29, 2023

Title: Organic contaminants and atmospheric nitrogen at the graphene–water interface: a simulation study
Ordered nanoscale patterns have been observed by atomic force microscopy at graphene–water and graphite–water interfaces. The two dominant explanations for these patterns are that (i) they consist of self-assembled organic contaminants or (ii) they are dense layers formed from atmospheric gases (especially nitrogen). Here we apply molecular dynamics simulations to study the behavior of dinitrogen and possible organic contaminants at the graphene–water interface. Despite the high concentration of N 2 in ambient air, we find that its expected occupancy at the graphene–water interface is quite low. Although dense (disordered) aggregates of dinitrogen have been observed in previous simulations, our results suggest that they are stable only in the presence of supersaturated aqueous N 2 solutions and dissipate rapidly when they coexist with nitrogen gas near atmospheric pressure. On the other hand, although heavy alkanes are present at only trace concentrations (micrograms per cubic meter) in typical indoor air, we predict that such concentrations can be sufficient to form ordered monolayers that cover the graphene–water interface. For octadecane, grand canonical Monte Carlo suggests nucleation and growth of monolayers above an ambient concentration near 6 μg m −3 , which is less than some literature values for indoor air. The thermodynamics of more » the formation of these alkane monolayers includes contributions from the hydration free-energy (unfavorable), the free-energy of adsorption to the graphene–water interface (highly favorable), and integration into the alkane monolayer phase (highly favorable). Furthermore, the peak-to-peak distances in AFM force profiles perpendicular to the interface (0.43–0.53 nm), agree with the distances calculated in simulations for overlayers of alkane-like molecules, but not for molecules such as N 2 , water, or aromatics. Taken together, these results suggest that ordered domains observed on graphene, graphite, and other hydrophobic materials in water are consistent with alkane-like molecules occupying the interface. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
Nanoscale Advances
Page Range or eLocation-ID:
1741 to 1757
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract
    <p>This data set for the manuscript entitled &#34;Design of Peptides that Fold and Self-Assemble on Graphite&#34; 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.</p> <p>Version: 2.0</p> <p>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.</p> <p><br /> Conventions Used in These Files<br /> &#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;</p> <p>Structure Files<br /> ----------------<br /> - graph_*.psf or sol_*.psf (original NAMD (XPLOR?) format psf file including atom details (type, charge, mass),More>>
  2. Sum frequency generation (SFG) * Equal contributors. spectroscopy was used to deduce the orientation of the terminal methyl (CH 3 ) group of self-assembled monolayers (SAMs) at the air–solid and air–liquid interfaces at surface concentrations as low as 1% protonated molecules in the presence of 99% deuterated molecules. The SFG spectra of octadecanethiol (ODT) and deuterated octadecanethiol (d 37 ODT) SAMs on gold were used for analysis at the air–solid interface. However, the eicosanoic acid (EA) and deuterated EA (d 39 EA) SAMs on the water were analyzed at the air–liquid interface. The tilt angle of the terminal CH 3 group was estimated to be ∼39 ° for a SAM of 1% ODT : 99% d 37 ODT, whereas the tilt angle of the terminal CH 3 group of the 1% EA : 99% d 39 EA monolayer was estimated to be ∼32 °. The reliability of the orientational analysis at low concentrations was validated by testing the sensitivity of the SFG spectroscopy. A signal-to-noise (S/N) ratio of ∼60 and ∼45 was obtained for the CH 3 symmetric stretch (SS) of 1% ODT : 99% d 37 ODT and 1% EA : 99% d 39 EA, respectively. The estimated increase in S/N ratio values, as a measure of the sensitivitymore »of the SFG spectroscopy, verified the capacity to acquire the SFG spectra at low concentrations of interfacial molecules under ambient conditions. Overall, the orientational analysis of CH 3 SS vibrational mode was feasible at low concentrations of protonated molecules due to increased S/N ratio. In support, the improved S/N ratio on varying incident power density of the visible beam was also experimentally demonstrated.« less
  3. Abstract
    <p>This data set provides files needed to run the simulations described in the manuscript entitled &#34;Organic contaminants and atmospheric nitrogen at the graphene–water interface: A simulation study&#34; using the molecular dynamics software NAMD and LAMMPS. The output of the simulations, as well as scripts used to analyze this output, are also included. The files are organized into directories corresponding to the figures of the main text and supplementary information. They include molecular model structure files (NAMD psf), force field parameter files (in CHARMM format), initial atomic coordinates (pdb format), NAMD or LAMMPS configuration files, Colvars configuration files, NAMD log files, and NAMD output including restart files (in binary NAMD format) and some trajectories in dcd format (downsampled). Analysis is controlled by shell scripts (Bash-compatible) that call VMD Tcl scripts. A modified LAMMPS C&#43;&#43; source file is also included.</p>
    This material is based upon work supported by the 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) andMore>>
  4. The formation of a suite of isoprene-derived hydroxy nitrate (IHN) isomers during the OH-initiated oxidation of isoprene affects both the concentration and distribution of nitrogen oxide free radicals (NOx). Experiments performed in an atmospheric simulation chamber suggest that the lifetime of the most abundant isomer, 1,2-IHN, is shortened significantly by a water-mediated process (leading to nitric acid formation), while the lifetime of a similar isomer, 4,3-IHN, is not. Consistent with these chamber studies, NMR kinetic experiments constrain the 1,2-IHN hydrolysis lifetime to less than 10 s in deuterium oxide (D2O) at 298 K, whereas the 4,3-IHN isomer has been observed to hydrolyze much less efficiently. These laboratory findings are used to interpret observations of the IHN isomer distribution in ambient air. The IHN isomer ratio (1,2-IHN to 4,3-IHN) in a high NOxenvironment decreases rapidly in the afternoon, which is not explained using known gas-phase chemistry. When simulated with an observationally constrained model, we find that an additional loss process for the 1,2-IHN isomer with a time constant of about 6 h best explains our atmospheric measurements. Using estimates for 1,2-IHN Henry’s law constant and atmospheric liquid water volume, we show that condensed-phase hydrolysis of 1,2-IHN can account for this lossmore »process. Simulations from a global chemistry transport model show that the hydrolysis of 1,2-IHN accounts for a substantial fraction of NOxlost (and HNO3produced), resulting in large impacts on oxidant formation, especially over forested regions.

    « less
  5. Abstract

    Spatially resolved tunable diode-laser absorption measurements of the absolute densities of He-I (23S1) metastables in a micro atmospheric pressure plasma jet operated in He/N2and driven by ‘peaks’- and ‘valleys’-type tailored voltage waveforms are presented. The measurements are performed at different nitrogen admixture concentrations and peak-to-peak voltages with waveforms that consist of up to four consecutive harmonics of the fundamental frequency of 13.56 MHz. Comparisons of the measured metastable densities with those obtained from particle-in-cell/Monte Carlo collision simulations show a good quantitative agreement. The density of helium metastables is found to be significantly enhanced by increasing the number of consecutive driving harmonics. Their generation can be further optimized by tuning the peak-to-peak voltage amplitude and the concentration of the reactive gas admixture. These findings are understood based on detailed fundamental insights into the spatio-temporal electron dynamics gained from the simulations, which show that voltage waveform tailoring allows to control the electron energy distribution function to optimize the metastable generation. A high degree of correlation between the metastable creation rate and the electron impact excitation rate from the helium ground state into the He-I ((3s)3S1) level is observed for some conditions which may facilitate an estimation of the metastable densities basedmore »on phase resolved optical emission spectroscopy measurements of the 706.5 nm He-I line originating from the above level and metastable density values at proper reference conditions.

    « less