Charge carriers have a tiny average drift inside the continuous wire. In a metal, negative electrons drift opposite the defined conventional-current arrow. They do not leap across the space between neighbouring windings; the wire, cryostat, and patient do not circulate.
Theme 02 · Fields and radio
Follow current into B₀, map and shim its variation, then use gradients and RF.
Manipulate 3D MRI magnet fields, perform a quantitative low-order B0 shim fit, then explore gradients, proton precession, RF transmit and receive hardware, I/Q detection, and a full-Bloch slice profile.
00 / WHERE THE MAIN FIELD COMES FROM
Cold wire carries current.
The current creates B₀.
A power supply first pushes electric current through many turns of superconducting wire. Those turns act together like one long solenoid: their magnetic fields add through the bore and return outside the magnet. After ramp-up, a closed superconducting path can keep that current circulating in persistent mode.
Drag to orbit · pinch, wheel, or use + / ↺ / − to zoom · tap any object · field lines are a map, not physical tubes
- 01Ramp the current
A controlled DC power supply applies voltage so current rises in the cold winding. “Current” means net electric charge crossing a wire section each second, measured in amperes. Yellow arrows follow the conductor in the conventional positive-charge direction; the wire does not move, and electrons do not jump between turns.
- 02Every turn contributes
Each loop makes a magnetic field. Inside a long coil, neighbouring loop fields point mainly the same way and add; more amperes or more turns per metre gives a larger central field.
- 03Close the persistent path
Once the target current is reached, an internal superconducting connection completes a very-low-resistance loop. The ramp supply can be disconnected while the established current continues.
- 04B₀ stays on
The stable main field exists between scans. RF and gradient pulses switch during imaging; the superconducting main-magnet current normally does not pulse on and off for each patient.
Bcenter ≈ μ0nI = μ0(N/L)I
4π × 10⁻⁷ T·m/A × 4,000 turns/m × 600 A ≈ 3.02 TAt 600 A and the selected winding density, the ideal long-solenoid estimate is 3.02 T. Raising current strengthens B₀ in direct proportion. Yellow arrowheads stay on the helical wire and mark conventional-current direction; they are not electrons and their display speed is not a measurement.
The winding length stays 1.50 m inside the 1.86 m package. The 3D winding shows 52 spaced turns to represent 6,000 effective turns. Move the slider: higher n visibly adds turns and closes their gaps; the drawing remains a compressed sample rather than literal manufacturer geometry.
This changes only how many direction-map curves the 3D lesson draws. It does not change current, B₀, tesla, gauss, fringe field, stored energy, or any clinical image.
More ampere-turns make a stronger main field.
At the center, increasing I or n increases the field almost linearly in this long-solenoid estimate. Clinically, a different B₀ changes the ¹H carrier frequency and many scanner/tissue behaviors; it is a magnet design and site platform, not a routine sequence knob.
Drawing lines connect the field direction at many points. They close into loops and crowd where the drawing represents stronger field. They are not strings, rays, current paths, or proton tracks.
Both measure magnetic flux density: 1 tesla = 10,000 gauss. Therefore 3 T = 30,000 G. Clinical MRI names B₀ in tesla; gauss is common in fringe-field discussions.
Ordinary resistance would turn sustained high current into continuous heat. Below its critical conditions, the magnet conductor can support persistent current with extremely small decay; cooling and protection remain essential.
3D FIELD BUILD-UP · MOVING CHARGE → CURRENT → LOOP → SOLENOID
One conductor makes a circular field.
Many loop fields add.
Blue moving markers now represent groups of negative charge crossing a counting plane; increasing current shows more marker crossings and a stronger calculated field. They are not a literal census of electrons in the wire. Because electrons are negative, their average drift is opposite the yellow conventional-current arrow. Build from a straight segment to one loop, then watch the field shape lengthen as aligned turns are added.
Drag to orbit · pinch, wheel, or use + / ↺ / − · tap an object · animation pace and line count are teaching choices
Count charge crossing a wire section—not electrons merely sitting in the wire.
The side view is kept two-dimensional because “crosses this plane during one second” is a measured rate. The rotatable model beside it shows the same conductor and circular field in 3D.
I = ΔQ / Δt = Ne|e| / Δt
Ampere measures charge crossing per second. Stationary charge gives no steady current-created B; current through the chosen geometry sets the calculated field.At 4.00 A, 4.00 coulombs of charge cross the chosen section each second—equivalent to about 2.50 × 10¹⁹ elementary electron charges per second. The conductor already contains vastly more carriers; this control changes net flow rate, not the mere existence of electrons.
More charge crossing per second strengthens every B contribution in direct proportion. The number of bright blue packets is a compressed rate display; their visible speed is deliberately not a microscopic drift-speed measurement.
Moving this control opens the many-loop stage. More equally oriented loops add more same-direction axial field at the center; it does not mean the 3D drawing represents a real MRI winding count.
Vector addition, not field-line counting
QUANTITATIVE AXIAL CROSS-SECTION · SAME CURRENT AND WINDING CONTROLS
Color tells how much field.
Arrows tell which way it points.
The 3D guides above are useful direction guides, but their spacing is not a measurement. This map calculates a field value at a grid of positions through the coil. Drag the red probe directly: moving it changes only the observation coordinate, never the magnet.
FINGER / POINTER: drag anywhere in the plot · KEYBOARD: focus the plot, then use arrow keys · the ordinary z and r sliders remain equivalent accessible controls
NEW QUANTITATIVE 3D LAB · MAPPED Δf → LOW-ORDER CORRECTION → RESIDUAL
Can broad shim fields flatten every B₀ error?
A field map describes where proton resonance is above or below a chosen reference. Fit the standard first- and second-order spatial shapes over one spherical region, then compare the original map, the added correction, and what the fit cannot remove.
Drag to orbit · pinch, wheel, or use + / ↺ / − · tap a field volume · needles plot signed Δf along the B₀ axis; they are not tissue displacement
Why can broad shim fields reduce—but not erase—a localized field defect?
All three 3D volumes show the same coordinates and the same fixed hertz scale. Change only the fitted basis, region size, or local feature, then ask which spatial patterns the residual keeps.
Increase model order without changing the measured map.
The fit uses 257 equally weighted Cartesian teaching samples inside a 22 cm sphere. Changing the region changes which physical positions contribute; it does not resize a patient or claim an optimal clinical shim volume.
This synthetic local feature is not anatomy or a measured susceptibility field. Raising it adds spatial detail that eight broad harmonic shapes cannot reproduce exactly.
The residual is the map plus the fitted correction.
Δfres(r) = Δfmap(r) + Σ cj Φj(r)
Δf is proton off-resonance in hertz and Δf = γ̄ΔBz. Φj are explicit dimensionless X, Y, Z, Z², ZX, ZY, X²−Y², and 2XY shapes normalized at a 12 cm reference radius. cj are correction-field coefficients in hertz—not shim-coil current.Which broad shapes oppose the map?
Signs describe the field added by this fit in the stated coordinate convention. A real scanner maps calibrated coil currents and constraints to fields; this lesson intentionally stops at field coefficients.
The low-order part shrinks; the local feature survives.
The least-squares correction lowers RMS off-resonance from 19.31 Hz to 4.77 Hz over the same samples. The residual is not a failure of arithmetic: the local Gaussian shape is deliberately absent from the allowed basis.
+Z is parallel to nominal B₀. Camera orbit changes only the view.
Center-frequency adjustment is separated from spatial shim coefficients.
Color repeats sign with a visibility floor plus square-root boost; shells, thickness, and labels are display geometry.
This is not a measured patient field map. The map is deterministic synthetic data sampled on an unweighted Cartesian lattice inside one sphere. The local Gaussian is an illustrative hard-to-fit feature, not modeled anatomy. The fit is unconstrained least squares over ideal normalized field shapes; it omits measured coil calibration, current and heating limits, cross-terms, concomitant/transverse fields, passive shims, higher-order or multi-coil arrays, regularization, motion, dynamic shimming, slice-wise optimization, noise, drift, and sequence-specific acceptance criteria. Near-zero residual for the low-order-only preset means the synthetic map was constructed from the same basis—not that a real subject can be shimmed perfectly.
Zoom from the patient to the nuclear physics.
Use the numbered journey in order. Every solid, sphere, arrow, field sheet, and ring can be tapped for its exact meaning. The jumps between body, tissue, atom, and nucleus are separate scale models—nothing literally swells inside the patient.
Drag to orbit · pinch, wheel, or use + / ↺ / − · tap an object to identify it · model sizes are not one continuous scale
First choose where inside the body the signal comes from.
The scanner surrounds the patient with B₀, gradient, and RF hardware. The yellow target marks one tiny tissue region inside the head. The next buttons magnify that region conceptually; the tissue does not move out of the body.
Patient, table, scanner bore, tissue water, fat, and many other materials.
A glowing target and connector make one future region easy to follow.
MRI never receives a ready-made picture from one location; it must excite and encode signals distributed through the body.
At 3.0 T, ¹H phase precesses at 127.73 million cycles per second. The 3D arrow is slowed enormously; its animation speed is only a trend cue.
The selected stage assigns this field slope to physical Z for slice selection, Y for phase encoding, or X for readout. Gradient current changes precession rate with position; it does not push atoms along that axis.
At 20.0 mT/m, a 2.0 kHz RF band selects an ideal 2.35 mm slab centered at z = 0 mm. Moving RF center frequency moves the slab; it does not move atoms.
f0 = γ̄B0
Δz = BWRF / (γ̄|Gslice|)
dk/dt = γ̄G(t)
vcoil(t) = −dΦ(t)/dt
A proton has intrinsic quantum angular momentum and a magnetic moment. Those are measured properties, but they do not require a little material surface spinning like Earth. The 3D arrow represents magnetic-moment direction or an ensemble-average magnetization. Its cone represents precession of that direction around B₀. Clinical MRI detects the collective changing field from vast ensembles—not one isolated proton.
The patient, tissue cube, water molecule, electron cloud, proton, spin-packet grid, receive loop, and voltage trace cannot share one honest geometric or numerical scale: body dimensions are metres, encoding positions centimetres, molecules nanometres, nuclei femtometres, and receive signals often microvolts. Each stage resets scale and camera. Water bonding, electron probability, nuclear spin, thermal polarization, quantum energy levels, Bloch-equation ensemble motion, slice profiles, relaxation, diffusion, susceptibility, chemical shift, exchange, coil sensitivity, noise, hardware fields, and clinical safety are simplified. The printed formulas state what is quantitative; colors, object counts, arrow sizes, cone angles, trace amplitude, and animation rates do not.
01 / FROM FIELD TO POSITION
A controlled tilt in the magnetic field.
A gradient coil adds a small, nearly linear variation to the much larger main field. Proton frequency now carries an address: where a spin sits determines how fast its phase turns.
Bz(x) = B0 + Gxx
Δf = γ̄Gxx → Δφ = 2πΔfΔt
kx = γ̄GxΔt → Δφ(x) = 2πkxx
PHASE ADDRESSAt x = +8.0 cm, Δf = +85.15 kHz accumulates +0.341 cycles (+122.6°) in 4.00 µs.
f0 = γ̄B0
42.58 MHz/T × 3.0 T = 127.73 MHzAt fixed G, position, B₁⁺, and pulse duration, changing B₀ retunes the MHz carrier. The gradient offset, ideal flip angle, and dk/dt stay unchanged.
Scope: this control retunes the field/RF frequency model. The TE/TR tissue constants remain the explicitly stated illustrative 3 T model; real coil fields, SAR, and relaxation are field-dependent.Bz = B₀ + Gₓx, Δf = γ̄Gₓx, kₓ = γ̄GₓΔt, and Δφ(x) = 2πkₓx are calculated exactly for ¹H under the displayed constant rectangular-gradient preview. Positive Δf is drawn counterclockwise about +Z in this declared convention. The transparent planes are equal-frequency locations; they are not RF sheets. Graph height and color magnify ΔB because its millitesla variation is visually negligible beside B₀ in tesla. The phase hands use a frame rotating at the isocenter carrier and show one declared time snapshot—not literal MHz animation, individual protons, gradient-coil geometry, nonlinear fields, concomitant fields, eddy currents, slew, diffusion, relaxation, or an acquired image.
The gradients are tiny beside B₀, but they are switched quickly and precisely. The audible knocking in MRI is gradient hardware moving under Lorentz force.
BEGINNER GAME / THREE HIDDEN WATER BUCKETS
You are the operator.
Find water with fields.
Three water buckets are hidden inside the bore. From the next room, choose a gradient vector, tune a narrow RF frequency offset, send a short probe, and listen for the echo. Find each bucket’s X, Y, and Z coordinate—then see why a frequency under one gradient selects a plane rather than a unique point.
Drag to orbit · solid ring = bore · blue plane = same-frequency locations · buckets reveal only after detection
Start with X. Sweep or tune across the spectrum and record three peaks; repeat for Y and Z. Peak height identifies the 0.6 L, 1.0 L, and 1.4 L buckets in this deliberately simplified game.
G adds a tiny, known B₀ slope. Positions on one perpendicular plane now share one resonance-frequency offset.
A short B₁⁺ pulse contains a center frequency plus a finite frequency band. Water inside that band is excited; the RF wave is not a locating beam.
After switching and sequence timing, excited water can contribute an echo. The loop detects the total induced voltage; software compares the result with the known field map.
WHAT THE FREQUENCY KNOB CHANGESIt changes RF carrier cycles per second around f₀, selecting another frequency band—and therefore another slab while G is on. It does not change RF amplitude, gradient strength, water amount, or directly select one 3D point.
No probe sent yet. Tuning changes the highlighted plane; press Probe & listen to make a measurement.
Match peak height across X, Y, and Z.
With Gx on, frequency tells X but says nothing about Y or Z. Rotate the gradient job and repeat. Real MRI performs a much richer, systematic version with many phase and frequency encodings.
No gradient: no position clue
With G = 0, all three identical-water resonances sit near the same f₀. Frequency can say “water signal exists,” but not where it came from.
One gradient: one projection
Δf = γ̄G·r. A measured frequency gives position along G, while an entire perpendicular plane has that same frequency.
More encodings: reconstruct location
Real imaging collects many known phase/frequency patterns. Reconstruction solves the mixed whole-object voltages for a spatial image; an operator does not hunt voxel by voxel.
02 / RF TRANSMIT, RELAXATION & RECEIVE
Excite broadly.
Listen locally.
The transmit coil creates the rotating B₁⁺ field that tips magnetization. After excitation, the scanner stops transmitting and receive coils detect the tiny voltage induced by precessing transverse magnetization. A common 1.5 T / 3 T workflow is body-coil transmit with a close local array receiving.
RF WAVEFORM DECODER · THREE DIFFERENT SHAPES, THREE DIFFERENT JOBS
Frequency is horizontal spacing. B₁⁺ amplitude is height. Pulse duration is width.
The plots share the live controls below. They are separated because “a faster wave,” “a stronger pulse,” and “a larger received signal” do not mean the same thing.WHY THIS PLOT IS 2DLeft–right is elapsed time. Up–down is one field component or received-voltage magnitude. Neither axis is a patient direction, so a third spatial axis would be invented and misleading.
The plotted field goes from a positive peak, through zero and a negative peak, back to the next positive peak: one complete 360° oscillation. It is a cycle of the RF magnetic field—not a proton traveling in a circle.
The slow outline connecting the peaks of thousands of fast carrier cycles. It is not a second radio wave and not a container; it summarizes how strongly the transmitter is driven over time.
The size of the complex voltage after the receiver mathematically removes the MHz carrier. Phase is stored too, even though this row shows only non-negative magnitude.
B₀ sets how many RF-field cycles occur each second. B₁⁺ and pulse duration set the ideal flip angle. Coil geometry changes the receive-sensitivity proxy.
Cycles pack closer or farther apart. The RF coil must retune, but ideal flip stays fixed here when B₁⁺ and duration are held.
The carrier and envelope grow taller. Flip angle changes; carrier spacing does not. Excess RF power can increase heating and SAR.
The envelope grows wider and contains more cycles. Flip angle changes; cycles per second do not.
Only the geometry-only receive curve changes here. The transmitter’s carrier and ideal flip do not.
How one fast RF coil voltage becomes I and Q.
Before demodulation there is one rapidly alternating voltage. The receiver compares it with two synchronized reference waves; low-pass filtering leaves two slower signed voltages that preserve the signal’s magnitude and phase.
WHY THIS EXPLANATION USES 2DThe wave rows plot voltage versus time; the round I/Q view plots two mathematical receiver-reference components. I and Q are not scanner X and Y. Use “See this split in the 3D receiver” to locate the calculation in the hardware chain.
One signed RF voltage
Changing magnetic flux induces a voltage that alternates above and below zero near the ¹H carrier. At the coil it can be microvolt-scale; gain changes its voltage scale but not the encoded phase.
Two references, 90° apart
The I mixer uses a cosine reference. The Q mixer uses a quarter-cycle-shifted reference. “In phase” and “quadrature” describe alignment with those references.
Fast ripple removed, phase retained
Multiplication creates slow and very fast terms. Low-pass filters discard the fast term, leaving continuous baseband I and Q voltages that change on the echo/offset time scale.
ADC stores two numbers
The ADC samples I and Q. Software writes the pair as S = I + iQ; lowercase i marks the perpendicular Q axis and is not electrical current.
WHAT DOES A MEAN HERE?A is a symbol for the received sample’s amplitude before it is split into I and Q. This teaching plot sets A = 1, so its I and Q numbers are unitless fractions of that reference. In an analog receiver the two outputs are voltages; after ADC they are signed digital counts or calibrated relative units. A does not mean ampere in these equations.
WHAT CHANGED?At +45°, I and Q are equal and positive. The vector rotates, but its magnitude stays 1.000 because only phase changed.
Drag to orbit · violet vectors are B₁⁺ · yellow arrow is net M · coral loops are local elements
α = γ ∫ B₁⁺(t) dt
vRx(t) ∝ −dΦM/dt
The large built-in body coil drives a broad, comparatively uniform B₁⁺ field. During transmit the nearby receive array is actively detuned; during reception the transmitter is isolated and the local elements feed low-noise preamplifiers.
The local sensitivity number uses the exact on-axis field falloff of an ideal 50 mm-radius circular loop, C(d)/C(3 cm) = [(1 + (d/5 cm)²)/(1 + (3/5)²)]−3/2. It is a normalized geometry lesson, not a scanner specification or a full spatial array map. The 3D ruler declares 1 world unit = 100 mm. The shared B₀ control retunes the displayed ¹H carrier; it does not rescale geometry, SNR, relaxation, or SAR. Actual values require coil- and patient-specific electromagnetic measurements.
LIVE 3D LOCAL-COIL RECEIVE JOURNEY · USES THE MODEL ABOVE
What a local coil physically captures—and what happens after.
Drag the stage scrubber with a finger or press Play. The 3D camera follows the signal from the patient to the final combined image.The body coil finishes the high-power B₁⁺ pulse. The transmitter is isolated, the receive-only local loops come out of their protected detuned state, and their low-noise paths are connected.
The programmed pulse ends and the transmit/receive switching network changes electrical connections after a short recovery interval.
Excited transverse magnetization exists in the patient. The local array is now electrically able to respond to it.
Switch timing is measured in seconds, commonly microseconds (µs). No image value has been measured yet.
The receive signal is tiny compared with transmit power. Isolation protects the preamplifier and prevents transmitter leakage from hiding early signal.
The RF pulse did not travel into the local coil as an image. It prepared magnetization; reception is a later electromagnetic induction event.
GdB = 20 log10(Vout/Vin) ↔ Vout = Vin·10GdB/20
Gain scales signal and existing input noise; it does not create localization or improve input SNR.CURRENT SETTINGSAt 40 dB, voltage is multiplied by 100 before digitization. A 128 kHz sample rate gives 7.81 µs between samples in this one-sample-per-dwell model.
Moving transverse magnetization distributed over position r.
One continuous signed RF voltage per element.
Two slower signed components per element.
Digital complex samples assigned to gradient-created addresses.
Aligned coil images become one relative magnitude image; tap for limits.
Transmit field
RF power at the Larmor frequency rotates magnetization. Amplitude and pulse duration set flip angle; spatial B₁⁺ variation makes flip angle nonuniform.
Receive sensitivity
Precessing transverse magnetization induces a tiny voltage. Close local elements couple strongly to nearby anatomy and admit less distant noise.
Isolation matters
Receive-only elements are detuned during the high-power transmit pulse. The transmit path is then isolated while low-noise preamplifiers listen for the echo.
ADVANCED 3D MODEL · FINITE-DURATION RF · FULL BLOCH ROTATION PRODUCT
A slice is not a plane.
It is a solved excitation profile.
A shaped RF pulse and a slice-select gradient act at the same time. Position becomes resonance offset; the pulse spectrum tips a band of those positions; a signed gradient moment then removes most of the phase ramp left by a symmetric pulse. Every vector, curve, and number below comes from the same finite-step Bloch calculation.
The gradient does not cut tissue. It makes resonance frequency vary continuously along the selected axis.
The sinc spectrum predicts the low-flip response; the full Bloch product remains valid when the tip is large.
Magnitude can look acceptable while residual phase still matters for coherent downstream sequences.
The selected band is transverse while off-band magnetization remains near +Mz. Vector direction retains the complex phase that a magnitude-only slice drawing hides.
The quantitative RF waveform and slice-profile plots remain available if 3D rendering is unavailable.
Arrow base = physical position z. Arrow direction = magnetization M.
DRAG TO ORBIT · + / − TO ZOOM · TAP AN OBJECT TO IDENTIFY IT
Mn+1(z) = R[−γ|Beff,n(z)|Δt] Mn(z)
Δf(z) = γ̄Gz + Δf₀ − ΔfRF · BW = TBW / T · Δznom = BW / (γ̄|G|) · Arephase = −GT/2Moving the RF carrier shifts the selected band; increasing |G| spreads the same RF bandwidth over a thinner physical slice.
Apodization reduces sidelobes but broadens the transition. The nominal thickness is not the exact FWHM of every pulse.
The dashed small-tip curve is center-scaled only to compare shape; the RMS metric uses its raw complex amplitude. Its error grows because Mz no longer remains approximately one.
A symmetric pulse accumulates approximately half the selection-gradient area after its effective excitation time; a negative half moment unwinds it.
This is a one-dimensional, single-band, nonadiabatic windowed-sinc pulse. It uses 256 midpoint RF samples and 181 spatial positions with exact norm-preserving rotations. Relaxation, flow, diffusion, chemical shift, multiband excitation, Shinnar–Le Roux pulse design, gradient ramps, concomitant fields, eddy currents, RF amplifier limits, and measured B₀/B₁ maps are not simulated. The post-RF rephaser is represented by its signed moment, so its duration and off-resonance evolution are omitted.