I'm performing multiscale calculations with mechanical embedding in ORCA and I have a question about how atomic charges are handled during geometry optimization and reaction coordinate scans. Are the electron density and resulting atomic charges recalculated and updated throughout each iteration of the optimization/scan process?