SciPost Submission Page
The selfconsistent quantumelectrostatic problem in strongly nonlinear regime
by P. Armagnat, A. Lacerda, B. Rossignol, C. Groth, X. Waintal
This is not the current version.
Submission summary
As Contributors:  Xavier Waintal 
Arxiv Link:  https://arxiv.org/abs/1905.01271v1 
Date submitted:  20190507 
Submitted by:  Waintal, Xavier 
Submitted to:  SciPost Physics 
Domain(s):  Computational 
Subject area:  Condensed Matter Physics  Computational 
Abstract
The selfconsistent quantumelectrostatic (also known as PoissonSchr\"odinger) problem is notoriously difficult in situations where the density of states varies rapidly with energy. At low temperatures, these fluctuations make the problem highly nonlinear which renders iterative schemes deeply unstable. We present a stable algorithm that provides a solution to this problem with controlled accuracy. The technique is intrinsically convergent including in highly nonlinear regimes. We illustrate our approach with (i) a calculation of the compressible and incompressible stripes in the integer quantum Hall regime and (ii) a calculation of the differential conductance of a quantum point contact geometry. Our technique provides a viable route for the predictive modeling of the transport properties of quantum nanoelectronics devices.
Current status:
Ontology / Topics
See full Ontology or Topics database.Submission & Refereeing History
Reports on this Submission
Anonymous Report 2 on 201965 Invited Report
 Cite as: Anonymous, Report on arXiv:1905.01271v1, delivered 20190605, doi: 10.21468/SciPost.Report.997
Report
In this manuscript, the authors present the computational scheme to perform SchrodingerPoisson simulations of nanoscale quantum systems. The main result of the manuscript is eliminating the convergence issues of the SchrodingerPoisson problem and reducing the necessity to perform computationally intensive steps of quantum simulation to only a few iterations.
The framework consists of several steps:
 Inverting Poisson equation as a capacitance equation for the density as a function of potential, thus rewriting the original problem to that of finding the intersection between solutions of inverted Poisson and Schrodinger equations. The authors show that the capacitance is a better behaved numerical quantity than its inverse.
 Developing a solver for inverted Poisson equation to obtain the capacitance and density for the given potential.
 Performing adiabatic approximations for Poisson and Schrodinger equations, allowing for linearization of electronic density and reducing the intersection problem to a set of local problems that can be converged rapidly.
 Performing a set of steps to get the adiabatic problem to converge to the full selfconsistent solution.
Additionally, the authors implement the finite volume method to discretize the Poisson equation, which allows them to enforce the charge conservation.
SchrodingerPoisson is a nonlinear problem and thus has convergence issues. In fact, as also described in this paper, the speed of convergence decreases with the system size  the potential needs to be converged on all lattice sites. This paper eliminates this unwanted prefactor by mapping it to multiple 1d selfconsistent equations. While the principal complexity of SchrodingerPoisson method is in the cubic scaling of the calculation of density from the Schrodinger equation, this manuscript eliminates the additional overhead, significantly helping to achieve a truly realistic computational modeling. I am sure this paper is useful as a handbook for various implementations of the SchrodingerPoisson method and thus should definitely be published in SciPost.
The paper is very clearly written, all the steps are explained in detail.
Additionally to the typos noted by the other referee, several edits might improve the manuscript:
1. Since the main contribution of the manuscript is the algorithm, it would be worth spelling it out completely with all the steps in a subsection.
2. Please provide more details on the parametrization of the problem. Section 2 says "We assume … a 2DEG density $2.11 x 10^{11} cm^{2]$". Since density is selfconsistently evaluated and changes in the presence of gates, what is that density the authors refer to? Is it the density of dopants or is the dopant density fitted such as to have this density in 2DEG in the absence of external gates? I believe that both this and the companion paper would also benefit from writing down all the parameters and boundary conditions.
Anonymous Report 1 on 201958 Invited Report
 Cite as: Anonymous, Report on arXiv:1905.01271v1, delivered 20190508, doi: 10.21468/SciPost.Report.939
Strengths
1. Clearly written
2. Well presented and documented
3. Detailed discussion of a difficult problem
4. Valuable information for researchers discussed in one place
5. New stable approach to a known problem
Weaknesses
No obvious weakness
Report
The authors propose an effective method to solve "the
the selfconsistent quantumelectrostatic problem" at
zero temperature. Commonly, this has been avoided by
performing the calculations at a finite, but low
temperature. Indeed, the suggested methodology is
well explained, documented, and supported by examples.
All this together makes the manuscript very valuable
to many researchers that have to address this problem
in some form in their modeling of physical phenomena.
Requested changes
For the authors to consider:
1.
The manuscript is well written, but I am note quite
sure if in the sentence:
"The technique is intrinsically convergent including
in highly nonlinear regimes."
the word "including" is the best choice.
2.
The authors do investigate the "selfconsistent
quantumelectrostatic problem", which they also
refer to as the "PoissonSchrödinger problem".
In view of a difference in use of the terms within
the physics and the mathematical communities I would
like the authors briefly to relate this problem to
the Hartree Approximation (HA) which in the physics
and the quantum chemistry communities is also the
mentioned simultaneous solution of the Schrödinger
and the Poisson equations with the condition that
the wavefunctions are to be orthogonal.
In the mathematics community the problem is often
confronted without this last condition.
The different stand point can be referred back to
the need of physicists to relate the HA to higher
order ones, the HartreeFock Approximation (HFA), or
higher order Green functions schemes. This issue
is briefly mentioned in the Introduction of:
The European Physical Journal B 84, 699 (2011).
We thank our referees for the report. In response to his/her comments,

we now mention the selfconsistent Hartree approximation as an equivalent name in the introduction.

We have replaced the word "including" by "even" in the abstract.
(in reply to Report 2 on 20190605)
We thank our referee for the report. In response to his/her comment,
We have added a new section 8 "Summary of the Algorithm". This section contains a chart that explains how the different parts of the algorithm are articulated.
We have added a few more details on the parameters used in the simulations.