Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set

Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

다중혈관 관상동맥 환자에서 y-문합을 이용하여 양쪽 내흉동맥만을 사용한 우회술의 조기 성적

QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

/pdf/quantum-espresso-a-modular-and-open-source-software-project-4t6vlupbeh.pdf

QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials

We present ab initio quantum-mechanical molecular-dynamics simulations of the liquid-metal--amorphous-semiconductor transition in Ge. Our simulations are based on (a) finite-temperature density-functional theory of the one-electron states, (b) exact energy minimization and hence calculation of the exact Hellmann-Feynman forces after each molecular-dynamics step using preconditioned conjugate-gradient techniques, (c) accurate nonlocal pseudopotentials, and (d) Nos\'e dynamics for generating a canonical ensemble. This method gives perfect control of the adiabaticity of the electron-ion ensemble and allows us to perform simulations over more than 30 ps. The computer-generated ensemble describes the structural, dynamic, and electronic properties of liquid and amorphous Ge in very good agreement with experiment. The simulation allows us to study in detail the changes in the structure-property relationship through the metal-semiconductor transition. We report a detailed analysis of the local structural properties and their changes induced by an annealing process. The geometrical, bonding, and spectral properties of defects in the disordered tetrahedral network are investigated and compared with experiment.

Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium.

We study adatom-covered single layers of ${\mathrm{CrSiTe}}_{3}$ and ${\mathrm{CrGeTe}}_{3}$ using first-principles calculations based on hybrid functionals. We find that the insulating ground state of a monolayer of La (Lu) deposited on single-layer ${\mathrm{CrSiTe}}_{3}$ (${\mathrm{CrGeTe}}_{3}$) carries spontaneously generated current loops around the Cr sites. These ``flux states'' induce antiferromagnetically ordered orbital moments on the Cr sites and are also associated with nontrivial topological properties. The calculated Chern numbers for these systems are predicted to be $\ifmmode\pm\else\textpm\fi{}1$ even in the absence of spin-orbit coupling, with sizable gaps on the order of 100 meV. The flux states and the associated topological phases result from spontaneous time-reversal symmetry breaking due to the presence of nonlocal Coulomb interactions.

Flux States and Topological Phases from Spontaneous Time-Reversal Symmetry Breaking in CrSi ( Ge ) Te 3 -Based Systems

We develop an exact formalism for performing first-principles calculations for insulators at fixed electric polarization. As shown by Sai, Rabe, and Vanderbilt (SRV) [Phys. Rev. B 66, 104108 (2002)], who designed an approximate method to tackle the same problem, such an approach allows one to map out the energy landscape as a function of polarization, providing a powerful tool for the theoretical investigation of polar materials. We apply our method to a system in which the ionic contribution to the polarization dominates (a broken-inversion-symmetry perovskite), one in which this is not the case (a III-V semiconductor), and one in which an additional degree of freedom plays an important role (a ferroelectric phase of KNO3). We find that while the SRV method gives rather accurate results in the first case, the present approach provides important improvements to the physical description in the latter cases.

/pdf/first-principles-calculations-for-insulators-at-constant-pykuttrfll.pdf

First-principles calculations for insulators at constant polarization.

We point out that, contrary to an assertion by Chen et al. [Phys. Rev. B 88, 125110 (2013)], the nonquantized part of the intrinsic anomalous Hall conductivity can indeed be expressed as a Fermi-surface property even when Weyl points are present in the band structure.

/pdf/comment-on-weyl-fermions-and-the-anomalous-hall-effect-in-1k0qyl6i6i.pdf

Comment on "weyl fermions and the anomalous Hall effect in metallic ferromagnets"

Total energies in Se. I. The trigonal crystal

noncollinear spins induces much larger magnetic charges than those that depend on spin-orbit coupling. The hexagonal manganites RMnO3 and ferrites RFeO3 (R = Sc, Y, In, Ho-Lu) exhibit strong couplings between electric, magnetic and structural degrees of freedom, with the transition-metal ions in the basal plane antiferromagnetically coupled through super-exchange so as to form a 120 noncollinear spin arrangement. Here we present a theoretical study of the magnetic charges, and of the spin-lattice and spin-electronic ME constants, in these hexagonal manganites and ferrites, clarifying the conditions under which exchange striction leads to an enhancedZ m values and anomalously large in-plane spin-lattice ME eects.

/pdf/magnetic-charges-and-magnetoelectricity-in-hexagonal-rare-4e7wpdww8z.pdf

David Vanderbilt

Papers

Flux States and Topological Phases from Spontaneous Time-Reversal Symmetry Breaking in CrSi ( Ge ) Te 3 -Based Systems

First-principles calculations for insulators at constant polarization.

Comment on "weyl fermions and the anomalous Hall effect in metallic ferromagnets"

Total energies in Se. I. The trigonal crystal

Magnetic charges and magnetoelectricity in hexagonal rare-earth manganites and ferrites