Significance: Jacobians, or spatially resolved sensitivity profiles, are central to image reconstruction in model-based optical tomography of biological tissue. Although Monte Carlo (MC) simulations are the gold standard for modeling light transport in turbid media, methodology for frequency- and time-domain Jacobians remains incomplete. Aim: This work extends MC to directly compute absorption and scattering Jacobians for frequency-domain (amplitude and phase) and time-domain (intensity and mean time-of-flight) measurements and prism-terminated optical fiber detectors. Approach: Jacobians are derived in the perturbation MC framework and implemented in the high-performance, open-source Monte Carlo eXtreme (MCX) simulator. Results are validated against the diffusion approximation (DA) solved using the finite element method in neonatal head models. MC with split voxels on curved surfaces is extended to Jacobian computation. The detector model is implemented in post-processing and compared with isotropic reception at surface. Results: MC- and DA-derived Jacobians show excellent agreement only in high-scattering regimes, highlighting the importance of MC for low-scattering domains. The detector model reduces surface sensitivity and marginally increases sensitivity to deeper tissues at short (< 2 cm) source-detector separations. Conclusion: A complete theoretical framework and MC software for computing frequency- and time-domain Jacobians is provided. Realistic detector modeling is encouraged for short-separation channels.

Simultaneous measurement of electrophysiological and hemodynamic brain signals imposes special requirements on the instrumentation. Here, we developed a high-density fiberoptic probe for concurrent diffuse optical tomography (DOT) and magnetoencephalography (MEG) recordings. Transparent two-component silicone was mixed with carbon black dye to achieve a black, flexible, non-magnetic support for the dense optode arrangement and low (5 mm) probe thickness. The probe was used to record somatosensory responses to electrical right median nerve stimulation at 0.5, 1, 2, and 4 Hz in 18 adult human subjects. Brain activity was simultaneously measured with a commercial whole-head MEG system and with the DOT optode arrangement covering approximately 40 cm 2 over the parietal region in the contralateral left hemisphere. Two correlation-based clustering methods were developed to find regions where the reconstructed time course of total hemoglobin concentration (HbT) changes correlated with the predicted hemodynamic activity based on time-course characteristics of the MEG sources and the canonical hemodynamic response model. Two statistically significant clusters were found based on the correlation between HbT around the postcentral gyrus and MEG primary somatosensory cortical activity at ~35 ms (P35m response). In addition, correlation between HbT and secondary somatosensory cortical activity suggested a statistically significant cluster in the postcentral gyrus and parietal operculum. These results illustrate an improvement in localization over previous DOT studies using sparse optode arrangements, and demonstrate the feasibility of the system for simultaneous HD-DOT-MEG experiments. Furthermore, the techniques described here pave the way for understanding the coupling between hemodynamic and electrophysiological responses. Further research is needed to reveal the neuronal circuits giving rise to the correlating MEG and DOT response features. Significant improvements in the technology are still expected via optimization of the detected light power in the instrumentation.

Diffuse optical tomography uses near-infrared light for 3D brain imaging. We simulated errors due to inaccurate baseline parameters and developed a deep learning technique that mitigated these errors, enhancing contrast and localization of activation changes. (c) 2025 The Author(s)

Background: Touch is an essential form of mother-child interaction, instigating better social bonding and emotional stability. Methods: We used diffuse optical tomography to explore the relationship between total haemoglobin (HbT) responses to affective touch in the child's brain at two years of age and maternal self-reported prenatal depressive symptoms (EPDS). Affective touch was implemented via slow brushing of the child's right forearm at 3 cm/s and non-affective touch via fast brushing at 30 cm/s and HbT responses were recorded on the left hemisphere. Results: We discovered a cluster in the postcentral gyrus exhibiting a negative correlation (Pearson's r = -0.84, p = 0.015 corrected for multiple comparisons) between child HbT response to affective touch and EPDS at gestational week 34. Based on region of interest (ROI) analysis, we found negative correlations between child responses to affective touch and maternal prenatal EPDS at gestational week 14 in the precentral gyrus, Rolandic operculum and secondary somatosensory cortex. The responses to non-affective touch did not correlate with EPDS in these regions. Limitations: The number of mother-child dyads was 16. However, by utilising high-density optode arrangements, individualised anatomical models, and video and accelerometry to monitor movement, we were able to minimize methodological sources of variability in the data. Conclusions: The results show that maternal depressive symptoms during pregnancy may be associated with reduced child responses to affective touch in the temporoparietal cortex. Responses to affective touch may be considered as potential biomarkers for psychosocial development in children. Early identification of and intervention in maternal depression may be important already during early pregnancy.

Diffuse optical tomography (DOT) uses near-infrared light to image spatially varying optical parameters in biological tissues. In functional brain imaging, DOT uses a perturbation model to estimate the changes in optical parameters, corresponding to changes in measured data due to brain activity. The perturbation model typically uses approximate baseline optical parameters of the different brain compartments, since the actual baseline optical parameters are unknown. We simulated the effects of these approximate baseline optical parameters using parameter variations earlier reported in literature, and brain atlases from four adult subjects. We report the errors in estimated activation contrast, localization, and area when incorrect baseline values were used. Further, we developed a post-processing technique based on deep learning methods that can reduce the effects due to inaccurate baseline optical parameters. The method improved imaging of brain activation changes in the presence of such errors.

Frequency-domain (FD) optical tomography instruments modulate the intensity of the light source at a radio frequency and measure the amplitude and phase shift of the detected photon density wave. The differing spatial sensitivities of amplitude and phase to the optical properties of tissue suggest that inclusion of phase data can improve the image reconstruction accuracy. This study describes our methodology for improved use of FD data in conjunction with a Monte Carlo (MC) forward solver (Monte Carlo eXtreme; MCX) and a voxel-based model of a two-year-old child’s head. The child participated our previous study where subjects were stimulated with affective (slow brushing) and non-affective touch (fast brushing) to their right forearm, and the responses were measured from the left hemisphere with our in-house 16-channel high-density FD system. We implemented the computation of the FD sensitivity profiles to the MCX photon simulation software, and validated the output against our in-house MC code. We used simulated and the real experimental touch response data to observe the effects of including both FD data types to the image reconstruction instead of amplitude data alone. For the simulated and experimental case, we observed that the inclusion of phase data increases the reconstructed contrast in the brain. The individual touch responses showed similarity to the group-level results in our original publication with 16 subjects and amplitude data alone, and other literature.

Objective.Diffuse optical tomography (DOT) provides a relatively convenient method for imaging haemodynamic changes related to neuronal activity on the cerebral cortex. Due to practical challenges in obtaining anatomical images of neonates, an anatomical framework is often created from an age-appropriate atlas model, which is individualized to the subject based on measurements of the head geometry. This work studies the approximation error arising from using an atlas instead of the neonate's own anatomical model.Approach.We consider numerical simulations of frequency-domain (FD) DOT using two approaches, Monte Carlo simulations and diffusion approximation via finite element method, and observe the variation in (1) the logarithm of amplitude and phase shift measurements, and (2) the corresponding inner head sensitivities (Jacobians), due to varying segmented anatomy. Varying segmentations are sampled by registering 165 atlas models from a neonatal database to the head geometry of one individual selected as the reference model. Prior to the registration, we refine the segmentation of the cerebrospinal fluid (CSF) by separating the CSF into two physiologically plausible layers.Main results.In absolute measurements, a considerable change in the grey matter or extracerebral tissue absorption coefficient was found detectable over the anatomical variation. In difference measurements, a small local 10%-increase in brain absorption was clearly detectable in the simulated measurements over the approximation error in the Jacobians, despite the wide range of brain maturation among the registered models.Significance.Individual-level atlas models could potentially be selected within several weeks in gestational age in DOT difference imaging, if an exactly age-appropriate atlas is not available. The approximation error method could potentially be implemented to improve the accuracy of atlas-based imaging. The presented CSF segmentation algorithm could be useful also in other model-based imaging modalities. The computation of FD Jacobians is now available in the widely-used Monte Carlo eXtreme software.

Diffuse optical tomography (DOT) utilises near-infrared light for imaging spatially distributed optical parameters, typically the absorption and scattering coefficients. The image reconstruction problem of DOT is an ill-posed inverse problem, due to the non-linear light propagation in tissues and limited boundary measurements. The ill-posedness means that the image reconstruction is sensitive to measurement and modelling errors. The Bayesian approach for the inverse problem of DOT offers the possibility of incorporating prior information about the unknowns, rendering the problem less ill-posed. It also allows marginalisation of modelling errors utilising the so-called Bayesian approximation error method. A more recent trend in image reconstruction techniques is the use of deep learning, which has shown promising results in various applications from image processing to tomographic reconstructions. In this work, we study the non-linear DOT inverse problem of estimating the (absolute) absorption and scattering coefficients utilising a 'model-based' learning approach, essentially intertwining learned components with the model equations of DOT. The proposed approach was validated with 2D simulations and 3D experimental data. We demonstrated improved absorption and scattering estimates for targets with a mix of smooth and sharp image features, implying that the proposed approach could learn image features that are difficult to model using standard Gaussian priors. Furthermore, it was shown that the approach can be utilised in compensating for modelling errors due to coarse discretisation enabling computationally efficient solutions. Overall, the approach provided improved computation times compared to a standard Gauss-Newton iteration.

Touch is an important component of early parent-child interaction and plays a critical role in the socio-emotional development of children. However, there are limited studies on touch processing amongst children in the age range from one to three years. The present study used frequency-domain diffuse optical tomography (DOT) to investigate the processing of affective and non-affective touch over left frontotemporal brain areas contralateral to the stimulated forearm in two-year-old children. Affective touch was administered by a single stroke with a soft brush over the child's right dorsal forearm at 3 cm/s, while non-affective touch was provided by multiple brush strokes at 30 cm/s. We found that in the insula, the total haemoglobin (HbT) response to slow brushing was significantly greater than the response to fast brushing (slow > fast). Additionally, a region in the postcentral gyrus, Rolandic operculum and superior temporal gyrus exhibited greater response to fast brushing than slow brushing (fast > slow). These findings confirm that an adult-like pattern of haemodynamic responses to affective and non-affective touch can be recorded in two-year-old subjects using DOT. To improve the accuracy of modelling light transport in the two-year-old subjects, we used a published age-appropriate atlas and deformed it to match the exterior shape of each subject's head. We estimated the combined scalp and skull, and grey matter (GM) optical properties by fitting simulated data to calibrated and coupling error corrected phase and amplitude measurements. By utilizing a two-compartment cerebrospinal fluid (CSF) model, the accuracy of estimation of GM optical properties and the localization of activation in the insula was improved. The techniques presented in this paper can be used to study neural development of children at different ages and illustrate that the technology is well-tolerated by most two-year-old children and not excessively sensitive to subject movement. The study points the way towards exciting possibilities in functional imaging of deeper functional areas near sulci in small children.

Common disc-type ultrasonic atomizers are characterized by a liquid column of millimeters to centimeters above the disc to successfully atomize the liquid. This results in to a fountain and splash formation, along with heating of the liquid, which typically are undesired. In recent studies, this limitation has been overcome by introducing large-sized droplet generation as a pre-stage to the ultrasonic atomization process. This, however, typically requires a pump to feed the atomizer. In the present study, we designed a novel single drop dispenser as a pre-stage which can operate using gravity with timing option, and thereby eliminating the need for a pump for the pre-stage dispensing.

Emotional speech is one of the principal forms of social communication in humans. In this study, we investigated neural processing of emotional speech (happy, angry, sad and neutral) in the left hemisphere of 21 two-month-old infants using diffuse optical tomography. Reconstructed total hemoglobin (HbT) images were analysed using adaptive voxel-based clustering and region-of-interest (ROI) analysis. We found a distributed happy > neutral response within the temporo-parietal cortex, peaking in the anterior temporal cortex; a negative HbT response to emotional speech (the average of the emotional speech conditions < baseline) in the temporo-parietal cortex, neutral > angry in the anterior superior temporal sulcus (STS), happy > angry in the superior temporal gyrus and posterior superior temporal sulcus, angry < baseline in the insula, superior temporal sulcus and superior temporal gyrus and happy < baseline in the anterior insula. These results suggest that left STS is more sensitive to happy speech as compared to angry speech, indicating that it might play an important role in processing positive emotions in two-month-old infants. Furthermore, happy speech (relative to neutral) seems to elicit more activation in the temporo-parietal cortex, thereby suggesting enhanced sensitivity of temporo-parietal cortex to positive emotional stimuli at this stage of infant development.

BACKGROUND:Maternal pregnancy-related anxiety (PRA) is reportedly related to neurodevelopmental outcomes of infants. However, the relationship between maternal PRA and the processing of emotions in the infant brain has not been extensively studied with neuroimaging. The objective of the present pilot study is to investigate the relationship between maternal PRA and infant hemodynamic responses to emotional speech at two months of age. METHODS:The study sample included 19 mother-infant dyads from a general sample of a population of Caucasian mothers. Self-reported Pregnancy-Related Anxiety Questionnaire (PRAQ-R2) data was collected from mothers during pregnancy at gestational weeks (gwks) 24 (N = 19) and 34 (N = 18). When their infants were two months old, the infants' brains functional responses to emotional speech in the left fronto-temporoparietal cortex were recorded using diffuse optical tomography (DOT). RESULTS:Maternal PRAQ-R2 scores at gwk 24 correlated negatively with the total hemoglobin (HbT) responses to sad speech on both sides of the temporoparietal junction (Spearman's rank correlation coefficient ρ = -0.87). The correlation was significantly greater at gwk 24 than gwk 34 (ρ = -0.42). LIMITATIONS:The field of view of the measurement did not include the right hemisphere or parts of the frontal cortex. The sample size is moderate and the mothers were relatively highly educated, thus there may be some differences between the study sample and the general population. CONCLUSIONS:Maternal pregnancy-related anxiety may affect child brain emotion processing development. Further research is needed to understand the functional and developmental significance of the findings.

ABSTRACT Emotional stimuli processing during childhood helps us to detect salient cues in our environment and prepares us for our social life. In early childhood, the emotional valences of auditory and visual input are salient and relevant cues of social aspects of the environment, and it is of special interest to understand how exactly the processing of emotional stimuli develops. Near‐infrared spectroscopy (NIRS) is a noninvasive neuroimaging tool that has proven valuable in studying emotional processing in children. After conducting a systematic search of PubMed, Web of Science, and Embase databases, we examined 50 NIRS studies performed to study emotional stimuli processing in children in the first 2 years of age. We found that the majority of these studies are done in infants and the most commonly used stimuli are visual and auditory. Many of the reviewed studies suggest the involvement of bilateral temporal areas in emotional processing of visual and auditory stimuli. It is unclear which neural activation patterns reflect maturation and at what age the emotional encoding reaches those typically seen in adults. Our review provides an overview of the database on emotional processing in children up to 2 years of age. Furthermore, it demonstrates the need to include the less‐studied age range of 1 to 2 years, and suggests the use of combined audio‐visual stimuli and longitudinal studies for future research on emotional processing in children. Thus, NIRS might be a vital tool to study the associations between the early pattern of neural responses and socioemotional development later in life.

Caressing touch is an effective way to communicate emotions and to create social bonds. It is also one of the key mediators of early parental bonding. The caresses are generally thought to represent a social form of touching and indeed, slow, gentle brushing is encoded in specialized peripheral nerve fibers, the C-tactile (CT) afferents. In adults, areas such as the posterior insula and superior temporal sulcus are activated by affective, slow stroking touch but not by fast stroking stimulation. However, whether these areas are activated in infants, after social tactile stimulation, is unknown. In this study, we compared the total hemoglobin responses measured with diffuse optical tomography (DOT) in the left hemisphere following slow and fast stroking touch stimulation in 16 2-month-old infants. We compared slow stroking (optimal CT afferent stimulation) to fast stroking (non-optimal CT stimulation). Activated regions were delineated using two methods: one based on contrast between the two conditions, and the other based on voxel-based statistical significance of the difference between the two conditions. The first method showed a single activation cluster in the temporal cortex with center of gravity in the middle temporal gyrus where the total hemoglobin increased after the slow stroking relative to the fast stroking (p = 0.04 uncorrected). The second method revealed a cluster in the insula with an increase in total hemoglobin in the insular cortex in response to slow stroking relative to fast stroking (p = 0.0005 uncorrected; p = 0.04 corrected for multiple comparisons). These activation clusters encompass areas that are involved in processing of affective, slow stroking touch in the adult brain. We conclude that the infant brain shows a pronounced and adult-like response to slow stroking touch compared to fast stroking touch in the insular cortex but the expected response in the primary somatosensory cortex was not found at this age. The results imply that emotionally valent touch is encoded in the brain in adult-like manner already soon after birth and this suggests a potential for involvement of touch in bonding with the caretaker.

Difference imaging aims at recovery of the change in the optical properties of a body based on measurements before and after the change.Conventionally, the image reconstruction is based on using difference of the measurements and a linear approximation of the observation model.One of the main benefits of the linearized difference reconstruction is that the approach has a good tolerance to modeling errors, which cancel out partially in the subtraction of the measurements.However, a drawback of the approach is that the difference images are usually only qualitative in nature and their spatial resolution can be weak because they rely on the global linearization of the nonlinear observation model.To overcome the limitations of the linear approach, we investigate a nonlinear approach for difference imaging where the images of the optical parameters before and after the change are reconstructed simultaneously based on the two datasets.We tested the feasibility of the method with simulations and experimental data from a phantom and studied how the approach tolerates modeling errors like domain truncation, optode coupling errors, and domain shape errors.

We used emotional speech to elicit brain activity in two-month-old infants and recorded emotion-modulated total hemoglobin changes on the frontotemporal cortex using diffuse optical imaging with a multi-distance source-detector array.

Near-infrared spectroscopy (NIRS) and diffuse optical tomography (DOT) detect changes in brain blood volume and oxygenation by measuring light that has passed through the head, including the scalp and the skull. Extracerebral and systemic circulation interfere with optical measurements of cerebral hemodynamics, especially when measuring brain responses to stimuli or tasks that evoke strong systemic circulatory changes. We studied the effect of changes in systemic circulation on NIRS responses and DOT reconstructions in thirteen subjects during a hand motor task that increased the heart rate. Both the NIRS responses and the DOT reconstructions depended on the change in the heart rate. The NIRS response amplitudes during epochs with a large change in heart rate (24.8±0.8 bpm; highest third) were significantly larger (p < 0.05) than during epochs with a smaller change in heart rate (5.8±0.5 bpm; lowest third). Accordingly, we propose that comparing epochs associated with large and small changes in heart rate serves as a method for estimating whether NIRS signals are affected by the systemic circulation, given that there is variability in the systemic circulation between epochs.

The effect of task-related extracerebral circulatory changes on diffuse optical tomography (DOT) of brain activation was evaluated using experimental data from 14 healthy human subjects and computer simulations. Total hemoglobin responses to weekday-recitation, verbal-fluency, and hand-motor tasks were measured with a high-density optode grid placed on the forehead. The tasks caused varying levels of mental and physical stress, eliciting extracerebral circulatory changes that the reconstruction algorithm was unable to fully distinguish from cerebral hemodynamic changes, resulting in artifacts in the brain activation images. Crosstalk between intra- and extracranial layers was confirmed by the simulations. The extracerebral effects were attenuated by superficial signal regression and depended to some extent on the heart rate, thus allowing identification of hemodynamic changes related to brain activation during the verbal-fluency task. During the hand-motor task, the extracerebral component was stronger, making the separation less clear. DOT provides a tool for distinguishing extracerebral components from signals of cerebral origin. Especially in the case of strong task-related extracerebral circulatory changes, however, sophisticated reconstruction methods are needed to eliminate crosstalk artifacts.