Numerical prediction of frequency dependent 3D maps of mechanical index thresholds in ultrasonic brain therapy