Materials database: Part D - Lua material scripts

This page explains OghmaNano's optional Lua material scripts. Rather than storing a single fixed number for each property, a script lets you define electrical, optical, and thermal parameters as functions of temperature, position (x, y, z), composition, or photon density. It covers how a script is recognised, where scripts can live on disk, how the model chooses between them, how each function is structured, and how to bring a script's values into a simulation.

1. What is a Lua material script?

A Lua material script is a plain-text .lua file in which each physical quantity is written as a small function. Every function is given a state table (containing the local temperature, position, and photon density) and returns two values: the calculated quantity and an enabled flag. The enabled flag is central to how scripts are used — only functions that are enabled will override the value otherwise stored in the simulation. A typical function looks like this:

function material.Xi(state)
	-- Electron affinity, units: eV
	local enabled = true
	local value = 4.45
	return value, enabled
end

Any file that OghmaNano finds is treated as a material script if the table declaration local material = {} appears within the first 1000 bytes of the file. This means a licence header or a block of explanatory comments can precede the declaration — the marker simply has to appear near the top of the file. You can view and edit a material's script directly from the Code editor tab in the Material editor ??.

A material script can define as much or as little as you like. It might define only the band gap and electron affinity, or a full set of band-structure, transport, recombination, and thermal parameters. Functions you do not need can simply be left out, or returned with enabled = false so that the model falls back to the value stored in the standard editors.

2. Where scripts live on disk

A material script can live in one of two places, and this determines how it appears in OghmaNano:

• Inside a material folder — for example inside the Gallium arsenide or MAPI-CH3NH3PbI3 folder. A script placed here belongs to that material and is displayed within the material when you open it, alongside its n/k data and other files.
• Outside a material folder — a stand-alone .lua file placed elsewhere in the database. In this case it is shown as a .lua file in its own right within the Materials database browser, rather than being attached to a single material.

Figure ?? shows a script living inside a material folder. The MAPI-CH3NH3PbI3 perovskite folder contains several optical (n/k) datasets — Ball, Leguy, Leguy-hydrated, and Phillips — represented by atom icons, together with a single electrical script, CH3NH3PbI3.lua, represented by a code-file icon.

3. One electrical script, several optical datasets

It is common to have several sets of optical constants for one material but only one sensible set of electrical parameters. Different research groups measure n and k over different wavelength ranges, using different samples and different methods, so it is useful to keep each of those optical datasets separately and choose between them per simulation. The electrical parameters (band gap, affinity, mobilities, and so on), by contrast, are usually best kept as a single, consistent definition for the material. That is exactly the arrangement shown in ??: four n/k datasets, one electrical script.

When the material is used, OghmaNano tries to select the most appropriate file automatically. If a script is not present at the level of the specific material being simulated, the model will search backwards through the material folder structure to find a suitable script to use. This means you can attach one electrical script at a sensible level and have it apply to the datasets grouped beneath it, without having to duplicate the same electrical parameters into every optical variant.

4. Anatomy of a material script

The example below is a complete electrical material script for AlGaInAs — the lattice-matched InP-based quaternary used for telecom lasers and detectors. It is shown in full so you can see how a real, fully documented script is put together; each function is then broken out and explained in Section 5 below. Click to expand.

-- See end of file for copyright, licensing and documentation links.

local material = {}

-- ---------------------------------------------------------------------------
-- NOTES ON AlGaInAs (read before using):
--
-- (Al_x Ga_{1-x})_{0.47} In_{0.53} As, LATTICE-MATCHED to InP (a = 5.869 A). The
-- workhorse quaternary for 1.3-1.55 um telecom lasers, modulators and
-- photodetectors. In is fixed at 0.53 for the InP match, so the single design
-- knob is the Al fraction x on the group-III sublattice.
--
-- DIRECT gap for all x (no crossover in the useful range): ~0.75 eV at x=0
-- (Ga0.47In0.53As, "InGaAs", ~1.65 um) rising to ~1.46 eV at x=1
-- (Al0.48In0.52As, ~0.85 um). Compared with the older InGaAsP telecom
-- quaternary, AlGaInAs offers a larger conduction-band offset (better electron
-- confinement), which is why it is favoured for high-temperature laser
-- operation - a device advantage, captured here only via Xi/Eg.
--
-- The lattice constant is FIXED by the InP match (independent of x). Eg is
-- computed from x; the remaining parameters are given at a representative
-- x = 0.30 (Eg ~ 0.92 eV, ~1.34 um) with InGaAs/AlInAs end points noted.
--
-- REFERENCES
--
-- [1]  I. Vurgaftman, J. R. Meyer, L. R. Ram-Mohan, J. Appl. Phys. 89, 5815
--      (2001). [constituent band parameters, bowing]
-- [2]  S. Adachi, "Properties of Semiconductor Alloys", Wiley (2009); and
--      S. Adachi, "Physical Properties of III-V Semiconductor Compounds"
--      (1992). [Eg(x), transport, dielectric, thermal]
-- [3]  O. Madelung, "Semiconductors: Data Handbook", Springer (2004).
-- ---------------------------------------------------------------------------

function material.name()
	local enabled = true

	return "AlGaInAs", enabled
end


function material.description()
	local enabled = true

	return "(AlxGa1-x)0.47In0.53As on InP, x = Al ~ 0.30 (direct gap)", enabled
end


function material.formula()
	local enabled = true

	return "AlGaInAs", enabled
end


function material.Eg(state)
	-- Units: eV
	-- Refs: [1],[2]
	--
	-- Direct gap, lattice-matched to InP, in Al fraction x. 300 K relation
	-- [2]:  Eg(x) = 0.76 + 0.49*x + 0.20*x^2
	-- (0.76 eV at x=0 -> 1.45 eV at x=1). Temperature handled by a
	-- Varshni-type shift about 300 K (alpha = 4.0e-4 eV/K, beta = 200 K).
	-- x = 0.30 -> ~0.925 eV (~1.34 um) at 300 K. Change al_fraction here
	-- and in the functions noted in the header.

	local enabled = true
	local al_fraction = 0.30
	local x = al_fraction
	local T = state.T
	local eg_300 = 0.76 + 0.49*x + 0.20*x*x
	local tshift = 4.0e-4*(90000.0/500.0 - T*T/(T + 200.0))
	local value = eg_300 + tshift

	return value, enabled
end


function material.Xi(state)
	-- Electron affinity
	-- Units: eV
	-- Refs: [2],[3]
	--
	-- ~4.45 eV at x = 0.30. Falls with Al (InGaAs ~4.6 -> AlInAs ~4.1 eV);
	-- the resulting large conduction-band offset to InP / across the alloy
	-- is the key confinement advantage over InGaAsP.

	local enabled = true
	local value = 4.45

	return value, enabled
end


function material.Nc(state)
	-- Effective conduction-band density of states
	-- Units: m^-3
	-- Refs: [2],[3]
	--
	-- ~4e23 m^-3 at 300 K for x = 0.30 (small direct-valley mass, m_e* ~
	-- 0.05-0.07 m0). Endpoints: InGaAs ~2.1e23, AlInAs ~5.4e23 m^-3.

	local enabled = true
	local T = state.T
	local value = 4.0e23*(T/300.0)^1.5

	return value, enabled
end


function material.Nv(state)
	-- Effective valence-band density of states
	-- Units: m^-3
	-- Refs: [2],[3]
	--
	-- ~8e24 m^-3 at 300 K for x = 0.30 (InGaAs ~7.7e24, AlInAs ~9e24 m^-3).

	local enabled = true
	local T = state.T
	local value = 8.0e24*(T/300.0)^1.5

	return value, enabled
end


function material.mu_e(state)
	-- Low-field electron mobility
	-- Units: m^2 V^-1 s^-1
	-- Refs: [2]
	--
	-- ~0.3 m^2/V/s (3000 cm^2/V/s) at x = 0.30 - high (small mass) but
	-- alloy-scattering-reduced from InGaAs (~10000-13000 cm^2/V/s) toward
	-- AlInAs (~1000-4000 cm^2/V/s). (300/T)^1.5 approximate.

	local enabled = true
	local T = state.T
	local value = 0.30*(300.0/T)^1.5

	return value, enabled
end

function material.mue_x(state)
	return material.mu_e(state)
end

function material.mue_y(state)
	return material.mu_e(state)
end

function material.mue_z(state)
	return material.mu_e(state)
end


function material.mu_h(state)
	-- Low-field hole mobility
	-- Units: m^2 V^-1 s^-1
	-- Refs: [2]
	--
	-- ~0.02 m^2/V/s (200 cm^2/V/s) at x = 0.30 (InGaAs ~300, AlInAs ~100
	-- cm^2/V/s). (300/T)^1.5 approximate.

	local enabled = true
	local T = state.T
	local value = 0.02*(300.0/T)^1.5

	return value, enabled
end

function material.muh_x(state)
	return material.mu_h(state)
end

function material.muh_y(state)
	return material.mu_h(state)
end

function material.muh_z(state)
	return material.mu_h(state)
end


function material.epsilonr(state)
	-- Relative static permittivity
	-- Dimensionless
	-- Refs: [2],[3]
	--
	-- ~13.4 at x = 0.30 (InGaAs ~13.9 -> AlInAs ~12.5).

	local enabled = true
	local value = 13.4

	return value, enabled
end


function material.free_to_free_recombination(state)
	-- Radiative recombination coefficient
	-- Units: m^3 s^-1
	-- Refs: representative (see note)
	--
	-- Direct-gap value ~1e-10 cm^3/s = 1e-16 m^3/s. Representative; refine
	-- against measured lifetimes.

	local enabled = true
	local value = 1.0e-16

	return value, enabled
end


function material.auger_Cn(state)
	-- Electron Auger recombination coefficient
	-- Units: m^6 s^-1
	-- Refs: representative (see note)
	--
	-- ~2e-28 cm^6/s = 2e-40 m^6/s at the InGaAs-rich (narrow-gap) end.
	-- Auger is a major loss mechanism in 1.55 um AlGaInAs/InGaAsP lasers
	-- and grows sharply as x decreases (gap narrows) - scale accordingly.

	local enabled = true
	local value = 2.0e-40

	return value, enabled
end


function material.auger_Cp(state)
	-- Hole Auger recombination coefficient
	-- Units: m^6 s^-1
	-- Refs: representative (see note)
	--
	-- As auger_Cn: representative 2e-28 cm^6/s = 2e-40 m^6/s.

	local enabled = true
	local value = 2.0e-40

	return value, enabled
end


function material.ss_srh_trap_energy(state)
	-- SRH trap energy relative to the middle of the band gap.
	-- Units: eV
	--
	-- Positive values are above mid-gap (towards the conduction band).
	-- Negative values are below mid-gap (towards the valence band).
	--
	-- Defect-dependent. Lattice-matched InP-based material can be very high
	-- quality; set level and density from your own data. Mid-gap default.

	local enabled = true
	local value = 0.0

	return value, enabled
end


function material.ss_srh_Nt(state)
	-- SRH trap density
	-- Units: m^-3
	--
	-- Defect-dependent placeholder; set from the intended material
	-- quality.

	local enabled = true
	local value = 1.0e21

	return value, enabled
end


function material.ss_srh_sigma_n(state)
	-- Electron capture cross section
	-- Units: m^2
	--
	-- Defect-dependent placeholder; set from measurement.

	local enabled = true
	local value = 1.0e-19

	return value, enabled
end


function material.ss_srh_sigma_p(state)
	-- Hole capture cross section
	-- Units: m^2
	--
	-- Defect-dependent placeholder; set from measurement.

	local enabled = true
	local value = 1.0e-19

	return value, enabled
end


function material.thermal_kl(state)
	-- Thermal conductivity
	-- Units: W m^-1 K^-1
	-- Refs: [2]
	--
	-- ~4.5 W/m/K at x = 0.30 - LOW, from strong mass-disorder phonon
	-- scattering in this quaternary (well below InP's ~68 or InAs's ~27
	-- W/m/K). Weakly dependent in the alloy regime; (300/T)^0.4 approximate.
	-- This low kappa contributes to self-heating in InP-based lasers.

	local enabled = true
	local T = state.T
	local value = 4.5*(300.0/T)^0.4

	return value, enabled
end


function material.heat_capacity(state)
	-- Specific heat capacity
	-- Units: J kg^-1 K^-1
	-- Refs: [3]
	--
	-- ~320 J/kg/K near 300 K (interpolated across the constituents).

	local enabled = true
	local value = 320.0

	return value, enabled
end


function material.density(state)
	-- Mass density
	-- Units: kg m^-3
	-- Refs: [3]
	--
	-- ~5300 kg/m^3 at x = 0.30 (InGaAs ~5500 -> AlInAs ~4650 kg/m^3).

	local enabled = true
	local value = 5300.0

	return value, enabled
end


function material.lattice_constant(state)
	-- Cubic lattice constant
	-- Units: m
	-- Refs: [1],[3]
	--
	-- FIXED at the InP match, a = 5.8697 A, INDEPENDENT of x (Al substitutes
	-- for Ga at constant lattice constant along the InP-matched line).
	-- Linear thermal expansion ~4.6e-6 /K (InP-like).

	local enabled = true
	local T = state.T
	local a300 = 5.8697e-10
	local expansion = 4.6e-6
	local value = a300*(1.0 + expansion*(T - 300.0))

	return value, enabled
end


function material.thermal_tau_e(state)
	-- Electron energy relaxation time towards the lattice temperature
	-- Units: s
	--
	-- Value basis: III-V family estimate
	-- Confidence: Medium
	--
	-- Reference:
	-- https://www.mdpi.com/2673-3978/3/2/16
	--
	-- Comments:
	-- Representative III-V carrier-to-lattice relaxation time. GaAs-like values
	-- are typically sub-ps to ps and field dependent.

	local enabled = true
	local value = 5.000000e-13

	return value, enabled
end


function material.thermal_tau_h(state)
	-- Hole energy relaxation time towards the lattice temperature
	-- Units: s
	--
	-- Value basis: III-V family estimate
	-- Confidence: Medium
	--
	-- Reference:
	-- https://www.mdpi.com/2673-3978/3/2/16
	--
	-- Comments:
	-- Representative III-V carrier-to-lattice relaxation time. GaAs-like values
	-- are typically sub-ps to ps and field dependent.

	local enabled = true
	local value = 5.000000e-13

	return value, enabled
end


function material.print()
	local state = {
		T = 300.0,
		x = 0.0,
		y = 0.0,
		z = 0.0,
		photon_density = 0.0,
	}

	print(string.format("Material:               %s", material.name()))
	print(string.format("Description:            %s", material.description()))
	print(string.format("Formula:                %s", material.formula()))
	print(string.format("Temperature:            %.2f K", state.T))
	print(string.format("Position:               %.6e, %.6e, %.6e m", state.x, state.y, state.z))
	print(string.format("Photon density:         %.6e m^-3", state.photon_density))

	print(string.format("Band gap:               %.6f eV", material.Eg(state)))
	print(string.format("Electron affinity:      %.6f eV", material.Xi(state)))
	print(string.format("Electron mobility:      %.6e m^2/V/s", material.mu_e(state)))
	print(string.format("Hole mobility:          %.6e m^2/V/s", material.mu_h(state)))
	print(string.format("Nc:                     %.6e m^-3", material.Nc(state)))
	print(string.format("Nv:                     %.6e m^-3", material.Nv(state)))
	print(string.format("Relative permittivity:  %.6f", material.epsilonr(state)))

	print(string.format("Radiative coeff.:       %.6e m^3/s", material.free_to_free_recombination(state)))
	print(string.format("Electron Auger coeff.:  %.6e m^6/s", material.auger_Cn(state)))
	print(string.format("Hole Auger coeff.:      %.6e m^6/s", material.auger_Cp(state)))

	print(string.format("SRH trap energy:        %.6f eV", material.ss_srh_trap_energy(state)))
	print(string.format("SRH trap density:       %.6e m^-3", material.ss_srh_Nt(state)))
	print(string.format("SRH sigma n:            %.6e m^2", material.ss_srh_sigma_n(state)))
	print(string.format("SRH sigma p:            %.6e m^2", material.ss_srh_sigma_p(state)))

	print(string.format("Electron energy relax.:  %.6e s", material.thermal_tau_e(state)))
	print(string.format("Hole energy relax.:      %.6e s", material.thermal_tau_h(state)))

	print(string.format("Thermal conductivity:   %.6e W/m/K", material.thermal_kl(state)))
	print(string.format("Heat capacity:          %.6e J/kg/K", material.heat_capacity(state)))
	print(string.format("Mass density:           %.6e kg/m^3", material.density(state)))
end

return material

5. The functions in detail

Every function in a material script follows the same pattern: it receives the state table and returns a value together with an enabled flag. The subsections below document each group of functions from the AlGaInAs example.

5.1 The state table and the enabled flag

Each function is passed a state table describing the local conditions at the point in the device where the value is being requested. The fields are:

• state.T — local temperature (K).
• state.x, state.y, state.z — position within the device (m).
• state.photon_density — local photon density (m^-3).

Because the model calls your function with the local state, any parameter can be made a function of these quantities — this is the whole point of using a script rather than a fixed value. The second return value, enabled, controls whether the model actually uses the script for that parameter. If enabled is true, the script value is used; if it is false, the parameter is left to the value stored in the standard editors. This is what lets a script define a handful of properties while leaving the rest to the ordinary fixed-value machinery.

5.2 Identity: name, description, formula

These three functions return descriptive strings rather than physical values. They identify the material, give a short human-readable description, and give its chemical formula.

function material.name()
	local enabled = true
	return "AlGaInAs", enabled
end

function material.description()
	local enabled = true
	return "(AlxGa1-x)0.47In0.53As on InP, x = Al ~ 0.30 (direct gap)", enabled
end

function material.formula()
	local enabled = true
	return "AlGaInAs", enabled
end

5.3 Band structure: Eg and Xi

Eg returns the band gap (eV) and Xi the electron affinity (eV). Together these two quantities set the positions of the conduction and valence band edges, and they are used throughout the model (see Section 7). The AlGaInAs band gap is the most instructive function in the file: it computes the 300 K gap from the Al fraction x and then applies a Varshni-type temperature shift, so the returned value depends on state.T.

function material.Eg(state)
	-- Units: eV.  Direct gap, lattice-matched to InP, in Al fraction x.
	local enabled = true
	local al_fraction = 0.30
	local x = al_fraction
	local T = state.T
	local eg_300 = 0.76 + 0.49*x + 0.20*x*x      -- 300 K value from composition
	local tshift = 4.0e-4*(90000.0/500.0 - T*T/(T + 200.0))  -- Varshni-type shift
	local value = eg_300 + tshift
	return value, enabled
end

function material.Xi(state)
	-- Electron affinity, units: eV
	local enabled = true
	local value = 4.45
	return value, enabled
end

5.4 Density of states: Nc and Nv

Nc and Nv return the effective conduction- and valence-band densities of states (m^-3). Both are quoted at 300 K and scaled with temperature as (T/300)^1.5, the standard parabolic-band result.

function material.Nc(state)
	local enabled = true
	local T = state.T
	local value = 4.0e23*(T/300.0)^1.5
	return value, enabled
end

function material.Nv(state)
	local enabled = true
	local T = state.T
	local value = 8.0e24*(T/300.0)^1.5
	return value, enabled
end

5.5 Mobilities: mu_e, mu_h and directional aliases

mu_e and mu_h return the low-field electron and hole mobilities (m² V^-1 s^-1), scaled with temperature as (300/T)^1.5. OghmaNano can use direction-dependent mobilities, so the file also provides mue_x, mue_y, mue_z and muh_x, muh_y, muh_z. In this isotropic example each simply returns the scalar value, but you could give each direction its own expression for an anisotropic material.

function material.mu_e(state)
	-- Low-field electron mobility, units: m^2 V^-1 s^-1
	local enabled = true
	local T = state.T
	local value = 0.30*(300.0/T)^1.5
	return value, enabled
end

function material.mue_x(state) return material.mu_e(state) end
function material.mue_y(state) return material.mu_e(state) end
function material.mue_z(state) return material.mu_e(state) end

function material.mu_h(state)
	-- Low-field hole mobility, units: m^2 V^-1 s^-1
	local enabled = true
	local T = state.T
	local value = 0.02*(300.0/T)^1.5
	return value, enabled
end

function material.muh_x(state) return material.mu_h(state) end
function material.muh_y(state) return material.mu_h(state) end
function material.muh_z(state) return material.mu_h(state) end

5.6 Permittivity: epsilonr

epsilonr returns the relative static permittivity (dimensionless), which sets the electrostatics in the drift–diffusion solver.

function material.epsilonr(state)
	-- Relative static permittivity, dimensionless
	local enabled = true
	local value = 13.4
	return value, enabled
end

5.7 Recombination: radiative and Auger

free_to_free_recombination returns the band-to-band (radiative) recombination coefficient (m³ s^-1), while auger_Cn and auger_Cp return the electron and hole Auger coefficients (m⁶ s^-1). Auger recombination is a major loss channel in narrow-gap III–V devices, and its strength grows sharply as the gap narrows, so these are natural candidates for a composition- or temperature-dependent expression in a more advanced script.

function material.free_to_free_recombination(state)
	-- Radiative recombination coefficient, units: m^3 s^-1
	local enabled = true
	local value = 1.0e-16
	return value, enabled
end

function material.auger_Cn(state)
	-- Electron Auger coefficient, units: m^6 s^-1
	local enabled = true
	local value = 2.0e-40
	return value, enabled
end

function material.auger_Cp(state)
	-- Hole Auger coefficient, units: m^6 s^-1
	local enabled = true
	local value = 2.0e-40
	return value, enabled
end

5.8 SRH traps

The four Shockley–Read–Hall functions describe a single equilibrium trap level: ss_srh_trap_energy (trap energy relative to mid-gap, eV; positive is towards the conduction band), ss_srh_Nt (trap density, m^-3), and ss_srh_sigma_n / ss_srh_sigma_p (electron and hole capture cross-sections, m²). These are strongly defect- and sample-dependent, so the values in the example are sensible placeholders to be replaced with your own measured data.

function material.ss_srh_trap_energy(state)
	-- Trap energy relative to mid-gap, units: eV
	local enabled = true
	local value = 0.0
	return value, enabled
end

function material.ss_srh_Nt(state)
	-- SRH trap density, units: m^-3
	local enabled = true
	local value = 1.0e21
	return value, enabled
end

function material.ss_srh_sigma_n(state)
	-- Electron capture cross section, units: m^2
	local enabled = true
	local value = 1.0e-19
	return value, enabled
end

function material.ss_srh_sigma_p(state)
	-- Hole capture cross section, units: m^2
	local enabled = true
	local value = 1.0e-19
	return value, enabled
end

5.9 Thermal and structural parameters

The remaining physical functions describe the material's thermal and structural behaviour: thermal_kl (lattice thermal conductivity, W m^-1 K^-1), heat_capacity (J kg^-1 K^-1), density (kg m^-3), lattice_constant (m, with linear thermal expansion applied about 300 K), and the carrier energy-relaxation times thermal_tau_e / thermal_tau_h (s). Note how thermal_kl and lattice_constant read state.T to return temperature-dependent values.

function material.thermal_kl(state)
	-- Thermal conductivity, units: W m^-1 K^-1
	local enabled = true
	local T = state.T
	local value = 4.5*(300.0/T)^0.4
	return value, enabled
end

function material.heat_capacity(state)
	-- Specific heat capacity, units: J kg^-1 K^-1
	local enabled = true
	local value = 320.0
	return value, enabled
end

function material.density(state)
	-- Mass density, units: kg m^-3
	local enabled = true
	local value = 5300.0
	return value, enabled
end

function material.lattice_constant(state)
	-- Cubic lattice constant, units: m (with thermal expansion)
	local enabled = true
	local T = state.T
	local a300 = 5.8697e-10
	local expansion = 4.6e-6
	local value = a300*(1.0 + expansion*(T - 300.0))
	return value, enabled
end

function material.thermal_tau_e(state)
	-- Electron energy relaxation time, units: s
	local enabled = true
	local value = 5.000000e-13
	return value, enabled
end

function material.thermal_tau_h(state)
	-- Hole energy relaxation time, units: s
	local enabled = true
	local value = 5.000000e-13
	return value, enabled
end

5.10 print() — a self-test helper

Finally, print() is a convenience function rather than a material property. It builds a representative state table (300 K, at the origin, zero photon density) and prints every quantity the script computes. This is useful for checking a script from the command line or in a standalone Lua interpreter before using it in a simulation — a quick way to confirm the numbers come out with the expected magnitudes.

6. Using a script in a simulation

Everything above concerns the Materials database. The screenshots in this section, however, come from the Electrical parameter editor inside a simulation — not the Materials database. This is where you choose whether a given layer uses simple fixed values or a Lua script.

By default, the electrical parameter editor holds a single fixed value for each property, stored in the simulation's sim.json file. In this mode the newer electrical script is not used at all: every box holds an ordinary number that you can type into directly ??. If you instead depress the Script button, the model switches to using the Lua script, and the parameter boxes now read “See script for value” ??. A box only shows “See script for value” if the corresponding function is enabled in the script; disabled functions continue to use the fixed value.

With the script active, a Code editor tab appears next to the Parameters tab, letting you view and edit the Lua source directly from within the simulation ??.

The small menu to the right of the Script button offers two actions ??:

• Import material script — opens the Materials database browser so you can select the Lua file you want the simulation to use.
• Populate json from script — runs the selected script and writes its computed values into the fixed values stored in sim.json, overwriting whatever was there.

The second action is what gives you the best of both worlds. You can either rely on the script fully — keeping the parameter boxes reading “See script for value” ?? — or you can use the script once to populate the simple edit boxes with its numbers ??, after which those values live independently in sim.json and can be tweaked by hand. There is no need to rely on the script permanently once you have pulled its values in. When you populate the JSON, OghmaNano shows a confirmation dialog listing exactly which values would be overwritten and what they would change to, so nothing is altered without your approval ??.

7. Where the script's values are used

Values produced by a material script are used throughout the model, not just in the electrical solver. Quantities such as Eg and Xi set the band-edge positions everywhere they are needed. For example, the energy levels drawn in the band-diagram / generation-rate viewer are taken from these values ??. Editing the band gap or affinity in a script therefore propagates automatically to the band diagram, the optical model, and the electrical solver.