A comparison of past and present computational methods for shape analysis of double-shell x-ray radiographs.

Implosion symmetry is a key requirement in achieving a robust burning plasma in inertial confinement fusion experiments. In double-shell capsule implosions, we are interested in the shape of the inner shell as it pushes on the fuel. Shape analysis is a popular technique for studying said symmetry during implosion. Combinations of filtering and contour-finding algorithms are studied for their promise in reliably recovering Legendre shape coefficients from synthetic radiographs of double-shell capsules with applied levels of noise. A radial lineout max(slope) method when used on an image pre-filtered with non-local means and a variant of the marching squares algorithm are able to recover p0, p2, and p4 maxslope Legendre shape coefficients with mean pixel discrepancy errors of 2.81 and 3.06, respectively, for the noisy synthetic radiographs we consider. This improves upon prior radial lineout methods paired with Gaussian filtering, which we show to be unreliable and whose performance is dependent on input parameters that are difficult to estimate.

Dynamics of the exponential integrate-and-fire model with slow currents and adaptation.

In order to properly capture spike-frequency adaptation with a simplified point-neuron model, we study approximations of Hodgkin-Huxley (HH) models including slow currents by exponential integrate-and-fire (EIF) models that incorporate the same types of currents. We optimize the parameters of the EIF models under the external drive consisting of AMPA-type conductance pulses using the current-voltage curves and the van Rossum metric to best capture the subthreshold membrane potential, firing rate, and jump size of the slow current at the neuron's spike times. Our numerical simulations demonstrate that, in addition to these quantities, the approximate EIF-type models faithfully reproduce bifurcation properties of the HH neurons with slow currents, which include spike-frequency adaptation, phase-response curves, critical exponents at the transition between a finite and infinite number of spikes with increasing constant external drive, and bifurcation diagrams of interspike intervals in time-periodically forced models. Dynamics of networks of HH neurons with slow currents can also be approximated by corresponding EIF-type networks, with the approximation being at least statistically accurate over a broad range of Poisson rates of the external drive. For the form of external drive resembling realistic, AMPA-like synaptic conductance response to incoming action potentials, the EIF model affords great savings of computation time as compared with the corresponding HH-type model. Our work shows that the EIF model with additional slow currents is well suited for use in large-scale, point-neuron models in which spike-frequency adaptation is important.