MRI image data acquisition and preprocessing
memory foam
All participants were screened for MRI safety prior to any scanning procedures. MRI scans were acquired, and clinical assessments were performed on the same day. Each participant underwent identical baseline and post-treatment MRI scans consisting of structural and resting-state functional MRI acquisitions. All MRI scans were acquired using a 3TGE Discovery MR750 scanner with a 32-channel head-neck imaging coil at the Center for Cognitive and Neurobiological Imaging at Stanford.High-resolution structural images using GE’s BRAVO sequence (three-dimensional, T1-weighted) were acquired for the whole brain (FOV = 256 × 256mm; matrix=256 × 256 voxel; slice thickness = 0.9 mm; TR = 2530 ms, TE = 2.98 ms, flip angle = 7°). During the 8-min resting state scan, participants were instructed to keep their eyes open, and their attention focused on a central fixation point, which consisted of a black screen with a white fixation cross.Participants were also instructed to let their minds wander freely and avoid repetitive thoughts. Whole-brain resting-state scans were collected with a 3X simultaneous multi-slice (i.e., multiband) acquisition echo planar imaging (EPI) sequence: TR = 2000ms, TE = 30 ms, flip angle = 77°, slice acceleration factor = 3, FOV=230 × 230mm, matrix = 128 × 128 voxel, 1.8 × 1.8 mm2 in-plane resolution, 87 contiguous axial slices parallel to the anterior commissure–posterior commissure line, yielding >1.4 M voxels every 2 s. The head motion of participants was effectively minimized using memory foam and inflatable padding. Participant alertness during the resting state task was monitored using in-scanner video cameras.As noted above, 5 subjects were excluded (see CONSORT diagram—Fig. MRI data were preprocessed using FMRIPREP version 20.2.025 [RRID:SCR 016216]. Each T1-weighted (T1w) volume was corrected for intensity non-uniformity using N4BiasFieldCorrection v2.1.0 [Resting state data were motion-corrected using mcflirt (FSL v5.0.9 [