md_trajectory = nlread('md_trajectory.hdf5', MDTrajectory)[0]
thermal_conductivity=ThermalConductivity(
    md_trajectory,
    start_time=10*picosecond,
    end_time=md_trajectory.lastTime(),
    time_resolution=md_trajectory.timeInterval(),
    heat_source_tag_name='heat_source',
    heat_sink_tag_name='heat_sink',
)