Cyclotron Info

Design Guide › Modeling

Modeling design rules

49 of the guide’s 1374 rules carry the modeling tag. Rules for computing before cutting: field maps, orbit codes and their initial conditions, scale-model limits, and the resonances that should be crossed on paper before they are crossed in beam. Each rule keeps its formula where the source gives one, a verbatim quote, a page-level citation, and a stable identifier (dg-NNNN) that resolves here and on the all-in-one guide.

To combine this tag with another (rules carrying both), use the filterable view: /design-guide/?domain=modeling and add a second chip. Related domains, by how often they share a rule with this one: Beam dynamics (16), Magnet (16), RF (13), Dee (5), Beam measurement (2).

Verify before use. Every rule here is a source extract in the vocabulary of the editorial methodology — faithful to its cited page, not an independently validated engineering requirement. Re-read any rule that drives a real design decision at the cited page before committing metal, money, or high voltage to it.

  1. For orbit-code initial conditions, model ions leaving a slit chimney from a flat plasma boundary at ~35,000 K plasma temperature (4.5 eV central starting energy); hole chimneys need a concave boundary. This recipe reproduced measured emittance well enough "that construction of actual cyclotrons can proceed".

    T_plasma ~ 35,000 K -> E_start ~ 4.5 eV; flat boundary (slit), concave (hole)

    ion-sourcebeam-dynamicsmodeling dg-625

    Source, quote & tabletop applicability
    the plasma temperature that provides the best match for experimental beams is approximately 35,000 K (resulting in a central starting energy of 4.5 eV)

    Forringer, Phase Space Characterization of an Internal Ion Source for Cyclotrons — MSU dissertation (2004) — p. 3

    Tabletop: Drop-in starting condition for the reference machine's central-region orbit models — start protons at 4.5 eV from a flat sheet across the slit, not from rest at a point.

  2. Design a regenerator by the seven-step procedure: nonlinear equations of motion -> amplitude- dependent radial/axial frequencies from the measured field (Krylov-Bogoliubov) -> matrix analysis of radial amplitude growth and required momentum kick -> field perturbation from the kick -> resulting axial kick -> axial growth -> pick position from the family of solutions. Gain per turn follows a = -sin(wr*th1)/sin(wr*(th2-th1)); a smaller crossing interval means more gain but a stronger regenerator.

    a = -sin(wr*th1)/sin(wr*(th2-th1)), wr = wr(r>R); delta(p') = -p0''*sin(wr*th2)/ sin(wr*(th2-th1))

    extractionbeam-dynamicsmodeling dg-774

    Source, quote & tabletop applicability
    The determination of the required perturbation for extracting the beam of a synchrocyclotron is made in seven steps.

    Stubbins, Extraction of Synchrocyclotron Beams Near the Maximum Energy — UCRL-3476 (1956) — p. 6

    Tabletop: The workflow (measured field -> amplitude-dependent tunes -> impulse-matrix tracking -> element strength) is exactly the CYCLOPS-lite pipeline planned for a next machine; the specific peeler-regenerator field shapes are synchrocyclotron machinery and need not transfer.

  3. Distrust scale-model magnet studies at excitation extremes: the 0.275-scale ORIC model set the Davis central-region design but could not be operated at the lowest planned field (3.5 kG); full-scale mapping then revealed a defocusing radial-profile depression at low fields that the model never showed, forcing an iron redesign.

    magnetmodeling dg-811

    Source, quote & tabletop applicability
    it was not possible to operate the model magnet at the extremely low (3.5 kilogauss) field levels at which we might like to operate the full scale machine.

    Jungerman, Kibbe & Peek, Central Region Studies for Incorporating an Axial Ion Source in the Davis 76-in. Cyclotron — UCD-CNL-49 (1966) — p. 5

    Tabletop: Modern translation for the FEMM pipeline: a model (physical or FEM) validated at one excitation does not certify another — iron saturation state changes the profile shape, so re-run the field solution at every planned operating point, especially the lowest, and verify the real magnet across its full excitation range.

  4. Get central-region starting conditions by backward tracking: Davis estimated ion starting conditions by placing ions on a known-good 12-in equilibrium orbit and de-accelerating them to the center, then used those conditions to launch forward acceleration runs — bypassing the ill-defined source-gap region on the first pass.

    integrate equations of motion with reversed energy gain from EO inward to r=0

    beam-dynamicsmodeling dg-814

    Source, quote & tabletop applicability
    The starting conditions for all cases were estimated by starting the ions on an equilibrium orbit of 12 inch radius and de-accelerating them to the center.

    Jungerman, Kibbe & Peek, Central Region Studies for Incorporating an Axial Ion Source in the Davis 76-in. Cyclotron — UCD-CNL-49 (1966) — p. 6

    Tabletop: Directly implementable in the Python orbit tracker — find the equilibrium orbit at modest radius (easy, well-conditioned), integrate backwards to the source region, and read off where the source slit and puller must be; cheaper and more robust than guessing forward launch conditions in the messy first gap.

  5. Cross a betatron resonance on paper before crossing it in beam: the Davis orbit code showed particles pass the 3/3 radial resonance at 6-7 in radius, and that the radial oscillation build-up "is not excessive and soon damps to 0.3 inch" — the resonance was accepted, quantitatively, rather than avoided.

    compute amplitude growth through resonance; accept if bounded and damping (here to 0.3 in)

    beam-dynamicsmodeling dg-816

    Source, quote & tabletop applicability
    The particles pass through the 3/3 resonance at a radius of 6-7 inches. The computer calculations show that the radial oscillation build-up at resonance is not excessive and soon damps to 0.3 inch.

    Jungerman, Kibbe & Peek, Central Region Studies for Incorporating an Axial Ion Source in the Davis 76-in. Cyclotron — UCD-CNL-49 (1966) — p. 6

    Tabletop: Method for the CYCLOPS-lite tracker: don't just plot nu_r(r) and forbid resonance lines — integrate through them and report amplitude growth in millimeters against the aperture; a fast-crossed resonance with bounded growth is a non-event even on a small machine.

  6. A computed field map validated by orbit code can produce first beam without empirical shimming iteration: Davis obtained a 21-MeV H2+ internal beam on the first attempt using the computed field, taken as confirmation of both the magnetic measurements and the orbit calculations.

    modelingmagnetbeam-dynamics dg-817

    Source, quote & tabletop applicability
    the validity of the calculations and magnetic field data is supported by the fact that we obtained an internal beam of 21 MeV H2+ ions using the computed field on the first attempt.

    Jungerman, Kibbe & Peek, Central Region Studies for Incorporating an Axial Ion Source in the Davis 76-in. Cyclotron — UCD-CNL-49 (1966) — p. 7

    Tabletop: The 1966 proof that a next machine's compute-first pipeline (field map -> tracker -> build) is sound — if the field is measured carefully and the tracker is honest, first beam on the first pump-down is a reasonable expectation, not luck. Also a period example of commissioning on H2+ rather than protons for shielding reasons.

  7. Develop cyclotron rf on an electrical model before metal is cut: the variable-energy oscillator study used an 8-ft section of the 63-inch dee-stem electrical model with capacitors simulating the dees, and selected a self-excited push-pull oscillator direct on the stems (two tuning controls, drive-insensitive frequency) from competing circuits tested on that model.

    rfmodeling dg-948

    Source, quote & tabletop applicability
    The push-pull circuit for this test was constructed by using an 8-ft section of the electrical model of the 63-in. cyclotron dee stems as the resonant system.

    Howard (ed.), Electronuclear Research Division Semiannual, period ending 20 September 1953 — ORNL-1663 (1954) — p. 19

    Tabletop: A bench-scale dee-stem mockup (pipe sections plus padding capacitors) lets the next machine's oscillator/coupling scheme be raced against alternatives for pocket change - the same measure-on-model philosophy as the deferred filament-coupling rule, one report earlier in hardware form.

  8. Build the model magnet for measurement access: the quarter-scale 114-inch model put the gap in a VERTICAL plane "to provide the greatest access for making field measurements" and made the pole tips removable "so that shims of any shape can be inserted readily" — a model magnet's geometry should serve the probe and the shim swap, not mimic the final machine's orientation.

    magnetmodelingfabrication dg-956

    Source, quote & tabletop applicability
    The pole tips are removable so that shims of any shape can be inserted readily.

    Howard (ed.), Electronuclear Research Division Semiannual, period ending 20 September 1954 — ORNL-1795 (1954) — p. 19

    Tabletop: For any next-machine shim-test rig (or a scaled FEMM-validation magnet), design for the measurement campaign - open sightlines for the Hall probe, pole tips that unbolt, gap oriented for jig access. ORNL judged these features worth 14.4 tons of model.

  9. Retire beam-dynamics risk with an electron-model machine: before committing to a 1-BeV proton AVF cyclotron, ORNL planned "an electron-model accelerator to be used in assessing the importance of imperfection resonances and the feasibility of their penetration" — electrons let you walk the same tune diagram at bench field, energy, and cost (1-BeV accelerator study).

    modelingbeam-dynamics dg-963

    Source, quote & tabletop applicability
    Plans are being made for an electron-model accelerator to be used in assessing the importance of imperfection resonances and the feasibility of their penetration.

    Howard (ed.), Electronuclear Research Division Semiannual, period ending 20 March 1955 — ORNL-1884 (1955) — p. 35

    Tabletop: Historical validation of an electron-model approach and of risk-ordered development - when the open question is orbit dynamics (resonance crossing, field tolerance), a tabletop electron machine answers it before proton iron is bought. (Same page — Thomas-1938 AVF theory, ORACLE computing, MURA spiral sectors - sector focusing arriving in this collection's timeline.)

  10. Retire RF-system risk with a scaled electrical model before cutting full-size metal: build the complete RF circuit at reduced scale (frequency scales inversely with size), verify tuning range, voltage distribution, and power on the bench, then commit to full-scale construction on the model dimensions. The 184-inch followed a three-stage chain: calculation (MacKenzie BP-140), half-scale model (this report), full-size bench test before installation.

    half-scale model resonates at 2x full-scale frequency; geometric ratios and line impedances are scale-invariant

    rfmodelingcyclotron-general dg-999

    Source, quote & tabletop applicability
    Performance of the model is considered sufficiently satisfactory to proceed with the full scale design and construction based on the model dimensions.

    Anderson, Half-Scale Model Tests on the Three Quarter Wave R.F. System — UCRL-31 (1947) — p. 16

    Tabletop: A next machine's dee/stem/tank is already benchtop-sized, so the transferable form is the mockup itself — a cheap RF-only copy (no vacuum) of the dee-liner geometry, swept with a VNA before the vacuum parts are machined. Same lineage as UCRL-64 and MDDC-1045 already in this collection.

  11. Scope the model to the physics it must answer: only the RF circuit was reproduced; the vacuum system, purely mechanical equipment, and the dee-bias insulation were omitted because they have no RF function. Known infidelities were listed, not ignored.

    rfmodeling dg-1000

    Source, quote & tabletop applicability
    Only the radio frequency circuit was simulated in the model, the vacuum system and purely mechanical equipment was not included.

    Anderson, Half-Scale Model Tests on the Three Quarter Wave R.F. System — UCRL-31 (1947) — p. 7

    Tabletop: License to mock up the next machine's RF cavity in bare copper/aluminum on a bench plate — no chamber, no pumps — provided every conducting surface that carries RF current (liner included) is reproduced.

  12. Extrapolate model power to full scale as P proportional to V^2 with a sqrt(2) shunt-impedance credit for the half-scale model (skin depth: doubled size at halved frequency raises Q and R_sh by sqrt(2)). The printed numbers obey it exactly: 520 W at 1.5 kV on the model becomes 146 kW at 30 kV full scale (x400/sqrt(2)); oscillator efficiency held at 59-64% across the band.

    P_full = P_model * (V_full/V_model)^2 * sqrt(s), s = model/full linear scale (=1/2 here, so divide by sqrt(2)); R_sh scales as s^(-1/2) at scaled frequency [scaling law implied by the printed 520 W -> 146 kW pair; verified against all four frequencies]

    rfmodeling dg-1001

    Source, quote & tabletop applicability
    For 30 kv on the full scale system the above power input figures become 146 kw at 9.25 mc, 132 kw at 12 mc, 146 kw at 15 mc, and 98 kw at 23 mc for continuous operation.

    Anderson, Half-Scale Model Tests on the Three Quarter Wave R.F. System — UCRL-31 (1947) — p. 15

    Tabletop: The V^2 term is the live part for a next machine's power budgeting — measured drive power at a safe low dee voltage extrapolates as (V_target/V_test)^2 on the same hardware, since Q is voltage-independent until multipactor/breakdown; a 1500-V measurement predicts the 5-13 kV LDMOS requirement.

  13. A scale model's known infidelities must be listed with the results: substitute 304-TL triodes have much larger internal inductance than the final 9C21s, the filament line's impedance changes where it enters the vacuum system, and mismatched plate capacities between the two tubes skewed early power measurements.

    rfmodeling dg-1011

    Source, quote & tabletop applicability
    the inductance inherent in the 304-TL triodes is large compared with that in the 9C21 triodes to be used in the final oscillator.

    Anderson, Half-Scale Model Tests on the Three Quarter Wave R.F. System — UCRL-31 (1947) — p. 12

    Tabletop: When bench-testing a next machine's RF with a stand-in amplifier or without the real chamber wall, write the fidelity caveats into the test log — the model predicts the cavity, not the parts that were substituted.

  14. Model cyclotron acceleration as kick-plus-coast: an impulsive energy change at each gap azimuth followed by a static (coasting) trajectory to the next gap. Accelerated behavior is then inferred from static phase plots at a few energies, interpolated — validated throughout this study against fully accelerated runs.

    per gap crossing dE = qV_gap; r and p_r unchanged at the kick (radial gaps); coast on the static map between gaps

    beam-dynamicsmodeling dg-1012

    Source, quote & tabletop applicability
    the acceleration can be considered to good approximation as being a simple impulsive change in the energy of the particle at the azimuth of the accelerating gap

    Blosser & Gordon, Computational Study of a Resonant Extraction System for a 3-Sector Cyclotron — MSUCP-9 (1961) — p. 4

    Tabletop: The core architecture for CYCLOPS-lite — thin-gap kicks alternating with magnetic coasting maps is 1961-validated practice, and MSUCP-12's analytic gap field upgrades the kick from a delta function to a distributed one when transit time matters.

  15. Build the orbit toolchain as two codes: a closed-orbit finder using a linear transfer-matrix (Newton-type) search, and a general tracker integrating median-plane-exact equations with the field supplied as tables of Fourier coefficients versus radius, with acceleration switchable on or off.

    B(r,theta) = B0(r) + sum_j [H_3j(r) cos(3j*theta) + G_3j(r) sin(3j*theta)] (field input format, Table I)

    modelingbeam-dynamics dg-1013

    Source, quote & tabletop applicability
    The Fixed Point Code locates closed orbits by means of a highly effective linear transfer matrix procedure, the General Orbit Code tracks arbitrary orbits as desired

    Blosser & Gordon, Computational Study of a Resonant Extraction System for a 3-Sector Cyclotron — MSUCP-9 (1961) — p. 5

    Tabletop: This is the CYCLOPS architecture in embryo (the lineage the planned "CYCLOPS-lite" copies); a next machine's tracker should likewise separate the equilibrium-orbit /tune solver from the general tracker, sharing one Fourier-vs-radius field representation fed by FEMM.

  16. To map the phase-space topology at an energy, locate the unstable fixed points first, then launch orbits displaced slightly from them and integrate both forward AND backward in time — the resulting trajectories trace the separatrices that bound every region of interest, far cheaper than blanketing the plane with orbits.

    beam-dynamicsmodeling dg-1014

    Source, quote & tabletop applicability
    the general orbit code is employed to trace forward and backward in time orbits with initial conditions displaced slightly from the unstable fixed points.

    Blosser & Gordon, Computational Study of a Resonant Extraction System for a 3-Sector Cyclotron — MSUCP-9 (1961) — p. 5

    Tabletop: Directly reusable in a Python tracker (integrate with negative dt for the backward branch); the efficient way to draw the r-pr stability picture of any field candidate for a next machine near resonances.

  17. Median-plane-only tracking is a justified economy: the small axial aperture holds the beam where the field's z-dependence is linear, so off-median aberrations stay small relative to median-plane effects — but the assumption was spot-checked with a few off-plane trial runs, not just asserted.

    beam-dynamicsmodeling dg-1015

    Source, quote & tabletop applicability
    the small axial beam space in a cyclotron constrains the particles to move in a region where the z dependence of the field is quite linear.

    Blosser & Gordon, Computational Study of a Resonant Extraction System for a 3-Sector Cyclotron — MSUCP-9 (1961) — p. 3

    Tabletop: Permission to build the next machine's first tracker 2-D (r, pr, E, phase) and add axial motion as a linearized afterthought — with the same obligation to verify by a handful of full-3D spot checks.

  18. Orbit studies can run on measured scale-model magnet fields long before the machine exists: the B26.1R field came from an 8.75-inch model magnet, radially scaled by 64/8.75 to the full machine, Fourier-analyzed assuming perfect 120-degree symmetry, flutter smoothed of measurement noise, average field isochronized, and harmonics above 99 dropped as negligible.

    r_machine = r_model * (64/8.75); field tabulated at radial increment 0.0080924 cyclotron units (1 cyc unit = E0/(q*B0*c))

    magnetmodeling dg-1016

    Source, quote & tabletop applicability
    The radial spacing of the table entrys is interpreted as increased by the factor 64/8.75 corresponding to the ratio of pole diameters

    Blosser & Gordon, Computational Study of a Resonant Extraction System for a 3-Sector Cyclotron — MSUCP-9 (1961) — p. 6

    Tabletop: The historical analog of the CadQuery->FEMM->field-map pipeline, plus two habits worth copying — symmetrize and smooth measured maps before tracking, and document every cleanup applied to the field the tracker actually ate.

  19. Simplify the accelerating waveform first, validate later: square-wave energy gain was used deliberately to decouple (E,t) from (r,pr) phase space; a closing check with sinusoidal voltage shifted the final beam position but left distortion essentially unchanged, adding only ~30 keV spread across a beam-sized area from differential phase slip.

    sinusoidal check after 8 turns: 62 keV total spread over 5 tracked particles (~30 keV across a beam-sized subarea), 55 deg mean phase drift

    modelingbeam-dynamicsrf dg-1022

    Source, quote & tabletop applicability
    The sinusoidal voltage, it is seen, shifts the final position of the beam spot but has almost no effect on the distortion.

    Blosser & Gordon, Computational Study of a Resonant Extraction System for a 3-Sector Cyclotron — MSUCP-9 (1961) — p. 40

    Tabletop: A permission slip for CYCLOPS-lite staging — start with constant energy gain per gap to get the radial dynamics right, then add cos(phi) gain and phase slip as a second-stage refinement, checking that conclusions survive.

  20. Validate the tracker against hand analytics at every opportunity: the gap-crossing-resonance amplitude (generated because each energy kick shifts the applicable equilibrium orbit while r, pr stay fixed) was computed by hand from tabulated orbit separations and linear mappings, and reproduced the tracked grid's amplitude and phase. The effect needs the two gap kicks to add coherently, which happens only when the field lacks 180-degree symmetry.

    amplitude generated per crossing = -(shift of E.O. between E and E+dE); example chain 0.00126 at 138 deg -> 0.00161 at 28 deg over one turn (Table III)

    modelingbeam-dynamics dg-1023

    Source, quote & tabletop applicability
    the result is seen to fairly accurately predict the actual amplitude and 0 of this point of the grid

    Blosser & Gordon, Computational Study of a Resonant Extraction System for a 3-Sector Cyclotron — MSUCP-9 (1961) — p. 44

    Tabletop: Two lessons — build point analytic cross-checks into a next machine's tracker test suite (transfer-matrix estimates vs tracked orbits), and note the physics is benign for a 180-degree-symmetric two-dee tabletop field where the paired kicks cancel.

  21. Let the computation overrule the folklore: tracking showed beams entering the extraction region centered on the equilibrium orbit behave as well or better than deliberately displaced beams — contradicting the group's own earlier qualitative proposal (MSUCP-2) that displacement would help.

    modelingextraction dg-1024

    Source, quote & tabletop applicability
    beams entering the extraction region approximately centered on the equilibrium orbit behave as well or better than beams entering with substantial displacement.

    Blosser & Gordon, Computational Study of a Resonant Extraction System for a 3-Sector Cyclotron — MSUCP-9 (1961) — p. 45

    Tabletop: The project-level lesson for a next machine — run the cheap simulation before committing hardware to any orbit-dynamics intuition, including intuitions published by people as good as Blosser and Gordon.

  22. A closed-form, size-independent solution exists for the cyclotron dee-gap field (Schwarz-Christoffel, per Murray & Ratner 1953 with corrections): for zero-thickness semi-infinite plate pairs at y = +/-h, tips at x = +/-k, potentials -/+V0, the median-plane field and potential are two-line formulas once one transcendental equation is solved. Geometry caution: k is the HALF-gap and h the HALF-aperture (plate tips map exactly to x = +/-k; re-derived from eq. 1 during extraction — the Fig. 1 scan invites misreading the full gap as k).

    median plane (eqs. 6-8): E_x(x,0) = (V0/h)*sech(X1)/(1 + alpha*sech^2(X1)); V(x,0) = sign(x)*(2*V0/pi)*arccos(sech(X1)); with X = pi*x/(2h) = X1 + alpha*tanh(X1); alpha = (1-a^2)/a^2; a from (pi/2)*(k/h) = arccosh(1/a) + sqrt(1-a^2)/a^2. E_y = 0 on the median plane; E_x even, V odd in x. (Report writes E = +dV/dx — fix sign on implementation.)

    deemodelingbeam-dynamics dg-1025

    Source, quote & tabletop applicability
    This paper presents in summary formulas for the computation of electric fields and potentials of an idealized cyclotron dee geometry.

    Beal, Computation of Electric Field and Potential of an Idealized Dee Geometry — MSUCP-12 (1961) — p. 6

    Tabletop: TRACKER SEED (flagged): this is directly implementable as the gap-field model in the tiny and the next machine's Python trackers — roughly ten lines plus a Newton solve — replacing or validating FEMM electrostatic maps. Identify 2h with the dee aperture, 2k with the dee-to-dummy-dee gap, 2V0 with the full dee-to-dummy-dee voltage (a grounded dummy dee is the same solution shifted by a constant, V0 = V_dee/2).

  23. The same solution gives the full off-median-plane E field — the ingredient needed for electric (gap) focusing models: E_x and E_y anywhere in the aperture follow from two coupled transcendental equations in (X1, Y1). Beal tabulated only the median plane, but eqs. 1-5 contain the whole 2-D field.

    general (eqs. 2-5): E_x = (V0/h)*Xv/(Xv^2+Xu^2); E_y = (V0/h)*Xu/(Xv^2+Xu^2); Xv = cosh(X1)*cos(Y1)*(1 + 1/(A*F^2)); Xu = -sinh(X1)*sin(Y1)*(1 - 1/(A*F^2)); F = sqrt(cosh^2(X1) - sin^2(Y1)); potential v = arccos(cos(Y1)/F), v = pi*V/(2*V0); solve X = X1 + sinh(X1)*cosh(X1)/(A*F^2) and -Y = Y1 + sin(Y1)*cos(Y1)/(A*F^2), with X = pi*x/(2h), Y = pi*y/(2h), A = a^2/(1-a^2).

    deemodelingbeam-dynamics dg-1026

    Source, quote & tabletop applicability
    Therefore, equations 2 and 4 coupled with equations 3 and 5 can be used to determine the electric field and potential at a point X, Y of the dee region.

    Beal, Computation of Electric Field and Potential of an Idealized Dee Geometry — MSUCP-12 (1961) — p. 5

    Tabletop: E_y(x,y) is exactly what a tracker needs for the Rose/Wilson electric gap-focusing term that dominates axial stability on the first turns of a sub-kV machine like the reference machine — available here analytically at any (x,y), no field map required.

  24. Solve the gap-field transcendental equation with Gordon's Newton iteration, not Murray-Ratner's original (which converges slowly for small alpha and fails for large alpha): linearize tanh(X1) about the current guess; convergence is quadratic and works for all gap ratios given the two-branch initial guess.

    eq. 9: X1_new = [X - alpha*tanh(X1*) + alpha*X1**sech^2(X1*)] / [1 + alpha*sech^2(X1*)]; initial guess X1 = X/(1+alpha) if (1+alpha) > X, else X1 = X - alpha

    modelingdee dg-1027

    Source, quote & tabletop applicability
    This iteration is Newtonian in character such that if a given X1 has an error of order e, then X1 given by equation 9 will have an error at order e squared.

    Beal, Computation of Electric Field and Potential of an Idealized Dee Geometry — MSUCP-12 (1961) — p. 7

    Tabletop: Copy the iteration and its initial-guess branch verbatim into the tracker's field routine; a handful of iterations reaches machine precision, cheap enough to call per integration step (or use once to build a spline).

  25. Table 1 is a ready-made verification dataset: E/(V0/h) and V/V0 at x/h = 0 to 5.0 in steps of 0.2, five significant figures, for nine gap ratios k/h = 0.1, 0.3, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5 — computed on MISTIC (D. A. Johnson's fixed-point program), the same machine as the MSUCP-9 orbit codes.

    benchmark anchors: E(0)/(V0/h) = 0.99388 (k/h=0.1), 0.87049 (0.5), 0.65448 (1.0), 0.48916 (1.5), 0.37823 (2.0), 0.21623 (3.5); V/V0 at x/h=1.0: 0.73760, 0.70319, 0.60615, 0.48545, 0.38270, 0.21862 respectively

    modelingdee dg-1031

    Source, quote & tabletop applicability
    Table 1 gives values of electric field and potential for a wide range of dee gap arrangements. [tables span PDF pages 11-19]

    Beal, Computation of Electric Field and Potential of an Idealized Dee Geometry — MSUCP-12 (1961) — p. 11

    Tabletop: Unit-test targets for the tracker's gap-field routine AND an independent check on FEMM electrostatic runs — model the same idealized geometry in FEMM once and match these 5-digit values before trusting FEMM on the real next-machine electrode shapes.

  26. Know the idealization's edges before leaning on it: the solution is 2-D (infinitely long straight edge — no dee-tip curvature, corners, or azimuthal variation), zero plate thickness, semi-infinite plates, electrostatic (quasi-static per RF cycle), no space charge, and symmetric +/-V0 drive. Beal's MISTIC computations covered alpha < 4, i.e. k/h < ~3.77, though the formulas themselves have no such limit.

    alpha < 4 corresponds to k/h < (2/pi)*(arccosh(sqrt(5)) + 4*sqrt(5)/5) = 3.77

    modelingdee dg-1032

    Source, quote & tabletop applicability
    A fixed point computer program written by D.A. Johnson for use on MISTIC was used to calculate the above equations for Ex and Vx at the point (X,0) with a<4.

    Beal, Computation of Electric Field and Potential of an Idealized Dee Geometry — MSUCP-12 (1961) — p. 8

    Tabletop: For use on a next machine, the real deviations to check against FEMM are finite dee thickness, the rounded/blunted tip, and the curved gap line near the source at small radius — expect the analytic solution to be excellent at mid-radius and approximate near center where the ion-source chimney dominates the field anyway.

  27. Shape sector iron by formula-guided empirical iteration, not a priori specification: pick "reasonable" <B>(r) choices, observe the flutter F(r) that results, and test the combination against tune formulae rather than demanding the iron fit pre-selected profiles exactly.

    iterate {<B>(r), F(r), tan(spiral)} -> Smith-Garren vz^2, vr -> accept/reject; do not fix profiles a priori

    magnetbeam-dynamicsmodeling dg-1089

    Source, quote & tabletop applicability
    The process is a trial and error search, with general guidelines and test criteria for success.

    Rainwater et al., The Columbia University Nevis Synchrocyclotron Major Modification — NEVIS-189 / R-774 / CU-295 (1971) — p. 5

    Tabletop: Directly transferable design-process pattern for any pole/shim work on a next machine — let FEMM play the role of the Nevis model magnets, with analytic tune formulae as the accept/reject criteria instead of a target field profile.

  28. Evaluate axial and radial tunes first with analytic formulae (Smith-Garren), then verify at critical places — especially large radius where derivative terms grow — by exact orbit-integration computer solutions.

    analytic vz,vr everywhere; exact orbit codes at critical radii (large r, extraction)

    beam-dynamicsmodeling dg-1090

    Source, quote & tabletop applicability
    first evaluated using the Smith-Garren formula, checked at critical places, especially at larger r, by exact orbit motion computer solutions.

    Rainwater et al., The Columbia University Nevis Synchrocyclotron Major Modification — NEVIS-189 / R-774 / CU-295 (1971) — p. 5

    Tabletop: Exactly the field-solver-plus-orbit-tracker pipeline an amateur design can run; the Nevis precedent says spend the expensive tracking only where the cheap formulae are least trustworthy.

  29. Establish all adjustable RF-system parameters (tuning range, mode spectrum, coupling ratios, voltage distribution) on a reduced-scale model plus computer calculation before committing to full-scale construction; also verify the beam-excited cross mode stays clear of harmonics of the main mode.

    1/2-scale RF model + computation -> full-scale build; cross mode kept well below 2x main mode over tuning range

    rfmodeling dg-1098

    Source, quote & tabletop applicability
    The design has used a 1/2 scale model, in conjunction with detailed computer calculations, to establish all parameters which can be varied

    Rainwater et al., The Columbia University Nevis Synchrocyclotron Major Modification — NEVIS-189 / R-774 / CU-295 (1971) — p. 9

    Tabletop: At tabletop size the "scale model" is the full-size mockup on the bench — cold-test a next machine's dee/stem with a VNA before power exists; the mode-spectrum audit transfers verbatim.

  30. When a calculation needs an empirical constant (here the effective image factor of saturated pole iron), measure it directly with a precisely known conductor configuration in the real field environment rather than taking a handbook value.

    septum fields = conductors + 5 image sets scaled by (mu-1)/(mu+1); measured fit gave mu = 5 to <1%

    magnetmodelingbeam-measurement dg-1103

    Source, quote & tabletop applicability
    The value of mu used was found experimentally by measuring the field from a precisely known configuration of conductors

    Rainwater et al., The Columbia University Nevis Synchrocyclotron Major Modification — NEVIS-189 / R-774 / CU-295 (1971) — p. 13

    Tabletop: A model-calibration pattern the FEMM-based pipeline should copy — one deliberate known-geometry measurement (a wire loop, a known coil) in the actual gap pins the permeability/saturation assumptions the whole field model rests on.

  31. Scale-model law for RF resonators: a 1/2-scale model runs at 2x frequency with L and C halved, its Q is 0.7x (1/sqrt(2)) the full-scale Q, and it therefore needs 1.4x the proportional power for a given dee voltage — coupling taps and loops in the model must pick up 1.4x the relative voltage.

    f_model = s*f_full, L,C /s, Q_model = Q_full/sqrt(s), P_model = sqrt(s)*P_full for equal V (s = scale factor 2 for half scale)

    rfmodeling dg-1125

    Source, quote & tabletop applicability
    The Q of the model will be 0.7 times the Q of the actual installation and so will require 1.4 times as much power for a given dee voltage.

    MacKenzie, Preliminary Report on the “Three Quarter Wave” R.F. System for Frequency Modulated Cyclotrons — AECD-1850, University of California (1947) — p. 6

    Tabletop: Bench-model a dee-stem or resonator geometry at reduced size before cutting full-size copper; apply the sqrt(scale) Q correction before comparing model power and coupling measurements to full-scale predictions.

  32. Pick the model scale so a real, available tube is a valid stand-in for the power tube: MacKenzie abandoned a 1/4-scale model because electron transit time at the scaled frequency falsified oscillator behavior, then chose 1/2 scale where an Eimac 304TL (1200 V, 600 mA) scales the 9C21 (12 kV, 6-8 A) faithfully.

    model frequency must stay low enough that tube transit-time effects remain negligible

    rfmodeling dg-1126

    Source, quote & tabletop applicability
    It was not excited satisfactorily due to the fact that at 100 megacycles the transit time effects were quite noticeable on the fundamental mode.

    MacKenzie, Preliminary Report on the “Three Quarter Wave” R.F. System for Frequency Modulated Cyclotrons — AECD-1850, University of California (1947) — p. 6

    Tabletop: Cold measurements (network analyzer) scale exactly, but any POWERED model test needs a driver whose parasitics and transit time scale with the geometry, or the model will exhibit modes the real system never sees and vice versa.

  33. Where geometry is complicated, trust model tests over calculation: MacKenzie preferred a mechanically awkward layout that could only be settled empirically, noting that dimensions calculable "fairly exactly" were the sole advantage of the calculable variant.

    rfmodeling dg-1127

    Source, quote & tabletop applicability
    dimensions can be calculated fairly exactly whereas in the system shown in Figure 5 one must depend on model tests (which are safer anyway).

    MacKenzie, Preliminary Report on the “Three Quarter Wave” R.F. System for Frequency Modulated Cyclotrons — AECD-1850, University of California (1947) — p. 6

    Tabletop: Transmission-line formulas ignore end effects, bends, and support hardware; for any resonator whose geometry is not a textbook line, a cheap model measurement outranks the calculation it checks.

  34. Predict full-scale RF power from model measurements and state both numbers: 400 W (48 Mc) and 600 W (18 Mc) of model input for 1500 V on the dee scaled - via V-squared and the sqrt(2) Q correction - to 28 and 42 kW for 15 kV; MacKenzie also flags that the model's bad joints and brass surfaces bias it pessimistic against the copper full-scale build.

    P scales as V^2 x (Q_full/Q_model)^-1; model bad joints/brass make prediction conservative

    rfmodeling dg-1137

    Source, quote & tabletop applicability
    The 1/2 scale model uses about 400 watts input to the oscillator at 48 megacycles and 600 watts input at 18 megacycles to produce 1500 volts on the dee.

    MacKenzie, Preliminary Report on the “Three Quarter Wave” R.F. System for Frequency Modulated Cyclotrons — AECD-1850, University of California (1947) — p. 16

    Tabletop: Dee power scales as voltage squared - measure watts-per-volt-squared on the bench, correct for Q, and the amplifier requirement for any target dee voltage falls out; joint quality and surface material shift Q enough to budget for.

  35. Match the probe to the field structure: sample the field at a point, not an average — the model-survey coils were 0.20 in. diameter x 0.15 in. high (2000 turns of No. 46, ~350 turn-cm^2) so that 0.5 in. resolution on the full-scale magnet (1/32 in. on the 1/16 model) was preserved, and coil placement errors were held below 1/32 in. Verify by repeat runs that the finite coil size does not degrade the map.

    probe dimension << field-structure scale / model scale factor

    beam-measurementmodeling dg-1307

    Source, quote & tabletop applicability
    It was desirable to know the magnetic field accurately to within 0.5 in. on the full-scale magnet. This corresponded to 1/32 in. on the models.

    Wakerling & Guthrie (eds.), Magnets and Magnetic Measuring Techniques — TID-5215, Radiation Laboratory, University of California (1949) — p. 73

    Tabletop: The active-area rule for Hall mapping a shim edge: a 1-2 mm sensor is marginal where the gradient scale is a few mm (pole edge, shim step); position repeatability of the mapper jig belongs in the same error budget as the sensor.

  36. The model-magnet scaling law: a linear scale model built from steel with the same magnetic properties, operated at the same field strength (same B everywhere, currents scaled to keep NI per gap-length), reproduces the prototype's field distribution exactly — magnetostatics has no intrinsic length scale until saturation properties differ. Leakage coefficients measured on the model apply directly to the full-scale magnet; forces follow with area (L^2) scaling.

    geometric scaling at fixed B and fixed material B-H curve; L_leakage(model) = L_leakage(full scale)

    modelingmagnet dg-1311

    Source, quote & tabletop applicability
    a linear scale model built from steel with the same magnetic properties as planned for the prototype magnet and operated at the same field strength will give results directly applicable to the prototype.

    Wakerling & Guthrie (eds.), Magnets and Magnetic Measuring Techniques — TID-5215, Radiation Laboratory, University of California (1949) — p. 25

    Tabletop: The physics that lets FEMM stand where models stood — and the terms of validity are the same for both: correct B-H data and correct geometry. Any cheap sub-scale mock-up of a planned magnet obeys it too, provided the steel matches and B is held, not NI.

  37. Match model scale to question precision: 1/16 scale answered every production design question (excitation, end effects, forces, flux allocation, stray field), but for the finest field-uniformity contour maps the team deferred to the 1/8-scale model of the same pole geometry as inherently more accurate. Bigger model only where the question demands it.

    modelingmagnet dg-1312

    Source, quote & tabletop applicability
    since the X Beta model was built to 1/8 scale it seemed true that the results would be more accurate than those which could be obtained on a 1/16-scale model.

    Wakerling & Guthrie (eds.), Magnets and Magnetic Measuring Techniques — TID-5215, Radiation Laboratory, University of California (1949) — p. 149

    Tabletop: The mesh-refinement decision in physical form — coarse resolution for excitation/force/leakage questions, fine resolution only for the ppm-level uniformity region; spending fine-model effort on questions the coarse model already answers is waste in either medium.

  38. The standard model-test suite, in reporting form: (1) saturation curve H_g vs NI/l_g; (2) efficiency vs NI/l_g; (3) leakage coefficient at various points (directly applicable to full scale); (4) uniformity contour maps of (H-H_g)/H_g on normalized pole coordinates; (5) stray-field map. Plus, in practice: gap-to-gap comparison, flux audit of every iron member, and magnetic forces. This is the complete characterization a predictive model owes the design.

    report H_g(NI/l_g), eta(NI/l_g), L(x), (H-Hg)/Hg contour map, stray map

    modelingmagnet dg-1313

    Source, quote & tabletop applicability
    The usual measurements made included a magnetization curve, a comparison of gap performance at different points in the magnet, uniformity contour maps, determination of magnetic forces, the density of flux through various parts of the magnet

    Wakerling & Guthrie (eds.), Magnets and Magnetic Measuring Techniques — TID-5215, Radiation Laboratory, University of California (1949) — p. 26

    Tabletop: A ready-made deliverables checklist for a simulation campaign on a new magnet — a FEMM study that produces these five plots plus a member-by-member flux audit has done what the 1944 model program did, in the same order.

  39. Track efficiency (gap mmf / total mmf) at TWO field levels as the saturation health check: the revised Alpha II model measured 95.6 +/- 2.0 per cent at 4600 Oe and 95.4 +/- 2.0 per cent at 3400 Oe — equality within error at both excitations demonstrated the iron was nowhere near saturation and the design conservative. Falling efficiency with rising field is the first global symptom of a saturating member.

    eta = 2.02 * H_avg[G] * l_gap[in] / (NI); compare at two field levels

    magnetmodeling dg-1314

    Source, quote & tabletop applicability
    With 3400 oersteds in the gaps the efficiency obtained was 95.4 +/- 2.0 per cent. Within experimental error they were the same at both field strengths. This indicates that the magnet design is rather conservative.

    Wakerling & Guthrie (eds.), Magnets and Magnetic Measuring Techniques — TID-5215, Radiation Laboratory, University of California (1949) — p. 140

    Tabletop: A two-point excitation scan (measured or simulated) separating "efficiency constant" from "efficiency dropping" localizes saturation onset without any interior probe — directly usable on an H-frame by comparing measured B vs I against the linear NI prediction at two currents.

  40. Audit the flux through EVERY iron member with wound loops and a ballistic integrator: loop-flux differences divided by enclosed-area differences give local leakage flux components (horizontal and vertical separately, by choosing loop pairs); dividing member flux by member steel area gives its working induction. Alpha II verdicts: core steel 12,000 G comfortable; core rim 17,460 G too high (fixed by the thicker full-scale rim, ~15,000 G, mu ~ 550); yokes at 10,000-15,000 G "good"; deliberately sacrificial cooling-tank walls saturated at ~30,000 G.

    B_member = (phi_loop difference)/(A_steel); leakage component = d-phi/d-A between loop pairs

    magnetmodeling dg-1315

    Source, quote & tabletop applicability
    In the neighborhood of 10,000 to 15,000 gauss the flux density is good, yet the yokes are not overloaded to the extent that the permeability of the steel drops excessively.

    Wakerling & Guthrie (eds.), Magnets and Magnetic Measuring Techniques — TID-5215, Radiation Laboratory, University of California (1949) — p. 141

    Tabletop: The 10-15 kG working band for structural mild steel is the same number modern small-magnet guidance gives (cf. Wouters; Zickler CAS) — and the loop-audit method is the measurement twin of integrating B over member cross-sections in a FEMM postprocessor: every member gets a number, every number gets a verdict.

  41. Correct model predictions for known model/prototype differences, with signs stated: the team measured permeability of BOTH the model steel and the full-scale steel (full-scale better -> full-scale performs slightly better), and tallied deliberate geometry differences (model rim proportionally 1 in. thinner -> model worse; model coil-tank iron 0.5 in. thicker -> model worse). Every known discrepancy got a direction, so the prediction became a bound, not a guess.

    modelingmaterials dg-1317

    Source, quote & tabletop applicability
    The permeability of the material of the full-scale unit is higher than that of the model. This indicates that a slightly better performance could be expected from the full-scale unit than from the model.

    Wakerling & Guthrie (eds.), Magnets and Magnetic Measuring Techniques — TID-5215, Radiation Laboratory, University of California (1949) — p. 146

    Tabletop: The sign-audit habit transfers to simulation directly — list every model-vs-hardware difference (B-H table provenance, fillets, packing factor, gaps at joints) with the direction it biases the prediction, so measured-vs-predicted discrepancies arrive pre-explained.

  42. The calutron model program's validation verdict — the benchmark for trusting scaled prediction: full-scale tests confirmed the 1/16-scale models as "dependable and accurate," with full-scale performance slightly BETTER than predicted (source-region field more uniform than model results, stray field weaker, track efficiency ~94% vs 95.6 +/- 2% model, core-to-tank field concentration 18% vs model 22%); 71 of 71 production tanks met the theoretical field criteria, worst case 2.8 cm against a 3.0 cm limit.

    modelingmagnet dg-1322

    Source, quote & tabletop applicability
    The magnetic performance of the track is better than predicted from the model experiment.

    Wakerling & Guthrie (eds.), Magnets and Magnetic Measuring Techniques — TID-5215, Radiation Laboratory, University of California (1949) — p. 191

    Tabletop: The historical calibration point for any predict-then-verify magnet pipeline: a faithful scaled model (same-steel/same-B then; validated FEM now) lands within a few per cent on global quantities and errs conservative when the prototype''s iron is better than the model''s — but the trust was EARNED by one full validation campaign, not assumed.

  43. After the first article validates the prediction chain, degrade acceptance testing to mechanical metrology: once track 1's 71 tanks all passed magnetic tests and shim-position measurements were shown sufficient to guarantee the field, track 6 was accepted on a dimensional check of shim positions alone, with magnetic spot checks only for special questions (end-tank asymmetry).

    modelingproject-management dg-1323

    Source, quote & tabletop applicability
    the excellent results obtained in testing track 1 showed that a dimensional check of the shim positions was entirely adequate.

    Wakerling & Guthrie (eds.), Magnets and Magnetic Measuring Techniques — TID-5215, Radiation Laboratory, University of California (1949) — p. 191

    Tabletop: The economic payoff of validation: once field-vs-geometry is established (by model/simulation plus one measured article), later shim changes can be accepted on caliper and indicator readings, reserving full field maps for genuinely new configurations.

  44. Expect field CORRECTION to be trial and error, and budget for it: the theoretically grounded shim-tilt correction scheme failed validation (predicted and measured tilt effects disagreed near the tilted shim, partly because a theory assumption — iron stuffing behind the tilted shim — was not implemented in hardware), and the team concluded the only feasible correction method for out-of-criteria fields was iterative cut-and-try.

    magnetmodeling dg-1324

    Source, quote & tabletop applicability
    It would appear that the only feasible method of making corrections when the necessity arises is by trial and error.

    Wakerling & Guthrie (eds.), Magnets and Magnetic Measuring Techniques — TID-5215, Radiation Laboratory, University of California (1949) — p. 181

    Tabletop: A 1944 warning that survives every FEMM run: analysis predicts the as-built field well but predicts CHANGES to an as-built field only if the change is modeled exactly as executed. Plan shimming as measure-cut-measure iterations, and keep shim stock adjustable.

  45. A scale-model magnet is a close call — list the reasons before building one. UW's five advantages of the 1/12 model: cheap/fast shim iteration (especially if special alloys needing heat treatment were tried), work proceeds before the full magnet is done, no interference with other construction, convenient minor measurements, and future re-studies while the cyclotron operates. The two difficulties: spatial resolution of field measurement, and coil heat (current density scales as 1/L, heat per unit volume as its square). Verdict after the fact: "advantages and disadvantages ... were fairly closely balanced" — partly because part geometry constrained shim options more than expected.

    Scaling at constant B: J ~ 1/L; heat/volume ~ J^2 ~ 1/L^2

    magnetmodeling dg-1332

    Source, quote & tabletop applicability
    the amount of testing required to arrive at a final shim design was much less than anticipated. This was due in part to the fact that the geometry of parts limited the possible variations more closely than was expected.

    The University of Washington 60-Inch Cyclotron: Progress and Status Report of Design and Construction — AECU-1951, University of Washington (c. 1950) — p. 17

    Tabletop: With FEMM the model-magnet role is filled by simulation (ucrl-31 showed the scale-model method itself; MacKenzie AECD-1850 the model-test discipline), but the balanced verdict is the lesson — physical iteration budget should go where the computable model is least trustworthy (saturation, real steel, mechanical tolerances).

  46. If you do build a model magnet, pour it from the same heat as the full core and make it a precise replica: UW's 1/12 model used Midvale forgings "poured from the same heat and it can be assumed magnetic properties are identical," a precise replica except bolts and carrying lugs (227 lb), with cover plates made from scraps of the actual cover plate stock.

    magnetmodelingmaterials dg-1333

    Source, quote & tabletop applicability
    The steel for both the cyclotron magnet and the model was poured from the same heat and it can be assumed magnetic properties are identical.

    The University of Washington 60-Inch Cyclotron: Progress and Status Report of Design and Construction — AECU-1951, University of Washington (c. 1950) — p. 18

    Tabletop: The transferable rule is identity of material between test article and final article — a next machine's FEMM model should use a B-H curve measured on the actual purchased steel, not a library curve for its nominal grade.

  47. Verify the model's prediction on the full magnet before committing to shims: UW measured the unshimmed full-scale radial dependence first, found "close agreement" with the model, and only then cut cyclotron shims to the model design — after which "the predicted radial dependence was verified and the results were considered satisfactory." Prediction, cross-check, commit.

    magnetmodeling dg-1339

    Source, quote & tabletop applicability
    showed that the model data could be used as a basis of prediction with confidence.

    The University of Washington 60-Inch Cyclotron: Progress and Status Report of Design and Construction — AECU-1951, University of Washington (c. 1950) — p. 33

    Tabletop: The FEMM-era version — survey the bare magnet, reconcile with the simulation, THEN machine shims from the reconciled model. Same model-then-verify discipline as MacKenzie's aecd-1850.

  48. Model the RF system at quarter scale (4x frequency) before building it: UW's quarter-scale model of the entire resonant system validated the calculated line lengths ("sufficiently accurate") and needed only minor adjustments; measured model Q ~3500 with no special joint precautions — about half the full-scale expectation, as scaling predicts (Q ~ sqrt(scale) at fixed geometry). Calculated equivalent-circuit constants were treated as guides, with the report noting it is "sometimes desirable" to add capacitance at the tube in parallel with the interelectrode capacitance to make a practical stub line.

    rfmodeling dg-1354

    Source, quote & tabletop applicability
    there are many uncertainties in the exact determination of the constants of the equivalent circuit ... the calculated values ... serve well as a guide

    The University of Washington 60-Inch Cyclotron: Progress and Status Report of Design and Construction — AECU-1951, University of Washington (c. 1950) — p. 82

    Tabletop: A tabletop resonator IS the scale model — build the dee/stem mockup on the bench, measure f and Q before committing to vacuum hardware, and trust lumped calculations as guides not gospel. Berkeley used the same quarter-scale method on the 88-inch (ucrl-9435), which also confirms the Q-degradation-at-joints lesson (ornl-2648).

  49. Kill parasitics on paper first: compute/measure the higher modes of the anode and grid circuits and ADJUST CIRCUIT ELEMENTS so that no mode coincides with a class-C plate-current harmonic anywhere in the tuning range. Harmonic content falls roughly as 1/n, so only the low harmonics (2nd, 3rd) can excite a mode destructively; the 88-inch verified mode placement on a quarter-scale RF model and reports the consequences of getting it wrong — destructive voltages at the grid vacuum insulator and reduced fundamental output.

    Design constraint: f_mode(k) != n * f_osc for n = 2, 3 over the whole tuning range (harmonic amplitude ~ 1/n)

    rfmodeling dg-1369

    Source, quote & tabletop applicability
    The circuit elements of the rf system were adjusted so that the first two higher modes would not be excited by an oscillator harmonic.

    Smith, The RCA 6949 as a Self-Excited Cyclotron Oscillator — UCRL-9435, Lawrence Radiation Laboratory (1960) — p. 5

    Tabletop: DIRECT for a fixed-frequency machine: sweep the dee system to a few hundred MHz on a VNA, list the modes, and check none sits at 2f or 3f of the drive — if one does, detune it with a strap or stub BEFORE blaming the amplifier for instability. Same mode-vs-harmonic discipline as the msucp-9 and ornl-2648 lines.