qojulia · ChristophHotter · Aug 20, 2026 · Aug 20, 2026
diff --git a/benchmarks/correlations.jl b/benchmarks/correlations.jl
@@ -23,7 +23,7 @@ function benchmark_correlations!(SUITE)
     G_cas = ▷(G_u, G_c, G_v)
 
     H_sym = hamiltonian(G_cas)
-    L_sym = lindblad(G_cas)[1]
+    L_sym = jump_operator(G_cas)[1]
 
     γ_ = 1.0
     σ_pulse = 1 / γ_

diff --git a/benchmarks/interaction_picture.jl b/benchmarks/interaction_picture.jl
@@ -67,7 +67,7 @@ function benchmark_interaction_picture!(SUITE)
     G_cas = ▷(G_u, G_s, G_v)
 
     H = hamiltonian(G_cas)
-    L = lindblad(G_cas)[1]
+    L = jump_operator(G_cas)[1]
     H_uv = hamiltonian(▷(G_u, G_v))
     H_int_ = simplify(H - H_uv)
 

diff --git a/benchmarks/slh_algebra.jl b/benchmarks/slh_algebra.jl
@@ -26,7 +26,7 @@ function benchmark_slh_algebra!(SUITE)
     SUITE["SLH Algebra"]["symbolic"]["3-cavity cascade"] = @benchmarkable begin
         G_cas = ▷($G_u, $G_c, $G_v)
         hamiltonian(G_cas)
-        lindblad(G_cas)
+        jump_operator(G_cas)
     end
 
     SUITE["SLH Algebra"]["symbolic"]["concatenate + cascade"] = @benchmarkable begin
@@ -87,7 +87,7 @@ function benchmark_slh_algebra!(SUITE)
         G_L_t = $G_L2 ▷ $G_ϕ_12 ▷ $G_L1
         G_t = G_R_t ⊞ G_L_t
         hamiltonian(G_t)
-        lindblad(G_t)
+        jump_operator(G_t)
     end
 
     ## --- Time-dependent closure evaluation (the ODE hot loop) ---
@@ -99,7 +99,7 @@ function benchmark_slh_algebra!(SUITE)
     G_t = G_R_t ⊞ G_L_t
 
     H_f = hamiltonian(G_t)
-    L_f = lindblad(G_t)
+    L_f = jump_operator(G_t)
 
     t_mid = T[length(T)÷2]
 

diff --git a/benchmarks/translation.jl b/benchmarks/translation.jl
@@ -30,7 +30,7 @@ function benchmark_translation!(SUITE)
     G_v = SLH(1, gv_sym * av_, 0)
     G_cas = ▷(G_u, G_c, G_v)
     H_sym = hamiltonian(G_cas)
-    L_sym = lindblad(G_cas)[1]
+    L_sym = jump_operator(G_cas)[1]
 
     ## --- Static translation (no time dependence) ---
 

diff --git a/docs/src/api.md b/docs/src/api.md
@@ -11,6 +11,7 @@ scattering
 ```
 
 ```@docs
+jump_operator
 lindblad
 ```
 

diff --git a/docs/src/implementation.md b/docs/src/implementation.md
@@ -21,13 +21,13 @@ av = Destroy(h, :a_v, 3)
 
 An SLH component is represented as `(S, L, H)` by the [`SLH`](@ref) type. The cascade [`▷`](@ref), concatenation [`⊞`](@ref), and feedback reduction [`feedback`](@ref) rules implement the standard network composition from the SLH framework.
 
-The resulting effective operators are accessed by [`hamiltonian`](@ref) and [`lindblad`](@ref) and remain symbolic until translation. This is especially useful when you want to further manipulate the expressions, e.g. to transform into the interaction picture. 
+The resulting effective operators are accessed by [`hamiltonian`](@ref) and [`jump_operator`](@ref) and remain symbolic until translation. This is especially useful when you want to further manipulate the expressions, e.g. to transform into the interaction picture.
 
 ```julia
 G_cas = ▷(G_u, G_s, G_v)
 
 H = hamiltonian(G_cas)
-L = lindblad(G_cas)
+L = jump_operator(G_cas)
 ```
 
 For networks with internal loops, the symbolic model can be reduced directly with [`feedback`](@ref), which applies the SLH feedback reduction rule before translation. This keeps the symbolic workflow consistent: build a network from cascades and concatenations, eliminate internal connections symbolically, and only then translate the reduced Hamiltonian and Lindblad operators to numerics.

diff --git a/docs/src/tutorial.md b/docs/src/tutorial.md
@@ -52,7 +52,7 @@ H = hamiltonian(G_cas)
 ```
 
 ```@example tutorial
-L = lindblad(G_cas)[1]
+L = jump_operator(G_cas)[1]
 ```
 
 ## 2. Numerical parameters and input pulse

diff --git a/examples/01-1_cavity-scattering__PRL2019_123-123604_fig2-fig3.jl b/examples/01-1_cavity-scattering__PRL2019_123-123604_fig2-fig3.jl
@@ -44,7 +44,7 @@ H = hamiltonian(G_cas)
 
 #
 
-L = lindblad(G_cas)[1] # only one Lindblad term in this example
+L = jump_operator(G_cas)[1] # only one Lindblad term in this example
 
 # To solve the dynamics of the system we translate the symbolic expressions into numeric operators (matrices) of QuantumOptics.jl. To do so, we define the numerical parameters and operator basis.
 

diff --git a/examples/01-2_stimulated-emission__PRL2019_123-123604_fig4.jl b/examples/01-2_stimulated-emission__PRL2019_123-123604_fig4.jl
@@ -43,7 +43,7 @@ H = G_cas.hamiltonian
 
 #
 
-L = G_cas.lindblad[1] # only one Lindblad in this example
+L = G_cas.jump_operator[1] # only one Lindblad in this example
 
 #
 

diff --git a/examples/02-1_cavity-phase-noise__PRA2020_102- 023717_fig2.jl b/examples/02-1_cavity-phase-noise__PRA2020_102- 023717_fig2.jl
@@ -41,7 +41,7 @@ H = hamiltonian(G_cas)
 
 #
 
-L = lindblad(G_cas)[1] # only one Lindblad in this example
+L = jump_operator(G_cas)[1] # only one jump operator in this example
 
 #
 

diff --git a/examples/02-3_mode-entanglement__PRA2020_102-023717_fig4.jl b/examples/02-3_mode-entanglement__PRA2020_102-023717_fig4.jl
@@ -48,7 +48,7 @@ G_v2 = SLH(1, g_v2' * av2, 0)
 G = cascade(G_s, G_v1, G_v2)
 
 H = hamiltonian(G)
-L = lindblad(G)[1]
+L = jump_operator(G)[1]
 nothing # hide
 
 # We use the parameters quoted in the paper and initialize the emitter in the excited state.

diff --git a/examples/03-1_beam-combiner__PRA2023_107-023715_fig2-fig3.jl b/examples/03-1_beam-combiner__PRA2023_107-023715_fig2-fig3.jl
@@ -49,7 +49,7 @@ H = hamiltonian(G_cas)
 
 #
 
-L = lindblad(G_cas)[1] # only one Lindblad in this example
+L = jump_operator(G_cas)[1] # only one Lindblad in this example
 
 # Next, the numerical parameters and functions of the system are defined.
 

diff --git a/examples/04-1_two-sided-cavity_with-atom_coh-drive.jl b/examples/04-1_two-sided-cavity_with-atom_coh-drive.jl
@@ -46,11 +46,11 @@ H1 = hamiltonian(G_cav_L_R_drive)
 
 # 
 
-L1_L = lindblad(G_cav_L_R_drive)[1]
+L1_L = jump_operator(G_cav_L_R_drive)[1]
 
 # 
 
-L1_R = lindblad(G_cav_L_R_drive)[2]
+L1_R = jump_operator(G_cav_L_R_drive)[2]
 
 # Here, the usual classical cavity drive-term $\sqrt{\kappa_L} E (a^\dagger + a)$ appears as a combination of Hamiltonian and Lindblad term. 
 # To solve the dynamics of the system we translate the symbolic expressions into numeric operators (matrices) of [QuantumOptics.jl](https://github.com/qojulia/QuantumOptics.jl). Since we do not want to include the basis of the atoms, we provide a dictionary of operators with the kwarg `operators` in the function [`to_numeric`](@ref). 
@@ -140,11 +140,11 @@ H2 = G_ac_drive.hamiltonian
 
 #
 
-L2_L = G_ac_drive.lindblad[1]
+L2_L = G_ac_drive.jump_operator[1]
 
 #
 
-L2_R = G_ac_drive.lindblad[2]
+L2_R = G_ac_drive.jump_operator[2]
 
 #
 

diff --git a/examples/04-2_two-sided-cavity_with-atom_coh-drive__cumulants.jl b/examples/04-2_two-sided-cavity_with-atom_coh-drive__cumulants.jl
@@ -45,11 +45,11 @@ H1 = hamiltonian(G_cav_L_R_drive)
 
 # 
 
-L1_L = lindblad(G_cav_L_R_drive)[1]
+L1_L = jump_operator(G_cav_L_R_drive)[1]
 
 # 
 
-L1_R = lindblad(G_cav_L_R_drive)[2]
+L1_R = jump_operator(G_cav_L_R_drive)[2]
 
 # The typical cavity drive-term $\sqrt{\kappa_L} E (a^\dagger + a)$ is a combination of Hamiltonian term and Lindblad. 
 # We use the function `meanfield` to obtain the equation for the intra-cavity field, which leads to a closed set of equations in this particular case.  
@@ -153,11 +153,11 @@ H2 = G_ac_drive.hamiltonian
 
 #
 
-L2_L = G_ac_drive.lindblad[1]
+L2_L = G_ac_drive.jump_operator[1]
 
 #
 
-L2_R = G_ac_drive.lindblad[2]
+L2_R = G_ac_drive.jump_operator[2]
 
 # We derive the equations of motion for system with a second-order mean-field approximation. 
 

diff --git a/examples/05-1_N-QDs_bidirectional-waveguide_coherent-pulse.jl b/examples/05-1_N-QDs_bidirectional-waveguide_coherent-pulse.jl
@@ -55,7 +55,7 @@ H = hamiltonian(G_t)
 
 #
 
-L = lindblad(G_t)
+L = jump_operator(G_t)
 L_R = L[1]
 
 #

diff --git a/examples/05-2_N-QDs_bidirectional-waveguide_quantum-pulse_qo.jl b/examples/05-2_N-QDs_bidirectional-waveguide_quantum-pulse_qo.jl
@@ -68,7 +68,7 @@ nothing # hide
 # The full Hamiltonian and Lindblad terms are extracted from the final SLH element. Note that as soon as one time-dependent function is involved in a cascade or concatenate, the returned $H$ and $L$ will also be time-dependent. 
 
 H = hamiltonian(G_t)
-L = lindblad(G_t)
+L = jump_operator(G_t)
 L_R = L[1]
 L_L = L[2]
 

diff --git a/examples/05-3_N-QDs_bidirectional-waveguide_feedback-reduction.jl b/examples/05-3_N-QDs_bidirectional-waveguide_feedback-reduction.jl
@@ -51,7 +51,7 @@ nothing # hide
 # second one to the transmitted right-moving output.
 
 H = hamiltonian(G_t)
-L = lindblad(G_t)
+L = jump_operator(G_t)
 L_L = L[1]
 L_R = L[2]
 nothing # hide

diff --git a/examples/06-1_interaction-picture__PRA2023_107-013706_fig2.jl b/examples/06-1_interaction-picture__PRA2023_107-013706_fig2.jl
@@ -42,7 +42,7 @@ H = hamiltonian(G_cas)
 
 #
 
-L = lindblad(G_cas)[1]
+L = jump_operator(G_cas)[1]
 
 # Usually we deal with the above derived Hamiltonian and Lindblad. In this example, however, we transform the system into the interaction picture of the virtual cavity-cavity interaction Hamiltonian $H_{uv}$. 
 

diff --git a/examples/07-1_beamsplitter_loss__quantum-pulse.jl b/examples/07-1_beamsplitter_loss__quantum-pulse.jl
@@ -51,7 +51,7 @@ H = hamiltonian(G)
 
 #
 
-L = lindblad(G)
+L = jump_operator(G)
 L[1]
 
 #

diff --git a/examples/07-2_hong-ou-mandel__quantum-pulse.jl b/examples/07-2_hong-ou-mandel__quantum-pulse.jl
@@ -50,7 +50,7 @@ H = hamiltonian(G)
 
 #
 
-L = lindblad(G)
+L = jump_operator(G)
 L[1]
 
 #

diff --git a/examples/08-1_pulse-delay__simple.jl b/examples/08-1_pulse-delay__simple.jl
@@ -49,7 +49,7 @@ G_v2 = concatenate(SLH(1, 0, 0), G_v)
 
 G_cas = cascade(G_u2, G_d, G_v2)
 H = hamiltonian(G_cas)
-L = lindblad(G_cas)
+L = jump_operator(G_cas)
 nothing # hide
 
 #

diff --git a/examples/09-1_coherent-feedback-squeezing__Gough-Wildfeuer-2009.jl b/examples/09-1_coherent-feedback-squeezing__Gough-Wildfeuer-2009.jl
@@ -36,7 +36,7 @@ nothing # hide
 #
 
 S_loop = scattering(G_loop)
-L_loop = lindblad(G_loop)[1]
+L_loop = jump_operator(G_loop)[1]
 
 #
 

diff --git a/examples/10-1_SUPER_excitation.jl b/examples/10-1_SUPER_excitation.jl
@@ -50,7 +50,7 @@ nothing # hide
 
 ## Hamiltonian and Lindbladian
 Hcas = hamiltonian(G_cas)
-Lcas = lindblad(G_cas)[1]
+Lcas = jump_operator(G_cas)[1]
 Lcasd = adjoint(Lcas)
 nothing # hide 
 

diff --git a/examples/drafts/02-2_traveling-cat-state__PRA2020_102-023717_fig3b.jl b/examples/drafts/02-2_traveling-cat-state__PRA2020_102-023717_fig3b.jl
@@ -38,7 +38,7 @@ H_s = p / 2 * (a'^2 + a^2) - K / 2 * (a'^2) * (a^2) + Δ * a' * a
 G_s = SLH(1, √(γ) * a, H_s)
 
 H = hamiltonian(G_s)
-L = lindblad(G_s)[1]
+L = jump_operator(G_s)[1]
 nothing # hide
 
 # Next, we define the numerical parameters from Sec. II.D. The classical pump
@@ -127,7 +127,7 @@ G_v = SLH(1, g_v * av, 0)
 G = G_s2 ▷ G_v
 
 H_2 = hamiltonian(G)
-L_2 = lindblad(G)[1]
+L_2 = jump_operator(G)[1]
 
 gv_t = coupling_output(v_mode, T)
 dict_p_t_2 = Dict(p => p_t, g_v => gv_t)

diff --git a/examples/drafts/08-2_pulse-propagation-delay__PRA2026_113-013730_fig4.jl b/examples/drafts/08-2_pulse-propagation-delay__PRA2026_113-013730_fig4.jl
@@ -66,7 +66,7 @@ G_atom = SLH(I4, L_atom, Δ*σ(2, 2))
 
 G_u_bs_d1_d2_atom = G_u_bs_d1_d2 ▷ G_atom
 H = hamiltonian(G_u_bs_d1_d2_atom)
-L = lindblad(G_u_bs_d1_d2_atom)
+L = jump_operator(G_u_bs_d1_d2_atom)
 
 #
 

diff --git a/examples/drafts/10-1_SUPER_excitation.jl b/examples/drafts/10-1_SUPER_excitation.jl
@@ -184,7 +184,7 @@ G_cas = ▷(G_u2, G_u1, G_2lvl, G_v1, G_v2) # cascade
 
 # Hamiltonian and Lindbladian
 Hcas = hamiltonian(G_cas)
-Lcas = lindblad(G_cas)[1]
+Lcas = jump_operator(G_cas)[1]
 Lcasd = adjoint(Lcas)
 
 # Time-dependent couplings
@@ -577,8 +577,8 @@ Gcon = ⊞(G_ch1, G_ch2) # concatenation
 
 # Hamiltonian and Lindblad operators (two channels)
 Hcon = hamiltonian(Gcon)
-Lch1g = lindblad(Gcon)[1]
-Lch2g = lindblad(Gcon)[2]
+Lch1g = jump_operator(Gcon)[1]
+Lch2g = jump_operator(Gcon)[2]
 Ldch1g = adjoint(Lch1g)
 Ldch2g = adjoint(Lch2g)
 
@@ -734,7 +734,7 @@ G_cas_1m = ▷(G_u, G_2lvls, G_v)
 
 # Hamiltonian and Lindbladian
 Hcas_1m = hamiltonian(G_cas_1m)
-Lcas_1m = lindblad(G_cas_1m)[1]
+Lcas_1m = jump_operator(G_cas_1m)[1]
 Lcas_1md = adjoint(Lcas_1m)
 
 # Time-dependent couplings
@@ -869,7 +869,7 @@ T = [dt:dt:Tend;]
 G_cas_1m0 = ▷(G_u, G_2lvls, G_v)
 
 Hcas_1m0 = hamiltonian(G_cas_1m0)
-Lcas_1m0 = lindblad(G_cas_1m0)[1]
+Lcas_1m0 = jump_operator(G_cas_1m0)[1]
 Lcas_1m0d = adjoint(Lcas_1m0)
 
 # Set v-coupling to zero (collect only input-system correlations)
@@ -1201,7 +1201,7 @@ G_cas_1m = ▷(G_u, G_2lvls, G_v)
 
 # Hamiltonian and Lindbladian
 Hcas_1m = hamiltonian(G_cas_1m)
-Lcas_1m = lindblad(G_cas_1m)[1]
+Lcas_1m = jump_operator(G_cas_1m)[1]
 Lcas_1md = adjoint(Lcas_1m)
 
 # Time-dependent couplings

diff --git a/examples/drafts/11-1_two-level-system_spontaneous-emission.jl b/examples/drafts/11-1_two-level-system_spontaneous-emission.jl
@@ -36,7 +36,7 @@ G_v = SLH(1, g_v' * av, 0)
 G = cascade(G_s, G_v)
 
 H = hamiltonian(G)
-L = lindblad(G)[1]
+L = jump_operator(G)[1]
 nothing # hide
 
 # We use the spontaneous-emission rate as the frequency unit. At

diff --git a/src/QuantumInputOutput.jl b/src/QuantumInputOutput.jl
@@ -18,6 +18,7 @@ const SQA = SecondQuantizedAlgebra
 export SLH,
     Gaussian,
     scattering,
+    jump_operator,
     lindblad,
     hamiltonian,
     # Composition
-Original file line number
+Diff line change
@@ Expand Up / @@ -11,6 +11,7 @@ scattering @@
     ```
     ```@docs
+    jump_operator
     lindblad
     ```
@@ Expand Down @@