dielectric_region = TubeRegion(
    value=4.0,
    start_point=[1.0, 2.0, 3.0]*Angstrom, end_point=[4.0, 5.0, 6.0]*Angstrom,
    inner_radius=3.5*Angstrom,
    thickness=1.2*Angstrom
)