Search

Although the nuclear fusion process has received a great deal of attention in recent years, the amount of mathematical analysis that supports the stability of the system seems to be relatively insufficient. This paper deals with the mathematical analysis of the magnetic confinement of the plasma via kinetic equations. We prove the global wellposedness of the Vlasov-Maxwell system in a two-dimensional annulus when a huge (but finite-in-time) external magnetic potential is imposed near the boundary. We assume that the solution is axisymmetric. The authors hope that this work is a step towards a more generalized work on the three-dimensional Tokamak structure. The highlight of this work is the physical assumptions on the external magnetic potential well which remains finite within a finite time interval and from that, we prove that the plasma never touches the boundary. In addition, we provide a sufficient condition on the magnitude of the external magnetic potential to guarantee that the plasma is confined in an annulus of the desired thickness which is slightly larger than the initial support. Our method uses the cylindrical coordinate forms of the Vlasov-Maxwell system., Comment: Accepted in a KRM special volume in the memory of Bob Glassey; 38 pages

Obtaining observational constraints on the role of turbulent effects for the solar dynamo is a difficult, yet crucial, task. Without such knowledge, the full picture of the operation mechanism of the solar dynamo cannot be formed. The magnetic helicity spectrum provides important information about the $\alpha$ effect. Here we demonstrate a formalism in spherical geometry to infer magnetic helicity spectra directly from observations of the magnetic field, taking into account the sign change of magnetic helicity across the Sun's equator. Using an angular correlation function of the magnetic field, we develop a method to infer spectra for magnetic energy and helicity. The retrieval of the latter relies on a fundamental definition of helicity in terms of linkage of magnetic flux. We apply the two-scale approach, previously used in Cartesian geometry, to spherical geometry for systems where a sign reversal of helicity is expected across the equator at both small and large scales. We test the method by applying it to an analytical model of a fully helical field, and to magneto-hydrodynamic simulations of a turbulent dynamo. The helicity spectra computed from the vector potential available in the models are in excellent agreement to the spectra computed solely from the magnetic field using our method. In a next test, we use our method to obtain the helicity spectrum from a synoptic magnetic field map corresponding to a Carrington rotation. We observe clear signs of a bihelical spectrum of magnetic helicity. Our formalism makes it possible to infer magnetic helicity in spherical geometry, without the necessity of computing the magnetic vector potential. This has the advantage of being gauge invariant. It has many applications in solar and stellar observations, but can also be used to analyze global magnetoconvection models of stars and compare them with observations., Comment: 10 pages, 6 figures, 3 appendices, submitted to Astronomy and Astrophysics

We study a class of fractional parabolic equations involving a time-dependent magnetic potential and formulate the corresponding inverse problem. We determine both the magnetic potential and the electric potential from the exterior partial measurements of the Dirichlet-to-Neumann map.

In this study we test 30 variants of 5 physical scaling laws that describe different aspects of solar flares. We express scaling laws in terms of the magnetic potential field energy $E_p$, the mean potential field strength $B_p$, the free energy $E_{free}$, the dissipated magnetic flare energy $E_{diss}$, the mean loop length scale $L$, the mean helically twisted flux tube radius $R$, the sunspot radius $r$, the emission measure-weighted flare temperature $T_w$, the electron density $n_e$, and the total emission measure $EM$, measured from a data set of $\lapprox 400$ GOES M- and X-class flare events. The 5 categories of physical scaling laws include (i) a scaling law of the potential-field energy, (ii) a scaling law for helical twisting, (iii) a scaling law for Petschek-type magnetic reconnection, (iv) the Rosner-Tucker-Vaiana scaling law, and (v) the Shibata-Yokoyama scaling law. We test the self-consistency of these theoretical scaling laws with observed parameters by requiring two conditions: a cross-corrleation coefficient of CCC$>$0.5 between the observed and theoretically predicted scaling laws, and a linear regression fit with a slope of $\alpha \approx 1$. With these two criteria we find that 10 out of the 30 tested scaling law variants are consistent with the observed data, which strongly corroborates the existence and validity of the tested flare scaling laws., Comment: 9 Figures

We study a class of fractional semilinear elliptic equations and formulate the corresponding Calder\'on problem. We determine the nonlinearity from the exterior partial measurements of the Dirichlet-to-Neumann map by using first order linearization and the Runge approximation property., Comment: 14 pages

Super-resolution of the Lie-Trotter splitting ($S_1$) and Strang splitting ($S_2$) is rigorously analyzed for the nonlinear Dirac equation without external magnetic potentials in the nonrelativistic regime with a small parameter $0, Comment: 26 pages, 2 figures and 8 tables

A program has been designed for accurately generating a potential magnetic field on a staggered grid by extrapolating the magnetic field normal to the photospheric surface. The code first calculates a magnetic potential using the Green's function method and then uses a finite differencing scheme to calculate the magnetic field from the potential. A new finite differencing formula was derived which accounts for grid staggering, it is shown that this formula gives a numerical approximation that is closest to the real potential field. It is also shown that extending the region over which normal photospheric field is specified can improve the accuracy of the potential field produced. The program is a FORTRAN 90 code that can be used to generate potential magnetic field inputs for Lare3d and other MHD solvers that use a staggered grid for magnetic field components. The program can be parallelised to run quickly over multiple computing cores. The code and supporting description are provided in the appendices and can also be found at https://github.com/calboo/Potential_Field_Calculator., Comment: Revision of some typos spotted whilst writing thesis

We study the problem of unique recovery of a non-smooth one-form $\mathcal A$ and a scalar function $q$ from the Dirichlet to Neumann map, $\Lambda_{\mathcal A,q}$, of a hyperbolic equation on a Riemannian manifold $(M,g)$. We prove uniqueness of the one-form $\mathcal A$ up to the natural gauge, under weak regularity conditions on $\mathcal A,q$ and under the assumption that $(M,g)$ is simple. Under an additional regularity assumption, we also derive uniqueness of the scalar function $q$. The proof is based on the geometric optic construction and inversion of the light ray transform extended as a Fourier Integral Operator to non-smooth parameters and functions.

This paper is devoted to the symmetry and symmetry breaking properties of a two-dimensional magnetic Schrödinger operator involving an Aharonov–Bohm magnetic vector potential. We investigate the symmetry properties of the optimal potential for the corresponding magnetic Keller–Lieb–Thirring inequality. We prove that this potential is radially symmetric if the intensity of the magnetic field is below an explicit threshold, while symmetry is broken above a second threshold corresponding to a higher magnetic field. The method relies on the study of the magnetic kinetic energy of the wave function and amounts to study the symmetry properties of the optimal functions in a magnetic Hardy–Sobolev interpolation inequality. We give a quantified range of symmetry by a non-perturbative method. To establish the symmetry breaking range, we exploit the coupling of the phase and of the modulus and also obtain a quantitative result., SCOPUS: ar.j, info:eu-repo/semantics/published

We construct black hole solutions to the leading order of string effective action in five dimensions with the source given by dilaton and magnetically charged antisymmetric gauge B-field. Presence of the considered B-field leads to the unusual asymptotic behavior of the solutions which are neither asymptotically flat nor asymptotically (A)dS. We consider the three-dimensional space part to correspond to the Bianchi classes and so the horizons of these topological black hole solutions are modeled by seven homogeneous Thurston’s geometries of $$E^3$$E3, $$S^3$$S3, $$H^3$$H3, $$ H^2 \times E^1$$H2×E1, $$\widetilde{{SL_2R}}$$SL2R~, nilgeometry, and solvegeometry. Calculating the quasi-local mass, temperature, entropy, dilaton charge, and magnetic potential, we show that the first law of black hole thermodynamics is satisfied by these quantities and the dilaton charge is not independent of mass and magnetic charge. Furthermore, for Bianchi type V, the T-dual black hole solution is obtained which carries no charge associated with B-field and the entropy turns to be invariant under the T-duality.

We propose a novel approach to reconstruct the surface magnetic helicity density on the Sun or sun-like stars. The magnetic vector potential is determined via decomposition of vector magnetic field measurements into toroidal and poloidal components. The method is verified using data from a non-axisymmetric dynamo model. We apply the method to vector field synoptic maps from Helioseismic and Magnetic Imager (HMI) onboard of Solar Dynamics Observatory (SDO) to study evolution of the magnetic helicity density during solar cycle 24. It is found that the mean helicity density of the non-axisymmetric magnetic field of the Sun evolves in a way which is similar to that reported for the current helicity density of the solar active regions. It has predominantly the negative sign in the northern hemisphere, and it is positive in the southern hemisphere. Also, the hemispheric helicity rule for the non-axisymmetric magnetic field showed the sign inversion at the end of cycle 24. Evolution of magnetic helicity density of large-scale axisymmetric magnetic field is different from that expected in dynamo theory. On one hand, the mean large- and small-scale components of magnetic helicity density display the hemispheric helicity rule of opposite sign at the beginning of cycle 24. However, later in the cycle, the two helicities exhibit the same sign in contrast with the theoretical expectations., Comment: 8 pages 5 figures

This paper is devoted to the study of fractional Schrödinger-Poisson type equations with magnetic field of the type $$\varepsilon^{2s}(-\Delta)_{A/\varepsilon}^{s}u + V(x)u + {\rm e}^{-2t}(\vert x \vert^{2t-3} \ast \vert u\vert ^{2})u = f(\vert u \vert^{2})u \quad \hbox{in} \ \open{R}^{3},$$ where ε > 0 is a parameter, s, t ∈ (0, 1) are such that 2s+2t>3, A:ℝ3 → ℝ3 is a smooth magnetic potential, (−Δ)As is the fractional magnetic Laplacian, V:ℝ3 → ℝ is a continuous electric potential and f:ℝ → ℝ is a C1 subcritical nonlinear term. Using variational methods, we obtain the existence, multiplicity and concentration of nontrivial solutions for e > 0 small enough.

This paper concerns mathematical theory of Meissner states of a bulk superconductor of type $I\!I$, which occupies a bounded domain $\Omega$ in $\Bbb R^3$ and is subjected to an applied magnetic field below the critical field $H_{S}$. A Meissner state is described by a solution $(f,\mathbf A)$ of a nonlinear partial differential system called Meissner system, where $f$ is a positive function on $\Omega$ which is equal to the modulus of the order parameter, and $\mathbf A$ is the magnetic potential defined on the entire space such that the inner trace of the normal component on the domain boundary $\partial\Omega$ vanishes. Such a solution is called a Meissner solution. Various properties of the Meissner solutions are examined, including regularity, classification and asymptotic behavior for large value of the Ginzburg-Landau parameter $\kappa$. It is shown that the Meissner solution is smooth in $\Omega$, however the regularity of the magnetic potential outside $\Omega$ can be rather poor. This observation leads to the ides of decomposition of the Meissner system into two problems, a boundary value problem in $\Omega$ and an exterior problem outside of $\Omega$. We show that the solutions of the boundary value problem with fixed boundary data converges uniformly on $\Omega$ as $\kappa$ tends to $\infty$, where the limit field of the magnetic potential is a solution of a nonlinear curl system. This indicates that, the magnetic potential part $\mathbf A$ of the solution $(f,\mathbf A)$ of the Meissner system, which has same tangential component of $curl \mathbf A$ on $\partial\Omega$, converges to a solution of the curl system as $\kappa$ increases to infinity, which verifies that the curl system is indeed the correct limit of the Meissner system in the case of three dimensions., Comment: Published in: J. Elliptic and Parabolic Equations, {\bf 4} (2) (2018), 441-523. https://link.springer.com/article/10.1007/s41808-018-0027-0 Springer Nature Sharedit initiative link https://rdcu.be/bbM8F

Clebsch potential gauge field theory for magnetohydrodynamics is developed based in part on the theory of Calkin (1963). It is shown how the polarization vector ${\bf P}$ in Calkin's approach, naturally arises from the Lagrange multiplier constraint equation for Faraday's equation for the magnetic induction ${\bf B}$, or alternatively from the magnetic vector potential form of Faraday's equation. Gauss's equation, (divergence of ${\bf B}$ is zero), is incorporated in the variational principle by means of a Lagrange multiplier constraint. Noether's theorem, coupled with the gauge symmetries is used to derive the conservation laws for (a)\ magnetic helicity (b)\ cross helicity, (c) fluid helicity for non-magnetized fluids, and (d) a class of conservation laws associated with curl and divergence equations, which applies to Faraday's equation and Gauss's equation. The magnetic helicity conservation law is due to a gauge symmetry in MHD and not due to a fluid relabelling symmetry. The analysis is carried out for the general case of a non-barotropic gas, in which the gas pressure and internal energy density depend on both the entropy $S$ and the gas density $\rho$. The cross helicity and fluid helicity conservation laws in the non-barotropic case, are nonlocal conservation laws, that reduce to local conservation laws for the case of a barotropic gas. The connections between gauge symmetries, Clebsch potentials and Casimirs are developed. It is shown that the gauge symmetry functionals in the work of Henyey (1982) satisfy the Casimir determining equations., Comment: 38 pages

Magnetic reconnection is a ubiquitous phenomenon in turbulent plasmas. It is an important part of the turbulent dynamics and heating of space and astrophysical plasmas. We examine the statistics of magnetic reconnection using a quantitative local analysis of the magnetic vector potential, previously used in magnetohydrodynamics simulations, and now employed to fully kinetic particle-in-cell (PIC) simulations. Different ways of reducing the particle noise for analysis purposes, including multiple smoothing techniques, are explored. We find that a Fourier filter applied at the Debye scale is an optimal choice for analyzing PIC data. Finally, we find a broader distribution of normalized reconnection rates compared to the MHD limit with rates as large as 0.5 but with an average of approximately 0.1.

Orbital advection is a significant bottleneck in disk simulations, and a particularly tricky one when used in connection with magnetohydrodynamics. We have developed an orbital advection algorithm suitable for the induction equation with magnetic potential. The electromotive force is split into advection and shear term, and we find that we do not need an advective gauge as solving the orbital advection implicitly precludes the shear term from canceling the advection term. We prove and demonstrate the third order in time accuracy of the scheme. The algorithm is also suited for non-magnetic problems. Benchmarked results of (hydrodynamical) planet-disk interaction and the of the magnetorotational instability are reproduced. We include detailed descriptions of the construction and selection of stabilizing dissipations (or high frequency filters) needed to generate practical results. The scheme is self-consistent, accurate, and elegant in its simplicity, making it particularly efficiently in straightforward finite-difference methods. As a result of the work, the algorithm is incorporated in the public version of the Pencil Code, where it can be used by the community., Comment: 11 pages + 4 appendix; 8+1 figures. Accepted by AJ; in press

A new class of black hole solutions of the five dimensional minimal gauged supergravity is presented. They are characterized by the mass, the electric charge, two equal magnitude angular momenta and the magnitude of the magnetic potential at infinity. These black holes possess a horizon of spherical topology; however, both the horizon and the sphere at infinity can be arbitrarily squashed, with nonextremal solutions interpolating between black strings and black branes. A particular set of extremal configurations corresponds to a new one-parameter family of supersymmetric black holes. While their conserved charges are determined by the squashing of the sphere at infinity, these supersymmetric solutions possess the same horizon geometry., Comment: 5 pages, 4 figures; v3: typos corrected, refs added; matches published version

In this paper we compute all the smooth solutions to the Hamilton-Jacobi equation associated with the horocycle flow. This can be seen as the Euler-Lagrange flow (restricted to the energy level set $E^{-1}(\frac 12)$) defined by the Tonelli Lagrangian $L:T\mathbb H\rightarrow \mathbb R$ given by (hyperbolic) kinetic energy plus the standard magnetic potential. The method we use is to look at Lagrangian graphs that are contained in the level set $\{H=\frac 12\}$, where $H:T^*\mathbb H\rightarrow \mathbb R$ denotes the Hamiltonian dual to $L$., Comment: Comments welcome

We study the motion of a particle in a particular magnetic field configuration both classically and quantum mechanically. For flux-free radially symmetric magnetic fields defined on circular regions, we establish that particle escape speeds depend, classically, on a gauge-fixed magnetic vector potential, and demonstrate some trajectories associated with this special type of magnetic field. Then we show that some of the geometric features of the classical trajectory (perpendicular exit from the field region, trapped and escape behavior) are reproduced quantum mechanically using a numerical method that extends the norm-preserving Crank-Nicolson method to problems involving magnetic fields. While there are similarities between the classical trajectory and the position expectation value of the quantum mechanical solution, there are also differences, and we demonstrate some of these., Comment: To appear, American Journal of Physics

This manuscript considers the main features of the nonlinear dynamics of multiple irregular celestial body systems. The gravitational potential, static electric potential, and magnetic potential are considered. Based on the three established potentials, we show that three conservative values exist for this system, including a Jacobi integral. The equilibrium conditions for the system are derived and their stability analyzed. The equilibrium conditions of a celestial system comprised of n irregular bodies are reduced to 12n minus 9 equations. The dynamical results are applied to simulate the motion of multiple-asteroid systems. The simulation is useful for the study of the stability of multiple irregular celestial body systems and for the design of spacecraft orbits to triple asteroid systems discovered in the solar system. The dynamical configurations of the five triple-asteroid systems 45 Eugenia, 87 Sylvia, 93 Minerva, 216 Kleopatra, and 136617 1994CC, and the six-body system 134340 Pluto are calculated and analyzed., Comment: 60 pages, 20 figures, The online version of this article (doi:10.1007/s10509-016-2884-z) contains supplementary material, which is available to authorized users. 6 videos, 5 rar files, 1 pdf file of figures

We combine experiment and theory to investigate the diffusive and subdiffusive dynamics of paramagnetic colloids driven above a two-state flashing potential. The magnetic potential was realized by periodically modulating the stray field of a magnetic bubble lattice in a uniaxial ferrite garnet film. At large amplitudes of the driving field, the dynamics of particles resembles an ordinary random walk with a frequency-dependent diffusion coefficient. However, subdiffusive and oscillatory dynamics at short time scales is observed when decreasing the amplitude. We present a persistent random walk model to elucidate the underlying mechanism of motion, and perform numerical simulations to demonstrate that the anomalous motion originates from the dynamic disorder in the structure of the magnetic lattice, induced by slightly irregular shape of bubbles., Comment: 9 pages, 8 figures

We consider the semi-classical Dirac operator coupled to a magnetic potential on a large class of manifolds including all metric contact manifolds. We prove a sharp local Weyl law and a bound on its eta invariant. In the absence of a Fourier integral parametrix, the method relies on the use of almost analytic continuations combined with the Birkhoff normal form and local index theory.

We propose a high-order finite difference weighted ENO (WENO) method for the ideal magnetohydrodynamics (MHD) equations. The proposed method is single-stage (i.e., it has no internal stages to store), single-step (i.e., it has no time history that needs to be stored), maintains a discrete divergence-free condition on the magnetic field, and has the capacity to preserve the positivity of the density and pressure. To accomplish this, we use a Taylor discretization of the Picard integral formulation (PIF) of the finite difference WENO method proposed in Christlieb et al. (2015) 23, where the focus is on a high-order discretization of the fluxes (as opposed to the conserved variables). We use the version where fluxes are expanded to third-order accuracy in time, and for the fluid variables space is discretized using the classical fifth-order finite difference WENO discretization. We use constrained transport in order to obtain divergence-free magnetic fields, which means that we simultaneously evolve the magnetohydrodynamic (that has an evolution equation for the magnetic field) and magnetic potential equations alongside each other, and set the magnetic field to be the (discrete) curl of the magnetic potential after each time step. In this work, we compute these derivatives to fourth-order accuracy. In order to retain a single-stage, single-step method, we develop a novel Lax-Wendroff discretization for the evolution of the magnetic potential, where we start with technology used for Hamilton-Jacobi equations in order to construct a non-oscillatory magnetic field. The end result is an algorithm that is similar to our previous work Christlieb et al. (2014) 8, but this time the time stepping is replaced through a Taylor method with the addition of a positivity-preserving limiter. Finally, positivity preservation is realized by introducing a parameterized flux limiter that considers a linear combination of high and low-order numerical fluxes. The choice of the free parameter is then given in such a way that the fluxes are limited towards the low-order solver until positivity is attained. Given the lack of additional degrees of freedom in the system, this positivity limiter lacks energy conservation where the limiter turns on. However, this ingredient can be dropped for problems where the pressure does not become negative. We present two and three dimensional numerical results for several standard test problems including a smooth Alfven wave (to verify formal order of accuracy), shock tube problems (to test the shock-capturing ability of the scheme), Orszag-Tang, and cloud shock interactions. These results assert the robustness and verify the high-order of accuracy of the proposed scheme.

We prove global existence and uniqueness of weak solutions for the time-dependent Ginzburg–Landau equations in a three-dimensional curved polyhedron which is not necessarily convex, where the gradient of the magnetic potential may not be square integrable. Preceding analyses all required the gradient of the solution to be square integrable, which is only true in convex or smooth domains.

By the introduction of a generalized magnetic vector potential, which contains the contributions of both the magnetic and electric parts, and the use of the Ampere's law within the quasistatic approximation as the state equation, the partial differential equations for governing the electromagnetic properties of superconductors as well as the surrounding coolant were established and numerically discretized by resorting to the finite-element technique and finite-difference scheme, respectively, in the spatial and temporal domain. In conjunction with an analytic method to calculate the magnetic field generated by permanent magnet, we compiled a numerical tool for performing an intricate study of the mutual effect among the superconducting constituents in a superconducting levitation system with translational symmetry. Taking a superconducting unit with three constituents inside as a practice, we simulated the electromagnetic responses of this unit while moving in the nonuniform magnetic field generated by permanent magnet guideway and, identified the influences of the mutual effect on the levitation force as well as on the distributions of the magnetic flux density, the supercurrent density, and the levitation force density by comparing to an envisaged reference, one constituent was simulated with all the rest absent to remove the mutual effect. The insights attained by the present study, mostly being inaccessible from the experiments, are aimed to provide useful implications for the design of a superconducting levitation system for the transit and analogous purposes, which usually employ multiple superconductors to achieve the desired capability., Comment: 22 pages, 13 figures

Magnetic fields and flows in coronal structures, for example, in gradual phases in flares, can be described by 2D and 3D magnetohydrostatic (MHS) and steady magnetohydrodynamic (MHD) equilibria. Within a physically simplified, but exact mathematical model, we study the electric currents and corresponding electric fields generated by shear flows. Starting from exact and analytically calculated magnetic potential fields, we solveid the nonlinear MHD equations self-consistently. By applying a magnetic shear flow and assuming a nonideal MHD environment, we calculated an electric field via Faraday's law. The formal solution for the electromagnetic field allowed us to compute an expression of an effective resistivity similar to the collisionless Speiser resistivity. We find that the electric field can be highly spatially structured, or in other words, filamented. The electric field component parallel to the magnetic field is the dominant component and is high where the resistivity has a maximum. The electric field is a potential field, therefore, the highest energy gain of the particles can be directly derived from the corresponding voltage. In our example of a coronal post-flare scenario we obtain electron energies of tens of keV, which are on the same order of magnitude as found observationally. This energy serves as a source for heating and acceleration of particles., Comment: 11 pages, 9 figures, accepted to Astronomy and Astrophysics

We introduce a general class of long-range magnetic potentials and derive high velocity limits for the corresponding scattering operators in quantum mechanics, in the case of two dimensions. We analyze the high velocity limits that we obtain in the presence of an obstacle and we uniquely reconstruct from them the electric potential and the magnetic field outside the obstacle, that are accessible to the particles. We additionally reconstruct the inaccessible fluxes (magnetic fluxes produced by fields inside the obstacle) modulo 2π, which give a proof of the Aharonov–Bohm effect. For every magnetic potential A in our class, we prove that its behavior at infinity [Formula: see text] can be characterized in a natural way; we call it the long-range part of the magnetic potential. Under very general assumptions, we prove that [Formula: see text] can be uniquely reconstructed for every [Formula: see text]. We characterize properties of the support of the magnetic field outside the obstacle that permit us to uniquely reconstruct [Formula: see text] either for all [Formula: see text] or for [Formula: see text] in a subset of 𝕊1. We also give a wide class of magnetic fields outside the obstacle allowing us to uniquely reconstruct the total magnetic flux (and [Formula: see text] for all [Formula: see text]). This is relevant because, as it is well-known, in general the scattering operator (even if it is known for all velocities or energies) does not define uniquely the total magnetic flux (and [Formula: see text]). We analyze additionally injectivity (i.e. uniqueness without giving a method for reconstruction) of the high velocity limits of the scattering operator with respect to [Formula: see text]. Assuming that the magnetic field outside the obstacle is not identically zero, we provide a class of magnetic potentials for which injectivity is valid.

Starting with a mathematical boundary value problem for the magnetic vector potential in an axisymmetric cylindrical coordinate system, we derive a general solution for any arbitrary current distribution using the method of Green's functions. We use this to derive an analytic form for an hourglass magnetic field pattern created by electrical currents that are concentrated near (but not confined within) the equatorial plane of a cylindrical coordinate system. Our solution is not characterized by a cusp at the equatorial plane, as in previous solutions based on a current sheet. The pattern we derive provides a very good fit to hourglass magnetic field patterns emerging from three-dimensional numerical simulations of core formation, and can in principle be used for source-fitting of observed magnetic hourglass patterns., Comment: 6 pages, 4 figures, to appear in ApJ (20 April 2013)

In this paper we study inverse boundary value problems with partial data for the bi-harmonic operator with first order perturbation. We consider two types of subsets of $\mathbb{R}^{n}(n\geq 3)$, one is an infinite slab, the other is a bounded domain. In the case of a slab, we show that, from Dirichlet and Neumann data given either on the different boundary hyperplanes of the slab or on the same boundary hyperplane, one can uniquely determine the magnetic potential and the electric potential. In the case of a bounded domain, we show the unique determination of the magnetic potential and the electric potential from partial Dirichlet and Neumann data under two different assumptions. The first assumption is that the magnetic and electric potentials are known in a neighborhood of the boundary, in this situation we obtain the uniqueness result when the Dirichlet and Neumann data are only given on two arbitrary open subsets of the boundary. The second assumption is that the Dirichlet and Neumann data are known on the same part of the boundary whose complement is a part of a hyperplane, we also establish the unique determination result in this local data case., Comment: 32 pages. Comments welcome. arXiv admin note: text overlap with arXiv:1104.0789 by other authors

Analytical models of axially symmetric thin disks of finite extension in presence of magnetic field are presented based on the well-known Morgan-Morgan solutions. The source of the magnetic field is cons\-tructed separating the equation corresponding to the Ampere's law of electrodynamics in spheroidal oblate coordinates. This produces two asso\-ciated Legendre equations of first order for the magnetic potential and hence that can be expressed as a series of associated Legendre functions of the same order. The discontinuity of its normal derivate across the disk allows us to interpret the source of the magnetic field as a ringlike current distribution extended on all the plane of the disk. We also study the circular speed curves or rotation curve for equatorial circular orbits of charged test particles both inside and outside the disk. The stability of the orbits is analyzed for radial perturbation using a extension of the Rayleigh criterion., Comment: 18 pages, 11 figures. Accepted for publication in Indian Journal of Physics

We report the first demonstration of an inductively coupled magnetic ring trap for cold atoms. A uniform, ac magnetic field is used to induce current in a copper ring, which creates an opposing magnetic field that is time-averaged to produce a smooth cylindrically symmetric ring trap of radius 5 mm. We use a laser-cooled atomic sample to characterise the loading efficiency and adiabaticity of the magnetic potential, achieving a vacuum-limited lifetime in the trap. This technique is suitable for creating scalable toroidal waveguides for applications in matterwave interferometry, offering long interaction times and large enclosed areas., Comment: 12 pages, 6 figures

A transfer matrix method is presented for solving the scattering problem for the quasi one-dimensional massless Dirac equation applied to graphene in the presence of an arbitrary inhomogeneous electric and perpendicular magnetic field. It is shown that parabolic cylindrical functions, which have previously been used in literature, become inaccurate at high incident energies and low magnetic fields. A series expansion technique is presented to circumvent this problem. An alternate method using asymptotic expressions is also discussed and the relative merits of the two methods are compared., Comment: Accepted for Publication in Modelling and Simulation in Materials Science and Engineering

We consider the semiclassical asymptotics of the sum of negative eigenvalues of the three-dimensional Pauli operator with an external potential and a self-generated magnetic field $B$. We also add the field energy $\beta \int B^2$ and we minimize over all magnetic fields. The parameter $\beta$ effectively determines the strength of the field. We consider the weak field regime with $\beta h^{2}\ge {const}>0$, where $h$ is the semiclassical parameter. For smooth potentials we prove that the semiclassical asymptotics of the total energy is given by the non-magnetic Weyl term to leading order with an error bound that is smaller by a factor $h^{1+\e}$, i.e. the subleading term vanishes. However, for potentials with a Coulomb singularity the subleading term does not vanish due to the non-semiclassical effect of the singularity. Combined with a multiscale technique, this refined estimate is used in the companion paper \cite{EFS3} to prove the second order Scott correction to the ground state energy of large atoms and molecules., Comment: Small typos corrected on Sep 24, 2011

We consider the unperturbed operator $H_0: = (-i \nabla - {\bf A})^2 + W$, self-adjoint in $L^2({\mathbb R}^2)$. Here ${\bf A}$ is a magnetic potential which generates a constant magnetic field $b>0$, and the edge potential $W = \bar{W}$ is a ${\mathcal T}$-periodic non constant bounded function depending only on the first coordinate $x \in {\mathbb R}$ of $(x,y) \in {\mathbb R}^2$. Then the spectrum $\sigma(H_0)$ of $H_0$ has a band structure, the band functions are $b {\mathcal T}$-periodic, and generically there are infinitely many open gaps in $\sigma(H_0)$. We establish explicit sufficient conditions which guarantee that a given band of $\sigma(H_0)$ has a positive length, and all the extremal points of the corresponding band function are non degenerate. Under these assumptions we consider the perturbed operators $H_{\pm} = H_0 \pm V$ where the electric potential $V \in L^{\infty}({\mathbb R}^2)$ is non-negative and decays at infinity. We investigate the asymptotic distribution of the discrete spectrum of $H_\pm$ in the spectral gaps of $H_0$. We introduce an effective Hamiltonian which governs the main asymptotic term; this Hamiltonian could be interpreted as a 1D Schroedinger operator with infinite-matrix-valued potential. Further, we restrict our attention on perturbations $V$ of compact support. We find that there are infinitely many discrete eigenvalues in any open gap in the spectrum of $\sigma(H_0)$, and the convergence of these eigenvalues to the corresponding spectral edge is asymptotically Gaussian., Comment: Lemma 2.1 added, the proof of Theorem 3.1 streamlined, typos corrected. 21 pages

International audience; We consider a Schrödinger operator H with a non-vanishing radial magnetic field B=dA and Dirichlet boundary conditions on the unit disk. We assume growth conditions on B near the boundary which guarantee in particular the compactness of the resolvent of this operator. Under some assumptions on an additional radial potential V the operator H + V has a discrete negative spectrum and we obtain an upper bound on the number of negative eigenvalues. As a consequence we get an upperbound of the number of eigenvalues of H smaller than any positive value, which involves the minimum of B and the square of the L^2 -norm of A(r)/r, where A(r) is the specific magnetic potential defined as the flux of the magnetic field through the disk of radius r centerde in the origin.

We prove optimal high-frequency resolvent estimates for perturbations of the Laplacian by large long-range magnetic and electric potentials in all dimensions $n\ge 3$. As an application, we prove dispersive estimates for the corresponding wave group in the case $n=3$., Comment: 20 pages, some improvements are made

Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e. when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential., Comment: 11 pages, 10 figures, submitted to Physics of Plasmas

This work examines the relation between the fractal properties of the photospheric magnetic patterns and those of the coronal magnetic fields in solar active regions.Our results show no correlation between the fractal dimensions of the observed 2d photospheric structures and the extrapolated unstable volumes in the corona, when nonlinear force-free extrapolation is used. This result is independent of efforts to (1) bring the photospheric magnetic fields closer to a nonlinear force-free equilibrium and (2) omit the lower part of the modeled magnetic field volume that is almost completely filled by unstable volumes. A significant correlation between the fractal dimensions of the photospheric and coronal magnetic features is only observed at the zero level (lower limit) of approximation of a current-free (potential) magnetic field extrapolation. We conclude that the complicated transition from photospheric non-force-free fields to coronal force-free ones hampers any direct correlation between the fractal dimensions of the 2d photospheric patterns and their 3d counterparts in the corona at the nonlinear force-free limit., Comment: 9 pages, 7 figures, 7 tables. Astronomy & Astrophysics (in press)

Using two- and three-dimensional hydromagnetic simulations for a range of different flows, including laminar and turbulent ones, it is shown that solutions expressing the field in terms of Euler potentials (EP) are in general incorrect if the EP are evolved with an artificial diffusion term. In three dimensions, standard methods using the magnetic vector potential are found to permit dynamo action when the EP give decaying solutions. With an imposed field, the EP method yields excessive power at small scales. This effect is more exaggerated in the dynamic case, suggesting an unrealistically reduced feedback from the Lorentz force. The EP approach agrees with standard methods only at early times when magnetic diffusivity did not have time to act. It is demonstrated that the usage of EP with even a small artificial magnetic diffusivity does not converge to a proper solution of hydromagnetic turbulence. The source of this disagreement is not connected with magnetic helicity or the three-dimensionality of the magnetic field, but is simply due to the fact that the nonlinear representation of the magnetic field in terms of EP that depend on the same coordinates is incompatible with the linear diffusion operator in the induction equation., Comment: 8 pages, 10 figures, MNRAS, accepted

The configuration-fixed deformation constrained relativistic mean field approach with time-odd component has been applied to investigate the ground-state properties of $^{33}$Mg with effective interaction PK1. The ground state of $^{33}$Mg has been found to be prolate deformed, $\beta_2=0.23$, with the odd neutron in $1/2[330]$ orbital and the energy -251.85 MeV which is close to the data -252.06 MeV. The magnetic moment $- 0.9134 \mu_\mathrm{N}$ is obtained with the effective electromagnetic current which well reproduces the data $- 0.7456 \mu_\mathrm{N}$ self-consistently without introducing any parameter. The energy splittings of time reversal conjugate states, the neutron current, the energy contribution from the nuclear magnetic potential, and the effect of core polarization are discussed in detail., Comment: 13 pages, 4 figures

In this paper we consider a magnetic analog of the double-Reissner-Nordstrom solution and construct the corresponding magnetic potential A_\varphi in the explicit form. The behavior of the resulting solution under the Harrison transformation then naturally singles out the asymmetric black diholes - configurations composed of two non-extreme black holes possessing unequal masses, and charges equal in magnitude but opposite in sign - as its most general subclass for which equilibrium of the black-hole constituents can be achieved with the aid of the external magnetic (or electric) field. We also generalize the double-Reissner-Nordstrom solution to the dilaton gravity with arbitrary dilaton coupling, yielding as the result the 4-dimensional double-Gibbons-Maeda spacetime. The study of some physical properties of the solutions obtained leads, in particular, to very simple formulas for the areas of the horizons and surface gravities., Comment: 18 pages, 1 figure; title changed, typos corrected; a considerably extended version which includes the discussion of the magnetostatic case and the explicit formula for the magnetic potential

This paper is concerned with small parameter asymptotics of magnetic quantum systems. In addition to a semiclassical parameter \eps, the case of small coupling $\lambda$ to the magnetic vector potential naturally occurs in this context. Magnetic Weyl calculus is adapted to incorporate both parameters, at least one of which needs to be small. Of particular interest is the expansion of the Weyl product which can be used to expand the product of operators in a small parameter, a technique which is prominent to obtain perturbation expansions. Three asymptotic expansions for the magnetic Weyl product of two H\"ormander class symbols are proven: (i) \eps \ll 1 and \lambda \ll 1, (ii) \eps \ll 1 and \lambda = 1 as well as (iii) \eps = 1 and \lambda \ll 1. Expansions (i) and (iii) are impossible to obtain with ordinary Weyl calculus. Furthermore, I relate results derived by ordinary Weyl calculus with those obtained with magnetic Weyl calculus by one- and two-parameter expansions. To show the power and versatility of magnetic Weyl calculus, I derive the semirelativistic Pauli equation as a scaling limit from the Dirac equation up to errors of 4th order in 1/c., Comment: 37 pages (shortened version, fixed typos)

Coronal magnetic field may be characterized by how its field lines interconnect regions of opposing photospheric flux -- its connectivity. Connectivity can be quantified as the net flux connecting pairs of opposing regions, once such regions are identified. One existing algorithm will partition a typical active region into a number of unipolar regions ranging from a few dozen to a few hundred, depending on algorithmic parameters. This work explores how the properties of the partitions depend on some algorithmic parameters, and how connectivity depends on the coarseness of partitioning for one particular active region magnetogram. We find the number of connections among them scales with the number of regions even as the number of possible connections scales with its square. There are several methods of generating a coronal field, even a potential field. The field may be computed inside conducting boundaries or over an infinite half-space. For computation of connectivity, the unipolar regions may be replaced by point sources or the exact magnetogram may be used as a lower boundary condition. Our investigation shows that the connectivities from these various fields differ only slightly -- no more than 15%. The greatest difference is between fields within conducting walls and those in the half-space. Their connectivities grow more different as finer partitioning creates more source regions. This also gives a quantitative means of establishing how far away conducting boundaries must be placed in order not to significantly affect the extrapolation. For identical outer boundaries, the use of point sources instead of the exact magnetogram makes a smaller difference in connectivity: typically 6% independent of the number of source regions.

Impulse formulations of Hall magnetohydrodynamic (MHD) equations are developed. The Lagrange invariance of a generalized ion magnetic helicity is established for Hall MHD. The physical implications of this Lagrange invariant are discussed. The discussion is then extended to compressible Hall MHD and a generalized ion magnetic potential helicity Lagrange invariant is established. The physical implications of the generalized ion magnetic potential helicity Lagrange invariant are shown to be the same, as to be expected, as those of the generalized ion magnetic helicity Lagrange invariant.

A proof is given of the vector identity proposed by Gubarev, Stodolsky and Zakarov that relates the volume integral of the square of a 3-vector field to non-local integrals of the curl and divergence of the field. The identity is applied to the case of the magnetic vector potential and magnetic field of a rotating charged shell. The latter provides a straightforward exercise in the use of the addition theorem of spherical harmonics., Comment: Based upon the paper in Am. J. Phys., but contains also a derivation of the central result that is valid for the scalar product of two different vector fields. Also, a proof is given that the general expression for the vector potential in the Coulomb gauge satisfies the central result. A discussion is given of the requirements on the gauge function. 12 pages pdf

We give a procedure for reconstructing a magnetic field and electric potential from boundary measurements given by the Dirichlet to Neumann map for the magnetic Schrodinger operator in R n , n ≥ 3. The magnetic potential is assumed to be continuous with L ∞ divergence and zero boundary values. The method is based on semiclassical pseudodifferential calculus and the construction of complex geometrical optics solutions in weighted Sobolev spaces.

The estimation of potential fields such as the gravitational or magnetic potential at the surface of a spherical planet from noisy observations taken at an altitude over an incomplete portion of the globe is a classic example of an ill-posed inverse problem. Here we show that the geodetic estimation problem has deep-seated connections to Slepian's spatiospectral localization problem on the sphere, which amounts to finding bandlimited spherical functions whose energy is optimally concentrated in some closed portion of the unit sphere. This allows us to formulate an alternative solution to the traditional damped least-squares spherical harmonic approach in geodesy, whereby the source field is now expanded in a truncated Slepian function basis set. We discuss the relative performance of both methods with regard to standard statistical measures as bias, variance and mean-square error, and pay special attention to the algorithmic efficiency of computing the Slepian functions on the region complementary to the axisymmetric polar gap characteristic of satellite surveys. The ease, speed, and accuracy of this new method makes the use of spherical Slepian functions in earth and planetary geodesy practical., Comment: 14 figures, submitted to the Geophysical Journal International

We write down scalar field theory and gauge theory on two-dimensional noncommutative spaces ${\cal M}$ with nonvanishing curvature and non-constant non-commutativity. Usual dynamics results upon taking the limit of ${\cal M}$ going to i) a commutative manifold ${\cal M}_0$ having nonvanishing curvature and ii) the noncommutative plane. Our procedure does not require introducing singular algebraic maps or frame fields. Rather, we exploit the K\"ahler structure in the limit i) and identify the symplectic two-form with the volume two-form. As an example, we take ${\cal M}$ to be the stereographically projected fuzzy sphere, and find magnetic monopole solutions to the noncommutative Maxwell equations. Although the magnetic charges are conserved, the classical theory does not require that they be quantized. The noncommutative gauge field strength transforms in the usual manner, but the same is not, in general, true for the associated potentials. We develop a perturbation scheme to obtain the expression for gauge transformations about limits i) and ii). We also obtain the lowest order Seiberg-Witten map to write down corrections to the commutative field equations and show that solutions to Maxwell theory on ${\cal M}_0$ are stable under inclusion of lowest order noncommutative corrections. The results are applied to the example of noncommutative AdS${}^2$., Comment: 27 pages

The metric ansatz is used to describe the gravitational field of a beam-pulse of spinning radiation (gyraton) in an arbitrary number of spacetime dimensions D. First we demonstrate that this metric belongs to the class of metrics for which all scalar invariants constructed from the curvature and its covariant derivatives vanish. Next, it is shown that the vacuum Einstein equations reduce to two linear problems in (D-2)-dimensional Euclidean space. The first is to find the static magnetic potential created by a point-like source. The second requires finding the electric potential created by a point-like source surrounded by given distribution of the electric charge. To obtain a generic gyraton-type solution of the vacuum Einstein equations it is sufficient to allow the coefficients in the corresponding harmonic decompositions of solutions of the linear problems to depend arbitrarily on retarded time and substitute the obtained expressions in the metric ansatz. We discuss properties of the solutions for relativistic gyratons and consider special examples., Comment: 11 pages

We provide an exact microscopic statistical treatment of particle and field correlations in a system of quantum charges in equilibrium with a classical radiation field. Using the Feynman-Kac-Ito representation of the Gibbs weight, the system of particles is mapped onto a collection of random charged wires. The field degrees of freedom can be integrated out, providing an effective pairwise magnetic potential. We then calculate the contribution of the transverse field coupling to the large-distance particle correlations. The asymptotics of the field correlations in the plasma are also exactly determined., Comment: 31 pages, 0 figures. PACS 05.30.-d, 05.40.-a, 11.10.Wx. Changes: Improved comparison with existing literature on field correlations. Added Concluding Remarks. References updated