Search

The heart beats due to the synchronized contraction of cardiomyocytes triggered by a periodic sequence of electrical signals called action potentials, which originate in the sinoatrial node and spread through the heart's electrical system. A large body of work is devoted to modeling the propagation of the action potential and to reproducing reliably its shape and duration. Connection of computational modeling of cells to macroscopic phenomenological curves such as the electrocardiogram has been also intense, due to its clinical importancce in analyzing cardiovascular diseases. In this work we simulate the dynamics of action potential propagation using the three-variable Fenton-Karma model that can account for both normal and damaged cells through spatially inhomogeneous voltage diffusion coefficient. We monitor the action potential propagation in the cardiac tissue and calculate the pseudo-electrocardiogram that reproduces the R and T waves. The R wave amplitude varies according to a double exponential law as a function of the (spatially homogeneous, for an isotropic tissue) diffusion coefficient. The addition of spatial inhomogeneity in the diffusion coefficient by means of a defected region representing damaged cardiac cells, may result in T-wave inversion in the calculated pseudo-electrocardiogram. The transition from positive to negative polarity of the T-wave is analyzed as a function of the length and the depth of the defected region., Comment: 12 pages, figures, 39 references

We investigate the use of Quantum Neural Networks for discovering and implementing quantum error-correcting codes. Our research showcases the efficacy of Quantum Neural Networks through the successful implementation of the Bit-Flip quantum error-correcting code using a Quantum Autoencoder, effectively correcting bit-flip errors in arbitrary logical qubit states. Additionally, we employ Quantum Neural Networks to restore states impacted by Amplitude Damping by utilizing an approximative 4-qubit error-correcting codeword. Our models required modification to the initially proposed Quantum Neural Network structure to avoid barren plateaus of the cost function and improve training time. Moreover, we propose a strategy that leverages Quantum Neural Networks to discover new encryption protocols tailored for specific quantum channels. This is exemplified by learning to generate logical qubits explicitly for the bit-flip channel. Our modified Quantum Neural Networks consistently outperformed the standard implementations across all tasks.

Artificial intelligence in the form of deep learning is now very popular and directly implemented in many areas of science and technology. In the present work we study time evolution of Discrete Breathers in one-dimensional nonlinear chains using the framework of Convolutional Neural Networks. We focus on differentiating discrete breathers which are localized nonlinear modes from linearized phonon modes. The breathers are localized in space and time-periodic solutions of non-linear discrete lattices while phonons are the linear collective oscillations of interacting atoms and molecules. We show that deep learning neural networks are indeed able not only to distinguish breather from phonon modes but also determine with high accuracy the underlying nonlinear on-site potentials that generate breathers. This work can have extensions to more complex natural systems.

In quantum targeted energy transfer, bosons are transferred from a certain crystal site to an alternative one, utilizing a nonlinear resonance configuration similar to the classical targeted energy transfer. We use a novel computational method based on machine learning algorithms in order to investigate selectivity as well as efficiency of the quantum transfer in the context of a dimer and a trimer system. We find that our method identifies resonant quantum transfer paths that allow boson transfer in unison. The method is readily extensible to larger lattice systems involving nonlinear resonances., Comment: REVTeX 4.2, 9 pages, 7 figures

We focus on chaotic dynamical systems and analyze their time series with the use of autoencoders, i.e., configurations of neural networks that map identical output to input. This analysis results in the determination of the latent space dimension of each system and thus determines the minimal number of nodes necessary to capture the essential information contained in the chaotic time series. The constructed chaotic autoencoders generate similar maximal Lyapunov exponents as the original chaotic systems and thus encompass their essential dynamical information.

The nonlinear dimer obtained through the nonlinear Schr{\"o}dinger equation has been a workhorse for the discovery the role nonlinearity plays in strongly interacting systems. While the analysis of the stationary states demonstrates the onset of a symmetry broken state for some degree of nonlinearity, the full dynamics maps the system into an effective $\phi^4$ model. In this latter context, the self-trapping transition is an initial condition dependent transfer of a classical particle over a barrier set by the nonlinear term. This transition has been investigated analytically and mathematically it is expressed through the hyperbolic limit of Jacobian elliptic functions. The aim of the present work is to recapture this transition through the use of methods of Artificial Intelligence (AI). Specifically, we used a physics motivated machine learning model that is shown to be able to capture the original dynamic self-trapping transition and its dependence on initial conditions. Exploitation of this result in the case of the non-degenerate nonlinear dimer gives additional information on the more general dynamics and helps delineate linear from nonlinear localization. This work shows how AI methods may be embedded in physics and provide useful tools for discovery.

For an ensemble of nonlinear systems that model, for instance, molecules or photonic systems, we propose a method that finds efficiently the configuration that has prescribed transfer properties. Specifically, we use physics-informed machine-learning (PIML) techniques to find the parameters for the efficient transfer of an electron (or photon) to a targeted state in a non-linear dimer. We create a machine learning model containing two variables, $\chi_D$, and $\chi_A$, representing the non-linear terms in the donor and acceptor target system states. We then introduce a data-free physics-informed loss function as $1.0 - P_j$, where $P_j$ is the probability, the electron being in the targeted state, $j$. By minimizing the loss function, we maximize the occupation probability to the targeted state. The method recovers known results in the Targeted Energy Transfer (TET) model, and it is then applied to a more complex system with an additional intermediate state. In this trimer configuration, the PIML approach discovers desired resonant paths from the donor to acceptor units. The proposed PIML method is general and may be used in the chemical design of molecular complexes or engineering design of quantum or photonic systems.

The spread of COVID-19 during the initial phase of the first half of 2020 was curtailed to a larger or lesser extent through measures of social distancing imposed by most countries. In this work, we link directly, through machine learning techniques, infection data at a country level to a single number that signifies social distancing effectiveness. We assume that the standard SIR model gives a reasonable description of the dynamics of spreading, and thus the social distancing aspect can be modeled through time-dependent infection rates that are imposed externally. We use an exponential ansatz to analyze the SIR model, find an exact solution for the time-independent infection rate, and derive a simple first-order differential equation for the time-dependent infection rate as a function of the infected population. Using infected number data from the "first wave" of the infection from eight countries, and through physics-informed machine learning, we extract the degree of linear dependence in social distancing that led to the specific infections. We find that in the two extremes are Greece, with the highest decay slope on one side, and the US on the other with a practically flat "decay". The hierarchy of slopes is compatible with the effectiveness of the pandemic containment in each country. Finally, we train our network with data after the end of the analyzed period, and we make week-long predictions for the current phase of the infection that appear to be very close to the actual infection values.

The COVID-19 pandemic has affected all countries of the world producing a substantial number of fatalities accompanied by a major disruption in their social, financial, and educational organization. The strict disciplinary measures implemented by China were very effective and thus were subsequently adopted by most world countries to various degrees. The infection duration and number of infected persons are of critical importance for the battle against the pandemic. We use the quantitative landscape of the disease spreading in China as a benchmark and utilize infection data from eight countries to estimate the complete evolution of the infection in each of these countries. The analysis predicts successfully both the expected number of daily infections per country and, perhaps more importantly, the duration of the epidemic in each country. Our quantitative approach is based on a Gaussian spreading hypothesis that is shown to arise as a result of imposed measures in a simple dynamical infection model. This may have consequences and shed light in the efficiency of policies once the phenomenon is over.

Chimeras and branching are two archetypical complex phenomena that appear in many physical systems; because of their different intrinsic dynamics, they delineate opposite non-trivial limits in the complexity of wave motion and present severe challenges in predicting chaotic and singular behavior in extended physical systems. We report on the long-term forecasting capability of Long Short-Term Memory (LSTM) and reservoir computing (RC) recurrent neural networks, when they are applied to the spatiotemporal evolution of turbulent chimeras in simulated arrays of coupled superconducting quantum interference devices (SQUIDs) or lasers, and branching in the electronic flow of two-dimensional graphene with random potential. We propose a new method in which we assign one LSTM network to each system node except for "observer" nodes which provide continual "ground truth" measurements as input; we refer to this method as "Observer LSTM" (OLSTM). We demonstrate that even a small number of observers greatly improves the data-driven (model-free) long-term forecasting capability of the LSTM networks and provide the framework for a consistent comparison between the RC and LSTM methods. We find that RC requires smaller training datasets than OLSTMs, but the latter requires fewer observers. Both methods are benchmarked against Feed-Forward neural networks (FNNs), also trained to make predictions with observers (OFNNs).