Search NASA⌕ Search

SEARCH · Search NASA

Results for “MD simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

Feasibility study of Mg storage in a bilayer silicene anode via application of an external electric field

With the goal of developing a Si-based anode for Mg-ion batteries (MIBs) that is both efficient and compatible with the current semiconductor industry, the current research utilized classical Molecular Dynamics (MD) simulation in investigating the intercalation of a Mg 2+ ion under an external electric field (E-field) in a 2D bilayer silicene anode (BSA). First principles density functional theory calculations were used to validate the implemented EDIP potentials. Our simulation shows that there exists an optimum E-field value in the range of 0.2–0.4 V Å –1 for Mg 2+ intercalation in BSA. To study the effect of the E-field on Mg 2+ ions, an exhaustive spread of investigations was carried out under different boundary conditions, including calculations of mean square displacement (MSD), interaction energy, radial distribution function (RDF), and trajectory of ions. Our results show that the Mg 2+ ions form a stable bond with Si in BSA. The effects of E-field direction and operating temperature were also investigated. In the X–Y plane in the 0°–45° range, 15° from the X-direction was found to be the optimum direction for intercalation. The results of this work also suggest that BSA does not undergo drastic structural changes during the charging cycles with the highest operating temperature being ~300 K

25 ENERGY STORAGE↗

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Macromolecular recognition: Structural aspects of the origin of the genetic system

Theoretical simulation of prebiotic chemical processes is an invaluable tool for probing the phenomenon of the evolution of life. Using computational and modeling techniques and guided by analogies from present day systems, we seek to understand the emergence of the genetic apparatus, enzymatic catalysis and protein synthesis under prebiotic conditions. Modeling of the ancestral aminoacyl-tRNA-synthetases (aRS) may provide important clues to the emergence of the genetic code and the protein synthetic machinery. The minimal structural requirements for the catalysis of tRNA aminoacylation are being explored. A formation of an aminoacyl adenylate was studied in the framework of ab initio molecular orbital theory. The role of individual residues in the vicinity of the TyrRS active site was examined, and the effect of all possible amino acids substitutions near the active site was examined. A formation of aminoacyl tRNA was studied by the molecular modeling system SYBYL with the high resolution crystallographic structures of the present day tRNA, aRS's complexes. The ultimate goal is to propose a simple RNA segment that is small enough to be build in the primordial chemical environment but maintains the specificity and catalytic activity of the contemporary RNA enzyme. To understand the mechanism of ribozyme catalyzed reactions, ab initio and semi-empirical (ZINDO) programs were used to investigate the reaction path of transphosphorylation. A special emphasis was placed on the possible catalytic and structural roles played by the coordinated magnesium cation. Both the inline and adjacent mechanisms of transphosphorylation were studied. The structural characteristics of the target helices, particularly a possible role for the G-T pair, is also studied by a molecular dynamics (MD) simulation technique.

Rein, Robert↗

Strong stretching theory for pH-responsive polyelectrolyte brushes in large salt concentrations

Here, we develop a theory for describing the thermodynamics, configuration, and electrostatics of strongly-stretched, pH-responsive polyelectrolyte (PE) brushes in the presence of large salt concentrations. The aim of the paper, therefore, is to study the properties of a PE brush in a salt concentration regime (namely, large concentrations of several molars) that has been hitherto unexplored theoretically in the context of PE brushes but can be routinely encountered in molecular scale simulations of the problem. The brushes are modelled using our recently developed augmented Strong Stretching Theory (SST), while the effect of the presence of the large salt concentration is accounted for by including the contributions of three different types of non-Poisson–Boltzmann (non-PB) effects to the free energy description of the PE brush induced electric double layer (EDL). These non-PB effects are ionic non-mean-field ion–ion correlations, solvent polarization, and the finite size effect of the ions and water dipoles. We study the individual influences of these different effects and show that the ion–ion correlations and solvent polarization effect reduce the brush height which consequently enhances the monomer density and leads to an electrostatic potential distribution of the brush induced EDL that has a larger magnitude at near-wall locations and becomes zero at shorter distances from the wall. The finite size effect, on the other hand, increases the brush height and therefore, weakens the monomer density and leads to a smaller near-wall magnitude of the EDL potential that becomes zero at larger distances from the wall. In time, we consider the impact of all the three non-PB effects simultaneously and show that the ion–ion correlations and solvent polarization effect dominate the size effects and dictate the overall brush configuration and the EDL electrostatics. We also point out that the influence of all the three non-PB effects becomes the largest for a larger salt concentration and a smaller bulk pH. Finally, we compare our theoretical predictions with those obtained from our recently developed all-atom MD simulation model and obtain an excellent match.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insight into the nanostructure of “water in salt” solutions: A SAXS/WAXS study on imide-based lithium salts aqueous solutions

"Water-in-salt " electrolyte (WISE) series have been studied for lithium-ion batteries and supercapacitor applications, but so far, most of the focus has been on the LiTFSI salt-based systems. Herein, we used small-angle X-ray scattering/wide-angle X-ray scattering (SAXS/WAXS) and molecular dynamics (MD) simulation to investigate the solvation structure of a series of lithium salts with four different symmetric anions: (bis(fluoro sulfonyl)imide (FSI), bis(trifluoromethane sulfonyl)imide (TFSI), bis(pentafluoroethane sulfonyl) imide (BETI) and bis(nonafluorobutane sulfonyl)imide (BNTI)), which have similar parent structures but different lengths of the fluorocarbon chains. Two competing structures were found in the four lithium salts aqueous solutions: anion solvated structure and anion network. The anion network plays a crucial role in obtaining a stable and enlarged electrochemical window. The d spacing of the anion solvated structure follows a linear correlation with the number of carbons in the fluorocarbon chains and an exponential correlation with the concentrations. Further, we also observed that the transition from one structure to another is predominantly controlled by the salt volume fraction for all imide-based aqueous solutions. This work provides a systematic study of the atomistic scale structures of a series of imide-based lithium salt aqueous solutions that could extend the knowledge of WIS electrolytes.

25 ENERGY STORAGE↗

Corrosion resistant hot melt adhesive to bind metals

Hot melt adhesives (HMAs) play an important role in many industries, and their demand is expected to grow. HMAs don't require any solvents, and their application results in the formation of strong bonds with the substrate upon cooling within seconds. These properties differentiate them from liquid glues and make them preferable for practical application. Currently, commercial HMAs are used in bonding lightweight materials such as paper, polymers, and cartons, and have limited usage in areas necessitating the bonding of heavier objects like metals. Here, in this study, we report a design and testing of versatile platform comprising an ion-coordinating polymer and ionic fillers for performance optimization and understanding of structure–property relationships, enabling the rational design of HMAs with improved adhesion to metal surfaces. All-atom Molecular dynamics (MD) simulations and various characterization methods are used to elucidate the adhesion mechanism in model composite system containing polyethylene oxide mixed with chemically diverse salt particles. The maximum adhesion strength is found in composites with Al(OH) 3 and FeCl 3 fillers. Interestingly, the presence of Al(OH) 3 also provides a multifunctional anticorrosion property as measured electrochemically using the Tafel method. The discovered path to formulations with improved adhesion to metal surfaces constitutes an important step toward advancing HMAs for use in the structural and semi-structural metal work domain.

Adhesive and cohesive forces↗

Machine-learning and first-principles investigation of lightweight medium-entropy alloys for hydrogen-storage applications

The transition to a low-carbon economy demands efficient and sustainable energy-storage solutions, with hydrogen emerging as a promising clean-energy carrier and with metal hydrides recognized for their hydrogen-storage capacity. Here, we leverage machine learning (ML) to predict hydrogen-to-metal (H/M) ratios and solution energy by incorporating thermodynamic parameters and local lattice distortion (LLD) as key features. Our best-performing ML model provides improvements to H/M ratios and solution energies over a broad class of medium-entripy alloys (easily extendable to multi-principal-element alloys), such as Ti–Nb-X (X = Mo, Cr, Hf, Ta, V, Zr) and Co–Ni-X (X = Al, Mg, V). Ti–Nb–Mo alloys reveal compositional effects in H-storage behavior, in particular Ti, Nb, and V enhance H-storage capacity, while Mo reduces H/M and hydrogen weight percent by 40–50 %. We attributed results in molybdenum-rich alloys to slow hydrogen kinetics, as validated by our pressure-composition-temperature (PCT) isotherm experiments on pure Ti and Ti 5 Mo 95 alloys. Density functional theory (DFT) and molecular dynamics (MD) simulations also confirm that Ti and Nb promote H diffusion, whereas Mo hinders it, highlighting the interplay between electronic structure, lattice distortions, and hydrogen uptake. Notably, our Gradient Boosting Regression model identifies LLD as a critical factor in H/M predictions. Here, to aid material selection, we present two periodic tables illustrating elemental effects on (a) H 2 wt% and (b) solution energy, derived from ML, and provide a reference for identifying alloying elements that enhance hydrogen solubility and storage.

08 HYDROGEN↗

Antielectrophoretic Response-Driven Bending–Tilting Deformation of Cationic Polyelectrolyte Brushes Drives Nonlinear Electroosmotic Transport in Brush-Grafted Nanochannels

In this paper, we use all-atom molecular dynamics (MD) simulations to describe a non-linearly enhanced electroosmotic (EOS) flow, where, in a nanochannel grafted with cationic PMETAC ([Poly(2-(Methacryloyloxy)Ethyl) Trimethylammonium Chloride]) brushes, a two-fold increase in the electric field strength leads to a several-fold (more than two-fold) increase in the EOS flow strength and volume flow rate. The electric field enforces the PMETAC brushes to undergo a bending-tilting driven deformation with a significant portion of the brush layer becoming parallel to the grafting surface. In response, a substantial fraction of the counterions leave the brush layer (hence become more mobile), but instead of going into the bulk, accumulate at the brush-bulk interface, i.e., stay in proximity of the brush segments aligned parallel to the grafting surface. This creates an interesting situation, where the counterions are not completely within the brush layer, yet they fully screen the brush charges. Such “freer” conditions enable the counterions to achieve very high velocity, thereby ensuring that the water solvating the counterions themselves move very fast triggering the significantly augmented EOS transport. Probing deeper we can identify that the bending-tilting driven brush deformation, enforcing the brushes to align parallel to the substrate, results from the anti-electrophoretic behavior of the brushes, where despite being positively charged, the brushes move against the electric field direction. Such an anti-electrophoretic behavior of the PE brushes, which has not been reported before, can be associated with the very fast velocities of the negatively charged counterions and the electrostatic and hydrodynamic coupling of the counterions with the positive functional groups of the brushes. Here, we anticipate that the findings of this paper will shed light on strategies for nanochannel flow, the anti-electrophoretic response of charged polymer chains, and the significance of capturing the detailed chemical architecture of polyelectrolytes in nanoscale science and engineering.

36 MATERIALS SCIENCE↗

CO 2 in Ionene–Ionic Liquid Composite Membranes

Abstract Ionene – ionic liquid (IL) composites are promising materials for CO 2 separation, yet a molecular‐level understanding of their structure and its impact on CO 2 speciation, solubility, rotation, and diffusivity remains unclear. Herein, using multimodal nuclear magnetic resonance (NMR), time‐of‐flight secondary ion mass spectrometry (ToF‐SIMS), atomic force microscopy (AFM), and molecular dynamics (MD) simulations, we reveal that the composites contain IL‐rich domains extending across hundreds of nanometres within the ionene matrix, and these bicontinuous domains span the entire membrane depth. CO 2 also absorbs into the ionene matrix, with the distribution between two CO 2 species varying with temperature and time. The rotational correlation times of these two species are on the timescale of 0.1 and 1 ns, respectively. As IL content increases, the ionic domains expand, resulting in higher CO 2 solubility due to enhanced molecular dynamics and increased free volume in both ionene backbones and IL‐rich regions. Although CO 2 diffusion in the membranes is an order of magnitude slower than in bulk IL, the activation energy for CO 2 diffusion remains comparable. Ionene‐IL composites represent a promising platform for designing CO 2 separation membranes, offering enhanced CO 2 diffusion and selectivity through IL‐rich domains, and increased CO 2 solubility and mechanical integrity from the ionene matrix.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic Monkey Bar Mechanism of Superionic Li‐ion Transport in LiTaCl 6

Abstract The LiTaCl 6 solid electrolyte has the lowest activation energy of ionic conduction at ambient conditions (0.165 eV), with a record high ionic conductivity for a ternary compound (11 mS cm −1 ). However, the mechanism has been unclear. We train machine‐learning force fields (MLFF) on ab initio molecular dynamics (AIMD) data on‐the‐fly and perform MLFF MD simulations of AIMD quality up to the nanosecond scale at the experimental temperatures, which allows us to predict accurate activation energy for Li‐ion diffusion (at 0.164 eV). Detailed analyses of trajectories and vibrational density of states show that the large‐amplitude vibrations of Cl − ions in TaCl 6 − enable the fast Li‐ion transport by allowing dynamic breaking and reforming of Li−Cl bonds across the space in between the TaCl 6 − octahedra. We term this process the dynamic‐monkey‐bar mechanism of superionic Li + transport which could aid the development of new solid electrolytes for all‐solid‐state lithium batteries.

Lei, Ming↗

Fluorophilic Sigma(σ)‐Lock Self‐Healable Copolymers

Abstract Although F‐Containing molecules and macromolecules are often used in molecular biology to increase the binding with Lewis acidic groups by introducing favorable C−F dipoles, there is virtually no experimental evidence and limited understanding of the nature of these interactions, especially their role in synthetic polymeric materials. These studies elucidate the molecular origin of inter‐ and intra‐Chain interactions responsible for self‐healing of F‐Containing copolymers composed of pentafluorostyrene and n‐butyl acrylate units (p(PFS/nBA). Guided by dynamic surface oscillating force (SOF) and spectroscopic measurements supported by molecular dynamics (MD) simulations, these studies show that the reformation of σ‐σ orbitals in −C−F of PFS and CH 3 CH 2 − of nBA units enables the recovery of entropic energy via fluorophilic‐σ‐lock van der Waals forces when PFS/nBA molar ratios are ~50/50. The strength of these interactions determined experimentally for self‐healable PFS/nBA compositions is in the order ~0.3 kcal/mol which primarily comes from fluorophilic‐σ‐lock (~70 %) contributions. These interactions are significantly diminished for non‐self‐healable counterparts. Strongly polarized −C−F σ orbitals create lateral dipolar forces enhancing the affinity towards −C−H orbitals, facilitating energetically favorable interactions. Entropic recovery driven by non‐Covalent bonding offers a valuable tool in designing materials with unique functionalities, particularly self‐healable batteries and energy storage devices.

Gaikwad, Samruddhi↗

Dynamic Monkey Bar Mechanism of Superionic Li–ion Transport in LiTaCl 6

Here, the LiTaCl 6 solid electrolyte has the lowest activation energy of ionic conduction at ambient conditions (0.165 eV), with a record high ionic conductivity for a ternary compound (11 mS cm –1 ). However, the mechanism has been unclear. We train machine-learning force fields (MLFF) on ab initio molecular dynamics (AIMD) data on-the-fly and perform MLFF MD simulations of AIMD quality up to the nanosecond scale at the experimental temperatures, which allows us to predict accurate activation energy for Li-ion diffusion (at 0.164 eV). Detailed analyses of trajectories and vibrational density of states show that the large-amplitude vibrations of Cl – ions in TaCl 6 – enable the fast Li-ion transport by allowing dynamic breaking and reforming of Li–Cl bonds across the space in between the TaCl 6 – octahedra. We term this process the dynamic-monkey-bar mechanism of superionic Li + transport which could aid the development of new solid electrolytes for all-solid-state lithium batteries.

25 ENERGY STORAGE↗

Fluorophilic Sigma(σ)‐Lock Self‐Healable Copolymers

Although F-Containing molecules and macromolecules are often used in molecular biology to increase the binding with Lewis acidic groups by introducing favorable C−F dipoles, there is virtually no experimental evidence and limited understanding of the nature of these interactions, especially their role in synthetic polymeric materials. These studies elucidate the molecular origin of inter- and intra-Chain interactions responsible for self-healing of F-Containing copolymers composed of pentafluorostyrene and n-butyl acrylate units (p(PFS/nBA). Guided by dynamic surface oscillating force (SOF) and spectroscopic measurements supported by molecular dynamics (MD) simulations, these studies show that the reformation of σ-σ orbitals in −C−F of PFS and CH 3 CH 2 − of nBA units enables the recovery of entropic energy via fluorophilic-σ-lock van der Waals forces when PFS/nBA molar ratios are ~50/50. The strength of these interactions determined experimentally for self-healable PFS/nBA compositions is in the order ~0.3 kcal/mol which primarily comes from fluorophilic-σ-lock (~70 %) contributions. These interactions are significantly diminished for non-self-healable counterparts. Strongly polarized −C−F σ orbitals create lateral dipolar forces enhancing the affinity towards −C−H orbitals, facilitating energetically favorable interactions. Entropic recovery driven by non-Covalent bonding offers a valuable tool in designing materials with unique functionalities, particularly self-healable batteries and energy storage devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis of Gold Hydride at High Pressure and High Temperature

Gold is an unreactive metal and its chemical interactions with hydrogen have only recently been explored. Here, in this study, we report the formation of gold hydride above 40 GPa and 2200 K in X-ray free electron laser heated diamond anvil cells using various hydrocarbons as hydrogen sources. Above 40 GPa, a hexagonal phase emerges close to the gold melting point, corresponding to a hydride with stoichiometry Au 2 H x , with x increasing from 0 to near 1 with pressure from 40 to 80 GPa. This is a high-temperature phase which reverts to face centered cubic gold on cooling to 295 K. Accompanying DFT-MD simulations are in excellent agreement with experiment and reveal the structure to consist of an hexagonal close packed gold lattice with atomic hydrogen disordered in the interstices. The hydrogen is superionic and exhibits high diffusivity through the crystalline gold lattice. Our results present the first solid-state binary compound of gold and hydrogen.

Ab initio calculations↗

Self‐Healable Poly(ionic liquid) Copolymers Driven by Polar and Dipolar Forces

Abstract Commodity aliphatic and aromatic acrylic‐based copolymers self‐heal due to ubiquitouskey‐and‐lock,ring‐and‐lock, andfluorophilic‐σ‐lockvan der Waals (vdW) interactions. However, the role of these interactions in the presence of covalently copolymerized ionic liquid (IL) is not known. This study is driven by the hypothesis that covalently incorporated cation–anion pairs to form poly(ionic liquid) copolymers (PILCs) can perturb inter‐ or intra‐chain vdW interactions reflected in mechanical and electrical responses. To test this hypothesis, we synthesized a series of PILCs comprising of pentafluorostyrene (PFS) and imidazolium‐based IL monomers with variable‐length aliphatic tails (methyl and butyl). Using a combination of 2D 1 H‐ 1 H and 19 F ‐ 19 F NOESY NMR and FTIR measurements supplemented by molecular dynamic (MD) simulations, these studies demonstrate that preferentially alternating/random PILCs topologies facilitate self‐healing. The introduction of cation–anion moieties modifies thefluorophilic‐σ‐lockinteractions and, along with longer aliphatic tails ─(CH 2 ) 3 CH 3 covalently attached to the imidazolium cation, enhances cation‐anion mobility, thus faster recovery from mechanical damage occurs. These findings underline how precise control over dipolar and ionic interactions through copolymer composition enables self‐healing in PILCs. These insights may open pathways for designing sustainable, mechanically resilient materials for applications in energy storage and energy harvesting.

Chemistry↗

Theoretical Investigation of the Oxygen Reduction Reaction over Platinum Catalysts Supported by Multi-Edged Vertically Aligned Carbon Nanofiber for Electrocatalyst Preparation

Vertically aligned carbon nanofibers (VACNFs) are promising supports for oxygen reduction reaction (ORR) electrocatalysts in fuel cells. Although experimentally these catalytic systems have shown great potential, there is lack of molecular understanding of the catalytic sites and reaction mechanisms. This work investigated the origin of the ORR reactivities of the platinum catalysts on multi-edged VACNFs (Pt/VACNF) using a multiscale modeling approach combining Density Functional Theory (DFT) and classical Molecular Dynamics (MD) simulations. Based on the ReaxFF potential, here all nanoscale Pt particles (Pt 55 , P 100 , and Pt 147 ) are stabilized by the open edges located axially along the VACNF walls. The calculated first-shell coordination numbers, CN, of surface Pt atoms are 6.63, 7.27, and 7.85, respectively, suggesting that the percentage of low coordination sites increases as the particle size decreases. The adsorption energies of OOH, O, and OH on Pt-55 were systematically probed using DFT calculations. These adsorption energies retain a linear correlation against the generalized coordination numbers ($\overline{CN}$). For Pt nanoparticles supported on VACNF, we found that the OOH and OH bind stronger than on Pt (111) by 0.14 and 0.17 eV, respectively, which can hinder the ORR activity with lower limiting potential than Pt (111). Our theoretical prediction is in good agreement with the linear sweeping voltammetry that revealed a left shift of the half-wave potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure- and Interaction-Based Design of Anti-SARS-CoV-2 Aptamers

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. In this study, we aim to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants.

60 APPLIED LIFE SCIENCES↗

Direct Synthesis of Layer-Tunable and Transfer-Free Graphene on Device-Compatible Substrates Using Ion Implantation Toward Versatile Applications

Direct synthesis of layer-tunable and transfer-free graphene on technologically important substrates is highly valued for various electronics and device applications. State of the art in the field is currently a two-step process: a high-quality graphene layer synthesis on metal substrate through chemical vapor deposition (CVD) followed by delicate layer transfer onto device-relevant substrates. Here, we report a novel synthesis approach combining ion implantation for a precise graphene layer control and dual-metal smart Janus substrate for a diffusion-limiting graphene formation to directly synthesize large area, high quality, and layer-tunable graphene films on arbitrary substrates without the post-synthesis layer transfer process. Carbon (C) ion implantation was performed on Cu–Ni film deposited on a variety of device-relevant substrates. A well-controlled number of layers of graphene, primarily monolayer and bilayer, is precisely controlled by the equivalent fluence of the implanted C-atoms (1 monolayer ~4 × 10 15 C-atoms/cm 2 ). Upon thermal annealing to promote Cu-Ni alloying, the pre-implanted C-atoms in the Ni layer are pushed toward the Ni/substrate interface by the top Cu layer due to the poor C-solubility in Cu. As a result, the expelled C-atoms precipitate into a graphene structure at the interface facilitated by the Cu-like alloy catalysis. After removing the alloyed Cu-like surface layer, the layer-tunable graphene on the desired substrate is directly realized. The layer-selectivity, high quality, and uniformity of the graphene films are not only confirmed with detailed characterizations using a suite of surface analysis techniques but more importantly are successfully demonstrated by the excellent properties and performance of several devices directly fabricated from these graphene films. Molecular dynamics (MD) simulations using the reactive force field (ReaxFF) were performed to elucidate the graphene formation mechanisms in this novel synthesis approach. With the wide use of ion implantation technology in the microelectronics industry, this novel graphene synthesis approach with precise layer-tunability and transfer-free processing has the promise to advance efficient graphene-device manufacturing and expedite their versatile applications in many fields.

36 MATERIALS SCIENCE↗