Related papers: A Simple and Efficient Lattice Summation Method fo…
Accounting for nuclear quantum effects (NQEs) can significantly alter material properties at finite temperatures. Atomic modeling using the path-integral molecular dynamics (PIMD) method can fully account for such effects, but requires…
A novel single MRT-featured lattice Boltzmann method (SmrtLBM) is introduced. The intricate algorithm of the multiple-relaxation-time (MRT) collision operator is reimagined and distilled into a streamlined, single-relaxation-time-like…
The application of molecular dynamics (MD) simulations to the interpretation of Raman scattering spectra is hindered by inability of atomistic simulations to account for the dynamic evolution of electronic polarizability, requiring the use…
We present a measurement-cost efficient implementation of Strongly-Contracted $N$-Electron Valence Perturbation Theory (SC-NEVPT2) for use on near-term quantum devices. At the heart of our algorithm we exploit the properties of adaptive…
Coulomb interaction, following an inverse-square force-law, quantifies the amount of force between two stationary and electrically charged particles. The long-range nature of Coulomb interactions poses a major challenge to molecular…
In many applications, such as plasma edge simulation of a nuclear fusion reactor, a coupled PDE/kinetic description is required, which is usually solved with a coupled finite-volume/Monte-Carlo method. Different procedures have been…
We introduce a Monte-Carlo algorithm for the simulation of charged particles moving in the continuum. Electrostatic interactions are not instantaneous as in conventional approaches, but are mediated by a constrained, diffusing electric…
We explore a novel simulation route for Plasma Wakefield Acceleration (PWFA) by using the computational method known as the Lattice Boltzmann Method (LBM). LBM is based on a discretization of the continuum kinetic theory while assuring the…
The development of next-generation molecular simulation models requires moving beyond pre-defined functional forms toward machine learning (ML) techniques that directly capture multiscale physics. Here, we demonstrate such an approach using…
We derive equations for the strongly coupled system of light and dense atomic ensembles. The formalism includes an arbitrary internal level structure for the atoms and is not restricted to weak excitation of atoms by light. In the low light…
Computational Fluid Dynamics (CFD) simulations are often constrained by the memory-bound nature of sparse matrix-vector operations, which eventually limits performance on modern high-performance computing (HPC) systems. This work introduces…
Quantum simulation holds the promise of improving the atomic simulations used at EDF to anticipate the ageing of materials of interest. One simulator in particular seems well suited to modeling interacting electrons: the Rydberg atoms…
We have extended the multilevel summation (MLS) method, originally developed to evaluate long-range Coulombic interactions in molecular dynamics (MD) simulations [Skeel et al., J. Comput. Chem., 23, 673 (2002)], to handle dispersion…
Most current machine learning interatomic potentials (MLIPs) rely on short-range approximations, without explicit treatment of long-range electrostatics. To address this, we recently developed the Latent Ewald Summation (LES) method, which…
The relativistic coupled-cluster (RCC) theory has been applied recently to a number of heavy molecules to determine their properties very accurately. Since it demands large computational resources, the method is often approximated to…
Memory-bound algorithms show complex performance and energy consumption behavior on multicore processors. We choose the lattice-Boltzmann method (LBM) on an Intel Sandy Bridge cluster as a prototype scenario to investigate if and how…
We develop an algorithm suitable for parallel molecular dynamics simulations in $d$ spatial dimensions and describe its implementation in C++. All routines work in arbitrary $d$; the maximum simulated $d$ is limited only by available…
The emergence of new nanoporous materials, based e.g. on 2D materials, offers new avenues for water filtration and energy. There is accordingly a need to investigate the molecular mechanisms at the root of the advanced performances of these…
For nearly the past 30 years, Centroid Molecular Dynamics (CMD) has proven to be a viable classical-like phase space formulation for the calculation of quantum dynamical properties. However, calculation of the centroid effective force…
Combined extracorporeal membrane oxygenation (ECMO) and continuous renal replacement therapy (CRRT) pose complex hemodynamic challenges in intensive care. In this study, a comprehensive lumped parameter model (LPM) is developed to simulate…