Search

While boundary plasmas in present day tokamaks generally fall in a fluid regime, neutral species near the boundary often require kinetic models due to long mean-free-paths compared to characteristic spatial scales in the region. Monte-Carlo (MC) methods provide a complete, high-fidelity approach to solving kinetic models, and must be coupled to fluid plasma models to simulate the full plasma-neutrals system. The statistical nature of MC methods, however, prevents convergence of coupled fluid-kinetic simulations to an exact self-consistent steady-state. Moreover, this forces the use of explicit methods that can suffer from numerical errors and require huge computational resources. Correlated Monte-Carlo (CMC) methods are expected to alleviate these issues, but have historically enjoyed only mixed success. Here, a fully implicit method for coupled plasma-neutral systems is demonstrated in 1D using the UEDGE plasma code and a homemade CMC code. In particular, it is shown that ensuring the CMC method is a differentiable function of the background plasma is sufficient to employ a Jacobian-Free Newton-Krylov solver for implicit time steps. The convergence of the implicit coupling method is explored and compared with explicit coupling and uncorrelated methods. It is shown that ensuring differentiability by controlling random seeds in the MC is sufficient to achieve convergence, and that the use of implicit time-stepping methods has the potential for improved stability and runtimes over explicit coupling methods., Comment: 6 pages, 6 figures. Comments welcome!

Given a pair of $\mathbb{Z}_2$-harmonic spinors (resp. 1-forms) on closed Riemannian 3-manifolds $(Y_1, g_1)$ and $(Y_2,g_2)$, we construct $\mathbb{Z}_2$-harmonic spinors (resp. 1-forms) on the connected sum $Y_1 \# Y_2$ and the torus sum $Y_1 \cup_{T^2} Y_2$ using a gluing argument. The main tool in the proof is a parameterized version of the Nash-Moser implicit function theorem established by Donaldson and the second author. We use these results to construct an abundance of new examples of $\mathbb Z_2$-harmonic spinors and 1-forms. In particular, we prove that for every closed 3-manifold $Y$, there exist infinitely many $\mathbb{Z}_2$-harmonic spinors with singular sets representing infinitely many distinct isotopy classes of embedded links, strengthening an existence theorem of Doan-Walpuski. Moreover, combining this with previous results, our construction implies that if $b_1(Y) > 0$, there exist infinitely many $\mathrm{spin}^c$ structures on $Y$ such that the moduli space of solutions to the two-spinor Seiberg-Witten equations is non-empty and non-compact., Comment: 38 pages, comments welcome!

Given a $\mathbb Z_2$-harmonic spinor satisfying some genericity assumptions, this article constructs a 1-parameter family of two-spinor Seiberg-Witten monopoles converging to it after renormalization. The proof is a gluing construction beginning with model solutions on a neighborhood of the $\mathbb Z_2$-harmonic spinor's singular set. The gluing is complicated by the presence of an infinite-dimensional obstruction bundle for the singular limiting linearized operator. This difficulty is overcome by introducing a generalization of Donaldson's alternating method in which a deformation of the $\mathbb Z_2$-harmonic spinor's singular set is chosen at each stage of the alternating iteration to cancel the obstruction components., Comment: 67 pages, 1 Figure. Comments welcome!

This article studies a class of Dirac operators of the form $D_\varepsilon= D+\varepsilon^{-1}\mathcal A$, where $\mathcal A$ is a zeroth order perturbation vanishing on a subbundle. When $\mathcal A$ satisfies certain additional assumptions, solutions of the Dirac equation have a concentration property in the limit $\varepsilon\to 0$: components of the solution orthogonal to $\ker(\mathcal A)$ decay exponentially away from the locus $\mathcal Z$ where the rank of $\ker(\mathcal A)$ jumps up. These results are extended to a class of non-linear Dirac equations. This framework is then applied to study the compactness properties of moduli spaces of solutions to generalized Seiberg-Witten equations. In particular, it is shown that for sequences of solutions which converge weakly to a $\mathbb Z_2$-harmonic spinor, certain components of the solutions concentrate exponentially around the singular set of the $\mathbb Z_2$-harmonic spinor. Using these results, the weak convergence to $\mathbb Z_2$-harmonic spinors proved in existing convergence theorems is improved to $C^\infty_{loc}$., Comment: Comments welcome!

A $\mathbb Z_2$-harmonic spinor on a 3-manifold $Y$ is a solution of the Dirac equation on a bundle that is twisted around a submanifold $\mathcal Z$ of codimension 2 called the singular set. This article investigates the local structure of the universal moduli space of $\mathbb Z_2$-harmonic spinors over the space of parameters $(g,B)$ consisting of a metric and perturbation to the spin connection. The main result states that near a $\mathbb Z_2$-harmonic spinor with $\mathcal Z$ smooth, the universal moduli space projects to a codimension 1 submanifold in the space of parameters. The analysis is complicated by the presence of an infinite-dimensional obstruction bundle and a loss of regularity in the first variation of the Dirac operator with respect to deformations of the singular set $\mathcal Z$, necessitating the use of the Nash-Moser Implicit Function Theorem., Comment: Comments welcome!

Given a compact 3-manifold $Y$ and a $\mathbb Z_2$-harmonic spinor $(\mathcal Z_0, A_0,\Phi_0)$ with singular set $\mathcal Z_0$, this article constructs a family of local solutions to the two-spinor Seiberg-Witten equations parameterized by $\epsilon \in (0,\epsilon_0)$ on tubular neighborhoods of $\mathcal Z_0$. These solutions concentrate in the sense that the $L^2$-norm of the curvature near $\mathcal Z_0$ diverges as $\epsilon\to 0$, and after renormalization they converge locally to the original $\mathbb Z_2$-harmonic spinor. In a sequel to this article, these model solutions are used in a gluing construction showing that any $\mathbb Z_2$-harmonic spinor satisfying some mild assumptions arises as the limit of a family of two-spinor Seiberg-Witten solutions on $Y$., Comment: 104 pages, 3 figures

A stochastic differential equation that describes the dynamics of single-domain magnetic particles at any temperature is derived using a classical formalism. The deterministic terms recover existing theory and the stochastic process takes the form of a mean-reverting random walk. In the ferromagnetic state diffusion is predominantly angular and the relevant diffusion coefficient increases linearly with temperature before saturating at the Curie point ($T_c$). Diffusion in the macrospin magnitude, while vanishingly small at room temperature, increases sharply as the system approaches $T_c$. Beyond $T_c$, in the paramagnetic state, diffusion becomes isotropic and independent of temperature. The stochastic macrospin model agrees well with atomistic simulations., Comment: 5 pages, 2 figures