A system for treating neurological conditions, the system including a processor, a transcranial magnetic stimulation coil and an application that is configured to direct the processor to stimulate a brain stimulation target using the transcranial magnetic stimulation coil in accordance with an aTBS protocol, wherein the application is configured to direct the processor to: obtain brain imaging data from a magnetic resonance imaging machine capable of obtaining functional magnetic resonance imaging, fMRI, image data of a patient's brain, where the brain imaging data describes neuronal activations within the patient's brain; map at least one region of interest to the patient's brain; locate functional subregions within the at least one region of interest based on the fMRI image data; determine functional relationships between at least two functional subregions; generate subregion parameters for each functional subregion; generate a target quality score for each functional subregion based on the subregion parameters by determining voxel size-weighted correlation coefficients for a set of subregions, the voxel size-weighted correlation coefficients for a given subregion being the sum of a two-dimensional matrix of Spearman's or Pearson's correlation coefficients derived from a hierarchical clustering algorithm describing the correlation coefficients between a surface region of interest subregion and all deep region of interest subregions; select as an aiTBS target a brain stimulation target based on its target quality score and the patient's neurological condition; direct the transcranial magnetic stimulation coil to apply magnetic stimulation according to an accelerated intermittent Theta-Burst Stimulation, aiTBS, protocol to the aiTBS target in the patient's brain during a plurality of sessions; separate each session by an intersession interval during which no magnetic stimulation is applied, wherein the intersession interval is between 25 and 120 minutes; and cause provision of, during each session in the plurality of sessions, the magnetic stimulation by generating a plurality of pulse trains that each comprise a set of pulses, wherein the frequency of pulses in each set of pulses is between 20 Hz and 70 Hz and the frequency of sets of pulses in each pulse train is between 3 Hz and 7 Hz, with pulse trains every 4 seconds to 10 seconds, wherein each pulse train is separated by an intertrain interval. obtain brain imaging data from a magnetic resonance imaging machine capable of obtaining functional magnetic resonance imaging, fMRI, image data of a patient's brain, where the brain imaging data describes neuronal activations within the patient's brain; map at least one region of interest to the patient's brain; locate functional subregions within the at least one region of interest based on the fMRI image data; determine functional relationships between at least two functional subregions; generate subregion parameters for each functional subregion; generate a target quality score for each functional subregion based on the subregion parameters by determining voxel size-weighted correlation coefficients for a set of subregions, the voxel size-weighted correlation coefficients for a given subregion being the sum of a two-dimensional matrix of Spearman's or Pearson's correlation coefficients derived from a hierarchical clustering algorithm describing the correlation coefficients between a surface region of interest subregion and all deep region of interest subregions; select as an aiTBS target a brain stimulation target based on its target quality score and the patient's neurological condition; direct the transcranial magnetic stimulation coil to apply magnetic stimulation according to an accelerated intermittent Theta-Burst Stimulation, aiTBS, protocol to the aiTBS target in the patient's brain during a plurality of sessions; separate each session by an intersession interval during which no magnetic stimulation is applied, wherein the intersession interval is between 25 and 120 minutes; and cause provision of, during each session in the plurality of sessions, the magnetic stimulation by generating a plurality of pulse trains that each comprise a set of pulses, wherein the frequency of pulses in each set of pulses is between 20 Hz and 70 Hz and the frequency of sets of pulses in each pulse train is between 3 Hz and 7 Hz, with pulse trains every 4 seconds to 10 seconds, wherein each pulse train is separated by an intertrain interval.