Quantum calcium-ion interactions with EEG

Previous papers have developed a statistical mechanics of neocortical interactions (SMNI) fit to short-term memory and EEG data. Adaptive Simulated Annealing (ASA) has been developed to perform fits to such nonlinear stochastic systems. An N-dimensional path-integral algorithm for quantum systems, qPATHINT, has been developed from classical PATHINT. Both fold short-time propagators (distributions or wave functions) over long times. Previous papers applied qPATHINT to two systems, in neocortical interactions and financial options. \textbf{Objective}: In this paper the quantum path-integral for Calcium ions is used to derive a closed-form analytic solution at arbitrary time that is used to calculate interactions with classical-physics SMNI interactions among scales. Using fits of this SMNI model to EEG data, including these effects, will help determine if this is a reasonable approach. \textbf{Method}: Methods of mathematical-physics for optimization and for path integrals in classical and quantum spaces are used for this project. Studies using supercomputer resources tested various dimensions for their scaling limits. In this paper the quantum path-integral is used to derive a closed-form analytic solution at arbitrary time that is used to calculate interactions with classical-physics SMNI interactions among scales. \textbf{Results}: The mathematical-physics and computer parts of the study are successful, in that there is modest improvement of cost/objective functions used to fit EEG data using these models. \textbf{Conclusion}: This project points to directions for more detailed calculations using more EEG data and qPATHINT at each time slice to propagate quantum calcium waves, synchronized with PATHINT propagation of classical SMNI.


Introduction
This project calculates quantum Ca 2+ interactions with EEG. In this paper, EEG is synonymous with large-scale neocortical firings during attentional tasks as measured by large-amplitude electroencephalographic (EEG) recordings. In this paper, only very specific calcium ions, Ca 2+ , are considered, those arising from regenerative calcium waves generated at tripartite neuron-astrocyteneuron synapses. Indeed, it is important to note that Ca 2+ ions, and specifically Ca 2+ waves, influence many processes in the brain, but this study focuses on free waves generated at tripartite synapses because of their calculated direct interactions with large synchronous neuronal firings.
Section 2 reviews the background of the main model used, Statistical Mechanics of Neocortical Interactions (SMNI).
Section 3 reviews the code Adaptive Simulated Annealing (ASA), used for optimization of many systems -fitting models to real data, e.g., fits to EEG data reported here.
Section 4 reviews the development of path-integral codes, PATHINT and qPATHINT, used for propagation of conditional probabilities and quantum-mechanical wave-functions, as reported here.
Section 5 gives new results of inclusion of quantum-mechanical interactions of Ca 2+ wave-packets with EEG.
Section 6 reviews some applications of this project. Section 7 gives the conclusion. The theory and codes for ASA and [q]PATHINT have been well tested across many disciplines by multiple users. This particular project most certainly is speculative, but it is testable. As reported here, fitting such models to EEG tests some aspects of this project. This is a somewhat indirect path, but not novel to many physics paradigms that are tested by experiment or computation. A detailed future path is described in the [q]PATHINT review Section.
While SMNI has been developed since 1981, and been confirmed by many tests, this evolving model including ionic scales has been part of multiple papers relatively recently, since 2012. Classical physics calculations support these extended SMNI models and are consistent with experimental data. Quantum physics calculations also support these extended SMNI models and, while they too are consistent with experimental data, it is quite speculative that they can persist in neocortex. Admittedly, it is surprising that detailed calculations continue to support this model, and so it is worth continued examination it until it is theoretically or experimentally proven to be false.

Synaptic Interactions
The short-time conditional probability distribution of firing of a given neuron firing given justprevious firings of other neurons is calculated from chemical and electrical intra-neuronal interactions (Ingber, , 1983. Given its previous interactions with k neurons within τ j of 5-10 msec, the conditional probability that neuron j fires (σ j = +1) or does not fire (σ j = −1) is The contribution to polarization achieved at an axon given activity at a synapse, taking into account averaging over different neurons, geometries, etc., is given by Γ, the "intra-neuronal" probability distribution. Ψ is the "inter-neuronal" probability distribution, of thousands of quanta of neurotransmitters released at one neuron's presynaptic site effecting a (hyper-)polarization at another neuron's postsynaptic site, taking into account interactions with neuromodulators, etc. This development holds for Γ Poisson, and for Ψ Poisson or Gaussian. V j is the depolarization threshold in the somatic-axonal region. v jk is the induced synaptic polarization of E or I type at the axon, and φ jk is its variance. The efficacy a jk is a sum of A jk from the connectivity between neurons, activated if the impinging k-neuron fires, and B jk from spontaneous background noise. The efficacy is related to the impedance across synaptic gaps.

Neuronal Interactions
Aggregation up to the mesoscopic scale from the microscopic synaptic scale uses mesoscopic probability P M represents a mesoscopic scale of columns of N neurons, with subsets E and I, represented by p q i . The "delta"-functions δ-constraint represents an aggregate of many neurons in a column. G is used to represent excitatory (E) and inhibitory (I) contributions.Ḡ designates contributions from both E and I. The path integral is derived in terms of mesoscopic Lagrangian L. The short-time distribution of firings in a minicolumn, given its just previous interactions with all other neurons in its macrocolumn, is thereby defined.

Columnar Interactions
In the prepoint (Ito) representation the SMNI Lagrangian L is The threshold factor F G is derived as where A G G is the columnar-averaged direct synaptic efficacy, B G G is the columnar-averaged backgroundnoise contribution to synaptic efficacy. The " ‡ " parameters arise from regional interactions across many macrocolumns.

SMNI Parameters From Experiments
All values of parameters and their bounds are taken from experimental data, not arbitrarily fit to specific phenomena.
N G = {N E = 160, N I = 60} was set for for visual neocortex, {N E = 80, N I = 30} was set for all other neocortical regions, M G and N G in F G are afferent macrocolumnar firings scaled to efferent minicolumnar firings by N/N * ≈ 10 −3 . N * is the number of neurons in a macrocolumn, about 10 5 . V includes nearest-neighbor mesocolumnar interactions. τ is usually considered to be on the order of 5-10 ms.
Nearest-neighbor interactions among columns give dispersion relations that were used to calculate speeds of mental visual rotation (Ingber, , 1983. The wave equation cited by EEG theorists, permitting fits of SMNI to EEG data (Ingber, 1995a), was derived using the variational principle applied to the SMNI Lagrangian.
This creates an audit trail from synaptic parameters to the statistically averaged regional Lagrangian.

Verification of basic SMNI Hypothesis
The core SMNI hypothesis first developed circa 1980 (Ingber, 1981(Ingber, , 1983 is that highly synchronous patterns of neuronal firings in fact process high-level information. Only since 2012 has this hypothesis been verified experimentally (Asher, 2012;Salazar et al., 2012).

Three Basic SMNI Models
Three basic models were developed by slightly changing the background firing component of the columnar-averaged efficacies B G G within experimental ranges, which modify F G threshold factors to yield in the conditional probability: (a) case EC, dominant excitation subsequent firings (b) case IC, inhibitory subsequent firings (c) case BC, balanced between EC and IC This is consistent with experimental evidence of shifts in background synaptic activity under conditions of selective attention (Briggs et al., 2013;Mountcastle et al., 1981), This enables a Centering Mechanism (CM) on case BC, giving BC , wherein the numerator of F G only has terms proportional to M E , M I and M ‡E , i.e., zeroing other constant terms by resetting the background parameters B G G , still within experimental ranges. This brings in a maximum number of minima into the physical firing M G -space, due to the minima of the new numerator in being in a parabolic trough defined by about which nonlinearities develop multiple minima identified with STM phenomena.
In current projects a Dynamic CM (DCM) model is used, resetting B G G every few epochs of τ . Such changes in background synaptic activity on such time scales are seen during attentional tasks (Briggs et al., 2013).
Spline-Laplacian transformations on the EEG potential Φ are proportional to the SMNI M G firing variables at each electrode site. The electric potential Φ is experimentally measured by EEG, not A, but both are due to the same currents I. Therefore, A is linearly proportional to Φ with a simple scaling factor included as a parameter in fits to data. Additional parameterization of background synaptic parameters, B G G and B ‡E E , modify previous work. The A model outperformed the no-A model, where the no-A model simply has used A-nondependent synaptic parameters. Cost functions with an |A| model were much worse than either the A model or the no-A model. Runs with different signs on the drift and on the absolute value of the drift also gave much higher cost functions than the A model.

PATHINT STM
The evolution of a Balanced Centered model (BC) after 500 foldings of ∆t = 0.01, 5 unit of relaxation time τ , exhibits the existence of ten well developed peaks. These peaks are identified with possible trappings of firing patterns. This describes the "7 ± 2" rule, as calculated by SMNI PATHINT in Fig. 2 (Ingber and Nunez, 1995). Figure 2 illustrates SMNI STM Model BC at the evolution at 5τ (Ingber and Nunez, 1995).

PATHINT STM Visual
The evolution of a Balanced Centered Visual model (BCV) after 1000 foldings of ∆t = 0.01, 10 unit of relaxation time τ , exhibits the existence of four well developed peaks. These peaks are identified with possible trappings of firing patterns. Other peaks at lower scales are clearly present, numbering on the same order as in the BC' model, as the strength in the original peaks dissipates throughout firing space, but these are much smaller and therefore much less probable to be accessed. This describes the "4 ± 2" rule for visual STM, as calculated by SMNI PATHINT in Fig. 3 (Ingber and Nunez, 1995). Figure 3 illustrates SMNI STM Model BCV at the evolution at 10τ (Ingber and Nunez, 1995).
http://www.astrocyte.info claims they are the most numerous cells in the human brain. Unlike the previous ideology of astrocytes being "filler" cells, they are very active in the central nervous system and greatly outnumber neurons, Glutamate release from astrocytes through a Ca 2+ -dependent mechanism can activate receptors at the presynaptic terminals. Intercellular calcium waves (ICWs) may travel over hundreds of astrocytes propagating over many neuronal synapses. ICWs contribute to control synaptic activity. Glutamate is released in a regenerative manner, with subsequent cells that are involved in the calcium wave releasing additional glutamate (Ross, 2012).
[Ca 2+ ] (concentrations of Ca 2+ ) affect increased release probabilities at synaptic sites, by enhancing the release of gliotransmitters. (Free Ca 2+ waves are considered here, not intracellular astrocyte calcium waves in situ which also increase neuronal firings.) These free regenerative Ca 2+ waves, arising from astrocyte-neuron interactions, couple to the magnetic vector potential A produced by highly synchronous collective firings, e.g., during selective attention tasks, as measured by EEG.

2.8.1
Canonical Momentum Π = p + qA As derived in the Feynman (midpoint) representation of the path integral, the canonical momentum, Π, defines the dynamics of a moving particle with momentum p in an electromagnetic field. In SI units, where q = −2e for Ca 2+ , e is the magnitude of the charge of an electron = 1.6×10 −19 C (Coulomb), and A is the electromagnetic vector potential. (In Gaussian units Π = p + qA/c, where c is the speed of light.) A represents three components of a 4-vector.

Vector Potential of Wire
A columnar firing state is modeled as a wire/neuron with current I measured in A = Amperes = C/s, along a length z observed from a perpendicular distance r from a line of thickness r 0 . If far-field retardation effects are neglected, this yields where µ is the magnetic permeability in vacuum = 4π10 −7 H/m (Henry/meter). Note the insensitive log dependence on distance. The contribution to A includes minicolumnar lines of current from hundreds to thousands of macrocolumns, within a region not so large to include many convolutions, but still contributing to large synchronous bursts of EEG.
Electric E and magnetic B fields, derivatives of A with respect to r, do not possess this logarithmic insensitivity to distance, and therefore they do not linearly accumulate strength within and across macrocolumns.
Estimates of contributions from synchronous firings to P300 measured on the scalp are tens of thousands of macrocolumns spanning 100 to 100's of cm 2 . Electric fields generated from a minicolumn may fall by half within 5-10 mm, the range of several macrocolumns.
There are other possible sources of magnetic vector potentials not described as wires with currents (Majhi and Ghosh, 2018). Their net effects plausibly would be included the vector magnetic potential of net synchrous firings, but not their functional forms as derived here.

Effects of Vector Potential on Momenta
The momentum p for a Ca 2+ ion with mass m = 6.6 × 10 −26 kg, speed on the order of 50 µm/s to 100 µm/s, is on the order of 10 −30 kg-m/s. Molar concentrations of Ca 2+ waves, comprised of tens of thousands of free ions representing about 1% of a released set, most being buffered, are within a range of about 100 µm to as much as 250 µm, with a duration of more than 500 ms, and with [Ca 2+ ] ranging from 0.1-5 µM (µM = 10 −3 mol/m 3 ).
The magnitude of the current is taken from experimental data on dipole moments Q = |I|ẑ whereẑ is the direction of the current I with the dipole spread over z. Q ranges from 1 pA-m = 10 −12 A-m for a pyramidal neuron (Murakami and Okada, 2006), to 10 −9 A-m for larger neocortical mass (Nunez and Srinivasan, 2006). These currents give rise to qA ≈ 10 −28 kg-m/s. The velocity of a Ca 2+ wave can be ≈20-50 µm/s. In neocortex, a typical Ca 2+ wave of 1000 ions, with total mass m = 6.655 × 10 −23 kg times a speed of ≈20-50 µm/s, gives p ≈ 10 −27 kg-m/s.

Reasonable Estimates
Estimates used here for Q come from experimental data. These include shielding and material effects. When coherent activity among many macrocolumns associated with STM is considered, |A| may be much larger. Since Ca 2+ waves influence synaptic activity, there is direct coherence between these waves and the activity of A. Classical physics calculates qA from macroscopic EEG to be on the order of 10 −28 kg-m/s, while the momentum p of a Ca 2+ ion is on the order of 10 −30 kg-m/s. This numerical comparison demonstrates the possible importance of the influence of A on p at classical scales.

Model of Models (MOM)
Deep Learning (DL) has invigorated AI approaches to parsing data in complex systems, often to develop control processes of these systems. A couple of decades ago, Neural Net AI approaches fell out of favor when concerns were apparent that such approaches offered little guidance to explain the "why" or "how" such algorithms worked to process data, e.g., contexts which were deemed important to deal with future events and outliers, etc.
The success of DL has overshadowed these concerns. However, that should not diminish their importance, especially if such systems are placed in positions to affect lives and human concerns; humans are ultimately responsible for structures they build.
An approach to dealing with these concerns can be called Model of Models (MOM). An argument in favor of MOM is that humans over thousands of years have developed models of reality across many disciplines, e.g., ranging over Physics, Biology, Mathematics, Economics, etc.
A good use of DL might be to process data for a given system in terms of a collection of models, then again use DL to process the models over the same data to determine a superior model of models (MOM). Eventually, large DL (quantum) machines could possess a database of hundreds or thousands of models across many disciplines, and directly find the best (hybrid) MOM for a given system.
In particular, SMNI offers a reasonable model upon which to further develop MOM, wherein multiple scales of observed interactions are developed. This is just one example of how physics modeling and computational physics can be used to better understand complex systems.

Ideas by Statistical Mechanics
A project sympathetic to this MOM context was proposed as Ideas by Statistical Mechanics (ISM) (Ingber, 2006a(Ingber, , 2007(Ingber, , 2008. using ASA (Ingber, 1989(Ingber, , 1993a(Ingber, , 2012c to fit parameters of a generic nonlinear multivariate colored-noise Gaussian-Markovian short-time conditional probability distribution to data, useful for many systems. Models developed using ASA have been applied in many contexts across many systems (Ingber, 1993b), including applications to neural networks (Atiya et al., 2003).
Many of these ASA applications have used Ordinal representations of features, to permit parameterization of their inclusion into models, quite similar in spirit to DL.
ASA can be used again in the expanded context of MOM. This is suggested as a first step in a new discipline to which MOM is to be applied, to help develop a range of parameters useful for DL, as DL by itself may get stuck in non-ideal local minima of the importance-sampled space. Then, after a reasonable range of models is found, DL can take over to permit much more efficient and accurate development of MOM for a given discipline/system.

Importance Sampling
Nonlinear and/or stochastic systems often require importance-sampling algorithms to scan or to fit parameters. Methods of simulated annealing (SA) are often used. Proper annealing (not "quenching") possesses a proof of finding the deepest minimum in searches.
The ASA code is open-source software, and can be downloaded and used without any cost or registration at https://www.ingber.com/#ASA (Ingber, 1993a(Ingber, , 2012c. This algorithm fits empirical data to a cost function over a D-dimensional parameter space, adapting for varying sensitivities of parameters during the fit. This ASA algorithm is faster than fast Cauchy annealing, which has schedule T i = T 0 /k, and much faster than Boltzmann annealing, which has schedule T i = T 0 / ln k (Ingber, 1989).

Outline of ASA Algorithm
For parameters in terms of parameter "temperatures" The default ASA uses the same type of annealing schedule for the acceptance function h as used for the generating function g. All default functions in ASA can be overridden with user-defined functions.

Path Integral in Stratonovich (Midpoint) Representation
The path integral in the Feynman (midpoint) representation is used to examine discretization issues in time-dependent nonlinear systems (Langouche et al., 1979(Langouche et al., , 1982Schulman, 1981). (N.b. g † in DM implies a prepoint evaluation.) Unless explicitly stated, the Einstein summation convention is used which implies repeated indices signify summation; bars | . . . | imply no summation.
Non-constant diffusions add terms to drifts, and a Riemannian-curvature potential R/6 is induced for dimension > 1 in the Stratonovich/Feynman discretization.

Path Integral in Ito (Prepoint) Representation
In the Ito (prepoint) representation: Here the diagonal diffusions are g |GG| and the drifts are g G .

Path-Integral Riemannian Geometry
The midpoint derivation derives a Riemannian geometry with metric defined by the inverse of the covariance matrix and where R is the Riemannian curvature An Ito prepoint discretization for the same probability distribution P gives a simpler algebraic form, but the Lagrangian L so specified does not satisfy a variational principle useful for moderate to large noise. Its variational principle is only useful in the weak-noise limit. This often means that finer meshes are required.

Three Approaches Are Mathematically Equivalent
Three basic different approaches are mathematically equivalent: (a) Fokker-Planck/Chapman-Kolmogorov partial-differential equations (b) Langevin coupled stochastic-differential equations (c) Lagrangian or Hamiltonian path-integrals All three are described here as many researchers are familiar with at least one of these approaches to complex systems. The path-integral approach is useful to define intuitive physical variables from the Lagrangian L in terms of underlying variables M G : Differentiation especially of noisy systems often introduces more noise. The path-integral often gives superior numerical performance because integration is a smoothing process.

Stochastic Differential Equation (SDE)
The Stratonovich (midpoint discretized) Langevin equations can be analyzed in terms of the Wiener process dW i . This can be developed with Gaussian noise η i = dW i /dt, with some care taken in the limit of small dt.
η i represents Gaussian white noise.

Partial Differential Equation (PDE)
The Fokker-Planck, sometimes defines as Chapman-Kolmogorov, partial differential equation is: g G replaces f G in the SDE if the Ito (prepoint discretized) calculus is used. If boundary conditions are added as Lagrange multipliers, these enter as a "potential" V creating a Schrodinger-type equation.

PATHINT/qPATHINT Code
qPATHINT is an N-dimensional code developed to calculate the propagation of quantum variables in the presence of shocks. Many real systems propagate in the presence of sudden changes of state dependent on time. qPATHINT is based on the classical-physics code, PATHINT, which has been useful in several systems across several disciplines. Applications have been made to SMNI and Statistical Mechanics of Financal Markets (SMFM) (Ingber, 2017a,b,c). To numerically calculate the path integral for serial changes in time, standard Monte Carlo techniques generally are not useful. PATHINT was originally developed for this purpose. The PATHINT C code of about 7500 lines of code using the GCC C-compiler was rewritten to use double complex variables instead of double variables, and further developed for arbitrary N dimensions, creating qPATHINT. The outline of the code is described here for classical or quantum systems, using generic coordinates q (Ingber, 2016b(Ingber, , 2017a. The distribution (probabilities for classical systems, wave-functions for quantum systems) can be numerically approximated to a high degree of accuracy using a histogram procedure, developing sums of rectangles of height P i and width ∆q i at points q i .

Shocks
Many real-world systems propagate in the presence of continual "shocks".
In SMNI, collisions occur via regenerative Ca 2+ waves. There also are interactions with changing A due to changing highly synchronous neuronal firings.
In SMFM applications, shocks occur due to future dividends, changes in interest rates, changes in asset distributions, etc.

PATHINT/qPATHINT Histograms
A one-dimensional path-integral in variable q in the prepoint Ito discretization is developed in terms of the kernel/propagator G, for each of its intermediate integrals, as P (q; t + ∆t) = dq [g 1/2 (2π∆t) −1/2 exp(−L ∆t)]P (q ; t) = dq G(q, q ; ∆t)P (q ; t) This yields T ij is a banded matrix representing the Gaussian nature of the short-time probability centered about the drift. Several projects have used this algorithm (Ingber and Nunez, 1995;Ingber and Wilson, 1999;Wolfer, 1983a,b, 1987). Special 2-dimensional codes were developed for specific projects in Statistical Mechanics of Combat (SMC), SMNI and SMFM Ingber and Nunez, 1995).

Meshes For [q]PATHINT
Explicit dependence of L on time t can be included. The mesh ∆q i is strongly dependent on diagonal elements of the diffusion matrix, e.g., This constrains the dependence of the covariance of each variable to be a (nonlinear) function of that variable to present a rectangular underlying mesh. Since integration is inherently a smoothing process (Ingber, 1990), coarser meshes are used relative to the corresponding stochastic differential equation(s) (Wehner and Wolfer, 1983a). By considering the contributions to the first and second moments, conditions on the time and variable meshes can be derived. Thus ∆t can be measured by the diffusion divided by the square of the drift.

Lessons Learned From SMFM and SMNI
SMNI qPATHINT has emphasized the requirement of broad-banded kernels for oscillatory quantum states.
SMFM PATHTREE, and its derived qPATHTREE, is a different options code, based on pathintegral error analyses, permitting a new very fast binary calculation, also applied to nonlinear timedependent systems (Ingber et al., 2001). However, in contrast to the present PATHINT/qPATHINT code that has been generalized to N dimensions, currently an SMFM [q]PATHTREE is only a binary tree with J = 1 and cannot be effectively applied to quantum oscillatory systems (Ingber, 2016b(Ingber, , 2017a.

4.7.1
Calculations At Each Node At Each Time Slice At each node of each time slice, a proposed algorithm is to calculate quantum-scale Ca 2+ wavepacket (2-way) interactions with macroscopic large-scale EEG/A. This entails algorithms: PATHINT using the Classical SMNI Lagrangian qPATHINT using the Quantum Ca 2+ wave-packet Lagrangian Sync in time during P300 attentional tasks. Time/phase relations between classical and quantum systems may be important. ASA-fit synchronized classical-quantum PATHINT-qPATHINT model to EEG data. A is determined experimentally from EEG, and includes all synaptic background B G G effects.

Results Including Quantum Scales
The wave function ψ e describing the interaction of A with p of Ca 2+ wave packets was derived in closed form from the Feynman representation of the path integral using path-integral techniques (Schulten, 1999), modified here to include A.
ψ e (t) = dr 0 ψ 0 ψ F = 1 − iht/(m∆r 2 ) 1 + iht/(m∆r 2 ) 1/4 where ψ 0 is the initial Gaussian packet, ψ F is the free-wave evolution operator,h is the Planck constant, q is the electronic charge of Ca 2+ ions, m is the mass of a wave-packet of 1000 Ca 2+ ions, ∆r 2 is the spatial variance of the wave-packet, the initial momentum is p 0 , and the evolving canonical momentum is Π = p + qA. Detailed calculations show that p of the Ca 2+ wave packet and qA of the EEG field make about equal contributions to Π (Ingber, 2015).
A changes slower than p, so static approximation of A used to derive ψ e and < p > ψ * ψ is reasonable to use within P300 EEG epochs, resetting t = 0 at the onset of each classical EEG measurement (1.953 ms apart), using the current A. This permits tests of interactions across scales in a classical context.

Supercomputer Resources
The XSEDE.org University of California San Diego (UCSD) supercomputer resource is Comet, described at https://portal.xsede.org/sdsc-comet . About 1000 hours of supercomputer CPUs are required for an ASA fit of SMNI to the same EEG data used previously, i.e., from http://physionet.nlm.nih.gov/pn4/erpbci (Citi et al., 2010;Goldberger et al., 2000), using mostly the same codes used previously (Ingber, 2016a). Many such sets of runs are required. Including quantum processes will take even longer.

Results
Using < p > ψ * ψ < p > ψ * ψ was used in classical-physics SMNI fits to EEG data using ASA. Runs using 1M or 100K generated states gave results not much different. Training with ASA used 100K generated states over 12 subjects with and without A, followed by 1000 generated states with the simplex local code contained with ASA. Training and Testing runs on XSEDE.org for this project has taken an equivalent of several months of CPU on the XSEDE.org UCSD platform Comet. These calculations use one additional parameter across all EEG regions to weight the contribution to synaptic background B G G . A is taken to be proportional to the currents measured by EEG, i.e., firings M G . Otherwise, the "zero-fit-parameter" SMNI philosophy was enforced, wherein parameters are picked from experimentally determined values or within experimentally determined ranges (Ingber, 1984).
As with previous studies using this data, results sometimes give Testing cost functions less than the Training cost functions. This reflects on great differences in data, likely from great differences in subjects' contexts, e.g., possibly due to subjects' STM strategies only sometimes including effects calculated here. Further tests of these multiple-scale models with more EEG data are required, and with the PATHINT-qPATHINT coupled algorithm described previously. Table 1 gives recent results on such tests. Cost functions are the effective Action, A ef f , which is L ∆t−log(prefactor), where the prefactor multiplier of the exponential arises from the normalization of the short-time conditional probability distribution and L ∆t is the argument of the exponential factor. Eq. (3) defines the Lagrangian L, and the normalization is defined in DM in Eq. (11).

Quantum Zeno Effects
The quantum-mechanical wave function of the wave packet was shown to "survive" overlaps after multiple collisions, due to their regenerative processes during the observed long durations of hundreds of ms. Thus, Ca 2+ waves may support a Zeno or "bang-bang" effect which may promote long coherence times (Burgarth et al., 2018;Facchi et al., 2004;Facchi and Pascazio, 2008;Giacosa Kozlowski et al., 2015;Muller et al., 2016;Patil et al., 2015;Wu et al., 2012;Zhang et al., 2014). Of course, the Zeno/"bang-bang" effect may exist only in special contexts, since decoherence among particles is known to be very fast, e.g., faster than phase-damping of macroscopic classical particles colliding with quantum particles (Preskill, 2015).
The wave may be perpetuated by the constant collisions of ions as they enter and leave the wave packet due to the regenerative collisions by the Zeno/"bang-bang" effect. qPATHINT can calculate the coherence stability of the wave due to serial shocks.

Survival of Wave Packet
In momentum space, the wave packet φ(p, t) is considered as being "kicked" from p to p + δp. Assume that random repeated kicks of δp result in < δp >≈ 0, and that each kick keeps the variance ∆(p + δp) 2 ≈ ∆(p) 2 . Then, the overlap integral at the moment t of a typical kick between the new and old state is where φ(p + δp, t) is the normalized wave function in p + δp momentum space. A crude estimate is obtained of the survival time amplitude A(t) and survival probability p(t) (Facchi and Pascazio, 2008), These numbers yield: Even many small repeated kicks do not appreciably affect the real part of φ, and these projections do not appreciably destroy the original wave packet, giving a survival probability per kick as p(t) ≈ exp(−2.5 × 10 −7 ) ≈ 1 − 2.5 × 10 −7 . The time-dependent phase terms are sensitive to times of tenths of a sec. These times are prominent in STM and in synchronous neural firings. Therefore, A effects on Ca 2+ wave functions may maximize their influence on STM at frequencies consistent with synchronous EEG during STM.
All these calculations support this model, in contrast to other models of quantum brain processes without such specific calculations and support (Hagan et al., 2002;Hameroff and Penrose, 2013;McKemmish et al., 2009). 6 Quantum Applications

Nano-Robotic Applications
There is the possibility of carrying pharmaceutical products in nanosystems that could affect unbuffered Ca 2+ waves in neocortex (Ingber, 2015). A Ca 2+ -wave momentum-sensor could act like a piezoelectric device.
At the onset of a Ca 2+ wave (on the order of 100's of ms), a change of momentum can be on the order of 10 −30 kg-m/s for a typical Ca 2+ ion. A Ca 2+ wave packet of 1000 ions with onset time of 1 ms, exerts a force on the order of 10 −24 N (1 N ≡ 1 Newton = 1 kg-m/s 2 ). A nano-robot would be attracted to this site, depositing chemicals/drugs that interact with the regenerative Ca 2+ -wave process.
An area of the receptor of the nanosystem of 1 nm 2 would require pressure sensitivity of 10 −6 Pa (1 Pa = 1 pascal = 1 N/m 2 ).
The nano-roboot could be switched on/off at a regional/columnar level by sensitivity to local electric/magnetic fields. Highly synchronous firings during STM processes can be affected by these piezoelectric nanosystems which affect background/noise efficacies via control of Ca 2+ waves. In turn, this would affect the influence of of Ca 2+ waves via the vector potential A, etc.

Free Will
There is interest in researching possible quantum influences on highly synchronous neuronal firings relevant to STM to understand connections to consciousness and "Free Will" (FW) (Ingber, 2016a,b).
If experimental evidence is gained of quantum-level processes of tripartite synaptic interactions with large-scale synchronous neuronal firings, then FW may be established using the Conway-Kochen quantum no-clone "Free Will Theorem" (FWT) Kochen, 2006, 2009).
The essence of FWT is that, since quantum states cannot be cloned, a Ca 2+ quantum wavepacket may not generate a state proven to have previously existed. As explained by the authors Kochen, 2006, 2009), experimenters have specific choices in selecting measurements, which are shared by (twinned) particles, including the choice of any random number generator that might be used to aid such choices. The authors maintain that their proof and description of quantum measurements used is general enough to rule out classical randomness, and that classical determinism cannot be supported by such processes as exist in the quantum world.

Conclusion
The SMNI model has demonstrated it is faithful to experimental data, for EEG recordings under STM experimental paradigms. qPATHINT permits an inclusion of quantum scales in the multiplescale SMNI model, by evolving Ca 2+ wave-packets with momentum p, including serial shocks, interacting with the magnetic vector potential A derived from EEG data, marching forward in time lock-step with experimental EEG data. This presents a time-dependent propagation of interacting quantum and classical scales.
This quantum path-integral algorithm with serial random shocks will be further studied as it can be used for many quantum systems.