Leader-Directed Differential Evolution for Adapter Training
Implementation & benchmark of arXiv:2607.29228 — Linear Proposal Operators and Stochastic Search Geometry in SOMA and Differential Evolution
Source paper: Linear Proposal Operators and Stochastic Search Geometry in SOMA and Differential Evolution arXiv:2607.29228 ⓘ · analyzed Aug 31, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Use differential evolution over adapter or prompt parameters, combining attraction to the current best parameter vector with a population-difference direction. Binomial crossover supplies coordinate-level exploration, while the operator-selection separation makes it possible to measure raw proposal geometry independently from parameter repair and noisy fitness selection.
Formulas
Mathematical statement
For target vector x_i in R^d, best vector x_best, two distinct population vectors x_r1 and x_r2, scale factor F, and diagonal crossover mask C = diag(C_1,...,C_d), the paper gives the DE/current-to-best/1/bin proposal u_i = x_i + C[F(x_best-x_i) + F(x_r1-x_r2)]. Each C_j is Bernoulli with crossover probability CR, except one forced coordinate is set to one so the offspring cannot equal x_i solely because of crossover. The first difference is leader attraction; the second is a population-difference direction. The paper's variation-repair-selection factorization requires analyzing raw u_i before feasibility repair because repair can change mean, covariance, support, and coordinate dependence. In neural adapter training, x_i is the flattened LoRA, prefix, or prompt vector and f is validation loss or a noisy minibatch loss. Optional per-layer normalization should be applied to population differences so large adapter tensors do not dominate the proposal.
Implementation notes
Integrate this at the optimizer level for a frozen pretrained network with trainable LoRA, prefix-tuning, or soft-prompt parameters. Flatten each candidate's trainable parameters into x_i, while retaining a layer map so the differential vector can optionally be normalized separately per tensor. For layer l, compute d_l = (x_r1,l - x_r2,l)/(RMS(x_r1,l-x_r2,l)+epsilon), and use the corresponding normalized leader difference if per-layer scale control is enabled. At every generation, identify x_best using a fixed validation minibatch or an averaged score over two minibatches. For each target, sample distinct r1 and r2, sample C_j from Bernoulli(CR), force one random coordinate or tensor block to be active, calculate u_i = x_i + C elementwise-multiplied by [F(x_best-x_i) + F(x_r1-x_r2)], project into LoRA norm bounds, evaluate the candidate, and retain elites plus the best offspring. Use population size 16 to 32, F in [0.3,0.8], and CR in [0.5,0.9] for the first sweep. Record raw proposal norms, changed-coordinate fraction, coordinate coverage, repaired-coordinate fraction, and loss improvement per forward evaluation. The first experiment should use a small ViT or 100M-parameter transformer with rank-4 or rank-8 adapters on CIFAR-100 or SST-2, comparing against random search, isotropic evolutionary strategies, and gradient-based LoRA under an equal forward-pass budget. The predicted signal is faster validation improvement than isotropic perturbations when population differences capture useful adapter directions; a negative result is equal or worse loss descent after accounting for population evaluations.
Verification
This idea has not been verified yet.
Verification happens in two stages: Stage 1 — a mechanism check on a toy system confirms the claimed mathematical phenomenon reproduces; Stage 2 — a benchmark implements the idea on a real (small) neural network task and compares it against a tuned baseline over 8 paired seeds with a permutation test.
Artifacts
Artifacts unavailable.