Pith. sign in

REVIEW 5 major objections 4 minor 1 cited by

The Bloch Equation Generator -- SimuF\'isica

T0 review · 5 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The Bloch Equation Generator derives and numerically solves the complete coupled Bloch equations for any configured N-level system with up to 30 levels, returning compilable C code.

desk verdict A genuine and useful tool, decently validated by three examples, but the printed derivation contradicts the shipped code and the generality claim is unproven. read the letter →

arxiv 2506.01108 v1 pith:OPN7Y4GE submitted 2025-06-01 quant-ph physics.atom-phphysics.comp-phphysics.optics

M. P. M. de Souza , G. H. H. Pavão , A. A. C. de Almeida , S. S. Vianna This is my paper
classification quant-phphysics.atom-phphysics.comp-phphysics.optics
keywords opticalBlochequationsmultilevelquantumsystemsdensitymatrixrotatingwaveapproximationcoherentpopulationtrappingelectromagneticallyinducedtransparencyrubidium-87Zeemansublevelsnumericalsimulationsoftware
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper presents the Bloch Equation Generator (BEG), a free browser-based tool that turns a user-drawn level diagram into the full set of optical Bloch equations and into C source code that integrates them. The claim is that for any system a user can configure — up to 30 levels, with chosen transitions, decay rates, and an optional rotating-wave approximation — the tool automatically produces the $N(N+1)/2$ independent equations and solves them both in the time domain and as a function of field detuning. The authors illustrate the claim with a two-level system, a $\Lambda$-type system showing coherent population trapping, electromagnetically induced transparency, and Autler–Townes splitting, and a realistic 12-level Zeeman-resolved $^{87}$Rb system with $\sigma$ and $\pi$ transitions. If the automated derivation is correct across configurations, the tool removes a laborious and error-prone manual step that currently limits multilevel Bloch-equation work.

What carries the argument

The central object is the generator itself: an interactive level-diagram editor that stores levels, allowed dipole transitions, decay pathways, and Rabi frequencies, then applies the slow-variable change $\rho_{ij}=\sigma_{ij}e^{i\omega_m t}$ for off-diagonal elements and drops $2\omega_m$ terms when the rotating-wave approximation is selected. From these rules it emits the independent density-matrix equations, treating $\sigma_{ij}$ and $\sigma_{ji}$ as conjugates, and produces C code whose numerical engine is fourth-order Runge–Kutta. Equation (9), $\gamma_{ij}=1/T_{ii}+1/T_{jj}$, fixes coherence relaxation from the configured decay rates, and the chosen field-mode frequencies fix each detuning $\delta_{ij}=\omega_m-\omega_{ij}$; the generator's ability to thread those detunings through multi-field configurations is what the paper asks users to verify manually.

What would settle it

Take a four-level diamond or other multi-field configuration with nonzero detunings, run the BEG, and compare the displayed equations and generated C code term-by-term against an independent hand-derived or symbolic derivation; any error in a detuning sum, a coherence coupling, or a decay bookkeeping term would falsify the general claim. A cheaper check is to run the emitted code on random parameters and require the population sum and density-matrix trace to stay at machine precision for the whole integration.

Watch

Extended reading notes

Core claim

The central claim is that the Bloch Equation Generator replaces the manual construction of optical Bloch equations for multilevel systems. For $N$ levels it generates $N(N+1)/2$ independent coupled differential equations consistent with Liouville–von Neumann density-matrix evolution, phenomenological decay terms, and the user's rotating-wave approximation choice, and it packages the equations as compilable C code using fourth-order Runge–Kutta integration. Time-domain evolution and steady-state spectra versus detuning are both supported. The three worked configurations reproduce known physics: Rabi oscillations and a power-broadened Lorentzian spectrum in two levels; the dark-state coherence $\rho_{13} \to -\Omega_{12}\Omega_{23}^*/(\Omega_{12}^2+\Omega_{23}^2)$, the EIT transparency dip, and the Autler–Townes doublet in the $\Lambda$ system; and, in the 12-level $^{87}$Rb case, population funneling into stretched Zeeman states for $\sigma^+$ excitation that matches an effective two-level description at strong Rabi frequency.

Load-bearing premise

The load-bearing premise is that the automatically generated equations are correct for every allowed configuration, since the paper validates only three hand-built cases and warns users to check the detuning sums manually.

Editorial extensions

If this is right

  • A user can go from a drawn level scheme to a working numerical solver without handwriting equations, for up to 30 levels (up to 465 coupled equations).
  • The $\Lambda$-system example shows that CPT, EIT, and Autler–Townes splitting can be obtained within one configurable tool, covering several nonlinear-optics regimes.
  • The 12-level $^{87}$Rb example shows the tool can handle realistic Zeeman-resolved transitions with 78 coupled equations and polarization-dependent selection rules.
  • Generated code supports both temporal evolution and detuning-domain spectra, so steady-state line shapes and power broadening follow from the same configuration.
  • Because the equations are displayed before code generation, users can audit the automated derivation before running a numerical experiment.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • General correctness would be better established by comparing generated equations against an independent symbolic derivation for a family of level schemes rather than only hand-built examples; the paper's own warning about $\delta_{31}=\delta_{21}+\delta_{32}$ suggests detuning bookkeeping is the likeliest failure mode.
  • The same architecture could plausibly extend to in-browser steady-state spectra and GPU acceleration, which the paper lists as future work; a demanding test would be a scheme with more than one field mode per transition, where the single-mode slow-variable ansatz in Eq. (6) is strained.
  • If the general claim holds, this style of automated equation generation could serve as a cross-check for hand-derived Bloch equations in EIT, laser-cooling, and quantum-memory models, where 12–30 level schemes are common.
  • An inexpensive additional validation would be running the emitted code on many random parameter sets and checking trace and population-sum conservation to machine precision, since the paper reports only a two-level accuracy proxy.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

5 major / 4 minor

Summary. The manuscript presents the Bloch Equation Generator (BEG), a free browser-based tool on the SimuFísica platform that automatically constructs and numerically integrates optical Bloch equations for N-level systems with 2≤N≤30. Users configure levels, allowed transitions, decay rates, Rabi frequencies, and the rotating-wave approximation; the tool displays the full set of N(N+1)/2 equations and provides downloadable C source code for time-domain and detuning-swept (frequency-domain) solutions. The paper validates the tool on three examples: a two-level system showing Rabi oscillations and power-broadened spectra, a Λ-type three-level system showing CPT, EIT, and Autler–Townes splitting, and a 12-level Zeeman-resolved model of the 87Rb 5S1/2(F=2)→5P3/2(F=3) transition for σ+ and π polarizations. The authors argue that the automation is important because manual formulation of Bloch equations becomes laborious and error-prone for large N.

Significance. If the tool works as claimed, this is a useful community resource: it is freely accessible, requires no local installation, and it ships reproducible example configurations and source-code links. The three benchmarks are well chosen and reproduce known physics: the two-level population spectrum matches the expected FWHM to excellent accuracy, the Λ-system dark-state coherence approaches the analytical value −0.5, and the 12-level σ+ example shows optical pumping into the extreme mF states. These successes give concrete evidence that the implemented numerical code solves physically meaningful equations for the tested cases. However, the printed theoretical derivation in Section 2 is inconsistent with the generated code in several places, so the paper as written does not establish that the automated equations produced by BEG are correct for arbitrary user configurations or for the full claimed parameter space.

major comments (5)
  1. [Section 2, Eq. (7b)] For N=2, Eq. (7b) reduces to dσ12/dt = iΩ12σ22 + (iδ−γ)σ12, with σ22=ρ22. This equation has no population-difference source term of the form proportional to (ρ22−ρ11), so an initial state with ρ11=1 would never develop coherence. Since Fig. 4(b) shows Rabi oscillations starting from that initial condition, the implemented equation cannot be Eq. (7b). The paper must either correct Eq. (7b) to include the population-difference driving term or explicitly state that the printed derivation is only schematic and not the exact specification of the generated code.
  2. [Section 2, Eq. (9)] Equation (9) gives γij = 1/Tii + 1/Tjj, which for the two-level case with T22=1/Γ21 and T11=∞ predicts γ12=Γ21. However, Listing 1 sets γ12=0.5*decay (i.e., Γ21/2), and Listing 3 sets γ12=γ23=0.5*decay. The generated code therefore uses γij=(Γi+Γj)/2, differing from Eq. (9) by a factor of two. The printed formula and the code must be reconciled, or the convention must be stated explicitly.
  3. [Section 2, Eq. (8c)] The decay term in Eq. (8c) is written as +Γjlρll(t), but the population of state |l⟩ should be fed by the decay from |j⟩, so the term should read +Γjlρjj(t). As printed, the equation describes a decay of ρll rather than an increase. This error, together with the issue in Eq. (7b), means the displayed equations are not a reliable reference for what the generator actually computes.
  4. [Section 3 and Section 2, Eq. (6)] The abstract and Section 3 claim that BEG handles systems under coherent excitation involving multiple frequencies, but Eq. (6) transforms each off-diagonal density-matrix element with a single frequency ωm, and the derivation does not cover the case of several field modes driving the same transition. The paper never states this restriction as a limitation, and none of the three examples tests multi-mode-per-transition configurations. The generality claim for arbitrary N≤30 systems is therefore not supported by the derivation or the validation.
  5. [Section 4.2, Listing 4] The text explicitly says that for the Λ-system spectrum the user must verify whether detunings such as δ31=δ21+δ32 are correctly written (see line 16 of Listing 4). This is a direct admission that the generator does not automatically enforce multi-photon detuning consistency. The paper should state this limitation prominently and, ideally, describe an automated consistency check in the tool, since manual verification is precisely the kind of error-prone step the tool claims to eliminate.
minor comments (4)
  1. [Section 4.1, Fig. 4(d)] The caption calls the plotted quantity 'absolute error' but the axis label shows ρ11+ρ22−1, which is the signed deviation from unity. Please clarify the terminology.
  2. [Section 4.1, FWHM expression] The expression FWHM = 2√(γ12² + 4Ω12²(γ12/Γ21)) is nonstandard in form; please define the power-broadening convention and give a textbook reference in addition to the thesis [12].
  3. [Listing 5] The branching factors for spontaneous decay from the excited F=3 sublevels (e.g., Gamma61=decay, Gamma71=decay/2, Gamma81=decay/3) appear ad hoc; please specify whether these come from Clebsch–Gordan coefficients or represent a simplified equal-branching assumption.
  4. [References] Reference [14] contains a malformed 'Preprint' field with a duplicated URL; please clean up the entry.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular derivation: the Bloch-equation generator is validated against independent analytic and literature results; the only self-citation, for the two-level FWHM, is not load-bearing.

full rationale

The paper's derivation chain is standard and self-contained: Section 2 obtains the optical Bloch equations from the Liouville-von Neumann equation with explicit dipole and RWA assumptions, and Section 3 describes the generator as an automated implementation of those equations. The three validation examples benchmark the generated C code against independent results: (i) the two-level Rabi oscillations and power-broadened spectrum are compared with the analytic FWHM = 2 sqrt(gamma12^2 + 4 Omega12^2 (gamma12/Gamma21)) expression cited to reference [12]; (ii) the Lambda system is checked against the analytic dark-state coherence rho13 -> -Omega12 Omega23*/(Omega12^2+Omega23^2) and the EIT/ATS literature; (iii) the 12-level Rb system is compared with the known two-level reduction for sigma polarization from reference [19]. In every case the expected values are computed from the same physical parameters but are not fitted to the tool's output, and the equations are not defined in terms of the quantities they are used to predict. The only self-citation is the FWHM formula from the first author's PhD thesis [12]; this is a standard analytic benchmark, independent of the BEG itself, and it does not support any uniqueness claim or restrict the derivation, so it is not load-bearing. Concerns raised by the manuscript's own caveat in Section 4.2, that users must manually verify two-photon detunings such as delta31 = delta21 + delta32, are validation-completeness issues rather than circularity: they show the generator's general correctness for arbitrary N<=30 configurations is assumed rather than exhaustively demonstrated, but they do not make any prediction equivalent to its input by construction. No fitted parameter is renamed as a prediction and no derivation step reduces to a self-citation chain. Hence the circularity score is low, reflecting only the minor non-load-bearing self-citation.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The paper makes no fitted-parameter claims; the examples use hand-chosen rates and Rabi frequencies, which do not affect the validity of the tool itself. No new physical entities are postulated.

free parameters (2)
  • Gamma21 (two-level example) = 2*pi*5 MHz
    Set in Section 4.1 to approximate the 87Rb D2 linewidth; an input choice for illustration, not fitted to data.
  • Rabi frequencies Omega_ij = Omega12 = Gamma21; Omega12 = Omega23 = 0.1 Gamma21; Omega_ij = Gamma
    Chosen by hand in each example to access different physical regimes; not fitted to data.
assumptions (5)
  • standard math Liouville-von Neumann equation (Eq. 1)
    Invoked as the starting point for density matrix evolution.
  • domain assumption Electric dipole approximation (Eq. 3)
    Assumes the atom-field interaction is dominated by the electric dipole term.
  • domain assumption Rotating-wave approximation
    Neglects rapidly oscillating 2*omega_m terms, valid for near-resonant driving.
  • domain assumption Phenomenological decay terms
    Population and coherence relaxation introduced by hand via rates Gamma_ij and gamma_ij.
  • domain assumption Non-interacting ensemble
    Equation (1) describes a single system; ensemble averaging is implicit.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The Bloch Equation Generator -- SimuF\'isica." pith.science (2026). https://pith.science/paper/OPN7Y4GE

@misc{pith2026250601108,
  author       = {Pith},
  title        = {Pith review of: The Bloch Equation Generator -- SimuF\'isica},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OPN7Y4GE}},
  note         = {Machine review of arXiv:2506.01108}
}
abstract

The interaction between multilevel quantum systems and coherent radiation underlies several phenomena in modern atomic optics. The formulation and solution of the Bloch equations, which describe the dynamics of such systems, become complex as the number of levels increases. In this work, we present the Bloch Equation Generator, a free, browser-based computational tool developed to automate the generation and numerical solution of Bloch equations for systems with up to 30 levels. Users can configure the level diagram, select allowed transitions, define decay rates, and choose whether or not to apply the rotating wave approximation. The software automatically generates the complete set of equations and provides C source code for numerical solutions in both the time and frequency domains. To illustrate its applicability, we present three examples: (i) a two-level system, (ii) a $\Lambda$-type system with analysis of CPT, EIT, and the Autler-Townes effect, and (iii) a realistic 12-level system based on the Zeeman-resolved $5S_{1/2} \to 5P_{3/2}$ transition of rubidium-87.

Figures

Figures reproduced from arXiv: 2506.01108 by the authors.

Figure 1
Figure 1. Bloch Equation Generator, version 1.4.2 (18/04/2025), configured for a four-level diamond [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Bloch equations generated for the configuration shown in Fig. 1. [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Code download panel in the BEG interface. [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (6 more)
Figure 4
Figure 4. Figure 4: (a) Two-level system. (b) Time evolution of populations. (c) Steady-state populations as a [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: (a) Λ-type three-level system. (b) Time evolution of the populations. (c) Time evo [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: Absorption of the weak probe field (Ω12/Γ21 = 0.1) in a Λ system as a function of δ21, for two values of Ω32. Source code: https://github.com/marcopolomoreno/bloch-equation-generator/ blob/main/three-level-%CE%9B/detuning.c. 10 double Gamma21 = decay /2.0; 11 double Ga…
Figure 7
Figure 7. Figure 7: (a) σ + transitions between Zeeman sublevels in the 5S1/2, F = 2 → 5P3/2, F = 3 line of 87Rb. (b) Time evolution of the ground state populations. (c) Time evolution of the excited state populations. Link to this configuration: https://simufisica.com/PPT3E. Source code:…
Figure 8
Figure 8. Figure 8: Comparison between the excited-state population in a two-level system ( [PITH_FULL_IMAGE:figures/full_fig_p009_8.png]
Figure 9
Figure 9. Figure 9: (a) π transitions between Zeeman sublevels in the 5S1/2, F = 2 → 5P3/2, F = 3 line of 87Rb. (b) Time evolution of the ground-state Zeeman populations. (c) Time evolution of the excited-state Zeeman populations. Link to this configuration: https://simufisica.com/ nwOow.…

Discussion (0). Sign in to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Exploring Heat Exchanges with the Calorimetry Simulator -- SimuF\'isica

    physics.ed-ph 2025-06 conditional novelty 4.0 of 10

    The Calorimetry SimuFísica simulator reproduces textbook calorimetry calculations in two worked examples, but the paper's stated thermal equilibrium equations omit masses and are inconsistent with the reported results.

Reference graph

Works this paper leans on

21 extracted references · 15 canonical work pages · cited by 1 Pith paper

  1. [1]

    KGaA) ISBN 9783527617197 URL https://onlinelibrary.wiley.com/doi/book/10.1002/9783527617197

    Cohen-Tannoudji C, Dupont-Roc J and Grynberg G 1998Atom—Photon Interactions: Basic Pro- cess and Appilcations(WILEY-VCH Verlag GmbH & Co. KGaA) ISBN 9783527617197 URL https://onlinelibrary.wiley.com/doi/book/10.1002/9783527617197

    doi:10.1002/9783527617197

  2. [2]

    Allen L and Eberly J 2012Optical Resonance and Two-Level AtomsDover Books on Physics (Dover Publications) ISBN 9780486136172 URLhttps://books.google.com.br/books? id=36DDAgAAQBAJ

    work page

  3. [3]

    Zigdon T, Wilson-Gordon A D and Friedmann H 2009Phys. Rev. A80(3) 033825 URLhttps: //link.aps.org/doi/10.1103/PhysRevA.80.033825

    work page doi:10.1103/physreva.80.033825

  4. [4]

    Express296495–6508 URLhttps://opg.optica.org/oe/ abstract.cfm?URI=oe-29-5-6495

    Noh H R and Moon H S 2021Opt. Express296495–6508 URLhttps://opg.optica.org/oe/ abstract.cfm?URI=oe-29-5-6495

    work page

  5. [5]

    Sakurai J and Napolitano J 2017Modern Quantum Mechanics(Cambridge University Press) ISBN 9781108422413 URLhttps://books.google.com.br/books?id=010yDwAAQBAJ 10

    work page

  6. [6]

    Souza M P M, Oliveira S P and Luiz V 2024arXiv:2410.18721 [physics.ed-ph]This is the English version of a peer-reviewed paper published in Portuguese: Rev. Bras. Ensino F ´ ıs. 46, e20230219 (2024). DOI:https://doi.org/10.1590/1806-9126-RBEF-2023-0219URLhttps://doi.org/ 10.48550/arXiv.2410.18721

    arXiv doi:10.48550/arxiv.2410.18721 2024

  7. [7]

    DOI:https://doi.org/10.59727/fne.v22i1.173URLhttps://doi.org/10

    Souza M P M, Oliveira C M and Ara´ ujo R P P 2024arXiv:2410.17951 [physics.ed-ph]This is the English version of a peer-reviewed paper published in Portuguese: A F ´ ısica na Escola 22, 240173 (2024). DOI:https://doi.org/10.59727/fne.v22i1.173URLhttps://doi.org/10. 48550/arXiv.2410.17951

    arXiv doi:10.48550/arxiv.2410.17951 2024

  8. [8]

    Press W H 2007Numerical Recipes in C: The Art of Scientific Computing(Cambridge University Press)

    work page

Show all 21 references
  1. [9]

    Ara´ ujo M O, Marinho L S and Felinto D 2022Phys. Rev. Lett.128(8) 083601 URLhttps://link. aps.org/doi/10.1103/PhysRevLett.128.083601

    doi

  2. [10]

    Jacques V, Hingant B, Allafort A, Pigeard M and Roch J F 2009European Journal of Physics30 921 URLhttps://dx.doi.org/10.1088/0143-0807/30/5/001

    doi

  3. [11]

    Steck D A Rubidium 87 d line data (revision 2.3.3, 28 may 2024) URLhttp://steck.us/ alkalidata

    2024

  4. [12]

    thesis Universidade Federal de Pernambuco

    Souza M P M 2012Excita¸ c˜ ao coerente de um vapor atˆ omico por trens de pulsos ultracurtos e lasers cont ´ ınuosPh.D. thesis Universidade Federal de Pernambuco

  5. [13]

    Arimondo E 1996 Coherent population trapping in laser spectroscopy (Progress in Opticsvol 35) ed Wolf E (Elsevier) pp 257–354 URLhttps://www.sciencedirect.com/science/article/pii/ S0079663808705316

    1996

  6. [14]

    Harris S E 1997Physics Today5036–42 ISSN 0031-9228 (Preprinthttp://arxiv.org/abs/https://pubs.aip.org/physicstoday/article- pdf/50/7/36/8312202/36 1 online.pdf:https://pubs.aip.org/physicstoday/article- pdf/50/7/36/8312202/36 1 online.pdf) URLhttps://doi.org/10.1063/1.881806

    doi

  7. [15]

    Fleischhauer M, Imamoglu A and Marangos J P 2005Rev. Mod. Phys.77(2) 633–673 URLhttps: //link.aps.org/doi/10.1103/RevModPhys.77.633

    doi

  8. [16]

    Hao L, Jiao Y, Xue Y, Han X, Zhao S B J and Raithel G 2018New Journal of Physics20073024 URLhttps://iopscience.iop.org/article/10.1088/1367-2630/aad153

    doi

  9. [17]

    Jadoon Z, Noh H and Kim J 2022Sci Rep12145 URLhttps://www.nature.com/articles/ s41598-021-04038-5

  10. [18]

    Foot C J 2005Atomic Physics(Oxford University Press)

  11. [19]

    Lezama A, Cardoso G C and Tabosa J W R 2000Phys. Rev. A63(1) 013805 URLhttps://link. aps.org/doi/10.1103/PhysRevA.63.013805

    doi

  12. [20]

    Zigdon T, Wilson-Gordon A D and Friedmann H 2008Phys. Rev. A77(3) 033836 URLhttps: //link.aps.org/doi/10.1103/PhysRevA.77.033836

    doi

  13. [21]

    Demeter G 2013Computer Physics Communications1841203–1210 ISSN 0010-4655 URLhttps: //www.sciencedirect.com/science/article/pii/S0010465512004250 11

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.