Nonlinear dynamics of wave packets in parity-time-symmetric optical lattices near the phase-transition point is analytically studied. A nonlinear Klein–Gordon equation is derived for the envelope of these wave packets. A variety of phenomena known to exist in this envelope equation are shown to also exist in the full equation, including wave blowup, periodic bound states, and solitary wave solutions.

Two-dimensional localized edge modes in optical honeycomb lattices are found and analyzed analytically and computationally. Weak nonlinearity and transverse modulation are found to introduce self-phase modulation in the phase and create internal nonlinear interactions as the electromagnetic field propagates through the lattice. Even with relatively strong nonlinearity localization and persistence of modes along the edge are found.

Nonlinear waves in deformed optical honeycomb lattices are investigated. Discrete couple mode equations are used to find higher order continuous nonlinear Dirac systems which are employed to describe key underlying phenomena. For weak deformation and nonlinearity the wave propagation is circular–ellliptical. At strong nonlinearity the diffraction pattern becomes triangular in structure which is traced to appropriate nonequal energy propagation in momentum space. At suitably large deformation the dispersion structure can have nearby Dirac points or small gaps. The effective dynamics of the wave packets is described by two maximally balanced nonlocal nonlinear Schrödinger type equations.

We observe vortex generation by selective excitation of two honeycomb sublattices at the vicinity of Dirac points. Such vortices arise from graphene “pseudospin”, suggesting that “pseudospin” could be observable and possess real angular momentum.

We propose a conservation-dissipation formalism (CDF) for coarse-grained descriptions of irreversible processes. This formalism is based on a stability criterion for non-equilibrium thermodynamics. The criterion ensures that non-equilibrium states tend to equilibrium in long time. As a systematic methodology, CDF provides a feasible procedure in choosing non-equilibrium state variables and determining their evolution equations. The equations derived in CDF have a unified elegant form. They are globally hyperbolic, allow a convenient definition of weak solutions, and are amenable to existing numerics. More importantly, CDF is a genuinely nonlinear formalism and works for systems far away from equilibrium. With this formalism, we formulate novel thermodynamics theories for heat conduction in rigid bodies and non-isothermal compressible Maxwell fluid flows as two typical examples. In these examples, the non-equilibrium variables are exactly the conjugate variables of the heat fluxes or stress tensors. The new theory generalizes Cattaneo’s law or Maxwell’s law in a regularized and nonlinear fashion.