Related papers: Minimal Basis Iterative Stockholder: Atoms in Mole…
The rigorous stability analysis of high-order implicit-explicit multistep (IEMS) methods for nonlinear parabolic equations by using discrete energy arguments is a long standing open issue due to their non-A-stable property. A novel…
A parameterization strategy for molecular models on the basis of force fields is proposed, which allows a rapid development of models for small molecules by using results from quantum mechanical (QM) ab initio calculations and thermodynamic…
The event-chain Monte Carlo (ECMC) method is an irreversible Markov process based on the factorized Metropolis filter and the concept of lifted Markov chains. Here, ECMC is applied to all-atom models of multi-particle interactions that…
Long ranged electrostatic interactions are time consuming to calculate in molecular dynamics and Monte-Carlo simulations. We introduce an algorithmic framework for simulating charged particles which modifies the dynamics so as to allow…
The quantum theory of atoms in molecules, QTAIM, is employed to identify AIM and quantify their interactions through the partitioning of molecule into atomic basins in the real space and it is confined only to the purely electronic systems…
The electric microfield distribution at charged particles is studied for two-component electron-ion plasmas using molecular dynamics simulation and theoretical models. The particles are treated within classical statistical mechanics using…
Electronic structure methods for accurate calculation of molecular properties have a high cost that grows steeply with the problem size, therefore, it is helpful to have the underlying atomic basis functions that are less in number but of…
Fundamental understanding of interatomic forces in molecules must emerge from quantum mechanics, yet widely used empirical force fields rely on simplified mechanistic approximations that often fail to capture the complexity of many-body…
To accurately model the electron density and polarization, a polarizable multipole (PM) model using the AMOEBA force field has been introduced \cite{Ren:2003, Shi:2013} recently. In the AMOEBA force field, the traditional point atomic…
We propose a simplified version of local molecular field (LMF) theory to treat Coulomb interactions in simulations of ionic fluids. LMF theory relies on splitting the Coulomb potential into a short-ranged part that combines with other…
We present an alternative, univocal characterization of the continuous transition from atomic to molecular shape in the Coulomb system constituted by two identical particles and a third particle with the opposite charge, as the mass ratio…
The starting point is the problem of finding the interaction energy of two coinciding homogeneous cubic charge distributions. The brute force method of subdividing the cube into $N^3$ sub-cubes and doing the sums results in slow convergence…
A double-atom partitioning of the molecular one-electron density matrix is used to describe atoms and bonds. All calculations are performed in Hilbert space. The concept of atomic weight functions (familiar from Hirshfeld analysis of the…
The recently developed energy conserving semi-implicit method (ECsim) for PIC simulation is applied to multiple scale problems where the electron-scale physics needs to be only partially retained and the interest is on the macroscopic or…
The multi-component quantum theory of atoms in molecules (MC-QTAIM) analysis is done on methane, ethylene, acetylene and benzene as selected basic hydrocarbons. This is the first report on applying the MC-QTAIM analysis on polyatomic…
We push the boundaries of electronic structure-based \textit{ab-initio} molecular dynamics (AIMD) beyond 100 million atoms. This scale is otherwise barely reachable with classical force-field methods or novel neural network and machine…
Electrochemical Impedance Spectroscopy (EIS) and Equivalent Circuit Models (ECMs) are widely used to characterize the impedance and estimate parameters of electrochemical systems such as batteries. We use a generic ECM with ten parameters…
Molecular dynamics (MD) simulations allow atomistic insights into chemical and biological processes. Accurate MD simulations require computationally demanding quantum-mechanical calculations, being practically limited to short timescales…
The evaluation of Coulomb forces is a difficult task. The summations that are involved converge only conditionally and care has to be taken in selecting the appropriate procedure to define the limits. The Ewald method is a standard method…
We propose a device for studying the Fermi-Hubbard model with long-range Coulomb interactions using an array of quantum dots defined in a semiconductor two-dimensional electron gas system. Bands with energies above the lowest energy band…