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.

We use the polynomial formulation of the holomorphic anomaly equations governing perturbative topological string theory to derive the free energies in a scaling limit to all orders in perturbation theory for any CalabiYau threefold. The partition function in this limit satisfies an Airy differential equation in a rescaled topological string coupling. One of the two solutions of this equation gives the perturbative expansion and the other solution provides geometric hints of the non-perturbative structure of topological string theory. Both solutions can be expanded naturally around strong coupling.

In recent years, fermionic topological phases of quantum matter has attracted a lot of attention. In a pioneer work by Gu, Wang and Wen, the concept of equivalence classes of fermionic local unitary(FLU) transformations was proposed to systematically understand non-chiral topological phases in 2D fermion systems and an incomplete classification was obtained. On the other hand, the physical picture of fermion condensation and its corresponding super pivotal categories give rise to a generic mathematical framework to describe fermionic topological phases of quantum matter. In particular, it has been pointed out that in certain fermionic topological phases, there exists the so-called q-type anyon excitations, which have no analogues in bosonic theories. In this paper, we generalize the Gu, Wang and Wen construction to include those fermionic topological phases with q-type anyon excitations. We argue that all non-chiral fermionic topological phases in 2+1D are characterized by a set of tensors (N^{ij}_k,F^{ij}_k,F^{ijm,αβ}_{kln,χδ},n_i,d_i), which satisfy a set of nonlinear algebraic equations parameterized by phase factors Ξ^{ijm,αβ}_{kl}, Ξ^{ij}_{kln,χδ}, Ω^{kim,αβ}_{jl} and Ω^{ki}_{jln,χδ}. Moreover, consistency conditions among algebraic equations give rise to additional constraints on these phase factors which allow us to construct a topological invariant partition for an arbitrary triangulation of 3D spin manifold. Finally, several examples with q-type anyon excitations are discussed, including the Fermionic topological phase from Tambara-Yamagami category for \mathbb{Z}_{2N}, which can be regarded as the \mathbb{Z}_{2N} parafermion generalization of Ising fermionic topological phase.

We derive a lower bound on the ground state energy of the Hubbard model for given value of the total spin. In combination with the upper bound derived previously by Giuliani, our result proves that in the low density limit, the leading order correction compared to the ground state energy of a non-interacting lattice Fermi gas is given by $8\pi a \rho_u \rho_d$, where $\rho_{u(d)}$ denotes the density of the spin-up (down) particles, and $a$ is the scattering length of the contact interaction potential. This result extends previous work on the corresponding continuum model to the lattice case.

In recent years, great success has been achieved on the classification of symmetry-protected topological (SPT) phases for interacting fermion systems by using generalized cohomology theory. However, the explicit calculation of generalized cohomology theory is extremely hard due to the difficulty of computing obstruction functions. In this paper, based on the physical picture of topological invariants and mathematical techniques in homotopy algebra, we develop an algorithm to resolve this hard problem. It is well known that cochains in the cohomology of the symmetry group, which are used to enumerate the SPT phases, can be expressed equivalently in different linear bases, known as the resolutions. By expressing the cochains in a reduced resolution containing much fewer basis than the choice commonly used in previous studies, the computational cost is drastically reduced. In particular, it reduces the computational cost for infinite discrete symmetry groups, like the wallpaper groups and space groups, from infinity to finity. As examples, we compute the classification of two-dimensional interacting fermionic SPT phases, for all 17 wallpaper symmetry groups.