Source author record

Huajie Chen

Huajie Chen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Stability of the Minimum Energy Path

The minimum energy path (MEP) is the most probable transition path that connects two equilibrium states of a potential energy landscape. It has been widely used to study transition mechanisms as well as transition rates in the fields of chemistry, physics, and materials science. In this paper, we derive a novel result establishing the stability of MEPs under perturbations of the energy landscape. The result also represents a crucial step towards studying the convergence of various numerical approximations of MEPs, such as the nudged elastic band and string methods.

preprint2021arXiv

Rigorous body-order approximations of an electronic structure potential energy landscape

We show that the local density of states (LDOS) of a wide class of tight-binding models has a weak body-order expansion. Specifically, we prove that the resulting body-order expansion for analytic observables such as the electron density or the energy has an exponential rate of convergence both at finite Fermi-temperature as well as for insulators at zero Fermi-temperature. We discuss potential consequences of this observation for modelling the potential energy landscape, as well as for solving the electronic structure problem.

preprint2020arXiv

A plane wave study on the localized-extended transitions in the one-dimensional incommensurate systems

Based on our recently proposed plane wave framework, we theoretically study the localized-extended transition in the one dimensional incommensurate systems with cosine type of potentials, which are in close connection to many recent experiments in the ultracold atom and photonic crystal. We formulate a propagator based scattering picture for the transition at the ground state and single particle mobility edge, in which the deeper connection between the incommensurate potentials, eigenstate compositions and transition mechanism is revealed. We further show that there exists a upper limit of localization length for all localized eigenstates, leading to an fundamental difference to the Anderson localization. Numerical calculations are presented alongside the analysis to justify our statements. The theoretical analysis and numerical methods can also be generalized to systems in higher dimensions, with different potentials or beyond the single particle regime, which would benefit the future studies in the related fields.

preprint2020arXiv

A Posteriori Error Estimates for Adaptive QM/MM Coupling Methods

Hybrid quantum/molecular mechanics models (QM/MM methods) are widely used in material and molecular simulations when MM models do not provide sufficient accuracy but pure QM models are computationally prohibitive. Adaptive QM/MM coupling methods feature on-the-fly classification of atoms during the simulation, allowing the QM and MM subsystems to be updated as needed. In this work, we propose such an adaptive QM/MM method for material defect simulations based on a new residual based it a posteriori error estimator, which provides both lower and upper bounds for the true error. We validate the analysis and illustrate the effectiveness of the new scheme on numerical simulations for material defects.

preprint2019arXiv

Locality of Interatomic Forces in Tight Binding Models for Insulators

The tight binding model is a minimalistic electronic structure model for predicting properties of materials and molecules. For insulators at zero Fermi-temperature we show that the potential energy surface of this model can be decomposed into exponentially localised site energy contributions, thus providing qualitatively sharp estimates on the interatomic interaction range which justifies a range of multi-scale models. For insulators at finite Fermi-temperature we obtain locality estimates that are uniform in the zero-temperature limit. A particular feature of all our results is that they depend only weakly on the point spectrum. Numerical tests confirm our analytical results. This work extends and strengthens (Chen, Ortner 2016) and (Chen, Lu, Ortner 2018) for finite temperature models.

preprint2015arXiv

Pair densities in density functional theory

The exact interaction energy of a many-electron system is determined by the electron pair density, which is not well-approximated in standard Kohn-Sham density functional models. Here we study the (complicated but well-defined) exact universal map from density to pair density. We survey how many common functionals, including the most basic version of the LDA (Dirac exchange with no correlation contribution), arise from particular approximations of this map. We develop an algorithm to compute the map numerically, and apply it to one-parameter families {a*rho(a*x)} of one-dimensional homogeneous and inhomogeneous single-particle densities. We observe that the pair density develops remarkable multiscale patterns which strongly depend on both the particle number and the "width" 1/a of the single-particle density. The simulation results are confirmed by rigorous asymptotic results in the limiting regimes a>>1 and a<<1. For one-dimensional homogeneous systems, we show that the whole spectrum of patterns is reproduced surprisingly well by a simple asymptotics-based ansatz which slowly smoothens out the "strictly correlated" a=0 pair density while slowly turning on the a=infty "exchange" terms as a increases. Our findings lend theoretical support to the celebrated semi-empirical idea [Becke93] to mix in a fractional amount of exchange, albeit not to assuming the mixing to be additive and taking the fraction to be a system independent constant.

preprint2015arXiv

QM/MM methods for crystalline defects. Part 1: Locality of the tight binding model

The tight binding model is a minimal electronic structure model for molecular modelling and simulation. We show that the total energy in this model can be decomposed into site energies, that is, into contributions from each atomic site whose influence on their environment decays exponentially. This result lays the foundation for a rigorous analysis of QM/MM coupling schemes.

preprint2015arXiv

QM/MM methods for crystalline defects. Part 2: Consistent energy and force-mixing

QM/MM hybrid methods employ accurate quantum (QM) models only in regions of interest (defects) and switch to computationally cheaper interatomic potential (MM) models to describe the crystalline bulk. We develop two QM/MM hybrid methods for crystalline defect simulations, an energy-based and a force-based formulation, employing a tight binding QM model. Both methods build on two principles: (i) locality of the QM model; and (ii) constructing the MM model as an explicit and controllable approximation of the QM model. This approach enables us to establish explicit convergence rates in terms of the size of QM region.

preprint2014arXiv

Numerical methods for a Kohn-Sham density functional model based on optimal transport

In this paper, we study numerical discretizations to solve density functional models in the "strictly correlated electrons" (SCE) framework. Unlike previous studies our work is not restricted to radially symmetric densities. In the SCE framework, the exchange-correlation functional encodes the effects of the strong correlation regime by minimizing the pairwise Coulomb repulsion, resulting in an optimal transport problem. We give a mathematical derivation of the self-consistent Kohn-Sham-SCE equations, construct an efficient numerical discretization for this type of problem for N = 2 electrons, and apply it to the H2 molecule in its dissociating limit. Moreover, we prove that the SCE density functional model is correct for the H2 molecule in its dissociating limit.

preprint2013arXiv

Adaptive Finite Element Approximations for Kohn-Sham Models

The Kohn-Sham equation is a powerful, widely used approach for computation of ground state electronic energies and densities in chemistry, materials science, biology, and nanosciences. In this paper, we study the adaptive finite element approximations for the Kohn-Sham model. Based on the residual type a posteriori error estimators proposed in this paper, we introduce an adaptive finite element algorithm with a quite general marking strategy and prove the convergence of the adaptive finite element approximations. Using D{\" o}rfler's marking strategy, we then get the convergence rate and quasi-optimal complexity. We also carry out several typical numerical experiments that not only support our theory,but also show the robustness and efficiency of the adaptive finite element computations in electronic structure calculations.

preprint2011arXiv

Numerical Analysis of Finite Dimensional Approximations of Kohn-Sham Models

In this paper, we study finite dimensional approximations of Kohn-Sham models, which are widely used in electronic structure calculations. We prove the convergence of the finite dimensional approximations and derive the a priori error estimates for ground state energies and solutions. We also provide numerical simulations for several molecular systems that support our theory.