Warn when mode solver pml covers a significant portion of the mode plane

Author: caseyflex

Diff for: CHANGELOG.md

@@ -17,6 +17,7 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 - Bug in `LayerRefinementSpec` that refines grids outside the layer region when one in-plane dimension is of size infinity.
 - Querying tasks was sometimes erroring unexpectedly.
 - Fixed automatic creation of missing output directories.
+- Warn when mode solver pml covers a significant portion of the mode plane.
 
 ## [2.8.0] - 2025-03-04
 

Diff for: tests/test_components/test_simulation.py

@@ -10,6 +10,7 @@
 from tidy3d.components.scene import MAX_GEOMETRY_COUNT, MAX_NUM_MEDIUMS
 from tidy3d.components.simulation import MAX_NUM_SOURCES
 from tidy3d.exceptions import SetupError, Tidy3dError, Tidy3dKeyError
+from tidy3d.plugins.mode import ModeSolver
 
 from ..utils import (
     SIM_FULL,
@@ -2543,7 +2544,11 @@ def test_sim_subsection(unstructured, nz):
     sim_red = SIM_FULL.subsection(
         region=region, boundary_spec=td.BoundarySpec.all_sides(td.Periodic())
     )
-    sim_red = SIM_FULL.subsection(region=region, sources=[], grid_spec=td.GridSpec.uniform(dl=20))
+    sim_red = SIM_FULL.subsection(
+        region=region,
+        sources=[],
+        grid_spec=td.GridSpec.uniform(dl=20),
+    )
     assert len(sim_red.sources) == 0
     sim_red = SIM_FULL.subsection(region=region, monitors=[])
     assert len(sim_red.monitors) == 0
@@ -3085,6 +3090,81 @@ def test_validate_sources_monitors_in_bounds():
         )
 
 
+def test_mode_pml_warning():
+    sim_size = (3, 3, 3)
+    lambda0 = 1.55
+    freq0 = td.C_0 / lambda0
+    si = td.material_library["cSi"]["Li1993_293K"]
+    sio2 = td.material_library["SiO2"]["Horiba"]
+    wg = td.Structure(geometry=td.Box(size=(0.22, 0.5, td.inf)), medium=si)
+    mode_plane = td.Box(size=(2, 2, 0))
+    mode_spec = td.ModeSpec(num_pml=(22, 22))
+    grid_spec = td.GridSpec.auto(wavelength=lambda0, min_steps_per_wvl=30)
+    symmetry = (0, 0, 0)
+    with AssertLogLevel(None):
+        sim = td.Simulation(
+            size=sim_size,
+            medium=sio2,
+            structures=[wg],
+            grid_spec=grid_spec,
+            run_time=1e-30,
+            monitors=[
+                td.ModeSolverMonitor(
+                    size=(2, 2, 0),
+                    name="mode",
+                    freqs=[freq0],
+                    mode_spec=mode_spec.updated_copy(num_pml=(10, 10)),
+                )
+            ],
+            symmetry=symmetry,
+        )
+    with AssertLogLevel("WARNING", contains_str="covers more than"):
+        sim = td.Simulation(
+            size=sim_size,
+            medium=sio2,
+            structures=[wg],
+            grid_spec=grid_spec,
+            run_time=1e-30,
+            monitors=[
+                td.ModeSolverMonitor(
+                    size=(2, 2, 0), name="mode", freqs=[freq0], mode_spec=mode_spec
+                )
+            ],
+            symmetry=symmetry,
+        )
+    with AssertLogLevel("WARNING", contains_str="covers more than"):
+        sim = td.Simulation(
+            size=sim_size,
+            medium=sio2,
+            structures=[wg],
+            grid_spec=grid_spec,
+            run_time=1e-30,
+            sources=[
+                td.ModeSource(
+                    size=(2, 2, 0),
+                    direction="+",
+                    source_time=td.GaussianPulse(freq0=freq0, fwidth=0.1 * freq0),
+                    mode_spec=mode_spec,
+                )
+            ],
+            symmetry=symmetry,
+        )
+    with AssertLogLevel("WARNING", contains_str="covers more than"):
+        mode_solver = ModeSolver(
+            simulation=sim, plane=mode_plane, mode_spec=mode_spec, freqs=[freq0]
+        )
+    with AssertLogLevel("WARNING", contains_str="covers more than"):
+        mode_sim = td.ModeSimulation(
+            size=sim_size,
+            medium=sio2,
+            structures=[wg],
+            grid_spec=grid_spec,
+            plane=mode_plane,
+            mode_spec=mode_spec,
+            freqs=[freq0],
+        )
+
+
 def test_fixed_angle_sim():
     wvl_um = 1.0
     freq0 = td.C_0 / wvl_um