Neurophotonics

Laser Interstitial Thermal Therapy (LITT) is a minimally invasive neurosurgical technique in which a stereotactically-implanted fiber delivers thermal energy to ablate intracranial lesions. Existing computer-assisted planning systems optimize trajectories against a one-dimensional line abstraction, then approximate the ablation zone as a fixed-radius cylinder post-hoc to estimate coverage. Trajectories selected as optimal under this model are not guaranteed to remain optimal once the cylindrical extent is applied, which introduces a mismatch between predicted and true ablation coverage. This may also underestimate spillover into surrounding healthy tissue. We present OptiLITT, a treatment planning system that represents the laser probe as a cylindrical ablation volume from the onset of optimization, jointly solving dual-fiber placement, lesion coverage, and healthy-tissue spillover as a single coupled problem. All planning parameters are exposed through a user-configurable graphical user interface supporting intraoperative refinement between planning stages.

Significance: Continuous, noninvasive, and depth-resolved monitoring of blood-flow-related tissue dynamics remains an important unmet need. Speckle-contrast optical spectroscopy (SCOS), including interferometric implementations such as iSCOS, provides a scalable optical route to blood-flow sensing, but conventional continuous-wave approaches lack intrinsic depth selectivity. Time-of-flight (TOF) gating offers a way to separate superficial and deeper dynamic contributions in layered tissues, such as skin-muscle or scalp-cortex, by resolving photon path lengths. Aim: We introduce a swept-source, single-channel implementation of interferometric speckle-contrast optical spectroscopy (iSCOS) to obtain TOF-resolved temporal speckle contrast, {kappa}^2, from the measured field autocorrelation g_1, and evaluate its feasibility for depth-resolved blood-flow sensing. Approach: A swept-source iNIRS system operating at 780 nm acquired interferometric signals, which were Fourier-transformed along the optical-frequency axis to recover complex TOF-resolved speckle fields. Temporal speckle contrast was then estimated at each TOF gate indirectly from g_1 using the speckle-visibility relation. Diffusion-based numerical simulations were first used to compare the direct variance-based estimator and the indirect g_1-based estimator under varying reduced scattering coefficient, diffusion coefficient, additive noise level, and bi-layer geometry. Because the simulations showed that the g_1-derived {kappa}^2 estimator was substantially less sensitive to additive noise than the direct estimator, this estimator was used for the main phantom and in vivo analyses, while the direct estimator served as a simulation comparator. The g_1-derived estimator was then applied to liquid and bi-layer phantoms, followed by proof-of-concept in vivo measurements on the human forearm during cuff occlusion and on the forehead during a Sudoku task. Results: TOF-resolved kappa2 curves recovered with the g_1-derived estimator matched DWS theory across scattering coefficients, photon path lengths, and exposure times. The estimator preserved theoretical accuracy for additive noise amplitudes up to 50% of the field amplitude, whereas the direct variance estimator showed substantial noise-induced bias and required correction. Bi-layer simulations and phantom experiments reproduced the predicted direction and onset of TOF-dependent decorrelation-rate trends in layered media. In vivo, the recovered blood-flow index tracked the expected TOF-dependent cuff-occlusion and reactive-hyperemia response in the forearm. During the single-subject Sudoku task, the left-forehead recording showed a TOF-dependent relative blood-flow-index increase of +0.8 {+/-} 1.9% at TOF = 400 ps, +9.8 {+/-} 2.2% at TOF = 600 ps, and +15.2 {+/-} 5.6% at TOF = 800 ps. This pattern is consistent with increased sensitivity to deeper tissue at longer photon path lengths, but requires cohort-level validation before quantitative interpretation as cognitive activation. Conclusions: Coupling temporal speckle-contrast analysis with swept-source iNIRS yields a proof-of-concept, depth-resolved platform for blood-flow sensing. By estimating TOF-resolved speckle contrast through the g_1-derived {kappa}^2 route, TOF-iSCOS suppresses additive-noise bias while preserving sensitivity to deeper dynamic tissue layers. The present single-channel results bridge continuous-wave iSCOS, interferometric NIRS and time-domain diffuse correlation spectroscopy (TD-DCS), and motivate future multi-channel and cohort studies for scalable cortical hemodynamic monitoring.

Quantifying blood-brain barrier (BBB) integrity from fluorescence microscopy remains limited by subjective scoring and categorical classification methods that lack reproducibility. For objective and consistent BBB phenotyping, we present two semi-automated image-analysis pipelines that replace manual scoring with quantitative, continuous-variable measurements. Our in vitro pipeline, implemented in Python, quantifies the connectivity of tight junction structures by measuring discrete ZO-1 fragment objects within manually traced junction regions. It outputs continuous metrics including average fragment area, total junctional area, and a junctional fragmentation ratio that captures degree of ZO-1 continuity versus discontinuity. In human brain microvascular endothelial cells subjected to glycocalyx component knockdown, the pipeline detected significantly reduced fragment area (37% decrease for both CD44 and syndecan-1 (SDC1) knockdown, p = 0.0148 and 0.0084) and junctional fragmentation ratio (p = 0.0061 and 0.0137). Our in vivo pipeline integrates ilastik-based pixel classification with FIJI macro automation to quantify vascular marker colocalization and to separate vessel signal from microglial contamination within a single fluorescence channel, eliminating the need for dedicated counterstains. Applied across four mouse cohorts [young, aged, Alzheimer's, traumatic brain injury (TBI)] and three brain regions [prefrontal cortex (PFC), hippocampus, midbrain], the pipeline revealed concurrent ZO-1 loss and ICAM-1 elevation in the PFC and hippocampus of aged and Alzheimer's mice, with Alzheimer's-specific doubling of eNOS occurring in the PFC (p = 0.0013). TBI mice showed persistent ZO-1 loss with transient ICAM-1 and eNOS changes. Both deterministic pipelines are available on GitHub and designed for adoption beyond the specific markers and systems analyzed here.

Accurate measurement of cellular oxygen levels is essential for understanding the balance between oxygen demand and supply in tissues. However, conventional methods only yield compromised results. We harnessed the oxygen dependence of bioluminescence to develop OxyBLI--a noninvasive optical method that directly monitors oxygen levels in specific cell populations of intact experimental animals. We characterized OxyBLI signals across various critical situations associated with common interventions. Hypoxic breathing and subsequent systemic tissue hypoxia caused blood to be redistributed in a way that prioritized brain oxygenation. In contrast, hyperoxic breathing sharply increased tissue oxygenation, which promptly returned to the target level owing to a vasoconstrictor response. These findings are expected to help resolve the long-standing clinical issue regarding the risks and benefits of administering supplemental oxygen to acutely ill patients. Our multifaceted approach, which presents multiple challenges to individual animals over time, will advance our understanding of the delicate interaction between hypoxia and hyperoxia.

Designing transcutaneous skeletal muscle oxygenation (SmO2) sensors requires jointly optimizing source--detector geometry and wavelength selection while guaranteeing performance across populations that vary in subcutaneous fat thickness and skin pigmentation. We present a multi-fidelity Bayesian optimization (MFBO) framework that couples Monte Carlo light-transport simulations at two photon-count fidelities to a distributionally robust design objective. An autoregressive Gaussian-process surrogate learns the correlation between inexpensive low-photon-count and accurate high-photon-count simulations, and a cost-aware acquisition function decides both where and at what fidelity to sample. Robustness across the population is enforced with Conditional Value-at-Risk (CVaR) and entropic-risk (ERM) objectives that target worst-case subjects rather than the population average. On a five-layer forearm tissue model with anthropometric variability we find (i) a fidelity regime that is favorable for MFBO where the low-fidelity surrogate is rank-informative (Spearman {rho} = 0.84) but biased, at 100x lower cost; (ii) MFBO attains 23% higher robust sensitivity than a strong high-fidelity single-fidelity baseline at equal budget (p = 0.035), and avoids the optimistic bias that causes low-fidelity-only optimization to collapse when its designs are validated at high fidelity; (iii) CVaR/ERM objectives improve worst-case tail performance by {approx}23% relative to a mean objective without sacrificing average sensitivity; and (iv) discovered designs improve robust tail sensitivity by roughly 3--6x over commercial and heuristic optode layouts, with the largest gains in the high-fat and high-melanin subpopulations. The methodology bridges stochastic light-transport physics with sample-efficient machine-learning optimization and generalizes to cerebral oximetry, photodynamic therapy planning, and wearable physiological monitors.

BackgroundConventional rodent models for the study of corneal pain commonly evoke eye blink reflex using methods that indiscriminately activate polymodal nociceptors, mechanoreceptors, and thermoreceptors at temporal resolutions that dont closely match the sub-second timescale of underlying neural dynamics. New methodWe introduce a novel automated behavioral paradigm for detecting blink reflexes in transgenic TRPV1-ChR2-EYFP mice, enabled by cell-type-specific, millisecond-precision optogenetic stimulation of corneal nociceptors (490 nm light). Using multi-feature quantification, we achieve robust automated detection using univariate and multivariate classifiers. ResultsTRPV1-ChR2-EYFP mice exhibited blink reflexes to high-intensity blue light (490 nm, 10 ms pulses) in a threshold-dependent manner (N=3). Blink probability was 77.1 {+/-} 17.1% at high intensity (2.77 mW/mm2) versus 4.2 {+/-} 4.2% at low intensity (0.46 mW/mm2). Red light (638 nm) produced no intensity-dependent change. Noxious air puff evoked blinks in >95% of trials under all conditions. DeepLabCut-based pose estimation extracted six features quantifying the blink reflex, enabling automated detection with [≥]98% accuracy using univariate and multivariate classifiers. Comparison with existing methodsUnlike conventional air puff paradigms, this optogenetic approach enables precise, cell-type-specific stimulation of corneal nociceptors, supporting automated analysis of blink responses at sub-second resolution. ConclusionsThis video tracking behavioral method using machine learning algorithms that accurately classify blink versus no-blink enables high-throughput and observer-independent empirical assessment of blink reflex, suggestive of corneal pain. Moreover, inducing blink reflex in TRPV1-ChR2 mice using high-intensity blue light also demonstrates nociceptive-specific behavioral responses analogous to somatosensory optogenetically-evoked hindpaw pain in the same animal genotype.

A test technique for measuring high-frequency transient current components in deep brain tissue is presented. The technique applies a voltage pulse with a high value in dV/dt, generating a corresponding current pulse with high dI/dt that can elicit measurable transient current responses from the electrode/tissue interface and adjacent brain tissue; responses are analyzed in the frequency domain by Fast Fourier Transform at a 200 kHz sampling frequency. The method was motivated by prior evidence that ferritin and neuromelanin in catecholaminergic tissue may support high-frequency conduction properties that have not previously been characterized in vivo. The protocol was applied in 277 measurements across five Sprague Dawley rats at cortical and basal ganglia locations in different locations in the brain. Preliminary spectral results show differences between catecholaminergic regions and cortical tissue that support further development and validation of the method.

Optical recording from large numbers of neurons is an indispensable technique for studying neuronal ensembles. We use optical sectioning through pinhole illumination to reduce the background fluorescence (F0) and increase the optical signal ({Delta}F/F0) in ex vivo brain slices densely labeled with GCaMP6f, allowing an ordinary fluorescence microscope to capture calcium transients from over 300 individual CA1 neurons - a marked increase compared to ordinary wide field fluorescence illumination. Multiple layers of overlapping neurons can be identified by their locations and the shape in space of their {Delta}F/F0 images. A single pinhole mask was placed at the field stop of a wide field illuminator, and the image of the pinhole was projected onto the tissue by a 20X NA 0.95 water immersion objective (Olympus). This created an illuminated disk with a diameter of [~]200 m and optical sections of hippocampal CA1 pyramidal layer tissue [~]100 m thick. This illumination blocked a large fraction of the F0, which in turn increased the {Delta}F/F0 5-10-fold compared to that of wide field illumination. When putative pyramidal neurons fire sparsely in the brain slice, up to 300 partially superimposed neurons can be identified by their shape and spatial location in the thick ([~]480 m) ex vivo slice in the CA1 area surrounding the pinhole image. The signal-to-noise ratio was adequate even at a low excitation light level of [~]20k photoelectrons per pixel well on the camera, allowing for 3,000 seconds of total recording time without significant bleaching. This pinhole "half confocal" method has created a useful way to sample calcium transient signals in thick tissue with a large population of neurons densely labeled with GCaMP-6f.

Fiber photometry (FP) has become a tool of choice for in vivo monitoring of genetically encoded biosensors. The ability to record and optogenetically manipulate circuits through the same fiber stub is powerful, but limited, as biosensors typically do not sample membrane voltage, leaving the experimenter blind to the direct effects of opsin photoactivation. Here we developed the Photometry Device with Optogenetics and Electrophysiology (P-DOpE) probe, fabricated via a new convergence taper-break (CTB) method that integrates industry-standard silica optical waveguides with low-impedance metal electrodes that can be arranged in experimenter-defined configurations. We demonstrate that chronically implanted P-DOpE probes provide months-long recordings of local field potential, single unit recording, and fiber photometry, with parallel optogenetic circuit perturbation. Conducting fiber photometry with same-site optogenetic stimulation, we identified a robust fluorescence signal that scaled with network activity and survived biosensor antagonism. As this confound could not be eliminated with standard isosbestic controls, we propose a simple correction strategy. As a first application, we used the probe to test a proposed mechanism for focal modulation of noradrenergic signaling in and by cortical circuits receiving afferents from the locus coeruleus. We found that increasing spiking activity in CA1 amplifies noradrenergic signaling evoked by contextual arousal by ~50%, but does not induce norepinephrine release in the absence of a phasic trigger - thus supporting the central prediction of the glutamate amplifies noradrenergic effects (GANE) hypothesis. The P-DOpE probe thus enables optogenetic manipulation and multimodal readout in a configurable low-cost, scalable, and robust format.

Transcranial photobiomodulation (PBM) is an emerging non-invasive brain stimulation method that is thought to increase neural metabolism by stimulating ATP production by the mitochondria. However, the mechanisms of action and the effects on the human brain are still unclear. In the present study, we investigated the potential of this method to enhance Blood Oxygen Level Dependent (BOLD) responses during the execution of a motor task in young and aged participants. Sixty young and aged participants were included in this single-blinded, sham-controlled, randomised, crossover study. They underwent an fMRI recording before and after 24-min stimulation with a 80-LEDs helmet emitting transcranially red and near infrared light. Post vs Pre BOLD signal was compared between PBM and SHAM, in each group. At baseline, aged participants showed reduced BOLD signal compared to young ones, in key regions of the sensorimotor processing, principally the left primary motor cortex and striatum. Transcranial PBM did not have a real impact in the young group. However in aged participants it increased BOLD signal in some regions that were underactivated compared to the young group at baseline. In particular, regional analysis showed increased BOLD response in the left primary motor cortex, and right dorsal and ventral premotor regions and striatum. These results suggest that transcranial PBM can increase fMRI BOLD activity in the task-related regions, particularly in aged subjects. Further research are needed to distinguish neural from vascular effects in transcranial PBM.

Functional ultrasound (fUS) maps cerebral blood volume (CBV) but lacks molecular and neuronal specificity. By simultaneously integrating fUS with optical imaging, we show that fUS-derived CBV correlates with both optically measured hemoglobin and neuronal calcium activity in awake mice. We further derive hemodynamic response functions linking calcium activity to CBV during spontaneous and sensory-evoked activity. Application to a mouse glioblastoma model demonstrates utility for studying neurovascular dysfunction in complex neuropathologies.

IntroductionChemogenetic tools such as Designer Receptors Exclusively Activated by Designer Drugs (DREADDs) provide a powerful means to causally manipulate defined neuronal populations in vivo. While chemogenetic fMRI studies have consistently demonstrated robust hemodynamic responses following circuit perturbation, considerably less is known about the accompanying metabolic consequences. Functional magnetic resonance spectroscopy (fMRS) offers the potential to probe these neurochemical processes, yet the relationship between hemodynamic and metabolic responses remains poorly understood. Here, we combined chemogenetics, pharmacological fMRI (ph-fMRI), and proton magnetic resonance spectroscopy (1H-MRS/fMRS) at 7 T to investigate the temporal evolution of metabolic and hemodynamic responses in the rat motor cortex. MethodsFemale Fischer rats received viral injections in the motor cortex to express either a pan-neuronal hM3D(Gq) DREADD construct (hSyn-hM3Dq) or an interneuron-targeted construct (hDlx-hM3Dq). Ph-fMRI, fMRS, and 1H-MRS measurements were performed before, during, and following systemic administration of clozapine-N-oxide (CNO, 1 mg/kg). Functional MRS was acquired during the acute response phase (0-60 min post-injection), while conventional 1H-MRS measurements were obtained at a delayed time point (70 min post-injection). ResultsChemogenetic modulation produced robust and opposing hemodynamic responses. Pan-neuronal activation elicited focal positive BOLD responses (+3.5 {+/-} 1.5%), whereas interneuron-targeted activation generated significant negative BOLD responses (-3.3 {+/-} 0.8%). In contrast, acute fMRS measurements revealed no significant changes in Glx or GABA concentrations during the first hour following CNO administration, despite the presence of strong hemodynamic effects. However, delayed metabolic alterations were detected 70 min after CNO administration. Animals expressing the pan-neuronal construct exhibited significant increases in GABA (+14.4%) and total choline compounds (+57.8%), whereas interneuron-targeted animals displayed reductions in several metabolites, including Glx (-15.6%), total NAA (-16.9%), glucose (-25.9%), and total creatine (-25.4%). ConclusionChemogenetic perturbation of cortical circuits produced robust hemodynamic responses but more subtle and temporally complex metabolic effects. The absence of detectable acute changes in Glx and GABA despite strong BOLD responses, together with the emergence of delayed neurochemical alterations, highlights the challenges of interpreting metabolic signals in relation to circuit activity.

Purpose: To investigate sex differences in cerebral blood flow through densely parcellated cortical and subcortical regions using explainable artificial intelligence methods and identify neurobiologically interpretable perfusion biomarkers. Methods: High-resolution pseudo-continuous arterial spin labelling (1.875 mm x 1.875 mm x 3 mm) and structural MRI data were curated from 215 healthy young adults (150 females, 95 males; age 18-30 years) from the publicly available I See your Brains (ISYB) dataset. Cerebral blood flow was quantified using atlas-based regional analysis with the Brainnetome Atlas (246 regions) and optimized registration procedures. Sex classification employed diverse machine learning paradigms including linear classifiers, ensemble methods, and kernel-based approaches for regional CBF features, with deep convolutional neural networks (CNN) applied to whole-brain 3D imaging data. Model interpretability was achieved using SHapley Additive exPlanations (SHAP), computed over an ensemble of 500 logistic regression models (100 iterations x 5-fold cross-validation). Regions appearing among the top 20% of discriminative features more than 289 times were considered statistically significant using binomial testing. GradCAM was used to obtain class-specific attribution maps from the CNN model. Results: Perfusion-based features demonstrated superior sex classification performance compared to structural morphometry. Regional CBF analysis using logistic regression achieved 91 +/- 2% balanced accuracy and 0.95 +/- 0.05 ROC-AUC, substantially outperforming morphometric features (85 +/- 8% balanced accuracy, 0.88 +/- 0.06 ROC-AUC). Deep learning classification of 3D CBF maps achieved a performance of 92 +/- 5% balanced accuracy, 0.92 +/- 0.05 ROC-AUC. SHAP analysis identified 30 statistically significant aggregation-agnostic CBF-based biomarker regions using regional CBF, predominantly involving frontoparietal control networks (27%) and default mode networks (17%). Grad-CAM revealed that the 3D CNN model primarily focused on regions within the frontal lobe. Morphometry-based analysis identified 28 discriminative regions with markedly different anatomical distribution (r = 0.21) emphasizing visual (32%) and default mode (14%) networks. Conclusion: Cerebral blood flow patterns provide highly sensitive and biologically interpretable markers of sex differences in young adult brain. The identification of robust perfusion biomarkers through explainable AI demonstrates the clinical potential of ASL imaging for precision medicine applications in neuroscience. We establish a methodological framework for investigating sex-specific brain physiology using non-invasive neuroimaging.

We present a non-invasive approach for continuous monitoring of lactate dynamics in-vivo using near-infrared (NIR) spectroscopy. Lactate-related spectral features were measured non-invasively within the overtone region (1600-1850 nm). Several anatomical measurement sites were evaluated, and the middle phalanx of the dorsal finger emerged as the most promising location due to its superior spectral quality and stable tissue perfusion, becoming the exclusive site for all further experiments. Across multiple exercise sessions, predictive models achieved high within-day accuracy (R2[≥] 0.8), while cross-day performance was affected by spectral drift and physiological variability. A dynamic offset-correction procedure effectively mitigated these baseline shifts, enabling stable prediction accuracy across days, weeks, and subjects. These findings demonstrate the feasibility of NIR-based lactate estimation and highlight the importance of adaptive correction strategies for reliable long-term, non-invasive monitoring.

Accurate assessment of cerebral hemodynamics impairment traditionally relies on arterial metrics, yet often overlooks venous drainage and arteriovenous dynamics, thereby limiting the evaluation of ischemia-induced microvascular dysfunction. To address this limitation, we implemented a signal-averaging framework, combined with co-fluctuation analysis, to extract predominantly arterial and venous hemodynamic signals and construct a dynamic arteriovenous co-fluctuation index that quantifies frame-by-frame coordination between arterial inflow and venous outflow activity. This time-resolved index enables spatial characterization of large-scale cortical arteriovenous coordination beyond conventional static correlation-based analyses. Comparative analyses between healthy controls and acute ischemic stroke mice demonstrated that the arteriovenous co-fluctuation index sensitively detects disruption of vascular coordination, revealing a slower state transition that occurs alongside distinct temporal abnormalities and regional heterogeneity between ischemic core and penumbral regions. These findings underscore the utility of arteriovenous coordination as a sensitive indicator of microcirculatory dysfunction, offering a practical analytical tool for assessing stroke-induced microvascular impairment.

Short-wave infrared (SWIR) sensors promise to expand the capabilities of optical sensing technologies but the lack of robust data characterising tissue-constituent optical properties in the SWIR makes instrument design challenging. We characterise and evaluate the optical properties of the dominant chromophores in tissue and tissue-mimicking phantoms, from visible to SWIR wavelengths. Using single-integrating sphere systems, we measured the optical properties of single-component chromophores (H2O, haemoglobin, corn oil, synthetic melanin) and multi-component tissues (whole blood, lard), to decouple contributions from optical scattering, H2O absorption and other contributing chromophores; we also characterised commonly-used phantom materials and investigated their potential to mimic soft tissues in the SWIR range using simulations. We provide a consistent dataset of absorption and reduced scattering coefficients that characterise the dominant tissue chromophores from 450 nm out to 1600 nm. These results were shown to be consistent with literature data, where available. We integrate these data into an open-source Python toolkit, SIMPA, for optical modelling and demonstrate soft tissue simulations that can be probed continuously from visible to SWIR wavelengths. Our findings are compared with tissue-mimicking phantoms, highlighting a need for additives for polymer-based phantoms that mimic SWIR water absorption. By providing this open-source dataset, we aim to enable future studies exploring SWIR light-tissue interactions that facilitate rapid assessment and prototyping of next-generation spectroscopy and imaging techniques.

Understanding brain-behavior relationships requires models capturing the distributed, interactive, and multiscale nature of neural systems. Traditional univariate approaches and single-biomarker models are inherently limited in this context, as they fail to represent dependencies across regions and the hierarchical organization of brain networks. In this study, we propose a graph-based multivariable framework for brain imaging analysis that integrates key organizational principles of brain function-including segregation, integration, modularity, and temporal dynamics-within a unified graph neural network architecture. The framework represents brain data as hierarchical graphs, where node features encode regional activation and temporal variability, and graph structure captures interactions within and between functional modules. The proposed approach is evaluated using functional near-infrared spectroscopy (fNIRS) data as a case study, where subject-specific brain graphs are constructed from task-based recordings acquired shortly after cochlear implant activation to predict speech understanding outcomes one year later. Under leave-one-subject-out validation, the model demonstrates strong predictive performance (R = 0.73, p < 0.001), outperforming previously reported single-biomarker approaches. Perturbation-based analyses further show that predictions are driven by distributed patterns of activity and interaction across regions and modalities, rather than isolated features. These results illustrate the capability of the proposed framework to capture complex brain organization and highlight its potential as a generalizable platform for multivariable analysis and prediction in neuroimaging applications beyond the specific clinical use case considered here.

Low-cost and open-source neural recording systems are increasingly important for expanding access to electrophysiological research. However, many existing platforms still rely on specialized hardware or limited modularity, restricting flexibility for laboratories seeking customizable solutions. Here, we developed and evaluated a modular neural recording platform constructed entirely from commercially available components. Recordings were compared against the ground truth. The platform successfully recovered local field potential (LFP)-like waveforms in most conditions and detected spike-like activity during direct connection recordings. Principal component analysis and k-means clustering further demonstrated the ability to distinguish multiple simulated spike waveforms. Signal quality varied across configurations, with saline recordings and preamplifier integration introducing increased noise and reduced detectability. These findings demonstrate the feasibility of building affordable and modular electrophysiology systems using widely accessible hardware. Although the current implementation has limitations in sampling rate, noise performance, and in vivo validation, the presented framework provides a practical foundation for future customizable open-source neural recording.

Clinical treatment planning of near-infrared (NIR) brain stimulation requires patient-specific light dosimetry to optimize fluence delivery to cortical targets. The gold-standard Monte Carlo (MC) photon transport forward solver is accurate but computationally expensive and non-differentiable for personalized inverse design across subjects. Here, we present a foundation-model (FM)-encoded, differentiable implicit-neural surrogate for the MC solver. A pretrained 3D MRI/CT foundation model, VISTA3D, is domain-adapted to head phantoms with known optical properties to encode the subject anatomy. Next, an implicit neural representation is used to predict light fluence at arbitrary continuous coordinates. This formulation enables off-grid queries and gradients with respect to illumination parameters. Trained with a physics-informed, decade-stratified loss, the surrogate attains R2 {approx} 0.90 on held-out subjects. Ablation results show that the FM benefit is contingent on domain adaptation. Benchmarked against standard learned surrogates, our model is the most accurate in the high-dose region and best on dose-fidelity metrics ({gamma}-index, treated-volume DICE). Finally, gradient-based optimization through the surrogate recovers MC-consistent illumination configurations 50-240 x faster.

Understanding movement encoding within human cortical circuits has been essential for advancing brain computer interfaces (BCIs). However, there are limited minimally invasive, high resolution neurorecording methods sensitive enough to detect single-trial movement-correlated neural activity. Functional ultrasound imaging (fUSI) provides submillimeter spatial resolution of deep cortical tissue with high sensitivity and, when paired with acoustically transparent skull implants, enables transcutaneous recording of human neurovascular changes. Prior studies have used fUSI in participants with acoustically transparent skull implants for on-off task mapping and decoding. Here, we demonstrate fUSI's ability to reliably resolve multi-body-part and single digit movement encoding within the primary sensorimotor cortex in a participant with an acoustically transparent skull implant. We obtained fine-grained mappings of individual effector representation that were consistent with classic somatotopy for both multi-body-part and single digit movement. We were able to resolve single-trial event-related activity, enabling single-trial decoding of both conditions. Analysis of voxels important for decoding suggested differential encoding of single digit movement information across the different Brodmann areas. Finally, we show that these patterns can be approximated across different sessions, allowing for cross session decoding. These results establish that fUSI can reliably delineate somatotopically organized motor representations at submillimeter resolution, bridging a critical gap between invasive electrophysiology and noninvasive hemodynamic imaging in a human subject.