Gauge transformation for pulse propagation and time ordered integrals

This paper introduces a gauge transformation that eliminates time-dependent onsite potentials by renormalizing hopping terms with phase factors, thereby simplifying time-ordered integrals and facilitating the simulation of pulse propagation in scattering systems.

The Gist

Imagine you are trying to predict how a crowd of people (electrons) moves through a complex building (a quantum system) when a sudden, loud announcement (a time-dependent pulse or voltage) is made over the intercom.

In the world of quantum physics, this is a nightmare to calculate. The standard math requires tracking every single person's movement while the announcement is happening, which involves incredibly messy, "time-ordered" calculations. It's like trying to calculate the path of every drop of water in a river while the wind is constantly changing direction.

This paper, by Adel Abbout, introduces a clever mathematical "magic trick" (called a Gauge Transformation) that simplifies this chaos. Here is how it works, broken down into simple concepts:

1. The Problem: The "Noisy" Room

Imagine a long hallway (a wire or lead) connected to a central room (the device you are studying). Suddenly, the lights in the hallway start flickering wildly (a time-dependent potential).

• The Old Way: To see how people move from the hallway into the room, you have to calculate the effect of those flickering lights on every single step taken by every person in the infinite hallway. It's computationally exhausting.
• The Math Struggle: The equations involve "time-ordered integrals." Think of this as trying to solve a puzzle where the pieces change shape every time you look at them. You can't just multiply them; you have to do it in a very specific, complex order.

2. The Solution: The "Silent" Trick

The author suggests a trick: Change your perspective.

Instead of watching the hallway lights flicker, imagine you put on a pair of special glasses (the Gauge Transformation).

• The Effect: Through these glasses, the flickering lights in the hallway disappear completely. The hallway looks perfectly calm and steady.
• The Catch: Where did the flickering go? It didn't vanish; it moved. The "flicker" is now concentrated entirely at the doorway (the interface) between the hallway and the central room.

3. How the Magic Works (The Phase Factor)

In physics terms, the author shows that you can mathematically "erase" the changing voltage at a specific spot.

• The Renormalization: When you erase the voltage at a spot, the "hopping" (movement) of electrons into and out of that spot gets a little "tag" or "phase factor" attached to it.
• The Analogy: Imagine people walking through a door. If the lights flicker inside the room, it's hard to track them. But if you say, "Okay, the room is dark and steady, but everyone entering or leaving the room has to do a little dance step (a phase shift) to compensate," the math becomes much easier.
• The Cancellation: If you do this for every spot in the infinite hallway, the "dance steps" people do when moving from one spot to the next cancel each other out. The hallway becomes invisible to the math, and the only place where the "dance" (the time-dependence) remains is right at the door where the hallway meets the room.

4. Why This Matters

This trick is a game-changer for two main reasons:

• Simplifying the Infinite: It allows scientists to simulate systems with infinite leads (like real-world wires) without having to do infinite calculations. They only need to calculate the messy part at the interface. It turns a problem involving an infinite number of variables into a problem involving just a few.
• Untangling the Knots: The "time-ordered integrals" mentioned earlier are like a tangled ball of yarn. This transformation cuts out the knots (the self-loops or onsite potentials) and leaves you with a straight, clean string. It makes the math for how the system evolves over time much faster and simpler to solve.

5. Real-World Applications

The paper gives a few examples of where this is useful:

• Pulse Propagation: Sending a signal down a wire. Instead of tracking the signal changing everywhere, you just track how it enters the device.
• Spin Pumping: Imagine a spinning top (a spin) that is wobbling. You can mathematically shift the "wobble" from the top itself to the connection point, making it much easier to calculate how much energy is being transferred.

The Bottom Line

Adel Abbout's paper is essentially saying: "Don't fight the noise everywhere. Move the noise to the edge, and let the rest of the system be quiet."

By shifting the complexity from the entire system to just the boundary, this method makes it possible to simulate complex quantum devices that were previously too difficult to calculate, opening the door to better designs for quantum computers and advanced electronic devices.

Technical Summary

Here is a detailed technical summary of the paper "Gauge transformation for pulse propagation and time ordered integrals" by Adel Abbout.

1. Problem Statement

The paper addresses two primary challenges in the simulation of mesoscopic quantum systems subject to time-dependent perturbations (such as electric pulses, periodic driving, or sudden potential changes):

• Complexity of Time-Dependent Schrödinger Equation (TDSE): Solving the TDSE for non-equilibrium states often requires handling infinite sums or integrals representing all possible absorption/emission processes, as energy is not conserved under external driving.
• Computational Burden of Time-Ordered Integrals: The evolution operator $U(t)$ involves time-ordered integrals ($T \exp[-i \int H(t') dt']$). When the Hamiltonian at different times does not commute ($[H(t), H(t')] \neq 0$), these integrals are difficult to evaluate analytically and often require expensive numerical methods.
• Handling Infinite Systems: Standard solvers (e.g., tkwant) struggle to apply time-dependent potentials directly to semi-infinite leads. They typically require gauge transformations on all sites of a lead, which can be computationally redundant if not optimized.

2. Methodology

The author proposes a successive gauge transformation technique based on the elimination of time-dependent onsite potentials at individual lattice sites.

• Hamiltonian Formulation: The system is described by a tight-binding Hamiltonian (not necessarily Hermitian):
  $$H = -\sum_{(k,k') \in \text{hop}} \gamma_{kk'} c^\dagger_k c_{k'} - e \sum_{k \in \text{list}} V_k(t) c^\dagger_k c_k$$
• The Transformation Operator: For a specific site $l$, the author defines a unitary operator $U_l$ to eliminate the potential $V_l(t)$:
  $$U_l = e^{-i\phi_l(t) c^\dagger_l c_l} = 1 + (e^{-i\phi_l(t)} - 1)c^\dagger_l c_l$$
  where the phase is defined as $\phi_l(t) = \frac{e}{\hbar} \int_{-\infty}^t V_l(t') dt'$.
• Renormalization of Hoppings: Applying the transformation $\tilde{H} = U^\dagger_l H U_l - i\hbar U^\dagger_l \partial_t U_l$ results in the removal of the onsite potential $V_l(t)$ but introduces phase factors into the hopping terms connected to site $l$:
  • Outward hopping (from $l$ to $k'$): Renormalized by $e^{-i\phi_l(t)}$.
  • Inward hopping (from $k'$ to $l$): Renormalized by $e^{+i\phi_l(t)}$.
  • Hoppings not involving site $l$ remain unchanged.
• Graph-Theoretic Approach: The time evolution operator is analyzed using the $\star$-product (star product) formalism and prime walks (cycles that do not visit a site more than once). The onsite potentials are treated as "self-loops" in the Hamiltonian graph. The gauge transformation effectively removes these self-loops, simplifying the path-sum expansion.

3. Key Contributions

• Generalized Gauge Transformation: The method is applicable to both finite and infinite systems and does not require the Hamiltonian to be Hermitian. It allows for the partial removal of potentials, not just the total elimination.
• Interface Localization of Time Dependence: The paper demonstrates that for a system with a central scattering region connected to semi-infinite leads, applying a uniform time-dependent potential to a lead can be transformed such that the time dependence vanishes from the lead's interior and concentrates solely at the interface between the lead and the central system.
• Simplification of Time-Ordered Integrals: By eliminating onsite potentials (self-loops) via gauge transformation, the number of complex integrals in the $\star$-product expansion is significantly reduced. This avoids the explicit computation of Neumann expansions involving non-commutative terms like $(1 - c)^{-1}$.
• Inverse Application (Phase to Potential): The method can be reversed: a time-dependent phase factor in a hopping term can be transformed into a time-dependent onsite potential. This is used to convert a time-dependent spin precession Hamiltonian into a time-independent one.

4. Results and Illustrative Examples

The paper validates the methodology through several scenarios:

1. Pulse Propagation in Scattering Systems:

   • In a setup with a central region and two semi-infinite leads, if a pulse $V(t)$ is applied to one lead, the gauge transformation renders the hopping terms inside the lead time-independent (phases cancel out between neighbors).
   • The time dependence is confined to the interface hopping (the bond connecting the lead to the central region), which acquires a phase factor $e^{\pm i\phi(t)}$. This drastically reduces the computational domain requiring time-dependent treatment.

2. Uniform Potential in Finite Systems:

   • For a ring-shaped system with a uniform time-dependent onsite potential, the transformation eliminates the potential everywhere inside the ring.
   • Internal hoppings remain unchanged because the phase factors from neighboring sites cancel exactly. Only the interface hoppings (if connected to external leads) retain the phase factors.

3. Spin Precession:

   • A single spin precessing at frequency $\omega$ is described by a time-dependent Hamiltonian. By applying the inverse gauge transformation (removing the phase from the hopping), the Hamiltonian is converted into a time-independent form with shifted energy levels ($\pm \omega/2$).
   • This simplifies the calculation of observables like spin density $\langle \sigma_z \rangle$ and facilitates the analysis of spin pumping in ferromagnetic chains adjacent to normal metals.

4. Reduction of Path-Sum Complexity:

   • In the graph representation of the Hamiltonian, removing self-loops (onsite potentials) reduces the number of cycles in the path-sum expansion. The resulting expression involves fewer renormalized hopping terms, making the evaluation of the time-evolution operator more efficient.

5. Significance

• Computational Efficiency: The method offers a powerful tool for simulating time-dependent quantum transport, particularly for pulse propagation and spin pumping, by reducing the spatial extent of time-dependent terms to the system boundaries (interfaces).
• Analytical Tractability: It transforms complex time-dependent problems into simpler forms, sometimes converting time-dependent Hamiltonians into time-independent ones, allowing for exact analytical solutions where numerical integration would be required.
• Versatility: The approach is robust for non-Hermitian systems and general graph structures, making it applicable to disordered systems (Random Matrix Theory) and complex mesoscopic devices.
• Theoretical Insight: It clarifies the physical equivalence between uniform potential shifts in leads and interface phase modulations, providing a deeper understanding of gauge invariance in non-equilibrium quantum transport.

In conclusion, Adel Abbout presents a rigorous gauge transformation framework that simplifies the mathematical structure of time-dependent quantum problems, offering significant advantages for both analytical derivations and numerical simulations of pulse propagation and spin dynamics.