dielectric_region = SphereRegion(
    value=4.0,
    radius=2.0*Angstrom,
    center=[1.0, 2.0 ,3.0]*Angstrom
)