BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Tapio Ala-Nissila

Tapio Ala-Nissila contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-sciphysics.comp-phcond-mat.softcond-mat.stat-mechcond-mat.mes-hallphysics.app-phphysics.opticsphysics.soc-phPopulations and Evolutionquant-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2026arXiv

Density-dependent sodium-storage mechanisms in hard carbon materials

Understanding the sodium-storage mechanism in hard carbon (HC) anodes is crucial for advancing sodium-ion battery (SIB) technology. However, the intrinsic complexity of HC microstructures and their interactions with sodium remain not fully elucidated. We present a multiscale methodology that integrates grand-canonical Monte Carlo (GCMC) simulations with a machine-learning interatomic potential based on the Gaussian approximation potential (GAP) framework to investigate sodium insertion mechanisms in hard carbons with different levels of porosity, achieved by simulating structural models with densities ranging from 0.7 to 1.9 g cm$^{-3}$. Structural and thermodynamic analyses reveal the interplay between pore size and accessibility and the relative contributions of adsorption, intercalation, and pore filling to the overall storage capacity. Low-density carbons favor pore-filling, achieving extremely high capacities at near-zero voltages, whereas high-density carbons primarily store sodium through adsorption and intercalation, leading to lower but more stable capacities. Intermediate-density carbons ($1.3-1.6$ g cm$^{-3}$) provide the most balanced performance, combining moderate capacity (480 and 310 mAh g$^{-1}$), safe operating voltages, and minimal volume expansion ($<10$\%). These findings establish a direct correlation between carbon density and electrochemical behavior, providing atomic-scale insight into how hard carbon morphology governs sodium-storage. The proposed framework offers a rational design principle for optimizing HC-based SIB anodes toward high energy density and long-term cycling stability.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Driven polymer translocation into a channel: Iso-flux tension propagation theory and Langevin dynamics simulations

Iso-flux tension propagation (IFTP) theory and Langevin dynamics (LD) simulations are employed to study the dynamics of channel-driven polymer translocation in which a polymer translocates into a narrow channel and the monomers in the channel experience a driving force $f_{\rm c}$. In the high driving force limit, regardless of the channel width, IFTP theory predicts $τ\propto f_{\textrm{c}}^β$ for the translocation time, where $β=-1$ is the force scaling exponent. Moreover, LD data show that for a very narrow channel fitting only a single file of monomers, the entropic force due to the subchain inside the channel does not play a significant role in the translocation dynamics, and the force exponent $β= -1$ regardless of the force magnitude. As the channel width increases the number of possible spatial configurations of the subchain inside the channel becomes significant, and the resulting entropic force causes the force exponent to drop below unity.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Modified Poisson-Boltzmann theory for polyelectrolytes in monovalent salt solutions with finite-size ions

We present a soft-potential-enhanced Poisson-Boltzmann (SPB) theory to efficiently capture ion distributions and electrostatic potential around rodlike charged macromolecules. The SPB model is calibrated with a coarse-grained particle-based model for polyelectrolytes (PEs) in monovalent salt solutions as well as compared to a full atomistic molecular dynamics simulations with explicit solvent. We demonstrate that our modification enables the SPB theory to accurately predict monovalent ion distributions around a rodlike PE in a wide range of ion and charge distribution conditions in the weak-coupling regime. These include excess salt concentrations up to 1 M, and ion sizes ranging from small ions, such as Na or Cl, to softer and larger ions with size comparable to the PE diameter. The work provides a simple way to implement an enhancement that effectively captures the influence of ion size and species into the PB theory in the context of PEs in aqueous salt solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport

We develop a neuroevolution-potential (NEP) framework for generating neural network based machine-learning potentials. They are trained using an evolutionary strategy for performing large-scale molecular dynamics (MD) simulations. A descriptor of the atomic environment is constructed based on Chebyshev and Legendre polynomials. The method is implemented in graphic processing units within the open-source GPUMD package, which can attain a computational speed over $10^7$ atom-step per second using one Nvidia Tesla V100. Furthermore, per-atom heat current is available in NEP, which paves the way for efficient and accurate MD simulations of heat transport in materials with strong phonon anharmonicity or spatial disorder, which usually cannot be accurately treated either with traditional empirical potentials or with perturbative methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Temperature-resilient anapole modes associated with TE polarization in semiconductor nanowire

Polarization-dependent scattering anisotropy of cylindrical nanowires has numerous potential applications in, for example, nanoantennas, photothermal therapy, thermophotovoltaics, catalysis, sensing, optical filters and switches. In all these applications, temperature-dependent material properties play an important role and often adversely impact performance depending on the dominance of either radiative or dissipative damping. Here, we employ numerical modeling based on Mie scattering theory to investigate and compare the temperature and polarization-dependent optical anisotropy of metallic (gold, Au) nanowires with indirect (silicon, Si) and direct (gallium arsenide, GaAs) bandgap semiconducting nanowires. Results indicate that plasmonic scattering resonances in semiconductors, within the absorption band, deteriorate with an increase in temperature whereas those occurring away from the absorption band strengthen as a result of the increase in phononic contribution. Indirect-bandgap thin ($20 \,\mathrm{nm}$) Si nanowires present low absorption efficiencies for both the transverse electric (TE, $E_{\perp}$) and magnetic (TM, $E_{\parallel}$) modes, and high scattering efficiencies for the TM mode at shorter wavelengths making them suitable as highly efficient scatterers. Temperature-resilient higher-order anapole modes with their characteristic high absorption and low scattering efficiencies are also observed in the semiconductor nanowires ($r \! = \! 125 \! - \! 130$ nm) for the TE polarization. Herein, the GaAs nanowires present $3 \! - \! 7$ times greater absorption efficiencies compared to the Si nanowires making them especially suitable for temperature-resilient applications such as scanning near-field optical microscopy (SNOM), localized heating, non-invasive sensing or detection that require strong localization of energy in the near field.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Adaptive and optimized COVID-19 vaccination strategies across geographical regions and age groups

We evaluate the efficiency of various heuristic strategies for allocating vaccines against COVID-19 and compare them to strategies found using optimal control theory. Our approach is based on a mathematical model which tracks the spread of disease among different age groups and across different geographical regions, and we introduce a method to combine age-specific contact data to geographical movement data. As a case study, we model the epidemic in the population of mainland Finland utilizing mobility data from a major telecom operator. Our approach allows to determine which geographical regions and age groups should be targeted first in order to minimize the number of deaths. In the scenarios that we test, we find that distributing vaccines demographically and in an age-descending order is not optimal for minimizing deaths and the burden of disease. Instead, more lives could potentially be saved by using strategies which emphasize high-incidence regions and distribute vaccines in parallel to multiple age groups. The level of emphasis that high-incidence regions should be given depends on the overall transmission rate in the population. This observation highlights the importance of updating the vaccination strategy when the effective reproduction number changes due to the general contact patterns changing and new virus variants entering.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Heat transport across graphene/hexagonal-BN tilted grain boundaries from phase-field crystal model and molecular dynamics simulations

We study the interfacial thermal conductance of grain boundaries (GBs) between monolayer graphene and hexagonal boron nitride (h-BN) sheets using a combined atomistic approach. First, realistic samples containing graphene/h-BN GBs with different tilt angles are generated using the phase-field crystal (PFC) model developed recently [P. Hirvonen \textit{et al.}, Phys. Rev. B \textbf{100}, 165412 (2019)] that captures slow diffusive relaxation inaccessible to molecular dynamics (MD) simulations. Then, large-scale MD simulations using the efficient GPUMD package are performed to assess heat transport and rectification properties across the GBs. We find that lattice mismatch between the graphene and h-BN sheets plays a less important role in determining the interfacial thermal conductance as compared to the tilt angle. In addition, we find no significant thermal rectification effects for these GBs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Directing Near-Infrared Photon Transport with Core@Shell Particles

Directing the propagation of near-infrared radiation is a major concern in improving the efficiency of solar cells and thermal insulators. A facile approach to scatter light in the near-infrared region without excessive heating is to embed compact layers with semiconductor particles. The directional scattering by semiconductor@oxide (core@shell) spherical particles (containing Si, InP, TiO$_2$, SiO$_2$, or ZrO$_2$) with a total radius varying from 0.1 to 4.0 μm and in an insulating medium at low volume fraction is investigated using Lorenz-Mie theory and multiscale modelling. The optical response of each layers is calculated under irradiation by the sun or a blackbody emitter at 1180 K. Reflectance efficiency factors of up to 83.7% and 63.9% are achieved for near-infrared solar and blackbody radiation in 200 μm thick compact layers with only 1% volume fraction of bare Si particles with a radius of 0.23 μm and 0.50 μm, respectively. The maximum solar and blackbody efficiency factors of layers containing InP particles was slightly less (80.2% and 60.7% for bare particles with a radius of 0.25 μm and 0.60 μm, respectively). The addition of an oxide coating modifies the surrounding dielectric environment, which improves the solar reflectance efficiency factor to over 90% provided it matches the scattering mode energies with the incident spectral density. The layers are spectrally-sensitive and can be applied as a back or front reflector for solar devices, high temperature thermal insulators, and optical filters in Gradient Heat Flux Sensors for fire safety applications.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Interpretation of apparent thermal conductivity in finite systems from equilibrium molecular dynamics simulations

We propose a way to properly interpret the apparent thermal conductivity obtained for finite systems using equilibrium molecular dynamics simulations (EMD) with fixed or open boundary conditions in the transport direction. In such systems the heat current autocorrelation function develops negative values after a correlation time which is proportional to the length of the simulation cell in the transport direction. Accordingly, the running thermal conductivity develops a maximum value at the same correlation time and eventually decays to zero. By comparing EMD with nonequilibrium molecular dynamics (NEMD) simulations, we conclude that the maximum thermal conductivity from EMD in a system with domain length 2L is equal to the thermal conductivity from NEMD in a system with domain length L. This facilitates the use of nonperiodic-boundary EMD for thermal transport in finite samples in close correspondence to NEMD.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kinetic roughening of the urban skyline

\begin{abstract} In this Letter we follow the asymptotic spatial correlation of buildings&#39; heights $G_{\infty}(r)$ in the whole of the Netherlands \cite{bag3d}, which comprises $\approx 10,000,000$ buildings, for the purpose of recovering its scaling with respect to space and time given respectively by the exponents $r^{2 α}$ and $t^{2β}$. This allows us to identify the universality class of the evolution of the urban skyline seen as a dynamically evolving interface. Two major classes of cities were identified based on the recovered value of $α=0.4$ and $α= 0$, which correspond respectively to the KPZ and the EW universality classes. Picking a discrete model from each of these classes and mapping it to physical rules for constructions in cities we conclude that imposed restrictions on buildings&#39; heights are reflected in the exponent $α$ and thus have implications on how the skyline evolves.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Polymer translocation through nanopore assisted by an environment of active rods

We use a combination of computer simulations and iso-flux tension propagation (IFTP) theory to investigate translocation dynamics of a flexible linear polymer through a nanopore into an environment composed of repulsive active rods in 2D. We demonstrate that the rod activity induces a crowding effect on the polymer, leading to a time-dependent net force that facilitates translocation into the active environment. Incorporating this force into the IFTP theory for pore-driven translocation allows us to characterise translocation dynamics in detail and derive a scaling form for the average translocation time as $\tildeτ \sim \tilde{L}_{\textrm{r}}^ν / \tilde{F}_{\textrm{SP}} $, where $\tilde{L}_{\textrm{r}}$ and $\tilde{F}_{\textrm{SP}}$ are the rod length and self-propelling force acting on the rods, respectively, and $ν$ is the Flory exponent.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pulling a DNA molecule through a nanopore embedded in an anionic membrane: tension propagation coupled to electrostatics

We consider the influence of electrostatic forces on driven translocation dynamics of a flexible polyelectrolyte being pulled through a nanopore by an external force on the head monomer. To this end, we augment the iso-flux tension propagation (IFTP) theory with electrostatics for a negatively charged biopolymer pulled through a nanopore embedded in a similarly charged anionic membrane. We show that for the realistic case such as a single-stranded DNA, the translocation dynamics at low salt where screening is weak and at finite negative membrane charge is unexpectedly accelerated despite the large repulsive electrostatic interactions between the polymer coil on the {\it cis} side and the charged membrane. This is due to the rapid release of the electrostatic potential energy of the coil during translocation, leading to an effectively attractive force that assists end-driven translocation. The speedup results in non-monotonic polymer length and membrane charge dependence of the exponent $α$ characterizing the translocation time $τ\propto N_0^α$ of the polymer with length $N_0$. In the regime of long polymers $N_0\gtrsim500$, the translocation exponent exceeds its upper limit $α=2$ previously observed for the same system without electrostatic interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

State leakage during fast decay and control of a superconducting transmon qubit

Superconducting Josephson junction qubits constitute the main current technology for many applications, including scalable quantum computers and thermal devices. Theoretical modeling of such systems is usually done within the two-level approximation. However, accurate theoretical modeling requires taking into account the influence of the higher excited states without limiting the system to the two-level qubit subspace. Here, we study the dynamics and control of a superconducting transmon using the numerically exact stochastic Liouville-von Neumann equation approach. We focus on the role of state leakage from the ideal two-level subspace for bath induced decay and single-qubit gate operations. We find significant short-time state leakage due to the strong coupling to the bath. We quantify the leakage errors in single-qubit gates and demonstrate their suppression with DRAG control for a five-level transmon in the presence of decoherence. Our results predict the limits of accuracy of the two-level approximation and possible intrinsic constraints in qubit dynamics and control for an experimentally relevant parameter set.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermal conductivity reduction in carbon nanotube by fullerene encapsulation: A molecular dynamics study

Single-walled carbon nanotubes (SWCNTs) in their pristine form have high thermal conductivity whose further improvement has attracted a lot of interest. Some theoretical studies have suggested that the thermal conductivity of a $(10,10)$ SWCNT is dramatically enhanced by C$_{60}$ fullerene encapsulation. However, recent experiments on SWCNT bundles show that fullerene encapsulation leads to a reduction rather than an increase in thermal conductivity. Here, we employ three different molecular dynamics methods to study the influence of C$_{60}$ encapsulation on heat transport in a $(10,10)$ SWCNT. All the three methods consistently predict a reduction of the thermal conductivity of $(10,10)$ SWCNT upon C$_{60}$ encapsulation by $20\%-30\%$, in agreement with experimental results on bundles of SWCNTs. We demonstrate that there is a simulation artifact in the Green-Kubo method which gives anomalously large thermal conductivity from artificial convection. Our results show that the C$_{60}$ molecules conduct little heat compared to the outer SWCNT and reduce the phonon mean free paths of the SWCNT by inducing extra phonon scattering. We also find that the thermal conductivity of a $(10,10)$ SWCNT monotonically decreases with increasing filling ratio of C$_{60}$ molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A minimal Tersoff potential for diamond silicon with improved descriptions of elastic and phonon transport properties

Silicon is an important material and many empirical interatomic potentials have been developed for atomistic simulations of it. Among them, the Tersoff potential and its variants are the most popular ones. However, all the existing Tersoff-like potentials fail to reproduce the experimentally measured thermal conductivity of diamond silicon. Here we propose a modified Tersoff potential and develop an efficient open source code called GPUGA (graphics processing units genetic algorithm) based on the genetic algorithm and use it to fit the potential parameters against energy, virial and force data from quantum density functional theory calculations. This potential, which is implemented in the efficient open source GPUMD (graphics processing units molecular dynamics) code, gives significantly improved descriptions of the thermal conductivity and phonon dispersion of diamond silicon as compared to previous Tersoff potentials and at the same time well reproduces the elastic constants. Furthermore, we find that quantum effects on the thermal conductivity of diamond silicon at room temperature are non-negligible but small: using classical statistics underestimates the thermal conductivity by about 10\% as compared to using quantum statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Influence of Boundaries and Thermostatting on Nonequilibrium Molecular Dynamics Simulations of Heat Conduction in Solids

Nonequilibrium molecular dynamics (NEMD) has been extensively used to study thermal transport at various length scales in many materials. In this method, two local thermostats at different temperatures are used to generate a nonequilibrium steady state with a constant heat flux. Conventionally, the thermal conductivity of a finite system is calculated as the ratio between the heat flux and the temperature gradient extracted from the linear part of the temperature profile away from the local thermostats. Here we show that, with a proper choice of the thermostat, the nonlinear part of the temperature profile should actually not be excluded in thermal transport calculations. We compare NEMD results against those from the atomistic Green&#39;s function method in the ballistic regime, and those from the homogeneous nonequilibrium molecular dynamics method in the ballistic-to-diffusive regime. These comparisons suggest that in all the transport regimes, one should directly calculate the thermal conductance from the temperature difference between the heat source and sink and, if needed, convert it to the thermal conductivity by multiplying it with the system length. Furthermore, we find that the Langevin thermostat outperforms the Nosé-Hoover (chain) thermostat in NEMD simulations because of its stochastic and local nature. We show that this is particularly important for studying asymmetric carbon-based nanostructures, for which the Nosé-Hoover thermostat can produce artifacts leading to unphysical thermal rectification. Our findings are important to obtain correct results from molecular dynamics simulations of nanoscale heat transport as the accuracy of the interatomic potentials is rapidly improving.

0 citations0 reviewsNo code linkedOpen access