A quasi-linear control theory analysis of timesharing skills

The compliance of the human ankle joint is measured by applying 0 to 50 Hz band-limited gaussian random torques to the foot of a seated human subject. These torques rotate the foot in a plantar-dorsal direction about a horizontal axis at a medial moleolus of the ankle. The applied torques and the resulting angular rotation of the foot are measured, digitized and recorded for off-line processing. Using such a best-fit, second-order model, the effective moment of inertia of the ankle joint, the angular viscosity and the stiffness are calculated. The ankle joint stiffness is shown to be a linear function of the level of tonic muscle contraction, increasing at a rate of 20 to 40 Nm/rad/Kg.m. of active torque. In terms of the muscle physiology, the more muscle fibers that are active, the greater the muscle stiffness. Joint viscosity also increases with activation. Joint stiffness is also a linear function of the joint angle, increasing at a rate of about 0.7 to 1.1 Nm/rad/deg from plantar flexion to dorsiflexion rotation

Further observations on the relationship of EMG and muscle force

Human skeletal muscle may be regarded as an electro-mechanical transducer. Its physiological input is a neural signal originating at the alpha motoneurons in the spinal cord and its output is force and muscle contraction, these both being dependent on the external load. Some experimental data taken during voluntary efforts around the ankle joint and by direct electrical stimulation of the nerve are described. Some of these experiments are simulated by an analog model, the input of which is recorded physiological soleus muscle EMG. The output is simulated foot torque. Limitations of a linear model and effect of some nonlinearities are discussed

Computational problems in autoregressive moving average (ARMA) models

The choice of the sampling interval and the selection of the order of the model in time series analysis are considered. Band limited (up to 15 Hz) random torque perturbations are applied to the human ankle joint. The applied torque input, the angular rotation output, and the electromyographic activity using surface electrodes from the extensor and flexor muscles of the ankle joint are recorded. Autoregressive moving average models are developed. A parameter constraining technique is applied to develop more reliable models. The asymptotic behavior of the system must be taken into account during parameter optimization to develop predictive models

A stochastic model of the electromyogram

The quantitative regularities of interference pattern formation by motor unit action potentials is investigated. The parameters of a single motor unit and how they relate to the Fourier transform analysis of an EMG are considered. The Fourier transform of the simulated electromyogram is compared with the Fourier transform of the actual EMG recorded from various human muscles using surface electrodes

Radiative Transitions in Charmonium from Lattice QCD

Radiative transitions between charmonium states offer an insight into the internal structure of heavy-quark bound states within QCD. We compute, for the first time within lattice QCD, the transition form-factors of various multipolarities between the lightest few charmonium states. In addition, we compute the experimentally unobservable, but physically interesting vector form-factors of the $\eta_c, J/\psi$ and $\chi_{c0}$. To this end we apply an ambitious combination of lattice techniques, computing three-point functions with heavy domain wall fermions on an anisotropic lattice within the quenched approximation. With an anisotropy $\xi=3$ at $a_s \sim 0.1 \mathrm{fm}$ we find a reasonable gross spectrum and a hyperfine splitting $\sim 90 \mathrm{MeV}$, which compares favourably with other improved actions. In general, after extrapolation of lattice data at non-zero $Q^2$ to the photopoint, our results agree within errors with all well measured experimental values. Furthermore, results are compared with the expectations of simple quark models where we find that many features are in agreement; beyond this we propose the possibility of constraining such models using our extracted values of physically unobservable quantities such as the $J/\psi$ quadrupole moment. We conclude that our methods are successful and propose to apply them to the problem of radiative transitions involving hybrid mesons, with the eventual goal of predicting hybrid meson photoproduction rates at the GlueX experiment.Comment: modified version as publishe

Two-Flavor Staggered Fermion Thermodynamics at N_t = 12

We present results of an ongoing study of the nature of the high temperature crossover in QCD with two light fermion flavors. These results are obtained with the conventional staggered fermion action at the smallest lattice spacing to date---approximately 0.1 fm. Of particular interest are a study of the temperature of the crossover a determination of the induced baryon charge and baryon susceptibility, the scalar susceptibility, and the chiral order parameter, used to test models of critical behavior associated with chiral symmetry restoration. From our new data and published results for N_t = 4, 6, and 8, we determine the QCD magnetic equation of state from the chiral order parameter using O(4) and mean field critical exponents and compare it with the corresponding equation of state obtained from an O(4) spin model and mean field theory. We also present a scaling analysis of the Polyakov loop, suggesting a temperature dependent ``constituent quark free energy.''Comment: LaTeX 25 pages, 15 Postscript figure

Quarkonium mass splittings in three-flavor lattice QCD

We report on calculations of the charmonium and bottomonium spectrum in lattice QCD. We use ensembles of gauge fields with three flavors of sea quarks, simulated with the asqtad improved action for staggered fermions. For the heavy quarks we employ the Fermilab interpretation of the clover action for Wilson fermions. These calculations provide a test of lattice QCD, including the theory of discretization errors for heavy quarks. We provide, therefore, a careful discussion of the results in light of the heavy-quark effective Lagrangian. By and large, we find that the computed results are in agreement with experiment, once parametric and discretization errors are taken into account.Comment: 21 pages, 17 figure

Probing the Region of Massless Quarks in Quenched Lattice QCD using Wilson Fermions

We study the spectrum of $H(m)=\gamma_5 W(-m)$ with $W(m)$ being the Wilson-Dirac operator on the lattice with bare mass equal to $m$. The background gauge fields are generated using the SU(3) Wilson action at $\beta=5.7$ on an $8^3\times 16$ lattice. We find evidence that the spectrum of $H(m)$ is gapless for $1.02 < m < 2.0$, implying that the physical quark is massless in this whole region.Comment: 22 pages, LaTeX file, uses elsart.sty, includes 11 figures A typographical error in one reference has been fixe

Color confinement and dual superconductivity in full QCD

We report on evidence that confinement is related to dual superconductivity of the vacuum in full QCD, as in quenched QCD. The vacuum is a dual superconductor in the confining phase, whilst the U(1) magnetic symmetry is realized a la Wigner in the deconfined phase.Comment: 4 pages, 4 eps figure

Determination of Inter-Phase Line Tension in Langmuir Films

A Langmuir film is a molecularly thin film on the surface of a fluid; we study the evolution of a Langmuir film with two co-existing fluid phases driven by an inter-phase line tension and damped by the viscous drag of the underlying subfluid. Experimentally, we study an 8CB Langmuir film via digitally-imaged Brewster Angle Microscopy (BAM) in a four-roll mill setup which applies a transient strain and images the response. When a compact domain is stretched by the imposed strain, it first assumes a bola shape with two tear-drop shaped reservoirs connected by a thin tether which then slowly relaxes to a circular domain which minimizes the interfacial energy of the system. We process the digital images of the experiment to extract the domain shapes. We then use one of these shapes as an initial condition for the numerical solution of a boundary-integral model of the underlying hydrodynamics and compare the subsequent images of the experiment to the numerical simulation. The numerical evolutions first verify that our hydrodynamical model can reproduce the observed dynamics. They also allow us to deduce the magnitude of the line tension in the system, often to within 1%. We find line tensions in the range of 200-600 pN; we hypothesize that this variation is due to differences in the layer depths of the 8CB fluid phases.Comment: See (http://www.math.hmc.edu/~ajb/bola/) for related movie