Computational Physics
See recent articles
Showing new listings for Friday, 7 August 2026
- [1] arXiv:2608.05398 [pdf, html, other]
-
Title: Multiphysics tritium transport modelling of the ARC breeding blanket with FESTIMSubjects: Computational Physics (physics.comp-ph); Mathematical Physics (math-ph)
Accurate prediction of tritium behaviour in molten salt breeding blankets is essential for the design and safe operation of ARC-class fusion reactors. This work presents a fully open-source, component-scale multiphysics framework for modelling tritium transport in an ARC liquid immersion blanket. Neutron transport, thermal hydraulics, and hydrogen isotope transport are coupled using OpenMC, OpenFOAM, and FESTIM, leveraging dedicated tools enabling direct transfer of spatially resolved fields between solvers. Assuming a zero inlet concentration, steady-state simulations predict a total tritium inventory of approximately 243 mg, with the blanket reaching steady-state tritium throughput within approximately 30 min, which is of a similar order to previous system-level estimates. The results show that tritium transport is dominated by turbulence-enhanced diffusion, with strong localisation in flow stagnation regions and reduced accumulation in highly turbulent zones. Sensitivity analyses indicate that predicted inventories are governed primarily by the numerical stabilisation scheme, with only a modest dependence on the turbulent Schmidt number. The proposed workflow provides a transparent and extensible basis for high-fidelity analysis of tritium transport in ARC-class breeding blankets.
- [2] arXiv:2608.05477 [pdf, html, other]
-
Title: Data-driven reconstruction of dynamical systems using Takens' Theorem, manifold learning, and universal function approximatorsSubjects: Computational Physics (physics.comp-ph)
Embedding theorems can be used to provide theoretical guarantees about the relation between low-dimensional observations of a system and its full-dimensional state and dynamics. Such theorems do not, however, provide guidance on observable choice, embedding construction, or methodologies to learn the mapping between the embedding and full-dimensional state. In this work, we develop an algorithmic framework, TAkens Reconstruction (TAR), to analyze and reconstruct arbitrary dynamical systems from low-dimensional time series using an integration of Takens' Delay Embedding Theorem, manifold learning techniques, and universal function approximators. We validate TAR in applications to a variety of simulated and observed dynamical systems and use it to investigate how delay vector structure impacts reconstruction accuracy. In an ecological system, we show that simple predator-prey dynamics can be reconstructed with observations taken over a wide variety of embedding time scales. In molecular dynamics simulations of the protein Villin, we demonstrate how including multiple time delays of the same observable series can be used to improve reconstruction of systems with multiple characteristic time scales. In the trade record of Vanguard S&P 500, we show how the approach exposes underlying dynamical phenomenologies in the data and accurate return predictions over short time horizons without access to full-dimensional market observations. We develop and release an open-source software package to enable the application of TAR to arbitrary dynamical systems.
New submissions (showing 2 of 2 entries)
- [3] arXiv:2608.05320 (cross-list from cond-mat.mtrl-sci) [pdf, other]
-
Title: Nitrogen Vacancy Centers in Hexagonal Diamond Exhibit Long Coherence TimesComments: 6 pages, 4 figures, data and codes uploaded on QRESP after final peer review completeSubjects: Materials Science (cond-mat.mtrl-sci); Computational Physics (physics.comp-ph); Quantum Physics (quant-ph)
We show that negatively charged nitrogen-vacancy (NV) centers in the hexagonal diamond polymorph lonsdaleite offer a route to spin qubits with enhanced coherence relative to their cubic-diamond counterparts. Using first-principles calculations, we examine two distinct defect configurations, AA, with the same symmetry as in cubic diamond and AB, with reduced symmetry. We find that the AB configuration of the NV center exhibits a finite transverse zero-field splitting, giving rise to an approximate fourfold enhancement of the Hahn-echo coherence time $T_2$ at zero magnetic field. The AA configuration, by contrast, closely reproduces the electronic structure and coherence properties of the cubic NV center. We further characterize the many-body electronic structure, vertical excitation energies, and photoluminescence spectra of both configurations, providing spectral fingerprints for their experimental identification. Our results establish symmetry-broken NV centers in lonsdaleite as promising candidates for quantum sensing and information science applications.
- [4] arXiv:2608.05555 (cross-list from physics.plasm-ph) [pdf, html, other]
-
Title: Millisecond-Scale Neural Operator Surrogates for Double-Null Free-Boundary Grad-Shafranov EquilibriaPlamen G. Krastev (Harvard University)Comments: 13 pages, 8 figures, 3 tablesSubjects: Plasma Physics (physics.plasm-ph); Nuclear Experiment (nucl-ex); Computational Physics (physics.comp-ph)
The Grad-Shafranov (GS) equation governs ideal magnetohydrodynamic equilibrium in tokamak plasmas. Free-boundary GS solvers are central to diverted-equilibrium modeling, but nonlinear Picard iteration introduces computational cost and sample-dependent latency that can become prohibitive in optimization, modeling, and control-oriented loops. Here we train a geometrically conditioned Fourier Neural Operator (FNO) to learn a constrained forward map from spatial coordinates, scalar operating parameters $(P_{\mathrm{axis}}, I_p, f_{\mathrm{vac}})$, and prescribed X-point locations to the poloidal-flux field $\psi(R,Z)$. The model is trained on a controlled family of constrained double-null free-boundary equilibria generated with \textsc{FreeGS} for a single fixed machine geometry and prescribed topology. The best model achieves a mean relative $L^2$ error of $0.05\%$, with test error following an empirical $N^{-0.68}$ power law over $N_{\mathrm{train}}\in\{500,1000,2000,5000\}$. It recovers both X-points to within $0.2$ cm and localizes the O-point to $0.03$ cm. As a physics-consistency diagnostic, the predicted fields satisfy an external finite-difference GS residual evaluation at the same level as the ground-truth fields, with mean normalized residual $2.29$, indistinguishable from the $2.29\pm0.06$ \textsc{FreeGS} baseline using the same diagnostic. The trained FNO evaluates one equilibrium in $2.77$ ms on GPU and $25.6$ ms on CPU, corresponding to speedups of ${\sim}640\times$ and ${\sim}69\times$ relative to \textsc{FreeGS} as configured here, with near-deterministic latency (p95/median $=1.01$). These results show that neural-operator surrogates can provide accurate, geometrically precise, millisecond-scale equilibrium evaluations for magnetic-confinement fusion workflows within a prescribed topology and machine geometry.
- [5] arXiv:2608.05589 (cross-list from cs.CE) [pdf, other]
-
Title: A hybrid s-version isogeometric strategy for dynamic crack propagation in 2D and 3D problemsSubjects: Computational Engineering, Finance, and Science (cs.CE); Computational Physics (physics.comp-ph)
A hybrid s-version of isogeometric analysis (hS-IGA) strategy is proposed for accurate and efficient evaluation of near-crack fracture quantities in dynamic crack propagation analysis. The strategy retains the global-local superposition framework of the conventional s-method, while introducing B-spline basis functions only into the global discretisation and preserving a Lagrange-based local mesh in the crack domain. This hybrid formulation is motivated by the continuity-related bottleneck in the global-local coupling integration of the conventional Lagrange-based s-method, and by the need to retain a Lagrange-based local mesh for crack representation and post-processing of the dynamic stress intensity factor (DSIF) and local stress. The resulting formulation removes discontinuities in the coupling integrands caused by the global approximation and enables accurate coupling integration by standard Gauss quadrature without recursive subdivision. The proposed strategy is verified using two-dimensional stationary and dynamic straight-crack problems against the standard finite element method and the conventional s-method, and is further assessed using three-dimensional stationary and dynamically propagating circular-crack problems against the conventional s-method. Results show that the proposed hS-IGA strategy accurately evaluates the DSIF and local stress while retaining the global-local modelling advantages of the s-method. It also substantially reduces the number of integration points required for coupling integration, by approximately 81% in the two-dimensional dynamic benchmark and 95.6% in the three-dimensional dynamic benchmark relative to the conventional s-method. These results demonstrate that the proposed hS-IGA framework provides an accurate and efficient global-local strategy for dynamic crack propagation analyses requiring reliable evaluation of near-crack fracture quantities.
- [6] arXiv:2608.05682 (cross-list from math.OC) [pdf, other]
-
Title: A Damped Subspace Splitting Algorithm for Constrained Density Functional TheoryComments: 30 pages, 8 figures, 2 tablesSubjects: Optimization and Control (math.OC); Chemical Physics (physics.chem-ph); Computational Physics (physics.comp-ph)
Constrained density functional theory (CDFT) provides a powerful framework for describing electronically excited and charge-localized states, which underlie a broad range of physical and chemical phenomena. However, the discretized optimization problems arising from CDFT calculations remain challenging, owing to the presence of both the Stiefel manifold constraint and additional nonconvex quadratic constraints. Existing algorithms either fail to enforce the quadratic constraints with high accuracy or face convergence issues due to double-loop iterative structures. In this paper, we first derive a subspace-splitting reformulation that decouples the two groups of constraints, by exploiting the inherent rotation invariance and introducing a nonlinear subspace alignment constraint. Based on this reformulation, we propose a single-loop damped alternating direction method of multipliers, called DASSP. To the best of our knowledge, DASSP is the first algorithm for CDFT calculations with rigorous convergence guarantees. Each iteration of DASSP comprises a spectral minimization step, a projected gradient step, and a damped dual ascent step, all of which admit efficient implementations. Numerical results on synthetic and realistic CDFT problems demonstrate that DASSP attains high feasibility accuracy and exhibits favorable efficiency without compromising robustness. We expect that this work will pave the way toward reliable and efficient large-scale CDFT applications.
- [7] arXiv:2608.05749 (cross-list from physics.optics) [pdf, other]
-
Title: Transverse quantum-state characterization of programmable electron opticsShengbo You, Paolo Rosi, Enzo Rotunno, Alberto Roncaglia, Luca Belsito, Amir H. Tavabi, Rafal E. Dunin-Borkowski, Vincenzo Grillo, Philipp M. PelzComments: 18 pages, 6 figuresSubjects: Optics (physics.optics); Materials Science (cond-mat.mtrl-sci); Computational Physics (physics.comp-ph); Quantum Physics (quant-ph)
Programmable electron optics -- electronically controlled phase plates -- underpin proposals from dose-efficient phase imaging to shaped-electron X-ray sources, nearly all assuming a pure, fully coherent delivered wave whose purity has never been measured. Here we reconstruct the transverse density matrix of a microelectromechanical electrostatic spiral phase plate by mixed-state ptychography, from one four-dimensional STEM scan per state and without added hardware. The delivered beam is substantially mixed: its purity falls from approximately 0.47 to approximately 0.24 as the applied bias grows, inconsistent with a fixed lateral source-blur model, while the real-space coherence width stays near 1 nm. The same scans calibrate the device in situ, allow virtual orbital-angular-momentum sorting and, through a partial-coherence-aware transfer theory, indicate that purifying the output could improve dose efficiency roughly threefold. One acquisition thus becomes a quantum-state acceptance test for programmable electron optics, supplying the purity and coherence that emerging phase-plate and diffractive-imaging schemes assume but leave unquantified.
- [8] arXiv:2608.05768 (cross-list from math.NA) [pdf, html, other]
-
Title: Minimization-based polynomial corrections for high-order curved boundaries on fixed and moving domains: assessment on finite volume and discontinuous Galerkin schemesSubjects: Numerical Analysis (math.NA); Computational Physics (physics.comp-ph)
In this work, we present two novel strategies to impose high-order boundary conditions on fixed and moving curved domains, approximated with piecewise affine triangulations. Achieving high-order accuracy on curved domains requires tackling both the PDE discretization error and the geometrical error simultaneously. While the former can be reduced by employing high-order numerical methods such as finite volume and discontinuous Galerkin, the latter demands either a high-order parametrization of the physical domain or a consistent approximation of the boundary conditions. Minimization-based approaches like the Reconstruction for Off-site Data (ROD) method allow one to skip the construction of high-order curvilinear meshes by defining high-order consistent boundary conditions on a computational boundary that does not match the physical one. The ROD approach mitigates the second-order geometrical error by retrieving a modified polynomial in each boundary cell, which enforces the boundary conditions exactly on the physical boundary. However, the standard ROD method requires the inversion of a local linear system, whose cost grows with the polynomial degree and mesh refinement. Inspired by a recent one-dimensional analysis, we show that ROD-type approaches can be recast as simple polynomial corrections, applicable without any linear system inversion. This greatly simplifies the development of minimization-based boundary treatments and reduces the associated computational cost. To prove the wide applicability of our strategy, we develop it within a Runge-Kutta discontinuous Galerkin framework and an ADER arbitrary-Lagrangian-Eulerian finite volume framework for compressible flows on fixed and moving domains. Several numerical experiments with Dirichlet and slip-wall boundary conditions are presented, with convergence analysis up to fifth order in both 2D and 3D.
- [9] arXiv:2608.05918 (cross-list from cond-mat.mtrl-sci) [pdf, html, other]
-
Title: Jsymm: A Python package for symmetry analysis of exchange tensors in magnetic HamiltoniansComments: 29 pages, 1 figureSubjects: Materials Science (cond-mat.mtrl-sci); Computational Physics (physics.comp-ph)
Symmetries of a crystal often restrict its physical properties. In particular, they determine possible forms of the tensors that describe interatomic exchange interaction, which governs a wide range of magnetic phenomena. Computationally demanding first-principles calculations of the exchange tensors can be greatly simplified by taking the symmetry constraints into account. Here, we present Jsymm, a Python package that derives the most general symmetry-compatible form of the exchange tensors directly from the crystallographic data. For any bond formed by magnetic ions, Jsymm produces the tensors of the Dzyaloshinskii-Moriya and anisotropic Heisenberg exchange interaction in symbolic form, as well as the tensors for all other bonds related to it by symmetry. This reduces the number of independent model parameters, dramatically lowering the computational cost of the ab initio calculations and preventing unphysical results arising from symmetry violations. The package accepts standard CIF files and provides a web interface in addition to an interactive text mode and a Python library. We demonstrate its utility on La$_2$CuO$_4$ and $\alpha$-Fe$_2$O$_3$, reproducing known symmetry constraints and revealing additional relations between components of the exchange tensors of different bonds.
- [10] arXiv:2608.06139 (cross-list from cs.SD) [pdf, html, other]
-
Title: Explicit and Stable Pseudospectral Time-Domain Method for the Föppl-von Kármán EquationsComments: To be presented at Forum Acusticum 2026, Graz, Austria, September 2026Subjects: Sound (cs.SD); Computational Physics (physics.comp-ph)
Modal synthesis is a widely-used technique for simulation of musical instrument dynamics. In the linear case, a modal decomposition leads to an uncoupled system of damped and forced harmonic oscillators which can be efficiently solved by standard time-stepping methods. However, extensions to nonlinear problems are challenging due to the presence of products of modal expansions in the governing equations. In the case of the Föppl-von Kármán plate, the nonlinear coupling between the modes is described by a fourth-order tensor and is prohibitively expensive to evaluate in the modal domain. In this work, we propose a pseudospectral method in which the products are evaluated on a grid in the spatial domain while spatial derivatives are computed exactly in the modal domain. Discrete sine and cosine transforms between the modal and spatial domains are used to impose simply supported boundary conditions for the plate. Finally, we prove non-negativity of the nonlinear potential energy of the system and employ a scalar auxiliary variable technique for explicit and stable time integration in the modal domain. As a result, we reduce the computational cost of modal synthesis while preserving its advantages like a precise control over the simulated frequency range. Sound examples are presented.
- [11] arXiv:2608.06175 (cross-list from physics.soc-ph) [pdf, html, other]
-
Title: The Structure of Spreading on Temporal NetworksComments: 8 pages, 5 figuresSubjects: Physics and Society (physics.soc-ph); Mathematical Physics (math-ph); Adaptation and Self-Organizing Systems (nlin.AO); Applied Physics (physics.app-ph); Computational Physics (physics.comp-ph)
The physics of spreading in static networks is well understood through mappings to percolation. We show that spreading dynamics on temporal networks can analogously be mapped to reachability in temporal event graphs. This provides a theoretical and computational framework for a class of processes, such as variants of the susceptible-infected-susceptible model. Without explicit simulations, through the component analysis of event graphs, we obtain epidemic prevalence and derive epidemic thresholds for temporal networks with arbitrary degree and inter-event time distributions, with significant computational advantages as compared to explicit simulations.
- [12] arXiv:2608.06274 (cross-list from astro-ph.IM) [pdf, html, other]
-
Title: NestyNet. II. Coherent Function-Space Posteriors from Scientific Neural Surrogates (or How to Avoid Expensive MCMC)Comments: 16 pages 4 figuresSubjects: Instrumentation and Methods for Astrophysics (astro-ph.IM); Computational Physics (physics.comp-ph)
Scientific analyses increasingly use flexible neural networks, but their thousands of correlated parameters make it challenging to interpret the associated uncertainties. Here we develop a low-dimensional posterior for the fitted function itself, for scientific neural surrogates trained with second-order optimization. Linearizing the fitting procedure with respect to the randomized residual rows gives a measurement-to-function transport, the linear map, assembled from the converged Jacobians and Gauss--Newton curvature, that carries measurement perturbations into the function perturbations that refitting would produce. Its leading singular functions define coherent deformation modes. Independent Gaussian coefficients then generate smooth function draws, so that any derived quantity, including those requiring derivatives or integrals of the draw, inherits the posterior. The construction distinguishes repeated-experiment covariance from the local Gauss--Newton/Laplace posterior and propagates both to correlated quantities of scientific interest. The result is conditional on the fit's declared choices (architecture, hyperparameters, active set, and optimization branch), and every fit is certified as converged by checking that a further optimization step would change the fitted predictions by less than a chosen small fraction of the measurement errors. Our primary example is an 800-parameter phase-space distribution function fit for a mock stellar disk. Four uncertainty coordinates, two orders of magnitude fewer than the fitted parameters and stable under refinement of the force basis, capture 99\% of the vertical-force posterior variance, and 4000 coherent draws propagate through the force, total-density, surface-density, and frequency calculations in 0.8s. The method provides a highly efficient route to uncertainty propagation for derivative-dependent scientific inference.
Cross submissions (showing 10 of 10 entries)
- [13] arXiv:2604.14083 (replaced) [pdf, html, other]
-
Title: Distributional Inverse HomogenizationSubjects: Computational Physics (physics.comp-ph); Materials Science (cond-mat.mtrl-sci); Computation (stat.CO)
For many materials, macroscopic mechanical behavior is determined by an intricate microstructure. Understanding the relation between these two scales helps scientists and engineers design better materials. The relation which maps microstructure to bulk material properties can be understood via the well-established theory of homogenization. However inverting the homogenization process, to recover microstructural information from measured macroscopic properties, is fraught with difficulties because of the averaging processes that underlie homogenization. Therefore, scientists and engineers usually need recourse to more invasive, often highly localized, investigations to estimate the microstructure. In this work, we develop a noninvasive methodology by which one can leverage large collections of measured bulk material properties to infer information about the statistics of microstructure at a global level. We call this, distributional inverse homogenization. We study this problem in one and two dimensions, considering both periodic and stochastic homogenization. We demonstrate the methodology in the context of 2D Voronoi constructions and underpin the observed empirical success with theory in 1D. We also show how the natural spatial variability of microstructure can be exploited to gather data that enables distributional inversion. And we concurrently learn a surrogate model, approximating the homogenization map, that accelerates the resulting computations in this setting. The work formulates a new class of inverse problems, bridging ideas from probability and homogenization to facilitate the inference of microstructural material variability from macroscopic measurements.
- [14] arXiv:2506.03014 (replaced) [pdf, html, other]
-
Title: Convergence and efficiency proof of quantum imaginary time evolution for bounded order systemsComments: 23 pagesSubjects: Quantum Physics (quant-ph); Computational Complexity (cs.CC); Computational Physics (physics.comp-ph)
Many current and near-future applications of quantum computing utilise parametric families of quantum circuits and variational methods that can suffer from obstacles including non-convergence to the global minimum due to local minima, critical slowing down, or exponential resource scaling. Here we show that quantum imaginary time evolution can overcome these obstacles if the underlying physical system satisfies a set of conditions. This includes many relevant applications such as ground state preparation for local theories in physics or chemistry, combinatorial optimisation problems, or quantum machine learning. In particular, we analyse the quantum imaginary time evolution showing convergence guarantees to the global minimum without critical slowing down and providing a priori estimates on the required evolution time which scale linearly in system size and inverse energy gap. Furthermore, a provided complexity analysis shows that quantum imaginary time evolution can be efficiently compiled into a parametric quantum circuit, finding the optimal parameters included, for a large class of physically relevant problems.
- [15] arXiv:2601.17459 (replaced) [pdf, other]
-
Title: Qhronology: A Python package for studying quantum models of closed timelike curvesComments: 93 pages, 4 figuresSubjects: Quantum Physics (quant-ph); Computational Physics (physics.comp-ph)
Qhronology is a novel scientific-computing package for studying quantum models of closed timelike curves (CTCs) and simulating general quantum information processing and computation. Written in Python, the program provides a comprehensive framework for analyzing quantum theories of antichronological time travel, including functionality to calculate quantum resolutions to temporal paradoxes. It also operates as a complete quantum circuit simulator, enabling the examination of quantum algorithms and protocols in both numerical and symbolic capacities. In this paper, we formally introduce Qhronology, beginning with discussion on aspects of its design philosophy and architecture. An overview of its basic usage is then presented, along with a collection of examples demonstrating its various capabilities within a variety of distinct contexts. Lastly, the performance of the package's circuit simulation component is characterized by way of some simple empirical benchmarking.
- [16] arXiv:2604.24138 (replaced) [pdf, other]
-
Title: Synchronized molecular dynamics method for thin-layer flows of complex fluidsSubjects: Fluid Dynamics (physics.flu-dyn); Soft Condensed Matter (cond-mat.soft); Computational Physics (physics.comp-ph)
We propose a multiscale computational method for thin-layer flows of complex fluids, termed the synchronized molecular dynamics (SMD) method, which directly couples local molecular dynamics (MD) simulations with a macroscopic lubrication description. In thin layers, the flow can be decomposed into cross-sectional dynamics that are strongly influenced by interfacial effects, and streamwise transport along the channel. The SMD method exploits this separation of scales by sparsely distributing local MD cells along the channel and synchronizing them through macroscopic conservation laws.
In this framework, the macroscopic continuity equation is enforced by iteratively updating the external forces applied to each MD cell, thereby allowing the cross-sectional velocity profiles and the streamwise pressure distribution to be obtained without prescribing constitutive relations or boundary conditions. The method is validated for pressure-driven and wall-driven flows of Lennard--Jones fluids in a wedge-shaped channel, demonstrating excellent agreement with a modified Reynolds equation that accounts for boundary slip.
The SMD method is further applied to polymeric lubrication flows modeled by the Kremer--Grest chain model. At large pressure differences, the present approach naturally captures pronounced shear-thinning behavior coupled with microscopic polymer conformation dynamics. The results demonstrate that the SMD method provides an efficient and physically consistent framework for the multiscale simulation of complex fluid thin-layer flows. - [17] arXiv:2607.01346 (replaced) [pdf, html, other]
-
Title: Hybrid Two-Level Transport Method with Solution Decomposition in Macro and Micro ComponentsComments: 9 pages, 6 figures, 2 tablesSubjects: Numerical Analysis (math.NA); Computational Physics (physics.comp-ph)
This paper presents a new hybrid MC/deterministic method for solving the one-group steady-state Boltzmann transport equation based on decomposition of solution in macro and micro components. The macro component captures the large-scale structure of the solution. It is represented by angular moments of the high-order transport solution. The $P_1$ approximation is applied to define the macro component. The first two angular moments are obtained as a solution of hybrid low-order moment equations with exact closures. The equation for the micro component is solved using a MC simulation. The hybrid two-level system of equations for macro and micro components is solved by fixed-point iteration scheme. Numerical results are presented to demonstrate variance reduction of stochastic numerical solution and improvement in computational efficiency.
- [18] arXiv:2608.03759 (replaced) [pdf, html, other]
-
Title: AFLOW-EMERALD: ElectroMagnetic modes EngineeRing in Advanced LayereD materialsComments: 16 pages, 7 figuresSubjects: Optics (physics.optics); Computational Physics (physics.comp-ph)
Layered and periodically patterned heterostructures underpin advanced optical, photonic, and plasmonic (meta)materials, whose rational design demands electromagnetic solvers that are both numerically robust and tightly linked to the underlying material properties. Here, we present AFLOW-EMERALD (ElectroMagnetic modes EngineeRing in Advanced LayereD materials), an open-source, modular, Python-based computational framework for simulating electromagnetic wave propagation in finite and periodic layered (meta)materials. Built around a unified object-oriented architecture, AFLOW-EMERALD combines a numerically stable scattering-matrix method with plane-wave expansion and extends to rigorous coupled-wave analysis for laterally patterned structures such as gratings. The software computes optical spectra, spatial field distributions, and photonic band structures, including complex-k dispersion in lossy, dispersive media. A streamlined YAML workflow allows users to seamlessly import dielectric function datasets from experimental, literature, or first-principles sources. Owing to its modular design, AFLOW-EMERALD is readily extensible and suitable for integration into computational materials-design pipelines, providing a practical platform for the coupled material-geometry engineering of dielectric photonic crystals, plasmonic multilayers, hyperbolic metamaterials, and more complex architectures supporting, e.g., surface and volume plasmon-polariton modes.