ترغب بنشر مسار تعليمي؟ اضغط هنا

Automated free energy calculation from atomistic simulations

90   0   0.0 ( 0 )
 نشر من قبل Sarath Menon
 تاريخ النشر 2021
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

We devise automated workflows for the calculation of Helmholtz and Gibbs free energies and their temperature and pressure dependence and provide the corresponding computational tools. We employ non-equilibrium thermodynamics for evaluating the free energy of solid and liquid phases at a given temperature and reversible scaling for computing free energies over a wide range of temperatures, including the direct integration of $P$-$T$ coexistence lines. By changing the chemistry and the interatomic potential, alchemical and upscaling free energy calculations are possible. Several examples illustrate the accuracy and efficiency of our implementation.

قيم البحث

اقرأ أيضاً

Diffuse scattering is a rich source of information about disorder in crystalline materials, which can be modelled using atomistic techniques such as Monte Carlo and molecular dynamics simulations. Modern X-ray and neutron scattering instruments can r apidly measure large volumes of diffuse-scattering data. Unfortunately, current algorithms for atomistic diffuse-scattering calculations are too slow to model large data sets completely, because the fast Fourier transform (FFT) algorithm has long been considered unsuitable for such calculations [Butler & Welberry, J. Appl. Cryst. 25, 391 (1992)]. Here, a new approach is presented for ultrafast calculation of atomistic diffuse-scattering patterns. It is shown that the FFT can actually be used to perform such calculations rapidly, and that a fast method based on sampling theory can be used to reduce high-frequency noise in the calculations. These algorithms are benchmarked using realistic examples of compositional, magnetic and displacive disorder. They accelerate the calculations by a factor of at least 100, making refinement of atomistic models to large diffuse-scattering volumes practical.
Atomistic effective Hamiltonian simulations are used to investigate electrocaloric (EC) effects in the lead-free Ba(Zr$_{0.5}$Ti$_{0.5}$)O$_{3}$ (BZT) relaxor ferroelectric. We find that the EC coefficient varies non-monotonically with the field at a ny temperature, presenting a maximum that can be traced back to the behavior of BZTs polar nanoregions. We also introduce a simple Landau-based model that reproduces the EC behavior of BZT as a function of field and temperature, and which is directly applicable to other compounds. Finally, we confirm that, for low temperatures (i.e., in non-ergodic conditions), the usual indirect approach to measure the EC response provides an estimate that differs quantitatively from a direct evaluation of the field-induced temperature change.
A coupled atomistic spin and lattice dynamics approach is developed which merges the dynamics of these two degrees of freedom into a single set of coupled equations of motion. The underlying microscopic model comprises local exchange interactions bet ween the electron spin and magnetic moment and the local couplings between the electronic charge and lattice displacements. An effective action for the spin and lattice variables is constructed in which the interactions among the spin and lattice components are determined by the underlying electronic structure. In this way, expressions are obtained for the electronically mediated couplings between the spin and lattice degrees of freedom, besides the well known inter-atomic force constants and spin-spin interactions. These former susceptibilities provide an atomistic ab initio description for the coupled spin and lattice dynamics. It is important to notice that this theory is strictly bilinear in the spin and lattice variables and provides a minimal model for the coupled dynamics of these subsystems and that the two subsystems are treated on the same footing. Questions concerning time-reversal and inversion symmetry are rigorously addressed and it is shown how these aspects are absorbed in the tensor structure of the interaction fields. By means of these results regarding the spin-lattice coupling, simple explanations of ionic dimerization in double anti-ferromagnetic materials, as well as, charge density waves induced by a non-uniform spin structure are given. In the final parts, a set of coupled equations of motion for the combined spin and lattice dynamics are constructed, which subsequently can be reduced to a form which is analogous to the Landau-Lifshitz-Gilbert equations for spin dynamics and damped driven mechanical oscillator for the ...
Statistical learning methods show great promise in providing an accurate prediction of materials and molecular properties, while minimizing the need for computationally demanding electronic structure calculations. The accuracy and transferability of these models are increased significantly by encoding into the learning procedure the fundamental symmetries of rotational and permutational invariance of scalar properties. However, the prediction of tensorial properties requires that the model respects the appropriate geometric transformations, rather than invariance, when the reference frame is rotated. We introduce a formalism that can be used to perform machine-learning of tensorial properties of arbitrary rank for general molecular geometries. To demonstrate it, we derive a tensor kernel adapted to rotational symmetry, which is the natural generalization of the smooth overlap of atomic positions (SOAP) kernel commonly used for the prediction of scalar properties at the atomic scale. The performance and generality of the approach is demonstrated by learning the instantaneous electrical response of water oligomers of increasing complexity, from the isolated molecule to the condensed phase.
The phase diagram of numerous materials of technological importance features high-symmetry high-temperature phases that exhibit phonon instabilities. Leading examples include shape-memory alloys, as well as ferroelectric, refractory, and structural m aterials. The thermodynamics of these phases have proven challenging to handle by atomistic computational thermodynamic techniques, due to the occurrence of constant anharmonicity-driven hopping between local low-symmetry distortions, while maintaining a high-symmetry time-averaged structure. To compute the free energy in such phases, we propose to explore the systems potential-energy surface by discrete sampling of local minima by means of a lattice gas Monte Carlo approach and by continuous sampling by means of a lattice dynamics approach in the vicinity of each local minimum. Given the proximity of the local minima, it is necessary to carefully partition phase space by using a Voronoi tessellation to constrain the domain of integration of the partition function, in order to avoid double-counting artifacts and enable an accurate harmonic treatment near each local minima. We consider the bcc phase of titanium as a prototypical examples to illustrate our approach.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا