skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Atomistic-Scale Simulations on Graphene Bending Near a Copper Surface
Molecular insights into graphene-catalyst surface interactions can provide useful information for the efficient design of copper current collectors with graphitic anode interfaces. As graphene bending can affect the local electron density, it should reflect its local reactivity as well. Using ReaxFF reactive molecular simulations, we have investigated the possible bending of graphene in vacuum and near copper surfaces. We describe the energy cost for graphene bending and the binding energy with hydrogen and copper with two different ReaxFF parameter sets, demonstrating the relevance of using the more recently developed ReaxFF parameter sets for graphene properties. Moreover, the draping angle at copper step edges obtained from our atomistic simulations is in good agreement with the draping angle determined from experimental measurements, thus validating the ReaxFF results.  more » « less
Award ID(s):
1808900
PAR ID:
10251270
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Catalysts
Volume:
11
Issue:
2
ISSN:
2073-4344
Page Range / eLocation ID:
208
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Fabric draping, which is referred to as the process of forming of textile reinforcements over a 3D mold, is a critical stage in composites manufacturing since it determines the fiber orientation that affects subsequent infusion and curing processes and the resulting structural performance. The goal of this study is to predict the fabric deformation during the draping process and develop in-depth understanding of fabric deformation through an architecture-based discrete Finite Element Analysis (FEA). A new, efficient discrete fabric modeling approach is proposed by representing textile architecture using virtual fiber tows modeled as Timoshenko beams and connected by the springs and dashpots at the intersections of the interlaced tows. Both picture frame and cantilever beam bending tests were carried out to characterize input model parameters. The predictive capability of the proposed modeling approach is demonstrated by predicting the deformation and shear angles of a fabric subject to hemisphere draping. Key deformation modes, including bending and shearing, are successfully captured using the proposed model. The development of the virtual fiber tow model provides an efficient method to illustrate individual tow deformation during draping while achieving computational efficiency in large-scale fabric draping simulations. Discrete fabric architecture and the inter-tow interactions are considered in the proposed model, promoting a deep understanding of fiber tow deformation modes and their contribution to the overall fabric deformation responses. 
    more » « less
  2. This is the simulation data set for the manuscript: Arvelo DM, Comer J, Schmit J, Garcia R (2024) Interfacial water is separated from a hydrophobic silica surface by a gap of 1.2 nm. ACS Nano 18:18683–18692. https://doi.org/10.1021/acsnano.4c05689 This data set 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. LAMMPS input files for the ReaxFF simulations are also included. 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 LAMMPS data), force field parameter files (in CHARMM format or ReaxFF format), initial atomic coordinates (pdb format), NAMD or LAMMPS configuration files, Colvars configuration files, NAMD or LAMMPS log files, and output including restart files (in binary NAMD format) and trajectories in dcd format (downsampled with a stride of 100 to 20 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: 1.0. Figure5AC: Simulation of pentadecane on a 5 chains/nm^2 OTS layer. Figure5B_FigureS7: Calculation of force profile for an SiO2 tip asperity model using adaptive biasing force. Systems: octane with 5 chains/nm^2 OTS, octane with 4 chains/nm^2 OTS, decane with 5 chains/nm^2 OTS, water with 5 chains/nm^2 OTS. FigureS6: Simulations showing the effect of octadecane on the structure of the OTS layer for 3 and 5 chains/nm^2 densities. FigureS8: Calculation of the adsorption free energy of tetracosane (C24) at the OTS–water interface using ABF. FigureS9: Python script for estimating the critical concentration to form an alkane layer at the OTS–water interface using the mean-field Ising model. FigureS10: ReaxFF simulation and modeling to create the silanol-terminated amorphous silica model of an AFM tip asperity. FigureS11: Molecular dynamics simulations showing spontaneous assembly of twelve or twenty-four tetracosane (C24) molecules at the interface between water and the alkyl groups of an OTS-conjugated silica surface. 
    more » « less
  3. Subramaniam, B. Executive Editor (Ed.)
    This research presents pioneering work on transforming a variety of waste plastic into synthetic graphite of high quality and purity. Six recycled plastics in various forms were obtained – including reprocessed polypropylene, high-density polyethylene flakes, shredded polyethylene films, reprocessed polyethylene (all obtained from Pennsylvania Recycling Markets Center), polystyrene foams and polyethylene terephthalate bottles (both sourced from a local recycling bin). The waste plastics were carbonized in sealed tubing reactors. The study shows that this versatile process can be used on a mix of waste plastics in a variety of recycled forms to obtain a uniform graphitic carbon phase, hence addressing the challenges of separation and transportation faced by the plastic recycling industry. The conversion yield to elemental carbon for recycled plastics was improved by up to 250% by using graphene oxide (GO) additives. Five different grades of GO and graphene were used to gain insights into the interaction mechanisms between plastics and GO during pyrolysis. The effect of GO additives on carbonization was analyzed using thermogravimetric analysis / differential scanning calorimetry and ReaxFF-based reactive molecular dynamics simulations. The obtained cokes were graphitized at 2500 ℃ and the graphitic quality of the synthetic graphites was analyzed using X-ray diffraction, transmission electron microscopy, and Raman spectroscopy. The plastic waste-derived synthetic graphites exhibit remarkable graphitic quality with crystallite sizes comparable with a model graphitizable material – anthracene coke. The thin, flake-like morphology and nanostructure featuring well-stacked contiguous lamellae make these graphitic carbons highly promising candidates for energy storage applications. Based on our experiments and atomistic-scale simulations we propose interaction mechanisms between the plastic polymers and the graphenic additives that explain the chemical conversion pathways for GO-assisted waste plastic carbonization and graphitization. 
    more » « less
  4. Transition metal dichalcogenides (TMDCs) are potential materials for future optoelectronic devices. Grain boundaries (GBs) can significantly influence the optoelectronic properties of TMDC materials. Here, we have investigated the mechanical characteristics of tungsten diselenide (WSe 2 ) monolayers and failure process with symmetric tilt GBs using ReaxFF molecular dynamics simulations. In particular, the effects of topological defects, loading rates, and temperatures are investigated. We considered nine different grain boundary structures of monolayer WSe 2 , of which six are armchair (AC) tilt structures, and the remaining three are zigzag (ZZ) tilt structures. Our results indicate that both tensile strength and fracture strain of WSe 2 with symmetric tilt GBs decrease as the temperature increases. We revealed an interfacial phase transition for high-angle GBs reduces the elastic strain energy within the interface at finite temperatures. Furthermore, brittle cracking is the dominant failure mode in the WSe 2 monolayer with tilted GBs. WSe 2 GB structures showed more strain rate sensitivity at high temperatures than at low temperatures. 
    more » « less
  5. Bhatele, A.; Hammond, J.; Baboulin, M.; Kruse, C. (Ed.)
    The reactive force field (ReaxFF) interatomic potential is a powerful tool for simulating the behavior of molecules in a wide range of chemical and physical systems at the atomic level. Unlike traditional classical force fields, ReaxFF employs dynamic bonding and polarizability to enable the study of reactive systems. Over the past couple decades, highly optimized parallel implementations have been developed for ReaxFF to efficiently utilize modern hardware such as multi-core processors and graphics processing units (GPUs). However, the complexity of the ReaxFF potential poses challenges in terms of portability to new architectures (AMD and Intel GPUs, RISC-V processors, etc.), and limits the ability of computational scientists to tailor its functional form to their target systems. In this regard, the convergence of cyber-infrastructure for high performance computing (HPC) and machine learning (ML) presents new opportunities for customization, programmer productivity and performance portability. In this paper, we explore the benefits and limitations of JAX, a modern ML library in Python representing a prime example of the convergence of HPC and ML software, for implementing ReaxFF. We demonstrate that by leveraging auto-differentiation, just-in-time compilation, and vectorization capabilities of JAX, one can attain a portable, performant, and easy to maintain ReaxFF software. Beyond enabling MD simulations, end-to-end differentiability of trajectories produced by ReaxFF implemented with JAX makes it possible to perform related tasks such as force field parameter optimization and meta-analysis without requiring any significant software developments. We also discuss scalability limitations using the current version of JAX for ReaxFF simulations. 
    more » « less