Sort by
End-to-End Differentiable Force Field Generator with Crystal Structure Differentiation and Matching

In this paper, we present a new force field (FF) parameterization method with direct matching of crystal structures and atomic charges optimization in end-to-end differentiable manner. The advancement of force field (FF) parameterization methods has been accelerated by differentiable programming. Automatic differentiation (AD) has facilitated energy and force matching by differentiating these quantities with respect to the FF parameters, which we mention as force differentiation and matching (FDM). Nevertheless, crystal structure matching with AD remains difficult due to the converged structures optimized by the iterative algorithm being non-differentiable with respect to the FF parameters. To overcome this limitation, we introduce structure differentiation and matching (SDM) technique for generating FFs of small organic molecules using reference data, including stable monomer structures, crystal structures, lattice energies, and potential energy surfaces (PESs) of dihedral angles. SDM employs implicit function differentiation (IFD) and differentiable Ewald techniques to optimize FF parameters and atomic charges correspondingly. Our case study of eight exemplified molecules demonstrates that SDM substantially outperforms the conventional FDM, with error factors reduced to less than one-quarter with the charge optimization method called SDM(q-opt). This strategy achieves remarkable precision in reproducing lattice constants, atomic configurations, lattice energies, and PESs. Furthermore, molecular dynamics simulations confirm the stability of the generated crystal structures. This method can be adapted to other FF categories, such as polarized FFs and those with explicit hydrogen bonding interactions. We foresee that SDM(q-opt) will emerge as a standard method for parameterizing FFs using crystal structures.

Open Access
Relevant
End-to-End Differentiable Force Field Generator with Crystal Structure Differentiation and Matching

Differentiable programming has accelerated the development of force-field (FF) parameterization techniques. Specifically, automatic differentiation (AD) facilitates energy and force matching by differentiating them with respect to the FF parameters; hereinafter, referred to as force differentiation and matching (FDM). Conversely, crystal structure matching with AD has persisted as a challenge because the converged structures optimized by the iterative algorithm cannot be differentiated with respect to the FF parameters. Therefore, in this paper, we propose a structure differentiation and matching (SDM) method, wherein the converged structures are directly differentiated using the parameters with implicit function differentiation and matched with the experimental crystal structures. Subsequently, with a case study, we compared the reproducibility of the crystal structures, internal atomic coordinates, and lattice energies on eight exemplary molecules with the differentiable Ewald method for long-range interactions. The results indicated that SDM outperformed FDM on all three criteria. The FFs generated by SDM reproduced the lattice constants with a mean error of 0.56 %, the internal atomic coordinates with an error of 0.16 Angstrom, and the lattice energies with an error of 0.14 kcal/mol. The corresponding accuracies obtained with FDM were 1.2 %, 0.22 Angstrom, and 2.40 kcal/mol, respectively. Furthermore, we performed molecular dynamic simulations on a supercell, containing more than 3000 atoms, to confirm if the crystal structures were preserved under temperature fluctuations at 300 K.Overall, this method is not limited to Amber-type FFs and can be easily applied to the other types of FFs. Thus, we believe that SDM will emerge as one of the new standards for parameterizing FFs with crystal structures.

Open Access
Relevant
The Group Egogram: The Method of Numeric Expression for Group’s Character in Discussions

This paper describes the following two points. (1) A proposal for the method of derivation for the group egogram, (2) a confirmation that the derived group egogram can be interpreted and used in the same way as the individual egogram. The group egogram is a method of numeric expression for group’s character. In this study, we combined two methods for quantification of group’s character: the “egograms” and the “group decision-making stress method.” The egogram is a method used to analyze an individual’s character, while the stress method is a type of group decision making support method. In the field of team sports research, the egogram has been used for team analysis, and the team’s egogram is derived from the arithmetic mean of the egograms of each individual team member. However, the shape of the graph produced by the five indicators is important for the egogram. Therefore, it is not appropriate to use simple arithmetic averages with other member’s egograms that differ in overall scale. For this reason, it is necessary to find a method to compare the graphs with the same scale without changing their shape. In mathematical views, the stress method can “align the scales of everyone’s egogram while keeping the meaning (shape) of each individual’s egogram, and then take the average.” In other words, we can derive the group’s egogram without contradiction mathematically. We use these methods to establish a method for deriving the “group egogram,” which gives a numeric expression for group’s character. In addition, we conduct evaluation experiments to confirm whether the mathematically derived group egogram really expresses the character of the group. The four evaluation items are “the risk taking tendency,” “the team orientation,” “the punctual tendency,” and “the harmony tendency.” As a result, we were able to confirm effectiveness of the group egogram on all evaluation items.

Open Access
Relevant