Papers

We show that chiral matters in the minimal supersymmetric standard model (MSSM) with semi-realistic flavor structures can be obtained from 7-brane configurations of magnetized super Yang-Mills (SYM) theory on a toroidal orbifold $R^{1,3} \times (T^2)^3/(Z_2 \times Z'_2)$, where background magnetic fluxes and Wilson-lines are turned on preserving four-dimensional ${\cal N}=1$ supersymmetry. The zero-mode spectrum of chiral multiplets in total is just MSSM ones, except the existence of those for three generations of right-handed neutrino and extra generations of MSSM Higgs pairs. Hierarchical Yukawa couplings can be obtained from the overlap integrals of wavefunctions localized in extra dimensions, allowing semi-realistic patterns of flavor structures for quarks and charged leptons. We also develop a systematic way to embed additional 7-branes into the configuration, those are sequestered from the visible sector toward a hidden sector model building.

We investigate Type IIB compactifications with spatially varying fluxes and dilaton profiles in the setting of dynamical cobordism. In particular, we analyze a G-theory motivated compactification in which the fluxes and the dilaton depend on coordinates of a complex two-dimensional plane. From the equations of motion, we deduce the existence of End of the World branes. In a cohomological interpretation, these branes appear precisely in order to trivialize the relevant cohomology class. Furthermore, we compute the associated bordism group and show that additional non-perturbative objects are needed to cancel the class, while retaining the cohomological contribution as a subgroup. This suggests a mathematical structure that connects energy scales with the emergence of perturbative and non-perturbative physics.

We compute the joint distributions of arbitrary numbers of eigenvectors of real and complex symmetric random tensors by the quantum field theoretical methods which were previously used to compute the mean distributions. We obtain the random matrix representations and the large-dimension asymptotics of the joint distributions. The latter can be expressed by a common function of tensor geometries, extending the universality found for the mean distributions to the joint distributions. Several crosschecks of our results are carried out by Monte Carlo computations.

We introduce the gravielectric (GE) and gravimagnetic (GM) fields in stationary spacetime using the Komar two-form and its dual. This opens the way to extend the Komar-Tomimatsu derivation of mass formulas to a more detailed picture in terms of the local lines of force. We show that Misner strings (MS) carry singular GE and GM fluxes connecting the horizon and the asymptotic zone. Moreover, MS are laterally transparent, so field lines can flow in and out of the bulk. This explains why the usual Komar mass integrals around the Misner strings in the Taub-NUT vacuum solution are negative: the pattern of field lines shows that they flow onto the string from the horizon, so it is necessary to calculate the incoming (positive) but not the outgoing Komar fluxes. This incoming flux is then turned back to the horizon through the Misner strings, realizing the closed circuit without sources. So Misner strings are massless empty tubes, but not rigid rods of negative mass. Similarly, GM field lines can connect positively and negatively charged regions of the horizon, generating, for example, the gravimagnetic dipole moment of the Kerr metric.

We investigate the divergent perturbative series of correlation functions for a massless, self-interacting scalar field in de Sitter space. First, we use our previously proposed method of autonomous equations to obtain finite time-dependent functions, and show that these functions approximate the time evolution of the correlation functions of the stochastic theory reasonably well. Second, we apply the technique of autonomous equations to the Borel-Le Roy transforms of correlation functions, and use solutions of these equations to perform Borel resummation. The results match the time evolution obtained in the stochastic picture substantially better. In addition, we propose an alternative method for extracting perturbative coefficients and provide a new derivation of our autonomous equation by truncating a system of Schwinger-Dyson-type differential equations.

Effective field theories provide a suitable framework for both particle physics and statistical physics. We delve deeper into the study of the effective three-dimensional scalar field theory for its application to statistical physics, especially considering the role of the sextic coupling in the tricritical-to-critical crossover. The three-loop renormalization of the mass and the two coupling constants that we perform allows us to obtain, for the first time, the complete renormalization group flow of the couplings in that order. We analyze what universality means in this problem and how we can recover non-universal terms from the renormalization group beta functions. The crossover is realized by the convergence of the renormalization group flow towards the line connecting the tricritical and critical fixed points.

Within modified gravity, we study static spherically and axially symmetric self-gravitating non-Abelian monopoles in $SU(2)$ Yang-Mills-Higgs theory. By comparing these monopoles with those obtained in Einstein-Yang-Mills-Higgs theory, we identify the differences introduced by the modification of gravity and show that they can be quite significant for systems with strong Higgs self-coupling.

We study the postinflationary dynamics of an Einstein-Cartan-Holst gravity-motivated inflationary scenario, known as Einstein-Cartan pseudoscalaron inflation, coupled to a type-I seesaw extension of the Standard Model with three heavy right-handed Majorana neutrinos. In particular, we show that nonthermal leptogenesis emerges as a necessary and self-consistent mechanism for generating the observed baryon asymmetry of the Universe, mainly because of the universal coupling of the inflaton to the additional heavy Majorana fermions. The resulting framework provides theoretical predictions that are fully compatible with the latest cosmological constraints from the Cosmic Microwave Background, Baryon Acoustic Oscillations, and Big Bang Nucleosynthesis, as well as with neutrino oscillation experiments, for a wide range of the fundamental Barbero-Immirzi model parameter $\gamma$, which controls the inflationary and postinflationary phases. In particular, for $\gamma \sim -1/100$ and a lightest Majorana-neutrino mass of order $10^{13}$ GeV, we find a scalar spectral index $n_s \sim 0.970$, a tensor-to-scalar ratio $r \sim 0.004$, for a number of e-folds before the end of inflation $N_e \lesssim 60$, and a baryon-to-entropy ratio $n_B/s \sim 8.7 \times 10^{-11}$.

We develop a generating-function formulation for the symbolic reduction of multi-loop Feynman integrals. In this framework, integration-by-parts identities are rewritten as differential equations for sector-wise generating functions, so the reduction problem can be studied in a non-commutative algebra of differential operators rather than only through relations among individual integrals. This viewpoint leads to an iterative algorithm that generates candidate equations, extracts symbolic reduction rules, updates the active rule set, and tests completeness on the lattice of integral indices. We illustrate the method with the sunset topology, planar and non-planar massless double-box topologies, representative subsectors, and a degenerate example in which the top sector contains no master integral. Together, these examples show how symbolic reduction rules, descendant equations, and completeness criteria can be organized within a single algebraic framework.

A constructive algorithm based on the theory of spectral pairs for constructing nonuniform wavelet basis in $L^2(\mathbb R)$ was considered by Gabardo and Nashed (J Funct. Anal. 158:209-241, 1998). In this setting, the associated translation set $\Lambda =\left\{ 0,r/N\right\}+2\,\mathbb Z$ is no longer a discrete subgroup of $\mathbb R$ but a spectrum associated with a certain one-dimensional spectral pair and the associated dilation is an even positive integer related to the given spectral pair. The main objective of this paper is to develop oblique and unitary extension principles for the construction nonuniform wavelet frames over non-Archimedean Local fields of positive characteristic. An example and some potential applications are also presented.

We present a method to quantize free fermions which eliminates the doublers when implemented on the lattice in any number of dimensions and in the $m=0$ limit. The elimination of doublers is achieved by combining a second-order description of fermions, with the tools associated with non-Hermitian Hamiltonians. We identify a new Pseudo-Hermitian symmetry of the second-order fermion equation, and we identify the associated $U(1)$ symmetry which will become charge when shifted to a local gauge theory. We validated the methods numerically.

We study finite-temperature $P$-form $\mathbb Z_N$ homological codes via an exact finite-Trotter quantum-to-classical map to a $(d+1)$-dimensional spacetime model with electric and magnetic topological background charges. The resulting background-resolved partition functions admit an exact reformulation in terms of closed magnetic and electric defect polymers, with opposite-species interactions governed by linking phases. By bounding this complex polymer gas by positive same-species hard-core majorant gases, we obtain an explicit low-activity criterion under which all background-dependent partition functions are uniformly controlled and homologically nontrivial polymers are exponentially suppressed on the scale of the spacetime systole. We also derive an exact higher-form Kramers-Wannier duality exchanging electric and magnetic backgrounds, Wilson and 't Hooft operators, and $P$-form and $(d-P)$-form theories. Finally, for prime $N$, we identify an exact source-free gauge-theory specialization coupled to the plaquette random-cluster model, which imports sharp phase-transition results on special geometries into the spacetime framework.

Machine learning has become a powerful tool in high-energy collider experiments, which enables the studies based on data-driven approaches to complex reconstruction and regression tasks. The study of identified hadron spectra in pseudorapidity regions beyond detector acceptance, which is limited to mid-rapidity regions, carries important information about particle production, yet remains unmeasured. In this work, we develop a physics-informed neural network, trained on PYTHIA8 $pp$ collisions at $\sqrt{s}=13.6$ TeV, to infer $p_{\rm T}$ spectra of $\pi^{\pm}$, $K^{\pm}$, $p/\bar{p}$, $\Lambda/\bar{\Lambda}$, and $K^{0}_{\mathrm{s}}$ in different rapidity regions. Physics-motivated constraints, including particle yield ratios, spectral shape, and smoothness, are incorporated into the loss function. A staged hyperparameter optimization strategy is used to ensure stability. The model achieves yield uncertainties of ${\sim}1.5\%$, $1.8\%$, and $5.83\%$ in the training, interpolation, and extrapolation regimes, respectively, outperforming XGBoost and LightGBM. It further reproduces key observables such as particle yield ratios, the multiplicity dependence of $\langle p_{\rm T} \rangle$, and kinetic freeze-out parameters, indicating that the model captures the underlying physics and provides reliable predictions beyond the measured phase space.

We analyze F-term hybrid inflation (FHI) within various grand unified theories (GUTs) in the presence of a Kähler potential for the inflaton field which parameterizes the Kähler manifolds $SU(1,1)/U(1)$ or $SU(2)/U(1)$. We take into account supergravity, radiative, and soft supersymmetry-breaking corrections to the tree-level potential and find that viable FHI can be realized without extrema along the inflationary trajectory for a broad region of the parameter space. For selected superpotential parameters the models' predictions are largely influenced by the curvature of the internal space and the magnitude of the tadpole parameter which are constrained so as to achieve compatibility with the current ACT and SPT data. We also discuss the formation of cosmic strings and their associated gravitational wave signals, potentially detectable by current and upcoming experiments.

We compare Picard–Lefschetz theory and resurgence in three basic one-dimensional exponential integrals: the Airy model, the Bessel model, and the Gamma model. On the Picard–Lefschetz side, we describe the Lefschetz thimbles and compute the connecting trajectories between critical points appearing at Stokes phases. On the resurgent side, we analyze the Borel singularities of the saddle expansions and use alien operators to recover the same Stokes coefficients. These examples serve as explicit finite-dimensional test cases for the dictionary between thimble wall-crossing and alien calculus.

We extend product Chern-Simons theory to develop several mixed $U(1)\times U(1)$ models where one gauge field is governed by a Chern-Simons term and the other by a Maxwell or Born-Infeld term. We show that, by choosing suitable potentials, the energy functional admits a topological lower bound saturated by first-order self-dual equations. The resulting dyonic systems can be divided into vortex-vortex and vortex-antivortex configurations, and the coexistence of vortices and antivortices in the latter extends the vortex-only result known in product Chern-Simons model. On a doubly periodic domain, we establish Bradlow-type bounds with distinct physical implications: for vortex-only systems, the vortex numbers stay below this bound and cannot be arbitrarily large; for vortex-antivortex systems, the bound is imposed on the difference between the vortex and antivortex numbers, while the individual numbers are arbitrary. This distinction results in a bounded energy spectrum for the former and an unbounded energy spectrum for the latter.

We develop a descriptive account of scientific reward in physics based on the concept of the time-dependent Polydoxon, defined as the structured set of empirically viable theories at a given time. We argue that highly rewarded contributions, such as those recognized by major prizes and professional honors, can be systematically understood as those that transform this space. These transformations take the form of expansion (adding viable theories), contraction (eliminating viable theories), reconfiguration (illuminating deeper structures and relations within and between theories), and enabling moves (methodological or technological advances that enable future transformations). The analysis is further refined by emphasizing that reward correlates with the transformation's magnitude, assessed along dimensions of scope, centrality, depth, and future leverage. This framework reframes the analysis of rewarded achievement away from isolated theoretical successes and toward the dynamics of a landscape of viable theories, providing a more unified descriptive interpretation of rewarded scientific activity in physics across its diverse set of theoretical and experimental discoveries.

In this work we derive novel analytical expressions for the mass and distance of a Schwarzschild black hole (BH), as well as for the orbital radius of test particles orbiting it, it terms of astrophysical observables measured throughout the entire orbit of the revolving particle. We use a general relativistic method to describe the frequency shifts of photons emitted in the vivinity of a BH by considering two emitters (or two positions of the same emitter) located symmetrically opposite to each other with respect to the observer's line of sight (LOS) when performing measurements along the orbit. Furthermore, the introduction of the redshift rapidity allows us to write independent expressions for the BH mass and its distance to Earth. We also extend our study to the case when astrophysical systems have a peculiar motion and derive the corresponding closed formulas.

We compute the logarithmic one-loop matter contribution to the thermodynamics and homogeneous twistless torsional Newton–Cartan scale response of a four-dimensional Lifshitz black brane. The background is the neutral planar member of the analytic Einstein–Maxwell–dilaton Lifshitz black-brane family, while the quantum fields are treated as probes: a real scalar with arbitrary mass and nonminimal curvature coupling, a four-component Dirac spinor, and an Abelian Maxwell field with its Faddeev–Popov ghosts. For each spin sector, the logarithmic coefficient separates into a smooth radial heat-kernel contribution, governed by $\calC_1$, and a horizon-localized conical contribution, governed by $\calW_h$. This separation identifies which part of the matter one-loop logarithm is visible in the boundary Lifshitz/TTNC Ward identity and which part instead belongs to horizon replica entropy. The smooth coefficient $\calC_1$ controls the source-induced homogeneous projection of the TTNC Weyl Ward identity, namely the smooth boundary-source contribution to the Ward combination $z\epsilon-2p$, whereas the thermodynamic logarithmic entropy is controlled by $\calC_{\rm therm}=\calC_1+zL^2\calW_h$. We give closed expressions for $\calC_1$, $\calW_h$, and $\calC_{\rm therm}$ for the scalar, Dirac, and Maxwell probe sectors, identify the gauge-field contact contribution in the conical entropy, and verify that the smooth source response vanishes in the relativistic planar $z=1$ limit while the standard horizon heat-kernel entropy coefficient remains.

Negative masses have long been explored, but their observational viability remains unclear. In this work, we develop a unified, observationally testable framework to constrain negative masses using both coupling level and dynamical probes. We establish that while dipole radiation bounds require universality of gravitational charge, the intrinsic dynamics of negative mass binaries generically lead to anomalous behaviors such as anti-chirps, dispersal and runaway motion. These signatures are absent in current gravitational wave observations, providing a robust exclusion channel independent of modified gravity assumptions.

We study the implementation of the island prescription in fuzzball-inspired models of black holes. As a simplified setup, we model a fuzzball by replacing the event horizon with a reflecting boundary (stretched horizon). In the framework of two-dimensional model with such boundary, we analyze the dynamics of entanglement entropy. We find that the presence of the boundary modifies the behavior of the island saddle, and for a range of parameter values we observe the effect of blinking island found in arXiv:2311.16244 which inevitably leads to the analogue of information paradox. We then extend the analysis to higher dimensions, incorporating both bulk and boundary contributions to the generalized entropy. The existence of island solutions is found to depend sensitively on the boundary conditions and the position of the stretched horizon, naturally leading to the absence of entanglement islands for a wide range of parameters. Finally, we consider more "realistic" stringy fuzzball geometries, including superstrata and bubbling solutions, and estimate whether island solutions can arise in these backgrounds. The results indicate that the existence of islands depends on the behavior of the geometric area near the cap, and is not guaranteed in general.

We investigate DC transport in a top-down construction of thermal QGP-like theories using a holographic $\mathcal{M}$-theoretic background, incorporating quartic curvature corrections. The DC and Hall conductivities are computed from the Dirac-Born-Infeld (DBI) action of the corresponding type-IIA probe $D6$ flavor branes via the reality condition method proposed in arXiv:0708.1994, arXiv:0705.3870. We further analyze pair-production contributions in the presence and absence of an external magnetic field and work out regimes where thermal pair production dominates over charge carrier transport. These findings extend earlier AdS/CFT results to non-conformal, higher-derivative settings relevant to thermal QGP-like theories.

The Bekenstein-Hawking entropy formula $\rho = e^{A/4G}$ receives significant corrections for charged black holes near extremality. Using standard results in JT gravity, the correction term can semiclassically be expressed as minus the exponential of the inner horizon area, $e^{A_{\text{inner}}/4G}$, and the cancellation between these two exponentials enforces a vanishing density of states towards extremality, when the two horizons collide. Building on arXiv:2402.10162, we argue that the correction term corresponds to a complex saddle geometry of the bulk gravitational path integral. The proposed geometry has a negative boundary length and caps off at the inner horizon; we refer to it as the inner horizon saddle. We discuss how the saddle, and its accompanying minus sign, contribute to the density of states through a Picard-Lefschetz analysis of the inverse Laplace contour, together with a stability analysis of the saddle. We also address the inner horizon saddle's violation of the Kontsevich-Segal-Witten (KSW) allowability criterion for the inclusion of complex metrics. Despite this violation, which is believed to cause unphysical divergences in path integral computations, one can describe one-loop effects on the inner horizon saddle by carefully treating wrong-sign modes. Motivated by this observation, we propose a weaker version of the KSW criterion, which we call the spectral KSW criterion. Its purpose is to characterize when one-loop corrections around complex gravitational saddles are well defined.

We provide the conditions for complete asymptotic freedom for chiral gauge theories including scalars, as motivated by grand unified models. These are generalised Georgi-Glashow and Bars-Yankielowicz theories that feature a scalar field transforming either in the fundamental or in the adjoint of the gauge group. In both scenarios, we consider the addition of multiple chiral fermion families. We systematically analyse the interplay between gauge, Yukawa, and quartic couplings required for all interactions to remain asymptotically free at short distances. We find that for both scalar representations, complete asymptotic free models can be obtained for a specific number of colours and multiplicity of vector-like and chiral families.

We study invariants of bosonic and fermionic (Grassmann-valued) matrices under the adjoint action of $U(N)$, weighted by the fermion number. Such models naturally appear as the supersymmetric indices of supersymmetric gauge theories and are captured by $U(N)$ matrix models. We discuss two features of the fermionic models that are qualitatively different from bosonic models. Firstly, the $2N^\text{th}$ power of a Grassmann matrix vanishes, which gives rise to many new trace relations. Secondly, trace relations in models involving fermions could cause an increase in the supersymmetric index as $N$ decreases, in contrast with purely bosonic models. We focus on a simple model involving one fermion and one derivative that corresponds to a $\frac14$-BPS supersymmetric index in $\mathcal{N}=4$ SYM theory, in which we find that the index is independent of $N$. We prove this rank-independence analytically, and experimentally study the cancellations between bosonic and fermionic trace relations that lead to it. Based on these observations, we make some conjectures on resulting algebraic structures, including the analogue of the polarized Cayley-Hamilton identities and the Second Fundamental Theorem of invariants in the presence of Grassmann matrices. Finally, we present various (smooth and singular) limits of the most general supersymmetric index in $\mathcal{N}=4$ SYM theory, and study some patterns in their behavior as a function of $N$.

The eigenstate thermalization hypothesis provides a framework for understanding thermalization in isolated quantum many-body systems by characterizing statistical properties of local observables in energy eigenstates. Here we demonstrate that distributions of matrix elements in macroscopic systems may depend not only on the macrostate parameters, such as the densities of local conserved charges, but generally also on the properties of ensembles used in sampling eigenstates. To this end, we depart from the conventional analysis of microcanonical windows and consider statistical ensembles with an adjustable scale parameter prescribing the magnitude of charge fluctuations. We specifically consider an integrable field theory that permits efficient numerical sampling of matrix elements and reliable extrapolation to the thermodynamic limit. Moreover, in this system, we identify a class of states that enables explicit closed-form computation of the suppression rate of matrix elements. Our findings reveal an underlying multiscale structure of matrix elements captured by a non-analytic fluctuation-scale dependence of algebraic exponents governing their statistical properties.

The chiral magnetic effect (CME), arising from the chiral anomaly and enabling a mutual conversion between magnetic topology and fermionic chirality, is a key mechanism in magnetar field evolution. Previous work by Dehman & Pons (2025) demonstrated that the CME can efficiently generate dipolar fields $\left(B_{\rm dip} \gtrsim 10^{14}\,\mathrm{G}\right)$, consistent with magnetar timing measurements, provided that the initial magnetic field carries net helicity. However, whether neutron stars are born with magnetic helicity remains uncertain. In this work, we investigate the CME across a range of initial helicity configurations, including non-helical initial conditions. We find that the CME efficiently generates magnetar-strength dipoles on timescales of decades, independently of the initial helicity content. The chiral instability is driven by localized helical structures that induce a residual chiral asymmetry and is primarily governed by the maximum chiral chemical potential, requiring $\mu_5^{\rm max} \gtrsim \mathrm{few}\times10^{-11}\,\mathrm{MeV}$ for onset in the magnetar regime. Our results further show that these dipoles may either remain stable and subsequently evolve through standard Ohmic decay, or become unstable if they acquire sufficient helicity, in which case they decay through the chiral anomaly, transferring energy to less helical modes. This outcome depends sensitively on the initial helicity distribution. These findings extend the applicability of the CME to more realistic magnetic-field configurations and underscore the importance of the helicity distribution at birth, a quantity that remains poorly constrained in newborn neutron stars, yet is crucial for determining their magnetic evolution and the emergence of magnetars.

We studied the constrained Hamiltonian formulation of a supersymmetric Korteweg-de Vries (KdV) equation, which is observed to be a constrained system similar to its classical version. We found a nontrivial Lagrangian description, where we select $a=2$ for the free parameter $a$ in the supersymmetric extension. The corresponding degenerate Lagrangian requires an exclusive consideration and the utilization of the Dirac-Bergmann algorithm. We explicitly determined the full set of primary and secondary constraints and constructed the total Hamiltonian governing the dynamics of the system. In this analysis, in addition to a nontrivial constraint involving the fermionic fields, the consistency conditions give rise to a nonlocal contribution to the Hamiltonian density. This highlights a distinctive feature of this supersymmetric extension. We showed that the resulting Hamilton equations of motion reproduce the supersymmetric KdV system in the component form. Finally, we derived a compact superspace representation of the Hamiltonian and demonstrated its consistency with the component-level formulation.

We consider three-dimensional Einstein gravity in Euclidean signature with a finite boundary of torus topology endowed with an induced metric of fixed conformal class and a constant trace of extrinsic curvature $K$. For vanishing, positive, and negative cosmological constant $\Lambda$, we analytically determine boundaries enclosing different patches of locally flat, de Sitter (dS$_3$), and Anti-de Sitter (AdS$_3$) spaces. We find solutions that depend non-trivially on either cycle of the torus, noting that some of them exhibit self-intersections. Adapting the Gibbons-Hawking prescription of interpreting the Euclidean gravitational path integral as a thermal partition function, we explore the rich semi-classical thermodynamic phase space of the problem. While most saddles are found to be either thermally unstable or metastable compared to those with uniform boundaries, we find inhomogeneous solutions that are thermodynamically favourable in the case of $\Lambda < 0$ and $2<K|\Lambda|^{-1/2}<3/\sqrt{2}$. Moreover, for all values of $\Lambda$, there exist patches of space with a non-contractible thermal circle and a macroscopic entropy. We comment on a recasting of the problem in terms of classical strings. We further analyse the problem in both the stretched dS$_3$ horizon limit and the AdS$_3$ boundary limit. Using the results for the latter, we find limitations to a proposed holographic dual that consists of a two-dimensional holographic conformal field theory coupled to timelike Liouville, deformed by a marginal $T\bar{T}$-type operator.

In his 1976 paper "More about the massive Schwinger model", Coleman introduced $1+1$-dimensional Quantum Electrodynamics coupled to two charged massive fermions. By applying Abelian bosonization, he elucidated much of the physics of this two-flavor Schwinger model, but he listed three puzzles at the end of his paper. We present new analytical and numerical calculations to solve Coleman's three puzzles and thereby deepen our understanding of this model. These puzzles pertain to the theory with equal fermion masses at $\theta = 0$ and at $\theta = \pi$, as well as the size of isospin-breaking effects when the fermion masses are unequal. For the puzzle at $\theta = \pi$, the solution is related to the structure of the zero-temperature phase diagram arXiv:2305.04437: for equal fermion masses $m$, the model exhibits spontaneous breaking of charge conjugation symmetry and absence of confinement for any value of the gauge coupling $g$, so that there is a smooth interpolation from weak to strong coupling. Using two-loop Renormalization Group and integrability methods, we show that the mass gap behaves as $\sim m e^{-0.111 g^2/m^2}$ in the strong coupling regime $m\ll g$. Our numerical results using the lattice Hamiltonian are in good agreement with this behavior. For the puzzle at $\theta = 0$, the solution is related to a level crossing between two isosinglet particles with different discrete quantum numbers; we demonstrate the necessity of such a crossing by comparing integrability and weak coupling calculations, and we also exhibit the crossing numerically. Finally, we provide a new estimate for the size of isospin-breaking effects caused by different fermion masses at strong coupling.