Interatomic Potentials Repository

Warning! Note that elemental potentials taken from alloy descriptions may not work well for the pure species. This is particularly true if the elements were fit for compounds instead of being optimized separately. As with all interatomic potentials, please check to make sure that the performance is adequate for your problem.

2010--Kim-H-K-Jung-W-S-Lee-B-J--Nb-C

Citation: H.-K. Kim, W.-S. Jung, and B.-J. Lee (2010), "Modified embedded-atom method interatomic potentials for the Nb-C, Nb-N, Fe-Nb-C, and Fe-Nb-N systems", Journal of Materials Research 25(7), 1288-1297. DOI: 10.1557/jmr.2010.0182.

Abstract: Modified embedded-atom method (MEAM) interatomic potentials for Nb-C, Nb-N, Fe-Nb-C, and Fe-Nb-N systems have been developed based on the previously developed MEAM potentials for lower order systems. The potentials reproduce various fundamental physical properties (structural properties, elastic properties, thermal properties, and surface properties) of NbC and NbN, and interfacial energy between bcc Fe and NbC or NbN, in generally good agreement with higher-level calculations or experimental information. The applicability of the present potentials to atomic-level investigations to the precipitation behavior of complex-carbonitrides (Nb,Ti)(C,N) as well as NbC and NbN, and their effects on the mechanical properties of steels are also discussed.

Related Models:
2005--Lee-B-J-Lee-J-W--C

2006--Lee-B-J--Fe-C

2008--Kim-Y-M-Lee-B-J--Ti-C

2009--Kim-H-K-Jung-W-S-Lee-B-J--Fe-Ti-C

2010--Kim-H-K-Jung-W-S-Lee-B-J--Fe-Nb-C

2020--Jeong-G-U-Lee-B-J--Pd-C

2020--Jeong-G-U-Lee-B-J--Pt-C
2010--Kim-H-K-Jung-W-S-Lee-B-J--Fe-Nb-C
2001--Lee-B-J-Baskes-M-I-Kim-H-Cho-Y-K--Nb

2008--Sa-I-Lee-B--Fe-Nb

2010--Kim-H-K-Jung-W-S-Lee-B-J--Fe-Nb-C

2022--Fereidonnejad-R-Moghaddam-A-O-Moaddeli-M--Al-Hf-Nb-Ta-Ti-Zr

2022--Fereidonnejad-R-Moghaddam-A-O-Moaddeli-M--Al-Nb
LAMMPS pair_style meam (2010--Kim-H-K--Nb-C--LAMMPS--ipr1)

See Computed Properties
Notes: These files are based on files obtained from http://cmse.postech.ac.kr/home_2nnmeam.
File(s):
library.meam
NbC.meam
OpenKIM (MO_145522277939)

See Computed Properties
Notes: Listing found at https://openkim.org.
Link(s):
KIM page MEAM_2NN_Fe_to_Ga__MO_145522277939_001

Implementation Information

This page displays computed properties for the 2010--Kim-H-K--Nb-C--LAMMPS--ipr1 implementation of the 2010--Kim-H-K-Jung-W-S-Lee-B-J--Nb-C potential. Computed values for other implementations can be seen by clicking on the links below:

Diatom Energy vs. Interatomic Spacing

Plots of the potential energy vs interatomic spacing, r, are shown below for all diatom sets associated with the interatomic potential. This calculation provides insights into the functional form of the potential's two-body interactions. A system consisting of only two atoms is created, and the potential energy is evaluated for the atoms separated by 0.02 Å <= r <= 6.0> Å in intervals of 0.02 Å. Two plots are shown: one for the "standard" interaction distance range, and one for small values of r. The small r plot is useful for determining whether the potential is suitable for radiation studies.

The calculation method used is available as the iprPy diatom_scan calculation method.

Clicking on the image of a plot will open an interactive version of it in a new tab. The underlying data for the plots can be downloaded by clicking on the links above each plot.

Notes and Disclaimers:

These values are meant to be guidelines for comparing potentials, not the absolute values for any potential's properties. Values listed here may change if the calculation methods are updated due to improvements/corrections. Variations in the values may occur for variations in calculation methods, simulation software and implementations of the interatomic potentials.
As this calculation only involves two atoms, it neglects any multi-body interactions that may be important in molecules, liquids and crystals.
NIST disclaimer

Version Information:

2019-11-14. Maximum value range on the shortrange plots are now limited to "expected" levels as details are otherwise lost.
2019-08-07. Plots added.

Download data

Click on plot to load interactive version

2010--Kim-H-K--Nb-C--LAMMPS--ipr1/diatom

Click on plot to load interactive version

2010--Kim-H-K--Nb-C--LAMMPS--ipr1/diatom_short

Cohesive Energy vs. Interatomic Spacing

Plots of potential energy vs interatomic spacing, r, are shown below for a number of crystal structures. The structures are generated based on the ideal atomic positions and b/a and c/a lattice parameter ratios for a given crystal prototype. The size of the system is then uniformly scaled, and the energy calculated without relaxing the system. To obtain these plots, values of r are evaluated every 0.02 Å up to 6 Å.

The calculation method used is available as the iprPy E_vs_r_scan calculation method.

Clicking on the image of a plot will open an interactive version of it in a new tab. The underlying data for the plots can be downloaded by clicking on the links above each plot.

Notes and Disclaimers:

These values are meant to be guidelines for comparing potentials, not the absolute values for any potential's properties. Values listed here may change if the calculation methods are updated due to improvements/corrections. Variations in the values may occur for variations in calculation methods, simulation software and implementations of the interatomic potentials.
The minima identified by this calculation do not guarantee that the associated crystal structures will be stable since no relaxation is performed.
NIST disclaimer

Version Information:

2020-12-18. Descriptions, tables and plots updated to reflect that the energy values are the measuredper atom potential energy rather than cohesive energy as some potentials have non-zero isolated atom energies.
2019-02-04. Values regenerated with even r spacings of 0.02 Å, and now include values less than 2 Å when possible. Updated calculation method and parameters enhance compatibility with more potential styles.
2019-04-26. Results for hcp, double hcp, α-As and L1₀ prototypes regenerated from different unit cell representations. Only α-As results show noticable (>1e-5 eV) difference due to using a different coordinate for Wykoff site c position.
2018-06-13. Values for MEAM potentials corrected. Dynamic versions of the plots moved to separate pages to improve page loading. Cosmetic changes to how data is shown and updates to the documentation.
2017-01-11. Replaced png pictures with interactive Bokeh plots. Data regenerated with 200 values of r instead of 300.
2016-09-28. Plots for binary structures added. Data and plots for elemental structures regenerated. Data values match the values of the previous version. Data table formatting slightly changed to increase precision and ensure spaces between large values. Composition added to plot title and structure names made longer.
2016-04-07. Plots for elemental structures added.

Select a composition:

Download data

Click on plot to load interactive version

2010--Kim-H-K--Nb-C--LAMMPS--ipr1/EvsR.C

Crystal Structure Predictions

Computed lattice constants and cohesive/potential energies are displayed for a variety of crystal structures. The values displayed here are obtained using the following process.

Initial crystal structure guesses are taken from:
1. The iprPy E_vs_r_scan calculation results (shown above) by identifying all energy minima along the measured curves for a given crystal prototype + composition.
2. Structures in the Materials Project and OQMD DFT databases.
All initial guesses are relaxed using three independent methods using a 10x10x10 supercell:
1. "box": The system's lattice constants are adjusted to zero pressure without internal relaxations using the iprPy relax_box calculation with a strainrange of 1e-6.
2. "static": The system's lattice and atomic positions are statically relaxed using the iprPy relax_static calculation with a minimization force tolerance of 1e-10 eV/Angstrom.
3. "dynamic": The system's lattice and atomic positions are dynamically relaxed for 10000 timesteps of 0.01 ps using the iprPy relax_dynamic calculation with an nph integration plus Langevin thermostat. The final configuration is then used as input in running an iprPy relax_static calculation with a minimization force tolerance of 1e-10 eV/Angstrom.
The relaxed structures obtained from #2 are then evaluated using the spglib package to identify an ideal crystal unit cell based on the results.
The space group information of the ideal unit cells is compared to the space group information of the corresponding reference structures to identify which structures transformed upon relaxation. The structures that did not transform to a different structure are listed in the table(s) below. The "method" field indicates the most rigorous relaxation method where the structure did not transform. The space group information is also used to match the DFT reference structures to the used prototype, where possible.
The cohesive energy, E_coh, is calculated from the measured potential energy per atom, E_pot$, by subtracting the isolated energy averaged across all atoms in the unit cell. The isolated atom energies of each species model is obtained either by evaluating a single atom atomic configuration, or by identifying the first energy plateau from the diatom scan calculations for r > 2 Å.

The calculation methods used are implemented into iprPy as the following calculation styles

Notes and Disclaimers:

These values are meant to be guidelines for comparing potentials, not the absolute values for any potential's properties. Values listed here may change if the calculation methods are updated due to improvements/corrections. Variations in the values may occur for variations in calculation methods, simulation software and implementations of the interatomic potentials.
The presence of any structures in this list does not guarantee that those structures are stable. Also, the lowest energy structure may not be included in this list.
Multiple values for the same crystal structure but different lattice constants are possible. This is because multiple energy minima are possible for a given structure and interatomic potential. Having multiple energy minima for a structure does not necessarily make the potential "bad" as unwanted configurations may be unstable or correspond to conditions that may not be relevant to the problem of interest (eg. very high strains).
NIST disclaimer

Version Information:

2025-07-02. All "mp-" reference structure calculations were re-relaxed using the updated Materials Project database rather than the original database structures. Also, a bug was fixed that caused the "static" relaxations to occasionally throw unnecessary errors. This was fixed and all affected calculations were reset and performed again.
2022-05-27. The "box" method results have all been redone with an updated methodology more suited for non-orthogonal systems.
2020-12-18. Cohesive energies have been corrected by making them relative to the energies of the isolated atoms. The previous cohesive energy values are now listed as the potential energies.
2019-06-07. Structures with positive or near zero cohesive energies removed from the display tables. All values still present in the raw data files.
2019-04-26. Calculations now computed for each implementation. Results for hcp, double hcp, α-As and L1₀ prototypes regenerated from different unit cell representations.
2018-06-14. Methodology completely changed affecting how the information is displayed. Calculations involving MEAM potentials corrected.
2016-09-28. Values for simple compounds added. All identified energy minima for each structure are listed. The existing elemental data was regenerated. Most values are consistent with before, but some differences have been noted. Specifically, variations are seen with some values for potentials where the elastic constants don't vary smoothly near the equilibrium state. Additionally, the inclusion of some high-energy structures has changed based on new criteria for identifying when structures have relaxed to another structure.
2016-04-07. Values for elemental crystal structures added. Only values for the global energy minimum of each unique structure given.

Select a composition:

Download raw data (including filtered results)

Reference structure matches:
A1--Cu--fcc = oqmd-1214512
A15--beta-W = oqmd-1214957
A2--W--bcc = oqmd-1215135
A3'--alpha-La--double-hcp = oqmd-1215403
A3--Mg--hcp = oqmd-1215313
A4--C--dc = mp-66, oqmd-675640
A5--beta-Sn = oqmd-1215581
A6--In--bct = mp-1181996, oqmd-1215670
A7--alpha-As = oqmd-594528
Ah--alpha-Po--sc = mp-998866, oqmd-636457, oqmd-1214601

prototype	method	E_coh (eV/atom)	E_pot (eV/atom)	a₀ (Å)	b₀ (Å)	c₀ (Å)	α (degrees)	β (degrees)	γ (degrees)
A7--alpha-As	dynamic	-7.4036	-7.4036	2.4546	2.4546	11.0545	90.0	90.0	120.0
mp-568806	dynamic	-7.3959	-7.3959	2.4556	2.4556	3.6684	90.0	90.0	120.0
oqmd-610551	box	-7.3958	-7.3958	2.4553	2.4553	3.6684	90.0	90.0	120.0
A7--alpha-As	box	-7.3902	-7.3902	2.4554	2.4554	10.9361	90.0	90.0	120.0
oqmd-683820	dynamic	-7.3887	-7.3887	2.4574	4.2469	7.3637	90.0	90.0	90.0
oqmd-683820	box	-7.3887	-7.3887	2.4568	4.248	7.3653	90.0	90.0	90.0
oqmd-683819	box	-7.3887	-7.3887	2.4567	4.248	7.3655	90.0	90.0	90.0
oqmd-1214690	box	-7.3886	-7.3886	2.4563	4.2487	7.3664	90.0	90.0	90.0
mp-937760	box	-7.3886	-7.3886	2.4572	4.2472	7.3696	90.0	90.0	90.0
mp-48	static	-7.3882	-7.3882	2.4554	2.4554	7.269	90.0	90.0	120.0
oqmd-590101	box	-7.3882	-7.3882	2.4554	2.4554	7.269	90.0	90.0	120.0
mp-568286	box	-7.3881	-7.3881	2.456	4.2492	7.3852	90.0	90.0	90.0
oqmd-1215937	box	-7.388	-7.388	2.4554	2.4554	7.2691	90.0	90.0	120.0
mp-568363	box	-7.3852	-7.3852	2.4565	4.2482	7.5068	90.0	90.0	90.0
mp-569416	dynamic	-7.3781	-7.3781	2.454	2.454	57.2026	90.0	90.0	120.0
oqmd-589492	static	-7.3781	-7.3781	2.454	2.454	57.8364	90.0	90.0	120.0
oqmd-589492	dynamic	-7.3781	-7.3781	2.454	2.454	54.4312	90.0	90.0	120.0
mp-569416	static	-7.3781	-7.3781	2.454	2.454	57.7586	90.0	90.0	120.0
mp-569416	box	-7.375	-7.375	2.4545	2.4545	57.5676	90.0	90.0	120.0
oqmd-589492	box	-7.375	-7.375	2.4545	2.4545	57.7433	90.0	90.0	120.0
mp-606949	box	-7.3706	-7.3706	2.4542	2.4542	30.9086	90.0	90.0	120.0
oqmd-589493	box	-7.3706	-7.3706	2.4541	2.4541	31.3601	90.0	90.0	120.0
A4--C--dc	dynamic	-7.37	-7.37	3.5565	3.5565	3.5565	90.0	90.0	90.0
mp-569304	dynamic	-7.3682	-7.3682	2.453	2.453	34.4928	90.0	90.0	120.0
oqmd-589491	dynamic	-7.3682	-7.3682	2.453	2.453	33.5472	90.0	90.0	120.0
oqmd-589491	static	-7.3682	-7.3682	2.453	2.453	36.378	90.0	90.0	120.0
mp-569304	static	-7.3682	-7.3682	2.453	2.453	35.5364	90.0	90.0	120.0
oqmd-589491	box	-7.3627	-7.3627	2.4535	2.4535	35.9984	90.0	90.0	120.0
mp-569304	box	-7.3618	-7.3618	2.4536	2.4536	35.08	90.0	90.0	120.0
mp-616440	dynamic	-7.3595	-7.3595	2.5146	2.5146	16.4546	90.0	90.0	120.0
oqmd-599494	dynamic	-7.3582	-7.3582	2.5146	2.5146	43.2046	90.0	90.0	120.0
oqmd-599494	static	-7.3581	-7.3581	2.514	2.514	43.2249	90.0	90.0	120.0
mp-611448	dynamic	-7.3563	-7.3563	2.5145	2.5145	12.3486	90.0	90.0	120.0
mp-569517	dynamic	-7.3543	-7.3543	2.5143	2.5143	30.8852	90.0	90.0	120.0
oqmd-599492	box	-7.3517	-7.3517	2.5044	2.5044	16.6112	90.0	90.0	120.0
mp-616440	box	-7.3514	-7.3514	2.5042	2.5042	16.6142	90.0	90.0	120.0
mp-611426	dynamic	-7.3513	-7.3513	2.5141	2.5141	8.2416	90.0	90.0	120.0
mp-569567	box	-7.3492	-7.3492	2.5032	2.5032	43.6638	90.0	90.0	120.0
mp-611448	box	-7.3463	-7.3463	2.5018	2.5018	12.4936	90.0	90.0	120.0
mp-47	dynamic	-7.3444	-7.3444	2.5135	2.5135	4.1264	90.0	90.0	120.0
mp-569517	box	-7.3426	-7.3426	2.5001	2.5001	31.2918	90.0	90.0	120.0
oqmd-599493	box	-7.3418	-7.3418	2.4998	2.4998	31.3055	90.0	90.0	120.0
mp-611426	box	-7.3375	-7.3375	2.498	2.498	8.3647	90.0	90.0	120.0
oqmd-637444	static	-7.3341	-7.3341	4.0031	4.0031	4.0031	90.0	90.0	90.0
mp-997182	dynamic	-7.3335	-7.3335	2.4513	2.4513	12.0458	90.0	90.0	120.0
mp-1040425	dynamic	-7.3335	-7.3335	2.4513	2.4513	19.9983	90.0	90.0	120.0
mp-990448	dynamic	-7.3335	-7.3335	2.4513	2.4513	8.5238	90.0	90.0	120.0
mp-47	box	-7.3197	-7.3197	2.4908	2.4908	4.2157	90.0	90.0	120.0
oqmd-599489	box	-7.3192	-7.3192	2.4906	2.4906	4.2167	90.0	90.0	120.0
oqmd-637443	static	-7.2495	-7.2495	5.1722	5.1722	5.1722	90.0	90.0	90.0
mp-1096869	dynamic	-7.2471	-7.2471	25.4572	14.7005	7.4494	90.0	90.0	90.0
mp-1078845	dynamic	-7.2222	-7.2222	4.2586	8.966	2.463	90.0	90.0	90.0
mp-1008374	box	-7.1418	-7.1418	2.6678	7.8119	2.5224	90.0	90.0	90.0
oqmd-637443	box	-7.0941	-7.0941	5.2116	5.2116	5.2116	90.0	90.0	90.0
mp-1188817	box	-7.087	-7.087	5.2126	5.2126	5.2126	90.0	90.0	90.0
mp-1078845	box	-7.0853	-7.0853	4.3158	8.9323	2.4829	90.0	90.0	90.0
oqmd-610557	box	-7.0564	-7.0564	4.7365	5.6117	4.3375	90.0	90.0	90.0
oqmd-610556	box	-7.0533	-7.0533	4.7355	5.6117	4.3397	90.0	90.0	90.0
mp-568410	box	-7.0497	-7.0497	4.7356	5.6133	4.341	90.0	90.0	90.0
oqmd-637353	box	-7.0125	-7.0125	4.4996	4.4996	2.4535	90.0	90.0	90.0
mp-1008395	box	-7.0095	-7.0095	4.5004	4.5004	2.4538	90.0	90.0	90.0
mp-683919	static	-6.9215	-6.9215	8.8935	17.5394	17.8315	90.0	90.0	90.0
oqmd-689315	static	-6.9215	-6.9215	8.8929	17.5385	17.8327	90.0	90.0	90.0
mp-1245190	dynamic	-6.9026	-6.9026	9.3464	9.4496	9.6926	89.5	89.8	87.9
oqmd-684315	dynamic	-6.8862	-6.8862	8.6623	13.8597	9.0025	90.0	90.0	90.0
mp-1196857	dynamic	-6.8711	-6.8711	7.8828	16.3193	3.9261	90.0	90.0	90.0
mp-680372	dynamic	-6.8688	-6.8688	9.1328	9.1328	24.7083	90.0	90.0	120.0
oqmd-614417	static	-6.8688	-6.8688	9.1328	9.1328	24.7079	90.0	90.0	120.0
mp-680372	static	-6.8688	-6.8688	9.1364	9.1364	24.8296	90.0	90.0	120.0
mp-683919	box	-6.8684	-6.8684	10.1639	17.5242	18.0335	90.0	90.0	90.0
oqmd-614417	dynamic	-6.8677	-6.8677	9.1329	9.1329	24.6311	90.0	90.0	120.0
mp-1182684	dynamic	-6.8616	-6.8616	6.6255	7.9062	13.6708	90.0	90.0	90.0
mp-630227	static	-6.8529	-6.8529	9.0541	9.1263	14.8231	90.0	90.0	90.0
mp-1147718	dynamic	-6.8513	-6.8513	10.3519	13.117	9.0774	90.0	90.0	90.0
mp-1147718	static	-6.8513	-6.8513	10.3516	13.117	9.0773	90.0	90.0	90.0
oqmd-684627	dynamic	-6.85	-6.85	9.0944	9.0944	14.503	90.0	90.0	90.0
oqmd-689315	box	-6.8468	-6.8468	10.0958	17.5905	18.0366	90.0	90.0	90.0
mp-1244964	dynamic	-6.8409	-6.8409	8.309	10.519	10.8624	86.5	87.7	79.6
oqmd-603085	static	-6.8366	-6.8366	13.6545	13.6545	13.6545	90.0	90.0	90.0
mp-568028	static	-6.8323	-6.8323	8.2306	14.0341	9.2616	90.0	90.0	90.0
mp-1205283	static	-6.831	-6.831	9.323	13.2657	9.0555	90.0	90.0	90.0
oqmd-614417	box	-6.829	-6.829	9.2421	9.2421	24.2629	90.0	90.0	120.0
mp-1080826	box	-6.8285	-6.8285	9.4175	2.5749	4.2133	90.0	98.2	90.0
oqmd-637352	box	-6.8269	-6.8269	9.4157	2.5753	4.2139	90.0	98.2	90.0
oqmd-595714	box	-6.8265	-6.8265	9.0922	9.2116	14.6616	90.0	90.0	90.0
mp-680372	box	-6.8262	-6.8262	9.2674	9.2674	25.2892	90.0	90.0	120.0
mp-667273	static	-6.8226	-6.8226	14.2353	14.2353	14.2353	90.0	90.0	90.0
oqmd-677227	dynamic	-6.8226	-6.8226	13.9442	13.9442	13.9442	90.0	90.0	90.0
mp-1196583	static	-6.8226	-6.8226	13.9445	13.9445	13.9445	90.0	90.0	90.0
mp-630227	box	-6.8127	-6.8127	9.0963	9.2025	15.6512	90.0	90.0	90.0
oqmd-603085	box	-6.8083	-6.8083	13.9525	13.9525	13.9525	90.0	90.0	90.0
mp-1147718	box	-6.8058	-6.8058	10.2719	14.7116	9.2256	90.0	90.0	90.0
mp-1190171	box	-6.805	-6.805	9.3042	2.58	4.2181	90.0	90.0	90.0
oqmd-637130	box	-6.8039	-6.8039	9.3024	2.5804	4.2184	90.0	90.0	90.0
mp-667273	box	-6.8027	-6.8027	14.4166	14.4166	14.4166	90.0	90.0	90.0
oqmd-677227	box	-6.7948	-6.7948	14.1406	14.1406	14.1406	90.0	90.0	90.0
mp-1196583	box	-6.7937	-6.7937	14.3577	14.3577	14.3577	90.0	90.0	90.0
oqmd-684627	box	-6.793	-6.793	9.176	9.176	15.276	90.0	90.0	90.0
mp-624889	box	-6.678	-6.678	4.7327	5.1652	5.668	90.0	90.0	90.0
oqmd-610552	box	-6.5849	-6.5849	4.7338	5.2685	6.5227	90.0	90.0	90.0
mp-568028	box	-6.5807	-6.5807	8.6235	14.2731	9.6082	90.0	90.0	90.0
mp-579909	box	-6.5756	-6.5756	4.733	5.267	6.9465	90.0	90.0	90.0
mp-1192619	dynamic	-6.5706	-6.5706	7.8231	7.8231	7.8231	90.0	90.0	90.0
mp-1244964	static	-6.5357	-6.5357	9.4897	9.774	10.8485	86.4	87.9	88.9
mp-1245190	static	-6.5092	-6.5092	9.5287	10.1674	10.1756	88.8	88.4	85.8
oqmd-684315	static	-6.4535	-6.4535	9.2672	15.7225	8.9224	90.0	90.0	90.0
mp-1244913	dynamic	-6.4004	-6.4004	9.2782	12.5595	12.6256	98.3	91.7	104.1
mp-1244913	static	-6.3916	-6.3916	9.528	10.2556	10.3809	82.7	82.9	88.8
mp-570002	static	-6.2284	-6.2284	5.0331	5.0331	5.0331	90.0	90.0	90.0
mp-24	box	-5.9897	-5.9897	4.689	4.689	4.689	90.0	90.0	90.0
mp-1196857	static	-5.7295	-5.7295	6.2546	18.1878	5.3434	90.0	90.0	90.0
mp-1205283	box	-5.7105	-5.7105	11.3543	13.4589	9.3143	90.0	90.0	90.0
mp-1056957	static	-5.6229	-5.6229	6.5479	3.1893	2.348	90.0	90.0	90.0
oqmd-585286	box	-5.4434	-5.4434	5.3564	5.3564	5.3564	90.0	90.0	90.0
mp-570002	box	-5.4423	-5.4423	5.3568	5.3568	5.3568	90.0	90.0	90.0
mp-1196857	box	-5.4349	-5.4349	8.715	19.2545	6.2116	90.0	90.0	90.0
mp-1197903	static	-5.4327	-5.4327	7.3096	15.5264	20.2166	76.2	88.5	76.8
oqmd-684315	box	-5.3916	-5.3916	10.6359	15.6987	9.6469	90.0	90.0	90.0
mp-1095633	dynamic	-4.7619	-4.7619	8.2206	8.2206	14.4908	90.0	90.0	90.0
mp-1095633	static	-4.7619	-4.7619	7.8933	7.8933	14.4509	90.0	90.0	90.0
mp-1245190	box	-4.6282	-4.6282	10.8584	12.0341	12.5891	84.5	90.0	88.5
A7--alpha-As	box	-4.3092	-4.3092	2.8301	2.8301	7.2453	90.0	90.0	120.0
A3'--alpha-La--double-hcp	box	-3.9973	-3.9973	2.8086	2.8086	4.9645	90.0	90.0	120.0
A5--beta-Sn	box	-3.953	-3.953	3.8633	3.8633	2.1405	90.0	90.0	90.0
Ah--alpha-Po--sc	static	-3.8984	-3.8984	2.0684	2.0684	2.0684	90.0	90.0	90.0
mp-1095534	box	-3.7946	-3.7946	13.0965	13.0965	2.6067	90.0	90.0	90.0
mp-1182684	box	-3.72	-3.72	10.3113	15.7487	20.5107	90.0	90.0	90.0
A6--In--bct	static	-3.6933	-3.6933	2.0015	2.0015	6.2095	90.0	90.0	90.0
mp-1205417	box	-3.5761	-3.5761	17.3176	17.3176	10.1475	90.0	90.0	90.0
oqmd-1214868	box	-3.4783	-3.4783	9.4301	9.4301	9.4301	90.0	90.0	90.0
mp-1192619	box	-3.4108	-3.4108	13.1054	13.1054	13.1054	90.0	90.0	90.0
mp-1097832	static	-3.4035	-3.4035	2.9722	2.9722	9.8799	90.0	90.0	90.0
A2--W--bcc	box	-2.9173	-2.9173	2.6077	2.6077	2.6077	90.0	90.0	90.0
A3--Mg--hcp	box	-2.8568	-2.8568	2.2337	2.2337	5.7752	90.0	90.0	120.0
mp-1095633	box	-2.6964	-2.6964	16.8903	16.8903	3.8226	90.0	90.0	90.0
A3--Mg--hcp	box	-2.529	-2.529	2.3834	2.3834	3.9133	90.0	90.0	120.0
A3'--alpha-La--double-hcp	box	-2.5271	-2.5271	2.3863	2.3863	7.8152	90.0	90.0	120.0
A1--Cu--fcc	static	-2.5252	-2.5252	3.3791	3.3791	3.3791	90.0	90.0	90.0
A15--beta-W	static	-2.4762	-2.4762	4.3042	4.3042	4.3042	90.0	90.0	90.0

prototype	method	E_coh (eV/atom)	E_pot (eV/atom)	a₀ (Å)	b₀ (Å)	c₀ (Å)	α (degrees)	β (degrees)	γ (degrees)
C1--CaF2--fluorite	static	-7.3205	-7.3205	5.1532	5.1532	5.1532	90.0	90.0	90.0
C1--CaF2--fluorite	static	-0.0	-0.0	10.6232	10.6232	10.6232	90.0	90.0	90.0
C1--CaF2--fluorite	box	0.1412	0.1412	9.1886	9.1886	9.1886	90.0	90.0	90.0

prototype	method	E_coh (eV/atom)	E_pot (eV/atom)	a₀ (Å)	b₀ (Å)	c₀ (Å)	α (degrees)	β (degrees)	γ (degrees)
L1_2--AuCu3	static	-5.831	-5.831	3.457	3.457	3.457	90.0	90.0	90.0
A15--Cr3Si	static	-5.7552	-5.7552	4.3892	4.3892	4.3892	90.0	90.0	90.0
D0_3--BiF3	box	-5.3264	-5.3264	5.3712	5.3712	5.3712	90.0	90.0	90.0
D0_3--BiF3	static	-0.0	-0.0	10.6232	10.6232	10.6232	90.0	90.0	90.0

prototype	method	E_coh (eV/atom)	E_pot (eV/atom)	a₀ (Å)	b₀ (Å)	c₀ (Å)	α (degrees)	β (degrees)	γ (degrees)
B1--NaCl--rock-salt	dynamic	-8.15	-8.15	4.56	4.56	4.56	90.0	90.0	90.0
B3--ZnS--cubic-zinc-blende	static	-7.5556	-7.5556	4.6861	4.6861	4.6861	90.0	90.0	90.0
B2--CsCl	static	-7.0422	-7.0422	2.8735	2.8735	2.8735	90.0	90.0	90.0
B3--ZnS--cubic-zinc-blende	static	-0.0	-0.0	10.6232	10.6232	10.6232	90.0	90.0	90.0

prototype	method	E_coh (eV/atom)	E_pot (eV/atom)	a₀ (Å)	b₀ (Å)	c₀ (Å)	α (degrees)	β (degrees)	γ (degrees)
C1--CaF2--fluorite	static	-7.1466	-7.1466	5.5282	5.5282	5.5282	90.0	90.0	90.0

prototype	method	E_coh (eV/atom)	E_pot (eV/atom)	a₀ (Å)	b₀ (Å)	c₀ (Å)	α (degrees)	β (degrees)	γ (degrees)
D0_3--BiF3	static	-7.1502	-7.1502	6.0622	6.0622	6.0622	90.0	90.0	90.0
A15--Cr3Si	box	-7.0625	-7.0625	4.9001	4.9001	4.9001	90.0	90.0	90.0
L1_2--AuCu3	static	-6.8237	-6.8237	3.8712	3.8712	3.8712	90.0	90.0	90.0
A15--Cr3Si	box	-1.6285	-1.6285	7.4786	7.4786	7.4786	90.0	90.0	90.0

prototype	method	E_coh (eV/atom)	E_pot (eV/atom)	a₀ (Å)	b₀ (Å)	c₀ (Å)	α (degrees)	β (degrees)	γ (degrees)
A15--beta-W	dynamic	-7.4754	-7.4754	5.2671	5.2671	5.2671	90.0	90.0	90.0
A2--W--bcc	dynamic	-7.47	-7.47	3.3024	3.3024	3.3024	90.0	90.0	90.0
mp-1104341	box	-7.4689	-7.4689	2.86	4.6701	9.7687	103.8	95.6	90.0
oqmd-1214814	dynamic	-7.4258	-7.4258	10.1993	10.1993	10.1993	90.0	90.0	90.0
oqmd-1214814	box	-7.4233	-7.4233	10.1975	10.1975	10.1975	90.0	90.0	90.0
mp-1094120	box	-7.3616	-7.3616	5.4115	6.6594	7.2901	107.9	110.0	101.9
oqmd-1215259	box	-7.3395	-7.3395	2.7783	5.4623	4.7614	90.0	90.0	90.0
A3--Mg--hcp	box	-7.3365	-7.3365	2.9666	2.9666	4.7495	90.0	90.0	120.0
A6--In--bct	static	-7.3358	-7.3358	2.7705	2.7705	4.7075	90.0	90.0	90.0
A3'--alpha-La--double-hcp	static	-7.3313	-7.3313	2.9564	2.9564	9.5657	90.0	90.0	120.0
oqmd-1214903	box	-7.3303	-7.3303	7.169	7.169	7.169	90.0	90.0	90.0
A1--Cu--fcc	box	-7.3263	-7.3263	4.1682	4.1682	4.1682	90.0	90.0	90.0
oqmd-1215081	box	-7.143	-7.143	3.1085	4.8442	9.6548	90.0	90.0	90.0
oqmd-1214725	box	-7.1386	-7.1386	3.1056	4.8613	9.6889	90.0	90.0	90.0
A5--beta-Sn	box	-7.0229	-7.0229	4.8378	4.8378	3.0666	90.0	90.0	90.0
oqmd-1215972	box	-6.6919	-6.6919	4.1446	4.1446	5.1185	90.0	90.0	120.0
Ah--alpha-Po--sc	static	-6.5806	-6.5806	2.6056	2.6056	2.6056	90.0	90.0	90.0
A7--alpha-As	box	-6.5284	-6.5284	3.6091	3.6091	9.6092	90.0	90.0	120.0
A4--C--dc	static	-6.0611	-6.0611	5.4209	5.4209	5.4209	90.0	90.0	90.0

Elastic Constants Predictions

Static elastic constants are displayed for the unique structures identified in Crystal Structure Predictions above. The values displayed here are obtained by measuring the change in virial stresses due to applying small strains to the relaxed crystals. The initial structure and the strained states are all relaxed using force minimization.

The calculation method used is available as the iprPy elastic_constants_static calculation method.

Notes and Disclaimers:

These values are meant to be guidelines for comparing potentials, not the absolute values for any potential's properties. Values listed here may change if the calculation methods are updated due to improvements/corrections. Variations in the values may occur for variations in calculation methods, simulation software and implementations of the interatomic potentials.
The presence of any structures in this list does not guarantee that those structures are stable.
The elastic constants have been computed for a variety of strains, and in some cases for slightly different lattice constant values. The static nature of this calculation can give poor predictions if the evaluated states straddle a functional discontinuity in the potential's third derivative. Be sure to compare the elastic constants for the different strains (positive and negative).
NIST disclaimer

Version Information:

2019-08-07. Data added.

Composition:

Prototype:

a₀:

strain:

Download raw data

Cij in GPa:

215.084	274.486	274.486	-0.000	0.000	0.000
274.486	215.084	274.486	0.000	-0.000	-0.000
274.486	274.486	215.084	-0.000	-0.000	-0.000
0.000	-0.000	-0.000	-79.796	-0.000	-0.000
0.000	0.000	-0.000	-0.000	-79.796	-0.000
0.000	-0.000	0.000	-0.000	-0.000	-79.796

215.084	274.486	274.486	0.000	0.000	0.000
274.486	215.084	274.486	-0.000	0.000	0.000
274.486	274.486	215.084	0.000	0.000	0.000
-0.000	-0.000	-0.000	-79.796	0.000	0.000
-0.000	-0.000	0.000	0.000	-79.796	0.000
-0.000	0.000	-0.000	0.000	0.000	-79.796

215.085	274.486	274.486	-0.000	0.000	0.000
274.486	215.085	274.486	0.000	0.000	-0.000
274.486	274.486	215.085	-0.000	-0.000	-0.000
0.000	0.000	0.000	-79.796	-0.000	-0.000
0.000	0.000	0.000	-0.000	-79.796	-0.000
0.000	0.000	0.000	-0.000	0.000	-79.796

215.084	274.486	274.486	0.000	0.000	0.000
274.486	215.084	274.486	0.000	0.000	0.000
274.486	274.486	215.084	0.000	0.000	0.000
-0.000	-0.000	-0.000	-79.796	0.000	0.000
-0.000	-0.000	-0.000	0.000	-79.796	-0.000
-0.000	-0.000	-0.000	0.000	-0.000	-79.796

215.086	274.487	274.487	-0.000	0.000	0.000
274.487	215.086	274.487	0.000	-0.000	-0.000
274.487	274.487	215.086	-0.000	-0.000	-0.000
0.000	0.000	0.000	-79.796	-0.000	-0.000
0.000	0.000	0.000	-0.000	-79.796	-0.000
0.000	0.000	0.000	0.000	0.000	-79.796

215.083	274.485	274.485	0.000	-0.000	-0.000
274.485	215.083	274.485	-0.000	0.000	0.000
274.485	274.485	215.083	0.000	0.000	0.000
-0.000	-0.000	-0.000	-79.796	0.000	0.000
-0.000	-0.000	-0.000	0.000	-79.796	-0.000
-0.000	-0.000	-0.000	0.000	-0.000	-79.796

59.279	312.236	312.236	-0.000	-0.000	-0.000
312.236	59.279	312.236	0.000	0.000	0.000
312.236	312.236	59.279	0.000	0.000	-0.000
-0.000	-0.000	-0.000	-340.729	0.000	0.000
-0.000	0.000	-0.000	0.000	-340.729	-0.000
-0.000	-0.000	-0.000	0.000	-0.000	-340.729

57.180	313.286	313.286	-0.000	-0.000	-0.000
313.286	57.180	313.286	-0.000	-0.000	-0.000
313.286	313.286	57.180	-0.000	0.000	-0.000
0.000	0.000	0.000	-340.730	-0.000	-0.000
0.000	0.000	0.000	-0.000	-340.730	0.000
0.000	0.000	0.000	0.000	-0.000	-340.729

59.279	312.236	312.236	-0.000	0.000	-0.000
312.236	59.279	312.236	0.000	0.000	-0.000
312.236	312.236	59.279	0.000	0.000	-0.000
-133272898.431	-132522967.282	-145308490.509	-8895809.641	-1202308.136	-5995124.730
-131849329.427	-128759213.333	-138361127.575	6292016.291	-11458649.470	-2605616.801
-138412999.703	-137770649.819	-136366776.248	1232351.337	-7080821.552	-4395871.766

59.279	312.236	312.236	0.000	-0.000	0.000
312.236	59.279	312.236	-0.000	-0.000	-0.000
312.236	312.236	59.279	-0.000	-0.000	-0.000
126146004.007	128687574.940	133932353.029	-11293978.254	-4010321.895	965214.681
133707998.110	123257384.976	130752304.472	-1059457.641	-10545033.637	990442.079
130095044.688	138537085.788	128778474.809	-3637754.681	-4056236.306	-10510950.355

59.280	312.237	312.237	0.000	0.000	0.000
312.237	59.280	312.237	0.000	0.000	0.000
312.237	312.237	59.280	0.000	0.000	0.000
-13061937.302	-13789029.537	-13642829.158	-518829.123	327421.989	-106029.858
-13878829.953	-12805869.118	-13120869.223	603915.464	-1133587.133	-64560.373
-13214798.062	-13629905.070	-13404079.528	480503.506	-131550.011	-702172.604

59.278	312.235	312.235	0.000	-0.000	0.000
312.235	59.278	312.235	0.000	-0.000	-0.000
312.235	312.235	59.278	-0.000	-0.000	0.000
12917513.829	13723308.841	14187314.727	-904065.145	188204.031	173812.534
13085747.890	12554590.809	13934478.086	86092.666	-1550430.799	715186.859
13116625.700	13427811.081	12949338.775	-380971.924	700274.748	-750593.574

1369.884	-242.672	-242.672	0.000	0.000	-0.000
-242.672	1369.884	-242.672	0.000	0.000	0.000
-242.672	-242.672	1369.884	0.000	0.000	-0.000
-0.000	0.000	0.000	-244.297	-0.000	0.000
0.000	-0.000	0.000	0.000	-244.297	0.000
0.000	0.000	-0.000	0.000	-0.000	-244.297

1369.884	-242.672	-242.672	-0.000	0.000	0.000
-242.672	1369.884	-242.672	-0.000	-0.000	-0.000
-242.672	-242.672	1369.884	-0.000	0.000	-0.000
0.000	-0.000	-0.000	-244.297	-0.000	0.000
-0.000	0.000	-0.000	-0.000	-244.297	0.000
-0.000	0.000	0.000	-0.000	-0.000	-244.297

1369.884	-242.672	-242.672	0.000	0.000	0.000
-242.672	1369.884	-242.672	-0.000	0.000	-0.000
-242.672	-242.672	1369.884	-0.000	0.000	-0.000
-0.000	0.000	0.000	-244.297	-0.000	-0.000
0.000	-0.000	0.000	0.000	-244.297	-0.000
0.000	0.000	-0.000	0.000	-0.000	-244.297

1369.883	-242.672	-242.672	-0.000	0.000	0.000
-242.672	1369.883	-242.672	-0.000	-0.000	-0.000
-242.672	-242.672	1369.883	-0.000	0.000	0.000
0.000	-0.000	-0.000	-244.297	0.000	0.000
-0.000	0.000	-0.000	-0.000	-244.297	0.000
-0.000	-0.000	0.000	-0.000	0.000	-244.297

1369.891	-242.671	-242.671	0.000	0.000	0.000
-242.671	1369.891	-242.671	0.000	0.000	0.000
-242.671	-242.671	1369.891	0.000	0.000	-0.000
-0.002	0.001	0.000	-244.297	-0.000	-0.000
0.001	-0.002	0.000	0.000	-244.297	0.000
0.001	0.000	-0.002	0.000	0.000	-244.297

1369.876	-242.674	-242.674	-0.000	0.000	0.000
-242.674	1369.876	-242.674	-0.000	-0.000	-0.000
-242.674	-242.674	1369.876	0.000	-0.000	0.000
0.002	-0.001	-0.000	-244.297	-0.000	0.000
-0.001	0.002	-0.000	-0.000	-244.297	0.000
-0.001	-0.000	0.002	0.000	-0.000	-244.297

164.884	297.323	307.378	0.000	-0.000	-0.000
316.731	145.476	307.378	-0.000	0.000	0.000
307.378	307.378	139.235	-0.000	0.000	0.000
-0.000	0.000	-0.000	-55.161	0.000	0.000
-0.000	-0.000	-0.000	-0.000	-55.161	-0.000
-1626792310.512	-1279890875.863	-1116922255.095	-33381828.882	18444204.366	-100539365.732

145.476	316.731	307.378	0.000	-0.000	-0.000
297.323	164.884	307.378	-0.000	0.000	-0.000
307.378	307.378	139.234	0.000	-0.000	-0.000
-0.000	0.000	0.000	-55.161	-0.000	-0.000
0.000	0.000	0.000	0.000	-55.161	-0.000
1545014351.474	1349603918.242	1108620734.279	-24216436.709	4704157.880	-93760861.125

-149618434.729	-126368587.007	-102455697.989	3906112.397	582006.924	2161567.716
297.323	164.884	307.378	-0.000	-0.000	0.000
307.378	307.378	139.235	-0.000	0.000	0.000
0.000	0.000	-0.000	-55.161	0.000	0.000
0.000	0.000	-0.000	-0.000	-55.161	0.000
0.000	0.000	0.000	-0.000	0.000	-85.627

168549915.293	169254120.802	105618549.347	-1023703.988	1622467.712	-7656956.165
316.731	145.476	307.378	0.000	-0.000	-0.000
307.378	307.378	139.234	-0.000	-0.000	-0.000
-0.000	-0.000	0.000	-55.161	-0.000	-0.000
-0.000	-0.000	0.000	-0.000	-55.161	-0.000
-0.000	-0.000	-0.000	-0.000	0.000	-85.627

164.885	297.324	307.379	0.000	-0.000	0.000
297.324	164.885	307.379	0.000	-0.000	0.000
307.379	307.379	139.236	0.000	0.000	0.000
0.000	0.000	-0.000	-55.161	-0.000	0.000
0.000	0.000	-0.000	0.000	-55.161	0.000
-17624545.491	-14088019.184	-10881392.305	-347174.722	121388.628	-616664.085

16867974.506	17628909.976	10526558.696	34802.041	460523.973	-481116.897
297.322	164.882	307.377	0.000	0.000	-0.000
307.377	307.377	139.233	0.000	-0.000	-0.000
-0.000	-0.000	0.000	-55.161	-0.000	-0.000
-0.000	-0.000	0.000	-0.000	-55.161	0.000
17231974.705	14392108.177	10951784.704	56510.089	115117.092	-1272590.446

581.389	380.434	-0.637	-0.000	-0.000	-0.000
429.073	532.750	-0.637	-0.000	-0.000	-0.000
-0.637	-0.637	1064.200	-0.000	0.000	-0.000
0.000	-0.000	-0.000	-34.025	-0.000	-0.000
0.000	0.000	-0.000	0.000	-34.025	-0.000
0.000	0.000	-0.000	0.000	-0.000	51.839

581.389	380.434	-0.637	0.000	0.000	0.000
429.073	532.750	-0.637	-0.000	-0.000	0.000
-0.637	-0.637	1064.200	-0.000	0.000	0.000
0.000	-0.000	0.000	-34.025	0.000	0.000
-0.000	-0.000	0.000	-0.000	-34.025	0.000
-0.000	-0.000	0.000	-0.000	-0.000	51.839

532.751	429.072	-0.637	-0.000	0.000	-0.000
429.073	532.751	-0.637	-0.000	0.000	-0.000
-0.637	-0.637	1064.201	0.000	0.000	0.000
0.000	0.000	-0.000	-34.025	-0.000	-0.000
0.000	0.000	-0.000	0.000	-34.025	-0.000
0.000	0.000	-0.000	-0.000	0.000	51.839

532.749	429.073	-0.637	-0.000	-0.000	0.000
429.072	532.750	-0.637	-0.000	-0.000	0.000
-0.637	-0.637	1064.199	-0.000	-0.000	0.000
-0.000	-0.000	0.000	-34.025	-0.000	0.000
-0.000	-0.000	0.000	0.000	-34.025	0.000
-0.000	-0.000	0.000	0.000	-0.000	51.839

532.762	429.069	-0.638	0.000	0.000	-0.000
429.073	532.758	-0.638	-0.000	0.000	-0.000
-0.639	-0.639	1064.209	-0.000	0.000	-0.000
0.000	0.001	-0.000	-34.025	-0.000	-0.000
0.001	0.000	-0.000	0.000	-34.025	-0.000
0.000	0.001	-0.000	-0.000	-0.000	100.478

532.739	429.076	-0.637	0.000	-0.000	0.000
429.072	532.743	-0.637	0.000	0.000	0.000
-0.636	-0.636	1064.191	0.000	-0.000	0.000
-0.000	-0.001	0.000	-34.025	-0.000	0.000
-0.001	-0.000	0.000	0.000	-34.025	0.000
-0.000	-0.001	0.000	0.000	0.000	100.478

247.581	358.616	-110.939	-0.000	-0.000	-0.000
364.969	241.228	-110.939	0.000	-0.000	-0.000
-110.938	-110.938	107.855	0.000	-0.000	0.000
0.000	0.000	-0.000	134.909	-0.000	0.000
0.000	0.000	-0.000	0.000	134.909	0.000
0.000	-0.000	-0.000	0.000	-0.000	-61.871

241.228	364.969	-110.939	-0.000	0.000	0.000
364.969	241.228	-110.939	-0.000	0.000	0.000
-110.938	-110.938	107.855	-0.000	-0.000	0.000
-0.000	-0.000	0.000	134.909	0.000	0.000
-0.000	-0.000	0.000	0.000	134.909	0.000
-0.000	0.000	-0.000	-0.000	-0.000	-61.871

241.228	364.969	-110.939	0.000	0.000	-0.000
364.969	241.228	-110.939	0.000	-0.000	-0.000
-110.938	-110.938	107.855	0.000	-0.000	0.000
0.000	0.000	-0.000	134.909	-0.000	-0.000
0.000	0.000	-0.000	0.000	134.909	-0.000
0.000	0.000	-0.000	0.000	-0.000	-61.871

241.227	364.969	-110.939	-0.000	0.000	0.000
364.969	241.227	-110.939	-0.000	-0.000	0.000
-110.938	-110.938	107.855	-0.000	0.000	0.000
-0.000	-0.000	0.000	134.909	-0.000	0.000
-0.000	-0.000	0.000	-0.000	134.909	-0.000
-0.000	-0.000	0.000	0.000	0.000	-61.871

247.582	358.617	-110.939	0.000	-0.000	-0.000
358.617	247.582	-110.939	-0.000	-0.000	-0.000
-110.938	-110.938	107.855	0.000	-0.000	-0.000
0.001	0.001	-0.001	134.909	0.000	-0.000
0.001	0.001	-0.001	0.000	134.909	-0.000
0.000	0.000	-0.000	0.000	-0.000	-61.871

241.226	364.968	-110.939	0.000	-0.000	0.000
364.968	241.226	-110.939	-0.000	0.000	0.000
-110.938	-110.938	107.854	-0.000	-0.000	-0.000
-0.001	-0.001	0.001	134.909	-0.000	-0.000
-0.001	-0.001	0.001	-0.000	134.909	0.000
-0.000	-0.000	0.000	-0.000	0.000	-61.871

122.514	345.614	306.855	0.000	0.000	0.000
345.614	122.514	306.855	0.000	0.000	0.000
306.855	306.855	130.402	0.000	0.000	0.000
-0.000	0.000	-0.000	-63.777	0.000	0.000
-0.000	0.000	-0.000	0.000	-63.777	0.000
0.000	-0.000	0.000	0.000	0.000	-69.305

164.759	303.369	306.855	-0.000	-0.000	0.000
303.369	164.759	306.855	-0.000	-0.000	0.000
306.855	306.855	130.402	-0.000	-0.000	-0.000
0.000	-0.000	0.000	-63.777	-0.000	-0.000
0.000	-0.000	0.000	-0.000	-63.777	-0.000
-0.000	-0.000	-0.000	0.000	0.000	-69.305

164.759	303.369	306.855	0.000	0.000	0.000
-155240118.572	-175738106.020	-98499459.119	5285842.209	212949.584	5099030.813
306.855	306.855	130.402	0.000	0.000	0.000
0.000	-0.000	-0.000	-63.777	0.000	0.000
0.000	0.000	-0.000	-0.000	-63.777	-0.000
0.000	-0.000	0.000	0.000	0.000	-69.305

164.759	303.369	306.855	-0.000	-0.000	0.000
303.369	164.759	306.855	0.000	-0.000	0.000
306.855	306.855	130.402	0.000	-0.000	0.000
-0.000	0.000	0.000	-63.777	-0.000	-0.000
-0.000	-0.000	0.000	0.000	-63.777	-0.000
-0.000	0.000	-0.000	-0.000	-0.000	-69.305

164.760	303.370	306.856	0.000	0.000	-0.000
-15064762.897	-16756276.849	-9555221.100	-117721.694	274271.663	435233.138
306.856	306.856	130.403	0.000	0.000	0.000
0.000	-0.000	-0.000	-63.777	0.000	0.000
0.000	0.000	-0.000	-0.000	-63.777	0.000
-16636978.215	-16412866.343	-9702057.759	140588.189	-188320.203	357201.590

15494015.074	16181217.957	9403332.644	-523432.005	-5503.738	50234.556
14796148.359	11821038.831	10027502.307	-92375.356	-752091.677	237183.251
306.854	306.854	130.401	-0.000	-0.000	-0.000
-0.000	-0.000	0.000	-63.777	-0.000	0.000
-0.000	-0.000	0.000	0.000	-63.777	0.000
16209582.818	16379865.739	9174375.413	48247.361	451204.133	-1378740.512

1079.358	127.240	127.240	-0.000	-0.000	-0.000
127.240	1079.358	127.240	-0.000	-0.000	-0.000
127.240	127.240	1079.358	-0.000	-0.000	-0.000
0.000	0.000	0.000	622.742	-0.000	-0.000
0.000	0.000	0.000	-0.000	622.742	-0.000
0.000	0.000	0.000	-0.000	-0.000	622.742

Phonon and Quasi-Harmonic Approximation Predictions

Phonon band structures and crystal properties estimated from quasi-harmonic approximation (QHA) calculations are displayed for select crystals. The calculations were performed using phonopy and LAMMPS. For the phonon calculations, 3x3 supercells of the potential-specific relaxed crystals were used. The QHA calculations were based on 11 strain states ranging from -0.05 to 0.05.

The calculation method used is available as the iprPy phonon calculation method.

Notes and Disclaimers:

The thermodynamic properties estimated from QHA are based on the assumption that only volumetric changes affect the phonon behaviors as the temperature changes. This tends to give good predictions at lower temperatures but ignores anharmonic effects such as phonon coupling and vacancy formation that can be important at higher temperatures.
Note that direct molecular dynamics (MD) simulations using the same potentials will disagree with the thermodynamics properties listed here for the lowest temperatures. The MD results are purely classical in nature and therefore lack a zero-point energy, whereas the phonon calculations inherently provide a zero-point energy.
The structures explored here are taken from the dynamically relaxed structures above. Despite the rigorous relaxation method used, some of these structures prove to be unstable once internal deformations are added. The phonon results may reflect this and give bad band gap predictions for these unstable crystals.
All QHA calculations performed here use the same set of strains which might not be ideal for all crystals. Be sure to check the bulk modulus and Helmholtz vs. volume plots to verify if the QHA strained states are reasonable for each crystal of interest.
QHA results may not be available for all crystal structures that have phonon results. Missing QHA results indicates an issue either with the strained states or with the QHA calculation itself.
NIST disclaimer

Version Information:

2023-03-14. Phonon and QHA plots added

Band Structures, Density of States, and QHA Verification Plots

Composition:

Prototype:

a₀:

plot:

Thermodynamic Predictions

Composition:

Plot:

Download data

Click on plot to load interactive version

2010--Kim-H-K--Nb-C--LAMMPS--ipr1/phonon.C.G.png

Free Surface Formation Energy Predictions

Static free surface formation energies are displayed for select crystals. The values displayed here are obtained by taking a perfect periodic bulk crystal, slicing along a crystallographic plane, and using force minimization to statically relax the surfaces. The free surface formation energy is computed by comparing the energy of the defect system to the bulk system and dividing by the total surface area created by the cut.

The calculation method used is available as the iprPy surface_energy_static calculation method.

Notes and Disclaimers:

These values are meant to be guidelines for comparing potentials, not the absolute values for any potential's properties. Values listed here may change if the calculation methods are updated due to improvements/corrections. Variations in the values may occur for variations in calculation methods, simulation software and implementations of the interatomic potentials.
The calculation only performs straight cuts along crystallographic planes and static relaxations. Lower energy configurations may exist that require atomic restructuring of the surfaces.
Multiple values may be listed for a given plane followed by a number indicating different unique atomic planar cuts for the same theoretical plane.
NIST disclaimer

Version Information:

2019-11-14. Calculations for the alternate #2 plane slices are now different from the #1 plane slices.
2019-08-07. Data added.

Composition:

Prototype:

a₀:

Download raw data

Surface	γ_fs (mJ/m²)
(111) 2	2069.22
(331) 2	2673.97
(110)	3843.80
(321)	3915.92
(332)	4172.62
(211)	4242.51
(320)	4719.00
(311) 2	4774.93
(322)	4789.84
(210)	5255.62
(311) 1	5351.34
(221)	5465.99
(310)	5629.16
(331) 1	5891.59
(111) 1	6320.92
(100)	7123.56

Point Defect Formation Energy Predictions

Static point defect formation energies, E_f, and elastic dipole tensors, p_ij, are displayed for select crystals. Relaxed defect configurations are obtained by taking a 12x12x12 supercell of a perfect periodic bulk crystal, inserting the point defect, and using force minimization to statically relax the atomic positions while keeping the system dimensions constant. E_f is computed by comparing the energy of the defect system to the same number of atoms in a perfect bulk crystal. p_ij is estimated as the difference in the system's global pressure with and without the defect multiplied by the system's volume.

Simple structural comparisons of the unrelaxed and relaxed defect configurations are used to help determine if the defect structure has relaxed to a different configuration. Relaxed structures that are identified as no longer consistent with the ideal defect definition are excluded from the table below. The only exception to this is if the lowest energy interstitial configuration does not coincide with a known ideal defect, its energies and pressure tensor are included under the listing "relaxed interstitial". The full list of calculation results including the transformed structures and the structural comparison values is included in the csv file available from the "Download raw data" link.

The calculation method used is available as the iprPy point_defect_static calculation method.

Notes and Disclaimers:

These values are meant to be guidelines for comparing potentials, not the absolute values for any potential's properties. Values listed here may change if the calculation methods are updated due to improvements/corrections. Variations in the values may occur for changes in calculation methods, simulation software and implementations of the interatomic potentials.
The computed formation energy and elastic dipole tensor values are sensitive to the size of the system used. The 12x12x12 supercell size was selected as it should provide a decent approximation of the true bulk values.
The tests for identifying configuration relaxations are not guaranteed to be comprehensive or suitable for all point defect types. The table only shows the point defects that are consistent with the ideal configurations for that defect type, and the lowest energy interstitial.
The "relaxed interstitial" label indicates that the structure is not consistent with any of the known ideal interstitial configurations. No information is provided as to what that relaxed strucutre is, whether it is an unknown structure, a "near"-ideal defect, or a collapse to amorphous.
NIST disclaimer

Composition:

Prototype:

a₀:

Download raw data

Point Defect	E_f (eV)	p₁₁ (eV)	p₂₂ (eV)	p₃₃ (eV)	p₁₂ (eV)	p₁₃ (eV)	p₂₃ (eV)
vacancy	3.349	4.649	4.649	4.649	-0.000	-0.000	-0.000
1nn divacancy	4.856	4.777	4.777	4.777	-5.493	-5.493	-5.493
100 dumbbell	9.383	17.141	17.141	74.886	-0.000	-0.000	-0.000

Dislocation Monopole Structures

Plots of the atomic structure of dislocation cores are shown here. The dislocations were constructed using dislocation monopole configurations in which a straight dislocation is created in the center of the system using anisotropic elasticity and the elastic constants computed above. A cylindrical region centered around the dislocation is allowed to relax while the remaining atoms are held fixed. The relaxation consists of 1 million NVT steps at 10 K followed by an energy/force minimization. The final plots show the atomic configuration near the center of the monopole system, and display differential displacement arrows and Nye tensor color maps that are both taken with respect to the Burgers vector direction.

These plots show a wide range of final structures, some corresponding to expected dislocation structures and others not. While the characterizations are dislocation-type dependent, some general classification guidelines are listed here.

Planar dislocations exist on a single slip plane, which results in a grouping of large differential displacement arrows between two atomic rows. All edge dislocations and some screw dislocations may be planar dislocations.
Non-planar dislocations are screw dislocations that are positioned at the intersection of more than one slip plane. This results in the differential displacement arrows circling the position of the dislocation core. These core structures often have rotational symmetry due to the intersecting planes.
The plotted atomic positions are also useful for characterizing the final configuration. For screw dislocations, the atomic projections should match projections in the bulk crystal if the crystal has not transformed. For edge dislocations, the projections should be the bulk projections with the extra half-plane of atoms directly above the dislocation core.
Images in which the atomic positions look correct but no large differential displacement arrows are observed are likely the result of the dislocation gliding out of view during relaxation.
Phase transformations are identified by large differential displacement arrows spread out in both plot directions, often across the whole image. Additionally, the atomic position projections will likely not correspond to what is expected.
Twins typically appear as a long line of large differential displacement arrows along a single plane. Note that twins can nucleate phase transformations so the distinction between the two is not always clear.

The calculation method used is available as the iprPy dislocation_monopole calculation method.

Links are provided above the images allowing for users to download the final defect configuration and a corresponding defect-free configuration if they wish to further analyze the structures or use them for subsequent simulations.

Notes and Disclaimers:

While the general dislocation monopole method is universal to all dislocations, these calculations rely on dislocation-type specific atomic configuration sizes and differential displacement/Nye tensor settings. For the system size, the specific super cell multipliers vary based on the crystal and crystal vector orientations, but the sizes are always selected to minimize the length along the dislocation line and attempt to make the two other box dimensions similar. For the structure analysis tools, both are sensitive to the cutoff distance used to identify neighbors, and both have some method-specific parameters. Where known, the settings used here correspond to the most common ones used in literature.
For bcc dislocations, the differential displacement vectors shown here are "wrapped" in that if their magnitudes exceed b/2 then a value of b is added/subtracted to wrap the vector back around. This convention reveals the dislocation's structure and location without showing the slip history.
NIST disclaimer

Version Information:

2025-08-07. Plots first added.

Composition:

Type:

Nye tensor + differential displacement plot

Reference system: base.dump
Dislocation system: disl.dump

MD Solid Property Predictions

Plots of lattice and elastic constants are shown as a function of temperature. The 0K points were taken from the Crystal Structure Predictions and the Elastic Constants Predictions sections above for the unique crystal structures relaxed with the "dynamic" method. Starting from the 0 K relaxed crystal unit cells, supercell systems are created by replicating all three dimensions by the same multiplier to achieve at least 4000 atoms. The systems are then relaxed at 50 K and zero pressure using 1 million NPT steps. Lattice constants are estimated by averaging the measured box dimensions. Temperatures are iteratively increased by 50 K, with each subsequent relaxation calculation starting from the final atomic configuration at the previous temperature and relaxing for another 1 million steps.

The elastic constants are calculated using the deformation–fluctuation hybrid method. Starting from the final atomic configurations of the dynamic relaxations, the system is allowed to evolve at constant volume with a Langevin thermostat. The Born matrix is computed during this run by evaluating how the atomic forces would vary due to applied linear strain fields. The elastic constants can then be estimated using the averaged Born matrix values and the averaged stresses on the system.

The calculation methods used are available as the iprPy relax_dynamic and elastic_constants_dynamic calculation methods.

Clicking on the image of a plot will open an interactive version of it in a new tab. The underlying data for the plots can be downloaded by clicking on the links above each plot.

Notes and Disclaimers:

The maximum temperature shown for a crystal either corresponds to the maximum value that has been computed so far, or the maximum value that the structure is believed to have remained untransformed. The unstable/transformation temperature identification is mostly automated and can miss transitions not associated with large discontinuities in property values with changing temperature.
For succinctness, the elastic constants displayed here are averaged according to the 0 K structure's crystal symmetry family. If the structures deviate from the expected symmetry at higher temperatures, the values may not be valid. Raw results are available for verification if requested.
NIST disclaimer

Version Information:

2025-08-07. Plots added.

Composition:

Prototype:

a₀:

Download data

Click on plot to load interactive version

2010--Kim-H-K--Nb-C--LAMMPS--ipr1/mdsolid.C.A4--C--dc.7d9409b9.a

Click on plot to load interactive version

2010--Kim-H-K--Nb-C--LAMMPS--ipr1//mdsolid.C.A4--C--dc.7d9409b9.Cij

MD Thermodynamic Predictions

Plots of internal energy, Gibbs free energy, entropy, heat capacity and volume are shown here as a function of temperature for various crystal structures and liquid phases. The included crystal structures correspond to those in the Solid Structures vs. Temperature section, and the liquid phases to those in the Liquid Properties section. Internal energy and volume are taken from the associated structure relaxations mentioned in those previous sections. Constant pressure heat capacity is estimated using 3-point numerical derivatives of enthalpy versus temperature. Note that since all simulations done here are at 0 pressure, internal energy and enthalpy are equivalent.

The Gibbs energies of the phases are estimated using thermodynamic integration between the interatomic potential in question and a simpler model with known Gibbs free energy values. For solids, the reference model is an Einstein solid, while for liquids it is the Uhlenbeck-Ford potential. Besides a short run at the start of the solid calculations to estimate Einstein model spring constants, the two calculations proceed similarly. Starting with the final relaxation configurations, the systems are stabilized for 25,000 steps. Then, over the next 50,000 steps the potential is swapped out for the reference potential. The system is then stabilized for another 25,000 steps with the reference model before a reverse swap of 50,000 steps is performed. The simulation ends with one final 25,000 step stabilization period. From the simulation, the irreversible work of transformation is estimated and used to compute the absolute Gibbs free energy of the target phase and potential relative to the reference potential. Entropy is estimated as the difference in enthalpy and Gibbs free energy and divided by temperature.

The calculation methods used are available as the iprPy relax_dynamic, relax_liquid_redo, free_energy, and free_energy_liquid calculation methods.

Clicking on the image of a plot will open an interactive version of it in a new tab. The underlying data for the plots can be downloaded by clicking on the links above each plot.

Notes and Disclaimers:

Good free energy estimates require that the initial and final phases before and after the thermodynamic integration remain the same. This is not an issue with the liquid calculations as the Uhlenbeck-Ford model always predicts a liquid, but may be an issue for the solid phase estimates. This is likely to be an issue for non-compact phases which the Einstein model is poor at representing.
NIST disclaimer

Version Information:

2025-08-07. Plots added.

Composition:

Property:

Download data

Click on plot to load interactive version

2010--Kim-H-K--Nb-C--LAMMPS--ipr1/mdthermo.C.U

1079.358	127.239	127.239	0.000	0.000	0.000
127.239	1079.358	127.239	0.000	0.000	0.000
127.239	127.239	1079.358	0.000	0.000	0.000
-0.000	-0.000	-0.000	622.742	0.000	0.000
-0.000	-0.000	-0.000	0.000	622.742	0.000
-0.000	-0.000	-0.000	0.000	0.000	622.742

1079.361	127.240	127.240	-0.000	-0.000	-0.000
127.240	1079.361	127.240	-0.000	-0.000	-0.000
127.240	127.240	1079.361	-0.000	-0.000	-0.000
0.000	0.001	0.001	622.742	-0.000	-0.000
0.001	0.000	0.001	-0.000	622.742	-0.000
0.001	0.001	0.000	-0.000	-0.000	622.742

1079.355	127.239	127.239	0.000	0.000	0.000
127.239	1079.355	127.239	0.000	0.000	0.000
127.239	127.239	1079.355	0.000	0.000	0.000
-0.000	-0.001	-0.001	622.742	0.000	0.000
-0.001	-0.000	-0.001	0.000	622.742	0.000
-0.001	-0.001	-0.000	0.000	0.000	622.742

463.488	342.745	269.870	-0.000	0.000	0.000
342.744	463.488	269.870	-0.000	0.000	0.000
269.868	269.872	496.197	-0.000	0.000	-0.000
-0.007	0.007	0.000	-1964.032	0.000	0.000
-0.006	0.006	0.000	-0.000	-1964.032	0.000
-0.002	0.002	0.000	-0.000	0.000	114.030

463.492	342.746	269.871	-0.000	0.000	0.000
342.746	463.492	269.871	-0.000	0.000	0.000
269.871	269.871	496.200	-0.000	0.000	0.000
-0.029	0.058	0.008	-1964.031	0.000	-0.000
0.058	-0.029	0.008	-0.000	-1964.031	-0.000
-0.000	0.000	-0.000	-0.000	0.000	114.030

463.484	342.743	269.869	0.000	-0.000	0.000
342.744	463.484	269.869	0.000	0.000	-0.000
269.869	269.869	496.193	0.000	-0.000	-0.000
0.029	-0.058	-0.008	-1964.031	-0.000	-0.000
-0.058	0.029	-0.008	0.000	-1964.031	-0.000
0.000	-0.000	0.000	0.000	-0.000	114.030

452.277	304.229	-72.915	-0.000	-0.000	-0.000
304.229	452.277	-72.915	-0.000	-0.000	0.000
-72.915	-72.915	126.682	-0.000	-0.000	0.000
0.000	0.000	-0.000	327.873	-0.000	-0.000
0.000	0.000	-0.000	-0.000	327.873	-0.000
-0.000	0.000	0.000	-0.000	-0.000	97.018

452.276	304.229	-72.915	0.000	0.000	0.000
304.229	452.276	-72.915	0.000	0.000	0.000
-72.915	-72.915	126.682	0.000	0.000	-0.000
-0.000	-0.000	0.000	327.873	0.000	0.000
-0.000	-0.000	0.000	0.000	327.873	0.000
-0.000	-0.000	0.000	0.000	0.000	97.018

452.281	304.230	-72.915	-0.000	-0.000	0.000
304.230	452.281	-72.915	-0.000	-0.000	0.000
-72.915	-72.915	126.682	-0.000	-0.000	0.000
0.002	0.002	-0.001	327.873	-0.000	-0.000
0.002	0.002	-0.001	-0.000	327.873	-0.000
0.000	0.000	-0.000	-0.000	-0.000	97.018

452.273	304.229	-72.915	-0.000	0.000	-0.000
304.229	452.273	-72.915	0.000	-0.000	0.000
-72.915	-72.915	126.681	0.000	0.000	0.000
-0.002	-0.002	0.001	327.873	0.000	0.000
-0.002	-0.002	0.001	0.000	327.873	0.000
-0.000	-0.000	0.000	0.000	0.000	97.018

975.560	9.514	9.225	22.441	0.000	0.000
9.514	975.560	9.225	-22.441	-0.000	0.000
9.857	9.857	51.050	0.000	0.000	0.000
20.347	-20.347	-0.000	283.269	-0.000	0.000
-0.005	-0.005	-0.007	-0.000	283.205	22.445
-0.000	0.000	0.000	-0.000	20.347	551.311

1044.297	-58.325	9.857	20.347	-0.000	-0.000
-58.325	1044.297	9.857	-20.347	-0.000	-0.000
9.265	9.265	50.215	0.000	-0.000	-0.000
22.450	-22.440	0.007	283.205	0.000	-0.000
0.005	0.005	0.007	0.000	283.205	22.445
0.000	-0.000	-0.000	-0.000	20.347	551.311

975.712	10.011	9.681	22.446	0.000	0.000
9.928	975.794	9.682	-22.443	0.000	0.000
9.247	9.247	50.189	-0.000	0.000	0.000
22.440	-22.441	-0.001	283.205	-0.000	0.000
-0.001	-0.001	-0.001	-0.000	283.205	22.441
-0.000	-0.000	-0.001	0.000	22.441	483.023

1044.296	-58.325	9.857	20.347	-0.000	-0.000
9.532	975.579	9.251	-22.441	-0.000	-0.000
9.857	9.857	51.050	0.000	-0.000	-0.000
22.441	-22.440	0.001	283.205	-0.000	-0.000
0.001	0.001	0.001	0.000	283.205	22.441
0.000	0.000	0.001	-0.000	22.441	483.023

975.749	9.977	9.682	22.445	-0.000	0.000
9.539	975.571	9.248	-22.441	-0.000	0.000
9.249	9.249	50.193	-0.000	0.000	0.000
22.440	-22.441	0.000	283.208	0.000	0.000
-0.000	-0.000	0.000	-0.003	283.205	22.441
-0.000	0.001	-0.000	0.000	22.441	483.023

975.289	9.820	9.252	22.450	-0.000	-0.000
9.522	975.581	9.248	-22.440	0.000	-0.000
9.857	9.857	51.050	0.000	-0.000	-0.000
20.347	-20.347	-0.000	283.266	0.000	-0.000
0.000	0.000	-0.000	0.003	283.205	22.441
0.000	-0.001	0.000	-0.000	22.441	483.023

1052.880	-58.542	11.147	20.982	-0.000	-0.000
-58.542	1052.880	11.147	-20.982	-0.000	0.000
10.380	10.380	57.124	-0.000	-0.000	0.000
20.982	-20.982	-0.000	279.173	-0.000	0.000
0.009	0.008	0.010	-0.000	279.117	22.958
0.009	0.008	0.010	-0.000	22.944	487.346

984.033	9.341	10.569	22.944	0.000	-0.000
9.791	984.276	10.984	-22.947	0.000	-0.000
9.807	9.810	56.439	0.000	0.000	-0.000
22.949	-22.966	-0.010	279.117	0.000	-0.000
-0.009	-0.008	-0.010	0.000	279.117	22.958
-0.009	-0.008	-0.010	0.000	22.944	487.346

984.300	9.786	10.995	22.947	-0.000	0.000
9.375	984.066	10.608	-22.944	-0.000	0.000
10.380	10.380	57.124	-0.000	-0.000	0.000
20.982	-20.982	0.000	279.174	-0.000	0.000
-0.000	-0.000	0.000	-0.000	279.173	20.982
-0.000	0.000	0.000	0.000	20.982	555.711

984.075	9.364	10.606	22.944	-0.000	0.000
9.795	984.292	10.994	-22.947	-0.000	0.000
9.839	9.839	56.475	-0.000	-0.000	0.000
22.944	-22.945	0.000	279.120	-0.000	0.000
-0.000	-0.000	0.000	-0.003	279.173	20.982
-0.000	0.001	-0.000	0.000	20.982	555.711

2010--Kim-H-K-Jung-W-S-Lee-B-J--Nb-C

C 7

C-Nb 1

Nb 5

Implementation Information

Diatom Energy vs. Interatomic Spacing

Cohesive Energy vs. Interatomic Spacing

Crystal Structure Predictions

Elastic Constants Predictions

Phonon and Quasi-Harmonic Approximation Predictions

Band Structures, Density of States, and QHA Verification Plots

Thermodynamic Predictions

Free Surface Formation Energy Predictions

Point Defect Formation Energy Predictions

Dislocation Monopole Structures

MD Solid Property Predictions

MD Thermodynamic Predictions

984.286	9.794	10.994	22.947	0.000	-0.000
9.362	984.068	10.604	-22.944	0.000	-0.000
9.838	9.838	56.473	0.000	0.000	-0.000
22.945	-22.944	-0.000	279.115	0.000	-0.000
0.000	0.000	-0.000	0.003	279.173	20.982
0.000	-0.001	0.000	-0.000	20.982	555.711

1055.911	-58.409	11.555	21.185	-0.000	-0.000
-58.409	1055.911	11.555	-21.185	-0.000	-0.000
10.625	10.625	58.971	0.000	-0.000	-0.000
23.088	-23.088	0.000	277.594	-0.000	0.001
0.000	-0.000	-0.000	-0.000	277.647	21.185
0.002	-0.001	0.000	-0.000	23.087	488.763

987.290	9.954	11.409	23.089	0.000	-0.001
9.527	987.054	11.034	-23.086	0.000	-0.001
10.625	10.625	58.971	0.000	-0.000	0.000
21.185	-21.185	0.000	277.647	-0.000	0.000
-0.000	0.000	0.000	0.000	277.647	21.185
-0.002	0.001	-0.000	0.000	23.087	488.762

987.291	9.954	11.409	23.089	0.000	-0.000
9.528	987.054	11.034	-23.087	0.000	-0.000
10.625	10.625	58.971	-0.000	0.000	0.000
21.185	-21.185	-0.000	277.647	0.000	0.000
-0.000	0.000	0.000	0.000	277.594	23.086
0.000	-0.000	0.000	-0.000	21.185	557.160

987.284	9.965	11.409	23.089	-0.000	0.000
9.537	987.048	11.034	-23.087	-0.000	0.000
10.104	10.104	58.382	-0.000	-0.000	0.000
23.086	-23.087	-0.000	277.597	-0.000	0.000
-0.000	-0.000	0.000	-0.003	277.647	21.185
-0.000	0.001	-0.000	0.000	23.087	488.762

987.281	9.960	11.409	23.089	0.000	-0.000
9.520	987.058	11.034	-23.086	0.000	-0.000
10.625	10.625	58.971	-0.000	-0.000	0.000
23.087	-23.086	0.000	277.592	0.000	-0.000
0.000	0.000	-0.000	0.003	277.647	21.185
0.000	-0.001	0.000	-0.000	23.087	488.762

392.769	-222.128	23.345	19.752	0.000	-0.000
-222.129	392.768	23.345	-19.752	0.000	0.000
-46.002	-46.002	83.044	-0.000	-0.000	0.000
57.810	-57.810	0.000	84.918	0.000	-0.000
0.000	0.000	0.000	0.000	84.918	57.810
0.000	0.000	0.000	0.000	57.810	201.076

392.768	-222.129	23.346	19.752	0.000	0.000
-128.939	273.213	122.163	-57.810	-0.000	-0.000
52.816	52.816	-657.660	0.000	-0.000	0.000
57.810	-57.810	-0.000	84.918	-0.000	0.000
0.000	-0.000	-0.000	-0.000	84.918	57.810
0.000	-0.000	-0.000	-0.000	57.810	201.076

279.033	-117.357	56.939	58.840	0.000	0.000
-117.358	279.034	56.939	-58.840	-0.000	0.000
52.816	52.816	-657.656	-0.000	0.000	0.000
57.810	-57.810	-0.000	84.918	-0.000	-0.000
0.000	0.000	-0.000	0.000	84.918	57.810
-0.000	-0.000	-0.000	-0.000	57.810	201.076

279.044	-117.364	56.936	58.838	0.000	-0.000
-128.941	273.218	122.162	-57.810	-0.000	-0.000
52.822	52.822	-657.696	0.000	-0.000	-0.000
57.810	-57.810	0.000	84.918	0.000	0.000
0.000	0.000	0.000	0.000	98.534	19.752
0.001	0.001	-0.000	0.000	57.810	201.076