We introduce a web of strongly correlated interacting 3+1D topological superconductors/ insulators of 10 particular global symmetry groups of Cartan classes, realizable in electronic condensed matter systems, and their generalizations. The symmetries include SU(N), SU(2), U(1), fermion parity, time reversal and relate to each other through symmetry embeddings. We overview the lattice Hamiltonian formalism. We complete the list of bulk symmetry-protected topological invariants (SPT invariants/partition functions that exhibit boundary 't Hooft anomalies) via cobordism calculations, matching their full classification. We also present explicit 4-manifolds that detect these SPTs. On the other hand, once we dynamically gauge part of their global symmetries, we arrive in various new phases of $SU(N)$ Yang-Mills (YM) theories, realizable as quantum spin liquids with emergent gauge fields. We discuss how coupling YM theories to time reversal-SPTs affects the strongly coupled theories at low energy. For example, we point out a possibility of having two deconfined gapless time-reversal symmetric $SU(2)$ YM theories at $\theta=\pi$ as two distinct conformal field theories, which although are secretly indistinguishable by gapped SPT states nor by correlators of local operators on oriented spacetimes, can be distinguished on non-orientable spacetimes or potentially by correlators of extended operators.

Topological Quantum Field Theories (TQFTs) pertinent to some emergent low energy phenomena of condensed matter lattice models in 2+1 and 3+1 dimensions are explored. Many of our TQFTs are highly-interacting without free quadratic analogs. Some of our bosonic TQFTs can be regarded as the continuum field theory formulation of Dijkgraaf-Witten twisted discrete gauge theories. Other bosonic TQFTs beyond the Dijkgraaf-Witten description and all fermionic TQFTs (namely the spin TQFTs) are either higher-form gauge theories where particles must have strings attached, or fermionic discrete gauge theories obtained by gauging the fermionic Symmetry-Protected Topological states (SPTs). We analytically calculate both the Abelian and non-Abelian braiding statistics data of anyonic particle and string excitations in these theories, where the statistics data can one-to-one characterize the underlying topological orders of TQFTs. Namely, we derive path integral expectation values of links formed by line and surface operators in these TQFTs. The acquired link invariants include not only the familiar Aharonov-Bohm linking number, but also Milnor triple linking number in 3 dimensions, triple and quadruple linking numbers of surfaces, and intersection number of surfaces in 4 dimensions. We also construct new spin TQFTs with the corresponding knot/link invariants of Arf(-Brown-Kervaire), Sato-Levine and others. We propose a new relation between the fermionic SPT partition function and the Rokhlin invariant. As an example, we can use these invariants and other physical observables, including ground state degeneracy, reduced modular $\mathcal{S}^{xy}$ and $\mathcal{T}^{xy}$ matrices, and the partition function on $\mathbb{RP}^3$ manifold, to identify all $\nu \in \mathbbm{Z}_8$ classes of 2+1 dimensional gauged $\mathbb{Z}_2$-Ising-symmetric $\mathbb{Z}_2^f$-fermionic Topological Superconductors (realized by stacking $\nu$ layers of a pair of chiral and anti-chiral $p$-wave superconductors [$p_x+ip_y$ and $p_x-ip_y$], where boundary supports non-chiral Majorana-Weyl modes) with continuum spin-TQFTs.

Solitary waves bifurcated from edges of Bloch bands in two-dimensional periodic media are determined both analytically and numerically in the context of a two-dimensional nonlinear Schrödinger equation with a periodic potential. Using multiscale perturbation methods, the envelope equations of solitary waves near Bloch bands are analytically derived. These envelope equations reveal that solitary waves can bifurcate from edges of Bloch bands under either focusing or defocusing nonlinearity, depending on the signs of the second-order dispersion coefficients at the edge points. Interestingly, at edge points with two linearly independent Bloch modes, the envelope equations lead to a host of solitary wave structures, including reduced-symmetry solitons, dipole-array solitons, vortex-cell solitons, and so on—many of which have not been reported before to our knowledge. It is also shown analytically that the centers of envelope solutions can be positioned at only four possible locations at or between potential peaks. Numerically, families of these solitary waves are directly computed both near and far away from the band edges. Near the band edges, the numerical solutions spread over many lattice sites, and they fully agree with the analytical solutions obtained from the envelope equations. Far away from the band edges, solitary waves are strongly localized, with intensity and phase profiles characteristic of individual families.

Previous work has shown that in a two-dimensional periodic medium under focusing or defocusing cubic nonlinearities, gap solitons in the form of low-amplitude and slowly modulated single-Bloch-wave packets can bifurcate out from the edges of Bloch bands. In this paper, linear stability properties of these gap solitons near band edges are determined both analytically and numerically. Through asymptotic analysis, it is shown that these gap solitons are linearly unstable if the slope of their power curve at the band edge has the opposite sign of nonlinearity
here focusing nonlinearity is said to have a positive sign, and defocusing nonlinearity to have a negative sign. An equivalent condition for linear instability is that the power of the gap solitons near the band edge is lower than the limit power value on the band edge. Through numerical computations of the power curves, it is found that this condition is always satisfied, thus two-dimensional gap solitons near band edges are linearly unstable. The analytical formula for the unstable eigenvalue of gap solitons near band edges is also asymptotically derived. It is shown that this unstable eigenvalue is proportional to the cubic power of the soliton’s amplitude, and it induces width instabilities of gap solitons. A comparison between this analytical eigenvalue formula and numerically computed eigenvalues shows excellent agreement.

In this paper, we present a novel pairwise-force smoothed particle hydrodynamics (PF-SPH) model to allow modeling of various interactions at interfaces in real time. Realistic capture of interactions at interfaces is a challenging problem for SPH-based simulations, especially for scenarios involving multiple interactions at different interfaces. Our PF-SPH model can readily handle multiple kinds of interactions simultaneously in a single simulation; its basis is to use a larger support radius than that used in standard SPH. We adopt a novel anisotropic filtering term to further improve the performance of interaction forces. The proposed model is stable; furthermore, it avoids the particle clustering problem which commonly occurs at the free surface. We show how our model can be used to capture various interactions. We also consider the close connection between droplets and bubbles, and show how to animate bubbles rising in liquid as well as bubbles in air. Our method is versatile, physically plausible and easy-to-implement. Examples are provided to demonstrate the capabilities and effectiveness of our approach.