No Arabic abstract
We give a detailed presentation of the theory and numerical implementation of an expression for the adiabatic energy flux in extended systems, derived from density-functional theory. This expression can be used to estimate the heat conductivity from equilibrium ab initio molecular dynamics, using the Green-Kubo linear response theory of transport coefficients. Our expression is implemented in an open-source component of the QE suite of computer codes for quantum mechanical materials modelling, which is being made publicly available.
We present a general and convenient first principle method to study near-field radiative heat transfer. We show that the Landauer-like expression of heat flux can be expressed in terms of a frequency and wave-vector dependent macroscopic dielectric function which can be obtained from the linear response density functional theory. A random phase approximation is used to calculate the response function. We computed the heat transfer in three systems -- graphene, molybdenum disulfide (MoS$_2$), and hexagonal boron nitride (h-BN). Our results show that the near-field heat flux exceeds the blackbody limit up to four orders of magnitude. With the increase of the distances between two parallel sheets, a $1/d^2$ dependence of heat flux is shown, consistent with Coulombs law. The heat transfer capacity is sensitive to the dielectric properties of materials. Influences from chemical potential and temperature are also discussed. Our method can be applied to a wide range of materials including systems with inhomogeneities.
Thermal and other transport coefficients were recently shown to be largely independent of the microscopic representation of the energy (current) densities or, more generally, of the relevant conserved densities/currents. In this paper we show how this gauge invariance, which is intimately related to the intrinsic indeterminacy of the energy of isolated atoms, can be exploited to optimize the statistical properties of the current time series from which the transport coefficients can be evaluated. To this end, we make use of a variational principle that relies on the metric properties of the conserved currents, treated as elements of an abstract linear space. Different metrics would result in different variational principles. We finally show how a recently proposed data-analysis methodology based on the theory of transport in multi-component systems can be recovered by a suitable choice of this metric.
Structural, electronic, vibrational and dielectric properties of LaBGeO$_5$ with the stillwellite structure are determined based on textit{ab initio} density functional theory. The theoretically relaxed structure is found to agree well with the existing experimental data with a deviation of less than $0.2%$. Both the density of states and the electronic band structure are calculated, showing five distinct groups of valence bands. Furthermore, the Born effective charge, the dielectric permittivity tensors, and the vibrational frequencies at the center of the Brillouin zone are all obtained. Compared to existing model calculations, the vibrational frequencies are found in much better agreement with the published experimental infrared and Raman data, with absolute and relative rms values of 6.04 cm$^{-1}$, and $1.81%$, respectively. Consequently, numerical values for both the parallel and perpendicular components of the permittivity tensor are established as 3.55 and 3.71 (10.34 and 12.28), respectively, for the high-(low-)frequency limit.
In this work we elaborate on two recently discovered invariance principles, according to which transport coefficients are, to a large extent, independent of the microscopic definition of the densities and currents of the conserved quantities being transported (energy, momentum, mass, charge). The first such principle, gauge invariance, allows one to define a quantum adiabatic energy current from density-functional theory, from which the heat conductivity can be uniquely defined and computed using equilibrium ab initio molecular dynamics. When combined with a novel topological definition of atomic oxidation states, gauge invariance also sheds new light onto the mechanisms of charge transport in ionic conductors. The second principle, convective invariance, allows one to extend the analysis to multi-component systems. These invariance principles can be combined with new spectral analysis methods for the current time series to be fed into the Green-Kubo formula to obtain accurate estimates of transport coefficients from relatively short molecular dynamics simulations.
Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Classical MD typically employ predefined functional forms which limits their ability to capture such complex size-dependent structural and dynamical transformation. Neural Network (NN) based potentials represent flexible alternatives and in principle, well-trained NN potentials can provide high level of flexibility, transferability and accuracy on-par with the reference model used for training. A major challenge, however, is that NN models are interpolative and requires large quantities of training data to ensure that the model adequately samples the energy landscape both near and far-from-equilibrium. Here, we introduce an active learning (AL) scheme that trains a NN model on-the-fly with minimal amount of first-principles based training data. Our AL workflow is initiated with a sparse training dataset (1 to 5 data points) and is updated on-the-fly via a Nested Ensemble Monte Carlo scheme that iteratively queries the energy landscape in regions of failure and updates the training pool to improve the network performance. Using a representative system of gold clusters, we demonstrate that our AL workflow can train a NN with ~500 total reference calculations. Our NN predictions are within 30 meV/atom and 40 meV/AA of the reference DFT calculations. Moreover, our AL-NN model also adequately captures the various size-dependent structural and dynamical properties of gold clusters in excellent agreement with DFT calculations and available experiments.