T
+34 946 567 842
F
+34 946 567 843
E
akhmatskaya@bcamath.org
Information of interest
 Orcid: 0000000251367991
My major research interest is in the development of numerical methods for modelling and simulation of complex systems. The main features of my research career have been: Work on the interface of Mathematics, Biology, Chemistry, Physics, Computational Statistics and Computer Science
Broad diapason of research activities from developing theoretical methods and implementing them on high performance, computers to application of software codes for study of important phenomena and processes.
Collaboration with multiinternational academic and industrial partners (Oxford University (UK), Imperial College London (UK), St Andrews University (UK), Berkeley Labs (USA), Accelrys Inc. (USA), NAG (UK), POLYMAT (Spain), energiGUNE (Spain), bioGUNE (Spain) among others).
Work experience:
 Academic: Cornell University, USA; University of Manchester, Imperial College, UK.
 Industrial: Fujitsu Centre for Information Technology, UK. Fujitsu Laboratories of Europe, UK.

Adaptive multistage integrators for optimal energy conservation in molecular simulations
(20161001)We introduce a new Adaptive Integration Approach (AIA) to be used in a wide range of molecular simulations. Given a simulation problem and a step size, the method automatically chooses the optimal scheme out of an available ...

Adaptive Splitting Integrators for Enhancing Sampling Efficiency of Modified Hamiltonian Monte Carlo Methods in Molecular Simulation
(20171024)The modified Hamiltonian Monte Carlo (MHMC) methods, i.e., importance sampling methods that use modified Hamiltonians within a Hybrid Monte Carlo (HMC) framework, often outperform in sampling efficiency standard techniques ...

Assessment of van der Waals inclusive density functional theory methods for layered electroactive materials
(20170101)Computationaldriven materials discovery requires efficient and accurate methods. Density functional theory (DFT) meets these two requirements for many classes of materials. However, DFTbased methods have limitations. One ...

Atomistic Insight into Ion Transport and Conductivity in Ga/AlSubstituted Li$_7$La$_3$Zr$_2$O$_{12}$ Solid Electrolytes
(20180109)Garnetstructured Li$_{7}$La$_{3}$Zr$_{2}$O$_{12}$ is a promising solid electrolyte for nextgeneration solidstate Li batteries. However, sufficiently fast Liion mobility required for battery applications only emerges ...

Causality, 'superluminality', and reshaping in undersized waveguides
(20121231)We analyse the reshaping mechanism leading to apparently 'superluminal' advancement of a pulse traversing an undersized section of a waveguide. For frequencies below the first inelastic threshold (cut off one), there are ...

Classification of resonance Regge trajectories and a modified Mulholland formula
(20111231)We employ a simple potential model to analyze the effects which a Regge trajectory, correlating with a bound or a metastable state at zero angular momentum, has on an integral cross section. A straightforward modification ...

Combining stochastic and deterministic approaches within high efficiency molecular simulations
(20131231)Generalized Shadow Hybrid Monte Carlo (GSHMC) is a method for molecular simulations that rigorously alternates Monte Carlo sampling from a canonical ensemble with integration of trajectories using Molecular Dynamics (MD). ...

Complex angular momentum theory of statetostate integral cross sections: Resonance effects in the $F + HD \to HF(v' = 3) + D$ reaction
(20151231)Statetostate reactive integral cross sections (ICSs) are often affected by quantum mechanical resonances, especially near a reactive threshold. An ICS is usually obtained by summing partial waves at a given value of ...

Computer software for understanding resonances and resonancerelated phenomena in chemical reactions
(20141231)In numerical modelling of chemical reactions one calculates the scattering matrix for the required values of energy and angular momentum. Having done so, one still faces the nontrivial task of extracting detailed information ...

Constant pressure hybrid Monte Carlo simulations in GROMACS
(20141231)Adaptation and implementation of the Generalized Shadow Hybrid Monte Carlo (GSHMC) method for molecular simulation at constant pressure in the NPT ensemble are discussed. The resulting method, termed NPTGSHMC, combines ...

Dynamic modeling of the morphology of latex particles with in situ formation of graft copolymer
(20121231)Modification of the polymerpolymer interfacial tension is a way to tailormake particle morphology of waterborne polymerpolymer hybrids. This allows achieving a broader spectrum of application properties and maximizing ...

Dynamic modeling of the morphology of multiphase waterborne polymer particles
(20131231)Multiphase waterborne polymer particles provide advantages in more demanding applications, and their performance depends on particle morphology. Currently, no dynamic model for the prediction of the development of the ...

Enhancing sampling in atomistic simulations of solid state materials for batteries: a focus on olivine NaFePO$_4$
(20170307)The study of ion transport in electrochemically active materials for energy storage systems requires simulations on quantum, atomistic and mesoscales. The methods accessing these scales not only have to be effective but ...

An even simpler understanding of quantum weak values
(201801)We explain the properties and clarify the meaning of quantum weak values using only the basic notions of elementary quantum mechanics.

Exploring Liion conductivity in cubic, tetragonal and mixedphase Alsubstituted Li7La3Zr2O12 using atomistic simulations and effective medium theory
(20190815)Garnet Li7La3Zr2O12 (LLZO) is a promising solid electrolyte candidate for solidstate Liion batteries, but at room temperature it crystallizes in a poorly Liion conductive tetragonal phase. To this end, partial substitution ...

Extracting Smatrix poles for resonances from numerical scattering data: TypeII Padé reconstruction
(20111231)We present a FORTRAN 77 code for evaluation of resonance pole positions and residues of a numerical scattering matrix element in the complex energy (CE) as well as in the complex angular momentum (CAM) planes. Analytical ...

Hartman effect and weak measurements that are not really weak
(20111231)We show that in wave packet tunneling, localization of the transmitted particle amounts to a quantum measurement of the delay it experiences in the barrier. With no external degree of freedom involved, the envelope of the ...

Impact of competitive processes on controlled radical polymerization
(20141231)The kinetics of radical polymerization have been systematically studied for nearly a century and in general are well understood. However, in light of recent developments in controlled radical polymerization many kinetic ...

Interference mechanism of seemingly superluminal tunneling
(20131231)Apparently "superluminal" transmission, e.g., in quantum tunneling and its variants, occurs via a subtle interference mechanism which allows reconstruction of the entire spatial shape of a wave packet from its front tail. ...

MesoGSHMC: A stochastic algorithm for mesoscale constant temperature simulations
(20111231)We consider the problem of timestepping/sampling for molecular and mesoscale particle dynamics. The aim of the work is to derive numerical timestepping methods that generate samples exactly from the desired target ...

Modified Hamiltonian Monte Carlo for Bayesian Inference
(20190722)The Hamiltonian Monte Carlo (HMC) method has been recognized as a powerful sampling tool in computational statistics. We show that performance of HMC can be significantly improved by incorporating importance sampling and ...

Molecular dynamics simulations of iron and aluminumloaded serum transferrin: Protonation of tyr188 is necessary to prompt metal release
(20121231)Serum transferrin (sTf) carries iron in blood serum and delivers it into cells by receptormediated endocytosis. The protein can also bind other metals, including aluminum. The crystal structures of the metalfree and ...

Multistage splitting integrators for sampling with modified Hamiltonian Monte Carlo methods
(20180717)Modified Hamiltonian Monte Carlo (MHMC) methods combine the ideas behind two popular sampling approaches: Hamiltonian Monte Carlo (HMC) and importance sampling. As in the HMC case, the bulk of the computational cost of ...

Multipletimestepping generalized hybrid Monte Carlo methods
(20141231)Performance of the generalized shadow hybrid Monte Carlo (GSHMC) method [1], which proved to be superior in sampling efficiency over its predecessors [24], molecular dynamics and hybrid Monte Carlo, can be further improved ...

Multiscale Modelling and Simulation of Advanced Battery Materials
(2021)Development of efficient strategies for the rational design of materials involved in the production and storage of renewable energy is essential for accelerating the transition to a lowcarbon economy. To contribute to ...

Multiscale simulations of the antimicrobial peptide maculatin 1.1: Water permeation through disordered aggregates
(20121231)The antimicrobial peptide maculatin 1.1 (M1.1) is an amphipathic αhelix that permeabilizes lipid bilayers. In coarsegrained molecular dynamics (CG MD) simulations, M1.1 has previously been shown to form membranespanning ...

No Time at the End of the Tunnel
(20180821)Modern attosecond experiments seek to provide an insight into a long standing question: “how much time does a tunnelling particle spend in the barrier?” Traditionally, quantum theory relates this duration to the delay ...

NonMarkovian models of the growth of a polymer chain
(20151231)Using simple exactly solvable models, we show that eventdependent time delays may lead to significant nonPoisson effects in the statistics of polymer chain growth. The results are confirmed by stochastic simulation of ...

Numerical Regge pole analysis of resonance structures in elastic, inelastic and reactive statetostate integral cross sections
(20141231)We present a detailed description of a FORTRAN code for evaluation of the resonance contribution a Regge trajectory makes to the integral statetostate cross section (ICS) within a specified range of energies. The ...

Numerical Regge pole analysis of resonance structures in statetostate reactive differential cross sections
(20220419)This is the third (and the last) code in a collection of three programs [Sokolovski et al. (2011), Akhmatskaya et al. (2014)] dedicated to the analysis of numerical data, obtained in an accurate simulation of an atomdiatom ...

Numerical simulation of extreme wave runup during storm events in Tramandaí Beach, Rio Grande do Sul, Brazil
(20151231)We present a high resolution analysis of the interaction of irregular waves with natural and urban structures leading to extreme wave runup. Horizontal runup data, instantaneous flooding maps, and wave propagation beyond ...

On the interfacial lithium dynamics in Li7La3Zr2O12:poly(ethylene oxide) (LiTFSI) composite polymerceramic solid electrolytes under strong polymer phase confinement
(20220527)A better molecularlevel understanding of Li+ diffusion through ceramic/polymer interfaces is key to designing highperformance composite solidstate electrolytes for allsolidstate batteries. By considering as a case ...

An optimal scaling to computationally tractable dimensionless models: Study of latex particles morphology formation
(202002)In modelling of chemical, physical or biological systems it may occur that the coefficients, multiplying various terms in the equation of interest, differ greatly in magnitude, if a particular system of units is used. Such ...

Particle Morphology
(20131231)This chapter discusses the morphology of latex particles obtained mainly by (mini)emulsion polymerisation. It describes some applications of these particles, and discusses the factors that influence the particle morphology. ...

Reexamination of continuous fuzzy measurement on twolevel systems
(20170410)Imposing restrictions on the Feynman paths of the monitored system has in the past been proposed as a universal modelfree approach to continuous quantum measurements. Here we revisit this proposition and demonstrate that ...

Relative frequencies of constrained events in stochastic processes: An analytical approach
(20151231)The stochastic simulation algorithm (SSA) and the corresponding Monte Carlo (MC) method are among the most common approaches for studying stochastic processes. They relies on knowledge of interevent probability density ...

Revealing the Mechanism of Sodium Diffusion in NaxFePO4 Using an Improved Force Field
(20180402)Olivine NaFePO4 is a promising cathode material for Naion batteries. Intermediate phases such as Na0.66FePO4 govern phase stability during intercalationdeintercalation processes, yet little is known about Na+ diffusion ...

'Superluminal paradox' in wave packet propagation and its quantum mechanical resolution
(20131231)We analyse in detail the reshaping mechanism leading to apparently 'superluminal' advancement of a wave packet traversing a classically forbidden region. In the coordinate representation, a barrier is shown to act as an ...

Tunnelling times, Larmor clock, and the elephant in the room
(20210511)A controversy surrounding the “tunnelling time problem” stems from the seeming inability of quantum mechanics to provide, in the usual way, a definition of the duration a particle is supposed to spend in a given region of ...

Unveiling Interfacial LiIon Dynamics in Li7La3Zr2O12/PEO(LiTFSI) Composite PolymerCeramic Solid Electrolytes for AllSolidState Lithium Batteries
(20210623)Unlocking the full potential of solidstate electrolytes (SSEs) is key to enabling safer and moreenergy dense technologies than today’s Liion batteries. In particular, composite materials comprising a conductive, flexible ...

Wigner's friends, tunnelling times and Feynman's "only mystery of quantum mechanics"
(20220113)Recent developments in elementary quantum mechanics have seen a number of extraordinary claims regarding quantum behaviour, and even questioning internal consistency of the theory. These are, we argue, different disguises ...
 An Introduction to Monte Carlo methods: Extensions, Variations and Applications of Monte Carlo Method  BCAM & UPV/EHU Courses 20122013 (BCAM)
 Workshop on Dynamical Systems and Applications, Bilbao, Basque Country  Member of Scientific Committee  1011December, 2013
 9th American Institute of Mathematical Sciences International Conference on Dynamical Systems and Differential Equations (AIMS 2012), Orlando, USA  Special Session Organizer  15 July, 2012
 Minisymposium on Nonlinear dynamics applied to biological modelling, Bilbao, Basque Country  Member of Organizing Committee  09 February, 2012
 The Second HPCGA Workshop, Bilbao, Basque Country  Member of Organizing Committee  1115 March, 2013
 Quantum Days in Bilbao IV: Mathematical Methods in Atomic and Molecular Physics, Bilbao, Spain  Organizer  1516 July, 2014
 10th American Institute of Mathematical Sciences Conference on Dynamical Systems, Differential Equations and Applications (AIMS 2014), Madrid, Spain  Special Session Organizer  0711 July, 2014
 The International Conference on Scientific Computation and Differential Equations, SciCADE 2015, Enhanced Sampling Methods, Minisymposium, Potsdam, Germany  Special Session Organizer  1418 September, 2015
 Quantum Days in Bilbao V: Mathematical Methods in Atomic and Molecular Physics  Satellite of QRS XIII, Bilbao, Spain  Organizer  1314 July, 2015
 BCAM  IMUVA Summer School on Uncertainty Quantification for Applied Problems, Bilbao, Basque Country, Spain  Organizer  47 July, 2016
 131st European Study Group with Industry (ESGI), Bilbao, Basque Country, Spain  Member of Scientific Committee  1519 May, 2017
 Minisymposium "Advances in Modelling Multiscale Systems and their Applications", the VIII International Conference on Coupled Problems in Science and Engineering, Coupled 2019, Sitges, Spain  Organizer  3 June, 2019  5 August, 2019
 International Summer School on Fractional and Other Nonlocal Models, Bilbao, Basque Country, Spain  Organizer  2831 May, 2018
 150th European Study Group with Industry (ESGI), Bilbao, Basque Country, Spain  Member of Scientific Committee  2125 October, 2019
 HPCEUROPA2 fellowship, 2012
 Ikerbasque Fellow, 20092010