Publications

Order by: publication date number of citations

• Codebase release 0.2 for Pychastic

  Radost Waszkiewicz, Maciej Bartczak, Kamil Kolasa, Maciej Lisicki

  SciPost Phys. Codebases 11-r0.2 (2023) · published 3 April 2023 | Toggle abstract · src

  In the last decade, Python-powered physics simulations ecosystem has been growing steadily, allowing greater interoperability, and becoming an important tool in numerical exploration of physical phenomena, particularly in soft matter systems. Driven by the need for fast and precise numerical integration in colloidal dynamics, here we formulate the problem of Brownian Dynamics (BD) in a mathematically consistent formalism of the Itō calculus, and develop a Python package to assist numerical computations. We show that, thanks to the automatic differentiation packages, the classical truncated Taylor-Itō integrators can be implemented without the burden of computing the derivatives of the coefficient functions beforehand. Furthermore, we show how to circumvent the difficulties of BD simulations such as calculations of the divergence of the mobility tensor in the diffusion equation and discontinuous trajectories encountered when working with dynamics on $S^2$ and $SO(3)$. The resulting Python package, Pychastic, is capable of performing BD simulations including hydrodynamic interactions at speeds comparable to dedicated implementations in lower-level programming languages, but with a much simpler end-user interface.

  1 citation

• Quantum rotor in a two-dimensional mesoscopic Bose gas

  Michał Suchorowski, Alina Badamshina, Mikhail Lemeshko, Michał Tomza, Artem G. Volosniev

  SciPost Phys. 18, 059 (2025) · published 19 February 2025 | Toggle abstract · pdf

  We investigate a molecular quantum rotor in a two-dimensional Bose-Einstein condensate. The focus is on studying the angulon quasiparticle concept in the crossover from few- to many-body physics. To this end, we formulate the problem in real space and solve it with a mean-field approach in the frame co-rotating with the impurity. We show that the system starts to feature angulon characteristics when the size of the bosonic cloud is large enough to screen the rotor. More importantly, we demonstrate the departure from the angulon picture for large system sizes or large angular momenta where the properties of the system are determined by collective excitations of the Bose gas.

  1 citation

• Topological defects in a double-mirror quadrupole insulator displace diverging charge

  Isidora Araya Day, Anton R. Akhmerov, Dániel Varjas

  SciPost Phys. Core 5, 053 (2022) · published 15 December 2022 | Toggle abstract · pdf

  We show that topological defects in quadrupole insulators do not host quantized fractional charges, contrary to what their Wannier representation indicates. In particular, we test the charge quantization hypothesis based on the Wannier representation of a disclination and a parametric defect. Since disclinations necessarily strain the lattice and parametric defects require closed curves in parameter space, both defects break four-fold rotation symmetry, even away from their origin. The Wannier representation of the defects is thus determined by local reflection symmetries. Contrary to the hypothesis, we find that the local charge density decays as $\sim 1/r^2$ with distance, leading to a diverging defect charge. Because topological defects are incompatible with four-fold rotation symmetry, we conclude that defect charge quantization is protected by sublattice symmetry, and not higher order topology.

  1 citation

• Codebase release 1.0 for SOLAX

  Louis Thirion, Philipp Hansmann, Pavlo Bilous

  SciPost Phys. Codebases 51-r1.0 (2025) · published 20 February 2025 | Toggle abstract · src

  Numerical modeling of fermionic many-body quantum systems presents similar challenges across various research domains, necessitating universal tools, including state-of-the-art machine learning techniques. Here, we introduce SOLAX, a Python library designed to compute and analyze fermionic quantum systems using the formalism of second quantization. SOLAX provides a modular framework for constructing and manipulating basis sets, quantum states, and operators, facilitating the simulation of electronic structures and determining many-body quantum states in finite-size Hilbert spaces. The library integrates machine learning capabilities to mitigate the exponential growth of Hilbert space dimensions in large quantum clusters. The core low-level functionalities are implemented using the recently developed Python library JAX. Demonstrated through its application to the Single Impurity Anderson Model, SOLAX offers a flexible and powerful tool for researchers addressing the challenges of many-body quantum systems across a broad spectrum of fields, including atomic physics, quantum chemistry, and condensed matter physics.

  1 citation

• Spin-lattice couplings and effect of displacements on magnetic interactions of a skyrmion system PdFe/Ir(111)

  Banasree Sadhukhan, Anders Bergman, Johan Hellsvik, Patrik Thunström, Anna Delin

  SciPost Phys. 18, 064 (2025) · published 21 February 2025 | Toggle abstract · pdf

  PdFe/Ir(111) has attracted tremendous attention for next-generation spintronics devices due to existence of magnetic skyrmions with the external magnetic field. Our density functional theoretical calculations in combination with spin dynamics simulation suggest that the spin spiral phase in fcc stacked PdFe/Ir(111) flips into the skyrmion lattice phase around B$_{ext} \sim$ 6 T. This leads to the microscopic understanding of the thermodynamic and kinetic behaviours affected by the intrinsic spin-lattice couplings (SLCs) in this skyrmion material for magneto-mechanical properties. Here we calculate fully relativistic SLC parameters from first principle computations and investigate the effect of SLC on dynamical magnetic interactions in skyrmion multilayers PdFe/Ir(111). The exchange interactions arising from next nearest-neighbors (NN) in this material are highly frustrated and responsible for enhancing skyrmion stability. We report the larger spin-lattice effect on both dynamical Heisenberg exchanges and Dzyaloshinskii-Moriya interactions for next NN compared to NN which is in contrast with recently observed spin-lattice effect in bulk bcc Fe and CrI$_3$ monolayer. Based on our analysis, we find that the effective measures of SLCs in fcc (hcp) stacking of PdFe/Ir(111) are $\sim 2.71 ( \sim 2.36)$ and $\sim 14.71 ( \sim21.89)$ times stronger for NN and next NN respectively, compared to bcc Fe. The linear regime of displacement for SLC parameters is $≤$ 0.02 Å which is 0.72% of the lattice constant for PdFe/Ir(111). The microscopic understanding of SLCs provided by our current study could help in designing spintronic devices based on thermodynamic properties of skyrmion multilayers.

  1 citation

• In-plane Fermi surface mapping of ZrSiS and HfSiS by de Haas-van Alphen oscillations

  Bruno Gudac, Mihovil Bosnar, Filip Orbanić, Trpimir Ivšić, Ivan Kokanović, Mario Novak

  SciPost Phys. Proc. 11, 019 (2023) · published 6 June 2023 | Toggle abstract · pdf

  We have conducted a comprehensive experimental investigation of the in-plane de Haas van Alphen quantum oscillations in ZrSiS and HfSiS and compared the results to the ab-initio calculations. We found a moderately good agreement for a set of frequencies located between 100 T and 200 T. These could be associated with extremal orbitals of the electron pockets located around the $k_z=\pm\pi/2c$ planes at the edge of the 1st Brillouin zone. In contrast, the experimentally detected low-frequency set of oscillations found around 20 T could not be reproduced by the calculations. Further refinement of the ab-initio calculation has to be conducted to fully capture all the observed features of the Fermi surface in ZrSiS and HfSiS.

  1 citation

• Identifying quantum phase transitions with minimal prior knowledge by unsupervised learning

  Mohamad Ali Marashli, Ho Lai Henry Lam, Hamam Mokayed, Fredrik Sandin, Marcus Liwicki, Ho-Kin Tang, Wing Chi Yu

  SciPost Phys. Core 8, 029 (2025) · published 6 March 2025 | Toggle abstract · pdf

  In this work, we proposed a novel approach for identifying quantum phase transitions in one-dimensional quantum many-body systems using AutoEncoder (AE), an unsupervised machine learning technique, with minimal prior knowledge. The training of the AEs is done with reduced density matrix (RDM) data obtained by Exact Diagonalization (ED) across the entire range of the driving parameter and thus no prior knowledge of the phase diagram is required. With this method, we successfully detect the phase transitions in a wide range of models with multiple phase transitions of different types, including the topological and the Berezinskii-Kosterlitz-Thouless transitions by tracking the changes in the reconstruction loss of the AE. The learned representation of the AE is used to characterize the physical phenomena underlying different quantum phases. Our methodology demonstrates a new approach to studying quantum phase transitions with minimal knowledge, small amount of needed data, and produces compressed representations of the quantum states.

  1 citation

• The Cytnx library for tensor networks

  Kai-Hsin Wu, Chang-Teng Lin, Ke Hsu, Hao-Ti Hung, Manuel Schneider, Chia-Min Chung, Ying-Jer Kao, Pochung Chen

  SciPost Phys. Codebases 53 (2025) · published 13 March 2025 | Toggle abstract · pdf

  We introduce a tensor network library designed for classical and quantum physics simulations called Cytnx (pronounced as sci-tens). This library provides an almost identical interface and syntax for both C++ and Python, allowing users to effortlessly switch between the two languages. Aiming at a quick learning process for new users of tensor network algorithms, the interfaces resemble the popular Python scientific libraries like NumPy, Scipy, and PyTorch. Not only multiple global Abelian symmetries can be easily defined and implemented, Cytnx also provides a new tool called Network that allows users to store large tensor networks and perform tensor network contractions in an optimal order automatically. With the integration of cuQuantum, tensor calculations can also be executed efficiently on GPUs. We present benchmark results for tensor operations on both devices, CPU and GPU. We also discuss features and higher-level interfaces to be added in the future.

  1 citation

• Codebase release 1.0 for Cytnx

  Kai-Hsin Wu, Chang-Teng Lin, Ke Hsu, Hao-Ti Hung, Manuel Schneider, Chia-Min Chung, Ying-Jer Kao, Pochung Chen

  SciPost Phys. Codebases 53-r1.0 (2025) · published 13 March 2025 | Toggle abstract · src

  We introduce a tensor network library designed for classical and quantum physics simulations called Cytnx (pronounced as sci-tens). This library provides an almost identical interface and syntax for both C++ and Python, allowing users to effortlessly switch between the two languages. Aiming at a quick learning process for new users of tensor network algorithms, the interfaces resemble the popular Python scientific libraries like NumPy, Scipy, and PyTorch. Not only multiple global Abelian symmetries can be easily defined and implemented, Cytnx also provides a new tool called Network that allows users to store large tensor networks and perform tensor network contractions in an optimal order automatically. With the integration of cuQuantum, tensor calculations can also be executed efficiently on GPUs. We present benchmark results for tensor operations on both devices, CPU and GPU. We also discuss features and higher-level interfaces to be added in the future.

  1 citation

• Codebase release 0.5 for HarmonicBalance.jl

  Jan Košata, Javier del Pino, Toni L. Heugel, Oded Zilberberg

  SciPost Phys. Codebases 6-r0.5 (2022) · published 24 August 2022 | Toggle abstract · src

  HarmonicBalance.jl is a publicly available Julia package designed to simplify and solve systems of periodic time-dependent nonlinear ordinary differential equations. Time dependence of the system parameters is treated with the harmonic balance method, which approximates the system's behaviour as a set of harmonic terms with slowly-varying amplitudes. Under this approximation, the set of all possible steady-state responses follows from the solution of a polynomial system. In HarmonicBalance.jl, we combine harmonic balance with contemporary implementations of symbolic algebra and the homotopy continuation method to numerically determine all steady-state solutions and their associated fluctuation dynamics. For the exploration of involved steady-state topologies, we provide a simple graphical user interface, allowing for arbitrary solution observables and phase diagrams. HarmonicBalance.jl is a free software available at https://github.com/NonlinearOscillations/HarmonicBalance.jl.

  1 citation