C–F Exchange Interactions in Quantum Materials

• C–F exchange interactions are quantum processes where carbon and fluorine atoms exchange electrons via direct and indirect mechanisms, impacting magnetic and reactivity properties.
• In functionalized graphene derivatives such as C₂F and C₂H, these interactions induce geometric frustration, stabilize quantum spin liquids, and support skyrmion formation.
• In CaF+Ca systems, barrierless atom-exchange and isotope-dependent effects enable tunable ultracold reactions and provide precise benchmarks for quantum chemical dynamics.

C-f exchange interactions refer to the quantum mechanical exchange processes involving carbon (C) and fluorine (F) atoms and ions. These interactions govern the magnetic coupling in functionalized graphene derivatives and the intermolecular atom-exchange pathways in CaF+Ca systems. Exchange interactions can be characterized by the competition between direct and indirect mechanisms, by their impact on geometric frustration and magnetism, and by their role in quantum dynamics such as skyrmion formation or ultracold chemical reactivity.

1. Fundamental Principles of C–F Exchange Interactions

At the quantum many-body level, exchange interactions between C and F orbitals arise from two primary mechanisms:

• Direct Exchange: This ferromagnetic coupling originates from the Coulomb integral overlapping neighboring atomic-like Wannier orbitals, most notably the $p_z$ functions on carbon sites in C$_2$F. The magnitude of the direct exchange $J^{\rm dir}_{ij}$ is calculated as

$$J^F_{ij} = \int d\mathbf r\,d\mathbf r'\; \frac{W_i^*(\mathbf r)\,W_j^*(\mathbf r')\,W_i(\mathbf r')\,W_j(\mathbf r)} {|\mathbf r-\mathbf r'|},$$

where $W_i$ denotes the Wannier function for site $i$.

• Indirect (Superexchange) Interaction: This antiferromagnetic contribution is mediated by virtual electron hops involving intermediate orbitals, for example through fluorine $2p$ states in graphene derivatives or metal–ligand bridges in transition metal compounds.

The competition between these mechanisms determines the net isotropic exchange coupling $J_{ij}$, which appears in effective Heisenberg Hamiltonians for the low-energy degrees of freedom:

$$\hat H_{\rm Heis} = -\sum_{i \neq j} J_{ij} \mathbf s_i \cdot \mathbf s_j.$$

2. C–F Exchange in Functionalized Graphene: C$_2$F and C$_2$0H

In single-side fluorinated graphene (C$_2$1F), the magnetism can be mapped onto a triangular lattice of unsaturated carbon $_2$2 orbitals. The key exchange channels and their numerical values (Rudenko et al., 2013, Mazurenko et al., 2016) are:

Coupling Type	Magnitude (meV)	Nature
$_2$3–$_2$4 (C–C, C$_2$5F)	$_2$6	Nearest-neighbor AFM (triangular frustration)
$_2$7–$_2$8 (C–F hybrid)	$_2$9	Ferromagnetic, negligible beyond first shell
$J^{\rm dir}_{ij}$0–$J^{\rm dir}_{ij}$1 (C$J^{\rm dir}_{ij}$2H)	$J^{\rm dir}_{ij}$3–40 (FM)	Strong leading coupling in semihydrogenated graphene
$J^{\rm dir}_{ij}$4–$J^{\rm dir}_{ij}$5 (C–C, C$J^{\rm dir}_{ij}$6H)	Weak AFM	Subdominant to FM $J^{\rm dir}_{ij}$7–$J^{\rm dir}_{ij}$8

Numerical evaluation of the exchange integrals in C$J^{\rm dir}_{ij}$9F (DFT with magnetic force theorem, $$J^F_{ij} = \int d\mathbf r\,d\mathbf r'\; \frac{W_i^*(\mathbf r)\,W_j^*(\mathbf r')\,W_i(\mathbf r')\,W_j(\mathbf r)} {|\mathbf r-\mathbf r'|},$$0 k-points, $$J^F_{ij} = \int d\mathbf r\,d\mathbf r'\; \frac{W_i^*(\mathbf r)\,W_j^*(\mathbf r')\,W_i(\mathbf r')\,W_j(\mathbf r)} {|\mathbf r-\mathbf r'|},$$1 eV) reveals rapid decay of AFM coupling beyond nearest neighbor, excluding long-range RKKY oscillations.

The triangular geometry of unsaturated C sites in C$$J^F_{ij} = \int d\mathbf r\,d\mathbf r'\; \frac{W_i^*(\mathbf r)\,W_j^*(\mathbf r')\,W_i(\mathbf r')\,W_j(\mathbf r)} {|\mathbf r-\mathbf r'|},$$2F, combined with the dominance of $$J^F_{ij} = \int d\mathbf r\,d\mathbf r'\; \frac{W_i^*(\mathbf r)\,W_j^*(\mathbf r')\,W_i(\mathbf r')\,W_j(\mathbf r)} {|\mathbf r-\mathbf r'|},$$3, leads to geometrical frustration. In the classical limit, the system prefers a $$J^F_{ij} = \int d\mathbf r\,d\mathbf r'\; \frac{W_i^*(\mathbf r)\,W_j^*(\mathbf r')\,W_i(\mathbf r')\,W_j(\mathbf r)} {|\mathbf r-\mathbf r'|},$$4 Néel arrangement of spins. Quantum mechanically, moderate itinerancy ($$J^F_{ij} = \int d\mathbf r\,d\mathbf r'\; \frac{W_i^*(\mathbf r)\,W_j^*(\mathbf r')\,W_i(\mathbf r')\,W_j(\mathbf r)} {|\mathbf r-\mathbf r'|},$$5), a small effective Hubbard $$J^F_{ij} = \int d\mathbf r\,d\mathbf r'\; \frac{W_i^*(\mathbf r)\,W_j^*(\mathbf r')\,W_i(\mathbf r')\,W_j(\mathbf r)} {|\mathbf r-\mathbf r'|},$$6, and a $$J^F_{ij} = \int d\mathbf r\,d\mathbf r'\; \frac{W_i^*(\mathbf r)\,W_j^*(\mathbf r')\,W_i(\mathbf r')\,W_j(\mathbf r)} {|\mathbf r-\mathbf r'|},$$7 local moment regime allow for quantum-spin-liquid ground states.

In C$$J^F_{ij} = \int d\mathbf r\,d\mathbf r'\; \frac{W_i^*(\mathbf r)\,W_j^*(\mathbf r')\,W_i(\mathbf r')\,W_j(\mathbf r)} {|\mathbf r-\mathbf r'|},$$8H, the strong direct $$J^F_{ij} = \int d\mathbf r\,d\mathbf r'\; \frac{W_i^*(\mathbf r)\,W_j^*(\mathbf r')\,W_i(\mathbf r')\,W_j(\mathbf r)} {|\mathbf r-\mathbf r'|},$$9–$W_i$0 ferromagnetic exchange suppresses long-range magnetic order at finite temperature, consistent with the Mermin-Wagner theorem.

3. Mechanisms of Direct and Indirect C–F Exchange

Direct exchange in C–F-related systems is quantified via constrained random phase approximation (cRPA) and real-space Coulomb integrals. The competition with kinetic (superexchange) terms is critical. For C$W_i$1F (Mazurenko et al., 2016):

Quantity	Value (meV)
$W_i$2 (AFM)	40
$W_i$3 (screened FM)	18
$W_i$4 (bare FM)	44
$W_i$5 (net, screened)	22
$W_i$6 (net, bare)	$W_i$7
$W_i$8 (DM vector)	$W_i$90.98

The net exchange $i$0 is antiferromagnetic unless the screened direct ferromagnetic term exceeds $i$1. Screening by environment or strain can tune $i$2 through zero, placing the system at the AFM–FM instability threshold.

Dzyaloshinskii–Moriya interaction (DMI) arises from spin–orbit coupling and superexchange processes. For C$i$3F, the DMI reaches magnitudes $i$41 meV, rendering $i$5 of order unity and enabling stabilization of Néel-type skyrmion lattices under moderate fields ($i$6 T, $i$7 K when $i$8 meV).

4. Atom–Exchange Pathways in CaF+Ca Systems

Ab-initio quantum chemistry methods (CCSD(T), MRCI) reveal the structure of ground and excited-state potential energy surfaces (PES) for CaF+Ca, which determine the C–F atom-exchange dynamics (Sardar et al., 27 Oct 2025). The nuclear Hamiltonian is

$i$9

with $2p$0 parametrized in Legendre polynomials. The ground-state X $2p$1 surface is deeply bound ($2p$2 cm$2p$3, bent geometry $2p$4), with a strong angular anisotropy from $2p$5. Excited (2) $2p$6 surface associated with CaF ($2p$7)+Ca ($2p$8) is even deeper ($2p$9 cm$J_{ij}$0 along linear $J_{ij}$1).

Atom–exchange is barrierless for the ground channel:

• No transition-state barrier is found in 2D $J_{ij}$2 scans.
• The process is exothermic or isoenergetic (isotope-dependent zero-point energy differences).
• Long-range van der Waals coefficients for $J_{ij}$3CaF+$J_{ij}$4Ca: $J_{ij}$5, $J_{ij}$6.

Excited-state atom-exchange is governed by deep PES wells and strong Ca($J_{ij}$7) spin–orbit coupling, supporting possible nonadiabatic transitions.

5. Isotope-Dependent Exchange and Ultracold Reaction Dynamics

Isotope-exchange reactions of the form

$J_{ij}$8

are controlled entirely by the zero‐point energy difference. The $J_{ij}$9-values, ranging 1–8 cm$$\hat H_{\rm Heis} = -\sum_{i \neq j} J_{ij} \mathbf s_i \cdot \mathbf s_j.$$0, are well below vibrational spacings ($$\hat H_{\rm Heis} = -\sum_{i \neq j} J_{ij} \mathbf s_i \cdot \mathbf s_j.$$1 cm$$\hat H_{\rm Heis} = -\sum_{i \neq j} J_{ij} \mathbf s_i \cdot \mathbf s_j.$$2) and above rotational spacings ($$\hat H_{\rm Heis} = -\sum_{i \neq j} J_{ij} \mathbf s_i \cdot \mathbf s_j.$$3 cm$$\hat H_{\rm Heis} = -\sum_{i \neq j} J_{ij} \mathbf s_i \cdot \mathbf s_j.$$4), ensuring product CaF in $$\hat H_{\rm Heis} = -\sum_{i \neq j} J_{ij} \mathbf s_i \cdot \mathbf s_j.$$5 across multiple $$\hat H_{\rm Heis} = -\sum_{i \neq j} J_{ij} \mathbf s_i \cdot \mathbf s_j.$$6 levels. Reaction is exothermic if $$\hat H_{\rm Heis} = -\sum_{i \neq j} J_{ij} \mathbf s_i \cdot \mathbf s_j.$$7, fully tunable via isotope selection. The absence of an activation barrier implies near-unit reaction probabilities at ultracold collision energies.

	$$\hat H_{\rm Heis} = -\sum_{i \neq j} J_{ij} \mathbf s_i \cdot \mathbf s_j.$$8CaF	$$\hat H_{\rm Heis} = -\sum_{i \neq j} J_{ij} \mathbf s_i \cdot \mathbf s_j.$$9CaF	$_2$0CaF	$_2$1CaF	$_2$2CaF	$_2$3CaF
$_2$4Ca	0	+2.247	+3.301	+4.305	+6.199	+7.947
$_2$5Ca	−2.247	0	+1.053	+2.059	+3.953	+5.701
...

6. Charge-Exchange Excitations and Continuum Effects

Charge-exchange excitations involving C–F pairs in nuclear and solid-state systems are modeled by self-consistent continuum RPA with finite-range Gogny-like interactions (Donno et al., 2016). The inclusion of both direct (Fock) and tensor–isospin channels ensures:

• Proper treatment of Fermi, Gamow–Teller, and spin–dipole operators.
• Preservation of energy-weighted sum rules.
• Strong dependence of SD($_2$6) modes on tensor contributions, evidenced by centroid shifts of 2–9 MeV.

Continuum coupling yields smooth, physical strength distributions above threshold, in contrast to discretized RPA approaches. Finite-range exchange (Fock) explicitly shapes the charge-exchange spectra, particularly for spin-dependent modes.

7. Implications and Tuning of C–F Exchange in Contemporary Research

C–F exchange interactions underpin several modern quantum phenomena:

• Geometric Frustration: Dominant AFM exchange on triangular lattices leads to nontrivial spin textures and suppresses conventional magnetic order.
• Quantum Spin Liquids: Small $_2$7 ratios and moderate moments in C$_2$8F shift ground states from classical order to highly entangled nonmagnetic quantum phases.
• Skyrmion Formation: The interplay of near-cancelled AFM/FM couplings and strong DMI in C$_2$9F enables stabilization of skyrmion crystals under experimentally accessible fields and temperatures.
• Ultracold Chemistry: Deep anisotropic PES and barrierless exchange in CaF+Ca systems allow for tunable synthesis, isotope-selective reactions, and precision benchmarking of quantum chemical dynamics.
• Spectroscopic and Dynamical Probes: The detailed mapping of exchange parameters, decay rates, and mode sensitivities informs experimental design in scanning probe, neutron scattering, and cold-molecule research.

A plausible implication is that environmental tuning (screening via substrates or strain, isotope engineering) can steer C–F exchange interactions across quantum phase boundaries, thereby enabling controlled exploration of frustrated magnetism and exotic reactivity regimes.