---

makeiv.cc

Go to the documentation of this file.

// makeiv.cc
// Generate a simulated SIS DC IV curve suitable for use with supermix.
// The output IV curve will be normalized such that Vgap = Rn = 1.
// The output will be written to the standard output stream.

#include "supermix.h"
#include "numerical.h"

// ---------------------------------------------------------------------------
// Function to print a usage prompt, with definitions of the parameters

const int cmd_line_params = 4;    // the expected number of command line params

void prompt(const char *const *const argv)
{
  cerr << argv[0] << ": construct a simulated, normalized SIS IV curve for"
       << endl << "use with SuperMix." << endl << endl;
  cerr << "Usage: " << argv[0] << "  R_subgap  Gap_width  I_leakage  I_defect"
       << endl << endl;

  cerr << "Where:" << endl
       << "       Parameter    Description                         Valid Range \n"
       << "       ---------    ---------------------------------   ----------- \n"
       << "       R_subgap     Subgap/Normal resistance ratio      > 1.0       \n"
       << "       Gap_width    Width of the gap relative to Vgap   0.0 - 1.0   \n"
       << "       I_leakage    Leakage current just above V = 0    0.0 - 1.0   \n"
       << "       I_defect     Current defect just above gap       0.0 - 1.0   \n"
       << endl
       << "       The output is normalized so that Vgap = Rn = 1.0.            \n"
       << "       Currents (I_leakage and I_defect) are normalized to Vgap/Rn. \n"
       << "       I_defect is the amount that the current just above Vgap      \n"
       << "       falls short of being equal to V/Rn. It rerpesents a constant \n"
       << "       offset below the I=V/Rn line for all currents with V>Vgap.   \n"
       << endl
       << "       Typical values for a fairly good SIS might be:  \n"
       << "          R_subgap  = 20          \n"
       << "          Gap_width = 0.02        \n"
       << "          I_leakage = 0.01        \n"
       << "          I_defect  = 0.1         \n"
       << endl;
}

// ---------------------------------------------------------------------------
// Global parameters:

double  G_subgap;      // Subgap/Normal conductance ratio  (0.0 - 1.0)
double  Gap_width;     // Relative width of the gap        (0.0 - 1.0)
double  I_leakage;     // Leakage current/(Vgap/Rn)        (0.0 - 1.0)
double  I_defect;      // Current defect above gap         (0.0 - 1.0)

const char *comment = "#";

// ---------------------------------------------------------------------------
// Functions used to simulate a simple SIS dc IV response:

// output a straight line fit: y = mx+b:
inline double line(double x, double m, double b) { return m*x + b; }

// a hyperbolic tangent step function with some parameter checking
double sym_step(double x, double w)
{
  if (w > 0.0) return tanh(x/w);
  else {
    if (x  < 0.0) return -1.0;
    else if (x == 0.0) return  0.0;
    else return  1.0;
  }
}

// output a finite-width generalization of the Heaviside step function:
inline double step(double x, double w) { return 0.5*(sym_step(x,w)+1.0); }

// the full simulation of the iv curve, using the curve parameters:
double f(double v)
{
  double av     = fabs(v);
  double i_sg   = line(av, G_subgap, I_leakage * sym_step(v, sqrt(Gap_width)));
  double i_ag   = line(av, 1.0, -I_defect) * sym_step(v, Gap_width);
  double weight = step(av - 1.0, Gap_width);
  return weight*i_ag + (1.0-weight)*i_sg;
}

// ---------------------------------------------------------------------------
int main(int argc, char ** argv)
{
  // Input command line parameters; print a usage prompt if needed
  bool bad_param = false;
  if (argc != 1+cmd_line_params) bad_param = true;
  else {
    Gap_width = atof(argv[2]);
    G_subgap  = 1/atof(argv[1]);
    I_leakage = atof(argv[3]);
    I_defect  = atof(argv[4]);

    if (Gap_width <= 0.0 || Gap_width >  1.0) bad_param = true;
    if (G_subgap  <  0.0 || G_subgap  >  1.0) bad_param = true;
    if (I_leakage <  0.0 || I_leakage >  1.0) bad_param = true;
    if (I_defect  <  0.0 || I_defect  >  1.0) bad_param = true;
  }

  if (bad_param) {
    prompt(argv);
    return 1;
  }


  // Adaptively fill an interpolator with simulated IV curve points. We use
  // an interpolator here, since the SuperMix class ivcurve does. We want to
  // ensure that the points we output will allow an interpolator to accurately
  // match the desired IV curve response. See interpolate.h and adaptive.h

  interpolator<double> idc;
  adaptive_fill(idc, f, 0, 2, .0001, .01);


  // Now output the (V,I) points in the interpolator.

  cout << fixed << setprecision(10);
  cout << comment << " Normalized simulated SIS DC IV characteristic curve" << endl;
  cout << comment << " Parameter values used in the simulation:" << endl;
  cout << comment << "  Gap Width                = " << Gap_width << endl;
  cout << comment << "  Subgap Resistance Ratio  = " << 1/G_subgap  << endl;
  cout << comment << "  Subgap Current Leakage   = " << I_leakage << endl;
  cout << comment << "  Normal Current Defect    = " << I_defect  << endl;
  cout << comment << endl;
  cout << comment << " V:         " << "\t" << "I:" << endl;

  for(unsigned i = 0; i < idc.size(); ++i)
    cout << idc.x(i) << "\t" << idc[i] << endl;
}

---

Please direct comments and corrections to supermix@submm.caltech.edu
Go to the supermix home page
Generated by 1.2.7