Thanks for the reply! This is very useful, particularly the paper you linked. With my research right now, I'm trying to come up with a visual representation of a "characteristic" atomic neighborhood around ions at different energy levels. More specifically, how the distribution of atoms surrounding a high vs low energy ion are different. But visually, there is no easily discernible difference, even though the radial distribution functions are very different for each energy level. That's why I'm studying other forms of distribution representations. Are you a member of one of the groups you listed, or just learning about it for your own research?
Yes, I'm part of the GAP research group (with Albert Bartok Partay and Gabor Csanyi). Representing environments is definitely still an ongoing research project. Other things I've played with in the past include identifying crystal structures at finite temperature (i.e. a classification rather than regression task), or differentiating between amorphous phases (since e.g. water has two amorphous solid phases with a 1st order transition in between, but there is no way in hell you would be able to tell one from the other visually.
We're currently working on a really ambitious new way to represent environments, but it's really preliminary at the moment.
Regarding your ion issue, what about the angular components? The radial functions really only tell you so much...
But more fundamentally, what do you mean by ion energy levels? I'm presuming you mean a metallic nucleus+core electrons, in a condensed phase at finite temperature. But of course that `atomic energy' -insofar as it exists- is a continuous function of position and not quantised, so I'm unsure what you mean by energy levels in this context.
That sounds like really exciting research. I guess the ultimate goal in representing environments would be to construct a function that captures all the statistically relevant information of a particular "type" of atomic neighborhood. It would reduce the degrees of freedom to the bare minimum necessary to recreate a similar environment that correctly reproduces any property of interest (kind of like data compression for materials). Would that be right?
The deal with the ions is that we are studying the transport and storage of lithium in a new type of carbon anode. The anode consists of small crystalline domains distributed throughout an amorphous carbon matrix. In order to capture all of the features of this material, we end up with a system of almost a million atoms. At the end of an equilibration run, the lithium ions can be found at different locations within the carbon matrix. It turns out that the potential energy of the lithium ions (as computed from the reactive potential the simulation was performed with) has a wide range of values. So we can sort these ions into bins of a histogram (this is what I meant by "energy levels"). And because there are so many samples, the radial distribution functions (RDFs) can be computed for all Li-C pairs in each separate energy bin. (The computed RDFs are useful because the results can be compared to the experimental RDFs obtained from neutron scattering.)
However, zooming in on ions of different potential energies reveals very little visual difference in the local environment. Yet we know there is definitely a difference in the structure because of the RDFs, but we cannot get a good understanding of it or provide a good representation of it. So that's when I discovered the 2010 paper by Bartók. As you mentioned, I want to figure out how the entire local atomic neighborhood affects the energy (as opposed to only the radial component), and I also want to create a 3D graphic that provides a good visualization of the differences between the atomic neighborhoods. In that sense, I need something that would compare the atomic neighborhoods without regard to translation, rotation, reflection, or permutation of identical atoms.
