https://quantumzeitgeist.com/tensorflow-for-quantumhits-first-birthday/ https://www.tensorflow.org/quantum Two previous XSEDE project-codes are merged for this project, "Electroencephalographic field influence on calcium momentum waves" and "Quantum pathintegral qPATHTREE and qPATHINT algorithms". These Classical codes on Classical computers define a Hybrid Classical-Quantum system. Previous fits to electroencephalographic (EEG) data have been published using quantum wave-packets defining a specific class of (re-)generated Ca2+ ions at tripartite neuron-astrocyte-neuron sites, which influence synaptic interactions [Ingber(2018)]. Since 2011 [Ingber(2011), Ingber(2012a)], During tasks requiring short-term memory (STM), classical as well as quantum calculations are consistent with interactions between momenta p of these wave-packets with a magnetic vector potential A generated by highly synchronous firings of many thousands of neocortical neurons. An anayltic (closed-form) path-integral calculation of the quantum process, in terms of a wavefunction with expectation value p interacting with A, defines quantum interactions coupled to a macroscopic system [Ingber (2017a)]. It is important to further address this system using these to enable realistic inclusion of shocks to the wave-packet. Due to the regenerative process of the wave-packet, e.g., due to collisions between Ca2+ ions in the wave-packet causing some ions to leave the wave-packet during its hundreds of msec lifetime, or as new ions enter from the astrocyte processes, the wavefunction is repeatedly projected onto quantum subspaces. These may be considered as random processes [Ross (2012)]. PATHTREE/qPATHTREE and PATHINT/qPATHINT can include shocks in the evolution of a short-time probability distribution over thousands of foldings [Ingber(2016a), Ingber(2017a), Ingber (2017b)]. # b) SMNI Statistical Mechanics of Neocortical Interactions (SMNI) was developed in the late 1970's [Ingber(1981), Ingber(1982), Ingber(1983)] and enhanced since by fitting experimental data from short-term memory (STM) and electroencephalography (EEG), including papers on fits to attention data [Ingber(2018)] and affective data [Alakus ~al.(2020) Alakus, Gonen, Turkoglu, Ingber(2021b)]. # c) PATHINT The path-integral defines code that numerically propagates short-time conditional probability distributions [Ingber ~al.(1991)Ingber, Fujio, Wehner, Wehner Wolfer(1983a), Wehner Wolfer(1983b)]. This was generalized to PATHINT in N dimensions, and applied in several disciplines [Ingber(2000a) # d) qPATHINT PATHTREE and PATHINT were generalized to quantum systems, qPATHTREE and qPATHINT [Ingber(2016a), Ingber(2017b), Ingber(2017c)]. A companion paper treats "Hybrid classicalquantum computing: Applications to statistical mechanics of financial markets" [Ingber(2021c)]. # e) Organization of Paper Section 2 further describes SMNI in the context of this project. Section 3 further describes Adaptive Simulated Annealing (ASA) in the context of this project. Section 4 further describes qPATHINT in the context of this project. Section 5 describes how the calculation proceeds between SMNI and qPATHINT. Section 6 describes performance and scaling issues. Section 7 gives motivation to this project. Calculations performed on the Ookami supercomputer at StonyBrook.edu test a closed-form derived solution with and without shocks. Shocks clearly lower cost functions. Section 8 is the Conclusion. # f) Caveat As stated previously in these projects [Ingber(2018)], "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." # II. SMNI Statistical Mechanics of Neocortical Interactions (SMNI) has been developed since 1981, in over 30+ papers, scaling aggregate synaptic interactions describing neuronal firings, scaling minicolumnarmacrocolumnar columns of neurons to mesocolumnar dynamics, and scaling columns of neuronal firings to regional (sensory) macroscopic sites identified in EEG studies [Ingber(1981), Ingber(1982), Ingber(1983), Ingber(1984), Ingber(1985a), Ingber(1994)]. Success of SMNI has been to discover agreement/fits with experimental data from modeled aspects of neocortical interactions, e.g., properties of short-term memory (STM) [Ingber(2012a)], including its capacity (auditory 72 and visual 42), duration, stability, primacy versus recency rule [Ingber(1984), Ingber (1985a), Ingber(1994), Ingber(1995), Ingber(2000b), Ingber(2012a), Ingber Nunez(1995)], EEG dispersion relations [Ingber(1985b)], as well other phenomenon, e.g., Hick's law [Hick(1952), Ingber(1999), Jensen (1987)], nearest-neighbor minicolumnar interactions within macrocolumns calculating rotation of images, etc [Ingber(1982), Ingber(1983), Ingber(1984), Ingber (1985a), Ingber(1994)]. SMNI was also scaled to include mesocolumns across neocortical regions to fit EEG data [Ingber(1997), Ingber(2012a), Ingber(2012b), Ingber (2018), Ingber Nunez(2010)]. # a) XSEDE EEG Project The Extreme Science and Engineering Discovery Environment (XSEDE.org) project since February 2013, "Electroencephalographic field influence on calcium momentum waves," has fit SMNI to EEG data, developing ionic Ca2+ momentum-wave effects among neuron-astrocyte-neuron tripartite synapses modified parameterization of background SMNI parameters. Both classical and quantum physics support the development of the vector magnetic potential of EEG from highly synchronous firings, e.g., as measured during selective attention, as directly interacting with these momentum-waves, creating feedback between these ionic/quantum and macroscopic scales [Ingber(2012a) # i. qPATHINT for SMNI qPATHINT includes quantum regenerative process to Ca2+ wave-packets as reasonable shocks to the waves that typically do not damage its coherence properties. A proof of principal has been published [Ingber(2017b)]. # b) SMNI With A A model of minicolumns as wires supporting neuronal firings, largely from large neocortical excitatory pyramidal cells in layer V (of I-VI), gives rise to currents, in turn giving rise to electric potentials measured as scalp EEG [Ingber(2011), Ingber(2012a), Nunez ~al.(2013)Nunez, Srinivasan, Ingber]. This gives rise to a magnetic vector potential. ?? = ?? 4?? ??log ? ?? ?? 0 ? (1) Note the log-insensitive dependence on distance. In neocortex, ?? ? ?? 0 , where 0 is the magnetic permeability in vacuum = 4??10 ?7 H/m (Henry/meter), where Henry has units of kg-m-C ?2 , the conversion factor from electrical to mechanical variables. For oscillatory waves, the magnetic field B = ? × A and the electric field E = (????/ð??"ð??")? × ? × A do not have this log dependence on distance. Thus, A fields can contribute collectively over large regions of neocortex [Ingber(2012a) The magnitude current is deermined by experimental data on dipole moments Q = |I|z ? where z ? 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 Okada(2006)], to 10 ?9 A-m for larger neocortical mass [Nunez Srinivasan(2006)]. Currents give rise to ??A ? 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 ?? = 6.655 × 10 ?23 kg times a speed of ?20-50 ??m/s, gives p ? 10 ?27 kgm/s. This yields p to be on the same order as ??A. # i. Results Including Quantum Scales The wave-function ð??"ð??" e of 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, not including shocks [Schulten(1999)], modified here to include A. ð??"ð??" e (??) = ? ???? 0 ð??"ð??" 0 ð??"ð??" ?? = ? 1 ? ?????/(??Î?"?? 2 ) 1 + ?????/(??Î?"?? 2 ) ? 1/4 [??Î?"?? 2 {1 + [???/(??Î?"?? 2 )] 2 }] ?1/4 × exp ?? [???(?? 0 +????)??/?? ] 2 2Î?"?? 2 1??????/(??Î?"?? 2 ) 1+[???/(??Î?"?? 2 )] 2 + ?? ?? 0 ??? ? ? ?? (?? 0 +????) 2 ?? 2??? ?(2) PATHINT also has been used with the SMNI Lagrangian ?? for STM for both auditory and visual memory [Ericsson Chase(1982), Zhang Simon(1985)], calculating the stability and duration of STM, the observed 7 ± 2 capacity rule of auditory memory, and the observed 4 ± 2 capacity rule of visual memory [Ingber(2000a), Ingber Nunez(1995)]. # ii. Results using < ?? > ð??"ð??" * ð??"ð??" The author previously used his derived analytic expression for < ?? > ð??"ð??" * ð??"ð??" in classical-physics SMNI fits to EEG data using ASA [Ingber(2016b) Calculations used one additional parameter across all EEG regions to weight the contribution to synaptic background ?? ?? ? ?? . A was proportional to the currents measured by EEG, i.e., firings ?? ?? . The "zero-fit-parameter" SMNI philosophy was enforced, wherein parameters are selected and enforced between experimentally determined ranges [Ingber(1984)]. Sopmetimes Testing cost functions were less than their Training cost functions, a result sometimes found in previous studies using this data. This likely is due to great differences in data, likely from great differences in subjects' contexts, e.g., possibly due to subjects' STM strategies including effects calculated here. ASA optimizations in this project always included "finishing" ASA importance-sampling with the modified Nelder-Mead simplex code included in the ASA code to ensure best precision. iii. Assumptions for quantum SMNI Some assumptions made for this quantum enhancement of SMNI can be determined by future experiments. The quantum wave-function of the Ca 2+ wavepacket was calculated, adding multiple collisions due to their regenerative processes, and it was demonstrated that overlaps with just-previous wave-functions during the observed long durations of hundreds of ms typical of Inclusion of repeated random shocks to the above wave-function ð??"ð??" ?? (??) demonstrated only small effects on the projections of the wave-packet after these shocks, i.e., the survival time was calculated [Facchi Pascazio(2008), Ingber(2018)]. Of course, the Zeno/"bang-bang" effect may exist only in special contexts, given that decoherence among particles is known to be very fast [Preskill(2015)]. Therefore, the constant collisions of Ca 2+ ions as they enter and leave the Ca 2+ wave-packet due to the regenerative process that maintains the wave may perpetuate at least part of the wave, permitting the Zeno/"bang-bang" effect. qPATHINT as used here provides an opportunity to explore the coherence stability of the wave due to serial shocks of this process. # iv. Nano-Robotic Applications It is possible that the above considerations could lead to pharmaceutical products contained in nanosystems that could affect unbuffered Ca 2+ waves in neocortex [Ingber(2015)]. A Ca 2+ -wave momentumsensor could act like a piezoelectric device [Ingber (2018)]. The nano-robot would be sensitive to local electric/magnetic fields. Highly synchronous firings during STM processes could be directed by piezoelectric nanosystems to affect background/noise efficacies via control of Ca 2+ waves. This could affect the influence of Ca 2+ waves via the vector potential A, etc. v. Free will Further qPATHINT calculations could give additional support to researching possible quantum influences on highly synchronous neuronal firings relevant to STM, yielding connections to consciousness and "Free Will" (FW). As described previously [Ingber(2016a), Ingber(2016b)], experimental feedback from quantumlevel processes of tripartite synaptic interactions with large-scale synchronous neuronal firings, recognized as being highly correlated with STM and states of attention, may be established using the quantum no-clone "Free Will Theorem" (FWT) [Conway Kochen(2006), Conway Kochen(2009)]. Ca 2+ quantum wave-packets may generate states proven to have not previously existed, since quantum states cannot be cloned. In this context, these quantum states may be influential in large-scale patterns of synchronous neuronal firings, rendering these patterns as truly new. The FWT may consider these patterns as new decisions not solely based on previous decisions. These considerations are quite different and independent of other philosophical considerations, e.g., as in https://plato.stanford.edu/ entries/qt-consciousness/. Note that only recently has the core SMNI hypothesis since circa 1980 [Ingber(1981), Ingber(1982), Ingber(1983)], that highly synchronous patterns of neuronal firings process high-level information, been verified experimentally [Asher(2012) # c) qPATHINT for SMNI A previous project tested applications of qPATHTREE and qPATHINT. The wave-function ð??"ð??" is numerically propagated from its initial state, growing into a tree of wave-function nodes. At each node, interaction of the of Ca 2+ wave-packet, via its momentum p, with highly synchronous EEG, via its collective magnetic vector potential A, determines changes of timedependent phenomena. Changes occur at microscopic scales, e.g., due to modifications of the regenerative wave-packet as ions leave and contribute to the wavepacket, determining the effect on tripartite contributions to neuron-astrocyte-neuron synaptic activity, affecting both p and A. Such changes also influence macroscopic scales, e.g., changes due to external and internal stimuli affecting synchronous firings, and thereby A. At every time slice, quantum effects on synaptic interactions are determined by expected values of the interactions over probabilities (ð??"ð??" * ð??"ð??") determined by the wave-functions at their nodes. Due to the form of the quantum Lagrangian/Hamiltonian, a multiplicative Gaussian form (with nonlinear drifts and diffusions) is propagated. This permits a straight-forward use of Gaussian quadratures for numerical integration of the expectation of the momenta of the wave-packet, i.e., of < ??(??) > ð??"ð??" * ð??"ð??" . E.g., see https://en.wikipedia.org/wiki/Gaussian_quadrature # d) Comparing EEG Testing Data with Training Data Using EEG data from [ due to the same currents I. Therefore, A is linearly proportional to ? with a simple scaling factor as a parameter in fits. Additional parameterization of background synaptic parameters also are included, ?? ?? ? ?? and ?? ?? ? ??? . # e) Investigation into Spline-Laplacian Transformation As is common practice, codes for the Spline-Laplacian transformations were applied to all electrodes measured on the scalp. However, the author thinks that the transformation should be applied to each Region of neocortex separately (e.g., visual, auditory, somatic, abstract, etc.), since each region typically participates in attention differently. This process is further tested in this project. # III. # ASA Algorithm For parameters ?? ?? ?? ? [?? ?? , ?? ?? ](3) sampling with the random variable ?? ?? ?? ?? ? [?1,1] ?? ??+1 ?? = ?? ?? ?? + ?? ?? (?? ?? ? ?? ?? )(4) the default generating function is ð??"ð??" ?? (??) = ? ?? ??=1 1 2ln (1+1/?? ?? )(|?? ?? |+?? ?? ) ? ? ?? ??=1 ð??"ð??" ?? ?? (?? ?? )(5) in terms of "temperatures" for parameters [Ingber(1989)] ?? ?? = ?? ??0 exp(??? ?? ?? 1/?? )(6) The default ASA uses the same distribution for the annealing schedule for the acceptance function ? used for the generating function ð??"ð??". The ASA default functions can be substituted with user-defined functions [Ingber(1993), Ingber (2012c)]. ASA has been applied to studies of COVID-19, fitting forms like ???? ?? , for variables ?? and parameters ?? and ??, in the drifts and covariances of conditional probability distributions [Ingber(2021d)]. With over 150 OPTIONS, ASA permits robust tuning over many classes of nonlinear stochastic systems. These many OPTIONS help ensure that ASA can be used robustly across many classes of systems. # IV. # Path-Integral Methodology a) Generic Applications Many systems are defined by (a) Fokker-Planck/Chapman -Kolmogorov partial -differential equations, (b) Langevin coupled stochastic-differential equations, and (c) Lagrangian or Hamiltonian pathintegrals. All three such systems of equations are equivalent mathematically, when limits of discretized variables are taken in the defined induced Riemannian geometry of the system due to nonlinear and time-dependent diffusions [Ingber(1982), Ingber(1983) i # . Path-Integral Algorithm In classical physics, the path integral of N variables indexed by ??,, at multiple times indexed by ??, is defined in terms of its Lagrangian ??: ??[?? ?? |?? ?? 0 ]????(??) = ? ? ? ????exp ??min ? ?? ?? 0 ???? ? ??? ??(??(?? 0 ) = ?? 0 )??(??(??) = ?? ?? ) ???? = lim ?? ?? ? ?? +1 ??=1 ð??"ð??" 1/2 ? ?? (2??Î?"??) ?1/2 ???? ?? ?? ??(?????, ?? ?? , ??) = 1 2 (????? ? ð??"ð??" ?? )ð??"ð??" ???? ? (????? ? ? ð??"ð??" ?? ? ) + ??/6 ð??"ð??" ???? ? = (ð??"ð??" ???? ? ) ?1 , ð??"ð??" = det(ð??"ð??" ???? ? )(7) The diagonal diffusion terms are ð??"ð??" ???? and the drift terms are ð??"ð??" ?? . If the diffusions terms are non-constant, there are additional terms in the drift, and in a Riemannian-curvature potential R/6 for dimension > 1 in The path-integral approach is useful to give mathematical support to physically intuitive variables in the Lagrangian ??, Canonical Momenta Indicators (CMI = ? ?? ) were used successfully in neuroscience [Ingber (1996), Ingber(1997), Ingber(1998)], combat analyses [Bowman Ingber(1997)], and financial markets [Ingber(1996), Ingber ]. The histogram algorithm numerically approximates sums of rectangles of height ?? ?? and width Î?"?? ?? at points ?? ?? . E.g., consider a one-dimensional system. In the prepoint Ito discretization, the pathintegral representation is given by the kernel ?? for each of its intermediate integrals, as ??(??; ?? + Î?"??) = ? ???? ? [ð??"ð??" 1/2 (2??Î?"??) ?1/2 exp(??? Î?"??)]??(?? ? ; ??) = ? ???? ? ??(??, ?? ? ; Î?"??)??(?? ? ; ??) ??(??; ??) = ? ?? ??=1 ??(?? ? ?? ?? )?? ?? (??) (t) ??(?? ? ?? ?? ) = 1, (?? ?? ? 1 2 Î?"?? ???1 ) ? ?? ? (?? ?? + 1 2 Î?"?? ?? ); 0, otherwise(9) This yields ?? ?? (?? + Î?"??) = ?? ???? (Î?"??)?? ?? (??)(10) ?? ???? is a banded matrix of the Gaussian short-time probability centered about the (possibly timedependent) drift. Explicit dependence of ?? on time ?? can be included. The mesh Î?"?? ?? is strongly dependent on diagonal elements of the diffusion matrix, e.g., Î?"?? ?? ? (Î?"??ð??"ð??" ???? ) 1/2(11) The covariance of each variable is a (nonlinear) function all variables, presenting a rectangular mesh. Given that integration is a smoothing process [ Ingber(1990)], fitting the data with integrals over the short-time probability distribution permits coarser meshes than the corresponding stochastic differential equation(s). The coarser resolution is appropriate for a numerical solution of the time-dependent path integral. Consideration of first and second moments yields conditions on the time and variable meshes [Wehner Wolfer(1983a)]. A scan of the time slice, Î?"?? ? ?? ? ?1 where ?? ? is the uniform/static Lagrangian, gives most important contributions to the probability distribution ??. # ii. Direct Kernel Evaluation Several projects have used this algorithm [Ingber Nunez(1995) [Ingber Nunez(1995)], and Statistical Mechanics of Financial Markets (SMFM) [Ingber(2000a)]. The 1-dimensional PATHINT code was generalized by the author to N dimensions. Also, a quantum generalization was made, changing all variables and functions to complex variables, encompassing about 7500 lines of PATHINT code. The N-dimensional code was developed for classical and quantum systems [Ingber(2016a), Ingber(2017a), Ingber(2017b)]. # iii. Monte Carlo vs Kernels Path-integral numerical applications often use Monte Carlo techniques [O'Callaghan Miller(2014)]. This includes the author's ASA code using ASA_SAMPLE OPTIONS [Ingber(1993)]. However, this project is concerned with time-sequential serial random shocks, which is not conveniently treated with Monte-Carlo/importance-sampling algorithms. T ij (?t) = 2 ?x i?1 + ?x i x i +?x i /2 x i ??x i?1 /2 dx x j +?x j /2 x j ??x j?1 /2 dx G(x, x ; ?t b) Quantum Path Integral Algorithms i. Scaling Issues qPATHINT has been tested with shocks to Ca 2+ waves [Ingber(2017b)], using the basic code also used for for quantum options on quantum money [Ingber(2017a)]. This has illustrated computational scaling issues, further described in the Performance and Scaling Section. # ii. Imaginary Time Imaginary-time Wick rotations permit imaginarytime to be transformed into real-time. Unfortunately, numerical calculations, after multiple foldings of the path integral, leaves no audit trail back to imaginary time to extract phase information (private communication with several authors of path-integral papers, including Larry Schulman on 18 Nov 2015) [Schulman(1981)]. V. # SMNI with Qpathint This defines a process fitting EEG using SMNI with qPATHINT numerically calculating the Quantum path-integral between EEG epochs. At the beginning of each EEG epoch time is reset (t=0); the wave-function is decohered ("collapsed") by any EEG measurement. Until the end of any EEG epoch, there are multiple calls to SMNI functions to calculate the evolution of the Classical distribution. This replaces the author's Quantum path-integral closed-form time-dependent analytic solution. # VI. # Performance and Scaling Code from a previous XSEDE grant "Electroencephalographic field influence on calcium momentum waves", is used for SMNI EEG fits. Code from a XSEDE previous XSEDE grant "Quantum pathintegral qPATHTREE and qPATHINT algorithm", is used for qPATHINT runs. # a) Scaling Estimates Estimates XSEDE.org's Expanse using 'gcc -O3'. Expanse is described in https://www.sdsc.edu/support/user_guides/expanse i. SMNI 100 ASA-iterations taking 7.12676s yields 0.0712676 sec/ASA-iteration over 2561 EEG epochs. With '-g' the total time is 29.9934s. # ii. qPATHINT The qPATHINT code uses a variable mesh covering 1121 points along the diagonal, with a maximum off-diagonal spread of 27. Corners require considerable CPU time to take care of boundaries. Oscillatory wave functions require a large off-diagonal spread [Ingber(2017b)]. dt=0.0002 secs requires 10 foldings of the distribution. This takes the code 0.0002 secs/qIteraction. With '-g' the code takes 0.004s to run. iii. Projected Hours/Service Units (SUs) for this Project nSubjects x 2 (switch Train/Test) yields a 24array set of 1-node jobs. ASA-iterations x (SMNI_time/ASA-iteration + nEpochs x qIterations x qPATHINT_time/qIteration) yields 100,000 x (0.07 + 2500 x 10 x 0.0002) = 507,000 sec = 140 hr/run = 6 day/run. Time for Gaussian quadratures calculations is not appreciable: https://en.wikipedia.org/wiki/Gaussian_quadrature Maximum duration of a normal job is 2 days. ASA has built in a simple way of ending jobs with printout required to restart, including sets of random numbers generated. # VII. # Closed-Form Calculations Calculations on the Ookami supercomputer at StonyBrook.edu tested the path-integral derived analytic (closed-form) expression [Ingber(2018)] for astrocyte Ca 2+ wave synaptic interactions < ?? > ð??"ð??" * ð??"ð??" = ?? ð??"ð??" * ð??"ð??" ?? ??? 0 = ??A +p 0 ?? 1/2 |Î?"r | ? (???) 2 +(??Î?"r 2 ) 2 ???+??Î?"r 2 ? 1/2 (12) Shocks were inserted into the mass ?? of the wave packet of 1000 ions, using a random number generator that contributed up to 1% of the synaptic contribution due to Ca 2+ wave contributions. I.e., the mass was perturbed as ?? = ??(1 ? ??) + ????, where ?? = 0.1 and ?? is a random number between 0 and 1. The results are given in Table 1. ASA was used with 200000 valid generated states for optimizations, then the modified Nelder-Mead code was used to sometimes gain extra precision with 5000 valid generated states. Training and Testing sets of data were used for 12 subjects, then the Training and Testing sets were exchanged. Runs were done with the added shocks and without these shocks. Each of the 48 runs took about 2 days on the Ookami supercomputer at StonyBrook.edu . As commented previously [Ingber(2018)], "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." Volume XXI Issue VI Version I The number 2561 of EEG epochs is a region of high amplitude of times from 17 to 22 secs after the tasks began, defining epochs to be about 0.002 sec. Table 1: Comparison using analytic derivation without and with shocks. The subject numbers are given as sNN, and a "-X" represents exchanging Training and Testing sets of EEG data. Under Results an "I" represents an improvement of a lower cost functions with shocks versus no-shocks. An "N" represents better results with no-shocks versus shocks. These results show 16 "I"s and 8 "N"s, clearly in favor of shocks. # Conclusion A numerical path-integral methodology is used by SMNI to fit EEG to test quantum evolution of astrocyte-(re-)generated wave-packets of Ca 2+ ions that suffer shocks due collisions and regeneration of free ions. SMNI is generalized with quantum variables using qPATHINT. SMNI has fit experimental data, e.g., STM and EEG recordings under STM experimental paradigms. qPATHINT includes quantum scales in the SMNI model, evolving Ca 2+ wave-packets with momentum p serial shocks, interacting with the magnetic vector potential A due to EEG, via the (p + ??A) interaction at each node at each time slice ??, in time with experimental EEG data. Published pilot studies have given rationales for developing this particular quantum path-integral algorithm, to study serial random shocks that occur in many systems. This quantum version can be used for many quantum systems, which is increasingly important as experimental data is increasing at a rapid pace for many quantum systems. https://docs.ocean.dwavesys.com/projects/hybrid/en/latest/index.htmlhttps://www.rigetti.com/whathttps://azure.microsoft.com/en-us/solutions/hybrid-cloud-app/#overviewhttps://www.ibm.com/it-infrastructure/z/capabilities/hybrid-cloudThese Classicalcomputersoftenrunoptimization program in systems that are described byquantum variables using these companies' Quantumcomputers [Introduction a) Hybrid Computing everal commercial Classical-Quantum computers now can be accessed via the Cloud, e.g., Rigetti, D-Wave, Microsoft, and IBM [Ingber(2021a)]; see ,Ingber(2012b),Ingber(2015),Ingber(2016b),Ingber(2017b),Ingber(2017c),Ingber(2018), 4Citi ~al. (2010)Citi, Poli, Cinel, Goldberger ~al. (2000)Goldberger, Amaral, Glass, Hausdorff, Ivanov, Mark, Mietus, Moody, Peng, Stanley] http://physionet.nlm.nih.gov/pn4/erpbci As was done previously, fitting SMNI to highly synchronous waves (P300) during attention tasks, for each of 12 subjects, it is possible to find 10 Training runs and 10 Testing runs [Ingber(2016b)]. A region of continuous high amplitude of 2561 lines represents times from 17 to 22 secs after the tasks began. Spline-Volume XXI Issue VI Version I ( E ) ( E ) Year 2021 © 2021 Global Journals Hybrid Classical-Quantum Computing: Applications to Statistical Mechanics of Neocortical Interactions © 2021 Global Journals Hybrid Classical-Quantum Computing: Applications to Statistical Mechanics of Neocortical Interactions ## Acknowledgments The author thanks the Extreme Science and Engineering Discovery Environment (XSEDE.org), for yearly supercomputer grants since February 2013. The author also thanks the Ookami supercomputer at StonyBrook.edu for its supercomputer resources for calculations. Note that in accordance with common practice, standalone website URLs are not necessarily referenced at the end, but often appear within the body of text. E.g., see https://libanswers.snhu.edu/faq/8627 https://www.scribbr.com/apa-examples/cite-a-website/ https://libguides.css.edu/APA7thEd/Webpages Computer Modelling 23 (3), 43-53. https:// www. ingber.com/path96_duffing.pdf. 23. [Ingber Wilson (1999)] Ingber, L., Wilson, J. (1999). Volatility of volatility of financial markets. Mathematical Computer Modelling 29 (5), 39-57. https://www.ingber.com/markets99_vol.pdf. 24. [Ingber Wilson(2000)] Ingber, L., Wilson, J. (2000). Statistical mechanics of financial markets: Exponential modifications to black-scholes. * GonenAlakus TTurkoglu] Alakus MGonen ITurkoglu References Références Referencias 1 2020. 2020 Alakus ~al. * Database for an emotion recognition system based on EEG signals andvarious computer games -GAMEEMO 10.1016/j.bspc.2020.101951 Biomedical Signal Processing and Control 60 101951 * Brain's code for visual working memory deciphered in monkeys NIHfunded study JAsher NIH Press Release 2012 NIH, Bethesda, MD Technical Report * Parameterized quantum circuits as machine learning models LloydBenedetti Sack MBenedetti ELloyd SSack MFiorentini 10.1088/2058-9565/ab5944 Quantum Science and Technology 4 4 2019. 2019 Benedetti ~al. * Canonical momenta of nonlinear combat MBowman LIngber Proceedings of the 1997 Simulation Multi-Conference the 1997 Simulation Multi-ConferenceAtlanta, GA; San Diego, CA Bowman Ingber 1997. 1997. 6-10 April 1997 Society for Computer Simulation * Quantum zeno dynamics from general quantum operations FacchiBurgarth Nakazato Pascazio DBurgarth PFacchi HNakazato SPascazio KYuasa arXiv:1809.09570v1 2018. 2018 Aberystwyth U Aberystwyth, UK Technical Report Burgarth ~al.. quant-ph * A threshold for quantum advantage in derivative pricing SChakrabarti RKrishnakumar GMazzola NStamatopoulos SWoerner WZeng arXiv:2012.03819v2 2020 Goldman Sachs New York, NY Technical Report quant-ph * PoliCiti LCiti RPoli CCinel 2010. 2010 Citi ~al. * Documenting, modelling and exploiting P300 amplitude changes due to variable target delays in Donchin's speller 10.1088/1741-2560/7/5/056006 Journal of Neural Engineering 7 056006 * The free will theorem [Conway Kochen ; Conway JKochen S arXiv:quant-ph/0604079 2006. 2006 Princeton U, Princeton, NJ Technical Report quant-ph * The strong free will theorem [Conway Kochen ; Conway JKochen S Notices of the American Mathematical Society 56 2 2009. 2009 * Exceptional memory KEricsson WChase American Scientist 70 1982. 1982 * Unification of dynamical decoupling and the quantum zeno effect LidarFacchi PPascazio] Facchi DLidar SPascazio Physical Review A 69 2004. 2004. 032314 Facchi ~al. * Quantum zeno dynamics: mathematical and physical aspects P[facchi Pascazio ; Facchi SPascazio Journal of Physics A 41 493001 2008. 2008 * Quantum zeno effect by general measurements GGiacosa GPagliara Giacosa Pagliara 2014. 2014 052107 * Physio Bank, PhysioToolkit, and Physio Net: components of a new research resource for complex physiologic signals AmaralGoldberger Glass Hausdorff Ivanov Mark Mietus Moody Stanley]Peng AGoldberger LAmaral LGlass JHausdorff PIvanov RMark JMietus GMoody CKPeng HStanley Circulation 101 23 2000. 2000 Goldberger ~al. * On the rate of gains of information WHick Quarterly Journal Experimental Psychology 34 4 1952 Hick (1952) * Mathematical comparison of combat computer models to exercise data ChenIngber Mondescu Muzzall LIngber CChen RMondescu DMuzzall MRenedo Ingber Fujio LWehner] Ingber HFujio MWehner Mathematical Computer Modelling 64 5 2001. 2001. 1991. 1991 Physical Review E * Optimization of trading physics models of markets LIngber RMondescu Neural Networks in Financial Engineering 12 4 2001 IEEE Transactions Neural Networks * Statistical mechanics of neocortical interactions: High resolution path-integral calculation of short-term memory LIngber PNunez Physical Review E 51 5 1995 Ingber Nunez (1995). Ingber Nunez (2010) * Neocortical dynamics at multiple scales: EEG standing waves, statistical mechanics, and physical analogs LIngber PNunez Ingber Pappalepore LStesiak] Ingber MPappalepore RStesiak Mathematical Biosciences 229 2010. 2014. 2014 Ingber ~al. * Electroencephalographic field influence on calcium momentum waves 10.1016/j.jtbi.2013.11.002 Journal of Theoretical Biology 343 * Pathintegral evolution of chaos embedded in noise: Duffing neocortical analog. Mathematical 44 Ingber Srinivasan LNunez] Ingber RSrinivasan PNunez LIngber Current Progress Journal 1 1 1996. 1996. 2012b Influence of macrocolumnar EEG on ca waves * Adaptive simulated annealing LIngber HOliveira APetraglia LIngber MMachado Stochastic global optimization and its applications with fuzzy adaptive simulated annealing MPetraglia New York Springer 2012c Ingber(2012c). Invited Paper * Calculating consciousness correlates at multiple scales of neocortical interactions LIngber Horizons in Neuroscience Research ACosta EVillalba Nova, Hauppauge, NY 2015 Ingber(2015). Invited paper. * Path-integral quantum PATHTREE and PATHINT algorithms LIngber International Journal of Innovative Research in Information Security 3 5 2016a Ingber(2016a) * Statistical mechanics of neocortical interactions: Large-scale EEG influences on molecular processes LIngber LIngber 10.1016/j.jtbi.2016.02.003.49 International Journal of Innovative Research in Information Security 395 2 2016b. 2017a Journal of Theoretical Biology * Evolution of regenerative ca-ion wave-packet in neuronal-firing fields: Quantum path-integral with serial shocks LIngber International Journal of Innovative Research in Information Security 4 2 2017b Ingber(2017b). Ingber(2017c) * Quantum path-integral qPATHINT algorithm LIngber LIngber 10.2174/1874110X01711010119.52 The Open Cybernetics Systemics Journal 11 7 2017c. 2018. 2018 EEG. Sci * Technological Breakthroughs and Future Business Opportunities in Education, Health, and Outer Space LIngber Hooke, A. 2021a IGI Global Hershey, PA Revisiting our quantum world * Quantum calcium-ion affective influences measured by EEG LIngber International Journal Physics and Space Science 2021b Ingber(2021b) * Hybrid classicalquantum computing: Applications to statistical mechanics of financial markets LIngber 2021: CQCA Physical Studies Institute 2021c Technical Report Ingber(2021c) * Forecasting COVID-19 with importance-sampling and pathintegrals LIngber WSEAS Transactions on Biology and Biomedicine 2021d. 1987 Ingber(2021d) * Non-hermitian dynamics in the quantum zeno limit AJensen CaballeroKozlowski Benitez ]Mekhov WKozlowski SCaballero-Benitez IMekhov arXiv1510.04857 Speed of Information-Processing and Intelligence PVernon Ablex, Norwood, NJ; Oxford, Oxford, UK 1987. 2015. 2015 58 Technical Report Individual differences in the hick paradigm. quant-ph * Functional Integration and Semiclassical Expansions RoekaertsLangouche FTirapegui] Langouche DRoekaerts ETirapegui Muller Gherardini MCaruso] Muller SGherardini FCaruso arXiv:1607.08871v1 Reidel 60 1982. 1982. 2016. 2016 Technical Report Quantum zeno dynamics through stochastic protocols. quant-ph * Contributions of principal neocortical neurons to magnetoencephalography and electroencephalography signals S[murakami Okad ; Murakami YOkada Journal of Physiology 575 3 2006. 2006 * Theoretical and experimental electrophysiology in human neocortex: Multiscale correlates of conscious experience Nunez Srinivasan PIngber] Nunez RSrinivasan LIngber 10.1002/9783527671632.ch06 Multiscale Analysis and Nonlinear Dynamics: From genes to the brain MPesenson New York Wiley 2013. 2013 Nunez ~al. * Electric Fields of the Brain: The Neurophysics of EEG P[nunez Srinivasan ; Nunez RSrinivasan 2006. 2006 Oxford University Press London 2nd Ed * Measurement-induced localization of an ultracold lattice gas [o'callaghanMiller ;O'callaghan MMiller BPatil Chakram YVengalattore] Patil SChakram MVengalattore 10.1103/PhysRevE.89.0421265 Physical Review Letters 89 2014. 2014. 042124. 2015. 2015. 140402 Physical Review E * Quantum mechanics JPreskill Lecture Notes 2015. 2015 Caltech, Pasadena, CA Technical Report * Understanding calcium waves and sparks in central neurons WRoss Nature Reviews Neuroscience 13 2012 Ross (2012 * DotsonSalazar Gray]Bressler RSalazar NDotson SBressler CGray 2012. 2012 Salazar ~al. * Content-specific fronto-parietal synchronization during visual working memory 10.1126/science Science 338 6110 * Techniques and Applications of Path Integration LSchulman 1981. 1981 J. Wiley & Sons New York * Quantum mechanics KSchulten Lecture Notes 1999. 1999 Technical Report * Illinois Urbana, IL * Numerical evaluation of path-integral solutions to fokker-planck equations. I M[wehner Wolfer ; Wehner WWolfer Physical Review A 27 1983a. 1983a * Numerical evaluation of path-integral solutions to fokker-planck equations. II. restricted stochastic processes M[wehner Wolfer ; Wehner WWolfer Physical Review A 28 1983b. 1983b * Numerical evaluation of path integral solutions to fokker-planck equations. III. time and functionally dependent coefficients M[wehner Wolfer ; Wehner WWolfer Physical Review A 35 1987. 1987 * Time-dependent decoherence-free subspace WangWu Yi]Wu SWang LYi X Journal of Physics A 405305 2012. 2012 Wu ~al. * ZhangSimon 1985 * STM capacity for chinese words and idioms: Chunking and acoustical loop hypotheses GZhang HSimon Memory & Cognition 13 1985