Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGauge optimization of time series for thermal-transport simulations
130
0
0.0
(
0
)
Added by
Aris Marcolongo Dr.
Publication date
2019
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
Applications of commodity polymers are often hindered by their low thermal conductivity. In these systems, going from the standard polymers dictated by weak van der Waals interactions to biocompatible hydrogen bonded smart polymers, the thermal transport coefficient k varies between 0.1 - 0.4 W/Km. Combining all-atom molecular dynamics simulations with some experiments, we study thermal transport and its link to the elastic response of commodity plastics. We find that there exists a maximum attainable stiffness (or sound wave velocity), thus providing an upper bound of k for these solid polymers. The specific chemical structure and the glass transition temperature play no role in controlling k, especially when the microscopic interactions are hydrogen bonding based. Our results are consistent with the minimum thermal conductivity model and existing experiments. The effect of polymer stretching on k is also discussed.
Simulations of systems with quenched disorder are extremely demanding, suffering from the combined effect of slow relaxation and the need of performing the disorder average. As a consequence, new algorithms, improved implementations, and alternative and even purpose-built hardware are often instrumental for conducting meaningful studies of such systems. The ensuing demands regarding hardware availability and code complexity are substantial and sometimes prohibitive. We demonstrate how with a moderate coding effort leaving the overall structure of the simulation code unaltered as compared to a CPU implementation, very significant speed-ups can be achieved from a parallel code on GPU by mainly exploiting the trivial parallelism of the disorder samples and the near-trivial parallelism of the parallel tempering replicas. A combination of this massively parallel implementation with a careful choice of the temperature protocol for parallel tempering as well as efficient cluster updates allows us to equilibrate comparatively large systems with moderate computational resources.
We implement the statistically sound G-JF thermostat for Langevin Dynamics simulations into the ESPREesSo molecular package for large-scale simulations of soft matter systems. The implemented integration method is tested against the integrator currently used by the molecular package in simulations of a fluid bilayer membrane. While the latter exhibits deviations in the sampling statistics that increase with the integration time step dt, the former reproduces near-correct configurational statistics for all dt within the stability range of the simulations. We conclude that, with very modest revisions to existing codes, one can significantly improve the performance of statistical sampling using Langevin thermostats.
The capabilities of CP2K, a density-functional theory package and OMEN, a nano-device simulator, are combined to study transport phenomena from first-principles in unprecedentedly large nanostructures. Based on the Hamiltonian and overlap matrices generated by CP2K for a given system, OMEN solves the Schroedinger equation with open boundary conditions (OBCs) for all possible electron momenta and energies. To accelerate this core operation a robust algorithm called SplitSolve has been developed. It allows to simultaneously treat the OBCs on CPUs and the Schroedinger equation on GPUs, taking advantage of hybrid nodes. Our key achievements on the Cray-XK7 Titan are (i) a reduction in time-to-solution by more than one order of magnitude as compared to standard methods, enabling the simulation of structures with more than 50000 atoms, (ii) a parallel efficiency of 97% when scaling from 756 up to 18564 nodes, and (iii) a sustained performance of 15 DP-PFlop/s.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
