SOF-8043: Gr/Ni(111) registry and separation simulation notebook - #364
Draft
VsevolodX wants to merge 4 commits into
Draft
SOF-8043: Gr/Ni(111) registry and separation simulation notebook#364VsevolodX wants to merge 4 commits into
VsevolodX wants to merge 4 commits into
Conversation
Reproduces the registry energetics of graphene on Ni(111) from Dahal & Batzill, Nanoscale 6, 2548 (2014): which high-symmetry registry is favourable, and how far the film sits above the surface. Two tiers. The film is placed at each of top-fcc, top-hcp, bridge-top and hollow — sites measured from the substrate's own top three Ni layers, and each registry labelled by where the second carbon sublattice lands — then scanned in z with MACE-MP + D3. A chemisorbing registry has two minima, so the comparison reads the chemisorbed branch and compares each registry at its own minimum; comparing at a shared height misranks them. The platform tier then computes one Total Energy job per registry at that geometry. The structure notebook additionally saves the base interface, which the simulation notebook loads by name: it previously saved only the empirically optimized variant. Verified in JupyterLite: top_fcc wins at 2.01 A (article: top-fcc at 2.1 A) and the hollow registry does not chemisorb. Co-Authored-By: Claude Fable 5 <noreply@anthropic.com>
|
Check out this pull request on See visual diffs & provide feedback on Jupyter Notebooks. Powered by ReviewNB |
The energy-vs-separation figure raised KeyError: 'e_min', a key removed when the scan was reworked into chemisorbed and dispersion-bound branches. Run All Cells continues past an error and the assertions were downstream, so it went unnoticed. Registries now carry the manuscript's own names and cover all four of its Fig. 1 configurations — hollow, atop/fcc, atop/hcp, bridge — with the figure itself embedded. Bridge is defined by its geometry rather than labelled by nearest site: one of its carbons is equidistant from two sites, so classifying it returned whichever the dict happened to list first. Claims match what the evidence supports. The two atop registries differ by a few meV per carbon, finer than this method resolves, so the check is on the atop family rather than on one of the two. The hollow registry's dispersion-bound distance is reported for context, not gated: MACE-MP + D3 places it near 4 A rather than graphite's 3.3 A. Two same-cell reference jobs (bare slab, free-standing film) now give an adsorption energy per carbon atom, with the cell, k-grid, cutoffs and smearing cancelling out of the difference. Also: the displaced variants are no longer written into uploads/, where load_material_from_folder's substring match over sorted filenames made them shadow the base material on a second run; degauss raised to 0.01 Ry for the metal; the scan-edge guard tests the sampled point rather than the interpolated minimum; dead label-mapping block removed; stray tildes in the introduction were rendering as strikethrough. Co-Authored-By: Claude Fable 5 <noreply@anthropic.com>
|
Important Draft PR not reviewedDraft PRs are not automatically reviewed by default.
To automatically review draft PRs, update your CodeRabbit configuration: reviews:
auto_review:
drafts: trueThanks for using CodeRabbit! It's free for OSS, and your support helps us grow. If you like it, consider giving us a shout-out. Comment |
The bridge registry did not match the manuscript's Fig. 1d. The figure puts a first-layer Ni under the midpoint of a C-C bond — the vertical bonds run through the centres of the surface atoms — while the code placed a carbon on the Ni-Ni midpoint, 1.9 A away, which also left that carbon equidistant from the fcc and hcp sites. The placement is now derived from the bond midpoint and verified rather than asserted, and it moves the bridge registry from 95 to 21 meV per carbon above atop/fcc, which is the shallow saddle it should be. starting_magnetization is indexed by position in ATOMIC_SPECIES, so the free-standing graphene reference would have started carbon with nickel's moment. The patch is now built per material by element, and a reference whose elements differ from the interface's gets its own workflow. The adsorption-energy references are off by default: they triple the job count of a run that is meant to finish one job unattended. Cutoffs drop to 40 Ry with an 8x density cutoff, per the GBRV guidelines already followed elsewhere in this repo. The scan-edge warning fired on every run, including where the minimum was properly bracketed by the point below it. It now fires only when the lowest chemisorbed sample is the first in the window, which is the case that actually means the well may lie outside it. Co-Authored-By: Claude Fable 5 <noreply@anthropic.com>
…r a default The density cutoff was 8x the wavefunction cutoff, a ratio taken from a sibling notebook that uses different pseudopotentials for a different system. GBRV publishes its ultrasoft set as a 40 / 200 Ry pair, which is also the platform default, so that is what this uses. Each remaining parameter now states which of the three it rests on. The k-point divisions are a multiple of three because K sits at (1/3, 1/3) and has to lie on the grid, and dense because a metal's Fermi surface needs it. The starting moment is Ni's bulk value. D3 is on because the hollow registry has no chemisorbed minimum at all and is held only by dispersion. The MACE model size is a measurement, not a preference: medium at float32 finds no chemisorbed minimum and inverts the result. The SCF settings are grounded in the failure they fix. A first job stopped at "convergence NOT achieved after 100 iterations" with the total energy oscillating in its fourth decimal — charge sloshing, not divergence. Cold smearing leaves the free energy insensitive to degauss where the gaussian default does not; local-TF mixing is built for the long-wavelength charge oscillation a slab supports; a smaller mixing fraction and more iterations let the magnetic moment settle. Co-Authored-By: Claude Fable 5 <noreply@anthropic.com>
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Specific Example reproducing the graphene/Ni(111) registry energetics of Dahal & Batzill, Graphene–nickel interfaces: a review, Nanoscale 6, 2548 (2014) — SOF-8043, epic SOF-7994.
What it reproduces
top_fccAbsolute adsorption energies are deliberately not compared — the manuscript's come from dispersion-corrected methods beyond semi-local DFT. The notebook compares the registry ordering and the separations, which transfer between methods.
How
Two tiers, so the notebook is useful in minutes and exact when it matters:
total_energyjob per registry at that geometry — PBE,nspin = 2with a starting moment on Ni, smearing, 12×12×1 (multiple of 3 so K is sampled),vdw_corr = 'grimme-d3', identical cutoffs across all four. A default run submits one job and names the registries to activate; the Cypress feature activates all four.Introduction.ipynb's Simulation column is filled for the Gr/Ni row.Change to the structure notebook
It saved only the empirically-optimized variant; it now also saves the base interface, which the simulation notebook loads by name. The simulation notebook raises if that name is missing rather than silently falling back to Standata.
Verification
Fast tier is green end-to-end in JupyterLite, driven by the Cypress feature in the companion web-app PR: provenance, all four registries, all three checks
ok.The DFT tier is not yet verified — every job on the local stack errors with
Account … did not sync to node …. That is an environment fault, not this notebook: mqueue-1'sceleryqueue holds 3249 messages with 0 consumers, the RabbitMQ brokers are not clustered, and the cluster's workers connect to a proxy address that is none of them, so publisher and consumer land on different brokers. The stock Silicontotal_energy.featurefails identically — its spec passes only because it asserts a label. Details inplan/log/platform.md.Draft until the DFT tier runs.
🤖 Generated with Claude Code