Job ID: {{ job_id }}
Composition: {{ composition }} | Potential: {{ potential }}
{% if tags %}{% for tag in tags %} {{ tag }} {% endfor %}
{% endif %}Glass structure after melt-quench procedure. Rotate and zoom with mouse. Hover over atoms for element info. Click on an atom to cut out its 5 Å neighbourhood.
{{ protocol_description }}
{% endif %} {% endif %}Pair distribution functions (RDF) computed up to rmax = 10 Å with 2000 bins. Bond cutoffs are automatically determined from the first RDF minimum (Savitzky-Golay smoothed). These cutoffs feed into coordination number distributions, Qn speciation & network connectivity, bond-angle distributions, and Guttman ring statistics (up to 40-membered rings). Click legend items to show/hide individual traces.
Viscosity computed via the Green-Kubo method by integrating off-diagonal stress
autocorrelation functions (Pxy, Pxz, Pyz).
{% if viscosity_params %}
Production runs of {{ "%.1f" | format(viscosity_params.production_ns) }} ns at
{{ viscosity_params.temperatures | join(', ') }} K
(max correlation lag {{ "%.1f" | format(viscosity_params.max_lag_ns) }} ns).
{% endif %}
The autocorrelation integration cutoff is determined automatically via a noise-threshold
criterion.
{% if vft %}
Dashed line: VFT extrapolation log10(η) = A + B / (T − T0)
with A = {{ "%.2f" | format(vft.A) }}, B = {{ "%.0f" | format(vft.B) }} K,
T0 = {{ "%.0f" | format(vft.T0) }} K ({{ "%.0f" | format(vft.T0 - 273.15) }} °C).
{% if vft.ref_temps %}
Reference points:
{% for label, temp_c in vft.ref_temps.items() %}
{{ label }} = {{ "%.0f" | format(temp_c) }} °C{{ "," if not loop.last else "." }}
{% endfor %}
{% endif %}
{% endif %}
{% if cte_plots.convergence %} Linear CTE (αL = αV / 3) from enthalpy–volume fluctuations via iterative NPT production runs {% if cte_params and cte_params.temperature %}at {{ cte_params.temperature | int }} K{% endif %} ({{ "%.2f" | format(cte_params.production_ns) }} ns each, up to {{ cte_params.max_production_runs }} runs). {% endif %} {% if cte_plots.volume_temperature %} {% if cte_plots.convergence %}Additionally, volume{% else %}Volume{% endif %} vs. temperature data collected from NPT simulations {% if cte_params and cte_params.temperatures %}at {{ cte_params.temperatures | join(', ') }} K{% endif %} ({{ "%.2f" | format(cte_params.production_ns) }} ns each), enabling a linear fit of volumetric thermal expansion. {% endif %}
Full 6×6 elastic stiffness tensor (Cij) computed via the central-difference stress-strain method. {% if elastic_params %} The simulation box is deformed by ±{{ "%.1f" | format(elastic_params.strain_pct) }}% strain along each of the three normal (xx, yy, zz) and three shear (yz, xz, xy) directions (12 deformations total, {{ elastic_params.temperature | int }} K, {{ "%.3f" | format(elastic_params.production_ns) }} ns production per direction). {% else %} The simulation box is deformed by ±0.1% strain along each of the three normal (xx, yy, zz) and three shear (yz, xz, xy) directions (12 deformations total). {% endif %} Stress tensors are averaged over the second half of each production run. Isotropic moduli (B, G, E, ν) obtained via Voigt-Reuss-Hill averaging.