Simulating Maxwell-Schr\"odinger Equations by High-Order Symplectic FDTD Algorithm

Guoda Xie; Ming Fang; Wei E.I. Sha; Zhixiang Huang

arxiv: 1906.09589 · v1 · pith:UDCFXJEXnew · submitted 2019-06-23 · ⚛️ physics.comp-ph · math.SG· physics.optics

Simulating Maxwell-Schr\"odinger Equations by High-Order Symplectic FDTD Algorithm

Guoda Xie , Zhixiang Huang , Ming Fang , Wei E.I. Sha This is my paper

Pith reviewed 2026-05-25 18:05 UTC · model grok-4.3

classification ⚛️ physics.comp-ph math.SGphysics.optics

keywords Maxwell-Schrödinger equationssymplectic integrationfinite-difference time-domainRabi oscillationlight-matter interactionnumerical methodsDirichlet boundary conditionshigh-order discretization

0 comments

The pith

A fourth-order symplectic FDTD algorithm simulates Maxwell-Schrödinger equations more accurately and efficiently than standard second-order methods for long-term Rabi oscillation studies.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper proposes a novel high-order symplectic finite-difference time-domain algorithm to solve the Maxwell-Schrödinger equations describing light-matter interactions. It discretizes the equations using fourth-order symplectic integration in time and fourth-order collocated finite differences in space. The Dirichlet boundary condition is enforced using an adapted image theory suitable for these high-order differences. Validation through three-dimensional simulations of population inversion in Rabi oscillations shows improved performance over the standard second-order FDTD approach in terms of accuracy and efficiency during extended time simulations. This establishes a method for energy-conserving modeling of coherent dynamics in such systems.

Core claim

The high-order SFDTD(4,4) algorithm, which applies fourth-order symplectic integration temporally and fourth-order collocated differences spatially to the Maxwell-Schrödinger equations, enables accurate and efficient simulation of coherent interactions between electromagnetic fields and artificial atoms under the semi-classical framework, as demonstrated by numerical studies of Rabi oscillations with population inversion.

What carries the argument

The symplectic finite-difference time-domain (SFDTD) algorithm using fourth-order symplectic integration and fourth-order collocated differences, with image theory adapted for Dirichlet boundaries.

Load-bearing premise

The adapted image theory successfully enforces the Dirichlet boundary condition for the fourth-order collocated differences without loss of accuracy or introduction of artifacts.

What would settle it

Running the Rabi oscillation test cases with the SFDTD(4,4) and FDTD(2,2) methods and observing that the high-order version does not show better accuracy or efficiency over long simulation times.

Figures

Figures reproduced from arXiv: 1906.09589 by Guoda Xie, Ming Fang, Wei E.I. Sha, Zhixiang Huang.

**Figure 3.** Figure 3: (b) The relative errors of the FDTD(2,2) approach and S [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

**Figure 4.** Figure 4: (b) The relative errors of the FDTD(2,2) approach and S [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗

**Figure 6.** Figure 6: (a) The calculation results for a large detuning Δ=0.3 [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗

**Figure 7.** Figure 7: The relative errors of the FDTD(2,2), FDTD(2,2)-DG and [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗

read the original abstract

A novel symplectic algorithm is proposed to solve the Maxwell-Schr\"odinger (M-S) system for investigating light-matter interaction. Using the fourth-order symplectic integration and fourth-order collocated differences, Maxwell-Schr\"odinger equations are discretized in temporal and spatial domains, respectively. The symplectic finite-difference time-domain (SFDTD) algorithm is developed for accurate and efficient study of coherent interaction between electromagnetic fields and artificial atoms. Particularly, the Dirichlet boundary condition is adopted for modeling the Rabi oscillation problems under the semi-classical framework. To implement the Dirichlet boundary condition, image theory is introduced, tailored to the high-order collocated differences. For validating the proposed SFDTD algorithm, three-dimensional numerical studies of the population inversion in the Rabi oscillation are presented. Numerical results show that the proposed high-order SFDTD(4,4) algorithm exhibits better numerical performance than the conventional FDTD(2,2) approach at the aspects of accuracy and efficiency for the long-term simulation. The proposed algorithm opens up a promising way towards a high-accurate energy-conservation modeling and simulation of complex dynamics in nanoscale light-matter interaction.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a concrete high-order SFDTD(4,4) discretization for the Maxwell-Schrödinger system with an adapted image method for Dirichlet boundaries, but the abstract supplies no error tables or convergence data to back the accuracy and efficiency claims.

read the letter

The main takeaway is a fourth-order symplectic time integrator combined with fourth-order collocated spatial differences for the Maxwell-Schrödinger equations, plus a custom image-theory treatment to enforce Dirichlet boundaries. That specific pairing for this coupled system is presented as the new piece. They apply it to three-dimensional Rabi oscillation tests and state that the scheme runs more accurately and efficiently than standard FDTD(2,2) over long times. The symplectic property is a sensible choice here because it should help keep the coherent dynamics stable without artificial energy drift. The boundary construction is described in enough detail that a reader could try to reproduce it. The soft spots are straightforward. The abstract asserts superior performance from the numerical studies yet shows none of the supporting numbers—no L2 errors, no observed convergence rates, no timing comparisons. Without those, the central claim stays unverified. The stress-test point about whether the image theory keeps the full fourth-order truncation error near the boundaries is also unresolved from what is given; if the stencil order drops at the boundary the long-term population-inversion results could be affected. This work is aimed at people who already run FDTD codes for nanophotonics or semi-classical quantum optics and want a higher-order option that preserves structure. A reader in that niche would get practical value from the boundary trick if the order is preserved. It deserves a serious referee because the method is explicit, the target equations are relevant, and the gaps are fixable with added tables and a short boundary-order analysis rather than a fundamental rewrite.

Referee Report

2 major / 0 minor

Summary. The manuscript proposes a novel high-order symplectic FDTD algorithm (SFDTD(4,4)) for the Maxwell-Schrödinger system, combining fourth-order symplectic time integration with fourth-order collocated spatial finite differences. Dirichlet boundary conditions for Rabi oscillation problems are enforced via an image theory construction tailored to the high-order stencils. Three-dimensional numerical experiments on population inversion are presented, with the central claim that SFDTD(4,4) outperforms conventional FDTD(2,2) in accuracy and efficiency for long-term simulations.

Significance. If the performance claims hold with the reported convergence and stability properties, the algorithm would provide a practical advance for energy-conserving, long-time modeling of coherent light-matter interactions at the nanoscale, where conventional low-order FDTD schemes suffer from dispersion and dissipation errors.

major comments (2)

[Section 3 (boundary conditions)] The adaptation of image theory to fourth-order collocated differences (described in the boundary-condition section, referenced as Section 3 in the stress-test note) must be shown to retain formal O(h^4) truncation error at the Dirichlet boundaries; if the construction only satisfies the boundary condition to lower order, the observed long-term accuracy gains in the Rabi population-inversion tests would not be generalizable.
[Numerical results / validation section] The numerical-results section supplies no quantitative error metrics (e.g., L2 norms of field or population errors), observed convergence rates, or tabulated comparisons against FDTD(2,2); without these data the headline claim of superior accuracy and efficiency cannot be verified from the manuscript itself.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments, which help clarify the presentation of our high-order symplectic FDTD algorithm. We address each major point below and will revise the manuscript accordingly to strengthen the validation and analysis.

read point-by-point responses

Referee: [Section 3 (boundary conditions)] The adaptation of image theory to fourth-order collocated differences (described in the boundary-condition section, referenced as Section 3 in the stress-test note) must be shown to retain formal O(h^4) truncation error at the Dirichlet boundaries; if the construction only satisfies the boundary condition to lower order, the observed long-term accuracy gains in the Rabi population-inversion tests would not be generalizable.

Authors: We agree that a formal demonstration of the truncation error order at the boundaries is essential for generalizability. The image theory construction enforces the Dirichlet condition exactly at boundary nodes and sets ghost values to preserve the fourth-order interior stencil consistency. However, the manuscript does not include an explicit local truncation error expansion near the boundaries. In the revision we will add this analysis in Section 3, together with supporting numerical convergence tests localized near the boundaries, to confirm that the O(h^4) order is retained. revision: yes
Referee: [Numerical results / validation section] The numerical-results section supplies no quantitative error metrics (e.g., L2 norms of field or population errors), observed convergence rates, or tabulated comparisons against FDTD(2,2); without these data the headline claim of superior accuracy and efficiency cannot be verified from the manuscript itself.

Authors: We acknowledge that the current numerical section relies on visual inspection of figures without tabulated quantitative metrics, making independent verification of the accuracy and efficiency claims difficult. In the revised manuscript we will add L2-norm error tables for both fields and populations, observed convergence rates under successive grid refinement, and direct side-by-side comparisons (error versus CPU time) between SFDTD(4,4) and FDTD(2,2) for the Rabi-oscillation tests. revision: yes

Circularity Check

0 steps flagged

No circularity; derivation is self-contained from standard numerical methods

full rationale

The paper constructs the SFDTD(4,4) scheme directly from fourth-order symplectic time integrators and fourth-order collocated spatial differences applied to the Maxwell-Schrödinger system, then adapts image theory for Dirichlet boundaries and validates via direct numerical experiments on Rabi oscillations. No equation, result, or performance claim reduces by construction to a fitted parameter, self-definition, or self-citation chain. The reported accuracy/efficiency gains are outputs of the implemented algorithm, not inputs renamed as predictions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The contribution is a numerical scheme built on standard mathematical properties of symplectic integrators and finite-difference approximations; the only notable domain assumption is the suitability of the semi-classical Maxwell-Schrödinger model for the chosen test problem.

axioms (1)

domain assumption The semi-classical Maxwell-Schrödinger framework is appropriate for modeling Rabi oscillations in artificial atoms under Dirichlet boundaries.
The paper adopts this framework for the validation studies described in the abstract.

pith-pipeline@v0.9.0 · 5743 in / 1260 out tokens · 39937 ms · 2026-05-25T18:05:09.188921+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/ArithmeticFromLogic.lean LogicNat initial Peano object; symplectic structure from Laws of Logic echoes

?

echoes
ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

the matrix U can be represented as ... U^T J + J U = 0 ... U is an infinitesimal real symplectic matrix ... the symplectic algorithm can be employed
IndisputableMonolith/Foundation/ArrowOfTime.lean Berry-phase monotonicity and Z-monotonicity for directed time echoes

?

echoes
ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

symplectic integrators ... energy-conserving and highly stable characteristics ... long-term simulation

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages

[1]

J, Theory and Computational of Electromagnetic Fields

J.-M. J, Theory and Computational of Electromagnetic Fields . Hoboken, NJ, USA: Wiley, 2010

work page 2010
[2]

Full-quantum numerical sche me of finite difference time domain method for high-order harmonic generation,

L. Zhang, H. Zeng, and R. S. Chen, “Full-quantum numerical sche me of finite difference time domain method for high-order harmonic generation,” IEEE J. Multiscale Multiphys. Comput. Techn ., vol. 3, pp. 74–79, Aug. 2018

work page 2018
[3]

A new 3-D transmis sion line matrix scheme for the combined Schrödinger–Maxwell problem in t he electronic/electromagnetic characterization of nanodevices,

L. Pierantoni, D. Mencarelli, and T. Rozzi, “A new 3-D transmis sion line matrix scheme for the combined Schrödinger–Maxwell problem in t he electronic/electromagnetic characterization of nanodevices,” IEEE Trans. Microw. Theory Techn., vol. 56, no. 3, pp. 654-662, Mar. 2008

work page 2008
[4]

Boundary immittanc e operators for the Schrödinger–Maxwell problem of carrier dynami cs in nanodevices,

L. Pierantoni, D. Mencarelli, and T. Rozzi, “Boundary immittanc e operators for the Schrödinger–Maxwell problem of carrier dynami cs in nanodevices,” IEEE Trans. Microw. Theory Techn ., vol. 57, no. 5, pp. 1147-1155, May 2009

work page 2009
[5]

Numerical solution of initial boundary value proble ms involving Maxwell’s equations in isotropic media,

K. S. Yee, “Numerical solution of initial boundary value proble ms involving Maxwell’s equations in isotropic media,” IEEE Trans. Antenna Propag., vol. 14, pp. 302-307, May 1966

work page 1966
[6]

Taflove and S

A. Taflove and S. C. Hagness, Computational electrodynamics the finite difference time-domain method (3rd edition), Artech House Boston London, 2005

work page 2005
[7]

Schrödinger electron s interacting with optical grat ings: quantum mechanical study of the inverse Smith–Purcell effect,

N. Talebi, “Schrödinger electron s interacting with optical grat ings: quantum mechanical study of the inverse Smith–Purcell effect,” New J. Phys., vol. 18, no.12, pp. 123006, Dec. 2016

work page 2016
[8]

Nonlinear nanopolaritonics: Finite-difference time-domain Maxwell–Schrödinger simulation of molecule-assisted plasmon transfer,

K. Lopata and D. Neuhauser, “Nonlinear nanopolaritonics: Finite-difference time-domain Maxwell–Schrödinger simulation of molecule-assisted plasmon transfer,” J. Chem Phys ., vol. 131, no. 1, pp. 014701, Jul. 2009

work page 2009
[9]

Derivation of co upled Maxwell–Schrödinger equatio ns describing matter-laser interac tion from first principles of qu antum electrodynamics,

M. Šindelka, “Derivation of co upled Maxwell–Schrödinger equatio ns describing matter-laser interac tion from first principles of qu antum electrodynamics,” Phys. Rev. A., vol. 81, no. 3, pp. 033833, Mar. 2010

work page 2010
[10]

D y n a m i c o p t i c a l r e s p o n s e o f a n e x c i t o nic quantum dot studied by solving the self consistent Maxwell-Schr ödinger equations nonperturbatively,

S . H e l l s t r ö m a n d Y . F u , “ D y n a m i c o p t i c a l r e s p o n s e o f a n e x c i t o nic quantum dot studied by solving the self consistent Maxwell-Schr ödinger equations nonperturbatively,” Phys. Rev. B, Condens. Matter., vol. 82, no. 24, pp. 245305, 2010

work page 2010
[11]

A Maxwell–Schrödinger solver for quantum optical few-level systems,

R. Fleischhaker and J. Evers, “A Maxwell–Schrödinger solver for quantum optical few-level systems,” Comp. Phys. Commun., vol. 182, no. 3, pp. 739-747, Mar. 2011

work page 2011
[12]

Maxwell–Schrödinger hybri d simulation for optically controlling quantum states: A scheme f or designing control pulses,

T. Takeuchi, S. Ohnuki, and T. Sake, “Maxwell–Schrödinger hybri d simulation for optically controlling quantum states: A scheme f or designing control pulses,” Phys. Rev. A., vol. 91, no. 3, pp. 033401, Mar. 2015

work page 2015
[13]

Hybrid simulation of Maxwell-Schrödinger equations for multi-physics problems characterized by anharmonic electrostatic potential,

T. Takeuchi, S. Ohnuki, and T. Sake, “Hybrid simulation of Maxwell-Schrödinger equations for multi-physics problems characterized by anharmonic electrostatic potential,” Prog. Electromagn. Res., vol. 148, pp. 73-82, 2014

work page 2014
[14]

Solution of the electric fiel d integral equation when it breaks down,

J. F. Zhu, S. Omar, and D. Jiao, “Solution of the electric fiel d integral equation when it breaks down,” IEEE Trans. Antenna Propag., vol. 62, no. 8, pp. 4122-4134, Aug. 2014

work page 2014
[15]

E F I E A n a l y s i s o f L o w - F r e q u e ncy Problems With Loop-Star Decomposition and Calderón Multiplicati ve Preconditioner,

S . Y a n , J . M . J i n , a n d Z . P . N i e , “ E F I E A n a l y s i s o f L o w - F r e q u e ncy Problems With Loop-Star Decomposition and Calderón Multiplicati ve Preconditioner,” IEEE Trans. Antenna Propag ., vol. 58, no. 3, pp. 857-867, Mar. 2010

work page 2010
[16]

Finite-Dif ference Time-Domain simulation of the Maxwell–Schrödinger system,

C. J. Ryu, A. Y. Liu, W. E. I. Sha, and W. C. Chew, “Finite-Dif ference Time-Domain simulation of the Maxwell–Schrödinger system,” IEEE J. Multiscale Multiphys. Computat. Techn., vol. 1, pp. 40-47, 2016

work page 2016
[17]

A high-order symplectic FDTD scheme for the Maxwell-Schrodinger system

X. Chen, F. M. Kong, K. Li, and M. Liu, “A high-order symplectic FDTD scheme for the Maxwell-Schrodinger system” IEEE J. Quantum Electron., vol. 54, no. 1, pp. 6100108, Feb. 2018

work page 2018
[18]

A unified Hamiltonian solution to Maxwell–Schrödinger equation s for modeling electromagnetic fie ld-particle interaction,

Y. P. Chen, W. E.I. Sha, L. Jiang, M. Meng, Y. M. Wu, and W. C. Chew, “A unified Hamiltonian solution to Maxwell–Schrödinger equation s for modeling electromagnetic fie ld-particle interaction,” Comput. Phys. Commuun., vol. 215, pp. 63-70, Jun. 2017

work page 2017
[19]

Krishnamurthy, K

V. Krishnamurthy, K. Ravi, and S.T. Ho, “Analytical framework f or the s t e a d y s t a t e a n a l y s i s o f w a v e l e ngth-dependent and intensity-dep endent interaction of multiple monochromatic beams in semiconductor ph otonic structures with multiple ac tive and passive sections,” IEEE J. Quantum Electron., vol. 48, no. 10, pp. 1282–1299, Oct. 2012

work page 2012
[20]

Designing optimal universal pulses using second-order, large-scale, non-linear optimization,

C. K. Anand, A. D. Bain, A. T. Curtis, and Z. H. Nie, “Designing optimal universal pulses using second-order, large-scale, non-linear optimization,” J. Magn. Reson., vol, 219, pp. 61-74, Apr. 2012

work page 2012
[21]

Application of the sympl ectic finite-difference time-domain scheme to electromagnetic simulat ion,

W. Sha, Z. Huang, X. Wu, and M. Chen, “Application of the sympl ectic finite-difference time-domain scheme to electromagnetic simulat ion,” J. Comput. Phys., vol. 225, no. 1, pp. 33-50. Jan. 2007

work page 2007
[22]

Dispersion relation equation preserving FDTD method for nonlinear cubic Schrödinger equation,

T. W. H. Sheu and L. Lin, “Dispersion relation equation preserving FDTD method for nonlinear cubic Schrödinger equation,” J. Comput. Phys., vol. 299, no. 4, pp. 1-21. Oct. 2015

work page 2015
[23]

Canonical symplectic structure and structure-preserving geomet ric algorithms for Schrödinger-Maxwell systems,

Q. Chen, H. Qin, J. Liu, J. Y. Xiao, R. L. Zhang, Y. He, and Y. L. Wang, “Canonical symplectic structure and structure-preserving geomet ric algorithms for Schrödinger-Maxwell systems,” J. Comput. Phys., vol. 349, pp. 441-452, 2017

work page 2017
[24]

Symmetric and symplectic ERKN for oscillatory Hamiltonian systems,

Z. Chen, X. You, W. Shi, and Z. Liu, “Symmetric and symplectic ERKN for oscillatory Hamiltonian systems,” Comput. Phys. Commun., vol. 183, no.1, pp. 86-98, 2012

work page 2012
[25]

Multiple image technique (MIT) a nd anisotropic perfectly matched layer (APML) in implementation of MRTD scheme for boundary truncations of microwave structures,

Q. Cao, Y. Chen, and R. Mittra, “Multiple image technique (MIT) a nd anisotropic perfectly matched layer (APML) in implementation of MRTD scheme for boundary truncations of microwave structures,” IEEE Trans. Microw. Theory Techn., vol. 50, no. 6, pp. 1578-1589, Jun. 2002

work page 2002
[26]

Kira and S

M. Kira and S. W. Koch, Semiconductor Quantum Optics . Cambridge Univ. Press, 2012

work page 2012
[27]

A new solution of Schrödinger equation based on symplectic algorithm,

Z. Huang, J. Xu, B. Sun, B. Wu, and X. Wu, “A new solution of Schrödinger equation based on symplectic algorithm,” Comput. Math. Applicat., vol. 69, no. 11, pp.1303-1312, Jun. 2015

work page 2015
[28]

A symplectic FDTD algorit hm for the simulation of lossy dispersive materials,

H. wang, B. Wu, Z. Huang, and X. Wu, “A symplectic FDTD algorit hm for the simulation of lossy dispersive materials,” Comput. Phys. Commun., vol. 185, no.3, pp. 862-872, 2014

work page 2014
[29]

Gerry and P

C. Gerry and P. Knight, Introduc tory Quantum Optics, Cambride University Press, New York, 2005

work page 2005
[30]

Marlon. O. Scully, and M. Suhail. Zubairy, Quantum Optics. Cambridge, U.K.:Cambridge Univ. Press, 1997

work page 1997

[1] [1]

J, Theory and Computational of Electromagnetic Fields

J.-M. J, Theory and Computational of Electromagnetic Fields . Hoboken, NJ, USA: Wiley, 2010

work page 2010

[2] [2]

Full-quantum numerical sche me of finite difference time domain method for high-order harmonic generation,

L. Zhang, H. Zeng, and R. S. Chen, “Full-quantum numerical sche me of finite difference time domain method for high-order harmonic generation,” IEEE J. Multiscale Multiphys. Comput. Techn ., vol. 3, pp. 74–79, Aug. 2018

work page 2018

[3] [3]

A new 3-D transmis sion line matrix scheme for the combined Schrödinger–Maxwell problem in t he electronic/electromagnetic characterization of nanodevices,

L. Pierantoni, D. Mencarelli, and T. Rozzi, “A new 3-D transmis sion line matrix scheme for the combined Schrödinger–Maxwell problem in t he electronic/electromagnetic characterization of nanodevices,” IEEE Trans. Microw. Theory Techn., vol. 56, no. 3, pp. 654-662, Mar. 2008

work page 2008

[4] [4]

Boundary immittanc e operators for the Schrödinger–Maxwell problem of carrier dynami cs in nanodevices,

L. Pierantoni, D. Mencarelli, and T. Rozzi, “Boundary immittanc e operators for the Schrödinger–Maxwell problem of carrier dynami cs in nanodevices,” IEEE Trans. Microw. Theory Techn ., vol. 57, no. 5, pp. 1147-1155, May 2009

work page 2009

[5] [5]

Numerical solution of initial boundary value proble ms involving Maxwell’s equations in isotropic media,

K. S. Yee, “Numerical solution of initial boundary value proble ms involving Maxwell’s equations in isotropic media,” IEEE Trans. Antenna Propag., vol. 14, pp. 302-307, May 1966

work page 1966

[6] [6]

Taflove and S

A. Taflove and S. C. Hagness, Computational electrodynamics the finite difference time-domain method (3rd edition), Artech House Boston London, 2005

work page 2005

[7] [7]

Schrödinger electron s interacting with optical grat ings: quantum mechanical study of the inverse Smith–Purcell effect,

N. Talebi, “Schrödinger electron s interacting with optical grat ings: quantum mechanical study of the inverse Smith–Purcell effect,” New J. Phys., vol. 18, no.12, pp. 123006, Dec. 2016

work page 2016

[8] [8]

Nonlinear nanopolaritonics: Finite-difference time-domain Maxwell–Schrödinger simulation of molecule-assisted plasmon transfer,

K. Lopata and D. Neuhauser, “Nonlinear nanopolaritonics: Finite-difference time-domain Maxwell–Schrödinger simulation of molecule-assisted plasmon transfer,” J. Chem Phys ., vol. 131, no. 1, pp. 014701, Jul. 2009

work page 2009

[9] [9]

Derivation of co upled Maxwell–Schrödinger equatio ns describing matter-laser interac tion from first principles of qu antum electrodynamics,

M. Šindelka, “Derivation of co upled Maxwell–Schrödinger equatio ns describing matter-laser interac tion from first principles of qu antum electrodynamics,” Phys. Rev. A., vol. 81, no. 3, pp. 033833, Mar. 2010

work page 2010

[10] [10]

D y n a m i c o p t i c a l r e s p o n s e o f a n e x c i t o nic quantum dot studied by solving the self consistent Maxwell-Schr ödinger equations nonperturbatively,

S . H e l l s t r ö m a n d Y . F u , “ D y n a m i c o p t i c a l r e s p o n s e o f a n e x c i t o nic quantum dot studied by solving the self consistent Maxwell-Schr ödinger equations nonperturbatively,” Phys. Rev. B, Condens. Matter., vol. 82, no. 24, pp. 245305, 2010

work page 2010

[11] [11]

A Maxwell–Schrödinger solver for quantum optical few-level systems,

R. Fleischhaker and J. Evers, “A Maxwell–Schrödinger solver for quantum optical few-level systems,” Comp. Phys. Commun., vol. 182, no. 3, pp. 739-747, Mar. 2011

work page 2011

[12] [12]

Maxwell–Schrödinger hybri d simulation for optically controlling quantum states: A scheme f or designing control pulses,

T. Takeuchi, S. Ohnuki, and T. Sake, “Maxwell–Schrödinger hybri d simulation for optically controlling quantum states: A scheme f or designing control pulses,” Phys. Rev. A., vol. 91, no. 3, pp. 033401, Mar. 2015

work page 2015

[13] [13]

Hybrid simulation of Maxwell-Schrödinger equations for multi-physics problems characterized by anharmonic electrostatic potential,

T. Takeuchi, S. Ohnuki, and T. Sake, “Hybrid simulation of Maxwell-Schrödinger equations for multi-physics problems characterized by anharmonic electrostatic potential,” Prog. Electromagn. Res., vol. 148, pp. 73-82, 2014

work page 2014

[14] [14]

Solution of the electric fiel d integral equation when it breaks down,

J. F. Zhu, S. Omar, and D. Jiao, “Solution of the electric fiel d integral equation when it breaks down,” IEEE Trans. Antenna Propag., vol. 62, no. 8, pp. 4122-4134, Aug. 2014

work page 2014

[15] [15]

E F I E A n a l y s i s o f L o w - F r e q u e ncy Problems With Loop-Star Decomposition and Calderón Multiplicati ve Preconditioner,

S . Y a n , J . M . J i n , a n d Z . P . N i e , “ E F I E A n a l y s i s o f L o w - F r e q u e ncy Problems With Loop-Star Decomposition and Calderón Multiplicati ve Preconditioner,” IEEE Trans. Antenna Propag ., vol. 58, no. 3, pp. 857-867, Mar. 2010

work page 2010

[16] [16]

Finite-Dif ference Time-Domain simulation of the Maxwell–Schrödinger system,

C. J. Ryu, A. Y. Liu, W. E. I. Sha, and W. C. Chew, “Finite-Dif ference Time-Domain simulation of the Maxwell–Schrödinger system,” IEEE J. Multiscale Multiphys. Computat. Techn., vol. 1, pp. 40-47, 2016

work page 2016

[17] [17]

A high-order symplectic FDTD scheme for the Maxwell-Schrodinger system

X. Chen, F. M. Kong, K. Li, and M. Liu, “A high-order symplectic FDTD scheme for the Maxwell-Schrodinger system” IEEE J. Quantum Electron., vol. 54, no. 1, pp. 6100108, Feb. 2018

work page 2018

[18] [18]

A unified Hamiltonian solution to Maxwell–Schrödinger equation s for modeling electromagnetic fie ld-particle interaction,

Y. P. Chen, W. E.I. Sha, L. Jiang, M. Meng, Y. M. Wu, and W. C. Chew, “A unified Hamiltonian solution to Maxwell–Schrödinger equation s for modeling electromagnetic fie ld-particle interaction,” Comput. Phys. Commuun., vol. 215, pp. 63-70, Jun. 2017

work page 2017

[19] [19]

Krishnamurthy, K

V. Krishnamurthy, K. Ravi, and S.T. Ho, “Analytical framework f or the s t e a d y s t a t e a n a l y s i s o f w a v e l e ngth-dependent and intensity-dep endent interaction of multiple monochromatic beams in semiconductor ph otonic structures with multiple ac tive and passive sections,” IEEE J. Quantum Electron., vol. 48, no. 10, pp. 1282–1299, Oct. 2012

work page 2012

[20] [20]

Designing optimal universal pulses using second-order, large-scale, non-linear optimization,

C. K. Anand, A. D. Bain, A. T. Curtis, and Z. H. Nie, “Designing optimal universal pulses using second-order, large-scale, non-linear optimization,” J. Magn. Reson., vol, 219, pp. 61-74, Apr. 2012

work page 2012

[21] [21]

Application of the sympl ectic finite-difference time-domain scheme to electromagnetic simulat ion,

W. Sha, Z. Huang, X. Wu, and M. Chen, “Application of the sympl ectic finite-difference time-domain scheme to electromagnetic simulat ion,” J. Comput. Phys., vol. 225, no. 1, pp. 33-50. Jan. 2007

work page 2007

[22] [22]

Dispersion relation equation preserving FDTD method for nonlinear cubic Schrödinger equation,

T. W. H. Sheu and L. Lin, “Dispersion relation equation preserving FDTD method for nonlinear cubic Schrödinger equation,” J. Comput. Phys., vol. 299, no. 4, pp. 1-21. Oct. 2015

work page 2015

[23] [23]

Canonical symplectic structure and structure-preserving geomet ric algorithms for Schrödinger-Maxwell systems,

Q. Chen, H. Qin, J. Liu, J. Y. Xiao, R. L. Zhang, Y. He, and Y. L. Wang, “Canonical symplectic structure and structure-preserving geomet ric algorithms for Schrödinger-Maxwell systems,” J. Comput. Phys., vol. 349, pp. 441-452, 2017

work page 2017

[24] [24]

Symmetric and symplectic ERKN for oscillatory Hamiltonian systems,

Z. Chen, X. You, W. Shi, and Z. Liu, “Symmetric and symplectic ERKN for oscillatory Hamiltonian systems,” Comput. Phys. Commun., vol. 183, no.1, pp. 86-98, 2012

work page 2012

[25] [25]

Multiple image technique (MIT) a nd anisotropic perfectly matched layer (APML) in implementation of MRTD scheme for boundary truncations of microwave structures,

Q. Cao, Y. Chen, and R. Mittra, “Multiple image technique (MIT) a nd anisotropic perfectly matched layer (APML) in implementation of MRTD scheme for boundary truncations of microwave structures,” IEEE Trans. Microw. Theory Techn., vol. 50, no. 6, pp. 1578-1589, Jun. 2002

work page 2002

[26] [26]

Kira and S

M. Kira and S. W. Koch, Semiconductor Quantum Optics . Cambridge Univ. Press, 2012

work page 2012

[27] [27]

A new solution of Schrödinger equation based on symplectic algorithm,

Z. Huang, J. Xu, B. Sun, B. Wu, and X. Wu, “A new solution of Schrödinger equation based on symplectic algorithm,” Comput. Math. Applicat., vol. 69, no. 11, pp.1303-1312, Jun. 2015

work page 2015

[28] [28]

A symplectic FDTD algorit hm for the simulation of lossy dispersive materials,

H. wang, B. Wu, Z. Huang, and X. Wu, “A symplectic FDTD algorit hm for the simulation of lossy dispersive materials,” Comput. Phys. Commun., vol. 185, no.3, pp. 862-872, 2014

work page 2014

[29] [29]

Gerry and P

C. Gerry and P. Knight, Introduc tory Quantum Optics, Cambride University Press, New York, 2005

work page 2005

[30] [30]

Marlon. O. Scully, and M. Suhail. Zubairy, Quantum Optics. Cambridge, U.K.:Cambridge Univ. Press, 1997

work page 1997