-
Thermalized and mixed meanfield ADP potentials for magnesium hydrides
Authors:
M. Molinos,
M. Ortiz,
M. P. Ariza
Abstract:
We develop meanfield approximation and numerical quadrature schemes for the evaluation of Angular-Dependent interatomic Potentials (ADPs) for magnesium and magnesium hydrides at finite temperature (thermalization) and arbitrary atomic molar fractions (mixing) within a non-equilibrium statistical mechanical framework and derive local equilibrium relations. We numerically verify and experimentally v…
▽ More
We develop meanfield approximation and numerical quadrature schemes for the evaluation of Angular-Dependent interatomic Potentials (ADPs) for magnesium and magnesium hydrides at finite temperature (thermalization) and arbitrary atomic molar fractions (mixing) within a non-equilibrium statistical mechanical framework and derive local equilibrium relations. We numerically verify and experimentally validate the accuracy and fidelity of the resulting thermalized/mixed ADPs (TADPs) by means of selected numerical tests including free entropy, heat capacity, thermal expansion, molar volumes, equation of state and elastic constants. We show that the local equilibrium properties predicted by TADPs agree closely with those computed directly from ADP by means of Molecular Dynamics (MD).
△ Less
Submitted 13 July, 2024;
originally announced July 2024.
-
Fracture and size effect in mechanical metamaterials
Authors:
J. Ulloa,
M. P. Ariza,
J. E. Andrade,
M. Ortiz
Abstract:
We resort to variational methods to evaluate the asymptotic behavior of fine metamaterials as a function of cell size. To zeroth order, the metamaterial behaves as a micropolar continuum with both displacement and rotation degrees of freedom, but exhibits linear-elastic fracture mechanics scaling and therefore no size effect. To higher order, the overall energetics of the metastructure can be char…
▽ More
We resort to variational methods to evaluate the asymptotic behavior of fine metamaterials as a function of cell size. To zeroth order, the metamaterial behaves as a micropolar continuum with both displacement and rotation degrees of freedom, but exhibits linear-elastic fracture mechanics scaling and therefore no size effect. To higher order, the overall energetics of the metastructure can be characterized explicitly in terms of the solution of the zeroth-order continuum problem by the method of Γ-expansion. We present explicit expressions of the second-order correction for octet frames. As an application, we evaluate the compliance of double-cantilever octet specimens to second order and use the result to elucidate the dependence of the apparent toughness of the specimen on cell size. The analysis predicts the discreteness of the metamaterial lattice to effectively shield the crack-tip, a mechanism that we term lattice shielding. The theory specifically predicts anti-shielding, i. e., coarser is weaker, in agreement with recent experimental observations.
△ Less
Submitted 26 June, 2024;
originally announced July 2024.
-
Interatomic-Potential-Free, Data-Driven Molecular Dynamics
Authors:
J. Bulin,
J. Hamaekers,
M. P. Ariza,
M. Ortiz
Abstract:
We present a Data-Driven (DD) paradigm that enables molecular dynamics calculations to be performed directly from sampled force-field data such as obtained, e.g., from ab initio calculations, thereby eschewing the conventional step of modeling the data by empirical interatomic potentials entirely. The data required by the DD solvers consists of local atomic configurations and corresponding atomic…
▽ More
We present a Data-Driven (DD) paradigm that enables molecular dynamics calculations to be performed directly from sampled force-field data such as obtained, e.g., from ab initio calculations, thereby eschewing the conventional step of modeling the data by empirical interatomic potentials entirely. The data required by the DD solvers consists of local atomic configurations and corresponding atomic forces and is, therefore, fundamental, i.e., it is not beholden to any particular model. The resulting DD solvers, including a fully explicit DD-Verlet algorithm, are provably convergent and exhibit robust convergence with respect to the data in selected test cases. We present an example of application to C60 buckminsterfullerenes that showcases the feasibility, range and scope of the DD molecular dynamics paradigm.
△ Less
Submitted 5 August, 2022;
originally announced August 2022.
-
Adiabatic lapse rate of real gases
Authors:
Bogar Díaz,
Miguel Ángel García Ariza,
J. E. Ramírez
Abstract:
We derive a formula for the dry adiabatic lapse rate of atmospheres composed of real gases. We restrict our study to those described by a family of two-parameter cubic equations of state and the recent Guevara non-cubic equation. Since our formula depends on the adiabatic curves, we compute them all at once, considering molecules that can move, rotate, and vibrate, for any equation of state. To il…
▽ More
We derive a formula for the dry adiabatic lapse rate of atmospheres composed of real gases. We restrict our study to those described by a family of two-parameter cubic equations of state and the recent Guevara non-cubic equation. Since our formula depends on the adiabatic curves, we compute them all at once, considering molecules that can move, rotate, and vibrate, for any equation of state. To illustrate our results, we estimate the lapse rate of the troposphere of Titan, obtaining a better approximation to the observed data in some instances, when compared to the estimation provided by the virial expansion up to third order.
△ Less
Submitted 19 July, 2022; v1 submitted 22 March, 2022;
originally announced March 2022.
-
A Line-free Method of Monopoles for 3D Dislocation Dynamics
Authors:
A. Deffo,
M. P. Ariza,
M. Ortiz
Abstract:
We develop an approximation scheme for three-dimensional dislocation dynamics in which the dislocation line density is concentrated at points, or monopoles. Every monopole carries a Burgers vector and an element of line. The monopoles move according to mobility kinetics driven by elastic and applied forces. The divergence constraint, expressing the requirement that the monopoles approximate a boun…
▽ More
We develop an approximation scheme for three-dimensional dislocation dynamics in which the dislocation line density is concentrated at points, or monopoles. Every monopole carries a Burgers vector and an element of line. The monopoles move according to mobility kinetics driven by elastic and applied forces. The divergence constraint, expressing the requirement that the monopoles approximate a boundary, is enforced weakly. The fundamental difference with traditional approximation schemes based on segments is that in the present approach an explicit linear connectivity, or 'sequence', between the monopoles need not be defined. Instead, the monopoles move as an unstructured point set subject to the weak divergence constraint. In this sense, the new paradigm is 'line-free', i. e., it sidesteps the need to track dislocation lines. This attribute offers significant computational advantages in terms of simplicity, robustness and efficiency, as demonstrated by means of selected numerical examples.
△ Less
Submitted 8 June, 2018;
originally announced June 2018.