Skip to content

dft.density

Orbital normalization and density construction.

import mlx_atomistic.dft.density

def density_from_orbitals(orbitals: mx.array, grid: RealSpaceGrid, *, occupations: Sequence[float] | None = None) -> mx.array

Build spin-unpolarized electron density from occupied orbitals.

With no explicit occupations, each orbital is treated as doubly occupied: ρ(r) = 2Σᵢ |ψᵢ(r)|².

Parameters

NameTypeDefaultDescription
orbitalsmx.array
gridRealSpaceGrid
occupationsSequence[float] | NoneNone

Returns

  • mx.array
def normalize_orbitals(orbitals: mx.array, grid: RealSpaceGrid) -> mx.array

Normalize each orbital so that ∫ |ψᵢ(r)|² dr = 1.

Parameters

NameTypeDefaultDescription
orbitalsmx.array
gridRealSpaceGrid

Returns

  • mx.array