2025-03-28 | | Total: 9
With recent advancements in machine learning for interatomic potentials, Python has become the go-to programming language for exploring new ideas. While machine-learning potentials are often developed in Python-based frameworks, existing molecular dynamics software is predominantly written in lower-level languages. This disparity complicates the integration of machine learning potentials into these molecular dynamics libraries. Additionally, machine learning potentials typically focus on local features, often neglecting long-range electrostatics due to computational complexities. This is a key limitation as applications can require long-range electrostatics and even flexible charges to achieve the desired accuracy. Recent charge equilibration models can address these issues, but they require iterative solvers to assign relaxed flexible charges to the atoms. Conventional implementations also demand very tight convergence to achieve long-term stability, further increasing computational cost. In this work, we present a scalable Python implementation of a recently proposed shadow molecular dynamics scheme based on a charge equilibration model, which avoids the convergence problem while maintaining long-term energy stability and accuracy of observable properties. To deliver a functional and user-friendly Python-based library, we implemented an efficient neighbor list algorithm, Particle Mesh Ewald, and traditional Ewald summation techniques, leveraging the GPU-accelerated power of Triton and PyTorch. We integrated these approaches with the Python-based shadow molecular dynamics scheme, enabling fast charge equilibration for scalable machine learning potentials involving systems with hundreds of thousands of atoms.
We present the python package DiPolMol-Py, which can be used to calculate the rotational and hyperfine structure of 2Σ molecules. The calculations can be performed in the presence of dc magnetic fields, dc electric fields and far off-resonant optical fields. We additionally include functions to calculate the polarisability of the molecule and the transition dipole moment between different energy eigenstates. The package is applicable to many of the molecules which can be laser cooled, specifically the alkaline earth fluorides. We provide a constants file which includes many of the required literature values for CaF, SrF and BaF. Additional species can easily be added by updating this file.
We investigate the behavior of heavy impurities in edge plasma turbulence by analyzing their trajectories using the Hasegawa-Wakatani model. Through direct numerical simulations, we track ensembles of charged impurity particles over hundreds of eddy turnover times within statistically steady turbulent flows. Assuming that heavy impurities lag behind the flow, a novel derivation of relaxation time of heavy impurities is proposed. Our results reveal that heavy impurities can cluster within turbulence. We provide multiscale geometrical Lagrangian statistics of heavy impurities trajectories. To quantify directional changes, we analyze the scale-dependent curvature angle, along with the influence of the Stokes number on the mean curvature angles and the probability distribution function of curvature angles.
The need for optimized structures with good mechanical performance for the minimum weight is common in industry. Solid Isotropic Material with Penalization (SIMP) is a Topology Optimization (TO) method offering a trade-off between minimum compliance (i.e., maximum stiffness) and a fixed material amount for a given set of boundary conditions. Since TO is a non-convex problem, its gradient can be tuned by filtering the topology's contour, creating sharper material profiles without necessarily compromising optimality. However, despite simplifying the layout, some filters fail to address manufacturability concerns such as capillarity (thin tweaks as struts) generated by uncertain loading, vibration or fatigue. A tailored density-based filtering strategy is offered to tackle this issue. Additionally, volume fraction is left unconstrained so material can be strategically replenished through a logarithmic rule acting on the updated compliance. In doing so, an interpolation space with three degrees of freedom (volume, compliance, minimum thickness) is created, yielding diverse topologies for the same boundary conditions and design values along different stages of evolving topological families with distinct features. The optimization process is further accelerated by introducing the volume-compliance iterative scheme as a physical loss function in a Double Distance Neural Network (D2NN), obtaining similar results to 2,000 steps worth of vanilla iteration within 500 training epochs. This proposal offers a novel topology optimization design space based on minimum strut thickness - via filtering - and topological families defined by minimum volume fraction and compliance. The methodology is tested on several examples with diverse loading and boundary conditions, obtaining similarly satisfactory results, and then boosted via Machine Learning, acting as a fast and cheap surrogate.
Silicon (Si) is the primary donor dopant in gallium nitride (GaN), introduced through epitaxial growth or ion implantation. However, precise control over Si diffusion remains a critical challenge for high-performance device applications. This study investigates Si diffusion mechanisms in bulk GaN using first-principles density functional theory (DFT) calculations, supported by ultra-high-pressure annealing (UHPA) experiments. Vacancy-mediated diffusion pathways were analyzed using the SIESTA code, with minimum energy paths (MEPs) and activation barriers determined via the nudged elastic band (NEB) method. The results indicate that Si diffusion barriers vary with crystallographic direction, with the lowest barrier of 3.2 eV along [11-20] and the highest barrier of ~9.9 eV along [1-100], rendering diffusion in this direction highly improbable. Alternative diffusion mechanisms, including direct exchange and ring-like migration, exhibit prohibitively high barriers (>12 eV). Phonon calculations confirm that temperature-induced reductions in effective diffusion barriers are minimal. Experimental validation using SIMS analysis on Si-implanted GaN samples subjected to UHPA (1450{\deg}C, 1 GPa) confirms negligible Si diffusion under these extreme conditions. These findings resolve inconsistencies in prior reports and establish that Si-doped GaN remains highly stable, ensuring reliable doping profiles for advanced electronic and optoelectronic applications.
Crystal structure prediction (CSP) is crucial for identifying stable crystal structures in given systems and is a prerequisite for computational atomistic simulations. Recent advances in neural network potentials (NNPs) have reduced the computational cost of CSP. However, searching for stable crystal structures across the entire composition space in multicomponent systems remains a significant challenge. Here, we propose a novel genetic algorithm (GA) -based CSP method using a universal NNP. Our GA-based methods are designed to efficiently expand convex hull volumes while preserving the diversity of crystal structures. This approach draws inspiration from the similarity between convex hull updates and Pareto front evolution in multi-objective optimization. Our evaluation shows that the present method outperforms the symmetry-aware random structure generation, achieving a larger convex hull with fewer trials. We demonstrated that our approach, combined with the developed universal NNP (PFP), can accurately reproduce and explore phase diagrams obtained through DFT calculations; this indicates the validity of PFP across a wide range of crystal structures and element combinations. This study, which integrates a universal NNP with a GA-based CSP method, highlights the promise of these methods in materials discovery.
The reaction dynamics for the H + HeH+ → He + H+2 reaction in its electronic ground state is investigated using two different representations of the potential energy surface (PES). The first uses a combined kernel and neural network representation of UCCSD(T) reference data whereas the second is a corrected PES (cR-PES) that eliminates an artificial barrier in the entrance channel appearing in its initial expansion based on full configuration interaction reference data. Despite the differences between the two PESs, both yield kv=0,j=0≈2×10−9 cm3/molecule/s at T=10 K which is consistent with a T−independent Langevin rate kL=2.1×10−9 cm3/molecule/s but considerably larger than the only experimentally reported value kICR=(9.1±2.5)×10−10 cm3/molecule/s from ion cyclotron resonance experiments. Similarly, branching ratios for the reaction outcomes are comparable for the two PESs. However, when analysing less averaged properties such as initial state-selected T−dependent rate coefficients and final vibrational states of the H+2 product for low temperatures, the differences in the two PESs manifest themselves in the observables. Thus, depending on the property analyzed, accurate and globally valid representations of the PES are required, whereas more approximate and empirical construction schemes can be followed for state-averaged observables.
We propose an efficient novel path sampling-based framework designed to accelerate the investigation of rare events in complex molecular systems. A key innovation is the shift from sampling restricted path ensemble distributions, as in transition path sampling, to directly sampling the distribution of shooting points. This allows for a rejection-free algorithm that samples the entire path ensemble efficiently. Optimal sampling is achieved by applying a selection bias that is the inverse of the free energy along a reaction coordinate. The optimal reaction coordinate, the committor, is iteratively constructed as a neural network using AI for Molecular Mechanism Discovery (AIMMD), concurrently with the free energy profile, which is obtained through reweighting the sampled path ensembles. We showcase our algorithm on theoretical and molecular bechnmarks, and demonstrate how it provides at the same time molecular mechanism, free energy, and rates at a moderate computational cost.
We present the AGATE simulation code, a Python-based framework developed primarily for solving the magnetohydrodynamics (MHD) equations while maintaining adaptability to other equation sets. The code employs a modular, object-oriented architecture that separates interface specifications from numerical implementations, allowing users to customize numerical methods and physics models. Built on a Godunov-type finite-volume scheme, AGATE currently supports the ideal, Hall, and Chew-Goldberger-Low (CGL) MHD equations, with multiple acceleration options ranging from Numpy to GPU-enabled computation via NVIDIA CUDA. Performance testing demonstrates that our GPU implementations achieve 40-60x speedups over CPU versions. Comprehensive validation through established benchmarks confirms accurate reproduction of both linear and nonlinear phenomena across different MHD regimes. This combination of modularity, performance, and extensibility makes AGATE suitable for multiple applications: from rapid prototyping to production simulations, and from numerical algorithm development to physics education.