2025-06-19 | | Total: 10
Twisted multilayer graphene, characterized by its moire patterns arising from inter-layer rotational misalignment, serves as a rich platform for exploring quantum phenomena. While first-principles calculations are computationally prohibitive and empirical interatomic potentials often lack accuracy, machine-learning interatomic potentials (MLIPs) present a promising alternative, offering (near-)DFT accuracy at a significantly reduced computational cost. Despite their success in two-dimensional monolayer materials, MLIPs remain under-explored in twisted multilayer graphene systems. In this work, we develop an Atomic Cluster Expansion (ACE) potential for simulating twisted multilayer graphene and test it on a range of simulation tasks. We propose an approach to construct training and test datasets that incorporate all possible twist angles and local stacking, including incommensurate ones. To achieve this, we generate configurations with periodic boundary conditions suitable for DFT calculations, and then introduce an internal twist and shift within those supercell structures. We further refine the dataset through active learning filtering, guided by Bayesian uncertainty quantification. Our model is validated for accuracy and robustness through a wide range of numerical tests.
Randomisation is widely used in quantum algorithms to reduce the number of quantum gates and ancillary qubits required. A range of randomised algorithms, including eigenstate property estimation by spectral filters, Hamiltonian simulation, and perturbative quantum simulation, though motivated and designed for different applications, share common features in the use of unitary decomposition and Hadamard-test-based implementation. In this work, we start by analysing the role of randomised linear-combination-of-unitaries (LCU) in quantum simulations, and present several quantum circuits that realise the randomised composite LCU. A caveat of randomisation, however, is that the resulting state cannot be deterministically prepared, which often takes an unphysical form UρV† with unitaries U and V. Therefore, randomised LCU algorithms are typically restricted to only estimating the expectation value of a single Pauli operator. To address this, we introduce a quantum instrument that can realise a non-completely-positive map, whose feature of frequent measurement and reset on the ancilla makes it particularly suitable in the fault-tolerant regime. We then show how to construct an unbiased estimator of the effective (unphysical) state U \rho V^\dagger and its generalisation. Moreover, we demonstrate how to effectively realise the state prepared by applying an operator that admits a composite LCU form. Our results reveal a natural connection between randomised LCU algorithms and shadow tomography, thereby allowing simultaneous estimation of many observables efficiently. As a concrete example, we construct the estimators and present the simulation complexity for three use cases of randomised LCU in Hamiltonian simulation and eigenstate preparation tasks.
When dense granular gases are continuously excited under microgravity conditions, spatial inhomogeneities of the particle number density can emerge. A significant share of particles may collect in strongly overpopulated regions, called clusters. This dynamical clustering, or gas-cluster transition, is caused by a complex interplay and balance between the energy influx and dissipation in particle collisions. Particle number density, container geometry, and excitation strength influence this transition. We perform Discrete Element Method (DEM) simulations for ensembles of frictional spheres in a cuboid container and apply the Kolmogorov Smirnov test and a caging criterion to the local packing fraction profiles to detect clusters. Machine learning can be used to study the gas-cluster transition, and can be a promising alternative to identify the state of the system for a given set of system parameters without time-consuming complex DEM simulations. We test various machine learning models and identify the best models to predict dynamical clustering of frictional spheres in a specific experimental geometry.
Physics-informed neural networks (PINNs) are extensively used to represent various physical systems across multiple scientific domains. The same can be said for cardiac electrophysiology, wherein fully-connected neural networks (FCNNs) have been employed to predict the evolution of an action potential in a 2D space following the two-parameter phenomenological Aliev-Panfilov (AP) model. In this paper, the training behaviour of PINNs is investigated to determine optimal hyperparameters to predict the electrophysiological activity of the myocardium in 3D according to the AP model, with the inclusion of boundary and material parameters. An FCNN architecture is employed with the governing partial differential equations in their strong form, which are scaled consistently with normalization of network inputs. The finite element (FE) method is used to generate training data for the network. Numerical examples with varying spatial dimensions and parameterizations are generated using the trained models. The network predicted fields for both the action potential and the recovery variable are compared with the respective FE simulations. Network losses are weighed with individual scalar values. Their effect on training and prediction is studied to arrive at a method of controlling losses during training.
Quantum neural networks (QNNs) are widely employed as ansätze for solving variational problems, where their expressivity directly impacts performance. Yet, accurately characterizing QNN expressivity remains an open challenge, impeding the optimal design of quantum circuits. In this work, we introduce the effective rank, denoted as \kappa, as a novel quantitative measure of expressivity. Specifically, \kappa captures the number of effectively independent parameters among all the variational parameters in a parameterized quantum circuit, thus reflecting the true degrees of freedom contributing to expressivity. Through a systematic analysis considering circuit architecture, input data distributions, and measurement protocols, we demonstrate that \kappa can saturate its theoretical upper bound, d_n=4^n-1, for an n-qubit system when each of the three factors is optimally expressive. This result provides a rigorous framework for assessing QNN expressivity and quantifying their functional capacity. Building on these theoretical insights, and motivated by the vast and highly structured nature of the circuit design space, we employ \kappa as a guiding metric for the automated design of highly expressive quantum circuit configurations. To this end, we develop a reinforcement learning framework featuring a self-attention transformer agent that autonomously explores and optimizes circuit architectures. By integrating theoretical characterization with practical optimization, our work establishes \kappa as a robust tool for quantifying QNN expressivity and demonstrates the effectiveness of reinforcement learning in designing high-performance quantum circuits. This study paves the way for building more expressive QNN architectures, ultimately enhancing the capabilities of quantum machine learning.
Neural network potentials (NNPs) offer a powerful alternative to traditional force fields for molecular dynamics (MD) simulations. Accurate and stable MD simulations, crucial for evaluating material properties, require training data encompassing both low-energy stable structures and high-energy structures. Conventional knowledge distillation (KD) methods fine-tune a pre-trained NNP as a teacher model to generate training data for a student model. However, in material-specific models, this fine-tuning process increases energy barriers, making it difficult to create training data containing high-energy structures. To address this, we propose a novel KD framework that leverages a non-fine-tuned, off-the-shelf pre-trained NNP as a teacher. Its gentler energy landscape facilitates the exploration of a wider range of structures, including the high-energy structures crucial for stable MD simulations. Our framework employs a two-stage training process: first, the student NNP is trained with a dataset generated by the off-the-shelf teacher; then, it is fine-tuned with a smaller, high-accuracy density functional theory (DFT) dataset. We demonstrate the effectiveness of our framework by applying it to both organic (polyethylene glycol) and inorganic (L_{10}GeP_{2}S_{12}) materials, achieving comparable or superior accuracy in reproducing physical properties compared to existing methods. Importantly, our method reduces the number of expensive DFT calculations by 10x compared to existing NNP generation methods, without sacrificing accuracy.
At room temperature, low frequency vibrations at far-infrared frequencies are thermally excited (k_B T > h \nu) and not restricted to harmonic fluctuations around a single potential energy minimum. For folded proteins, these intrinsically anharmonic vibrations can contain information on slow conformational transitions. Recently, we have developed FREquency-SElective ANharmonic (FRESEAN) mode analysis, a method based on time correlation functions that isolates low-frequency vibrational motions from molecular dynamics simulation trajectories without relying on harmonic approximations. We recently showed that low-frequency vibrations obtained from FRESEAN mode analysis are effective collective variables in enhanced sampling simulations of conformational ensembles. However, FRESEAN mode analysis is based on velocity time correlations between all degrees of freedom, which creates computational challenges for large biomolecules. To facilitate future applications, we demonstrate here how coarse-graining of all-atom simulation trajectories can be combined with FRESEAN mode analysis to extract information on low-frequency vibrations at minimal computational cost.
There exist elegant methods of aligning point clouds in \mathbb R^3. Unfortunately, these methods rely on the positive definite property of the Euclidean metric, and do not easily extend to the indefinite Minkowski metric. In this paper, we propose two solutions to the following problem: given inertial reference frames A and B, and given (possibly noisy) measurements of a set of 4-vectors \{v_i\} made in those reference frames with components \{v_{A,i}\} and \{v_{B,i}\}, find the optimal Lorentz transformation \Lambda such that \Lambda v_{A,i}=v_{B,i}. The method we outline is conceptually simple and easily extends to alignment problems in other matrix Lie groups.
Effectively leveraging data from multiple quantum-chemical methods is essential for building machine-learned force fields (MLFFs) that are applicable to a wide range of chemical systems. This study systematically investigates two multi-fidelity training strategies, pre-training/fine-tuning and multi-headed training, to elucidate the mechanisms underpinning their success. We identify key factors driving the efficacy of pre-training followed by fine-tuning, but find that internal representations learned during pre-training are inherently method-specific, requiring adaptation of the model backbone during fine-tuning. Multi-headed models offer an extensible alternative, enabling simultaneous training on multiple fidelities. We demonstrate that a multi-headed model learns method-agnostic representations that allow for accurate predictions across multiple label sources. While this approach introduces a slight accuracy compromise compared to sequential fine-tuning, it unlocks new cost-efficient data generation strategies and paves the way towards developing universal MLFFs.
The computational discovery and design of new crystalline materials, particularly metal-organic frameworks (MOFs), heavily relies on high-quality, computation-ready structural data. However, recent studies have revealed significant error rates within existing MOF databases, posing a critical data problem that hinders efficient high-throughput computational screening. While rule-based algorithms like MOSAEC, MOFChecker, and the Chen and Manz method (Chen-Manz) have been developed to address this, they often suffer from inherent limitations and misclassification of structures. To overcome this challenge, we introduce MOFClassifier, a novel machine learning approach built upon a positive-unlabeled crystal graph convolutional neural network (PU-CGCNN) model. MOFClassifier learns intricate patterns from perfect crys-tal structures to predict a crystal-likeness score (CLscore), effectively classifying MOFs as computation-ready. Our model achieves a ROC value of 0.979 (previous best 0.912) and, importantly, can identify subtle structural and chemical errors that are fundamentally undetectable by current rule-based methods. By accurately recovering previously misclassified false-negative structures, MOFClassifier reduces the risk of overlooking promising material candidates in large-scale computational screening efforts. This user-friendly tool is freely available and has been integrated into the preparation workflow for the updated CoRE MOF DB 2025 v1, contributing to accelerated computational discovery of MOF materials.