MathSciDoc: An Archive for Mathematician ∫

In this article, we propose a novel approach to construct macroscopic balance equations and constitutive equations describing various irreversible phenomena. It is based on the general principles of non-equilibrium thermodynamics and consists of four basic steps: picking suitable state variables, choosing a strictly concave entropy function, properly separating entropy fuxes and production rates, and determining a dissipation matrix. Our approach takes advantage of both extended irreversible thermodynamics and GENERIC formalisms and shows a direct correspondence with Levermore’s moment-closure hierarchies for the Boltzmann equation. As a direct application, a new ten-moment model beyond the classical hierarchies is constructed and is shown to recover the Euler equations in the equilibrium state. These interesting results may put various macroscopic modeling approaches, starting from the general principles of non-equilibrium thermodynamics, on a solid microscopic foundation based on the Boltzmann equation.

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.

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.

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.

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 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.

The linear spectrum and corresponding Bloch modes of shallow honeycomb lattices near Dirac points are investigated. Via perturbation theory, the dispersion relation is found to have threefold degeneracy at leading order with eigenvalue splitting at the following two orders; i.e., the threefold eigenvalue splits into single and double values. Multiscale perturbation methods are employed to describe the nonlinear dynamics of the associated wave envelopes. The dynamics of the envelope depends on different asymptotic balances whereupon a three-level nonlinear Dirac-type equation or a two-level nonlinear Dirac equation is derived. The analysis agrees well with direct numerical simulations.

Wave propagation in two-dimensional generalized honeycomb lattices is studied. By employing the tight-binding (TB) approximation, the linear dispersion relation and associated discrete envelope equations are derived for the lowest band. In the TB limit, the Bloch modes are localized at the minima of the potential wells and can analytically be constructed in terms of local orbitals. Bloch-mode relations are converted into integrals over orbitals. With this methodology, the linear dispersion relation is derived analytically in the TB limit. The nonlinear envelope dynamics are found to be governed by a unified nonlinear discrete wave system. The lowest Bloch band has two branches that touch at the Dirac points. In the neighborhood of these points, the unified system leads to a coupled nonlinear discrete Dirac system. In the continuous limit, the leading-order evolution is governed by a continuous nonlinear Dirac system. This system exhibits conical diffraction, a phenomenon observed in experiments. Coupled nonlinear Dirac systems are also obtained. Away from the Dirac points, the continuous limit of the discrete equation leads to coupled nonlinear Schrödinger equations when the underlying group velocities are nearly zero. With semiclassical approximations, all the parameters are estimated analytically.

Conical diffraction in honeycomb lattices is analyzed. This phenomenon arises in nonlinear Schrödinger equations with honeycomb lattice potentials. In the tight-binding approximation the wave envelope is governed by a nonlinear classical Dirac equation. Numerical simulations show that the Dirac equation and the lattice equation have the same conical diffraction properties. Similar conical diffraction occurs in both the linear and nonlinear regimes. The Dirac system reveals the underlying mechanism for the existence of conical diffraction in honeycomb lattices.

Pseudospin, an additional degree of freedom inherent in graphene, plays a key role in understanding many fundamental phenomena such as the anomalous quantum Hall effect, electron chirality and Klein paradox. Unlike the electron spin, the pseudospin was traditionally considered as an unmeasurable quantity, immune to Stern-Gerlach-type experiments. Recently, however, it has been suggested that graphene pseudospin is a real angular momentum that might manifest itself as an observable quantity, but so far direct tests of such a momentum remained unfruitful. Here, by selective excitation of two sublattices of an artificial photonic graphene, we demonstrate pseudospin-mediated vortex generation and topological charge flipping in otherwise uniform optical beams with Bloch momentum traversing through the Dirac points. Corroborated by numerical solutions of the linear massless Dirac-Weyl equation, we show that pseudospin can turn into orbital angular momentum completely, thus upholding the belief that pseudospin is not merely for theoretical elegance but rather physically measurable.