Nonlinear Two-Photon Quantum Battery

• The paper demonstrates that nonlinear two-photon interactions enable ultrafast charging dynamics and enhanced energy storage compared to traditional linear protocols.
• The methodology relies on engineered quadratic couplings in platforms like circuit-QED and photonic cavities, which induce squeezed states and collective Rabi oscillations.
• Key results reveal exponential energy scaling and maximal ergotropy from pure Gaussian states, paving the way for advanced hybrid quantum technologies.

A nonlinear two-photon driving quantum battery model utilizes nonlinear photon-photon interactions to implement energy transfer and storage protocols with substantial advantages over single-photon counterparts. These models leverage engineered quadratic (two-photon) couplings in cavity-QED, circuit-QED, and related bosonic systems, as well as in ensembles of two-level systems (TLS) coupled to a resonator. Central features include ultrafast charging dynamics, generation of nonclassical (squeezed) states, and enhanced scaling of stored energy and available work (ergotropy). The following sections synthesize definitions, formalism, dynamical analysis, performance metrics, and implementation considerations, with direct reference to the established literature (Crescente et al., 2020, Downing et al., 2023, Downing et al., 2024, Downing et al., 24 Oct 2025, Xu et al., 15 Nov 2025).

1. Theoretical Framework and Model Hamiltonians

Nonlinear two-photon quantum battery models generalize first-order (dipolar, “linear”) charging by incorporating quadratic coupling—driving or interaction terms proportional to $a^2 + a^{\dagger 2}$ (bosonic systems) or $[(a^\dagger)^2 + a^2]$ coupled to collective spin operators (Dicke-type models). Two principal classes arise:

(a) Collective Two-Photon Dicke Model

The system consists of $N$ two-level systems (TLS), described by collective spin operators $$J_\alpha = \frac{1}{2} \sum_{i=1}^N \sigma_\alpha^{(i)}$$, coupled to a single quantized mode (resonator) of frequency $\omega_c$: $$H = \omega_c a^\dagger a + \omega_a J_z + g_1 J_x (a^\dagger + a) + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$ Here, $g_1$ and $g_2$ parameterize single-photon and two-photon coupling strengths. Suppressing $g_1$ yields the pure two-photon Dicke Hamiltonian: $$H_{2ph} = \omega_c a^\dagger a + \omega_a J_z + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$ Typically, $[(a^\dagger)^2 + a^2]$0 is engineered via circuit-QED architectures (e.g., symmetric SQUIDs, bichromatic driving) to maximize nonlinear exchange (Crescente et al., 2020).

(b) Bosonic Two-Mode Quadratic Driven Model

A bosonic "charger" ($[(a^\dagger)^2 + a^2]$1) is coupled (either linearly or nonlinearly) to a "battery" ($[(a^\dagger)^2 + a^2]$2), both of frequency $[(a^\dagger)^2 + a^2]$3, with a time-dependent or static quadratic drive: $[(a^\dagger)^2 + a^2]$4 $[(a^\dagger)^2 + a^2]$5 realizes a two-photon (parametric) drive; $[(a^\dagger)^2 + a^2]$6 enables energy transfer between charger and battery (Downing et al., 2023, Downing et al., 24 Oct 2025). In some cases, direct quadratic ("downconversion") coupling $[(a^\dagger)^2 + a^2]$7 between charger and battery is considered (Downing et al., 24 Oct 2025).

2. Dynamical Evolution, Squeezing, and Charging Dynamics

Squeezing Generation and State Evolution

The quadratic drive executes a squeezing operation on the recipient mode (battery), with the propagator

$[(a^\dagger)^2 + a^2]$8

where $[(a^\dagger)^2 + a^2]$9 is the two-photon drive amplitude (Downing et al., 2024). For a bosonic battery initially in vacuum, squeezing populates higher Fock states, leading to

$N$0

i.e., the stored energy grows exponentially with the integrated pulse area $N$1. For time-dependent protocols (e.g., Gaussian pulse envelopes), both the time-maximum stored energy and maximum charging power display exponential scaling in the drive parameters.

Spin-Boson Collective Dynamics

In the two-photon Dicke model, with all TLSs initially in their ground state and the cavity loaded with $N$2 photons, the system undergoes collective Rabi oscillations at the frequency $N$3. The average stored energy evolves as

$N$4

with fluctuations $N$5 (Crescente et al., 2020).

3. Performance Metrics: Energy, Power, and Ergotropy

Key figures of merit for quantum batteries include the maximum stored energy ($N$6), maximum and average charging power ($N$7, $N$8), and ergotropy (extractable work, $N$9).

• Scaling Laws:
  • In the two-photon Dicke model:
  • $$J_\alpha = \frac{1}{2} \sum_{i=1}^N \sigma_\alpha^{(i)}$$0
  • $$J_\alpha = \frac{1}{2} \sum_{i=1}^N \sigma_\alpha^{(i)}$$1
  • Charging time $$J_\alpha = \frac{1}{2} \sum_{i=1}^N \sigma_\alpha^{(i)}$$2
  • For linear models, $$J_\alpha = \frac{1}{2} \sum_{i=1}^N \sigma_\alpha^{(i)}$$3 and $$J_\alpha = \frac{1}{2} \sum_{i=1}^N \sigma_\alpha^{(i)}$$4 (Crescente et al., 2020).
• Ergotropy and Minimum-Uncertainty States:
  • In continuous-variable two-photon models, squeezing induces pure minimum-uncertainty (Gaussian) states of the battery mode, saturating Heisenberg's relation, and thus all stored energy is extractable:

  $$J_\alpha = \frac{1}{2} \sum_{i=1}^N \sigma_\alpha^{(i)}$$5 - The determinant of the covariance matrix $$J_\alpha = \frac{1}{2} \sum_{i=1}^N \sigma_\alpha^{(i)}$$6 indicates absence of passive (non-extractable) energy (Downing et al., 24 Oct 2025). - For general Gaussian states, ergotropy is

  $$J_\alpha = \frac{1}{2} \sum_{i=1}^N \sigma_\alpha^{(i)}$$7

  (Downing et al., 2023).

• Nonreciprocal and Dissipative Effects:

  • Optimizing dissipation asymmetry (suppressing battery loss) and balancing collective dissipation channels enhances both $$J_\alpha = \frac{1}{2} \sum_{i=1}^N \sigma_\alpha^{(i)}$$8 and $$J_\alpha = \frac{1}{2} \sum_{i=1}^N \sigma_\alpha^{(i)}$$9 (Xu et al., 15 Nov 2025).

4. Comparison with Single-Photon Driving and Other Charging Mechanisms

Direct benchmarking against linear (single-photon) protocols reveals qualitative and quantitative superiority of two-photon models in several metrics:

• The two-photon Dicke protocol delivers an additional $\omega_c$0 speedup in charging and $\omega_c$1 enhancement in power, on top of the usual Dicke collective scaling (Crescente et al., 2020).
• In bosonic Gaussian batteries, both stored energy and ergotropy can become parametrically larger near critical points of the driven-dissipative system, while some energy becomes passive and unextractable near but not at minimum-uncertainty (Downing et al., 2023).
• Under strong two-photon driving in nonreciprocal models, the absolute energy storage and ergotropy advantages persist, though conversion efficiency per unit supplied energy may be lower than in linear cases (Xu et al., 15 Nov 2025).

5. Open-System and Nonequilibrium Considerations

Open quantum system dynamics, encompassing Markovian loss, engineered dissipation, and environmental coupling, define available performance regimes and physical limits:

• Master equations of Lindblad type describe decay (rates: $\omega_c$2), cross-dissipation (rate: $\omega_c$3), and engineered nonreciprocal energy flow (Xu et al., 15 Nov 2025).
• Dynamical equilibration times, steady-state existence, and critical points are set by interplay of two-photon drive strength and decay rates. Stability requires $\omega_c$4, where $\omega_c$5 is the relevant total loss rate (Xu et al., 15 Nov 2025).
• At strong driving and near dissipative exceptional points, parametrically large energy storage is feasible, while real devices require technical regularization to avoid unphysical divergences (Downing et al., 2023).

6. Experimental Implementation and Platform-Specific Parameters

Realization of nonlinear two-photon quantum batteries requires platforms supporting significant two-photon (parametric) interactions:

• Circuit QED: Superconducting circuits with flux qubits, symmetric SQUIDs, or Josephson parametric amplifiers; frequencies $\omega_c$6–$\omega_c$7 GHz; $\omega_c$8 up to $\omega_c$9–$$H = \omega_c a^\dagger a + \omega_a J_z + g_1 J_x (a^\dagger + a) + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$0; drive amplitudes up to $$H = \omega_c a^\dagger a + \omega_a J_z + g_1 J_x (a^\dagger + a) + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$1 MHz; coherence times $$H = \omega_c a^\dagger a + \omega_a J_z + g_1 J_x (a^\dagger + a) + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$2s–ms (Crescente et al., 2020, Downing et al., 2023, Xu et al., 15 Nov 2025).
• Photonic Microcavities: Second-order nonlinear ($$H = \omega_c a^\dagger a + \omega_a J_z + g_1 J_x (a^\dagger + a) + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$3) crystals for optical frequency implementations; parametric downconversion for two-photon driving (Downing et al., 2024).
• Trapped Ions and Magnonic Resonators: Bichromatic driving and engineered couplings enable implementation of the quadratic terms (Crescente et al., 2020, Xu et al., 15 Nov 2025).
• Optomechanical/Nanomechanical: Two-phonon modulation mechanisms support equivalent quadratic interactions (Downing et al., 2024).

A table summarizing representative parameter regimes:

Platform	Two-Photon Drive (G, $$H = \omega_c a^\dagger a + \omega_a J_z + g_1 J_x (a^\dagger + a) + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$4)	Loss Rate ($$H = \omega_c a^\dagger a + \omega_a J_z + g_1 J_x (a^\dagger + a) + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$5, $$H = \omega_c a^\dagger a + \omega_a J_z + g_1 J_x (a^\dagger + a) + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$6)	Coupling ($$H = \omega_c a^\dagger a + \omega_a J_z + g_1 J_x (a^\dagger + a) + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$7, $$H = \omega_c a^\dagger a + \omega_a J_z + g_1 J_x (a^\dagger + a) + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$8)
Circuit QED	$$H = \omega_c a^\dagger a + \omega_a J_z + g_1 J_x (a^\dagger + a) + g_2 J_x \left[ (a^\dagger)^2 + a^2 \right]$$9–$g_1$0 MHz	$g_1$1–$g_1$2 MHz	$g_1$3–$g_1$4 MHz
Photonic Cavity	$g_1$5–$g_1$6 GHz	$g_1$7–$g_1$8 MHz	platform dependent
Magnonic Resonators	$g_1$9–$g_2$0 MHz	$g_2$1–$g_2$2 MHz	$g_2$3 $g_2$4
Optomechanics	variable	$g_2$5 resonance frequency	system specific

7. Outlook and Implications

Nonlinear two-photon quantum battery models provide a platform for ultrafast, high-capacity quantum energy storage, with superextensive scaling properties and the intrinsic generation of nonclassical states (notably, squeezed minimum-uncertainty Gaussian states). These features yield maximal extractable work and open prospects for hybrid quantum technologies leveraging stored energy and quantum correlations for sensing or computation. Fine control of parametric driving, environmental couplings, and two-photon interactions is central to optimizing device performance and realizing experimental prototypes consistent with theoretical predictions (Crescente et al., 2020, Downing et al., 2023, Downing et al., 24 Oct 2025, Xu et al., 15 Nov 2025, Downing et al., 2024).