This work has been supported by Deutsches Elektronen-Synchtrotron DESY, a member of the Helmholtz Association (HGF). Molecular dynamics calculations were enabled through the Maxwell computational resources operated at Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany. YuVV is grateful to Deutsche Forschungsgemeinschaft (grant VI 713/1-3) for support. A part of computations has been done at HPC facilities of the Universität zu Köln.