Singularity correction at metallic edges

Author: weiliangjin2021

CHANGELOG.md

@@ -14,6 +14,7 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 - 'ModeSimulation' now supports 'PermittivityMonitor'.
 - Classmethod `from_frequency_range` in `GaussianPulse` for generating a pulse whose amplitude in the frequency_range [fmin, fmax] is maximized, which is particularly useful for running broadband simulations.
 - Differentiable function `td.plugins.autograd.interpolate_spline` for 1D linear, quadratic, and cubic spline interpolation, supporting differentiation with respect to the interpolated values (`y_points`) and optional endpoint derivative constraints.
+- Enable singularity correction at PEC and LossyMetalMedium edges.
 
 ### Changed
 - Performance enhancement for adjoint gradient calculations by optimizing field interpolation.

tidy3d/components/subpixel_spec.py

@@ -120,6 +120,14 @@ class PECConformal(AbstractSubpixelAveragingMethod):
         ge=0,
     )
 
+    edge_singularity_correction: bool = pd.Field(
+        True,
+        title="Apply Singularity Model At Metal Edges",
+        description="Apply field correction model at metallic edges where field singularity occurs. "
+        "The correction is applied to straight edges aligned with the primal grids, and the wedge angle is either "
+        "0 or 90 degree.",
+    )
+
     @cached_property
     def courant_ratio(self) -> float:
         """The scaling ratio applied to Courant number so that the courant number