W pracy zaproponowano nową metodę numerycznego otrzymywania strukturporowatych metodą dynamiki molekularnej. Metoda polega na rozpoczęciu symulacji z ładunkami ekranowanymi i stopniowym zwiększaniu ich do pełnych ładunków jonowych. Metodę zastosowano do układu PbSiO3.