sqzcomb

Squeezed light in Kerr microcombs, computed end to end: from the classical Lugiato-Lefever steady state, through the linearized quantum fluctuations around it, to the output quadrature-noise spectrum a homodyne detector would report. The package exists because intracavity squeezing is not the observable; what leaves the extraction port is, and the coupling that maximizes one does not maximize the other.

Entanglement of the comb (new in v0.6)

The entangle module quantifies twin-beam entanglement directly from the covariance matrices the package produces: the partial-transpose symplectic eigenvalue through the Simon invariants (necessary and sufficient for two-mode Gaussian states), the logarithmic negativity, and the symmetric Duan-Simon EPR sum with its separability bound. The invariant formula is cross-checked against an independent explicit-partial-transpose computation, against the closed-form two-mode squeezed vacuum (E_N = 2r exactly), and on the driven photonic molecule, whose two rings turn out to be PPT-entangled below threshold while the symmetric Duan sum misses it, a working demonstration of why the sharper criterion matters.

from sqzcomb import (photonic_molecule, intracavity_covariance,
                     covariance_xxpp, entanglement_report)
M, gammas = photonic_molecule(mu=0.8, J=1.0)
sigma = covariance_xxpp(intracavity_covariance(M, gammas))
print(entanglement_report(sigma, 0, 1))

Status

v0.6.0 (alpha). Implemented and tested:

• Lugiato-Lefever solver (Strang splitting; the Kerr step and the linear-plus-pump step are each exact)
• homogeneous steady states against the exact cubic
• linearized fluctuation (Bogoliubov) matrix around an arbitrary steady state, with a stability guard that refuses above-threshold states
• input-output quadrature spectra with an extraction port and intrinsic loss; single-mode and joint two-mode quadratures
• photonic molecule (new in v0.2): the two-ring coupled-mode model (photonic_molecule), per-mode-coupling output spectra (output_variance_ports), and the exact instability threshold (molecule_threshold)
• multimode comb molecule (new in v0.3): every retained comb line of the LLE fluctuation matrix coupled to a matching auxiliary-ring mode (molecule_fluctuation_matrix), with multi-line bus detection and joint twin-beam quadratures through the auxiliary ring
• Gaussian-state interop (new in v0.4): steady-state covariance matrices of the intracavity Gaussian state (intracavity_covariance, numpy-only Lyapunov solve), export in the standard xxpp quadrature ordering with an explicit hbar convention (covariance_xxpp, vacuum exactly the identity at hbar = 2), Williamson symplectic spectra (symplectic_eigenvalues), and a QuTiP adapter (drift_from_qutip) that turns any quadratic QuTiP Hamiltonian into a drift matrix this package's spectra machinery accepts, refusing non-quadratic Hamiltonians rather than silently linearizing them. Install sqzcomb[interop] for the adapter; its test asserts that the released two-ring molecule is reproduced exactly from a QuTiP Hamiltonian.

Verified against closed forms in the test-suite: vacuum passes a passive cavity, and a passive molecule, unchanged for every coupling, port and frequency; the degenerate parametric oscillator output spectrum is reproduced to 1e-10; and the textbook result that detectable squeezing saturates at 3 dB at critical coupling, while full extraction breaks that limit, emerges from the machinery rather than being asserted.

For the molecule, the test-suite additionally asserts: exact reduction to the single ring at zero coupling; passive supermodes split by exactly 2J; the resonant threshold mu = 1 + J^2/gamma (static branch, Hopf branch 1 + gamma beyond J = gamma); the quarter-turn quadrature rotation of the -iJ hop; the exact zero-frequency equivalence of the auxiliary-ring port to an effective single mode with escape efficiency (J^2/gamma)/(1 + J^2/gamma); and a J^2/gamma = 3 molecule reaching 6 dB detected squeezing through the auxiliary port although the Kerr ring itself has no extraction port, which is the molecule extraction mechanism in its simplest form.

For the multimode molecule, the asserts continue in the same spirit: at one retained line the builder equals the released two-ring matrix to machine precision; at zero coupling it reduces exactly to the plain fluctuation matrix and the v0.1 spectra, single-line and twin-beam; a passive multimode molecule returns exact vacuum through any bus; the resonant auxiliary ring at zero frequency is exactly the single ring with J^2/gamma_b extra loss per line and the quarter-turn rotation; and when J^2/gamma_b exceeds one, twin-beam squeezing detected through the auxiliary bus is strictly deeper than through the main bus of the same device. The stability guard is also asserted to refuse a flat state that is above a pair's modulational-instability threshold once the molecule's added loss is removed.

Later releases are documented in their own sections below: imperfect detection (v0.5) and two-mode Gaussian entanglement of the comb (v0.6).

Not yet implemented (the roadmap, in order): soliton-crystal steady-state continuation, supermode decomposition of the multimode covariance, and thermal input noise.

Install and use

pip install sqzcomb

For development, clone the repository and pip install -e .[test].

import numpy as np
from sqzcomb import (lle_evolve, fluctuation_matrix,
                     output_quadrature_variance, squeezing_db)

# steady state at pump F and detuning alpha, anomalous dispersion d2
psi = lle_evolve(np.full(256, 0.05 + 0j), F=1.2, alpha=0.8,
                 dispersion=(-0.02,), t_end=300.0)

M, modes = fluctuation_matrix(psi, alpha=0.8, dispersion=(-0.02,))
i0 = int(np.where(modes == 0)[0][0])
v = output_quadrature_variance(M, eta=0.5, omega=0.0,
                               mode_index=i0, n_modes=modes.size)
print(squeezing_db(v), "dB relative to vacuum")

Units are the standard normalized LLE units: time in photon lifetimes, eta = kappa_ex / kappa, vacuum variance 1/2.

Imperfect detection (new in v0.5)

Every spectrum above is the noise at the extraction port; the photodiodes report less. detection applies the standard beamsplitter model of optical loss and quantum efficiency plus additive electronic noise, in both languages the package speaks: as a scalar map on quadrature variances (detected_variance, detected_squeezing_db, with dark_from_clearance_db converting a receiver's dark clearance in dB into a variance) and as the lossy Gaussian channel on xxpp covariance matrices (lossy_channel_xxpp, per-mode efficiencies allowed). required_efficiency inverts the loss model into the number an experiment plans around: the minimum efficiency that still delivers a target squeezing from a given source. Loss stages compose by multiplying efficiencies, and the test suite asserts that composition exactly, along with the vacuum fixed point, physicality of the channel (symplectic eigenvalues never fall below hbar/2), and agreement of the scalar and matrix forms on a squeezed mode.

References: the beamsplitter model of detector inefficiency, U. Leonhardt, Measuring the Quantum State of Light (Cambridge, 1997); the Gaussian lossy channel, C. Weedbrook et al., Rev. Mod. Phys. 84, 621 (2012).

Methodological basis

T. M. Mahim, M. M. Rahman and A. S. M. Mohsin, "Overcoming the 3 dB squeezing extraction limit in silicon carbide microcombs with a photonic molecule" (under review); code for the paper: https://github.com/Tanvir-Mahmud-Mahim/sic-molecule-squeezer

This package is the general-purpose engine; the paper repository reproduces the specific published study. v0.2 adds the two-mode photonic molecule, the extraction mechanism in its simplest form; the paper's full multimode comb molecule remains in the paper repository.

Support and governance

The package is written and maintained by Tanvir Mahmud Mahim (Department of Electrical and Electronic Engineering, BRAC University), who reviews every change and takes the final decision on scope and releases. There is no separate governance body; design questions are discussed in the open in issues and pull requests, and the standing rule of CONTRIBUTING.md binds the maintainer exactly as it binds contributors: a change that touches physics arrives with a test, and a constant arrives with its source.

Support runs through the issue tracker at https://github.com/TaN-MM-Org/sqzcomb/issues. Usage questions are welcome there alongside bug reports; a docstring that left a unit or a sign convention unclear is treated as a documentation bug, not as user error. The maintainer aims to respond within a week.

While the version is below 1.0 the API may still move between minor versions; such changes are called out in the release notes. The normalized-unit conventions stated above are stable: any change to them would be a breaking change named in the release notes, never a quiet renormalization.

License

Apache-2.0