sqzcomb 0.1.1__tar.gz

sqzcomb-0.1.1/LICENSE

@@ -0,0 +1,202 @@
+
+                                 Apache License
+                           Version 2.0, January 2004
+                        http://www.apache.org/licenses/
+
+   TERMS AND CONDITIONS FOR USE, REPRODUCTION, AND DISTRIBUTION
+
+   1. Definitions.
+
+      "License" shall mean the terms and conditions for use, reproduction,
+      and distribution as defined by Sections 1 through 9 of this document.
+
+      "Licensor" shall mean the copyright owner or entity authorized by
+      the copyright owner that is granting the License.
+
+      "Legal Entity" shall mean the union of the acting entity and all
+      other entities that control, are controlled by, or are under common
+      control with that entity. For the purposes of this definition,
+      "control" means (i) the power, direct or indirect, to cause the
+      direction or management of such entity, whether by contract or
+      otherwise, or (ii) ownership of fifty percent (50%) or more of the
+      outstanding shares, or (iii) beneficial ownership of such entity.
+
+      "You" (or "Your") shall mean an individual or Legal Entity
+      exercising permissions granted by this License.
+
+      "Source" form shall mean the preferred form for making modifications,
+      including but not limited to software source code, documentation
+      source, and configuration files.
+
+      "Object" form shall mean any form resulting from mechanical
+      transformation or translation of a Source form, including but
+      not limited to compiled object code, generated documentation,
+      and conversions to other media types.
+
+      "Work" shall mean the work of authorship, whether in Source or
+      Object form, made available under the License, as indicated by a
+      copyright notice that is included in or attached to the work
+      (an example is provided in the Appendix below).
+
+      "Derivative Works" shall mean any work, whether in Source or Object
+      form, that is based on (or derived from) the Work and for which the
+      editorial revisions, annotations, elaborations, or other modifications
+      represent, as a whole, an original work of authorship. For the purposes
+      of this License, Derivative Works shall not include works that remain
+      separable from, or merely link (or bind by name) to the interfaces of,
+      the Work and Derivative Works thereof.
+
+      "Contribution" shall mean any work of authorship, including
+      the original version of the Work and any modifications or additions
+      to that Work or Derivative Works thereof, that is intentionally
+      submitted to Licensor for inclusion in the Work by the copyright owner
+      or by an individual or Legal Entity authorized to submit on behalf of
+      the copyright owner. For the purposes of this definition, "submitted"
+      means any form of electronic, verbal, or written communication sent
+      to the Licensor or its representatives, including but not limited to
+      communication on electronic mailing lists, source code control systems,
+      and issue tracking systems that are managed by, or on behalf of, the
+      Licensor for the purpose of discussing and improving the Work, but
+      excluding communication that is conspicuously marked or otherwise
+      designated in writing by the copyright owner as "Not a Contribution."
+
+      "Contributor" shall mean Licensor and any individual or Legal Entity
+      on behalf of whom a Contribution has been received by Licensor and
+      subsequently incorporated within the Work.
+
+   2. Grant of Copyright License. Subject to the terms and conditions of
+      this License, each Contributor hereby grants to You a perpetual,
+      worldwide, non-exclusive, no-charge, royalty-free, irrevocable
+      copyright license to reproduce, prepare Derivative Works of,
+      publicly display, publicly perform, sublicense, and distribute the
+      Work and such Derivative Works in Source or Object form.
+
+   3. Grant of Patent License. Subject to the terms and conditions of
+      this License, each Contributor hereby grants to You a perpetual,
+      worldwide, non-exclusive, no-charge, royalty-free, irrevocable
+      (except as stated in this section) patent license to make, have made,
+      use, offer to sell, sell, import, and otherwise transfer the Work,
+      where such license applies only to those patent claims licensable
+      by such Contributor that are necessarily infringed by their
+      Contribution(s) alone or by combination of their Contribution(s)
+      with the Work to which such Contribution(s) was submitted. If You
+      institute patent litigation against any entity (including a
+      cross-claim or counterclaim in a lawsuit) alleging that the Work
+      or a Contribution incorporated within the Work constitutes direct
+      or contributory patent infringement, then any patent licenses
+      granted to You under this License for that Work shall terminate
+      as of the date such litigation is filed.
+
+   4. Redistribution. You may reproduce and distribute copies of the
+      Work or Derivative Works thereof in any medium, with or without
+      modifications, and in Source or Object form, provided that You
+      meet the following conditions:
+
+      (a) You must give any other recipients of the Work or
+          Derivative Works a copy of this License; and
+
+      (b) You must cause any modified files to carry prominent notices
+          stating that You changed the files; and
+
+      (c) You must retain, in the Source form of any Derivative Works
+          that You distribute, all copyright, patent, trademark, and
+          attribution notices from the Source form of the Work,
+          excluding those notices that do not pertain to any part of
+          the Derivative Works; and
+
+      (d) If the Work includes a "NOTICE" text file as part of its
+          distribution, then any Derivative Works that You distribute must
+          include a readable copy of the attribution notices contained
+          within such NOTICE file, excluding those notices that do not
+          pertain to any part of the Derivative Works, in at least one
+          of the following places: within a NOTICE text file distributed
+          as part of the Derivative Works; within the Source form or
+          documentation, if provided along with the Derivative Works; or,
+          within a display generated by the Derivative Works, if and
+          wherever such third-party notices normally appear. The contents
+          of the NOTICE file are for informational purposes only and
+          do not modify the License. You may add Your own attribution
+          notices within Derivative Works that You distribute, alongside
+          or as an addendum to the NOTICE text from the Work, provided
+          that such additional attribution notices cannot be construed
+          as modifying the License.
+
+      You may add Your own copyright statement to Your modifications and
+      may provide additional or different license terms and conditions
+      for use, reproduction, or distribution of Your modifications, or
+      for any such Derivative Works as a whole, provided Your use,
+      reproduction, and distribution of the Work otherwise complies with
+      the conditions stated in this License.
+
+   5. Submission of Contributions. Unless You explicitly state otherwise,
+      any Contribution intentionally submitted for inclusion in the Work
+      by You to the Licensor shall be under the terms and conditions of
+      this License, without any additional terms or conditions.
+      Notwithstanding the above, nothing herein shall supersede or modify
+      the terms of any separate license agreement you may have executed
+      with Licensor regarding such Contributions.
+
+   6. Trademarks. This License does not grant permission to use the trade
+      names, trademarks, service marks, or product names of the Licensor,
+      except as required for reasonable and customary use in describing the
+      origin of the Work and reproducing the content of the NOTICE file.
+
+   7. Disclaimer of Warranty. Unless required by applicable law or
+      agreed to in writing, Licensor provides the Work (and each
+      Contributor provides its Contributions) on an "AS IS" BASIS,
+      WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or
+      implied, including, without limitation, any warranties or conditions
+      of TITLE, NON-INFRINGEMENT, MERCHANTABILITY, or FITNESS FOR A
+      PARTICULAR PURPOSE. You are solely responsible for determining the
+      appropriateness of using or redistributing the Work and assume any
+      risks associated with Your exercise of permissions under this License.
+
+   8. Limitation of Liability. In no event and under no legal theory,
+      whether in tort (including negligence), contract, or otherwise,
+      unless required by applicable law (such as deliberate and grossly
+      negligent acts) or agreed to in writing, shall any Contributor be
+      liable to You for damages, including any direct, indirect, special,
+      incidental, or consequential damages of any character arising as a
+      result of this License or out of the use or inability to use the
+      Work (including but not limited to damages for loss of goodwill,
+      work stoppage, computer failure or malfunction, or any and all
+      other commercial damages or losses), even if such Contributor
+      has been advised of the possibility of such damages.
+
+   9. Accepting Warranty or Additional Liability. While redistributing
+      the Work or Derivative Works thereof, You may choose to offer,
+      and charge a fee for, acceptance of support, warranty, indemnity,
+      or other liability obligations and/or rights consistent with this
+      License. However, in accepting such obligations, You may act only
+      on Your own behalf and on Your sole responsibility, not on behalf
+      of any other Contributor, and only if You agree to indemnify,
+      defend, and hold each Contributor harmless for any liability
+      incurred by, or claims asserted against, such Contributor by reason
+      of your accepting any such warranty or additional liability.
+
+   END OF TERMS AND CONDITIONS
+
+   APPENDIX: How to apply the Apache License to your work.
+
+      To apply the Apache License to your work, attach the following
+      boilerplate notice, with the fields enclosed by brackets "[]"
+      replaced with your own identifying information. (Don't include
+      the brackets!)  The text should be enclosed in the appropriate
+      comment syntax for the file format. We also recommend that a
+      file or class name and description of purpose be included on the
+      same "printed page" as the copyright notice for easier
+      identification within third-party archives.
+
+   Copyright [yyyy] [name of copyright owner]
+
+   Licensed under the Apache License, Version 2.0 (the "License");
+   you may not use this file except in compliance with the License.
+   You may obtain a copy of the License at
+
+       http://www.apache.org/licenses/LICENSE-2.0
+
+   Unless required by applicable law or agreed to in writing, software
+   distributed under the License is distributed on an "AS IS" BASIS,
+   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+   See the License for the specific language governing permissions and
+   limitations under the License.

sqzcomb-0.1.1/PKG-INFO

@@ -0,0 +1,98 @@
+Metadata-Version: 2.4
+Name: sqzcomb
+Version: 0.1.1
+Summary: Squeezed light in Kerr microcombs: Lugiato-Lefever steady states, linearized quantum noise, and detected output spectra
+Author-email: "Tanvir M. Mahim" <tanvir.mahim@bracu.ac.bd>
+License: Apache-2.0
+Project-URL: Homepage, https://tanvir-mahmud-mahim.github.io/software/
+Project-URL: Repository, https://github.com/TaN-MM-Org/sqzcomb
+Project-URL: Issues, https://github.com/TaN-MM-Org/sqzcomb/issues
+Project-URL: Changelog, https://github.com/TaN-MM-Org/sqzcomb/releases
+Keywords: squeezed light,microcomb,Kerr,Lugiato-Lefever,quantum optics,homodyne
+Classifier: Development Status :: 3 - Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: Apache Software License
+Classifier: Programming Language :: Python :: 3
+Classifier: Topic :: Scientific/Engineering :: Physics
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy>=1.22
+Provides-Extra: test
+Requires-Dist: pytest; extra == "test"
+Dynamic: license-file
+
+# sqzcomb
+
+[![DOI](https://img.shields.io/badge/DOI-10.5281%2Fzenodo.22015375-blue)](https://doi.org/10.5281/zenodo.22015375) [![tests](https://github.com/TaN-MM-Org/sqzcomb/actions/workflows/ci.yml/badge.svg)](https://github.com/TaN-MM-Org/sqzcomb/actions)
+
+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.
+
+## Status
+
+v0.1.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
+
+Verified against closed forms in the test-suite: vacuum passes a passive
+cavity unchanged for every coupling 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.
+
+Not yet implemented (the v0.2+ roadmap, in order): the auxiliary-ring
+photonic molecule and its Purcell extraction of below-threshold modes,
+soliton-crystal steady-state continuation, supermode decomposition of the
+multimode covariance, and thermal input noise.
+
+## Install and use
+
+```
+pip install -e .
+```
+
+```python
+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.
+
+## 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, including the photonic-molecule
+extraction that this package's v0.1 does not yet contain.
+
+## License
+
+Apache-2.0

sqzcomb-0.1.1/README.md

@@ -0,0 +1,74 @@
+# sqzcomb
+
+[![DOI](https://img.shields.io/badge/DOI-10.5281%2Fzenodo.22015375-blue)](https://doi.org/10.5281/zenodo.22015375) [![tests](https://github.com/TaN-MM-Org/sqzcomb/actions/workflows/ci.yml/badge.svg)](https://github.com/TaN-MM-Org/sqzcomb/actions)
+
+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.
+
+## Status
+
+v0.1.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
+
+Verified against closed forms in the test-suite: vacuum passes a passive
+cavity unchanged for every coupling 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.
+
+Not yet implemented (the v0.2+ roadmap, in order): the auxiliary-ring
+photonic molecule and its Purcell extraction of below-threshold modes,
+soliton-crystal steady-state continuation, supermode decomposition of the
+multimode covariance, and thermal input noise.
+
+## Install and use
+
+```
+pip install -e .
+```
+
+```python
+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.
+
+## 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, including the photonic-molecule
+extraction that this package's v0.1 does not yet contain.
+
+## License
+
+Apache-2.0

sqzcomb-0.1.1/pyproject.toml

@@ -0,0 +1,33 @@
+[build-system]
+requires = ["setuptools>=68"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "sqzcomb"
+version = "0.1.1"
+description = "Squeezed light in Kerr microcombs: Lugiato-Lefever steady states, linearized quantum noise, and detected output spectra"
+readme = "README.md"
+license = { text = "Apache-2.0" }
+authors = [{ name = "Tanvir M. Mahim", email = "tanvir.mahim@bracu.ac.bd" }]
+requires-python = ">=3.9"
+keywords = ["squeezed light", "microcomb", "Kerr", "Lugiato-Lefever", "quantum optics", "homodyne"]
+dependencies = ["numpy>=1.22"]
+classifiers = [
+  "Development Status :: 3 - Alpha",
+  "Intended Audience :: Science/Research",
+  "License :: OSI Approved :: Apache Software License",
+  "Programming Language :: Python :: 3",
+  "Topic :: Scientific/Engineering :: Physics",
+]
+
+[project.optional-dependencies]
+test = ["pytest"]
+
+[project.urls]
+Homepage = "https://tanvir-mahmud-mahim.github.io/software/"
+Repository = "https://github.com/TaN-MM-Org/sqzcomb"
+Issues = "https://github.com/TaN-MM-Org/sqzcomb/issues"
+Changelog = "https://github.com/TaN-MM-Org/sqzcomb/releases"
+
+[tool.setuptools.packages.find]
+where = ["src"]

sqzcomb-0.1.1/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

sqzcomb-0.1.1/src/sqzcomb/__init__.py

@@ -0,0 +1,22 @@
+"""sqzcomb: squeezed light in Kerr microcombs, from steady state to detected noise.
+
+Solves the Lugiato-Lefever equation for the classical intracavity field,
+linearizes the quantum fluctuations around that steady state, and computes
+output quadrature-noise spectra through the standard input-output relations,
+so that a resonator design can be judged by the number a homodyne detector
+would actually report.
+
+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).
+"""
+from .lle import lle_evolve, homogeneous_steady_states
+from .linearize import fluctuation_matrix, single_mode_parametric
+from .spectra import output_quadrature_variance, squeezing_db
+
+__version__ = "0.1.1"
+__all__ = [
+    "lle_evolve", "homogeneous_steady_states",
+    "fluctuation_matrix", "single_mode_parametric",
+    "output_quadrature_variance", "squeezing_db",
+]

sqzcomb-0.1.1/src/sqzcomb/linearize.py

@@ -0,0 +1,68 @@
+"""Linearized fluctuation dynamics around an LLE steady state.
+
+Writing psi = psi_s + delta and keeping first order, the doubled vector
+z = (delta_a(k1..kM), delta_a*(k1..kM)) obeys  dz/dt = M z + inputs, with
+
+    M = [[ A,        B       ],
+         [ conj(B),  conj(A) ]]
+
+    A_{kk'} = [-(1 + i alpha) + i D(k)] delta_{kk'} + 2 i (|psi_s|^2)^(k - k')
+    B_{kk'} = i (psi_s^2)^(k + k')  (hat = Fourier component)
+
+Below threshold every eigenvalue of M has a negative real part; the module
+checks this and refuses to build spectra on an unstable state.
+"""
+from __future__ import annotations
+
+import numpy as np
+
+from .lle import _linear_symbol
+
+
+def fluctuation_matrix(psi_s, alpha, dispersion=(0.0,), modes=None):
+    """Doubled-basis drift matrix M around steady state psi_s.
+
+    modes : iterable of integer mode numbers to keep (default: all grid
+        modes). Returns (M, modes_array).
+    """
+    psi_s = np.asarray(psi_s, dtype=complex)
+    n = psi_s.size
+    kgrid = np.fft.fftfreq(n, d=1.0 / n).astype(int)
+    if modes is None:
+        modes = kgrid.copy()
+    modes = np.asarray(list(modes), dtype=int)
+    m = modes.size
+
+    f_abs2 = np.fft.fft(np.abs(psi_s) ** 2) / n     # components of |psi_s|^2
+    f_sq = np.fft.fft(psi_s ** 2) / n               # components of psi_s^2
+    comp_abs2 = dict(zip(kgrid, f_abs2))
+    comp_sq = dict(zip(kgrid, f_sq))
+
+    L = _linear_symbol(alpha, dispersion, modes.astype(float))
+    A = np.zeros((m, m), dtype=complex)
+    B = np.zeros((m, m), dtype=complex)
+    for i, ki in enumerate(modes):
+        for j, kj in enumerate(modes):
+            if i == j:
+                A[i, j] += L[i]
+            A[i, j] += 2j * comp_abs2.get(ki - kj, 0.0)
+            B[i, j] += 1j * comp_sq.get(ki + kj, 0.0)
+    M = np.block([[A, B], [np.conj(B), np.conj(A)]])
+    return M, modes
+
+
+def single_mode_parametric(mu, delta=0.0):
+    """Drift matrix of a single-mode below-threshold parametric process.
+
+    dz/dt = M z with M = [[-1 + i delta, mu], [conj(mu), -1 - i delta]].
+    Threshold at |mu| = 1 (for delta = 0). For real mu > 0 the phi = pi/2
+    output quadrature is squeezed and phi = 0 anti-squeezed. Used for
+    validation against the closed-form spectra of the degenerate
+    parametric oscillator.
+    """
+    return np.array([[-1.0 + 1j * delta, mu],
+                     [np.conj(mu), -1.0 - 1j * delta]], dtype=complex)
+
+
+def is_stable(M):
+    return bool(np.max(np.linalg.eigvals(M).real) < 0.0)

sqzcomb-0.1.1/src/sqzcomb/lle.py

@@ -0,0 +1,63 @@
+"""Lugiato-Lefever solver (normalized units).
+
+    d psi / dt = -(1 + i alpha) psi + i |psi|^2 psi + i sum_m (d_m/m!) (-i d/dtheta)^m psi + F
+
+with time in units of the photon lifetime 2/kappa, alpha the normalized
+detuning, F the normalized pump, and d_m the normalized dispersion
+coefficients (d_2 < 0 anomalous in this sign convention: D(k) = d2/2 k^2 + ...).
+
+Integrator: Strang-split step. The Kerr phase step is exact (|psi|^2 is
+invariant under it); the linear-plus-pump step is exact in the Fourier
+domain (first-order linear ODE with constant inhomogeneity).
+"""
+from __future__ import annotations
+
+import math
+
+import numpy as np
+
+
+def _linear_symbol(alpha, dispersion, k):
+    D = np.zeros_like(k, dtype=float)
+    for m, dm in enumerate(dispersion, start=2):
+        D += dm / math.factorial(m) * k ** m
+    return -(1.0 + 1j * alpha) + 1j * D
+
+
+def lle_evolve(psi0, F, alpha, dispersion=(0.0,), t_end=50.0, dt=0.01):
+    """Evolve the LLE from psi0 to t_end; returns the final field.
+
+    psi0 : complex array over theta in [0, 2 pi) (uniform grid).
+    F : scalar pump amplitude (homogeneous pump).
+    dispersion : (d2, d3, ...) normalized coefficients.
+    """
+    psi = np.asarray(psi0, dtype=complex).copy()
+    n = psi.size
+    k = np.fft.fftfreq(n, d=1.0 / n)  # integer mode numbers
+    L = _linear_symbol(alpha, dispersion, k)
+    eL = np.exp(L * dt)
+    Fk = np.zeros(n, dtype=complex)
+    Fk[0] = F * n  # FFT convention: homogeneous term lives in k = 0
+    # exact linear-with-pump propagator: psi_k -> eL psi_k + (eL - 1)/L * Fk
+    with np.errstate(divide="ignore", invalid="ignore"):
+        pump_prop = np.where(np.abs(L) > 1e-14, (eL - 1.0) / L, dt)
+    steps = int(round(t_end / dt))
+    for _ in range(steps):
+        psi *= np.exp(1j * np.abs(psi) ** 2 * (0.5 * dt))   # half Kerr
+        pk = np.fft.fft(psi)
+        pk = eL * pk + pump_prop * Fk                        # exact linear
+        psi = np.fft.ifft(pk)
+        psi *= np.exp(1j * np.abs(psi) ** 2 * (0.5 * dt))   # half Kerr
+    return psi
+
+
+def homogeneous_steady_states(F, alpha):
+    """Real intensities rho = |psi|^2 of the flat steady states.
+
+    Roots of  rho * (1 + (alpha - rho)^2) = F^2.  Returns the real,
+    non-negative roots in ascending order.
+    """
+    # rho^3 - 2 alpha rho^2 + (1 + alpha^2) rho - F^2 = 0
+    roots = np.roots([1.0, -2.0 * alpha, 1.0 + alpha * alpha, -float(F) ** 2])
+    real = roots[np.abs(roots.imag) < 1e-9].real
+    return np.sort(real[real >= 0.0])

sqzcomb-0.1.1/src/sqzcomb/spectra.py

@@ -0,0 +1,70 @@
+"""Output quadrature-noise spectra via input-output theory.
+
+The cavity couples to the extraction port with rate fraction eta = kappa_ex
+/ kappa and to intrinsic loss with 1 - eta (total damping normalized to 1
+in the doubled drift matrix M, i.e. Re eigenvalues -> -1 far from any
+parametric process). With vacuum at both inputs,
+
+    z(Omega)      = G(Omega) [ sqrt(2 eta) z_ex + sqrt(2 (1-eta)) z_0 ]
+    z_out(Omega)  = sqrt(2 eta) z(Omega) - z_ex(Omega)
+    G(Omega)      = (-i Omega I - M)^{-1}
+
+and the symmetrized output covariance follows from the vacuum input
+covariance. Quadrature X_phi = (a e^{-i phi} + a^dagger e^{i phi}) /
+sqrt(2), vacuum variance 1/2.
+"""
+from __future__ import annotations
+
+import numpy as np
+
+from .linearize import is_stable
+
+
+def _output_covariance(M, eta, omega):
+    m2 = M.shape[0]
+    ident = np.eye(m2, dtype=complex)
+    G = np.linalg.inv(-1j * omega * ident - M)
+    T_ex = 2.0 * eta * G - ident
+    T_0 = 2.0 * np.sqrt(eta * (1.0 - eta)) * G
+    # vacuum covariance <z_in z_in^dagger> in the (a, a*) doubled ordering:
+    # <a a^dagger> = 1, all other blocks zero
+    half = m2 // 2
+    N = np.zeros((m2, m2), dtype=complex)
+    N[:half, :half] = np.eye(half)
+    S = T_ex @ N @ T_ex.conj().T + T_0 @ N @ T_0.conj().T
+    return S
+
+
+def output_quadrature_variance(M, eta, omega, mode_index, n_modes,
+                               phi=0.0, mode_index_b=None):
+    """Symmetrized variance of an output quadrature at frequency omega.
+
+    mode_index : index of the mode (within the retained mode list) whose
+        quadrature is detected. If mode_index_b is given, the joint
+        two-mode quadrature (a + b)/sqrt(2) rotated by phi is used, the
+        natural variable for twin-beam squeezing.
+    Vacuum level is 0.5. Requires a stable M.
+    """
+    if not is_stable(M):
+        raise ValueError("drift matrix is unstable (above threshold); "
+                         "linearized spectra are meaningless there")
+    S = _output_covariance(M, eta, omega)
+    u = np.zeros(2 * n_modes, dtype=complex)
+    if mode_index_b is None:
+        u[mode_index] = np.exp(-1j * phi) / np.sqrt(2.0)
+        u[n_modes + mode_index] = np.exp(1j * phi) / np.sqrt(2.0)
+    else:
+        for idx, w in ((mode_index, 0.5), (mode_index_b, 0.5)):
+            u[idx] += np.exp(-1j * phi) * w
+            u[n_modes + idx] += np.exp(1j * phi) * w
+    # With N giving <a a^dagger> = 1 and zero elsewhere, u^dag S u is the
+    # output quadrature spectrum with vacuum level exactly 1/2; the
+    # passive-cavity identity (T_ex N T_ex^dag + T_0 N T_0^dag = N for any
+    # eta, Omega) is verified in the test-suite.
+    var = np.real(u.conj() @ S @ u)
+    return float(var)
+
+
+def squeezing_db(variance):
+    """Convert a quadrature variance to dB relative to vacuum (0.5)."""
+    return 10.0 * np.log10(variance / 0.5)

sqzcomb-0.1.1/src/sqzcomb.egg-info/PKG-INFO

@@ -0,0 +1,98 @@
+Metadata-Version: 2.4
+Name: sqzcomb
+Version: 0.1.1
+Summary: Squeezed light in Kerr microcombs: Lugiato-Lefever steady states, linearized quantum noise, and detected output spectra
+Author-email: "Tanvir M. Mahim" <tanvir.mahim@bracu.ac.bd>
+License: Apache-2.0
+Project-URL: Homepage, https://tanvir-mahmud-mahim.github.io/software/
+Project-URL: Repository, https://github.com/TaN-MM-Org/sqzcomb
+Project-URL: Issues, https://github.com/TaN-MM-Org/sqzcomb/issues
+Project-URL: Changelog, https://github.com/TaN-MM-Org/sqzcomb/releases
+Keywords: squeezed light,microcomb,Kerr,Lugiato-Lefever,quantum optics,homodyne
+Classifier: Development Status :: 3 - Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: Apache Software License
+Classifier: Programming Language :: Python :: 3
+Classifier: Topic :: Scientific/Engineering :: Physics
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy>=1.22
+Provides-Extra: test
+Requires-Dist: pytest; extra == "test"
+Dynamic: license-file
+
+# sqzcomb
+
+[![DOI](https://img.shields.io/badge/DOI-10.5281%2Fzenodo.22015375-blue)](https://doi.org/10.5281/zenodo.22015375) [![tests](https://github.com/TaN-MM-Org/sqzcomb/actions/workflows/ci.yml/badge.svg)](https://github.com/TaN-MM-Org/sqzcomb/actions)
+
+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.
+
+## Status
+
+v0.1.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
+
+Verified against closed forms in the test-suite: vacuum passes a passive
+cavity unchanged for every coupling 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.
+
+Not yet implemented (the v0.2+ roadmap, in order): the auxiliary-ring
+photonic molecule and its Purcell extraction of below-threshold modes,
+soliton-crystal steady-state continuation, supermode decomposition of the
+multimode covariance, and thermal input noise.
+
+## Install and use
+
+```
+pip install -e .
+```
+
+```python
+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.
+
+## 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, including the photonic-molecule
+extraction that this package's v0.1 does not yet contain.
+
+## License
+
+Apache-2.0

sqzcomb-0.1.1/src/sqzcomb.egg-info/SOURCES.txt

@@ -0,0 +1,13 @@
+LICENSE
+README.md
+pyproject.toml
+src/sqzcomb/__init__.py
+src/sqzcomb/linearize.py
+src/sqzcomb/lle.py
+src/sqzcomb/spectra.py
+src/sqzcomb.egg-info/PKG-INFO
+src/sqzcomb.egg-info/SOURCES.txt
+src/sqzcomb.egg-info/dependency_links.txt
+src/sqzcomb.egg-info/requires.txt
+src/sqzcomb.egg-info/top_level.txt
+tests/test_sqzcomb.py

sqzcomb-0.1.1/src/sqzcomb.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

sqzcomb-0.1.1/src/sqzcomb.egg-info/requires.txt

@@ -0,0 +1,4 @@
+numpy>=1.22
+
+[test]
+pytest

sqzcomb-0.1.1/src/sqzcomb.egg-info/top_level.txt

@@ -0,0 +1 @@
+sqzcomb

sqzcomb-0.1.1/tests/test_sqzcomb.py

@@ -0,0 +1,81 @@
+import numpy as np
+import pytest
+
+from sqzcomb import (lle_evolve, homogeneous_steady_states,
+                     fluctuation_matrix, single_mode_parametric,
+                     output_quadrature_variance, squeezing_db)
+
+
+def test_free_decay():
+    n = 64
+    psi0 = np.ones(n, dtype=complex)
+    psi = lle_evolve(psi0, F=0.0, alpha=0.0, t_end=3.0, dt=0.005)
+    assert np.allclose(np.abs(psi), np.exp(-3.0), rtol=1e-3)
+
+
+def test_flat_state_matches_cubic():
+    F, alpha = 1.1, 0.5
+    rho = homogeneous_steady_states(F, alpha)
+    # evolve from small seed to the (single) flat state
+    psi = lle_evolve(np.full(32, 0.1 + 0j), F=F, alpha=alpha,
+                     t_end=200.0, dt=0.01)
+    assert np.allclose(np.abs(psi) ** 2, rho[0], rtol=1e-4)
+    # residual of the cubic at the solver's intensity
+    r = np.abs(psi[0]) ** 2
+    assert np.isclose(r * (1 + (alpha - r) ** 2), F * F, rtol=1e-3)
+
+
+def test_vacuum_passes_through_passive_cavity():
+    M = np.diag([-1.0 + 0.3j, -1.0 - 0.3j]).astype(complex)
+    rng = np.random.default_rng(0)
+    for _ in range(10):
+        eta = rng.uniform(0.05, 1.0)
+        omega = rng.uniform(-3, 3)
+        phi = rng.uniform(0, np.pi)
+        v = output_quadrature_variance(M, eta, omega, mode_index=0,
+                                       n_modes=1, phi=phi)
+        assert np.isclose(v, 0.5, atol=1e-12)
+
+
+def test_matches_parametric_oscillator_closed_form():
+    # S_Y(Omega) = 1 - 4 eta mu / ((1 + mu)^2 + Omega^2), vacuum = 1
+    for mu in (0.3, 0.9):
+        M = single_mode_parametric(mu)
+        for eta in (0.3, 0.5, 1.0):
+            for omega in (0.0, 0.7, 2.5):
+                v = output_quadrature_variance(M, eta, omega, mode_index=0,
+                                               n_modes=1, phi=np.pi / 2)
+                analytic = 0.5 * (1 - 4 * eta * mu / ((1 + mu) ** 2 + omega ** 2))
+                assert np.isclose(v, analytic, rtol=1e-10), (mu, eta, omega)
+
+
+def test_three_db_extraction_limit_at_critical_coupling():
+    # eta = 1/2: as mu -> 1 the detectable squeezing saturates at 3 dB
+    M = single_mode_parametric(0.9999)
+    v = output_quadrature_variance(M, 0.5, 0.0, 0, 1, phi=np.pi / 2)
+    assert squeezing_db(v) > -3.011
+    assert squeezing_db(v) < -2.99
+    # full extraction breaks the limit (checked away from the mu -> 1
+    # numerical degeneracy: at mu = 0.99, eta = 1 the closed form gives
+    # about -46 dB)
+    M99 = single_mode_parametric(0.99)
+    v_full = output_quadrature_variance(M99, 1.0, 0.0, 0, 1, phi=np.pi / 2)
+    assert squeezing_db(v_full) < -20.0
+
+
+def test_above_threshold_rejected():
+    M = single_mode_parametric(1.2)
+    with pytest.raises(ValueError):
+        output_quadrature_variance(M, 0.5, 0.0, 0, 1)
+
+
+def test_lle_linearization_stable_below_threshold():
+    F, alpha = 0.8, 0.2
+    psi = lle_evolve(np.full(32, 0.05 + 0j), F=F, alpha=alpha,
+                     t_end=200.0, dt=0.01)
+    M, modes = fluctuation_matrix(psi, alpha)
+    assert np.max(np.linalg.eigvals(M).real) < 0.0
+    # and the machinery produces a finite spectrum on the pumped mode
+    i0 = int(np.where(modes == 0)[0][0])
+    v = output_quadrature_variance(M, 0.5, 0.0, i0, modes.size)
+    assert np.isfinite(v) and v > 0.0