Publications

Order by: publication date number of citations

Python-JAX-based fast Stokesian dynamics

Kim William Torre, Raoul D. Schram, Joost de Graaf

SciPost Phys. Codebases 56 (2025) · published 30 June 2025 | · pdf

Stokesian Dynamics (SD) is a powerful computational framework for simulating the motion of particles in a viscous Newtonian fluid under Stokes-flow conditions. Traditional SD implementations can be computationally expensive as they rely on the inversion of large mobility matrices to determine hydrodynamic interactions. Recently, however, the simulation of thermalized systems with large numbers of particles has become feasible [Fiore and Swan, J. Fluid Mech. 878, 544 (2019)]. Their "fast Stokesian dynamics" (FSD) method leverages a saddle-point formulation to ensure overall scaling of the algorithm that is linear in the number of particles $\mathcal{O}(N)$; performance relies on dedicated graphics-processing-unit computing. Here, we present a different route toward implementing FSD, which instead leverages the Just-in-Time (JIT) compilation capabilities of Google JAX. We refer to this implementation as JFSD and perform benchmarks on it to verify that it has the right scaling and is sufficiently fast by the standards of modern computational physics. In addition, we provide a series of physical test cases that help ensure accuracy and robustness, as the code undergoes further development. Thus, JFSD is ready to facilitate the study of hydrodynamic effects in particle suspensions across the domains of soft, active, and granular matter.

2 citations
Codebase release 0.2 for JFSD

Kim William Torre, Raoul D. Schram, Joost de Graaf

SciPost Phys. Codebases 56-r0.2 (2025) · published 30 June 2025 | · src

Stokesian Dynamics (SD) is a powerful computational framework for simulating the motion of particles in a viscous Newtonian fluid under Stokes-flow conditions. Traditional SD implementations can be computationally expensive as they rely on the inversion of large mobility matrices to determine hydrodynamic interactions. Recently, however, the simulation of thermalized systems with large numbers of particles has become feasible [Fiore and Swan, J. Fluid Mech. 878, 544 (2019)]. Their "fast Stokesian dynamics" (FSD) method leverages a saddle-point formulation to ensure overall scaling of the algorithm that is linear in the number of particles $\mathcal{O}(N)$; performance relies on dedicated graphics-processing-unit computing. Here, we present a different route toward implementing FSD, which instead leverages the Just-in-Time (JIT) compilation capabilities of Google JAX. We refer to this implementation as JFSD and perform benchmarks on it to verify that it has the right scaling and is sufficiently fast by the standards of modern computational physics. In addition, we provide a series of physical test cases that help ensure accuracy and robustness, as the code undergoes further development. Thus, JFSD is ready to facilitate the study of hydrodynamic effects in particle suspensions across the domains of soft, active, and granular matter.

1 citation
Knizhnik-Zamolodchikov equations and integrable hyperbolic Landau-Zener models

Suvendu Barik, Lieuwe Bakker, Vladimir Gritsev, Emil A. Yuzbashyan

SciPost Phys. 18, 212 (2025) · published 27 June 2025 | · pdf

We study the relationship between integrable Landau-Zener (LZ) models and Knizhnik-Zamolodchikov (KZ) equations. The latter are originally equations for the correlation functions of two-dimensional conformal field theories, but can also be interpreted as multi-time Schrödinger equations. The general LZ problem is to find probabilities of tunneling from eigenstates at $t=t_\text{in}$ to eigenstates at $t\to+∞$ for an $N× N$ time-dependent Hamiltonian $\hat H(t)$. A number of such problems are exactly solvable in the sense that their tunneling probabilities are elementary functions of Hamiltonian parameters. Recently, it has been proposed that exactly solvable LZ models of this type map to KZ equations. Here we use this connection to identify and solve a class of integrable LZ models with hyperbolic time dependence, $\hat H(t)=\hat A+\hat B/t$, for $N=2, 3$, and $4$, where $\hat A$ and $\hat B$ are time-independent matrices.

1 citation
On symmetry-resolved generalized entropies

Fei Yan, Sara Murciano, Pasquale Calabrese, Robert Konik

SciPost Phys. 18, 211 (2025) · published 27 June 2025 | · pdf

Symmetry-resolved entanglement, capturing the refined structure of quantum entanglement in systems with global symmetries, has attracted a lot of attention recently. In this manuscript, introducing the notion of symmetry-resolved generalized entropies, we aim to develop a computational framework suitable for the study of excited state symmetry-resolved entanglement as well as the dynamical evolution of symmetry-resolved entanglement in symmetry-preserving out-of-equilibrium settings. We illustrate our framework using the example of (1+1)-d free massless compact boson theory, and benchmark our results using lattice computation in the XX chain. As a byproduct, our computational framework also provides access to the probability distribution of the symmetry charge contained within a subsystem and the corresponding full counting statistics.

1 citation
Boundary operator expansion and extraordinary phase transition in the tricritical O(N) model

Xinyu Sun, Shao-Kai Jian

SciPost Phys. 18, 210 (2025) · published 27 June 2025 | · pdf

We study the boundary extraordinary transition of a three-dimensional (3D) tricritical $O(N)$ model. We first compute the mean-field Green's function with a general coupling of $|\vec \phi|^{2n}$ (with $n=3$ corresponding to the tricritical model) at the extraordinary phase transition. Then, using layer susceptibility, we obtain the boundary operator expansion for the transverse and longitudinal modes within the $\epsilon=3 - d$ expansion. Based on these results, we demonstrate that the tricritical point exhibits an extraordinary transition characterized by an ordered boundary for any $N$. This provides the first nontrivial example of continuous symmetry breaking in 2D in the context of boundary criticality.

3 citations
Hybrid symmetry class topological insulators

Sanjib Kumar Das, Bitan Roy

SciPost Phys. 18, 209 (2025) · published 27 June 2025 | · pdf

Traditional topological materials belong to different Altland-Zirnbauer symmetry classes (AZSCs) depending on their non-spatial symmetries. Here we introduce the notion of hybrid symmetry class topological insulators (HSCTIs): A fusion of two different AZSC topological insulators (TIs) such that they occupy orthogonal Cartesian hyperplanes and their universal massive Dirac Hamiltonian mutually anticommute, a mathematical procedure we name hybridization. The boundaries of HSCTIs can also harbor TIs, typically affiliated with an AZSC that is different from the ones for the parent two TIs. As such, a fusion or hybridization between planar class AII quantum spin Hall and vertical class BDI Su-Schrieffer-Heeger insulators gives birth to a three-dimensional class A HSCTI, accommodating quantum anomalous Hall insulators (class A) of opposite Chern numbers and quantized Hall conductivity of opposite signs on the top and bottom surfaces. Such a response is shown to be stable against weak disorder. We extend this construction to encompass crystalline HSCTI and topological superconductors (featuring half-quantized thermal Hall conductivity of opposite sings on the top and bottom surfaces), and beyond three spatial dimensions. Non-trivial responses of three-dimensional HSCTIs to crystal defects (namely edge dislocations) in terms of mid-gap bound states at zero energy around its core only on the top and bottom surfaces are presented. Possible (meta)material platforms to harness and engineer HSCTIs are discussed.

1 citation
Quasiperiodic quadrupole insulators

Raul Liquito, Miguel Gonçalves, Eduardo V. Castro

SciPost Phys. 18, 208 (2025) · published 26 June 2025 | · pdf

Higher-order topological insulators are an intriguing new family of topological states that host lower-dimensional boundary states. Concurrently, quasiperiodic systems have garnered significant interest due to their complex localization and topological properties. In this work we study the impact of chiral symmetry preserving quasiperiodic modulations on the paradigmatic Benalcazar-Bernevig-Hughes model, which hosts topological insulating phases with zero-energy sublattice-polarized modes. We find that the topological properties are not only robust to the quasiperiodic modulation, but can even be enriched. In particular, we unveil the first instance of a quasiperiodic induced second-order topological insulating phase. Furthermore, in contrast with disorder, we find that quasiperiodic modulations can induce multiple reentrant topological transitions, showing an intricate sequence of localization properties. Our results open a promising avenue for exploring the rich interplay between higher-order topology and quasiperiodicity.
Theorems for the lightcone bootstrap

Balt C. van Rees

SciPost Phys. 18, 207 (2025) · published 26 June 2025 | · pdf

Consider a conformally covariant four-point function of identical scalar operators with a discrete spectrum, a twist gap, and compatible with the unitarity conditions. We give a mathematical proof confirming that the spectrum and OPE coefficients at large spin and fixed twist always become that of a generalized free field theory.

2 citations
The fermionic double smeared null energy condition

Duarte Fragoso, Lihan Guo

SciPost Phys. Core 8, 045 (2025) · published 26 June 2025 | · pdf

Energy conditions are crucial for understanding why exotic phenomena such as traversable wormholes and closed timelike curves remain elusive. In this paper, we prove the Double Smeared Null Energy Condition (DSNEC) for the fermionic free theory in 4-dimensional flat Minkowski spacetime, extending previous work on the same energy condition for the bosonic case [J. R. Fliss and B. Freivogel, SciPost Phys. 12, 084 (2022); J. R. Fliss et al., SciPost Phys. 14, 024 (2023)] by adapting Fewster and Mistry's method [C. J. Fewster and B. Mistry, Phys. Rev. D 68, 105010 (2003)] to the energy-momentum tensor $T_{++}$. A notable difference from previous works lies in the presence of the $\gamma_0 \gamma_+$ matrix in $T_{++}$, causing a loss of symmetry. This challenge is addressed by making use of its square-root matrix. We provide explicit analytic results for the massless case as well as numerical insights for the mass-dependence of the bound in the case of Gaussian smearing.
Interaction-induced strong zero modes in short quantum dot chains with time-reversal symmetry

A. Mert Bozkurt, Sebastian Miles, Sebastiaan L. D. ten Haaf, Chun-Xiao Liu, Fabian Hassler, Michael Wimmer

SciPost Phys. 18, 206 (2025) · published 25 June 2025 | · pdf

We theoretically explore the emergence of strong zero modes in a two-site chain consisting of two quantum dots coupled due to a central dot that mediates electron hopping and singlet superconducting pairing. In the presence of time-reversal symmetry, the on-site Coulomb interaction leads to a three-fold ground-state degeneracy when tuning the system to a sweet spot as a function of the inter-dot couplings. This degeneracy is protected against changes of the dot energies in the same way as "poor man's" Majorana bound states in short Kitaev chains. In the limit of strong interactions, this protection is maximal and the entire spectrum becomes triply degenerate, indicating the emergence of a "poor man's" version of a strong zero mode. We explain the degeneracy and protection by constructing corresponding Majorana Kramers-pair operators and $\mathbb{Z}_3$-parafermion operators. The strong zero modes share many properties of Majorana bound states in short Kitaev chains, including the stability of zero-bias peaks in the conductance and the behavior upon coupling to an additional quantum dot. However, they can be distinguished through finite-bias spectroscopy and the exhibit a different behavior when scaling to longer chains.

1 citation

Previous 1 ... 24 25 26 27 28 29 30 31 32 33 34 ... 363 Next

Publications

2 citations

1 citation

1 citation

1 citation

3 citations

1 citation

2 citations

1 citation