Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStochastic Constrained Extended System Dynamics for Solving Charge Equilibration Models
50
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a new stochastic extended Lagrangian solution to charge equilibration that eliminates self-consistent field (SCF) calculations, eliminating the computational bottleneck in solving the many-body solution with standard SCF solvers. By formulating both charges and chemical potential as latent variables, and introducing a holonomic constraint that satisfies charge conservation, the SC-XLMD method accurately reproduces structural, thermodynamic, and dynamics properties using ReaxFF, and shows excellent weak- and strong-scaling performance in the LAMMPS molecular simulation package.
rate research
Read More
We extend the Eckart theorem, from the ground state to excited statew, which introduces an energy augmentation to the variation criterion for excited states. It is shown that the energy of a very good excited state trial function can be slightly lower than the exact eigenvalue. Further, the energy calculated by the trial excited state wave function, which is the closest to the exact eigenstate through Gram-Schmidt orthonormalization to a ground state approximant, is lower than the exact eigenvalue as well. In order to avoid the variation restrictions inherent in the upper bound variation theory based on Hylleraas, Undheim, and McDonald [HUM] and Eckart Theorem, we have proposed a new variation functional Omega-n and proved that it has a local minimum at the eigenstates, which allows approaching the eigenstate unlimitedly by variation of the trial wave function. As an example, we calculated the energy and the radial expectation values of Triplet-S(even) Helium atom by the new variation functional, and by HUM and Eckart theorem, respectively, for comparison. Our preliminary numerical results reveal that the energy of the calculated excited states 3rd Triplet-S(even) and 4th Triplet-S(even) may be slightly lower than the exact eigenvalue (inaccessible by HUM theory) according to the General Eckart Theorem proved here, while the approximate wave function is better than HUM.
We demonstrate an efficient algorithm for inverse problems in time-dependent quantum dynamics based on feedback loops between Hamiltonian parameters and the solutions of the Schr{o}dinger equation. Our approach formulates the inverse problem as a target vector estimation problem and uses Bayesian surrogate models of the Schr{o}dinger equation solutions to direct the optimization of feedback loops. For the surrogate models, we use Gaussian processes with vector outputs and composite kernels built by an iterative algorithm with Bayesian information criterion (BIC) as a kernel selection metric. The outputs of the Gaussian processes are designed to model an observable simultaneously at different time instances. We show that the use of Gaussian processes with vector outputs and the BIC-directed kernel construction reduce the number of iterations in the feedback loops by, at least, a factor of 3. We also demonstrate an application of Bayesian optimization for inverse problems with noisy data. To demonstrate the algorithm, we consider the orientation and alignment of polyatomic molecules SO$_2$ and chiral propylene oxide (PPO) induced by strong laser pulses. We use simulated time evolutions of the orientation or alignment signals to determine the relevant components of the molecular polarizability tensors to within 1% accuracy. We show that, for the five independent components of the polarizability tensor of PPO, this can be achieved with as few as 30 quantum dynamics calculations.
We propose a fast method for the calculation of short-range interactions in molecular dynamics simulations. The so-called random-batch list method is a stochastic version of the classical neighbor-list method to avoid the construction of a full Verlet list, which introduces two-level neighbor lists for each particle such that the neighboring particles are located in core and shell regions, respectively. Direct interactions are performed in the core region. For the shell zone, we employ a random batch of interacting particles to reduce the number of interaction pairs. The error estimate of the algorithm is provided. We investigate the Lennard-Jones fluid by molecular dynamics simulations, and show that this novel method can significantly accelerate the simulations with a factor of several fold without loss of the accuracy. This method is simple to implement, can be well combined with any linked cell methods to further speed up and scale up the simulation systems, and can be straightforwardly extended to other interactions such as Ewald short-range part, and thus it is promising for large-scale molecular dynamics simulations.
The weighted ensemble method, introduced by Huber and Kim, [G. A. Huber and S. Kim, Biophys. J. 70, 97 (1996)], is one of a handful of rigorous approaches to path sampling of rare events. Expanding earlier discussions, we show that the technique is statistically exact for a wide class of Markovian and non-Markovian dynamics. The derivation is based on standard path-integral (path probability) ideas, but recasts the weighted-ensemble approach as simple resampling in path space. Similar reasoning indicates that arbitrary nonstatic binning procedures, which merely guide the resampling process, are also valid. Numerical examples confirm the claims, including the use of bins which can adaptively find the target state in a simple model.
We present an implementation of range separated functionals utilizing the Slater-function on grids in real space in the projector augmented waves method. The screened Poisson equation is solved to evaluate the necessary screened exchange integrals on Cartesian grids. The implementation is verified against existing literature and applied to the description of charge transfer excitations. We find very slow convergence for calculations within linear response time-dependent density functional theory and unoccupied orbitals of the canonical Fock operator. Convergence can be severely improved by using Huzinagas virtual orbitals instead. This combination furthermore enables an accurate determination of long-range charge transfer excitations by means of ground-state calculations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
