Skip to content
Choose an application

03 / Matter & spectra

One atom. Four ways
to examine the model.

Move from energy levels to probability, then vary the scales and inspect how one proposed level correction affects several transitions.

Established Coulomb benchmarks + declared hypothetical corrections

Spectral transitions

Choose the upper and lower level. Compare the gross Coulomb energy interval, frequency and vacuum wavelength.

Explore gross energy levels and spectral intervals

Select two principal levels. The tool calculates their energy difference in the reduced-mass Coulomb model and converts it into a frequency and a vacuum wavelength. The initial example is the n = 3 to n = 2 Balmer interval.

Choose the interval

The upper level must exceed the lower level. If they cross, the other selector adjusts to the nearest valid level.

Balmer: n = 3 → 2

Energy interval ΔE
1.888651 eV
Frequency ν = ΔE/h
4.566738e+14 Hz
Vacuum wavelength λ = c/ν
656.4696 nm

E₂ = -3.399572 eV; E₃ = -1.510921 eV.

Gross hydrogen energy levels: selected Balmer intervalThe interval from n = 3 to n = 2 is 1.888651 electronvolts. Levels are shown on a linear energy scale above n = 2. Full numerical tables follow.Energy / eVIonization: 0n = 2n = 3ΔE = 1.888651 eV
The vertical axis is a linear energy scale, zoomed to the selected lower level and all higher levels. The zero-energy threshold is ionization. Closely spaced high levels are resolved numerically in the tables.

Balmer interval 3 to 2: 1.888651 eV; 656.4696 nm in vacuum.

Intervals in the selected seriesScroll sideways to see every column.
Upper → lower nΔE / eVν / HzVacuum λ / nm
3 → 2 (selected)1.8886514.566738e+14656.4696
4 → 22.5496796.165096e+14486.2738
5 → 22.855646.904908e+14434.173
6 → 23.0218427.306781e+14410.2935
7 → 23.1220567.549097e+14397.1236
8 → 23.1870997.706370e+14389.019
9 → 23.2316927.814196e+14383.6511
10 → 23.2635897.891323e+14379.9014
11 → 23.2871897.948388e+14377.1739
12 → 23.3051397.991791e+14375.1255

Series limit: 3.399572 eV; 364.7053 nm in vacuum (upper n → ∞).

All principal levels and degeneracies

The n² count is the number of orbital states in the spinless Coulomb model. Multiplying by two includes electron spin before spin-dependent interactions are added. These counts do not assert exact degeneracy after fine, hyperfine, or external-field effects.

Principal-level energies relative to ionizationScroll sideways to see every column.
LevelEnergy / eVOrbital statesIncluding electron spin
n = 1-13.598287312
n = 2 (lower)-3.3995718248
n = 3 (upper)-1.51092081918
n = 4-0.8498929541632
n = 5-0.5439314912550
n = 6-0.3777302023672
n = 7-0.2775160674998
n = 8-0.21247323964128
n = 9-0.1678800981162
n = 10-0.135982873100200
n = 11-0.112382539121242
n = 12-0.0944325504144288

What the interval does and does not specify

An energy difference supplies a photon frequency if a transition with that energy occurs. An electric-dipole transition also needs suitable angular quantum numbers and a nonzero dipole matrix element. For example, the familiar one-photon 2p → 1s line and the two-photon 2s → 1s decay require different transition physics even though their principal levels are the same in this gross model. This tool does not calculate transition probabilities, intensities, or lifetimes.

The displayed wavelengths are in vacuum. Comparisons with visible reference lines must account for any air-wavelength convention and for resolved fine-structure components. See the NIST hydrogen line tables.

μ = memp/(me + mp);   En = −μc²α²/(2n²);   ΔE = Eupper − Elower;   ν = ΔE/h;   λvac = c/ν.

Scope: point nucleus, nonrelativistic Coulomb interaction, reduced mass, and isolated atom. Fine structure, radiative corrections, hyperfine splitting, finite nuclear size, external fields, and recoil terms beyond reduced mass are outside this calculation. Reproducing these formulas with measured constants is a benchmark calculation; it does not derive the Coulomb Hamiltonian from SMM/DTT.

Constants: NIST 2022 CODATA. The calculation uses c = 299792458 m/s, h = 6.62607015 × 10−34 J s, e = 1.602176634 × 10−19 C, α = 7.2973525643 × 10−3, me = 9.1093837139 × 10−31 kg, and mp = 1.67262192595 × 10−27 kg. Displayed digits are rounded model outputs; this tool does not propagate measurement uncertainties.

Try this investigationCompare the 3→2 and 2→1 transitions. Record the approximation before comparing these gross levels with precision spectroscopy.

Where is the electron likely to be found?

Distinguish volume density, radial probability density and cumulative probability. The supplied hydrogenic states are stationary benchmark states.

Explore atomic probability and radius

Compare volume density, radial probability, and the probability contained inside a sphere. The initial 1s example places the radius at one reduced-mass Bohr length, x = r/aμ = 1.

Choose the state and the quantity to display

Radius is in the proton–electron relative coordinate. Use the number field for an exact value. The plotted range changes to 24 aμ when you select an n = 2 state.

Probability inside selected radius
32.3324%
Probability outside selected radius
67.6676%
Selected radius in physical units
52.94654 pm
Radial probability at selected x
0.5413411
Angular-average density × aμ³
0.04307856
1s radial probability and selected radiusThe radial probability peaks at x = 1. The sphere with that radius contains about 32.33 percent of the normalized probability. Numerical values follow the chart.Radial probability per unit x02468x = r / aμ
For the radial view, the shaded area is the probability inside the selected sphere. The volume-density curve has a different integration measure. Every state is normalized over the full range 0 ≤ r < ∞.

1s. The radius 1 aμ contains 32.3324% probability.

The 1s volume density is largest at the origin. Its radial probability is largest at x = 1 because spherical shells have volume proportional to r².

The chart ends at x = 8; 0.001632% of the normalized probability remains beyond that radius.

Probability values in dimensionless radiusScroll sideways to see every column.
x = r/aμaμ³ × average densityP(x), per unit xProbability inside x
00.318309900%
0.50.11709970.36787948.03014%
1 (selected)0.043078560.541341132.3324%
20.0058300490.293050276.1897%
40.00010678110.0214696198.6246%
60.0000019557640.000884766699.9478%
80.000000035821060.00002880999.9984%

Several useful meanings of atomic radius

Mean radius ⟨r⟩
1.5 aμ
Median radius: 50% inside
1.33703 aμ
Radius containing 90%
2.66116 aμ
Root-mean-square radius √⟨r²⟩
1.73205 aμ

A radial probability is a density per unit x; its value at one point is not itself a probability for a finite shell. Integrate it over an interval to obtain that probability. The cumulative curve performs that integration from the origin. The most probable radius, the mean radius, and a probability-containing radius answer different questions.

Normalized functions used by the tool

Write x = r/aμ and let P(x) dx be the probability in a spherical shell. For normalized Coulomb states:

1s: P(x) = 4x²e−2x.
2s: P(x) = x²(2 − x)²e−x/8.
2p: P(x) = x⁴e−x/24.

The angular-average volume density is P(x)/(4πx²), with its continuous limiting value used at the origin. The cumulative probability is F(x) = ∫0xP(u)du. The code evaluates the resulting elementary expressions with a small-radius series to avoid cancellation.

For 1s, F(x) = 1 − e−2x(1 + 2x + 2x²). For 2s, F(x) = 1 − e−x(1 + x + x²/2 + x⁴/8). For 2p, F(x) = 1 − e−x(1 + x + x²/2 + x³/6 + x⁴/24).

The 2p radial function is independent of the magnetic quantum number. Its spatial density is not spherically symmetric for a chosen magnetic substate; the density view deliberately shows the average over directions.

The physical length uses aμ = ħ/(μcα), including reduced mass. These are stationary wavefunction distributions in the effective Coulomb model. The chart does not depict a classical electron orbit, an actual sequence of measurement outcomes, or a simulation of the SMM/DTT inner generative levels.

Constants: NIST 2022 CODATA. The calculation uses c = 299792458 m/s, h = 6.62607015 × 10−34 J s, e = 1.602176634 × 10−19 C, α = 7.2973525643 × 10−3, me = 9.1093837139 × 10−31 kg, and mp = 1.67262192595 × 10−27 kg. Displayed digits are rounded model outputs; this tool does not propagate measurement uncertainties.

Try this investigationCompare the position of the radial peak with the radius containing half of the probability. Explain why they answer different questions.

Change a supplied scale

Inspect how a declared change in coupling or mass ratio changes atomic size and binding energy. A sensitivity calculation does not select the input constant.

Explore linked atomic scales

Change a declared coupling or nuclear mass and see several atomic quantities respond together. Electron mass, Planck’s constant, and the speed of light remain fixed. The initial state uses the reference constants.

Choose the effective model inputs

Increasing α⁻¹ decreases α. The number field retains the full input precision; the slider moves in steps of 0.001.

The logarithmic slider covers R = 100 to 10000. Reference hydrogen uses the proton/electron mass ratio.

Reference inputs

Current coupling α = 0.0072973525643.

Reduced-mass Bohr length aμ
52.94654 pm
Ground-state binding magnitude |E₁|
13.59829 eV
Relative ground-state RMS speed
2187691 m/s
Electron RMS speed in the centre-of-mass frame
2186500 m/s
Atomic scales relative to reference hydrogenAll four ratios initially equal one. Adjusting the inputs changes the radius, binding, and RMS speed ratios together. Full values and percentage changes follow.Ratio to reference hydrogenRadius1Binding1Relative speed1Electron speed11 = reference
Every bar is divided by its own reference hydrogen value. The dashed line marks a ratio of one. The horizontal scale adapts so that all current values remain visible.

Reference inputs. Radius 52.94654 pm; binding energy 13.59829 eV.

Current values compared with reference hydrogenScroll sideways to see every column.
QuantityCurrentReferenceChange
Radius aμ / pm52.9465409552.946540950%
Binding |E₁| / eV13.5982872613.598287260%
Relative RMS speed / c0.0072973525640.0072973525640%
Electron COM RMS speed / c0.0072933804650.0072933804650%
Reduced-mass fraction μ / mₑ0.99945567940.99945567940%
2 → 1 vacuum interval / nm121.5684456121.56844560%

Read the response as a linked calculation

μ/me = R/(1 + R);   aμ = ħ/(μcα);   |E₁| = μc²α²/2;
vrel,rms/c = α;   ve,COM,rms/c = (μ/me)α.

At fixed mass ratio, a 1% increase in α produces approximately a 1% decrease in radius, a 2% increase in binding, and a 1% increase in either RMS speed. These are local fractional sensitivities; the full nonlinear formulas determine the displayed finite changes.

The local mass sensitivity is q = 1/(1 + R) = 0.0005443206. A 1% increase in R changes the binding by approximately q%, the radius by −q%, and the electron COM RMS speed by q%. The relative RMS speed does not change when only R changes.

Why there are two RMS speeds

The Coulomb wavefunction describes relative separation. Its momentum gives a relative RMS speed αc in the ground state. In the centre-of-mass frame, the electron’s momentum is the relative momentum, but its velocity is obtained by dividing by me, so ve,COM,rms = (μ/me)vrel,rms. This factor approaches one for an infinitely heavy nucleus.

RMS speed measures momentum variance. It does not imply a classical circular trajectory. In particular, a real stationary 1s wavefunction can have zero probability current and nonzero momentum variance.

Meaning of the 137 control: the button enters α = 1/137 into an effective hydrogen calculation. It displays the consequences of that chosen input. Identifying a structural count with a physical coupling requires an independent matching derivation and any physical correction; moving this control does not supply that derivation. The reference mode uses measured constants, and all changed settings are marked illustrative.

Constants: NIST 2022 CODATA. The calculation uses c = 299792458 m/s, h = 6.62607015 × 10−34 J s, e = 1.602176634 × 10−19 C, α = 7.2973525643 × 10−3, me = 9.1093837139 × 10−31 kg, and mp = 1.67262192595 × 10−27 kg. Displayed digits are rounded model outputs; this tool does not propagate measurement uncertainties.

Try this investigationChange the coupling slightly, then compare the fractional changes in radius and binding energy.

One correction, several spectral consequences

Apply one hypothetical level-shift rule consistently across transitions. This is a correlated model prediction, not an experimental residual dataset.

Explore shared level shifts and differential intervals

Synthetic illustration. Assign one correction to each state, then obtain every interval correction by subtracting the two state corrections. This shows why independent per-line adjustments can fail and why a shared response pattern may still be physically ambiguous.

δEnℓ/h = C + (A + B)δℓ0/n³, in kHz.
δνa→b = (δEb − δEa)/h.

The contact profile acts on S states (ℓ = 0). P and D states have zero contact weight in this simplified model. A and B are separately adjustable labels for exactly the same profile. C is a common offset. No measured spectral residuals are loaded.

Choose illustrative coefficients

A is an illustrative coefficient to be determined by a specified model. It is not a calculated DTT value.

B deliberately has the same state dependence as A, illustrating an indistinguishable competing contribution.

Changing a common energy reference shifts all levels equally and cancels from every interval.

The transfer control keeps A + B fixed. Close to a coefficient limit it reverses the transfer direction.

Total contact coefficient A + B
+10 kHz
Signed correction: 2S → 8D
−1.25 kHz
Synthetic interval corrections from shared level shiftsAt A = 10 kHz and B = C = 0, the 2S to 8D interval correction is minus 1.25 kHz because the lower 2S state shifts while the upper 8D contact contribution is zero. All values follow in tables.Synthetic interval correction / kHz1S → 2S2S → 8S1S → 8S2S → 8D2S → 2P2P → 8DNegative ← signed shift → Positive
Bars show signed corrections to specified energy intervals. A positive correction increases Eb − Ea. No claim about a single-photon transition, intensity, or experimental detectability follows from a bar.

Synthetic model. Total contact coefficient +10 kHz. The 2S to 8D correction is −1.25 kHz. The shared-level cycle closes.

Signed interval corrections from the same state modelScroll sideways to see every column.
Interval a → bWeight differenceA contribution / kHzB contribution / kHzTotal / kHz
1S → 2S-0.875-8.75-0-8.75
2S → 8S-0.123046875-1.23046875-0-1.23046875
1S → 8S-0.998046875-9.98046875-0-9.98046875
2S → 8D-0.125-1.25-0-1.25
2S → 2P-0.125-1.25-0-1.25
2P → 8D0000
State corrections used in every interval
Common state-level inputScroll sideways to see every column.
StateContact weightA / kHzB / kHzC / kHzδE/h / kHz
1S1+1000+10
2S0.125+1.2500+1.25
2P00000
3S0.037037037+0.3703703700+0.37037037
8S0.001953125+0.0195312500+0.01953125
8D00000

1. Verify the cycle relation

For three levels a, b, and c, a shared energy model gives δνa→c = δνa→b + δνb→c. The cancellation is algebraic, independent of the size of A, B, or C.

-8.75 + (-1.23046875) = -9.98046875 kHz
Cycle residual: 0 kHz (to numerical rounding).

2. Keep the lower S-state correction

For 2S → 8D, the upper D state has zero contact weight. The lower 2S state has weight 1/8, so the interval correction is −(A + B)/8. The contact contribution therefore remains present. Only an interval with equal contact weights cancels this profile.

3. Check whether the mechanism can be identified

Every interval depends on A and B through their sum. Moving 5 kHz from one coefficient to the other leaves all interval predictions unchanged. Even arbitrarily precise measurements of these same intervals could not separate A from B without additional physical information or observables with a different response pattern.

This is an illustration of a possible degeneracy with another contact contribution, such as a leading S-state finite-size term in an appropriate effective model. A real analysis must derive the coefficients, their dependence on mass and nuclear properties, their uncertainties, and the applicable perturbation operators. The equality of the two profiles here is deliberately stipulated to make the identifiability issue visible.

4. Specify the next theoretical calculation

A physical DTT correction needs a defined operator or state-energy map, a fixed normalization, and a coefficient derived from the source construction. A holonomy or structural quantity cannot be converted to kilohertz by selecting an arbitrary scale. Once that calculation is available, the same level-shift machinery can propagate it consistently across many intervals.

This tool contains synthetic coefficients and no measured residual dataset. Background and derivation requirements: differential comparisons, research status, and sources. The common formula is deliberately elementary so that every displayed value can be checked by hand.

Try this investigationUse one set of level-shift parameters for every transition. Identify which ratios or differences would provide an independent comparison.

From a benchmark to a test

The atomic Hamiltonian and reference constants are supplied. The SMM/DTT task is to derive a physical state space, interaction and matching law before claiming a new spectral prediction.

The gross-level tools exclude the detailed fine structure, Lamb shift and other corrections needed for precision hydrogen work. Use the published atomic references and the site’s differential-test guide when choosing a comparator.

SMONAD

Share this page

Share this page with its own title, summary and preview image.

More platforms