Source author record

C. G. Gray

C. G. Gray appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR Biological Physics cond-mat.stat-mech Neurons and Cognition physics.chem-ph physics.comp-ph

Catalog footprint

What is connected

4works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

Improved Hodgkin & Huxley-type model for action potentials in squid

By extending the crude Goldman-Hodgkin-Katz electrodiffusion model for resting-state membrane potentials in perfused giant axons of squid, we reformulate the Hodgkin-Huxley (HH) phenomenological quantitative model to create a new model which is simpler and based more fundamentally on electrodiffusion principles. Our dynamical system, like that of HH, behaves as a 4-dimensional resonator exhibiting subthreshold oscillations. The predicted speed of propagating action potentials at 20 degrees Celsius is in good agreement with the HH experimental value at 18.5 degrees Celsius. After the external concentration of calcium ions is reduced, the generation of repetitive rebound action potentials is predicted by our model, in agreement with experiment, when the membrane is stimulated by a brief (0.1 ms) depolarizing current. Unlike the HH model, our model predicts, in agreement with experiment, that prolonged constant-current stimulation does not generate spike trains in perfused axons. Our resonator model predicts rebound spiking following prolonged hyperpolarizing stimulation, observed at 18.5 degrees Celsius by HH but not predicted at this temperature by their quantitative model. Spiking promoted by brief hyperpolarization is also predicted, at room temperature, by our electrodiffusion model, but only at much lower temperatures (ca. 6 degrees Celsius) by the HH model. We discuss qualitatively, more completely than do HH, temperature dependences of the various physical effects which determine resting and action potentials.

preprint2016arXiv

Stiff-spring approximation revisited: inertial effects in non-equilibrium trajectories

Use of harmonic guiding potentials is the most common method for implementing steered molecular dynamics (SMD) simulations, performed to obtain potentials of mean force (PMFs) of molecular systems using non-equilibrium work (NEW) theorems. Harmonic guiding potentials are also the natural choice in single molecule force spectroscopy experiments. The stiff spring approximation (SSA) of Schulten and coworkers enables to use the work performed along SMD trajectories to obtain the PMF. We discuss and demonstrate how a high spring constant, k, required for the validity of the SSA can violate another requirement of this theory, i.e., the validity of Brownian dynamics of the system. Violation of the Brownian condition results in the introduction of kinetic energy contributions to the external work, performed during SMD simulations. These inertial effects result in skewed work distributions, rather than the Gaussian distributions predicted by SSA. The inertial effects also result in broader work distributions, which worsen the effect of the skewness when calculating work averages. Remarkably, our results strongly suggest that the skew and width of work distributions are independent of the average drift velocity and physical asymmetries. The skew and broadening of work distributions result in biased estimation of the PMF. The bias manifests itself in the form of a systematic error that increases with simulation time. We discuss the proper upper limit for k, such that the inertial effects are avoided. This limit, used together with the relation for the lower limit of k, enables to conduct accurate steering while satisfying the Brownian dynamics. Furthermore, we argue and demonstrate that using the peak-value (rather than the statistical mean) of the work distributions vastly reduces the bias in the calculated PMFs and improves the accuracy.

preprint2014arXiv

McMillan-Mayer Theory of Solutions Revisited: Simplifications and Extensions

McMillan and Mayer (MM) proved two remarkable theorems in their paper on the equilibrium statistical mechanics of liquid solutions. They first showed that the grand canonical partition function for a solution can be reduced to a one with an effectively solute-only form, by integrating out the solvent degrees of freedom. The total effective solute potential in the effective solute grand partition function can be decomposed into components which are potentials of mean force for isolated groups of one, two, three, etc, solute molecules. Secondly, from the first result, now assuming low solute concentration, MM derived an expansion for the osmotic pressure in powers of the solute concentration, in complete analogy with the virial expansion of gas pressure in powers of the density at low density. The molecular expressions found for the osmotic virial coefficients have exactly the same form as the corresponding gas virial coefficients, with potentials of mean force replacing vacuum potentials. In this paper we restrict ourselves to binary liquid solutions with solute species $A$ and solvent species $B$ and do three things: (a) By working with a semi-grand canonical ensemble (grand with respect to solvent only) instead of the grand canonical ensemble used by MM, and avoiding graphical methods, we have greatly simplified the derivation of the first MM result,(b) by using a simple nongraphical method developed by van Kampen for gases, we have greatly simplified the derivation of the second MM result,i.e.,the osmotic pressure virial expansion; as a by-product, we show the precise relation between MM theory and Widom potential distribution theory, and (c) we have extended MM theory by deriving virial expansions for other solution properties such as the enthalpy of mixing. The latter expansion, with changed independent variables corresponding to current experiments, is proving useful.

preprint2010arXiv

The magnetic fields of forming solar-like stars

Magnetic fields play a crucial role at all stages of the formation of low mass stars and planetary systems. In the final stages, in particular, they control the kinematics of in-falling gas from circumstellar discs, and the launching and collimation of spectacular outflows. The magnetic coupling with the disc is thought to influence the rotational evolution of the star, while magnetised stellar winds control the braking of more evolved stars and may influence the migration of planets. Magnetic reconnection events trigger energetic flares which irradiate circumstellar discs with high energy particles that influence the disc chemistry and set the initial conditions for planet formation. However, it is only in the past few years that the current generation of optical spectropolarimeters have allowed the magnetic fields of forming solar-like stars to be probed in unprecedented detail. In order to do justice to the recent extensive observational programs new theoretical models are being developed that incorporate magnetic fields with an observed degree of complexity. In this review we draw together disparate results from the classical electromagnetism, molecular physics/chemistry, and the geophysics literature, and demonstrate how they can be adapted to construct models of the large scale magnetospheres of stars and planets. We conclude by examining how the incorporation of multipolar magnetic fields into new theoretical models will drive future progress in the field through the elucidation of several observational conundrums.