dielectric_region = BoxRegion(
    4.0,
    xmin=0*Angstrom, xmax=41.439*Angstrom,
    ymin=0.5*Angstrom, ymax=3.57914*Angstrom,
    zmin=5*Angstrom, zmax=50.4135*Angstrom
    )