Consider electromagnetic waves in two-dimensional honeycomb structured media. The properties of transverse electric (TE) polarized waves are determined by the spectral properties of the elliptic operator L^A=-\nabla_x·A(x) \nabla_x, where A(x) is Λ_h− periodic (Λ_h denotes the equilateral triangular lattice), and such that with respect to some origin of coordinates, A(x) is PC− invariant (A(x)=\overline{A(-x)}) and 120° rotationally invariant (A(R*x)=R*A(x)R, where R is a 120° rotation in the plane). We first obtain results on the existence, stability and instability of Dirac points, conical intersections between two adjacent Floquet-Bloch dispersion surfaces. We then show that the introduction through small and slow variations of a domain wall across a line-defect gives rise to the bifurcation from Dirac points of highly robust (topologically protected) edge states. These are time-harmonic solutions of Maxwell's equations which are propagating parallel to the line-defect and spatially localized transverse to it. The transverse localization and strong robustness to perturbation of these edge states is rooted in the protected zero mode of a one-dimensional effective Dirac operator with spatially varying mass term. These results imply the existence of uni-directional propagating edge states for two classes of time-reversal invariant media in which C symmetry is broken: magneto-optic media and bi-anisotropic media. Our analysis applies and extends the tools previously developed in the context of honeycomb Schrödinger operators.

Photonic graphene, a photonic crystal with honeycomb structures, has been intensively studied in both theoretical and applied fields. Similar to graphene which admits Dirac Fermions and topological edge states, photonic graphene supports novel and subtle propagating modes (edge modes) of electromagnetic waves. These modes have wide applications in many optical systems. In this paper, we propose a novel gradient recovery method based on Bloch theory for the computation of topological edge modes in the honeycomb structure. Compared to standard finite element methods, this method provides higher order accuracy with the help of gradient recovery technique. This high order accuracy is highly desired for constructing the propagating electromagnetic modes in applications. We analyze the accuracy and prove the superconvergence of this method. Numerical examples are presented to show the efficiency by computing the edge mode for the P-symmetry and C-symmetry breaking cases in honeycomb structures.

The Onsager's reciprocal relation plays a fundamental role in the nonequilibrium thermodynamics. However, unfortunately, its classical version is valid only within a narrow region near equilibrium due to the linear regression hypothesis, which largely restricts its usage. In this paper, based on the conservation-dissipation formalism, a generalized version of Onsager's relations for the master equations and Fokker-Planck equations was derived. Nonlinear constitutive relations with nonsymmetric and positively stable operators, which become symmetric under the detailed balance condition, constitute key features of this new generalization. Similar conclusions also hold for many other classical models in physics and chemistry, which in turn make the current study as a benchmark for the application of generalized Onsager's relations in nonequilibrium thermodynamics.

The seek for a new universal formulation for describing various non-equilibrium processes is a central task of modern non-equilibrium thermodynamics. In this paper, a novel steady-state thermodynamic formalism was established for general Markov processes described by the Chapman-Kolmogorov equation. Furthermore, corresponding formalisms of steady-state thermodynamics for master equation and Fokker-Planck equation could be rigorously derived in mathematics. To be concrete, we proved that: 1) in the limit of continuous time, the steady-state thermodynamic formalism for the Chapman-Kolmogorov equation fully agrees with that for the master equation; 2) a similar one-to-one correspondence could be established rigorously between the master equation and Fokker-Planck equation in the limit of large system size; 3) when a Markov process is restrained to one-step jump, the steady-state thermodynamic formalism for the Fokker-Planck equation with discrete state variables also goes to that for master equations, as the discretization step gets smaller and smaller. Our analysis indicated that, with respect to the steady state, general Markov processes admit a unified and self-consistent non-equilibrium thermodynamic formulation, regardless of underlying detailed models.

We report on transport properties of the super-honeycomb lattice, the band structure of which possesses a flat band and Dirac cones, according to the tight-binding approximation. This super-honeycomb model combines the honeycomb lattice and the Lieb lattice and displays the properties of both. The super-honeycomb lattice also represents a hybrid fermionic and bosonic system, which is rarely seen in nature. By choosing the phases of input beams properly, the flat-band mode of the super-honeycomb will be excited and the input beams will exhibit strong localization during propagation. On the other hand, if the modes of Dirac cones of the super-honeycomb lattice are excited, one will observe conical diffraction. Furthermore, if the input beam is properly chosen to excite a sublattice of the super-honeycomb lattice and the modes of Dirac cones with different pseudospins, e.g., the three-beam interference pattern, the pseudospin-mediated vortices will be observed.