BiopowerLab

1. True Magenta — NV⁻ + SiV⁻ Co-Doped Diamond

Target emission: Dual peaks at 637 nm (NV⁻) and 738 nm (SiV⁻) → additive color mixing perceived as magenta.

Minimum carrier scale: Single nanodiamond ≥5 nm hosts one NV⁻; co-locating SiV⁻ requires ≥15 nm. Practical co-doped particles ≥20 nm.

Thermodynamic feasibility:

• N substitution in diamond: formation energy E_f(N_s) ≈ 3.0 eV — favorable at HPHT conditions where kT ≈ 0.14 eV (1400°C) is compensated by the ~6 GPa pressure driving C(graphite)→C(diamond) conversion (ΔG ≈ −2.9 kJ/mol at 6 GPa).
• Si split-vacancy: E_f(SiV) ≈ 5.5 eV in neutral state — requires vacancy supply from irradiation; formation is kinetically limited, not thermodynamically prohibited.
• NV⁻ charge state: stabilized when Fermi level E_F > NV(−/0) transition level at E_V + 2.6 eV — achieved with O-terminated surfaces (electron affinity ~1.7 eV) or with neighboring N donors.

Route A — HPHT Co-Doped Growth:

1. Carbon source preparation: Solid High-purity graphite (99.99% C, <1 ppm metallic impurity) crushed to 200-mesh powder. Dried at 200°C in Ar for 12h to remove adsorbed moisture.
   C(graphite, 200-mesh) — stored under Ar(g) at 1 atm, 25°C
2. Catalyst alloy preparation: Solid Fe-Ni-Co catalyst (64:28:8 wt%) arc-melted under Ar, with Si powder (99.999%) blended at 0.05–0.15 wt%. Compressed into pellets at 500 MPa.
   Fe(s) + Ni(s) + Co(s) → [Fe₆₄Ni₂₈Co₈](alloy, s)
   Si(s, 99.999%) → mechanically blended into catalyst pellet
3. N₂ atmosphere loading: Gas HPHT capsule (Re-lined Mo) sealed with N₂ partial pressure 0.3–0.8 atm to control N incorporation at 50–200 ppm in the grown crystal.
   N₂(g, 0.5 atm) ⇌ 2N(dissolved in Fe-Ni-Co melt, l)
   Sieverts' law: [N] ∝ √(P_N₂); at 0.5 atm and 1400°C, [N]_melt ≈ 0.02 wt%.
4. HPHT crystallization: Solid + Liquid Belt press or multi-anvil: ramp to 5.5–6.0 GPa over 30 min, heat to 1350–1500°C, hold 24–72h. Carbon dissolves into molten catalyst, supersaturates, and precipitates as diamond on seed crystal. Si and N incorporate substitutionally during growth.
   C(graphite, s) →[dissolves into Fe-Ni-Co(l) at P≥5.5 GPa]→ C(diamond, s)
   N(dissolved, l) → N_s(substitutional in diamond, s)
   Si(dissolved, l) → Si_s(substitutional in diamond, s)
   Crystallization rate: ~1–5 mg/h depending on ΔT between dissolution and growth zones (~30–50°C gradient). {111} faces grow fastest.
5. Catalyst removal — acid dissolution: Liquid Boil in concentrated HCl:HNO₃ (3:1, aqua regia) at 120°C for 48h, then H₂SO₄:HClO₄ (3:1) at 250°C for 24h.
   Fe(s) + 4HCl(aq) + HNO₃(aq) → FeCl₃(aq) + NO(g)↑ + 2H₂O(l)
   Ni(s) + 2HCl(aq) → NiCl₂(aq) + H₂(g)↑
   Co(s) + 2HCl(aq) → CoCl₂(aq) + H₂(g)↑
   Graphitic carbon: C(sp²) + HClO₄(aq) →[250°C]→ CO₂(g)↑ + HCl(aq)
   Double replacement: metal enters solution as chloride, acid anion replaces lattice.
6. Vacancy creation — electron irradiation: Solid 2 MeV electron beam, fluence 1×10¹⁸ e⁻/cm² (beam current ~50 µA, 6h exposure). Creates Frenkel pairs uniformly throughout bulk.
   C(lattice, s) + e⁻(2 MeV) → V(vacancy) + C_i(interstitial)
   Displacement threshold E_d ≈ 37–47 eV for C in diamond
   Each 2 MeV electron displaces ~10 C atoms along its track; at 10¹⁸ e⁻/cm², vacancy concentration ≈ 10¹⁹ cm⁻³.
7. Vacancy migration anneal: Solid in Vacuum (<10⁻⁵ mbar). Ramp: RT→400°C at 5°C/min (interstitials recombine), hold 1h; ramp 400→800°C at 2°C/min, hold 2h (vacancies mobile, E_a(V) ≈ 2.3 eV, diffusion length ~50 nm at 800°C/2h).
   V(mobile at 800°C) + N_s(stationary) → NV (nitrogen-vacancy center)
   V(mobile at 800°C) + Si_s(stationary) → SiV (silicon split-vacancy, D₃d)
   Competing reaction: V + V → V₂ (divacancy, undesired) — minimized by keeping [V] < [N]+[Si].
8. Charge-state stabilization: Plasma + Solid Oxygen-terminate surface via O₂ plasma (300 W, 5 min, 0.5 Torr) to create negative electron affinity, stabilizing NV⁻ over NV⁰. SiV⁻ is intrinsically stable when N donors are present (electron transfer N→SiV).
   Diamond-H(surface) + O₂(plasma) → Diamond-O(surface) + H₂O(g)↑
   NV⁰ + e⁻(from N donor or O-surface band bending) → NV⁻
   SiV⁰ + e⁻(from N donor) → SiV⁻

Route B — CVD Growth with Dual Precursors (alternative):

1. Gas Microwave plasma CVD: CH₄(2%)/H₂ + N₂(200 ppm) + SiH₄(50 ppm) at 850°C, 40 Torr, 1.5 kW. Growth rate ~1 µm/h on (100) single-crystal seed.
   CH₄(g) + H₂(g) →[plasma, 850°C]→ C·(radical) + H·(radical) + H₂(g)
   C·(radical) → C(diamond surface, s) — step-flow growth on (100)
   N₂(g) →[plasma dissociation]→ 2N· → N_s(in growing lattice)
   SiH₄(g) →[plasma]→ Si· + 2H₂(g) → Si_s(in growing lattice)
2. Solid Irradiation and anneal as Route A, steps 6–8.

Route C — Detonation Nanodiamond + Ion Implantation (nano-scale):

1. Solid + Gas Detonation nanodiamond (DND, 4–5 nm) from TNT/RDX detonation in inert atmosphere. Purify in boiling HClO₄/H₂SO₄ (1:3) for 48h.
   C₇H₅N₃O₆(TNT, s) →[detonation, ~3000°C, ~30 GPa, µs]→ C(diamond, s, 4–5 nm) + CO₂(g) + H₂O(g) + N₂(g)
2. Solid Spin-coat DND monolayer on SiO₂ substrate. Implant ²⁸Si⁺ at 30 keV (10¹³/cm²) — range ~15 nm, matching ND core. Implant ¹⁴N⁺ at 10 keV (10¹³/cm²) — range ~8 nm.
   ²⁸Si⁺(30 keV) + C(ND lattice) → Si_s + V + C_i (implant damage)
   ¹⁴N⁺(10 keV) + C(ND lattice) → N_s + V + C_i
3. Solid Anneal at 800°C, 2h, in forming gas → V migrates to Si and N. O₂ plasma for charge stability. Lift off substrate by sonication in DI water.

Environmental conditions matrix:

Verification protocol: Confocal PL at 532 nm excitation; expect peaks at 575 (NV⁰, weak), 637 (NV⁻, strong), 738 (SiV⁻, narrow). Hanbury Brown–Twiss g²(0) measurement on single particles to confirm single-photon emission from each center independently.

2. Broadband White — Multi-Center Ensemble Diamond

Target emission: Simultaneous N3 (415 nm blue), H3 (503 nm green), NV⁰ (575 nm yellow-orange), NV⁻ (637 nm red) → additive mixing to perceived white across CIE 1931 chromaticity.

Minimum carrier scale: ≥50 nm nanodiamond for sufficient defect diversity; ≥100 nm for balanced multi-center population.

Thermodynamic feasibility:

• N aggregation sequence: N_s(isolated, C-center) →[>700°C]→ N₂(A-aggregate) →[>1400°C, ~Gyr or HPHT hours]→ N₃V(B-aggregate) + N₃(N3 center). A→B transition activation energy E_a ≈ 5.0–5.5 eV; requires HPHT to complete in practical timescales.
• H3 formation: N₂(A-agg.) + V →[anneal 600–800°C]→ NVN(H3). Exothermic: ΔE ≈ −1.5 eV once V is mobile.
• Challenge: all four centers at balanced intensities requires spatial partitioning — high-N zones for N3/H3, low-N zones for NV — hence graded doping.

Route A — Graded-N CVD + Sequential Irradiation/Anneal:

1. Substrate preparation: Solid Type IIa HPHT seed (100)-oriented, polished Ra < 5 nm. Acid clean: H₂SO₄:H₂O₂ (4:1) boil 1h, then HF dip 30s.
   Surface contaminants → dissolved by H₂SO₄/H₂O₂ (piranha)
   SiO₂(native) + 6HF(aq) → H₂SiF₆(aq) + 2H₂O(l)
2. Layer 1 — High-N CVD (N3/H3 precursor zone): Gas CH₄(3%)/H₂ + N₂(2000 ppm), 900°C, 50 Torr, microwave 2.0 kW. Grow 30 µm. Incorporates [N] ≈ 100–200 ppm → forms A-aggregates during growth.
   CH₄(g) →[plasma]→ C(diamond, s) + 2H₂(g)
   N₂(g) →[plasma]→ 2N· → N_s(in lattice) →[growth T 900°C]→ partial N₂ aggregation
3. First irradiation — proton beam: Solid 300 keV H⁺, 5×10¹⁶ H⁺/cm². Bragg peak at ~2 µm depth. Creates vacancy-rich zone in top of Layer 1.
   H⁺(300 keV) + C(lattice) → V + C_i (primary knock-on)
   V + N-N(A-agg.) →[600°C anneal]→ NVN (H3 center, 503 nm green)
4. First anneal: Solid in Vacuum 600°C, 2h → V mobile (E_a(V) ≈ 2.3 eV, just becoming mobile at 600°C). V captured by A-aggregates → H3. Interstitials recombine or cluster.
   V(mobile) + N₂(A-aggregate) → H3 (NVN, s)
5. Layer 2 — Low-N CVD (NV precursor zone): Gas CH₄(1.5%)/H₂ + N₂(50 ppm), 850°C, 40 Torr, 1.5 kW. Grow 20 µm. Low N ensures isolated substitutional N (C-centers), not aggregates.
   CH₄(g) + H₂(g) + N₂(trace) → C(diamond):N_s(isolated, ~20 ppm)
6. Second irradiation — electron beam: Solid 2 MeV e⁻, 1×10¹⁸ e⁻/cm². Uniform vacancy creation through both layers.
   e⁻(2 MeV) + C(lattice) → V + C_i (uniform through depth)
7. Second anneal: Solid in Vacuum 800°C, 2h → V migrates to isolated N in Layer 2 → NV centers. In Layer 1, additional V are captured by remaining A-aggregates → more H3, or form NV at isolated N sites → NV⁰ (weaker).
   V + N_s(isolated, Layer 2) → NV (575/637 nm)
   V + N₂(remaining A-agg., Layer 1) → H3 (additional, 503 nm)
8. HPHT aggregation treatment for N3: Solid 1600°C, 6.0 GPa, 30 min in Ar-sealed capsule. Drives partial A→B aggregation in the high-N Layer 1. Creates N3 (N₃V) centers.
   3N_s(isolated or A-agg.) + V →[1600°C, 6 GPa]→ N₃V (N3 center, 415 nm blue)
   Single replacement: V displaces C at N-cluster site → N₃V
9. Surface treatment: Plasma O₂ plasma 300 W, 5 min → stabilizes NV⁻ charge state. NV⁰ persists in N-poor regions (desired for 575 nm yellow contribution).
   Diamond-H(s) + O·(plasma) → Diamond-OH(s) → Diamond=O(s) + H₂O(g)

Route B — Single-Crystal with Controlled N Gradient (alternative):

1. Solid HPHT growth with N₂ pressure modulated during growth: start at 1.0 atm N₂ (high-N zone), reduce to 0.05 atm over 48h (low-N zone). Creates radial N gradient in one crystal.
2. Solid Electron irradiation + 800°C anneal → NV in low-N zone, H3 in high-N zone.
3. Solid HPHT re-treatment 1600°C 6 GPa 20 min → N3 in high-N zone.

Crystallization kinetics: CVD diamond growth at 900°C on (100) face proceeds by step-flow mechanism. Growth rate R = k·[CH₃·]·exp(−E_a/kT) where E_a ≈ 0.7 eV for H-abstraction rate-limiting step. At 3% CH₄, R ≈ 3–5 µm/h. N incorporation efficiency η_N ≈ 10⁻³ (1 N per 1000 C atoms in gas → 1 N per 10⁶ C in lattice at these conditions).

Environmental conditions matrix:

Verification: PL mapping with 405 nm (excites N3, H3), 532 nm (excites NV), and 660 nm (selectively excites SiV if present) excitation lasers. CIE chromaticity analysis of total emission spectrum should fall within 0.01 of D65 white point (x=0.3127, y=0.3290).

3. Turquoise / Teal — Vacancy-Boron-Nitrogen Complex (VBN)

Target emission: ~490–510 nm from a ternary defect complex with simultaneous B and N near-neighbor substitution and an adjacent vacancy.

Minimum carrier scale: Single defect occupies ~3 lattice sites → ~1 nm. Host crystal ≥10 nm for quantum confinement not to shift energy levels.

Thermodynamic feasibility:

• B substitution: E_f(B_s) ≈ 1.1 eV — readily incorporated. Creates acceptor level 0.37 eV above VB.
• N substitution: E_f(N_s) ≈ 3.0 eV — creates donor level 1.7 eV below CB.
• B-N compensation: when B and N are nearest-neighbors, they form a donor-acceptor pair (DAP) with net charge compensation. Predicted emission from VBN ternary: ~2.5 eV (496 nm) based on DAP energy minus vacancy relaxation (~0.3 eV).
• Challenge: B and N compete for same substitutional sites; at HPHT, B preferentially enters lattice on {111} growth sectors, N on {100}. Must force co-location.

Route A — HPHT with h-BN Decomposition:

1. Precursor preparation: Solid Hexagonal boron nitride (h-BN) powder (99.5%) ball-milled with graphite (99.99%) at 1:500 mass ratio. This ensures atomic-scale mixing of B and N with C.
   BN(hexagonal, s) + C(graphite, s) →[ball mill, 24h, WC vial, Ar atm]→ C:BN(intimately mixed powder, s)
2. HPHT crystallization: Solid + Liquid Fe-Ni catalyst (70:30) + BN-doped graphite. 6.0 GPa, 1400°C, 48h. At these conditions, h-BN decomposes:
   BN(s) →[6 GPa, 1400°C]→ B(dissolved in Fe-Ni melt, l) + N(dissolved in Fe-Ni melt, l)
   C(graphite, s) →[dissolves into melt]→ C(diamond, s) incorporating B_s and N_s
   Net: C(graphite) + BN(s) →[catalyst, HPHT]→ diamond:(B_s,N_s)(s)
   Key: BN decomposition at 6 GPa releases both atoms into melt simultaneously → increases probability of adjacent-site incorporation.
3. Vacancy creation: Solid He⁺ irradiation at 150 keV, 10¹⁵ He⁺/cm². He implants at ~500 nm depth; creates ~5 vacancies per He ion. Post-implant, He diffuses out at >600°C (interstitial He in diamond, E_a(He diffusion) ≈ 0.3 eV).
   He⁺(150 keV) + C(lattice) → V + C_i + He(interstitial)
   At 700°C: He(interstitial) → He(g) (diffuses to surface and escapes)
4. VBN formation anneal: Solid in Gas (95% Ar / 5% H₂ forming gas). 700°C, 1h. Vacancies mobile; captured by B-N pairs to form VBN ternary.
   V(mobile) + B_s(near) + N_s(adjacent to B) → V-B-N (ternary complex, s)
   Forming gas: H₂ prevents surface oxidation; maintains stable surface Fermi level

Route B — Sequential Ion Implantation into Type IIa CVD Diamond:

1. Solid Start with high-purity Type IIa CVD diamond ([N] < 5 ppb, [B] < 1 ppb). Implant ¹¹B⁺ at 30 keV (range ~55 nm via SRIM) at 10¹³ ions/cm².
   ¹¹B⁺(30 keV) → stops at 55±15 nm depth in diamond
   Creates ~3 V per B ion along track
2. Solid Implant ¹⁴N⁺ at 35 keV (range ~55 nm, matched to B depth) at 10¹³ ions/cm².
   ¹⁴N⁺(35 keV) → stops at 55±18 nm depth
   Overlapping B and N implant profiles maximize B-N nearest-neighbor probability
3. Solid Anneal 700°C 1h forming gas → V captured by B-N pairs. Competing reactions: V+N→NV, V+B→BV (both undesired).
   Desired: V + B_s-N_s(pair) → VBN
   Competing: V + N_s(isolated) → NV (~637 nm, red — impurity signal)
   Competing: V + V → V₂ (neutral, no useful emission)
   Probability of VBN formation increases when B and N concentrations are comparable and co-located — which Route B achieves through matched implant energies.

Route C — CVD with Simultaneous B₂H₆ and N₂ (gas-phase co-doping):

1. Gas CH₄(1%)/H₂ + B₂H₆(0.5 ppm) + N₂(100 ppm). 800°C, 30 Torr. B₂H₆ thermal decomposition provides atomic B.
   B₂H₆(g) →[plasma]→ 2B· + 3H₂(g)
   N₂(g) →[plasma]→ 2N·
   B· + N· + C(diamond surface) → diamond:(B_s,N_s near-neighbor)
   Challenge: B₂H₆ is extremely toxic (TLV 0.1 ppm) — requires fully contained gas handling with double-containment and toxic gas monitors.
2. Solid Irradiate and anneal per Route A steps 3–4.

Environmental conditions matrix:

Verification: PL at 405 nm excitation at 10 K and 300 K; scan 450–600 nm for new ZPL not matching known N3/H3/NV lines. Electron paramagnetic resonance (EPR) to detect B-N-V coupling signature distinct from isolated NV or BV.

4. Pure Violet — Germanium-Vacancy Neutral (GeV⁰)

Target emission: ~400–430 nm violet from GeV in neutral charge state, or from a Ni-N complex emitting in the violet.

Minimum carrier scale: Single GeV: ≥8 nm nanodiamond. Ni-N complex: ≥10 nm.

Thermodynamic feasibility:

• GeV⁻ is well-characterized at 602 nm (ZPL). GeV⁰ is predicted to have a higher-energy transition due to one fewer electron in the defect orbital → estimated ~400–450 nm (blue-violet).
• GeV⁰ charge state: stabilized when Fermi level is below GeV(−/0) transition level at E_V + 1.8 eV (estimated). Requires p-type environment or H-terminated surface.
• Alternative Ni-N route: NE1 center (Ni-N complex) absorbs at ~430 nm; associated emission may appear in violet-blue if oscillator strength is sufficient.

Route A — CVD GeV + Charge Neutralization:

1. Isotopically pure CVD growth: Gas ¹²CH₄(99.99%)/H₂ + GeH₄(0.05%) at 750°C, 40 Torr, microwave 2.45 GHz. ¹²C lattice eliminates ¹³C nuclear spin noise for narrow linewidths.
   ¹²CH₄(g) + H₂(g) →[plasma]→ C(diamond, s) + 2H₂(g)
   GeH₄(g) →[plasma]→ Ge· + 2H₂(g) → Ge_s(in lattice)
   Growth creates ~1 intrinsic V per 100 Ge incorporations (CVD growth defects)
2. Electron irradiation for additional vacancies: Solid 2 MeV e⁻, 5×10¹⁷ e⁻/cm².
   e⁻(2 MeV) + C(lattice) → V + C_i
   V yield: ~5×10¹⁸ cm⁻³
3. GeV formation anneal: Solid in Vacuum 900°C, 2h. V migrates to Ge → forms split-vacancy GeV.
   V(mobile) + Ge_s → GeV (D₃d split-vacancy geometry)
4. Charge neutralization via H₂ plasma: Plasma Pure H₂ plasma, 800 W, 600°C, 10 Torr, 30 min. Hydrogenation of diamond surface creates positive surface charge layer (negative electron affinity with H-termination reverses to positive EA).
   Diamond-O(s) + H·(plasma) → Diamond-H(s) + OH·(g)
   H-termination: surface band bending → hole accumulation layer → p-type
   GeV⁻ + h⁺(surface-induced hole) → GeV⁰
5. Electrostatic gating (alternative charge control): Solid + Liquid Fabricate on-chip gate electrode: deposit 5 nm Ti / 100 nm Al₂O₃ (ALD) gate dielectric on polished diamond surface. Apply +2–5 V gate voltage to deplete electrons → stabilize GeV⁰.
   Al(CH₃)₃(g) + H₂O(g) →[ALD, 200°C]→ Al₂O₃(s) + CH₄(g) (gate dielectric)
   V_gate = +3 V → E_F shifts below GeV(−/0) → GeV⁰ stabilized

Route B — Ni-N Complex for Violet Emission:

1. Solid + Liquid HPHT growth using Ni-Mn catalyst (no Co/Fe) with 0.5 atm N₂. Ni enters lattice at specific sites.
   C(graphite) + Ni(catalyst) + N₂(g) →[5.5 GPa, 1350°C]→ diamond:(Ni,N) complexes
2. Solid Anneal at 1500°C, 2h, vacuum → Ni-N complexes aggregate into NE1 configuration.
   Ni_s + N_s →[1500°C, diffusion]→ Ni-N complex (NE1 type, ~430 nm absorption)
3. Note: Ni-N complexes are weak emitters (Φ_F ≈ 0.01–0.03). Multiple centers per particle needed for usable intensity.

Environmental conditions matrix:

Verification: Low-T PL (10 K) with 375 nm laser excitation → scan 390–460 nm for GeV⁰ ZPL. Photon correlation (g²(0)) to confirm single-emitter character. Compare with GeV⁻ at 602 nm under same excitation to confirm charge-state switching.

5. Bright Amber/Orange — Neutral Tin-Vacancy (SnV⁰)

Target emission: ~580–610 nm from SnV in neutral charge state. SnV⁻ emits at 619 nm; SnV⁰ predicted to blue-shift by ~20–40 nm due to altered orbital filling → ~580–600 nm (amber-orange).

Minimum carrier scale: ≥15 nm. Sn (covalent radius 1.39 Å vs C 0.77 Å) creates substantial lattice strain; host must accommodate ~4% local volume expansion.

Thermodynamic feasibility:

• SnV formation energy in diamond: E_f ≈ 7–8 eV (high due to size mismatch). Compensated by irradiation-supplied vacancies and favorable Sn-V binding energy (~3 eV).
• SnV⁰ vs SnV⁻: charge transition level (−/0) estimated at E_V + 2.0 eV. Requires Fermi level suppression (p-type environment) to depopulate SnV⁻ → SnV⁰.
• Spin-orbit splitting of SnV⁻ ground state: ~850 GHz → largest of Group-IV series. SnV⁰ splitting is unknown experimentally.

Route A — HPHT with Sn Metal Additive:

1. Catalyst preparation: Solid Fe-Co alloy (70:30) + metallic Sn powder (99.99%, −325 mesh) at 2 wt%. Pressed into pellet with graphite at 500 MPa.
   Fe(s) + Co(s) + Sn(s) + C(graphite, s) →[cold press, 500 MPa]→ pellet (s)
2. HPHT growth: Solid + Liquid 5.5 GPa, 1350°C, 48h. Sn dissolves into Fe-Co melt, co-precipitates into diamond substitutionally. [Sn] in diamond ≈ 1–10 ppm.
   Sn(dissolved in Fe-Co melt, l) + C(diamond growth front, s) → Sn_s(in diamond lattice)
   Single replacement: Sn atom substitutes for C during crystallization
3. Acid removal of catalyst + metallic Sn: Liquid Aqua regia 120°C, 48h (dissolves Fe, Co). Then concentrated HCl at 60°C, 12h (dissolves residual Sn).
   Fe(s) + 4HCl + HNO₃ → FeCl₃(aq) + NO(g) + 2H₂O
   Sn(s) + 2HCl(aq) → SnCl₂(aq) + H₂(g)↑ (single replacement)
   Co(s) + 2HCl(aq) → CoCl₂(aq) + H₂(g)↑
4. Electron irradiation + rapid thermal anneal (RTA): Solid 2 MeV e⁻ at 5×10¹⁷/cm². Then RTA: ramp 50°C/s to 1200°C, hold 5 min, quench at 200°C/s in N₂ gas jet.
   e⁻(2 MeV) + C(lattice) → V + C_i
   V(mobile at 1200°C) + Sn_s → SnV (split-vacancy)
   Rapid quench: freezes charge state before thermodynamic equilibrium → preserves SnV⁰
   Why RTA: At equilibrium, SnV⁻ is more stable. Rapid quench traps the kinetic product SnV⁰ before electrons from distant N donors can transfer to SnV.

Route B — FIB Implantation + Anneal:

1. Solid Type IIa CVD substrate. Focused ¹²⁰Sn²⁺ beam at 400 keV, 10¹² ions/cm², spot size 30 nm for localized single-center creation.
   ¹²⁰Sn²⁺(400 keV) → range ~80 nm in diamond (SRIM)
   Co-creates ~50 V per Sn ion (displacement cascade)
2. Solid RTA 1200°C, 5 min, 200°C/s quench. V captured by Sn. Most V recombine with C_i or form V₂.
   V + Sn_s → SnV (at implant depth ~80 nm)
3. Plasma H₂ surface termination to push Fermi level below SnV(−/0) → stabilize SnV⁰.
   Diamond-O(s) + H·(plasma) → Diamond-H(s) → surface p-type → SnV⁰

Route C — Electrolyte Charge Tuning (nanodiamond suspension):

1. Liquid Suspend SnV-containing nanodiamonds (from Route A + milling) in pH 3 buffered electrolyte (citrate buffer). Apply +1.5 V vs Ag/AgCl with Pt working electrode.
   SnV⁻(in ND) →[electrochemical oxidation, +1.5 V]→ SnV⁰(in ND) + e⁻(to electrode)
   pH dependence: lower pH → more positive surface charge → stabilizes SnV⁰. Above pH 7, SnV⁻ dominates.

Crystallization kinetics: HPHT diamond growth rate with Sn additive is reduced ~30% vs undoped (Sn lattice strain creates local growth barriers). Expected: R ≈ 0.7–3.5 mg/h. Sn incorporation is growth-sector dependent: [Sn] on {111} faces is ~3× higher than {100}.

Environmental conditions matrix:

Verification: PL at 532 nm excitation; expect SnV⁻ at 619 nm. Under H₂-terminated or electrochemical bias conditions, monitor for new peak at 580–600 nm (SnV⁰). Temperature-dependent PL from 10 K to 300 K to map SnV⁰ thermal quenching behavior.

6. Full Visible Spectrum — Multi-Layer Nanodiamond

Target emission: 400–700 nm continuous broadband emission from concentric doped/irradiated shells, each hosting a distinct color center.

Minimum carrier scale: ~150–200 nm total particle diameter (core + 4 shells of 20–40 nm each).

Thermodynamic feasibility:

• Each shell is independently optimized for its color center. The key challenge is that each irradiation/anneal cycle must not destroy centers created in previous shells.
• NV⁻ anneal temperature (800°C) is below the onset of N aggregation (~900°C in short timescales), so prior NV centers survive. H3 centers are stable to ~1000°C. N3 centers are stable to >1200°C.
• SiV⁻ is stable to ~1200°C. GeV⁻ to ~1000°C. Ordering: create the most stable centers first (N3 innermost), least stable last (NV⁻ outermost).

Process Chain (layer-by-layer build-up):

1. Core seed — detonation nanodiamond (5 nm): Solid + Gas Detonation synthesis: TNT/RDX (60:40) charge in steel chamber. Detonation at ~3000°C, ~30 GPa, microsecond timescale → carbon condenses as sp³ nanodiamond in cooling wave.
   2C₇H₅N₃O₆(TNT) + C₃H₆N₆O₆(RDX) →[detonation]→
   21C(diamond, s, 4–5 nm) + 6N₂(g) + 8H₂O(g) + 9CO₂(g)
   Note: incomplete combustion yields diamond in carbon-rich core of detonation wave
   Purify: boiling HClO₄/H₂SO₄ (1:3) 72h → removes sp² carbon, metals.
   C(sp², graphite, s) + 2HClO₄(aq) →[250°C]→ CO₂(g) + 2HCl(aq) + O₂(g)
   Fe/Cr/Ni(from chamber walls, s) + HCl(aq) → metal chlorides(aq)
2. Shell 1 — Blue (N3 centers, 415 nm): Gas CVD overgrowth on DND core in fluidized bed or rotating substrate reactor. CH₄(3%)/H₂ + N₂(2000 ppm), 900°C, 50 Torr. Grow 30 nm shell. High N → A-aggregates form during growth.
   ND-core + CH₄/H₂/N₂(high) →[CVD, 900°C]→ ND-core@Shell-1:N(A-agg.)
   Then e⁻ irradiation 10¹⁷/cm² + HPHT anneal (1600°C, 6 GPa, 20 min) → A→B + N3 formation.
   N₂(A-agg.) →[HPHT 1600°C]→ N₃V (N3, 415 nm blue) + N₄V₂ (B-agg.)
3. Shell 2 — Green (H3 centers, 503 nm): Gas CVD overgrowth CH₄(2%)/H₂ + N₂(500 ppm), 850°C. Grow 25 nm. Moderate N → A-aggregates.
   Shell-1 + CVD → Shell-2:N(A-agg., moderate)
   Then proton irradiation 300 keV (range ~2 µm, penetrates entire particle), 10¹⁶/cm² + anneal 600°C 2h.
   V + N₂(A-agg.) →[600°C]→ NVN (H3, 503 nm green)
   N3 centers in Shell 1 are stable at 600°C — no damage.
4. Shell 3 — Orange/Yellow (NV⁰ at 575 nm + SiV⁻ at 738 nm): Gas CH₄(1%)/H₂ + N₂(50 ppm) + SiH₄(30 ppm), 800°C. Grow 30 nm.
   Shell-2 + CH₄/H₂/N₂(low)/SiH₄(trace) →[CVD, 800°C]→ Shell-3:(N_s isolated, Si_s)
   e⁻ irradiation 5×10¹⁷/cm² + anneal 800°C 2h.
   V + N_s(isolated) → NV (575 nm as NV⁰ due to low N concentration)
   V + Si_s → SiV⁻ (738 nm, stabilized by N donors)
   Prior H3 in Shell 2 stable at 800°C. Prior N3 in Shell 1 stable.
5. Shell 4 — Red cap (NV⁻ at 637 nm): Gas CH₄(2%)/H₂ + N₂(500 ppm), 850°C. Grow 25 nm.
   Shell-3 + CH₄/H₂/N₂(moderate) →[CVD]→ Shell-4:N_s(isolated, ~50 ppm)
   e⁻ irradiation 10¹⁸/cm² + anneal 800°C 2h → high NV⁻ density.
   V + N_s → NV ; O₂ plasma terminates surface → NV⁻ stable
6. Final surface treatment: Plasma O₂ plasma 300 W 5 min → stabilizes NV⁻ in outermost shell. Inner shells unaffected (surface treatment penetration <2 nm).
   Surface C-H → Surface C=O (O₂ plasma)
   NV⁰(outer shell) + e⁻(from O-surface band bending) → NV⁻

Environmental conditions matrix:

Verification: Single-particle PL spectroscopy with 405 nm excitation → full 400–750 nm emission spectrum. Measure CIE coordinates per particle. Transmission electron microscopy (TEM) to confirm core-shell morphology. Raman spectroscopy to verify sp³ quality in each shell (1332 cm⁻¹ peak, FWHM < 10 cm⁻¹).

7. Deep UV Fluorescence — Free-Exciton Recombination in Ultra-Pure Diamond

Target emission: ~225–235 nm (5.27 eV) from band-edge free-exciton radiative recombination. This is diamond's intrinsic emission, suppressed at room temperature by phonon-assisted non-radiative decay.

Minimum carrier scale: Bulk single-crystal ≥50 µm thick (exciton mean free path at 10 K). Not achievable in nanodiamonds due to surface quenching of excitons.

Thermodynamic feasibility:

• Diamond indirect band gap: 5.47 eV → direct gap at Γ point: ~7.3 eV. Free-exciton binding energy: ~80 meV → exciton recombination at 5.47 − 0.08 = 5.39 eV (~230 nm).
• At room temperature: exciton-phonon scattering rate exceeds radiative rate by ~10⁴× → Φ_F ≈ 10⁻⁵. At 10 K: phonon population freezes → Φ_F ≈ 10⁻³.
• Any impurity (>1 ppb N or B) creates defect states that trap excitons before radiative recombination → requires extreme purity.

Route A — Ultra-Pure Homoepitaxial CVD:

1. Substrate: Solid Type IIa HPHT seed, (100)-oriented, [N] < 1 ppm, [B] < 50 ppb. Polish both faces to Ra < 1 nm with scaife (cast-iron wheel + diamond paste). Final reactive ion etch (RIE) in Ar/O₂ to remove subsurface damage.
   Diamond(surface damage) + O₂(plasma, RIE) →[Ar⁺ bombardment]→ CO₂(g) (removes ~200 nm)
2. Ultra-high-purity CVD chamber: Gas All-metal sealed chamber, base pressure <10⁻⁹ Torr (turbomolecular + ion pump). Bake at 200°C 48h before growth. Gas purity: ¹²CH₄ (isotopically enriched 99.999%), H₂ (99.99999%, "seven-nines"). N₂ < 0.1 ppb (getter-purified). B < 0.01 ppb.
   ¹²CH₄(99.999%) + H₂(99.99999%) →[microwave plasma, 2.45 GHz, 3 kW]→
   C(diamond, homoepitaxial, s) + 2H₂(g)
   Growth rate: 0.5–1 µm/h on (100) at 900°C, 150 Torr
   [N] in grown layer < 1 ppb, [B] < 0.1 ppb
3. Surface polish: Solid Chemical-mechanical polish (CMP) with colloidal SiO₂ on (100) face → Ra < 0.3 nm. Then O₂ plasma clean.
   Diamond(surface, s) + SiO₂(colloidal, slurry) →[mechanical + chemical]→ atomically flat (100)
4. Cryogenic excitation: Gas + Solid Mount sample in He closed-cycle cryostat, cool to 10 K. Excite with ArF excimer laser (193 nm, 6.42 eV — above band gap) or synchrotron radiation (tunable 190–220 nm). Collect emission through MgF₂ window (transparent to 120 nm).
   hν(193 nm, 6.42 eV) + Diamond →[above-gap excitation]→ e⁻(CB) + h⁺(VB)
   e⁻ + h⁺ →[Coulomb attraction, 10 K]→ free exciton (FE)
   FE → hν(230 nm, 5.39 eV) + phonon(s) (radiative recombination)
   Phonon replicas: 230 nm + n×TA (165 meV) or TO (141 meV)

Route B — Electron Beam Excitation (cathodoluminescence):

1. Solid Same ultra-pure sample in SEM with cryostage (10 K). 10 keV electron beam → generates e-h pairs throughout excitation volume.
   e⁻(10 keV) → ~2800 e-h pairs (E_gap/3 rule) per primary electron
   e-h → FE → hν(230 nm) at 10 K
   Advantage: no UV-transparent window needed. Collection via spectrometer mounted on SEM.

Environmental conditions matrix:

Verification: UV spectrometer (VUV-capable, 180–300 nm range). Expect free-exciton emission at ~230 nm with phonon replicas at ~237, ~244, ~251 nm (TA, TO, LO phonon sidebands). Compare with 10 K vs 77 K vs 300 K → should quench dramatically above ~50 K. Absence of 235 nm peak indicates N or B impurity exceeds threshold.

8. Extended IR Fluorescence (>1000 nm) — Divacancy Chains and Cluster States

Target emission: 1000–1600 nm from extended divacancy chains and cluster defect states for O-band (1260–1360 nm) and C-band (1530–1565 nm) telecom wavelengths.

Minimum carrier scale: ~30 nm (chain of ≥5 linked V₂ divacancies along ⟨110⟩). Bulk single-crystal preferred for highest quality.

Thermodynamic feasibility:

• Isolated divacancy V₂: formation energy ~7 eV from thermal equilibrium, but easily created by irradiation (each 2 MeV e⁻ creates ~10 Frenkel pairs).
• V₂ aggregation: V₂ is mobile above ~700°C (E_a ≈ 3.5 eV). V₂ + V₂ → V₄ (chain extension) is exothermic by ~1 eV per addition.
• Chain electronic states: as chain length increases, in-gap states form mini-bands → emission red-shifts progressively from 741 nm (GR1, single V) → 986 nm (H2, NVN⁻) → >1000 nm (V₂ chains).

Route A — High-Dose Neutron Irradiation + Staged Anneal:

1. Starting material: Solid Type IIa CVD single-crystal, [N] < 5 ppb. Low N ensures vacancies are not captured by N (which would form NV instead of V₂ chains). Polish both faces Ra < 5 nm.
2. Neutron irradiation: Solid Nuclear research reactor, thermal + fast neutron spectrum. Fluence 10¹⁹ n/cm² (requires ~1 month in a 10¹⁴ n/cm²/s flux reactor). Each fast neutron creates a displacement cascade of ~100 Frenkel pairs → [V] ≈ 10²¹ cm⁻³ (heavily damaged).
   n(fast, >0.1 MeV) + C(lattice) →[knock-on cascade]→ ~100 V + ~100 C_i
   n(thermal) + ¹²C → ¹³C (neutron capture, minor) or ¹⁴C (double capture, trace)
   ¹²C(n,γ)¹³C: σ = 3.5 mb — negligible isotopic shift at 10¹⁹ fluence
   Activation products: ¹²C(n,γ)¹³C (stable), no significant radioactivation of pure diamond. However, any metallic inclusions activate → acid-clean before irradiation.
3. Staged thermal anneal — 3-step vacancy aggregation: Solid in Gas (Ar, 1 atm or vacuum)
   Step 1 (400°C, 4h): Interstitials (C_i) mobile at E_a ≈ 1.7 eV. C_i recombines with nearby V → reduces total V count by ~50%. Surviving V are isolated.
   C_i(mobile at 400°C) + V(nearby) → C(lattice restored) — Frenkel recombination
   Step 2 (700°C, 8h): Single vacancies V become mobile (E_a ≈ 2.3 eV). V encounters other V → forms V₂ divacancy (binding energy ~2 eV). V₂ is stable and immobile at 700°C.
   V(mobile) + V(stationary or mobile) → V₂ (divacancy, oriented along ⟨110⟩)
   V + N_s → NV (competing — minimized in Type IIa with [N]<5 ppb)
   Step 3 (1000°C, 4h): V₂ becomes mobile (E_a ≈ 3.5 eV). V₂ encounters V₂ → V₄. V₄ + V₂ → V₆. Chain growth along ⟨110⟩ crystallographic direction (lowest energy configuration).
   V₂(mobile at 1000°C) + V₂ → V₄ (chain along ⟨110⟩)
   V₄ + V₂ → V₆
   V₆ + V₂ → V₈
   General: V_n + V₂ → V_n+2 (ΔE ≈ −1 eV per V₂ addition)
   Quench at 1000°C by turning off heater, cool under Ar → preserves chain configuration.
4. Optional HPHT stabilization: Solid 6 GPa, 800°C, 1h → compressive stress prevents chain dissociation; stabilizes extended defects against thermal fluctuation.
   V_n chain →[6 GPa compression]→ V_n chain (compressed, stable)

Route B — He⁺ Implantation for Localized V₂ Chains:

1. Solid He⁺ at 350 keV, 10¹⁶ ions/cm² → Bragg peak at ~750 nm depth. Creates ~20 V per He ion in end-of-range damage zone.
   He⁺(350 keV) + C(lattice) → V + C_i (along track) + He(interstitial at 750 nm depth)
   [V] at Bragg peak ≈ 10²⁰ cm⁻³ → high enough for chain formation
2. Solid Anneal 700°C 4h → He diffuses out (E_a ≈ 0.3 eV) + V migrates → V₂ formation in damage zone.
   He(interstitial) →[700°C]→ He(g, escapes via surface)
   V + V → V₂ (localized at 750 nm depth)
3. Solid Anneal 1000°C 2h → V₂ aggregates into chains. Confined to narrow damage zone → higher local chain density.

Route C — Carbon Ion Implantation for Self-Interstitial-Free Damage:

1. Solid ¹²C⁺ at 200 keV into Type IIa. Self-ion implantation creates Frenkel pairs but adds no foreign atoms.
   ¹²C⁺(200 keV) + ¹²C(lattice) → V + ¹²C_i + ¹²C(implanted, becomes C_i or C_s)
   No foreign atom residue — purest vacancy source
2. Solid Staged anneal as Route A steps.

Crystallization kinetics of V₂ chain growth: At 1000°C, V₂ diffusivity D_V₂ ≈ 10⁻¹⁴ cm²/s. In 4h, diffusion length L = √(6Dt) ≈ 9 nm. For chain formation, V₂ must encounter another V₂ within this range → requires [V₂] > (1/L³) ≈ 1.4×10¹⁸ cm⁻³. At initial [V] ≈ 10²¹ and 50% recombination → [V] ≈ 5×10²⁰; after V+V→V₂ → [V₂] ≈ 2.5×10²⁰ cm⁻³. This greatly exceeds the threshold → chain formation is efficient.

Environmental conditions matrix:

Verification: Near-IR PL at 77 K using InGaAs detector (spectral range 900–1700 nm). Excite with 785 nm laser (below GR1 at 741 nm, avoids GR1 fluorescence). Expect broad emission 1000–1600 nm from V₂ chain states. High-resolution PL to resolve individual ZPL features of V₃, V₄, V₅ … V_n (predicted spacing ~15–25 nm between successive chain ZPLs). Positron annihilation spectroscopy (PAS) to confirm vacancy cluster size distribution.

P1. Enhanced Blue Phosphorescence — High-Purity Boron-Doped (Type IIb-Analog)

Target: Blue phosphorescence at ~500 nm with τ > 20 s. The diagnostic signature of natural Type IIb (e.g., Hope Diamond).

Minimum scale: Bulk single-crystal ≥200 µm (DAP recombination requires macroscopic path lengths for donor-acceptor separation).

Thermodynamic feasibility:

• Boron acceptor level: E_V + 0.37 eV. Nitrogen donor level: E_C − 1.7 eV. DAP recombination energy: E_gap − E_A − E_D − e²/(4πε₀εr) where r = donor-acceptor separation. For r ~ 5–50 nm, emission peaks at ~2.5 eV (~500 nm).
• Phosphorescence lifetime scales with DAP separation: τ ∝ exp(2r/a_B) where a_B is the Bohr radius of the shallower impurity. Larger r → longer τ → bluer emission. This is why Type IIb phosphorescence is uniquely long-lived.
• Requires [B] > [N] for p-type behavior (net uncompensated boron). If [N] > [B], donors compensate acceptors and no phosphorescence occurs.

Route A — HPHT Growth with B Doping:

1. Nitrogen-free carbon source: Solid Graphite (99.999%) + Fe-Co catalyst (no Ni — Ni introduces competing defects). Add elemental boron (99.9%) at 0.5–5 ppm by weight. Load capsule under N₂-free conditions (glove box, [N₂] < 1 ppm).
   C(graphite, 5N) + B(s, 0.5–5 ppm) + Fe-Co(catalyst) →[assembled in N₂-free glove box]→ capsule
2. HPHT growth: Solid + Liquid 6.0 GPa, 1400°C, 72h. Gettering: add Ti or Zr powder (0.1 wt%) to catalyst as nitrogen getter — Ti forms TiN, preventing N from entering diamond.
   C(graphite) →[dissolves in Fe-Co melt, 6 GPa]→ C(diamond, s)
   B(dissolved in melt) → B_s(substitutional, acceptor, p-type)
   N(trace in melt) + Ti(s) → TiN(s) (removed from melt as stable nitride)
   Ti + ½N₂ → TiN (ΔG = −308 kJ/mol at 1400°C; thermodynamically favorable)
3. Acid clean + surface preparation: Liquid Aqua regia 120°C 48h + H₂SO₄/HClO₄ 250°C 24h. Polish to Ra < 5 nm. O₂ plasma 5 min for clean surface.
   Fe-Co(s) + HCl/HNO₃ → metal chlorides(aq)
   Ti/TiN inclusions: dissolved by HF(aq) + HNO₃(aq) if needed
4. Phosphorescence verification: Solid Excite with SWUV 254 nm (6 W Hg lamp) for 30 s. Remove excitation → photograph at 1 s intervals. Blue afterglow at ~500 nm should persist >20 s. Measure decay with Si photodiode + oscilloscope → fit I(t) = I₀·exp(−t/τ) + background.

Route B — CVD B-Doped Diamond:

1. Gas CH₄(0.5%)/H₂ + B₂H₆(0.1–1 ppm). 800°C, 30 Torr, microwave. Low CH₄ for high purity. No N₂ in gas feed ([N₂] < 0.01 ppm via getter purification).
   CH₄(g) + H₂(g) + B₂H₆(trace) →[plasma, 800°C]→ diamond:B_s(p-type)
   B₂H₆(g) →[plasma]→ 2B· + 3H₂(g)
   [B] in diamond ≈ 0.1–10 ppm depending on B₂H₆ flow
   Safety: B₂H₆ is pyrophoric and toxic (TLV 0.1 ppm). Double-containment gas lines, toxic gas monitoring, automatic shutdown systems required.
2. Solid No irradiation needed — phosphorescence is intrinsic to boron-doped diamond. Verify as Route A step 4.

Environmental conditions matrix:

P2. Orange Long-Persistent Phosphorescence — Deep Trap production in Type IaB

Target: Orange phosphorescence at ~590 nm with τ > 60 s. Replicates the rare natural phenomenon seen in specific Argyle-type stones.

Minimum scale: Bulk crystal ≥100 µm. Deep-trap phosphorescence requires macroscopic trap density for visible intensity.

Thermodynamic feasibility:

• B-aggregate (N₄V₂) creates deep trap state ~0.8 eV above VB (from thermoluminescence glow-curve data). Electrons excited by SWUV are captured in these traps. Thermal detrapping rate: k = ν₀·exp(−E_trap/kT) where ν₀ ~ 10¹² s⁻¹ (phonon attempt frequency). At RT (kT=0.026 eV), k ~ 10¹²·exp(−0.8/0.026) ≈ 10⁻¹ s⁻¹ → τ ≈ 10 s. For τ > 60 s, need E_trap ≈ 0.9 eV → achieved with specific B-aggregate configurations.
• Requires HPHT annealing of Type Ia diamond to drive A→B aggregation. A→B transition requires ~10⁹ years at 1150°C (mantle) but only ~hours at 2000°C under HPHT conditions.

Process Chain:

1. Starting material: Solid Type IaA natural diamond with [N] > 500 ppm (strong A-aggregate). Or synthetic Type Ib HPHT diamond with high isolated N.
2. A→B aggregation by HPHT treatment: Solid 2200°C, 8 GPa, 10h (requires large-volume press or toroidal cell). At these conditions, N diffusion coefficient D_N ≈ 10⁻¹² cm²/s → diffusion length ~6 nm/h → sufficient for A→B aggregate restructuring.
   2(N-N)(A-aggregate) →[2200°C, 8 GPa]→ N₄V₂(B-aggregate) + C_i
   A→B transition: 4N_s + 2C(lattice) → N₄V₂ + 2C_i
   ΔH ≈ +2 eV (endothermic but driven by high T and entropy)
3. Controlled irradiation for trap activation: Solid Low-dose electron irradiation: 2 MeV, 10¹⁶ e⁻/cm² (10× less than for NV creation). Creates dilute vacancies that are captured by B-aggregates, modifying trap depth.
   V + N₄V₂(B-agg.) → N₄V₃ (modified B-aggregate with deeper trap)
   E_trap increases from ~0.8 eV to ~0.9–1.0 eV → longer phosphorescence lifetime
4. Mild anneal to stabilize: Solid in Gas 500°C, 2h, Ar. Recovers interstitial damage without mobilizing vacancies captured by B-aggregates.
   C_i(mobile at 500°C) + V(isolated) → C(lattice restored)
   N₄V₃ is immobile at 500°C → preserved

Environmental conditions matrix:

Verification: SWUV 254 nm excitation for 60 s, then time-resolved photography in dark. Orange afterglow at ~590 nm should persist ≥60 s. Thermoluminescence (TL) measurement: heat at 2°C/s from RT to 400°C → expect glow peak at ~180°C corresponding to E_trap ≈ 0.9 eV.

P3. Bright Green Persistent Phosphorescence (>60 s) — SiV⁰ Triplet or B-N Trap

Target: Green phosphorescence at ~510–540 nm with τ > 60 s. No known natural or synthetic diamond exhibits this.

Minimum scale: ~50 nm (nanodiamond with embedded trap centers).

Thermodynamic feasibility:

• Approach 1 — SiV⁰ triplet: SiV⁰ has S=1 ground state (triplet). If excited state has singlet character, ISC could populate a metastable state emitting at lower energy than the 946 nm ZPL. Predicted: 946 nm fluorescence + possible ~520 nm triplet emission from higher crystal-field state. Entirely theoretical — never measured.
• Approach 2 — B-N donor-acceptor trap in co-doped diamond: B acceptor + N donor + produced vacancy trap → electron captured at deep level (~1.0 eV) → slow detrapping → recombination at ~2.4 eV (520 nm). Lifetime controlled by trap depth production.

Process Chain (B-N trap approach):

1. CVD co-doped growth: Gas CH₄(1%)/H₂ + B₂H₆(0.5 ppm) + N₂(50 ppm), 800°C, 30 Torr. Produces diamond with both B and N at ~1–10 ppm each.
   CH₄ + H₂ + B₂H₆ + N₂ →[plasma]→ diamond:(B_s, N_s) co-doped
2. Vacancy creation: Solid He⁺ at 100 keV, 5×10¹⁵/cm². Creates vacancy-rich zone at ~300 nm depth with B-N pairs.
   He⁺(100 keV) → V + C_i at co-doped depth
3. Low-temperature anneal for deep trap formation: Solid in Gas 500°C, 4h, Ar. Vacancies partially mobile → some captured by B-N pairs → creates B-V-N deep trap complex. Temperature kept below full V mobilization to prevent NV or BV formation (which consume V without creating traps).
   V(partially mobile, 500°C) + B_s-N_s(pair) → B-V-N (deep trap, E_trap ≈ 1.0 eV)
   Competing: V + N_s(isolated) → NV (consumed at 500°C only if V very close)
4. Surface optimization: Plasma O₂ plasma to remove surface damage, stabilize charge states.

Alternative — persistent phosphor coating on nanodiamond:

1. Solid + Liquid Coat fluorescent nanodiamonds (NV⁻, 637 nm) with SrAl₂O₄:Eu²⁺,Dy³⁺ persistent phosphor shell via sol-gel.
   Sr(NO₃)₂ + 2Al(NO₃)₃ + Eu(NO₃)₃(trace) + Dy(NO₃)₃(trace) →[sol-gel, citric acid]→
   SrAl₂O₄:Eu²⁺,Dy³⁺ precursor gel → calcine 1200°C 2h → green persistent phosphor
   ND@SrAl₂O₄:Eu,Dy (core-shell)
   SrAl₂O₄:Eu²⁺,Dy³⁺ emits at 520 nm with τ > 10 hours. As a shell on nanodiamond, provides persistent green glow while the NV⁻ core provides red fluorescence under active excitation → fusion carrier.

Environmental conditions matrix:

Verification: Excite with LWUV 365 nm for 60 s, measure decay at 520 nm in dark. For B-V-N approach: expect τ ~ 10–120 s if trap depth is correct. For phosphor-shell approach: expect τ > 10 h (SrAl₂O₄ is an industrial-grade persistent phosphor). TEM to confirm core-shell morphology.

P4. Violet Persistent Phosphorescence — N9 Triplet Pathway

Target: Violet delayed emission at ~400–420 nm with τ ~ 1–5 s from the N9 center triplet state.

Minimum scale: Bulk crystal ≥50 µm (N9 is a UV absorption center; phosphorescence requires macroscopic path for sufficient absorption).

Process Chain:

1. Starting material: Solid Type IaA or IaAB diamond with high [N] (>500 ppm). N9 center is common in nitrogen-rich diamonds as a UV absorption feature at 236 nm.
2. Deep-UV excitation protocol: Gas Deuterium lamp (160–400 nm continuum) with 230 nm bandpass filter (10 nm FWHM), or KrF excimer laser at 248 nm. Excite for 120 s at 10 mW/cm².
   hν(236 nm, 5.25 eV) + N9(S₀) → N9*(S₁) →[ISC, k_ISC]→ N9*(T₁)
   N9*(T₁) → N9(S₀) + hν(~410 nm, phosphorescence, τ ~ 1–5 s)
   Note: N9 ISC efficiency is unknown — this is the fundamental unknown. If N9 has strong spin-orbit coupling (due to multiple N atoms), ISC may be efficient enough for detectable phosphorescence.
3. Cryogenic enhancement: Solid Cool to 77 K (liquid N₂). Lower temperature suppresses non-radiative decay → enhances phosphorescence quantum yield. Measure T-dependent phosphorescence: 10 K, 77 K, 150 K, 300 K.
   At 77 K: k_nr reduced by ~10×; τ(phos.) expected to increase proportionally
   If τ(300K) ~ 1 s, then τ(77K) ~ 10 s (predicted)

Environmental conditions: Deep-UV excitation requires VUV-transparent optics (MgF₂ lens, CaF₂ windows). Measurement in complete darkness (<10⁻⁸ W/cm² stray light). PMT or EMCCD detector with 400–430 nm bandpass filter for phosphorescence detection. Signal may be extremely weak — photon counting mode may be needed.

Verification: Time-gated PL: excite at 236 nm, gate detector to open 10 ms after excitation ceases, collect 390–440 nm for 100 s. If violet phosphorescence exists, it will appear as a decaying signal distinguishable from instrumental background by its exponential decay form.

P5. Near-IR Persistent Phosphorescence — Vacancy Chain Trap States

Target: Near-IR delayed emission at ~800–1000 nm with τ ~ ms to seconds from vacancy chain defect traps.

Minimum scale: Bulk crystal ≥100 µm (vacancy chains require extended crystal for chain formation).

Process Chain:

1. Create V₂ chains per synthesis pathway 8 (IR fluorescence): Solid High-dose neutron irradiation (10¹⁹ n/cm²) + staged anneal (400→700→1000°C) in Type IIa.
   V₂ chains along ⟨110⟩ (from fluorescence synthesis pathway 8)
2. Additional shallow-trap doping: Solid Low-dose N⁺ implantation (10¹² ions/cm², 200 keV) into the chain-rich region. N creates shallow traps (~1.7 eV below CB) adjacent to V₂ chains.
   N⁺(200 keV) → N_s at ~150 nm depth (overlapping V₂ chain zone)
   N_s(shallow trap) captures electron → slow release to V₂ chain → delayed IR emission
3. Anneal 600°C 1h: Solid in Gas Ar. N becomes substitutional. V₂ chains not mobilized (need >700°C). N traps are now proximal to V₂ emitters.
   N(implant damage) →[600°C]→ N_s(substitutional) adjacent to V₂ chains
   Trap-mediated emission: hν(excitation) → e⁻ to CB → captured by N_s trap → thermally released → recombines at V₂ chain → hν(800–1000 nm, delayed)

Environmental conditions: Neutron irradiation per pathway 8. N⁺ implantation in standard ion implanter. Anneal in tube furnace under Ar flow. IR detection: InGaAs photodiode array (900–1700 nm) + lock-in amplifier for time-resolved phosphorescence at ms timescales.

Verification: Excite with 785 nm laser (above V₂ chain absorption), switch off, measure decay at 900–1100 nm with InGaAs detector. Expected: multi-exponential decay with components at ~10 ms (direct), ~100 ms (N-trap mediated), and possibly ~1 s (deep N-V₂ coupling). Compare sample with V₂ chains alone (no N) to sample with N+V₂ chains — N-doped should show longer-lived component.

FJH-1. Flash Joule Diamond (2 ms) — Nanodiamond to Lonsdaleite-Diamond Composites

Target: Nanodiamond (4–50 nm), lonsdaleite inclusions, and cubic-hexagonal composite particles. Fluorescence: NV⁻ (637–640 nm), stacking fault luminescence (450–550 nm).

Minimum scale: ~5 nm single particles. Maximum scale: Continuous roll-to-roll FJH → kg/h throughput → km-scale film deposition achievable.

Thermodynamic basis:

• FJH peak temperature: 3000–4000 K in ~1–2 ms. At these temperatures, carbon is above the graphite→diamond equilibrium line even at ambient pressure due to kinetic trapping — the cooling rate (~10⁴ K/s) quenches the high-T metastable sp³ phase before it can revert to graphite.
• Lonsdaleite forms preferentially under rapid shock-like conditions because wurtzite stacking (ABAB) has lower nucleation barrier than zinc-blende (ABCABC) at high undercooling. FJH cooling rate favors wurtzite nucleation → lonsdaleite-diamond composite.
• N incorporation from precursor (e.g., melamine, urea) creates NV centers directly during FJH — no post-irradiation needed.

Nanometer-scale pathway (lab, ~mg/pulse):

1. Precursor preparation: Solid Carbon black (Vulcan XC-72, 50 nm primary particles) + 5 wt% melamine (C₃H₆N₆, nitrogen source) + 2 wt% Fe₂O₃ nanoparticles (optional catalyst). Mix by ball-milling 1h. Pack 100 mg into quartz tube (4 mm ID) between copper electrodes.
   Carbon black(s) + melamine(s) + Fe₂O₃(s, optional) →[ball mill]→ packed precursor (ρ ~ 0.5 g/cm³)
2. FJH pulse — diamond formation (2 ms): Plasma Discharge capacitor bank (0.5 F, charged to 150 V) through sample. Peak current ~500 A, energy density ~5–10 kJ/g. Sample reaches ~3000–3500 K in <1 ms.
   C(amorphous/graphitic, s) →[3000–3500 K, 2 ms pulse]→ C(sp³, nanodiamond+lonsdaleite, s) + C(sp², residual graphene shells)
   N(from melamine) → N_s(substitutional in diamond lattice) at ~100–500 ppm
   V(thermal vacancies at 3500 K) + N_s →[cooling at 10⁴ K/s]→ NV centers (formed during quench)
   Key: at 3500 K, thermal vacancy concentration [V] ~ exp(−E_f/kT) ~ exp(−7/0.3) ~ 10⁻¹⁰ per site → insufficient for NV. However, the amorphous→crystalline transition creates copious structural vacancies during recrystallization → [V] ≫ equilibrium → NV forms during rapid crystallization.
3. Rapid quench: Solid Sample cools from 3500 K to RT in ~50 ms (natural radiative + conductive cooling in quartz tube). Cooling rate ~7×10⁴ K/s quenches sp³ phase. Product: mixed-phase powder: ~30–60% nanodiamond (4–20 nm), ~10–30% lonsdaleite, ~20–40% graphene/amorphous carbon shells.
4. Purification: Liquid Boiling HNO₃ (65%, 120°C, 12h) → dissolves Fe₂O₃ + amorphous carbon. Then H₂SO₄/HClO₄ (3:1, 250°C, 6h) → dissolves graphene shells. Centrifuge → collect nanodiamond + lonsdaleite fraction (>95% sp³).
   C(sp², amorphous) + HClO₄(aq, 250°C) → CO₂(g) + HCl(aq)
   Fe₂O₃(s) + 6HCl(aq) → 2FeCl₃(aq) + 3H₂O(l)
   C(sp³, nanodiamond) — resistant to acid oxidation at 250°C
5. Post-FJH NV activation (optional — enhances Φ_F): Solid Mild e⁻ irradiation (2 MeV, 10¹⁷/cm²) + 800°C anneal 2h → creates additional NV centers from N already in lattice. Increases [NV] by ~5–10× over FJH-only.
   e⁻(2 MeV) + C(ND lattice) → V + C_i
   V(mobile at 800°C) + N_s(from melamine, already in lattice) → NV

Large-scale pathway (industrial, kg/h → km-scale film deposition):

1. Continuous FJH reactor: Solid + Plasma Conveyor belt feeds carbon precursor pellets (50 g each) between roller electrodes. Pulse rate: 1 Hz (1 pellet/s). Each pulse: 2 ms, 150 V, 500 A. Throughput: 50 g × 0.5 (diamond yield) × 3600 s/h = ~90 kg/h of raw nanodiamond-lonsdaleite mixture.
   Precursor pellet(50 g) →[2 ms FJH pulse, conveyor]→ ND-Lons. mixture(~25 g) + volatiles(CO₂,N₂, ~25 g)
2. Inline acid wash: Liquid Continuous counter-current acid wash (HNO₃ → H₂SO₄/HClO₄) in PTFE-lined flow reactor. Residence time 4h. Output: purified nanodiamond-lonsdaleite suspension in DI water.
3. Film deposition: Liquid Spray-coat nanodiamond suspension onto substrate (glass, polymer film, Si wafer) via roll-to-roll slot-die coating. 1 µm thick film at 10 m/min line speed → 600 m/h → 14.4 km/day of fluorescent nanodiamond film.
   ND suspension(aq) →[slot-die coating, 10 m/min]→ ND thin film(1 µm, on flexible substrate)
   Post-bake: 150°C 10 min → evaporates water → adherent ND film

Fluorescence by wavelength — FJH diamond products:

Environmental conditions:

FJH-2. Flash Joule Graphene (3 ms) — GQDs, GO, and Fluorescent Carbon Variants

Target: Graphene (few-layer, turbostratic), graphene quantum dots (2–30 nm), and graphene oxide with tunable fluorescence from UV-blue (430 nm) through visible to red (650 nm). Φ_F range: 0.02–0.80 depending on size, doping, and functionalization.

Minimum scale: Individual GQDs ~2 nm. Maximum scale: Continuous FJH → tonnes/day of graphene; GQD extraction by sonication/oxidation of FJH graphene → kg/h GQDs.

Thermodynamic basis:

• FJH at 3 ms pulse reaches ~3000 K. Carbon precursors (waste plastics, biomass, coal, carbon black) undergo rapid graphitization: amorphous C → turbostratic few-layer graphene. Longer pulse (3 ms vs 2 ms for diamond) allows graphene sheet crystallization rather than sp³ nucleation.
• GQDs are produced by controlled fragmentation of FJH graphene sheets: oxidative cutting (HNO₃/H₂SO₄), sonication, or secondary low-energy FJH pulses that fragment sheets into quantum-confined dots.

Nanometer-scale pathway — Fluorescent GQDs from FJH graphene:

1. FJH graphene production (3 ms): Solid + Plasma Carbon black (200 mg) between Cu electrodes. Capacitor: 0.5 F, 120 V. 3 ms discharge. Product: turbostratic graphene (~90% few-layer, <10 layer).
   C(amorphous/carbon black, s) →[3000 K, 3 ms FJH pulse]→ C(turbostratic graphene, s, few-layer)
   Yield: ~90% sp² by Raman (I_2D/I_G > 1)
2. Oxidative cutting to GQDs: Liquid Disperse FJH graphene (50 mg) in conc. H₂SO₄/HNO₃ (3:1, 50 mL). Sonicate 2h at 40°C (ultrasonic bath, 40 kHz). Then reflux 100°C 12h. Cool, dilute 10× with DI water, dialyze (3.5 kDa MWCO) against DI water 48h.
   Graphene sheet(s) + H₂SO₄/HNO₃(aq) →[sonication + reflux]→
   GQD(2–10 nm) + small organic fragments + CO₂(g)
   Edge groups: -COOH, -OH, -C=O (from oxidative cutting)
   GQD size controlled by sonication time: 2h → 5–10 nm (green); 6h → 2–5 nm (blue). More sonication = smaller dots = bluer emission.
3. Size-selective separation: Liquid Centrifugal filtration through 10 kDa and 3 kDa membranes:
   >10 kDa retentate: GQDs 10–30 nm → red emission (600–650 nm)
   3–10 kDa fraction: GQDs 5–10 nm → green emission (510–530 nm)
   <3 kDa permeate: GQDs 2–5 nm → blue emission (430–460 nm)
4. N-doping for enhanced Φ_F: Solid + Gas Hydrothermal treatment: GQDs in DI water + urea (1:5 mass ratio) in Teflon-lined autoclave, 200°C, 10h. N atoms incorporate at edges as pyridinic-N and graphitic-N.
   GQD(-COOH, -OH) + CO(NH₂)₂(urea) →[200°C, 10h, hydrothermal]→
   N-GQD(-NH₂, pyridinic-N, graphitic-N) + CO₂(g) + H₂O(l)
   N/C ratio: ~5–10% → Φ_F increases from 0.10→0.50 (blue) or 0.20→0.60 (green)
5. S,N co-doping for orange emission: Liquid Hydrothermal: GQDs + thiourea (SC(NH₂)₂, 1:3 mass ratio), 180°C, 8h.
   GQD + SC(NH₂)₂ →[180°C, hydrothermal]→ S,N-GQD(thiophenic-S, pyridinic-N) + H₂S(g)
   Emission: 570–600 nm orange; Φ_F ≈ 0.15–0.40

Large-scale pathway — Industrial GQD production:

1. Continuous FJH graphene: Conveyor-fed FJH (as FJH-1) with 3 ms pulse and plastic or biomass precursor → tonnes/day of graphene flake.
   Waste PE/PP plastic(s) →[3 ms FJH, continuous]→ turbostratic graphene(s) + H₂(g) + light hydrocarbons(g)
   Throughput: ~100 kg/h per production line
2. Continuous oxidative cutting: Liquid Flow reactor: graphene suspension pumped through H₂SO₄/HNO₃ at 100°C with inline sonication (20 kHz horn). Residence time 6h. Output: GQD-rich solution.
   Graphene(aq suspension) + acid →[flow reactor, sonication]→ GQDs(aq) at ~5 g/L
   Throughput: ~10 kg/h GQDs per flow line
3. Membrane separation + spray drying: Liquid + Gas Cross-flow filtration for size selection. Spray-dry into powder (inlet 200°C, outlet 80°C). Package as dry fluorescent GQD powder.
   GQD solution(aq) →[spray dry, 200°C inlet]→ GQD powder(s) + H₂O(g)
   Output: color-sorted GQD powders (blue, green, yellow, red, orange)
4. km-scale fluorescent film: Redissolve GQD powder in water or ethanol. Slot-die coat onto PET film at 20 m/min. Dry at 80°C inline. → 28.8 km/day of fluorescent GQD film.
   GQD(aq or EtOH solution) →[slot-die, 20 m/min, 80°C dry]→ fluorescent GQD film on PET

Fluorescence by wavelength — GQD/GO products: