Introduction: Recent investigations into the biomechanics of the brain have unveiled alteration in tissue stiffness triggered by external stimuli. For instance, visual stimulation effects can be measured in elasticity images of the cortex generated by functional magnetic resonance elastography (MRE). Such a mechanical characterization method combined with non-invasive brain stimulation (NIBS), a technique that seeks to selectively modulate particular parts of the brain using weak electrical currents, has the potential to influence research on various neurological disorders. In this in silico study, we aimed to elucidate individual and interdependent aspects related to a synchronized biomechanical imaging and non-invasive brain stimulation methodology. Magnetic resonance electrical impedance tomography (MREIT) was incorporated to the pipeline, providing a promising way of evaluating NIBS-induced electrical current patterns in the brain while leveraging MRE and transcranial alternating current stimulation (tACS) experimental settings.Methods: A mouse head model was assembled using open-access atlases to include five anatomical structures: skin/subcutaneous tissue, skull, cerebrospinal fluid (CSF), brain white and grey matters. MRE, tACS, and MREIT experiments were simulated using Comsol Multiphysics with Matlab Livelink. Synthetic MRE and MREIT data were processed using the subzone non-linear inversion and harmonic Bz algorithm, respectively, to reconstruct images of the distributed complex shear modulus and electrical conductivity.Results and Discussion: Lorentz body forces arising from simultaneous MRE and tACS elicited elastic waves of negligible amplitude compared with the extrinsic actuation levels reported in the literature, which allowed accurate reconstructions of the complex shear modulus. Qualitative electrical conductivity maps retrieved by MREIT accurately delineated anatomical regions of the brain model and could be used to recover reasonably accurate distributions of tACS-induced currents. This multi-physics approach has potential for translation to human brain imaging, and may provide more possibilities for the characterization of brain function together than in isolation.