Lasko Basnarkov
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Lasko Basnarkov has a total of 45 co-authors. Here are the top co-authors by number of shared publications:
- Florin Avram
University of Pau and Pays de l'Adour, Pau, France, France
8 shared publications - Ljupco Kocarev
Research Center for Computer Science and Information Technologies, Macedonian Academy of Sciences and Arts 1 , Bul. Krste Misirkov 2, 1000 Skopje
6 shared publications - Rim Adenane
Laboratoire d’Analyse, Géométrie et Applications, Département des Mathématiques, Université Ibn-Tofail, Kenitra 14000, Morocco
6 shared publications Ljupco Kocarev
Faculty of Computer Science and Engineering SS. Cyril and Methodius University P.O. Box 393 1000 Skopje Macedonia
5 shared publications- Trifce Sandev
Research Center for Computer Science and Information Technologies, Macedonian Academy of Sciences and Arts, Bul. Krste Misirkov 2, 1000 Skopje, Macedonia
4 shared publications - Igor Tomovski
Macedonian Academy of Sciences and Arts, Skopje, Macedonia
3 shared publications
Latest Publications
- 1 June 2026
- Published by AIP Publishing in Chaos: An Interdisciplinary Journal of Nonlinear Science
- Vol. 36 (6)
The implicit regularizing effect of stochastic resetting in deep learning analysis of anomalous diffusion
Here, we investigate whether stochastic resetting, a technique that periodically reverts training to beneficial checkpoints, recently explored in the context of deep learning with noisy labels, can improve the decoding of anomalous diffusion trajectories, a task made challenging by noise and limited trajectory lengths, which often render subtle differences between diffusion processes indistinguishable. First, we find that incorporating stochastic resets of neural network parameters improves the validation loss across different hyperparameters and noise levels, confirming its applicability in trajectory decoding tasks. Second, we find that the relative benefit of resetting increases with trajectory length, and we offer a mechanistic explanation supported by minibatch-gradient ensemble diagnostics that links this observation to the underlying optimization dynamics. Furthermore, we observe that there exists an optimal resetting probability that yields the best performance, highlighting the importance of tuning this hyperparameter. Building on this insight, we introduce time-varying resetting mechanisms that dynamically adjust the resetting probability during training. Our results show that these mechanisms often match or surpass the performance of the best fixed-resetting probabilities and offer a solid basis for designing effective dynamic resetting strategies for regularization.
- 1 June 2026
- Published by AIP Publishing in Chaos: An Interdisciplinary Journal of Nonlinear Science
- Vol. 36 (6)
The implicit regularizing effect of stochastic resetting in deep learning analysis of anomalous diffusion
Here, we investigate whether stochastic resetting, a technique that periodically reverts training to beneficial checkpoints, recently explored in the context of deep learning with noisy labels, can improve the decoding of anomalous diffusion trajectories, a task made challenging by noise and limited trajectory lengths, which often render subtle differences between diffusion processes indistinguishable. First, we find that incorporating stochastic resets of neural network parameters improves the validation loss across different hyperparameters and noise levels, confirming its applicability in trajectory decoding tasks. Second, we find that the relative benefit of resetting increases with trajectory length, and we offer a mechanistic explanation supported by minibatch-gradient ensemble diagnostics that links this observation to the underlying optimization dynamics. Furthermore, we observe that there exists an optimal resetting probability that yields the best performance, highlighting the importance of tuning this hyperparameter. Building on this insight, we introduce time-varying resetting mechanisms that dynamically adjust the resetting probability during training. Our results show that these mechanisms often match or surpass the performance of the best fixed-resetting probabilities and offer a solid basis for designing effective dynamic resetting strategies for regularization.
- Open Access
- 14 April 2026
- Published by Springer Nature in Systematic Reviews
- Vol. 15 (1)
Optimizing document retrieval using massive text embeddings and LLM prompt engineering
Background: The rapid expansion of digital data poses a unique challenge for retrieving relevant and insightful information efficiently. In particular, the increasing volume of scientific publications has made literature reviews time-consuming. The emergence of large language models (LLMs) offers new opportunities to streamline this process. Methods: This paper explores the use of generative artificial intelligence (GenAI) for query reformulation and evaluates the performance of nine massive text embedding models, varying in size and fine-tuning strategies, in the context of document retrieval. We apply multiple prompt engineering techniques to evaluate the ability of LLMs to generate effective queries, comparing them with human-crafted queries. These are used to retrieve documents utilizing nine embedding models. The evaluation is across five datasets using metrics such as recall, average precision, and rank-based measures. Results: Results show that embedding models fine-tuned for semantic similarity consistently outperform general-purpose models, with UAE Large proving most robust across diverse domains. Furthermore, queries generated using zero-shot and few-shot prompting techniques often surpass the performance of human-formulated queries. Conclusion: These findings highlight the value of integrating LLMs and massive text embeddings to reduce manual effort in literature reviews. GenAI provides a reliable starting point for query formulation, with human input reserved for refinement when needed.
- 14 April 2026
- Published by Springer Nature in Systematic Reviews
- Vol. 15 (1)
Open AccessOptimizing document retrieval using massive text embeddings and LLM prompt engineering
Background: The rapid expansion of digital data poses a unique challenge for retrieving relevant and insightful information efficiently. In particular, the increasing volume of scientific publications has made literature reviews time-consuming. The emergence of large language models (LLMs) offers new opportunities to streamline this process. Methods: This paper explores the use of generative artificial intelligence (GenAI) for query reformulation and evaluates the performance of nine massive text embedding models, varying in size and fine-tuning strategies, in the context of document retrieval. We apply multiple prompt engineering techniques to evaluate the ability of LLMs to generate effective queries, comparing them with human-crafted queries. These are used to retrieve documents utilizing nine embedding models. The evaluation is across five datasets using metrics such as recall, average precision, and rank-based measures. Results: Results show that embedding models fine-tuned for semantic similarity consistently outperform general-purpose models, with UAE Large proving most robust across diverse domains. Furthermore, queries generated using zero-shot and few-shot prompting techniques often surpass the performance of human-formulated queries. Conclusion: These findings highlight the value of integrating LLMs and massive text embeddings to reduce manual effort in literature reviews. GenAI provides a reliable starting point for query formulation, with human input reserved for refinement when needed.
- 31 December 2025
- Published by World Scientific Pub Co Pte Ltd in International Journal of Modern Physics C
- Vol. 38 (04)
Calculation of the effective reproduction number
The basic reproduction number [ ] is a key epidemiological quantity, representing the expected number of new cases obtained from receiving a pathogen from one infected individual in totally susceptible population. Here is studied the effective reproduction number [ ] as its generalized version, by considering that in the population might be present other infected persons, besides the one introduced, which are absent in the definition of [ ]. To calculate the theoretical value of the effective reproduction number, we use an integral formula based on the age of infection. The approach is applied to the classical SIR model and the two-strain SIIS model. It is obtained that such theoretical values deviate from the classical values which assume constant susceptible population, if the susceptible compartment size changes substantially during the typical period of infectiousness. Numerical simulations of the stochastic version of the models give values that match very well to those from the integral formula. The findings suggest that in scenarios where the population has nonnegligible number of infected individuals and the change of the compartment sizes is considerable, one has to be cautious in using the classical results. The numerical simulations also show that the shape of the generation interval distribution depends on whether the susceptible population is changing. Thus, the estimates of the generation interval distribution should be considered carefully as well.
- 31 December 2025
- Published by World Scientific Pub Co Pte Ltd in International Journal of Modern Physics C
- Vol. 38 (04)
Calculation of the effective reproduction number
The basic reproduction number [ ] is a key epidemiological quantity, representing the expected number of new cases obtained from receiving a pathogen from one infected individual in totally susceptible population. Here is studied the effective reproduction number [ ] as its generalized version, by considering that in the population might be present other infected persons, besides the one introduced, which are absent in the definition of [ ]. To calculate the theoretical value of the effective reproduction number, we use an integral formula based on the age of infection. The approach is applied to the classical SIR model and the two-strain SIIS model. It is obtained that such theoretical values deviate from the classical values which assume constant susceptible population, if the susceptible compartment size changes substantially during the typical period of infectiousness. Numerical simulations of the stochastic version of the models give values that match very well to those from the integral formula. The findings suggest that in scenarios where the population has nonnegligible number of infected individuals and the change of the compartment sizes is considerable, one has to be cautious in using the classical results. The numerical simulations also show that the shape of the generation interval distribution depends on whether the susceptible population is changing. Thus, the estimates of the generation interval distribution should be considered carefully as well.
- Open Access
- 21 December 2025
- Published by MDPI AG in Mathematics
- Vol. 14 (1) , 23
The Similarity Between Epidemiologic Strains, Minimal Self-Replicable Siphons, and Autocatalytic Cores in (Chemical) Reaction Networks: Towards a Unifying Framework
Motivation: We aim to study the boundary stability and persistence of positive odes in mathematical epidemiology models by importing structural tools from chemical reaction networks. This is largely a review work, which attempts to congregate the fields of mathematical epidemiology (ME), and chemical reaction networks (CRNs), based on several observations. We started by observing that epidemiologic strains, defined as disjoint blocks in either the Jacobian on the infected variables, or as blocks in the next generating matrix (NGM), coincide in most of the examples we studied, with either the set of critical minimal siphons or with the set of minimal autocatalytic sets (cores) in an underlying CRN. We leveraged this to provide a definition of the disease-free equilibrium (DFE) face/infected set as the union of either all minimal siphons, or of all cores (they always coincide in our examples). Next, we provide a proposed definition of ME models, as models which have a unique boundary fixed point on the DFE face, and for which the Jacobian of the infected subnetwork admits a regular splitting, which allows defining the famous next generating matrix. We then define the interaction graph on minimal siphons (IGMS), whose vertices are minimal siphons, and whose edges indicate the existence of reactions producing species in one siphon from species in another. When this graph is acyclic, we say the model exhibits an Acyclic Minimal Siphon Decomposition (AMSD). For AMSD models whose minimal siphons partition the infection species, we show that the NGM is block triangular after permutation, which implies the classical max structure of the reproduction number for multi-strain models. In conclusion, using irreversible reaction networks, minimal siphons and acyclic siphon decompositions, we provide a natural bridge from CRN to ME. We implement algorithms to compute IGMS and detect AMSD in our Epid-CRN Mathematica package (which already contain modules to identify minimal siphons, criticality, drainability, self-replicability, etc.). Finally, we illustrate on several multi-strain ME examples how the block structure induced by AMSD, and the ME reproduction functions, allow expressing boundary stability and persistence conditions by comparing growth numbers to 1, as customary in ME. Note that while not addressing the general Persistence Conjecture mentioned in the title, our work provides a systematic method for deriving boundary instability conditions for a significant class of structured models.- 21 December 2025
- Published by MDPI AG in Mathematics
- Vol. 14 (1) , 23
Open AccessThe Similarity Between Epidemiologic Strains, Minimal Self-Replicable Siphons, and Autocatalytic Cores in (Chemical) Reaction Networks: Towards a Unifying Framework
Motivation: We aim to study the boundary stability and persistence of positive odes in mathematical epidemiology models by importing structural tools from chemical reaction networks. This is largely a review work, which attempts to congregate the fields of mathematical epidemiology (ME), and chemical reaction networks (CRNs), based on several observations. We started by observing that epidemiologic strains, defined as disjoint blocks in either the Jacobian on the infected variables, or as blocks in the next generating matrix (NGM), coincide in most of the examples we studied, with either the set of critical minimal siphons or with the set of minimal autocatalytic sets (cores) in an underlying CRN. We leveraged this to provide a definition of the disease-free equilibrium (DFE) face/infected set as the union of either all minimal siphons, or of all cores (they always coincide in our examples). Next, we provide a proposed definition of ME models, as models which have a unique boundary fixed point on the DFE face, and for which the Jacobian of the infected subnetwork admits a regular splitting, which allows defining the famous next generating matrix. We then define the interaction graph on minimal siphons (IGMS), whose vertices are minimal siphons, and whose edges indicate the existence of reactions producing species in one siphon from species in another. When this graph is acyclic, we say the model exhibits an Acyclic Minimal Siphon Decomposition (AMSD). For AMSD models whose minimal siphons partition the infection species, we show that the NGM is block triangular after permutation, which implies the classical max structure of the reproduction number for multi-strain models. In conclusion, using irreversible reaction networks, minimal siphons and acyclic siphon decompositions, we provide a natural bridge from CRN to ME. We implement algorithms to compute IGMS and detect AMSD in our Epid-CRN Mathematica package (which already contain modules to identify minimal siphons, criticality, drainability, self-replicability, etc.). Finally, we illustrate on several multi-strain ME examples how the block structure induced by AMSD, and the ME reproduction functions, allow expressing boundary stability and persistence conditions by comparing growth numbers to 1, as customary in ME. Note that while not addressing the general Persistence Conjecture mentioned in the title, our work provides a systematic method for deriving boundary instability conditions for a significant class of structured models. - 19 November 2025
- Published by Elsevier in Chaos, Solitons & Fractals
- Vol. 202
Semi-Markov random walk on complex networks
- 19 November 2025
- Published by Elsevier in Chaos, Solitons & Fractals
- Vol. 202
Semi-Markov random walk on complex networks
- Open Access
- 18 August 2025
- Published by American Physical Society (APS) in Physical Review E
- Vol. 112 (2) , 024120
Response to an external field of a generalized Langevin equation with stochastic resetting of the memory kernel
We study a generalized Langevin equation framework that incorporates stochastic resetting of a truncation power-law memory kernel. The inclusion of stochastic resetting enables the emergence of resonance phenomena even in parameter regimes where conventional settings (without resetting) do not exhibit such behavior. Specifically, we explore the response of the system to an external field under three scenarios: (i) a free particle, (ii) a particle in a harmonic potential, and (iii) the effect of truncation in the memory kernel. In each case, the primary focus is on understanding how the resetting mechanism interacts with standard parameters to induce stochastic resonance. In addition, we explore the effect of resetting on the dielectric loss.
- 18 August 2025
- Published by American Physical Society (APS) in Physical Review E
- Vol. 112 (2) , 024120
Open AccessResponse to an external field of a generalized Langevin equation with stochastic resetting of the memory kernel
We study a generalized Langevin equation framework that incorporates stochastic resetting of a truncation power-law memory kernel. The inclusion of stochastic resetting enables the emergence of resonance phenomena even in parameter regimes where conventional settings (without resetting) do not exhibit such behavior. Specifically, we explore the response of the system to an external field under three scenarios: (i) a free particle, (ii) a particle in a harmonic potential, and (iii) the effect of truncation in the memory kernel. In each case, the primary focus is on understanding how the resetting mechanism interacts with standard parameters to induce stochastic resonance. In addition, we explore the effect of resetting on the dielectric loss.
- Open Access
- 4 August 2024
- Published by MDPI AG in Mathematics
- Vol. 12 (15) , 2425
Some Probabilistic Interpretations Related to the Next-Generation Matrix Theory: A Review with Examples
The fact that the famous basic reproduction number , i.e., the largest eigenvalue of the next generation matrix , sometimes has a probabilistic interpretation is not as well known as it deserves to be. It is well understood that half of this formula, , is a Markovian generating matrix of a continuous-time Markov chain (CTMC) modeling the evolution of one individual on the compartments. It has also been noted that the not well-enough-known rank-one formula for of Arino et al. (2007) may be interpreted as an expected final reward of a CTMC, whose initial distribution is specified by the rank-one factorization of F. Here, we show that for a large class of ODE epidemic models introduced in Avram et al. (2023), besides the rank-one formula, we may also provide an integral renewal representation of with respect to explicit “age kernels” , which have a matrix exponential form.This latter formula may be also interpreted as an expected reward of a probabilistic continuous Markov chain (CTMC) model. Besides the rather extensively studied rank one case, we also provide an extension to a case with several susceptible classes.- 4 August 2024
- Published by MDPI AG in Mathematics
- Vol. 12 (15) , 2425
Open AccessSome Probabilistic Interpretations Related to the Next-Generation Matrix Theory: A Review with Examples
The fact that the famous basic reproduction number , i.e., the largest eigenvalue of the next generation matrix , sometimes has a probabilistic interpretation is not as well known as it deserves to be. It is well understood that half of this formula, , is a Markovian generating matrix of a continuous-time Markov chain (CTMC) modeling the evolution of one individual on the compartments. It has also been noted that the not well-enough-known rank-one formula for of Arino et al. (2007) may be interpreted as an expected final reward of a CTMC, whose initial distribution is specified by the rank-one factorization of F. Here, we show that for a large class of ODE epidemic models introduced in Avram et al. (2023), besides the rank-one formula, we may also provide an integral renewal representation of with respect to explicit “age kernels” , which have a matrix exponential form.This latter formula may be also interpreted as an expected reward of a probabilistic continuous Markov chain (CTMC) model. Besides the rather extensively studied rank one case, we also provide an extension to a case with several susceptible classes. - Open Access
- 21 December 2023
- Published by MDPI AG in Mathematics
- Vol. 12 (1) , 27
Algorithmic Approach for a Unique Definition of the Next-Generation Matrix
The basic reproduction number is a concept which originated in population dynamics, mathematical epidemiology, and ecology and is closely related to the mean number of children in branching processes (reflecting the fact that the phenomena of interest are well approximated via branching processes, at their inception). Despite the very extensive literature around for deterministic epidemic models, we believe there are still aspects which are not fully understood. Foremost is the fact that is not a function of the original ODE model, unless we also include in it a certain gradient decomposition, which is not unique. This is related to the specification of the “infected compartments”, which is also not unique. A second interesting question is whether the extinction probabilities of the natural continuous time Markovian chain approximation of an ODE model around boundary points (disease-free equilibrium and invasion points) are also related to the gradient decomposition. We offer below several new contributions to the literature: (1) A universal algorithmic definition of a gradient decomposition (and hence of the resulting ). (2) A fixed point equation for the extinction probabilities of a stochastic model associated to a deterministic ODE model, which may be expressed in terms of the decomposition. Last but not least, we offer Mathematica scripts and implement them for a large variety of examples, which illustrate that our recipe offers always reasonable results, but that sometimes other reasonable decompositions are available as well.- 21 December 2023
- Published by MDPI AG in Mathematics
- Vol. 12 (1) , 27
Open AccessAlgorithmic Approach for a Unique Definition of the Next-Generation Matrix
The basic reproduction number is a concept which originated in population dynamics, mathematical epidemiology, and ecology and is closely related to the mean number of children in branching processes (reflecting the fact that the phenomena of interest are well approximated via branching processes, at their inception). Despite the very extensive literature around for deterministic epidemic models, we believe there are still aspects which are not fully understood. Foremost is the fact that is not a function of the original ODE model, unless we also include in it a certain gradient decomposition, which is not unique. This is related to the specification of the “infected compartments”, which is also not unique. A second interesting question is whether the extinction probabilities of the natural continuous time Markovian chain approximation of an ODE model around boundary points (disease-free equilibrium and invasion points) are also related to the gradient decomposition. We offer below several new contributions to the literature: (1) A universal algorithmic definition of a gradient decomposition (and hence of the resulting ). (2) A fixed point equation for the extinction probabilities of a stochastic model associated to a deterministic ODE model, which may be expressed in terms of the decomposition. Last but not least, we offer Mathematica scripts and implement them for a large variety of examples, which illustrate that our recipe offers always reasonable results, but that sometimes other reasonable decompositions are available as well. - Open Access
Algorithmic approach for an unique definition of the next generation matrix
The basic reproduction number R0 is a concept which originated in population dynamics, mathematical epidemiology, and ecology and is closely related to the mean number of children in branching processes.We offer below three new contributions to the literature: 1) We order a universal algorithmic definition of a (F, V) gradient decomposition (and hence of the resulting R0), which requires a minimal input from the user, namely the specification of an admissible set of disease/infection variables. We also present examples where other choices may be more reasonable, with more terms in F, or more terms in V . 2) We glean out from the works of Bacaer a fixed point equation (8) for the extinction probabilities of a stochastic model associated to a deterministic ODE model, which may be expressed in terms of the (F, V ) decomposition. The fact that both R0 and the extinction probabilities are functions of (F, V ) underlines the centrality of this pair, which may be viewed as more fundamental than the famous next generation matrix FV^{-1}. 3) We suggest introducing a new concept of sufficient/minimal disease/infection set (sufficient for determining R0). More precisely, our universal recipe of choosing "new infections" once the "infections" are specified suggests focusing on the choice of the latter, which is also not unique. The maximal choice of choosing all compartments which become 0 at the given boundary point seems to always work, but is the least useful for analytic computations, therefore we propose to investigate the minimal one. As a bonus, this idea seems useful for understanding the Jacobian factorization approach for computing R0 . Last but not least, we offer Mathematica scripts and implement them for a large variety of examples, which illustrate that our recipe others always reasonable results, but that sometimes other reasonable (F, V ) decompositions are available as well.
Open AccessAlgorithmic approach for an unique definition of the next generation matrix
The basic reproduction number R0 is a concept which originated in population dynamics, mathematical epidemiology, and ecology and is closely related to the mean number of children in branching processes.We offer below three new contributions to the literature: 1) We order a universal algorithmic definition of a (F, V) gradient decomposition (and hence of the resulting R0), which requires a minimal input from the user, namely the specification of an admissible set of disease/infection variables. We also present examples where other choices may be more reasonable, with more terms in F, or more terms in V . 2) We glean out from the works of Bacaer a fixed point equation (8) for the extinction probabilities of a stochastic model associated to a deterministic ODE model, which may be expressed in terms of the (F, V ) decomposition. The fact that both R0 and the extinction probabilities are functions of (F, V ) underlines the centrality of this pair, which may be viewed as more fundamental than the famous next generation matrix FV^{-1}. 3) We suggest introducing a new concept of sufficient/minimal disease/infection set (sufficient for determining R0). More precisely, our universal recipe of choosing "new infections" once the "infections" are specified suggests focusing on the choice of the latter, which is also not unique. The maximal choice of choosing all compartments which become 0 at the given boundary point seems to always work, but is the least useful for analytic computations, therefore we propose to investigate the minimal one. As a bonus, this idea seems useful for understanding the Jacobian factorization approach for computing R0 . Last but not least, we offer Mathematica scripts and implement them for a large variety of examples, which illustrate that our recipe others always reasonable results, but that sometimes other reasonable (F, V ) decompositions are available as well.
- 30 March 2023
- Published by Springer Nature
- p. 14-25
Comparative Study of Random Walks with One-Step Memory on Complex Networks
- 30 March 2023
- Published by Springer Nature
- p. 14-25
Comparative Study of Random Walks with One-Step Memory on Complex Networks