Crossover from hc/e to hc/2e current oscillations in rings of s-wave superconductors

We analyze the crossover from an hc/e-periodicity of the persistent current in flux threaded clean metallic rings towards an hc/2e-flux periodicity of the supercurrent upon entering the superconducting state. On the basis of a model calculation for a one-dimensional ring we identify the underlying mechanism, which balances the hc/e versus the hc/2e periodic components of the current density. When the ring circumference exceeds the coherence length of the superconductor, the flux dependence is strictly hc/2e periodic. Further, we develop a multi-channel model which reduces the Bogoliubov - de Gennes equations to a one-dimensional differential equation for the radial component of the wave function. The discretization of this differential equation introduces transverse channels, whose number scales with the thickness of the ring. The periodicity crossover is analyzed close the critical temperature

Electron-beam-induced shift in the apparent position of a pinned vortex in a thin superconducting film

When an electron beam strikes a superconducting thin film near a pinned vortex, it locally increases the temperature-dependent London penetration depth and perturbs the circulating supercurrent, thereby distorting the vortex's magnetic field toward the heated spot. This phenomenon has been used to visualize vortices pinned in SQUIDs using low-temperature scanning electron microscopy. In this paper I develop a quantitative theory to calculate the displacement of the vortex-generated magnetic-flux distribution as a function of the distance of the beam spot from the vortex core. The results are calculated using four different models for the spatial distribution of the thermal power deposited by the electron beam.Comment: 9 pages, 6 figures, resubmitted to PRB with referee-suggested revisions, includes new paragraph on numerical evaluatio

Growth of a vortex polycrystal in type II superconductors

We discuss the formation of a vortex polycrystal in type II superconductors from the competition between pinning and elastic forces. We compute the elastic energy of a deformed grain boundary, that is strongly non-local, and obtain the depinning stress for weak and strong pinning. Our estimates for the grain size dependence on the magnetic field strength are in good agreement with previous experiments on NbMo. Finally, we discuss the effect of thermal noise on grain growth.Comment: 4 pages, 2 figure

Local Simulation Algorithms for Coulombic Interactions

We consider dynamically constrained Monte-Carlo dynamics and show that this leads to the generation of long ranged effective interactions. This allows us to construct a local algorithm for the simulation of charged systems without ever having to evaluate pair potentials or solve the Poisson equation. We discuss a simple implementation of a charged lattice gas as well as more elaborate off-lattice versions of the algorithm. There are analogies between our formulation of electrostatics and the bosonic Hubbard model in the phase approximation. Cluster methods developed for this model further improve the efficiency of the electrostatics algorithm.Comment: Proceedings Statphys22 10 page

Local Simulation Algorithms for Coulomb Interaction

Long ranged electrostatic interactions are time consuming to calculate in molecular dynamics and Monte-Carlo simulations. We introduce an algorithmic framework for simulating charged particles which modifies the dynamics so as to allow equilibration using a local Hamiltonian. The method introduces an auxiliary field with constrained dynamics so that the equilibrium distribution is determined by the Coulomb interaction. We demonstrate the efficiency of the method by simulating a simple, charged lattice gas.Comment: Last figure changed to improve demonstration of numerical efficienc

Fourier Acceleration of Langevin Molecular Dynamics

Fourier acceleration has been successfully applied to the simulation of lattice field theories for more than a decade. In this paper, we extend the method to the dynamics of discrete particles moving in continuum. Although our method is based on a mapping of the particles' dynamics to a regular grid so that discrete Fourier transforms may be taken, it should be emphasized that the introduction of the grid is a purely algorithmic device and that no smoothing, coarse-graining or mean-field approximations are made. The method thus can be applied to the equations of motion of molecular dynamics (MD), or its Langevin or Brownian variants. For example, in Langevin MD simulations our acceleration technique permits a straightforward spectral decomposition of forces so that the long-wavelength modes are integrated with a longer time step, thereby reducing the time required to reach equilibrium or to decorrelate the system in equilibrium. Speedup factors of up to 30 are observed relative to pure (unaccelerated) Langevin MD. As with acceleration of critical lattice models, even further gains relative to the unaccelerated method are expected for larger systems. Preliminary results for Fourier-accelerated molecular dynamics are presented in order to illustrate the basic concepts. Possible extensions of the method and further lines of research are discussed.Comment: 11 pages, two illustrations included using graphic

Maximum Flux Transition Paths of Conformational Change

Given two metastable states A and B of a biomolecular system, the problem is to calculate the likely paths of the transition from A to B. Such a calculation is more informative and more manageable if done for a reduced set of collective variables chosen so that paths cluster in collective variable space. The computational task becomes that of computing the "center" of such a cluster. A good way to define the center employs the concept of a committor, whose value at a point in collective variable space is the probability that a trajectory at that point will reach B before A. The committor "foliates" the transition region into a set of isocommittors. The maximum flux transition path is defined as a path that crosses each isocommittor at a point which (locally) has the highest crossing rate of distinct reactive trajectories. (This path is different from that of the MaxFlux method of Huo and Straub.) It is argued that such a path is nearer to an ideal path than others that have been proposed with the possible exception of the finite-temperature string method path. To make the calculation tractable, three approximations are introduced, yielding a path that is the solution of a nonsingular two-point boundary-value problem. For such a problem, one can construct a simple and robust algorithm. One such algorithm and its performance is discussed.Comment: 7 figure

Elasticity-driven interaction between vortices in type-II superconductors

The contribution to the vortex lattice energy which is due to the vortex-induced strains is calculated covering all the magnetic field range which defines the vortex state. This contribution is compared with previously reported ones what shows that, in the most part of the vortex state, it has been notably underestimated until now. The reason of such underestimation is the assumption that only the vortex cores induce strains. In contrast to what is generally assumed, both core and non-core regions are important sources of strains in high-$\kappa$ superconductors.Comment: 10 pages, 1 figure, revtex

Quantitative nanoscale vortex-imaging using a cryogenic quantum magnetometer

Microscopic studies of superconductors and their vortices play a pivotal role in our understanding of the mechanisms underlying superconductivity. Local measurements of penetration depths or magnetic stray-fields enable access to fundamental aspects of superconductors such as nanoscale variations of superfluid densities or the symmetry of their order parameter. However, experimental tools, which offer quantitative, nanoscale magnetometry and operate over the large range of temperature and magnetic fields relevant to address many outstanding questions in superconductivity, are still missing. Here, we demonstrate quantitative, nanoscale magnetic imaging of Pearl vortices in the cuprate superconductor YBCO, using a scanning quantum sensor in form of a single Nitrogen-Vacancy (NV) electronic spin in diamond. The sensor-to-sample distance of ~10nm we achieve allows us to observe striking deviations from the prevalent monopole approximation in our vortex stray-field images, while we find excellent quantitative agreement with Pearl's analytic model. Our experiments yield a non-invasive and unambiguous determination of the system's local London penetration depth, and are readily extended to higher temperatures and magnetic fields. These results demonstrate the potential of quantitative quantum sensors in benchmarking microscopic models of complex electronic systems and open the door for further exploration of strongly correlated electron physics using scanning NV magnetometry.Comment: Main text (5 pages, 4 figures) plus supplementary material (5 pages, 6 figures). Comments welcome. Further information under http://www.quantum-sensing.c

The regulatory subunit of PKA-I remains partially structured and undergoes β-aggregation upon thermal denaturation

Background: The regulatory subunit (R) of cAMP-dependent protein kinase (PKA) is a modular flexible protein that responds with large conformational changes to the binding of the effector cAMP. Considering its highly dynamic nature, the protein is rather stable. We studied the thermal denaturation of full-length RIα and a truncated RIα(92-381) that contains the tandem cyclic nucleotide binding (CNB) domains A and B. Methodology/Principal Findings: As revealed by circular dichroism (CD) and differential scanning calorimetry, both RIα proteins contain significant residual structure in the heat-denatured state. As evidenced by CD, the predominantly α-helical spectrum at 25°C with double negative peaks at 209 and 222 nm changes to a spectrum with a single negative peak at 212-216 nm, characteristic of β-structure. A similar α→β transition occurs at higher temperature in the presence of cAMP. Thioflavin T fluorescence and atomic force microscopy studies support the notion that the structural transition is associated with cross-β-intermolecular aggregation and formation of non-fibrillar oligomers. Conclusions/Significance: Thermal denaturation of RIα leads to partial loss of native packing with exposure of aggregation-prone motifs, such as the B' helices in the phosphate-binding cassettes of both CNB domains. The topology of the β-sandwiches in these domains favors inter-molecular β-aggregation, which is suppressed in the ligand-bound states of RIα under physiological conditions. Moreover, our results reveal that the CNB domains persist as structural cores through heat-denaturation. © 2011 Dao et al