Dataset and testbench: 4H-SiC photonic-molecule squeezed-light source

Complete raw simulation data and validation testbench underlying every figure and number in the article "Overcoming the 3 dB squeezing extraction limit in silicon carbide microcombs with a photonic molecule" (T. M. Mahim, M. M. Rahman, and A. S. M. Mohsin, 2026). The archive contains: (1) data/ - FEM effective-index sweeps and the selected-geometry solution with mesh-convergence table (femwell/scikit-fem, materials from the CC0 refractiveindex.info database), Lugiato-Lefever soliton-crystal states across the pump-detuning window and a third-order-dispersion scaling family, and the full quantum-model outputs: squeezing spectra, supermode decompositions, photonic-molecule Purcell-rate sweep, and the drop-port logarithmic-negativity matrix of the odd-mode lattice. (2) testbench/ - analytic validation of the quantum solver (vacuum output, Bogoliubov eigenvalues, 3 dB critical-coupling bound, Heisenberg uncertainty) and numerical convergence checks, with a frozen copy of the source code and reference output logs. All arrays are NumPy .npz/.json files; conventions and pass criteria are documented in the README inside the archive. Source code is maintained at https://github.com/Tanvir-Mahmud-Mahim/sic-molecule-squeezer (Apache-2.0). Key results: 8.5 dB detectable squeezing over 1.74 GHz at 8.3 mW pump; entangled 30-mode lattice with E_N up to 0.12; soliton crystal robust to 3.5x the design third-order dispersion. Version 1.1 (revision 1 of the article): adds the explicit two-ring photonic-molecule model and the third-order-dispersion stability-boundary study. New data: aux_ring.npz (two-ring resonance mismatch, heater alignment, mode-dependent Purcell rate for auxiliary-radius errors of 0/5/20 nm), q_tworing.npz (full non-adiabatic two-ring quantum model: extraction-rate and auxiliary-linewidth scan, adiabatic-convergence series, design-point spectra, supermodes, entanglement matrix), q_tworing_family.npz (full-model squeezing across the detuning family), d3_boundary.npz and d3_mech.npz (D3 boundary scan with step 0.25 and breathing-state diagnostics), plus the corresponding scripts and reference logs in the testbench. Revised key results (full two-ring model at kappa_aux/2pi = 8 GHz, J/2pi = 1.6 GHz): 7.9 dB detectable squeezing over 1.81 GHz at 8.3 mW pump, 8.5 dB near the annihilation boundary, E_N up to 0.13; the adiabatic limit remains 8.5-9.6 dB. Crystal survives 3.25x the design D3 and is lost at 3.5x via a dispersive-wave-driven breathing instability.