Erratum: Azimuthal Anisotropy at the Relativistic Heavy Ion Collider: The First and Fourth Harmonics [Phys. Rev. Lett. 92, 062301 (2004)]

This corrects the article DOI: 10.1103/PhysRevLett.92.062301.

Measurement of e^{+}e^{-} Momentum and Angular Distributions from Linearly Polarized Photon Collisions

The Breit-Wheeler process which produces matter and antimatter from photon collisions is experimentally investigated through the observation of 6085 exclusive electron-positron pairs in ultraperipheral Au+Au collisions at sqrt[s_{NN}]=200  GeV. The measurements reveal a large fourth-order angular modulation of cos4Δϕ=(16.8±2.5)% and smooth invariant mass distribution absent of vector mesons (ϕ, ω, and ρ) at the experimental limit of ≤0.2% of the observed yields. The differential cross section as a function of e^{+}e^{-} pair transverse momentum P_{⊥} peaks at low value with sqrt[⟨P_{⊥}^{2}⟩]=38.1±0.9  MeV and displays a significant centrality dependence. These features are consistent with QED calculations for the collision of linearly polarized photons quantized from the extremely strong electromagnetic fields generated by the highly charged Au nuclei at ultrarelativistic speed. The experimental results have implications for vacuum birefringence and for mapping the magnetic field which is important for emergent QCD phenomena.

Global Polarization of Ξ and Ω Hyperons in Au+Au Collisions at sqrt[s_{NN}]=200 GeV

Global polarization of Ξ and Ω hyperons has been measured for the first time in Au+Au collisions at sqrt[s_{NN}]=200  GeV. The measurements of the Ξ^{-} and Ξ[over ¯]^{+} hyperon polarization have been performed by two independent methods, via analysis of the angular distribution of the daughter particles in the parity violating weak decay Ξ→Λ+π, as well as by measuring the polarization of the daughter Λ hyperon, polarized via polarization transfer from its parent. The polarization, obtained by combining the results from the two methods and averaged over Ξ^{-} and Ξ[over ¯]^{+}, is measured to be ⟨P_{Ξ}⟩=0.47±0.10(stat)±0.23(syst)% for the collision centrality 20%-80%. The ⟨P_{Ξ}⟩ is found to be slightly larger than the inclusive Λ polarization and in reasonable agreement with a multiphase transport model. The ⟨P_{Ξ}⟩ is found to follow the centrality dependence of the vorticity predicted in the model, increasing toward more peripheral collisions. The global polarization of Ω, ⟨P_{Ω}⟩=1.11±0.87(stat)±1.97(syst)% was obtained by measuring the polarization of daughter Λ in the decay Ω→Λ+K, assuming the polarization transfer factor C_{ΩΛ}=1.

Brain fMRI during orientation selective epidural spinal cord stimulation

Epidural spinal cord stimulation (ESCS) is widely used for chronic pain treatment, and is also a promising tool for restoring motor function after spinal cord injury. Despite significant positive impact of ESCS, currently available protocols provide limited specificity and efficiency partially due to the limited number of contacts of the leads and to the limited flexibility to vary the spatial distribution of the stimulation field in respect to the spinal cord. Recently, we introduced Orientation Selective (OS) stimulation strategies for deep brain stimulation, and demonstrated their selectivity in rats using functional MRI (fMRI). The method achieves orientation selectivity by controlling the main direction of the electric field gradients using individually driven channels. Here, we introduced a similar OS approach for ESCS, and demonstrated orientation dependent brain activations as detected by brain fMRI. The fMRI activation patterns during spinal cord stimulation demonstrated the complexity of brain networks stimulated by OS-ESCS paradigms, involving brain areas responsible for the transmission of the motor and sensory information. The OS approach may allow targeting ESCS to spinal fibers of different orientations, ultimately making stimulation less dependent on the precision of the electrode implantation.

Nonmonotonic Energy Dependence of Net-Proton Number Fluctuations

Nonmonotonic variation with collision energy (sqrt[s_{NN}]) of the moments of the net-baryon number distribution in heavy-ion collisions, related to the correlation length and the susceptibilities of the system, is suggested as a signature for the quantum chromodynamics critical point. We report the first evidence of a nonmonotonic variation in the kurtosis times variance of the net-proton number (proxy for net-baryon number) distribution as a function of sqrt[s_{NN}] with 3.1  σ significance for head-on (central) gold-on-gold (Au+Au) collisions measured solenoidal tracker at Relativistic Heavy Ion Collider. Data in noncentral Au+Au collisions and models of heavy-ion collisions without a critical point show a monotonic variation as a function of sqrt[s_{NN}].

Directed, Elliptic, and Higher Order Flow Harmonics of Protons, Deuterons, and Tritons in Au+Au Collisions at sqrt[s_{NN}]=2.4 GeV

Flow coefficients v_{n} of the orders n=1-6 are measured with the High-Acceptance DiElectron Spectrometer (HADES) at GSI for protons, deuterons, and tritons as a function of centrality, transverse momentum, and rapidity in Au+Au collisions at sqrt[s_{NN}]=2.4  GeV. Combining the information from the flow coefficients of all orders allows us to construct for the first time, at collision energies of a few GeV, a multidifferential picture of the angular emission pattern of these particles. It reflects the complicated interplay between the effect of the central fireball pressure on the emission of particles and their subsequent interaction with spectator matter. The high precision information on higher order flow coefficients is a major step forward in constraining the equation of state of dense baryonic matter.

Quantitative but Not Qualitative Performance Changes in Predictive Motor Timing as a Result of Overtraining

The possibilities of substantial long-term improvement of predictive timing might be sometimes seen as limited, with scanty information of neural substrates underlying the potential learning process. To address this issue, we have investigated the performance of 21 baseball professionals and 21 matched controls in a predictive motor timing task previously shown to engage the cerebellum. Baseball players, hypothesized as a model of overtraining of the prediction of future state of the surroundings, showed significantly higher quantitative performance than nonathletic controls, with a substantial part of the baseball players reaching levels far beyond the range observed in common population. Furthermore, the qualitative performance profile of baseball players under various conditions as target speed and acceleration modes did not differ from the profile of healthy controls. Our results suggest that regular exigent training has the potential to vastly improve predictive motor timing. Moreover, the quantitative but not qualitative difference in the performance profile allows us to hypothesize that the selective honing of the same cerebellar processes and networks as in non-trained individuals is the substrate for the quantitative performance improvement, without substantial engagement of further neural nodes.

Differential diagnosis of tremor syndromes using MRI relaxometry

Differential diagnosis of the most common tremor syndromes - essential tremor (ET) and Parkinson's disease (PD) is burdened with high error rate. However, diagnostic MRI biomarkers applicable in this clinically highly relevant scenario remain an unfulfilled objective. The presented study was designed in search for possible candidate MRI protocols relevant for differential diagnostic process in tremor syndromes.10 non-advanced tremor-dominant PD patients meeting diagnostic criteria for clinically established PD, 12 isolated ET patients and 16 healthy controls were enrolled into this study. The study focused on relaxation MRI protocols - T1, T2, adiabatic T1ρ and adiabatic T2ρ due to their relatively low post-processing requirements enabling implementation into routine clinical practice. Compared to ET, PD patients had significantly longer T2 relaxation times in striata with dominant findings in the putamen contralateral to the clinically more affected body side. This difference was driven by alterations in the PD group as confirmed in the complementary comparison with healthy controls. According to the receiver operating characteristic analysis, this region provided a reasonable sensitivity of 0.91 and specificity of 0.89 in the differential diagnosis of PD and ET. In PD patients, we further found prolonged T1ρ in the substantia nigra compared to ET and healthy controls, and shorter T2 and T2ρ in the cerebellum compared to healthy controls. T2 relaxation time in the putamen contralateral to the clinically more affected body side is a plausible candidate diagnostic marker for the differentiation of PD and ET.

Structural Alterations in Deep Brain Structures in Type 1 Diabetes

Even though well known in type 2 diabetes, the existence of brain changes in type 1 diabetes (T1D) and both their neuroanatomical and clinical features are less well characterized. To fill the void in the current understanding of this disease, we sought to determine the possible neural correlate in long-duration T1D at several levels, including macrostructural, microstructural cerebral damage, and blood flow alterations. In this cross-sectional study, we compared a cohort of 61 patients with T1D with an average disease duration of 21 years with 54 well-matched control subjects without diabetes in a multimodal MRI protocol providing macrostructural metrics (cortical thickness and structural volumes), microstructural measures (T1-weighted/T2-weighted [T1w/T2w] ratio as a marker of myelin content, inflammation, and edema), and cerebral blood flow. Patients with T1D had higher T1w/T2w ratios in the right parahippocampal gyrus, the executive part of both putamina, both thalami, and the cerebellum. These alterations were reflected in lower putaminal and thalamic volume bilaterally. No cerebral blood flow differences between groups were found in any of these structures, suggesting nonvascular etiologies of these changes. Our findings implicate a marked nonvascular disruption in T1D of several essential neural nodes engaged in both cognitive and motor processing.

Orientation-selective and directional deep brain stimulation in swine assessed by functional MRI at 3T

Functional MRI (fMRI) has become an important tool for probing network-level effects of deep brain stimulation (DBS). Previous DBS-fMRI studies have shown that electrical stimulation of the ventrolateral (VL) thalamus can modulate sensorimotor cortices in a frequency and amplitude dependent manner. Here, we investigated, using a swine animal model, how the direction and orientation of the electric field, induced by VL-thalamus DBS, affects activity in the sensorimotor cortex. Adult swine underwent implantation of a novel 16-electrode (4 rows × 4 columns) directional DBS lead in the VL thalamus. A within-subject design was used to compare fMRI responses for (1) directional stimulation consisting of monopolar stimulation in four radial directions around the DBS lead, and (2) orientation-selective stimulation where an electric field dipole was rotated 0°−360° around a quadrangle of electrodes. Functional responses were quantified in the premotor, primary motor, and somatosensory cortices. High frequency electrical stimulation through leads implanted in the VL thalamus induced directional tuning in cortical response patterns to varying degrees depending on DBS lead position. Orientation-selective stimulation showed maximal functional response when the electric field was oriented approximately parallel to the DBS lead, which is consistent with known axonal orientations of the cortico-thalamocortical pathway. These results demonstrate that directional and orientation-selective stimulation paradigms in the VL thalamus can tune network-level modulation patterns in the sensorimotor cortex, which may have translational utility in improving functional outcomes of DBS therapy.

First Measurement of Λ_{c} Baryon Production in Au+Au Collisions at sqrt[s_{NN}]=200 GeV

We report on the first measurement of the charmed baryon Λ_{c}^{±} production at midrapidity (|y|<1) in Au+Au collisions at sqrt[s_{NN}]=200  GeV collected by the STAR experiment at the Relativistic Heavy Ion Collider. The Λ_{c}/D^{0} [denoting (Λ_{c}^{+}+Λ_{c}^{-})/(D^{0}+D[over ¯]^{0})] yield ratio is measured to be 1.08±0.16  (stat)±0.26  (sys) in the 0%-20% most central Au+Au collisions for the transverse momentum (p_{T}) range 3<p_{T}<6  GeV/c. This is significantly larger than the pythia model calculations for p+p collisions. The measured Λ_{c}/D^{0} ratio, as a function of p_{T} and collision centrality, is comparable to the baryon-to-meson ratios for light and strange hadrons in Au+Au collisions. Model calculations including coalescence hadronization for charmed baryon and meson formation reproduce the features of our measured Λ_{c}/D^{0} ratio.

Orientation selective deep brain stimulation of the subthalamic nucleus in rats

Deep brain stimulation (DBS) has become an important tool in the management of a wide spectrum of diseases in neurology and psychiatry. Target selection is a vital aspect of DBS so that only the desired areas are stimulated. Segmented leads and current steering have been shown to be promising additions to DBS technology enabling better control of the stimulating electric field. Recently introduced orientation selective DBS (OS-DBS) is a related development permitting sensitization of the stimulus to axonal pathways with different orientations by freely controlling the primary direction of the electric field using multiple contacts. Here, we used OS-DBS to stimulate the subthalamic nucleus (STN) in healthy rats while simultaneously monitoring the induced brain activity with fMRI. Maximal activation of the sensorimotor and basal ganglia-thalamocortical networks was observed when the electric field was aligned mediolaterally in the STN pointing in the lateral direction, while no cortical activation was observed with the electric field pointing medially to the opposite direction. Such findings are consistent with mediolateral main direction of the STN fibers, as seen with high resolution diffusion imaging and histology. The asymmetry of the OS-DBS dipolar field distribution using three contacts along with the potential stimulation of the internal capsule, are also discussed. We conclude that OS-DBS offers an additional degree of flexibility for optimization of DBS of the STN which may enable a better treatment response.

Rotating frame MRI relaxations as markers of diffuse white matter abnormalities in multiple sclerosis

•

T1ρ and RAFF4 - MRI protocols sensitive to slow motional regimes – are able to detect changes in T2w-defined normally appearing white matter of patients with multiple sclerosis.

•

RAFF4, T1ρ and T2ρ showed differences in hippocampus in patients with multiple sclerosis despite the absence of alterations in resting-state functional MRI metrics.

•

Hence, RAFF4, T1ρ and T2ρ hold promise as potential non-invasive tools for monitoring MS activity and eventually for the evaluation of therapeutic effects.

Even though MRI visualization of white matter lesions is pivotal for the diagnosis and management of multiple sclerosis (MS), the issue of detecting diffuse brain tissue damage beyond the apparent T2-hyperintense lesions continues to spark considerable interest. Motivated by the notion that rotating frame MRI methods are sensitive to slow motional regimes critical for tissue characterization, here we utilized novel imaging protocols of rotating frame MRI on a clinical 3 Tesla platform, including adiabatic longitudinal, T1ρ, and transverse, T2ρ, relaxation methods, and Relaxation Along a Fictitious Field (RAFF) in the rotating frame of rank 4 (RAFF4), in 10 relapsing-remitting multiple sclerosis patients and 10 sex- and age-matched healthy controls. T1ρ, T2ρ and RAFF4 relaxograms extracted from the whole white matter exhibited a significant shift towards longer relaxation time constants in MS patients as compared to controls. T1ρ and RAFF4 detected alterations even when considering only regions of normally appearing white matter (NAWM), while other MRI metrics such as T1w/T2w ratio and diffusion tensor imaging measures failed to find group differences. In addition, RAFF4, T2ρ and, to a lesser extent, T1ρ showed differences in subcortical grey matter structures, mainly hippocampus, whereas no functional changes in this region were detected in resting-state functional MRI metrics. We conclude that rotating frame MRI techniques are exceptionally sensitive methods for the detection of subtle abnormalities not only in NAWM, but also in deep grey matter in MS, where they surpass even highly sensitive measures of functional changes, which are often suggested to precede detectable structural alterations. Such abnormalities are consistent with a wide spectrum of different, but interconnected pathological features of MS, including the loss of neuronal cells and their axons, decreased levels of myelin even in NAWM, and altered iron content.

Behavioral and Neuroanatomical Account of Impulsivity in Parkinson's Disease

Impulse control disorder (ICD) is a major non-motor complication of Parkinson's disease (PD) with often devastating consequences for patients' quality of life. In this study, we aimed to characterize the phenotype of impulsivity in PD and its neuroanatomical correlates.

Methods: Thirty-seven PD patients (15 patients with ICD, 22 patients without ICD) and 36 healthy controls underwent a neuropsychological battery. The test battery consisted of anxiety and depression scales, self-report measures of impulsivity (Barratt scale and UPPS-P), behavioral measures of impulsive action (Go/No-Go task, Stop signal task) and impulsive choice (Delay discounting, Iowa gambling task), and measures of cognitive abilities (working memory, attention, executive function). Patients and controls underwent structural MRI scanning.

Results: Patients with ICD had significantly higher levels of self-reported impulsivity (Barratt scale and Lack of perseverance from UPPS-P) in comparison with healthy controls and non-impulsive PD patients, but they performed similarly in behavioral tasks, except for the Iowa gambling task. In this task, patients with ICD made significantly less risky decisions than patients without ICD and healthy controls. Patients without ICD did not differ from healthy controls in self-reported impulsivity or behavioral measurements. Both patient groups were more anxious and depressive than healthy controls. MRI scanning revealed structural differences in cortical areas related to impulse control in both patient groups. Patients without ICD had lower volumes and cortical thickness of bilateral inferior frontal gyrus. Patients with ICD had higher volumes of right caudal anterior cingulate and rostral middle frontal cortex.

Conclusions: Despite the presence of ICD as confirmed by both clinical follow-up and self-reported impulsivity scales and supported by structural differences in various neural nodes related to inhibitory control and reward processing, patients with ICD performed no worse than healthy controls in various behavioral tasks previously hypothesized as robust impulsivity measures. These results call for caution against impetuous interpretation of behavioral tests, since various factors may and will influence the ultimate outcomes, be it the lack of sensitivity in specific, limited ICD subtypes, excessive caution of ICD patients during testing due to previous negative experience rendering simplistic tasks insufficient, or other, as of now unknown aspects, calling for further research.

First Observation of the Directed Flow of D^{0} and D^{0}[over ¯] in Au+Au Collisions at sqrt[s_{NN}]=200 GeV

We report the first measurement of rapidity-odd directed flow (v_{1}) for D^{0} and D^{0}[over ¯] mesons at midrapidity (|y|<0.8) in Au+Au collisions at sqrt[s_{NN}]=200  GeV using the STAR detector at the Relativistic Heavy Ion Collider. In 10-80% Au+Au collisions, the slope of the v_{1} rapidity dependence (dv_{1}/dy), averaged over D^{0} and D^{0}[over ¯] mesons, is -0.080±0.017(stat)±0.016(syst) for transverse momentum p_{T} above 1.5  GeV/c. The absolute value of D^{0} meson dv_{1}/dy is about 25 times larger than that for charged kaons, with 3.4σ significance. These data give a unique insight into the initial tilt of the produced matter, and offer constraints on the geometric and transport parameters of the hot QCD medium created in relativistic heavy-ion collisions.

First measurement of neutrino oscillation parameters using neutrinos and antineutrinos by NOvA

The NOvA experiment has seen a 4.4σ signal of ν[over ¯]_{e} appearance in a 2 GeV ν[over ¯]_{μ} beam at a distance of 810 km. Using 12.33×10^{20} protons on target delivered to the Fermilab NuMI neutrino beamline, the experiment recorded 27 ν[over ¯]_{μ}→ν[over ¯]_{e} candidates with a background of 10.3 and 102 ν[over ¯]_{μ}→ν[over ¯]_{μ} candidates. This new antineutrino data are combined with neutrino data to measure the parameters |Δm_{32}^{2}|=2.48_{-0.06}^{+0.11}×10^{-3}  eV^{2}/c^{4} and sin^{2}θ_{23} in the ranges from (0.53-0.60) and (0.45-0.48) in the normal neutrino mass hierarchy. The data exclude most values near δ_{CP}=π/2 for the inverted mass hierarchy by more than 3σ and favor the normal neutrino mass hierarchy by 1.9σ and θ_{23} values in the upper octant by 1.6σ.

Polarization of Λ (Λ[over ¯]) Hyperons along the Beam Direction in Au+Au Collisions at sqrt[s_{NN}]=200 GeV

The Λ (Λ[over ¯]) hyperon polarization along the beam direction has been measured in Au+Au collisions at sqrt[s_{NN}]=200  GeV, for the first time in heavy-ion collisions. The polarization dependence on the hyperons' emission angle relative to the elliptic flow plane exhibits a second harmonic sine modulation, indicating a quadrupole pattern of the vorticity component along the beam direction, expected due to elliptic flow. The polarization is found to increase in more peripheral collisions, and shows no strong transverse momentum (p_{T}) dependence at p_{T} greater than 1  GeV/c. The magnitude of the signal is about 5 times smaller than those predicted by hydrodynamic and multiphase transport models; the observed phase of the emission angle dependence is also opposite to these model predictions. In contrast, the kinematic vorticity calculations in the blast-wave model tuned to reproduce particle spectra, elliptic flow, and the azimuthal dependence of the Gaussian source radii measured with the Hanbury Brown-Twiss intensity interferometry technique reproduce well the modulation phase measured in the data and capture the centrality and transverse momentum dependence of the polarization signal.

Observation of Excess J/ψ Yield at Very Low Transverse Momenta in Au+Au Collisions at sqrt[s_{NN}]=200 GeV and U+U Collisions at sqrt[s_{NN}]=193 GeV

We report on the first measurements of J/ψ production at very low transverse momentum (p_{T}<0.2  GeV/c) in hadronic Au+Au collisions at sqrt[s_{NN}]=200  GeV and U+U collisions at sqrt[s_{NN}]=193  GeV. Remarkably, the inferred nuclear modification factor of J/ψ at midrapidity in Au+Au (U+U) collisions reaches about 24 (52) for p_{T}<0.05  GeV/c in the 60%-80% collision centrality class. This noteworthy enhancement cannot be explained by hadronic production accompanied by cold and hot medium effects. In addition, the dN/dt distribution of J/ψ for the very low p_{T} range is presented for the first time. The distribution is consistent with that expected from the Au nucleus and shows a hint of interference. Comparison of the measurements to theoretical calculations of coherent production shows that the excess yield can be described reasonably well and reveals a partial disruption of coherent production in semicentral collisions, perhaps due to the violent hadronic interactions. Incorporating theoretical calculations, the results strongly suggest that the dramatic enhancement of J/ψ yield observed at extremely low p_{T} originates from coherent photon-nucleus interactions. In particular, coherently produced J/ψ's in violent hadronic collisions may provide a novel probe of the quark-gluon plasma.

Azimuthal Harmonics in Small and Large Collision Systems at RHIC Top Energies

The first (v_{1}^{fluc}), second (v_{2}), and third (v_{3}) harmonic coefficients of the azimuthal particle distribution at midrapidity are extracted for charged hadrons and studied as a function of transverse momentum (p_{T}) and mean charged particle multiplicity density ⟨N_{ch}⟩ in U+U (sqrt[s_{NN}]=193  GeV), Au+Au, Cu+Au, Cu+Cu, d+Au, and p+Au collisions at sqrt[s_{NN}]=200  GeV with the STAR detector. For the same ⟨N_{ch}⟩, the v_{1}^{fluc} and v_{3} coefficients are observed to be independent of the collision system, while v_{2} exhibits such a scaling only when normalized by the initial-state eccentricity (ϵ_{2}). The data also show that ln(v_{2}/ϵ_{2}) scales linearly with ⟨N_{ch}⟩^{-1/3}. These measurements provide insight into initial-geometry fluctuations and the role of viscous hydrodynamic attenuation on v_{n} from small to large collision systems.

Consensus paper: Decoding the Contributions of the Cerebellum as a Time Machine. From Neurons to Clinical Applications

Time perception is an essential element of conscious and subconscious experience, coordinating our perception and interaction with the surrounding environment. In recent years, major technological advances in the field of neuroscience have helped foster new insights into the processing of temporal information, including extending our knowledge of the role of the cerebellum as one of the key nodes in the brain for this function. This consensus paper provides a state-of-the-art picture from the experts in the field of the cerebellar research on a variety of crucial issues related to temporal processing, drawing on recent anatomical, neurophysiological, behavioral, and clinical research.The cerebellar granular layer appears especially well-suited for timing operations required to confer millisecond precision for cerebellar computations. This may be most evident in the manner the cerebellum controls the duration of the timing of agonist-antagonist EMG bursts associated with fast goal-directed voluntary movements. In concert with adaptive processes, interactions within the cerebellar cortex are sufficient to support sub-second timing. However, supra-second timing seems to require cortical and basal ganglia networks, perhaps operating in concert with cerebellum. Additionally, sensory information such as an unexpected stimulus can be forwarded to the cerebellum via the climbing fiber system, providing a temporally constrained mechanism to adjust ongoing behavior and modify future processing. Patients with cerebellar disorders exhibit impairments on a range of tasks that require precise timing, and recent evidence suggest that timing problems observed in other neurological conditions such as Parkinson's disease, essential tremor, and dystonia may reflect disrupted interactions between the basal ganglia and cerebellum.The complex concepts emerging from this consensus paper should provide a foundation for further discussion, helping identify basic research questions required to understand how the brain represents and utilizes time, as well as delineating ways in which this knowledge can help improve the lives of those with neurological conditions that disrupt this most elemental sense. The panel of experts agrees that timing control in the brain is a complex concept in whom cerebellar circuitry is deeply involved. The concept of a timing machine has now expanded to clinical disorders.

Tuning Neuromodulation Effects by Orientation Selective Deep Brain Stimulation in the Rat Medial Frontal Cortex

Previous studies that focused on treating major depressive disorder with conventional deep brain stimulation (DBS) paradigms produced inconsistent results. In this proof-of-concept preclinical study in rats (n = 8), we used novel paradigms of orientation selective DBS for stimulating the complex circuitry crossing the infralimbic cortex, an area considered analogous to human subgenual cingulate cortex. Using functional MRI at 9.4 T, we monitored whole brain responses to varying the electrical field orientation of DBS within the infralimbic cortex. Substantial alterations of functional MRI responses in the amygdala, a major node connected to the infralimbic cortex implicated in the pathophysiology of depression, were observed. As expected, the activation cluster near the electrode was insensitive to the changes of the stimulation orientation. Hence, our findings substantiate the ability of orientation selective stimulation (OSS) to recruit neuronal pathways of distinct orientations relative to the position of the electrode, even in complex circuits such as those involved in major depressive disorder. We conclude that OSS is a promising approach for stimulating brain areas that inherently require individualisation of the treatment approach.