Abstract:We develop a method that can constrain any local exchange-correlation potential to preserve basic exact conditions. Using the method of Lagrange multipliers, we calculate for each set of given Kohn-Sham orbitals a constraint-preserving potential which is closest to the given exchange-correlation potential. The method is applicable to both the time-dependent (TD) and independent cases. The exact conditions that are enforced for the time-independent case are Galilean covariance, zero net force and torque, and Levy-Perdew virial theorem. For the time-dependent case we enforce translational covariance, zero net force, Levy-Perdew virial theorem, and energy balance. We test our method on the exchange (only) Krieger-Li-Iafrate (xKLI) approximate-optimized effective potential for both cases. For the time-independent case, we calculated the ground state properties of some hydrogen chains and small sodium clusters for some constrained xKLI potentials and Hartree-Fock (HF) exchange. The results (total energy, Kohn-Sham eigenvalues, polarizability, and hyperpolarizability) indicate that enforcing the exact conditions is not important for these cases. On the other hand, in the time-dependent case, constraining both energy balance and zero net force yields improved results relative to TDHF calculations. We explored the electron dynamics in small sodium clusters driven by cw laser pulses. For each laser pulse we compared calculations from TD constrained xKLI, TD partially constrained xKLI, and TDHF. We found that electron dynamics such as electron ionization and moment of inertia dynamics for the constrained xKLI are most similar to the TDHF results. Also, energy conservation is better by at least one order of magnitude with respect to the unconstrained xKLI. We also discuss the problems that arise in satisfying constraints in the TD case with a non-cw driving force.

Optimizing working parameters of the twin-range cutoff method in terms of accuracy and efficiency

Tuning the range separation parameter in periodic systems

Accelerating scientific codes by performance and accuracy modeling

A Posteriori Error Estimation and Adaptive Algorithm for the Atomistic/Continuum Coupling in 2D

Adaptive Multiscale Coupling Methods of Molecular Mechanics based on a Unified Framework of a Posteriori Error Estimates

The HI Content and Extent of Low Surface Brightness Galaxies - Could LSB Galaxies be Responsible for Damped Ly-alpha Absorption?

Accuracy Assessment of Atomistic Neural Network Potentials: The Impact of Cutoff Radius and Message Passing

Efficient and Precise Force Field Optimization for Biomolecules Using DPA-2

Towards chemical accuracy using a multi-mesh adaptive finite element method in all-electron density functional theory

Improving the precision of forces in real-space pseudopotential density functional theory

Parameter Optimization of a Dual-Comb Ranging System by Using a Numerical Simulation Method.

Adaptive Finite Element Approximations for Kohn-Sham Models

Truncation Effects of Shift Function Methods in Bulk Water Systems

Efficient Hartree-Fock Exchange Algorithm with Coulomb Range Separation and Long-Range Density Fitting

Assessing the accuracy of tailored coupled cluster methods corrected by electronic wave functions of polynomial cost

The optimal particle-mesh interpolation basis

Improving approximate-optimized effective potentials by imposing exact conditions: Theory and applications to electronic statics and dynamics

Database optimization for empirical interatomic potential models

A Posteriori Error Estimate and Adaptive Mesh Refinement Algorithm for Atomistic/Continuum Coupling with Finite Range Interactions in Two Dimensions

Embedding parameters in ab initio theory to develop approximations based on molecular similarity

Spin parameter optimization for spin-polarized extended tight-binding methods