Skip to content

Anthracene inverse-design cavity: stopped at the Q-observability gate

Outcome

The Q≥25,000 target was not achieved. The inverse-design stage produced a nearly binary, feedthrough-connected Si3N4 structure in the requested stack, and both of its differentiable surrogate objectives improved. However, four independent direct-Q scans covering 16 gray or contour handoffs found zero certified resonances. Because there was no observable pole to track, the free-form boundary optimizer was not launched and no post-boundary geometry exists.

This is not a RAM failure or a hidden high-Q claim. Each final ringdown used a 35 nm fixed Cartesian mesh, double precision, 15 ps runtime, 12 probes, no early shutoff, and 6,800,640 full-domain cells before x/y symmetry reduction. The Tidy3D jobs completed normally. The failure is in cavity creation and objective alignment: the adjoint surrogates learned wavelength-selective scattering and local field enhancement without creating measurable stored energy.

The optimization ledger settled at 9.7851 of 12 FlexCredits, leaving 2.2149 FC unused. The separately authorized 2.5-FC final-evaluation ledger was never created because there was no fully Q-optimized frozen geometry to evaluate.

The cumulative cost stayed below 12 FC, while all four direct-Q gates certified zero modes

Physical stack and source

Final retained XZ and YZ cross-sections showing SiO2, patterned Si3N4 and air, 100 nm anthracene, PVA, and the embedded emitter

The simulated contract was:

  • a SiO2 substrate with n=1.45;
  • a 300 nm patterned Si3N4 layer with n=2.0 and air in its gaps and lateral voids;
  • a uniform 100 nm anthracene film with n=1.8;
  • an Ey dipole at (x, y, z) = (0, 0, 0.20) µm, exactly at the anthracene mid-plane;
  • PVA upper cladding with n=1.5;
  • no vacuum aperture; and
  • a 10 × 3 µm design region with 700 nm feedthrough guides and x/y symmetry.

The preflight sampled epsilon=3.24 at the source, exactly 1.8²; the dipole was not in air.

What inverse design produced

Matched-scale topology before optimization and after the final beta-64 local-field refinement

The weak seed was only a feedthrough guide and central pad. Projection continuation at beta=3, 8, 16, 32, and 64 formed a cavity-like distributed pattern and reduced the fraction of pixels between 0.1 and 0.9 from 84.67% to 6.57%. The final threshold-0.50 conversion yielded 84 exact contours with zero reported density-to-contour round-trip mismatch.

Adjoint objective histories, including the rejected high-learning-rate projection step and final emitter-field refinement

The initial mode-power contrast surrogate rose from 768.8 to 29,662.6 over the accepted continuation. A beta-32 learning rate of 0.02 collapsed the surrogate and was rejected; repeating at 0.002 recovered and improved it. The subsequent emitter-field contrast objective rose monotonically by 12.35%. These are real adjoint improvements, but neither metric is Q.

Paid stage Evaluated states / candidates Retained result Gray fraction Actual FC
beta=3 mode-power adjoint 24 768.8 → 10,670.8 84.67% 1.2000
beta=8 mode-power adjoint 20 650.7 → 8,893.9 44.87% 0.9750
beta=16 mode-power adjoint 16 5,447.2 → 19,855.4 25.33% 0.8000
beta=16 Q scan 5 0 certified modes 1.2074
beta=32 high-rate diagnostic 12 best state was index 0; later value 852.9 13.10% 0.6000
beta=32 corrected adjoint 16 16,710.4 → 26,562.4 13.40% 0.7250
beta=32 Q scan 5 0 certified modes 1.2074
beta=64 corrected adjoint 16 24,374.0 → 29,662.6 6.60% 0.7896
beta=64 Q scan 3 0 certified modes 0.7244
beta=64 emitter-field adjoint 8 1.2880e10 → 1.4471e10 6.57% 0.8320
final Q scan 3 0 certified modes 0.7244

Boundary handoff and the missing post-boundary result

Final adjoint density, its exact binary contour handoff, and the explicitly absent post-boundary geometry

The contour machinery worked: the final density became a signed-distance field and exact polygon set suitable for smooth normal motion. But a boundary-Q step requires a baseline pole whose frequency and decay can be tracked under positive and negative perturbations. With zero certified modes, any finite difference or QNM boundary direction would have been a fit-to-noise direction. Stopping here protected both the scientific claim and the remaining budget.

Field and ringdown evidence

Twelve anthracene-plane Ey probes show prompt excitation followed by a stationary numerical floor

The final point dipole launched substantial field: the largest sampled |Ey| was 742.97. After 2 ps, the collective RMS was 0.109626; after 10 ps it was still 0.109626. That unchanged late amplitude is a numerical floor, not an exponential cavity tail. Multi-window, multi-probe Tidy3D ResonanceFinder therefore returned zero certified modes at every threshold. A resonant spatial field map is intentionally absent because no resonant frequency was available at which such a map would be meaningful.

Diagnosis

The evidence rules out several tempting explanations:

  • Not the old source-placement error: the emitter is inside anthracene.
  • Not insufficient RAM: all 16 fixed-grid Q jobs completed, with no memory failure.
  • Not merely binarization loss: beta=16, beta=32, beta=64, and the final 6.57%-gray design all failed before and after thresholding.
  • Not a dead source or monitor: every one of the 12 probes received the prompt field.

The supported explanation is objective mismatch. Target-to-sideband output mode power can improve by interference, routing, and selective scattering. Even point-field contrast at three real frequencies can increase through a prompt near-field response. Neither objective directly rewards electromagnetic energy stored after the source is gone, so both can improve while ringdown Q remains unobservable.

The next paid campaign should not extend this seed with more boundary motion. It should first make the cavity-creation objective pass an observability gate:

  1. Validate the exact 35 nm, 12-probe ringdown protocol on a known resonant positive-control nanobeam in this same substrate/anthracene/PVA stack.
  2. Replace pure spectral contrast with a differentiable stored-energy/loss objective, such as omega * U_design / P_out, while retaining a minimum TE0 feedthrough-coupling constraint. Use a closed flux surface or explicit radiation/output channels so prompt field cannot masquerade as storage.
  3. Add a weak fixed Bragg/stop-band scaffold outside the inverse-design defect region. The adjoint still creates the cavity and coupling section, but it no longer has to discover both a bandgap and a defect from a uniform guide in one short budget.
  4. Run a short direct ringdown every continuation stage and retune the adjoint frequency to the measured pole. Do not continue projection unless at least one mode is consistent across probes and fit windows.
  5. Begin contour optimization only after a conservative baseline Q (for example Q>1,000) is independently observable. Then require direct positive and negative boundary controls before trusting each QNM-gradient step.

This turns ringdown from a final surprise into a campaign invariant.

Reproducibility

The retained campaign, exact simulation inputs, result summaries, raw SimulationData, credit decisions, and ledger are under engines/atom_fishbone_coupler/campaigns/anthracene_cavity_q25k_v1/. The visual manifest hashes every plotted source and output. The machine-readable campaign outcome is retained as campaign_outcome.json in the campaign directory.