What problem does this paper attempt to address?

This paper attempts to solve the accuracy problems encountered in calculating the high - temperature Coulomb density matrix under periodic boundary conditions. Specifically, the author proposes a new method to handle the Coulomb interaction in systems with periodic boundary conditions to ensure that all interactions are processed at the same level of accuracy. This helps to improve the accuracy of the solution of the Bloch equation in path - integral calculations. ### Main Problems and Solutions 1. **Background and Challenges**: - In many quantum Monte Carlo (QMC) simulations in the fields of physics and chemistry, especially in the path - integral Monte Carlo (PIMC) method, it is necessary to accurately handle the Coulomb interaction. - Most previous methods rely on decomposing the Ewald potential into short - range and long - range parts. This decomposition introduces additional approximations, leading to a decrease in simulation efficiency and accuracy. - For the Coulomb propagator under periodic boundary conditions, there is no generally accepted calculation method, resulting in unnecessary approximations, which in turn affect the efficiency and accuracy of many - body simulations. 2. **Proposed Method**: - The author proposes a method to directly calculate the Coulomb pair - density matrix, avoiding the decomposition of the potential into short - range and long - range parts. - The Coulomb pair - density matrix is calculated by summing over eigenstates, avoiding grid errors and being able to efficiently provide off - diagonal elements. - A new Ewald action is introduced, combining the original approximation and quantum correction terms, ensuring convergence to the original approximation at large separation distances while maintaining the accuracy of short - range behavior. 3. **Specific Steps**: - **Calculating the Density Matrix of Isolated Particle Pairs**: Use the eigenstate - summing method to calculate the density matrix of isolated particle pairs, avoiding grid errors. - **Handling Periodic Image Charges**: Apply the pair - approximation method to handle periodic image charges, assuming that path diffusion can be independently averaged, reducing approximation errors. - **Introducing Background Terms**: By introducing the background term \(u_{BG}\), correct the errors caused by neutralizing background charges, ensuring the consistency of diagonal and off - diagonal elements. 4. **Verification and Application**: - Verify the accuracy of the proposed method by testing the potential and kinetic energies of isolated hydrogen atoms and molecules. - In periodic systems, evaluate the accuracy of the high - temperature density matrix by calculating the residuals of the Bloch equation. - Apply it to PIMC simulations of hydrogen, helium and their mixtures, and study the quantum effects in one - component plasmas, demonstrating the effectiveness and superiority of this method. ### Conclusion This method has higher accuracy and fewer numerical errors compared to previous decomposition methods, and is especially suitable for path - integral calculations of high - temperature many - body systems. Through quantitative test results, the superior performance of this method in isolated and periodic systems is proven, providing a benchmark for future density - matrix construction. ### Formula Summary - High - temperature many - body density matrix expression: \[ \rho_{B/F}(R, R'; \beta)=\frac{1}{N!}\sum_P(\pm 1)^P\int_{R\rightarrow PR'}dR_i\rho(R, R_1; \tau)\rho(R_1, R_2; \tau)\cdots\rho(R_{M - 1}, PR'; \tau) \] - Free - particle density matrix: \[ \rho_0(r, r'; \tau)=(4\pi\lambda\tau)^{-D/2}\exp\left(-\frac{(r - r')^2}{4\lambda\tau}\right) \] - Ewald action: \[ u_{EW}(r, r'; \tau)\approx\frac{\tau}{2}Q_iQ_j[V_{EW}(r)+V_{EW}(r')]+\sum_L\Delta u(r + L, r' + L; \tau)+u_{BG} \] These formulas and methods provide more accurate and efficient solutions for path - integral calculations of high - temperature many - body systems.

Computation of the High Temperature Coulomb Density Matrix in Periodic Boundary Conditions

Force and Heat Current Formulas for Many-Body Potentials in Molecular Dynamics Simulations with Applications to Thermal Conductivity Calculations

The high-density limit of two-electron systems: Results from the extended Overhauser approach

Extended Hartree-Fock method based on pair density functional theory

Stochastic density functional theory combined with Langevin dynamics for warm dense matter

Bypassing the malfunction junction in warm dense matter simulations

Density functional calculation of many-electron systems in cartesian coordinate grid

Density fitting in periodic systems: application to TDHF in diamond and oxides

A note on accurate pressure calculations of Coulomb systems with periodic boundary conditions

A new tunnel technique with acellular dermal matrix for soft tissue preparation prior to symphyseal block graft--a description of technique and case report.

Bound state breaking and the importance of thermal exchange-correlation effects in warm dense hydrogen

The BCS Critical Temperature at High Density

A well-tempered density functional theory of electrons in molecules

Multicomponent Density-Functional Theory for Electrons and Nuclei

Linear scaling computation of the Fock matrix VII. Periodic Density Functional Theory at the $Γ$-point

Density functional theory in one-dimension for contact-interacting fermions

Wave-function Inspired Density Functional Applied to the H_2/H_2^+ Challenge

Interaction Picture Density Matrix Quantum Monte Carlo

Self-consistent calculation of the coupling constant in the Gross-Pitaevskii equation