Reactive force field study of Li/C systems for electrical energy storage

Evidence and attribution¶

Authority of statements

Prose sections below (Summary, Methods, Findings, etc.) are curated summaries of the publication identified by doi, title, and pdf_path in the front matter above. They are not new primary claims by this wiki.

For definitive numerical values, reaction schemes, and interpretations, use the peer-reviewed article (and optional records under normalized/papers/ when present)—not this page alone.

Summary¶

The work develops a ReaxFF description for Li–C interactions (trained against vdW-corrected DFT) and apply it to graphitic and nanostructured carbon anodes for Li-ion storage. Grand canonical Monte Carlo together with reactive sampling is used to study Li intercalation stages, in-plane ordering, and interlayer spacing in graphite, including vacancy-bearing “defective graphite” where increased Li/C and voltage shifts mirror trends associated with Li plating in recent experiments. Additional geometries (for example onion-like and nanorod-like carbon models) illustrate how dimensionality and defects change rate-capable Li populations at the atomistic level.

Methods¶

Force-field training¶

ReaxFF parameters for Li–C are fit to van der Waals–corrected DFT reference data (abstract), targeting Li–graphite staging motifs and related intercalation geometries summarized in JCTC Methods/tables.

MD / Monte Carlo application¶

Grand canonical Monte Carlo (GCMC) with ReaxFF explores Li loading in 3D periodic graphite supercells of order ~10³ atoms (abstract), including vacancy-bearing “defective graphite” and 0D/1D carbon constructs discussed in the results. After each accepted GCMC trial the authors perform a conjugate-gradient relaxation (0.5 kcal/mol convergence between CG steps) so the host can expand on lithiation. Temperature scans fix an external Li chemical potential (referenced to a BCC Li atom energy) at 300 K until maximum loading. Hybrid MD segments: 12.5 ps NVT molecular dynamics are run for every 50 successful GCMC moves to let C and Li relax between insertions. Metadynamics (used in separate Li diffusion free-energy calculations in the article) adds Gaussian hills (0.5 kcal/mol height, 0.35 Å width) every 100 MD timesteps as stated in the JCTC text. Barostat / NPT: N/A — lithiation is driven by GCMC at fixed Li chemical potential, not hydrostatic pressure control. Replica / applied electric fields: N/A for the summarized bulk intercalation campaign. MD package: the introduction cites LAMMPS and ADF as examples of scalable engines; the Computational Methods sections used here do not attach a single program name to every production run—N/A for a one-line engine label without SI confirmation. MD timestep / thermostat: N/A on this page for a consolidated value—the text quotes 12.5 ps NVT segments and GCMC cadence but scatters other integrator settings.

Findings¶

The parametrized Li/C ReaxFF reproduces intercalation voltage profiles and in-plane ordering / interlayer spacing signatures for stage I and II graphite consistent with experiment and selected DFT references in the paper. Vacancy defects increase effective Li/C and shift voltages in a direction discussed alongside Li plating seen in recent experiments. Nanostructures: a 0D defective onion-like carbon facilitates fast surface adsorption kinetics, whereas a 1D nanorod highlights edge-nucleated intercalation requiring a critical Li density. Sensitivity: staging and rate behavior shift with defect concentration and geometry.

Limitations¶

Electrolyte, SEI, and explicit electron transfer are outside the model. Transferability of the parametrization to very different chemical environments needs separate validation. Detailed GCMC/MD schedules and numerical tables are only in the JCTC article body and SI, not in this short wiki summary.

Relevance to group¶

Foundational ReaxFF + Monte Carlo coupling for battery-relevant carbon lithiation, co-authored with ORNL collaborators; demonstrates the group’s emphasis on scalable reactive models for energy materials.

Citations and evidence anchors¶

Abstract and Sec. I in papers/Raju_LiC_2015_JCTC.pdf summarize parameterization strategy, GCMC/ReaxFF intercalation results, and defective-graphite plating discussion; DOI: 10.1021/ct501027v.

reaxff-family
Pair with paper:2015raju-venue-research (duplicate PDF in corpus; same scientific content).