No Arabic abstract
We compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way, that a DNN model, trained on meta-GGA (SCAN) data, reduce the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.
Using first principles simulations we have investigated the structural and bonding properties of dense fluid oxygen up to 180 GPa. We have found that band gap closure occurs in the molecular liquid, with a slow transition from a semi-conducting to a poor metallic state occurring over a wide pressure range. At approximately 80 GPa, molecular dissociation is observed in the metallic fluid. Spin fluctuations play a key role in determining the electronic structure of the low pressure fluid, while they are suppressed at high pressure.
An interatomic potential for Al-Tb alloy around the composition of Al90Tb10 was developed using the deep neural network (DNN) learning method. The atomic configurations and the corresponding total potential energies and forces on each atom obtained from ab initio molecular dynamics (AIMD) simulations are collected to train a DNN model to construct the interatomic potential for Al-Tb alloy. We show the obtained DNN model can well reproduce the energies and forces calculated by AIMD. Molecular dynamics (MD) simulations using the DNN interatomic potential also accurately describe the structural properties of Al90Tb10 liquid, such as the partial pair correlation functions (PPCFs) and the bond angle distributions, in comparison with the results from AIMD. Furthermore, the developed DNN interatomic potential predicts the formation energies of crystalline phases of Al-Tb system with the accuracy comparable to ab initio calculations. The structure factor of Al90Tb10 metallic glass obtained by MD simulation using the developed DNN interatomic potential is also in good agreement with the experimental X-ray diffraction data.
The marriage of density functional theory (DFT) and deep learning methods has the potential to revolutionize modern research of material science. Here we study the crucial problem of representing DFT Hamiltonian for crystalline materials of arbitrary configurations via deep neural network. A general framework is proposed to deal with the infinite dimensionality and covariance transformation of DFT Hamiltonian matrix in virtue of locality and use message passing neural network together with graph representation for deep learning. Our example study on graphene-based systems demonstrates that high accuracy ($sim$meV) and good transferability can be obtained for DFT Hamiltonian, ensuring accurate predictions of materials properties without DFT. The Deep Hamiltonian method provides a solution to the accuracy-efficiency dilemma of DFT and opens new opportunities to explore large-scale materials and physics.
The ab-initio theory of low-field electronic transport properties such as carrier mobility in semiconductors is well-established. However, an equivalent treatment of electronic fluctuations about a non-equilibrium steady state, which are readily probed experimentally, remains less explored. Here, we report a first-principles theory of electronic noise for warm electrons in semiconductors. In contrast with typical numerical methods used for electronic noise, no adjustable parameters are required in the present formalism, with the electronic band structure and scattering rates calculated from first-principles. We demonstrate the utility of our approach by applying it to GaAs and show that spectral features in AC transport properties and noise originate from the disparate time scales of momentum and energy relaxation, despite the dominance of optical phonon scattering. Our formalism enables a parameter-free approach to probe the microscopic transport processes that give rise to electronic noise in semiconductors.
Transport properties of dense liquid helium under the conditions of planets core and cool atmosphere of white dwarfs have been investigated by using the improved centroid path-integral simulations combined with density functional theory. The self-diffusion is largely higher and the shear viscosity is notably lower predicted with the quantum mechanical description of the nuclear motion compared with the description by Newton equation. The results show that nuclear quantum effects (NQEs), which depends on the temperature and density of the matter via the thermal de Broglie wavelength and the ionization of electrons, are essential for the transport properties of dense liquid helium at certain astrophysical conditions. The Stokes-Einstein relation between diffusion and viscosity in strongly coupled regime is also examined to display the influences of NQEs.