Structure and energy of a line vortex in a three-dimensional ordered Josephson medium

Two possible equilibrium configurations of line vortices in a three-dimensional ordered Josephson medium for any value of structural factor b are considered: the center of the vortex coincides with the center of one of the cells and the center of the vortex is on one of the contacts. Infinite sets of equations describing these configurations are derived. The infinite set can be made finite if currents away from the center are neglected. The assumption b = 0 is shown to be valid if pinning parameter I is less than 0.25. For I > 0.25, the structures and energies of both configurations of line isolated vortices are calculated throughout the range of structural factor b. As structural factor b increases, phase jumps at the contacts, currents in the central part of the vortex, and the total energies of the vortices decrease in both configurations. This leads to a decrease in critical field H _c1. For all values of I and b, the energy of the vortex centered on the contact is higher than that of the vortex centered in the middle of the cell.     

Interaction and pinning of plane vortices in a three-dimensional Josephson medium and possible distances between two isolated vortices

Within a continuous vortex model, exact expressions are obtained for the Josephson and magnetic energies of plane (laminar) vortices, as well as for the energy and force of pinning by cells in a three-dimensional Josephson medium. If the porosity of the medium is taken into account, the Josephson and magnetic energies of the vortex differ from those for the continuum case. The contributions to the pinning energy from the Josephson and magnetic energies have opposite signs. An algorithm for numerically solving a system of difference equations is proposed in order to find the shape and the energy of the vortex in its stable and unstable states. The continuous vortex model is shown to fail in predicting correct values of the Josephson and magnetic energy of the vortex, as well as of the pinning energy components. Expressions for the least possible distances between two isolated vortices are obtained for a small pinning parameter. Analytical results are in close agreement with computer simulation. An algorithm for numerically solving a system of difference equations is proposed in order to find the least possible distances between two isolated vortices when the pinning parameter I is not small. The minimal value of I at which the center-to-center distance N of the vortices equals three cells is 1.428; for N=2, I _min=1.947. At I>2.907, the vortices can be centered in adjacent cells.     

Pinning of vortices and penetration of the magnetic field into a periodically modulated long Josephson contact

The upper field of the Meissner regime, H _up, and overheat field H′_c1, above which vortices start penetrating into a Josephson contact, are calculated throughout the range of pinning parameter I. The stability of likely configurations is investigated. It is shown that H _up = H′_c1 at any I. The existence of a single vortex centered at the extreme cell in the contact is demonstrated to be a possibility. At I > 3.69, such a vortex may exist even in a zero magnetic field. At 1.48 < I < 3.69, this vortex can exist in an external field in the range from some H _v to H _up. At I < 1.48, the vortex cannot exist under any conditions. From the equality of H _up and H′_c1 at any I, the conclusion is drawn that penetration of vortices into any Josephson medium is conditioned by the need to satisfy flux quantization conditions. Here, not the forces of vortex pinning at defects in the medium but quantization requirements are of major importance, which are satisfied in specific quantum ways rather than by meeting equilibrium conditions for vortices, forces, etc.     

Penetration of the magnetic field into a 3D ordered Josephson medium at very small values of the pinning parameter

The Meissner state of a 3D Josephson medium is analyzed for stability against small fluctuations of phase discontinuities at contacts. For any form of fluctuations, there exists value I ₀ of pinning parameter I such that the Meissner configuration remains stable if I < I ₀. Reasons why the configuration remains stable at small I are considered. Instability arises when the quadratic form of the second variation of Gibbs potential G is not a positively definite quantity. At small I, the contribution of the Josephson energy to G is small. The second variation of the magnetic energy, the other component of G, is always a positively definite quadratic form. Therefore, instability may arise only if I has a finite value. This statement holds true not only for the Meissner but also for any equilibrium configuration. At I < I ₀, stability persists up to the boundary of the Meissner state. Then, a sequence of plane vortices parallel to the boundary appears throughout the sample. Thus, vortices appearing at I < I ₀ are plane vortices rather than linear. The configurations of currents and the magnetic field profile inside the sample are calculated for I < I ₀. Calculation is based on analyzing the continuous variation of the current configuration toward a decrease in the Gibbs potential.     

Equilibrium states of planar vortices in a three-dimensional josephson medium and meaning of the pinning-energy concept

Various ways of specifying the pinning-energy concept for planar vortices in a three-dimensional cellular Josephson medium are analyzed. It is shown that, for values of the pinning parameter I that are not small, a universal characteristic of vortex interaction with the lattice cannot be found, since the displacement of a vortex distorts its shape. At small values of I, the maximum pinning force can be chosen for such a characteristic. Two equilibrium states of a vortex are analyzed for stability. It is revealed that the state of higher energy is not inevitably unstable. A correct analysis of stability must be based on exploring a quadratic form that describes the energy of a current configuration. Such an investigation is performed for the equilibrium state of a vortex. At small values of the pinning parameter, the vortex state of higher energy is quasistable.     

Finite vortex numbers and symmetric vortex structures in a rotating trapped Fermi gas in the BCS-BEC crossover

The ground state of a three-dimensional (3D) rotating trapped superfluid Fermi gas in the BCS-BEC crossover is mapped to finite N _v -body vortex states by a simple ansatz. The total vortex energy is measured from the ground-state energy of the system in the absence of the vortices. The vortex state is stable since the vortex potential and rotation energies are attractive while the vortex kinetic energy and interaction between vortices are repulsive. By combining the analytical and numerical works for the minimal vortex energy, the 2D configurations of N _v vortices are studied by taking into account of the finite size effects both on xy-plane and on z-direction. The calculated vortex numbers as a function of the interaction strength are appropriate to the renew experimental results by Zwierlein in [High-temperature superfluidity in a ultracold Fermi gas, Ph.D. thesis, Massachusetts Institute of Technology, 2006]. The numerical results show that there exist two types of vortex structures: the trap center is occupied and unoccupied by a vortex, even in the case of N _v < 10 with regular polygon and in the case of N _v ≥ 10 with finite triangle lattice. The rotation frequency dependent vortex numbers with different interaction strengths are also discussed.     

Line vortices in a three-dimensional ordered Josephson medium at a small pinning parameter

The structure and energy of a line vortex whose axis is aligned with the symmetry axis of a finite-thickness slab indefinitely long in two directions is calculated by solving a set of linear finite-difference equations. Fluxoid quantization conditions in cells near the center of the vortex serve as boundary conditions. An exact solution is approached by iterations in phase stepwise discontinuities that cannot be considered small. A close similarity between the configuration under study and a periodic sequence (chain) of vortices makes it possible to allow for the effect of the domain boundary on the structure and energy of the vortex. It is shown that, at any width of the slab, one can find a pinning parameter value so small that the vortex cannot be viewed as solitary and contributions from other vortices should be taken into account in calculation. Proceeding in this way, one can find the structure and energy of the vortex however small the pinning parameter is. The total energy of the vortex is its intrinsic energy plus the sum of its energies of interaction with other members of the chain. In turn, the intrinsic energy is the sum of the energies of the small discrete core and quasi-continuous outer shell. It is demonstrated that the energy of the core is a linear function of the pinning parameter and is comparable to the energy of the shell.     

Topological self-dual configurations in a Maxwell–Higgs model with a CPT-odd and Lorentz-violating nonminimal coupling

We have studied the existence of topological self-dual configurations in a nonminimal CPT-odd and Lorentz-violating (LV) Maxwell–Higgs model, where the LV interaction is introduced by modifying the minimal covariant derivative. The Bogomol’nyi–Prasad–Sommerfield formalism has been implemented, revealing that the scalar self-interaction implying self-dual equations contains a derivative coupling. The CPT-odd self-dual equations describe electrically neutral configurations with finite total energy proportional to the total magnetic flux, which differ from the charged solutions of other CPT-odd and LV models previously studied. In particular, we have investigated the axially symmetrical self-dual vortex solutions altered by the LV parameter. For large distances, the profiles possess general behavior similar to the vortices of Abrikosov–Nielsen–Olesen. However, within the vortex core, the profiles of the magnetic field and energy can differ substantially from ones of the Maxwell–Higgs model depending if the LV parameter is negative or positive.     

Gauge cosmic chiral strings in general relativity

Cylindrically symmetric self-gravitating configurations of string (vortex) type are considered within the framework of the chiral SU(2) model with the inclusion of the Yang-Mills proper gauge field. In the approximation of the large topological charge n the solutions to the field equations are found, with the magnetic field of the vortex being longitudinal. The linear energy density of the vortex configuration is estimated. The text was submitted by the author in English.     

Penetration of magnetic field into a long periodically modulated Josephson contact

A new approach to magnetic field profiling inside a Josephson contact is suggested. Its essence consists in analyzing continuous variation of a current configuration leading to a decrease in the Gibbs potential. With this approach, one can find a configuration into which the Meissner state turns when an external field slightly exceeds the upper boundary of the Meissner regime and trace the evolution of this configuration with increasing field. Calculations show that there exists critical value I _c of the pinning parameter in the range 0.95–1.00. This critical value separates two possible conditions of magnetic field penetration into the contact. At I > I _c, a near-boundary current configuration completely compensating for the external field inside the contact arises irrespective of the external field strength. At I < I _c, such a situation is observed only until the external field strength exceeds certain value H _max. Higher fields penetrate into the contact indefinitely deep. In nearboundary configurations, the magnetic field drops with increasing depth almost linearly. Its slope k has rational values, which remain constant within finite intervals of I. As I goes beyond a given interval, k rises stepwise and takes on another rational value. When an external magnetic field is switched on adiabatically, configurations with a maximal growth rate of the magnetic field are observed.     

Vortex rings moving into turbulence

We report on an experimental study of turbulent vortex rings injected with velocity U _v0 into a grid-generated turbulent flow (with RMS streamwise velocity u _*) and followed relative to the mean flow. The initial Reynolds number of the vortices varies from 4500 to 11,500. The turbulence was characterised by an intensity I_t =u _*/U _v0, which varied over the range 0<I_t <0.03. A mathematical model based on a stochastic model of the vortex core is developed to explain and interpret the results. The vortex radius grows diffusively in time with the rate of increase of the square of the vortex radius increasing linearly with I_t . As the vortices grow, they slow down sufficiently rapidly in a manner that they penetrate a finite distance into the turbulence. The vortex velocity, averaged over many experiments, showed an initial t ^?1 decay, consistent with Maxworthy’s experiments. The analysis and experiments show that such vortices ultimately only move a finite distance from their point of generation and this distance varies inversely with I_t .     

Finite lifetime broadening of calculated X‐ray absorption spectra: possible artefacts close to the edge

X‐ray absorption spectra calculated within an effective one‐electron approach have to be broadened to account for the finite lifetime of the core hole. For methods based on Green's function this can be achieved either by adding a small imaginary part to the energy or by convoluting the spectra on the real axis with a Lorentzian. By analyzing the Fe K‐ and L_2,3‐edge spectra it is demonstrated that these procedures lead to identical results only for energies higher than a few core‐level widths above the absorption edge. For energies close to the edge, spurious spectral features may appear if too much weight is put on broadening via the imaginary energy component. Special care should be taken for dichroic spectra at edges which comprise several exchange‐split core levels, such as the L₃‐edge of 3d transition metals.     

Entropy of self-gravitating radiation

We examine the entropy of self-gravitating radiation confined to a spherical box of radiusR in the context of general relativity. We expect that configurations (i.e., initial data) which extremize total entropy will be spherically symmetric, time symmetric distributions of radiation in local thermodynamic equilibrium. Assuming this is the case, we prove that extrema ofS coincide precisely with static equilibrium configurations of the radiation fluid. Furthermore, dynamically stable equilibrium configurations are shown to coincide with local maxima ofS. The equilibrium configurations and their entropies are calculated and their properties are discussed. However, it is shown that entropies higher than these local extrema can be achieved and, indeed, arbitrarily high entropies can be attained by configurations inside of or outside but arbitrarily near their own Schwarzschild radius. However, if we limit consideration to configurations which are outside their own Schwarzschild radius by at least one radiation wavelength, then the entropy is bounded and we find S_max MR, whereM is the total mass. This supports the validity for self-gravitating systems of the Bekenstein upper limit on the entropy to energy ratio of material bodies.     

Penetration of the magnetic field into a 3D ordered Josephson medium at small values of the pinning parameter

The current configurations and the profile of the magnetic field penetrating into a 3D ordered Josephson medium are calculated for I < I _C . The calculation algorithm (modified for finite-length samples) is based on analyzing the continuous variation of the configuration toward a decrease in the Gibbs potential. This algorithm makes it possible to find a configuration into which the Meissner state passes when I < I _C and an external field slightly exceeds H _max and trace the evolution of this configuration with a further rise in the field. At H > H _max, the magnetic field penetrates into the sample as a quasi-uniform sequence of plane vortices. When H is roughly equal to H ₀/2, where H ₀ is the outer field at which one fluxoid Φ₀ passes through each cell, the plane vortices disintegrate into linear ones centered in cells neighboring along the diagonal. As the field grows, the vortex pattern condenses: zero-fluxoid cells are gradually “filled” starting from the boundary. When the field approaches H ₀, a sequence of plane vortices centered in adjacent rows arises near the boundary. With a further increase in the field, sequences of linear vortices with a double fluxoid form at the boundary. Then, such a scenario is periodically repeated with a period H ₀ in the external field.     

Finite number of vortices and bending of finite vortex lines in a confined rotating Bose-Einstein condensate

The minimal energy configurations of finite N_v-body vortices in a rotating trapped Bose-Einstein condensate is studied analytically by extending the previous work [Y. Castin, R. Dum, Eur. Phys. J. D 7, 399 (1999)], and taking into account the finite size effects on z-direction and the bending of finite vortex lines. The calculation of the energy of the vortices as a function of the rotation frequency of the trap gives number, curvature, configuration of vortices and width of vortex cores self-consistently. The numerical results show that (1) the simplest regular polynomial of the several vortex configurations is energetically favored; while the hexagonal vortex lattice is more stable than square lattice; (2) bending is more stable then straight vortex line along the z-axis for λ<1; (3) the boundary effect is obvious: compared with the estimation made under infinite boundary, the finite size effect leads to a lower vortex density, while the adding vortex bending results in a less higher density because of the expansion. The results are in well agreement with the other authors' ones.     

Different methods of calculating the pinning energy of plane vortices in a 3D Josephson medium: A comparative study

The pinning energy of plane (laminar) vortices in a 3D Josephson medium is calculated within a continuous vortex model considering functions of two types: V=1−cosϕ and V= 2/π⁴ϕ²(2π−ϕ)². The shape and energy of the stable and unstable vortices are found with an algorithm for the exact numerical solution of a set of difference equations. The vortex magnetic and Josephson energies diverge. The magnetic and Josephson components of the pinning energy are close in magnitude but differ in sign; as a result, the total pinning energy is smaller than its components by one order of magnitude. This result is confirmed analytically. An analytical computing method within the continuous vortex model is suggested. This method preserves the difference terms in the energy expression. The magnetic energy found by this method differs from the Josephson energy in magnitude, and the magnetic component of the pinning energy is opposite in sign to the Josephson component. Comparative analysis of the approximate approaches to energy calculation within the continuous vortex model when the difference terms are retained and when they are replaced by derivatives is performed. It is shown that the continuous vortex model gives incorrect values of the Josephson and magnetic components of the pinning energy. The actual values are several tens or several hundreds of times higher than those obtained with the continuous vortex model. Yet, since the Josephson and magnetic components of the pinning energy have different signs, the exact value of the total pinning energy and the approximate value obtained within the continuous vortex model differ insignificantly.     

Self-interstitial configuration in molybdenum studied by modified analytical embedded atom method

The stability of various atomic configurations containing a self-interstitial atom (SIA) in a model representing Mo has been investigated using the modified analytical embedded atom method (MAEAM). The lattice relaxations are treated with the molecular dynamics (MD) simulation at absolute zero of temperature. Six relatively stable self-interstitial configurations and formation energies have been described and calculated. The results indicate that the [111] dumbbell interstitial S₁₁₁ has the lowest formation energy, and in ascending order, the sequence of the configurations is predicted to be S₁₁₁, C, S₁₁₀, T, S₀₀₁ and O. From relaxed displacement field up to the fifth-NN atoms of six configurations, we know that the relaxed displacements depend not only on separation distances of the NN atoms from the defect centre but also strongly on the direction of the connected line between the NN atoms and the defect centre. The equilibrium distances between two nearest atoms in the core of the S₁₁₁, C, S₁₁₀, T, S₀₀₁ and O configurations are 0.72a, 0.72a, 0.71a, 0.72a, 0.70a and 0.70a, respectively.      

Estimate of the core size of a vortex in an extreme type II superconductor

The core size of a single straight vortex in an extreme type II superconductor is estimated at zero temperature. The free energy is expressed as a non-local functional of the order parameter and vector potential associated with the vortex. Profiles for the order parameter and vector potential that depend on the length scales ξ and λ, respectively, are assumed. For an order parameter profile that achieves its bulk value after the finite distance, ξ, the free energy approximately separates into a sum of two contributions, one depending on the order parameter and the other on the magnetic field. Using two different approximation schemes the two parts of the free energy are minimised. For an order parameter profile that, in its region of change, varies linearly with distance from the vortex axis, the length scale ξ is found to be comparable with the zero temperature coherence length, in agreement with the findings of Ullah et al. (Phys. Rev. B 42 (1990), 9950) for a one-parameter fit to the order parameter. The length scale associated with the magnetic field is also estimated.     

Vortex structures in planar lattices of magnetic dipoles in the presence of exchange coupling

Vortex equilibrium states of planar square lattices of magnetic dipoles in the presence of the exchange interaction have been studied. It has been shown that the vortex equilibrium configurations differ in the position of the vortex center and, correspondingly, in the magnitude and direction of the total magnetic moment of the system. In the case of the position of the vortex center in the center of the array, the total magnetic moment of the system is zero. The vortex center moves in the direction perpendicular to the field under the action of the external planar magnetic field on the system. Thus, the transitions between different equilibrium vortex configurations are implemented and the magnetic moment of the system of dipoles is controlled.     

Macroscopic potential-energy surfaces for symmetric fission and heavy-ion reactions

We calculate the macroscopic potential energy of deformation for symmetric configurations of interest in fission and heavy-ion reactions. The shape of the system is characterized in terms of two moments of the matter distribution. These moments correspond to the distance between the centers of mass of the two halves of the system and to the elongation of each half about its center of mass. The configurations studied include a continuous sequence of shapes from the sphere to two-, three-, and four-fragment scission lines. Beyond the scission lines and prior to the line of first contact in heavy-ion reactions we represent the system in terms of separated oblate and prolate spheroids. The macroscopic energy is calculated as the sum of a Coulomb energy and a nuclear macroscopic energy that takes into account the finite range of the nuclear force. For systems throughout the periodic table we display the calculated energy as a function of the two moments in the form of contour maps. Some important features of the contour maps are the binary, ternary, and quaternary saddle points, the fission and fusion (or two-fragment) valleys, and the three- and four-fragment valleys. The maps illustrate how the topography of the potential energy changes as a function of the nuclear system considered. For example, as we move from lighter to heavier nuclear systems the binary saddle point moves from outside the point of first contact in heavy-ion reactions to inside the contact point. Because of this, the formation of a heavy compound nucleus requires additional energy relative to the maximum in a one-dimensional interaction barrier. The maps also illustrate for moderately heavy systems the presence of separate valleys for binary fission and fusion. For still heavier systems the ternary and quaternary saddle points are no longer present. This means that the ternary and quaternary valleys are accessible by paths that decrease monotonically in energy beyond the binary saddle point. Finally, for nuclear systems heavier than about ³⁰⁰120, the binary saddle point itself disappears, which in the absence of single-particle effects precludes altogether the formation of a compound system.