Circuit Simulator including Josephson Junctions in Python

# PyJJASim

Circuit Simulator including Josephson Junctions in Python

# Installation

PyJJASim requires numpy, scipy and matplotlib


pip install pyjjasim



# Introduction

PyJJASim is a circuit simulator including Josephson Junctions as components,

intended to be used on large Josephson Junction Arrays (JJAs).

PyJJASim is specialized in keeping track of Josephson vortices in the circuit.

It can also compute static configurations that have vortices at desired

locations in the circuit.

This requires that the circuit is a planar embedding (in 2D), such that

one can unambiguously refer to faces of the circuit, and vortices reside

at faces. This imposes that nodes in the circuit must be placed

at 2D coordinates, and that no junctions can cross. This also means no

hierarchical structure is supported.

# Features

• supports basic components (inductors, resistors, capacitors and current- and voltage sources)

• keep track of (and place) Josephson vortices in the circuit

• compute static configurations

• determine dynamic stability of static configurations

• maximize parameters that have stable static configurations

• compute time evolutions

• define external magnetic flux through each face

• thermal fluctuations modeling nonzero temperature

• visualization and animation of simulation results

# Documentation

API \

whitepaper (in progress)

# Example Usage

from pyjjasim import *

sq_array = SquareArray(10, 10)

problem = StaticProblem(sq_array, frustration=0.1)

config, status, info = problem.compute()

config.plot(node_quantity="phi")


Program output:

n = np.zeros(sq_array.face_count())

n[sq_array.locate_faces(x=[2.5,6.5], y=[2.5,6.5])] = 1

config, status, info = problem.new_problem(vortex_configuration=n).compute()

config.plot(node_quantity="phi")

plt.show()


Program output:

More examples:

figure 1: Example of a frustrated honeycomb array with horizontal current bias.

Snapshot of a time evolution; the vortices are drifting downward. The electric potential

gradually increases in the direction of the external current

(see examples/images/biassed_honeycomb.py).

figure 2: Example of giant Shapiro steps in a square array. The external

current has a DC and AC component causing resonance in the form of

voltage plateaus (see examples/images/shapiro_steps.py).

figure 3: Example of a disordered array with different types of Josephson vortices

(see examples/images/disordered.py).

If you have any questions, comments, complaints,

bug reports, feature requests, etc.

## Release history Release notifications | RSS feed

Uploaded source
Uploaded py3