Urban Marhl
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Urban Marhl has a total of 18 co-authors. Here are the top co-authors by number of shared publications:
- Vojko Jazbinšek
Institute of Mathematics, Physics and Mechanics, 1000 Ljubljana, Slovenia
7 shared publications - Tilmann Sander
Physikalisch-Technische Bundesanstalt, 10587 Berlin, Germany
6 shared publications - Anna Jodko-Władzińska
Warsaw University of Technology
3 shared publications - Rok Hren
Syreon Research Institute, 1142 Budapest, Hungary
3 shared publications Vojko Jazbinšek
Affiliation Institute of Mathematics, Physics and Mechanics, Ljubljana, Slovenia
3 shared publications- Adam Liebert
Nalecz Institute of Biocybernetics and Biomedical Engineering, Polish Academy of Sciences, Warsaw, Poland
2 shared publications
Latest Publications
- Open Access
- 2 September 2026
- Published by MDPI AG in Biosensors
- Vol. 16 (9) , 487
Fetal Magnetocardiography Using Optically Pumped Magnetometers: A Literature Review
Fetal magnetocardiography (fMCG) provides direct non-invasive assessment of fetal cardiac electrophysiology, enabling detailed evaluation of cardiac rhythm, conduction, and repolarization. However, the clinical adoption of conventional fMCG has been limited by its reliance on superconducting quantum interference device (SQUID) systems, which require cryogenic cooling and specialized infrastructure. Optically pumped magnetometers (OPMs) have emerged as a promising cryogen-free alternative with the potential to broaden access to fetal electrophysiological assessment. This review summarizes the technological evolution and early clinical evaluation of OPM-based fMCG through an analysis of original in vivo human studies published up to June 2026. Twelve eligible studies were identified and synthesized narratively. Advances in sensor design, magnetic shielding, acquisition strategies, and signal-processing algorithms have enabled SQUID-comparable signal quality and cardiac interval measurements while substantially reducing cryogenic and infrastructure requirements. OPM-fMCG has demonstrated the potential to assess fetal cardiac time intervals, heart rate variability, fetal movement, and clinically important arrhythmias, including congenital long QT syndrome, atrioventricular block, and supraventricular and ventricular tachyarrhythmias. However, the available evidence remains dominated by small, single-centre studies, with relatively few fetuses affected by clinically significant arrhythmias. Prospective multicenter clinical validation, protocol standardization, independent replication, and regulatory evaluation are therefore required before OPM-fMCG can be integrated into routine diagnostic pathways for pregnancies requiring advanced fetal electrophysiological assessment.
- 2 September 2026
- Published by MDPI AG in Biosensors
- Vol. 16 (9) , 487
Open AccessFetal Magnetocardiography Using Optically Pumped Magnetometers: A Literature Review
Fetal magnetocardiography (fMCG) provides direct non-invasive assessment of fetal cardiac electrophysiology, enabling detailed evaluation of cardiac rhythm, conduction, and repolarization. However, the clinical adoption of conventional fMCG has been limited by its reliance on superconducting quantum interference device (SQUID) systems, which require cryogenic cooling and specialized infrastructure. Optically pumped magnetometers (OPMs) have emerged as a promising cryogen-free alternative with the potential to broaden access to fetal electrophysiological assessment. This review summarizes the technological evolution and early clinical evaluation of OPM-based fMCG through an analysis of original in vivo human studies published up to June 2026. Twelve eligible studies were identified and synthesized narratively. Advances in sensor design, magnetic shielding, acquisition strategies, and signal-processing algorithms have enabled SQUID-comparable signal quality and cardiac interval measurements while substantially reducing cryogenic and infrastructure requirements. OPM-fMCG has demonstrated the potential to assess fetal cardiac time intervals, heart rate variability, fetal movement, and clinically important arrhythmias, including congenital long QT syndrome, atrioventricular block, and supraventricular and ventricular tachyarrhythmias. However, the available evidence remains dominated by small, single-centre studies, with relatively few fetuses affected by clinically significant arrhythmias. Prospective multicenter clinical validation, protocol standardization, independent replication, and regulatory evaluation are therefore required before OPM-fMCG can be integrated into routine diagnostic pathways for pregnancies requiring advanced fetal electrophysiological assessment.
- Open Access
Optimizing OPM-MEG Sensor Layouts Using the Sequential Selection Algorithm with Simulated Sources and Individual Anatomy
Magnetoencephalography (MEG) based on optically pumped magnetometers (OPMs) offers the flexibility to position sensors closer to the scalp, which improves the signal-to-noise ratio compared to conventional superconducting quantum interference device (SQUID) systems. However, the spatial resolution of OPM-MEG critically depends on sensor placement, especially when the number of sensors is limited. In this study, we present a methodology for optimizing OPM-MEG sensor layouts using each subject’s anatomical information derived from individual magnetic resonance imaging (MRI). We generated realistic forward models from reconstructed head surfaces and simulated magnetic fields produced by equivalent current dipoles (ECDs). We compared multiple simulation strategies, including ECDs randomly distributed across the cortical surface and ECDs constrained to regions of interest. For each simulated magnetic field map (MFM) database, we applied the sequential selection algorithm (SSA) to identify sensor positions that maximized information capture. Unlike previous approaches relying on large measurement databases, this simulation-driven strategy eliminates the need for extensive pre-existing recordings. We benchmarked the performance of the personalized layouts using OPM-MEG datasets of auditory evoked fields (AEFs) derived from real whole-head SQUID-MEG measurements. Our results show that simulation-based SSA optimization improves the coverage of cortical regions of interest, reduces the number of sensors required for accurate source reconstruction, and yields sensor configurations that perform comparably to layouts optimized using measured data. To evaluate the quality of estimated MFMs, we applied metrics such as the correlation coefficient (CC), root-mean-square error, and relative error. Our results show that the first 15 to 20 optimally selected sensors (CC > 0.95) capture most of the information contained in full-head MFMs. Additionally, we performed source localization for the highest auditory response (M100) by fitting equivalent current dipoles and found that localization errors were < 5 mm. The results further indicate that SSA performance is insensitive to individualized head geometry, supporting the feasibility of using representative anatomical models and highlighting the potential of this approach for clinical OPM-MEG applications.
Open AccessOptimizing OPM-MEG Sensor Layouts Using the Sequential Selection Algorithm with Simulated Sources and Individual Anatomy
Magnetoencephalography (MEG) based on optically pumped magnetometers (OPMs) offers the flexibility to position sensors closer to the scalp, which improves the signal-to-noise ratio compared to conventional superconducting quantum interference device (SQUID) systems. However, the spatial resolution of OPM-MEG critically depends on sensor placement, especially when the number of sensors is limited. In this study, we present a methodology for optimizing OPM-MEG sensor layouts using each subject’s anatomical information derived from individual magnetic resonance imaging (MRI). We generated realistic forward models from reconstructed head surfaces and simulated magnetic fields produced by equivalent current dipoles (ECDs). We compared multiple simulation strategies, including ECDs randomly distributed across the cortical surface and ECDs constrained to regions of interest. For each simulated magnetic field map (MFM) database, we applied the sequential selection algorithm (SSA) to identify sensor positions that maximized information capture. Unlike previous approaches relying on large measurement databases, this simulation-driven strategy eliminates the need for extensive pre-existing recordings. We benchmarked the performance of the personalized layouts using OPM-MEG datasets of auditory evoked fields (AEFs) derived from real whole-head SQUID-MEG measurements. Our results show that simulation-based SSA optimization improves the coverage of cortical regions of interest, reduces the number of sensors required for accurate source reconstruction, and yields sensor configurations that perform comparably to layouts optimized using measured data. To evaluate the quality of estimated MFMs, we applied metrics such as the correlation coefficient (CC), root-mean-square error, and relative error. Our results show that the first 15 to 20 optimally selected sensors (CC > 0.95) capture most of the information contained in full-head MFMs. Additionally, we performed source localization for the highest auditory response (M100) by fitting equivalent current dipoles and found that localization errors were < 5 mm. The results further indicate that SSA performance is insensitive to individualized head geometry, supporting the feasibility of using representative anatomical models and highlighting the potential of this approach for clinical OPM-MEG applications.
- Open Access
- 10 December 2025
- Published by Elsevier in NeuroImage
- Vol. 325, 121644
Evaluating task-evoked neurovascular coupling using integrated OPM-MEG and fNIRS imaging
- 10 December 2025
- Published by Elsevier in NeuroImage
- Vol. 325, 121644
Open AccessEvaluating task-evoked neurovascular coupling using integrated OPM-MEG and fNIRS imaging
- Open Access
Optimization of OPM-MEG Layouts with a Limited Number of Sensors
Magnetoencephalography (MEG) is a non-invasive neuroimaging technique that measures weak magnetic fields generated by neural electrical activity in the brain. Traditional MEG systems use superconducting quantum interference device (SQUID) sensors, which require cryogenic cooling and employ a dense array of sensors to capture magnetic field maps (MFMs) around the head. Recent advancements have introduced optically pumped magnetometers (OPMs) as a promising alternative. Unlike SQUIDs, OPMs do not require cooling and can be placed closer to regions of interest (ROIs). This study aims to optimize the layout of OPM-MEG sensors, maximizing information capture with a limited number of sensors. We applied a sequential selection algorithm (SSA), originally developed for body surface potential mapping in electrocardiography, which requires a large database of full-head MFMs. While modern OPM-MEG systems offer full-head coverage, expected future clinical use will benefit from simplified procedures, where handling a lower number of sensors is easier and more efficient. To explore this, we converted full-head SQUID-MEG measurements of auditory-evoked fields (AEFs) into OPM-MEG layouts with 80 sensor sites. System conversion was done by calculating a current distribution on the brain surface using minimum norm estimation (MNE). We evaluated the SSA’s performance under different protocols, for example, using measurements of single or combined OPM components. We assessed the quality of estimated MFMs using metrics, such as the correlation coefficient (CC), root-mean-square error, and relative error. Additionally, we performed source localization for the highest auditory response (M100) by fitting equivalent current dipoles. Our results show that the first 15 to 20 optimally selected sensors (CC > 0.95, localization error < 1 mm) capture most of the information contained in full-head MFMs. Our main finding is that for event-related fields, such as AEFs, which primarily originate from focal sources, a significantly smaller number of sensors than currently used in conventional MEG systems is sufficient to extract relevant information.
Open AccessOptimization of OPM-MEG Layouts with a Limited Number of Sensors
Magnetoencephalography (MEG) is a non-invasive neuroimaging technique that measures weak magnetic fields generated by neural electrical activity in the brain. Traditional MEG systems use superconducting quantum interference device (SQUID) sensors, which require cryogenic cooling and employ a dense array of sensors to capture magnetic field maps (MFMs) around the head. Recent advancements have introduced optically pumped magnetometers (OPMs) as a promising alternative. Unlike SQUIDs, OPMs do not require cooling and can be placed closer to regions of interest (ROIs). This study aims to optimize the layout of OPM-MEG sensors, maximizing information capture with a limited number of sensors. We applied a sequential selection algorithm (SSA), originally developed for body surface potential mapping in electrocardiography, which requires a large database of full-head MFMs. While modern OPM-MEG systems offer full-head coverage, expected future clinical use will benefit from simplified procedures, where handling a lower number of sensors is easier and more efficient. To explore this, we converted full-head SQUID-MEG measurements of auditory-evoked fields (AEFs) into OPM-MEG layouts with 80 sensor sites. System conversion was done by calculating a current distribution on the brain surface using minimum norm estimation (MNE). We evaluated the SSA’s performance under different protocols, for example, using measurements of single or combined OPM components. We assessed the quality of estimated MFMs using metrics, such as the correlation coefficient (CC), root-mean-square error, and relative error. Additionally, we performed source localization for the highest auditory response (M100) by fitting equivalent current dipoles. Our results show that the first 15 to 20 optimally selected sensors (CC > 0.95, localization error < 1 mm) capture most of the information contained in full-head MFMs. Our main finding is that for event-related fields, such as AEFs, which primarily originate from focal sources, a significantly smaller number of sensors than currently used in conventional MEG systems is sufficient to extract relevant information.
- 1 June 2024
- Published by Springer Nature
- p. 170-177
Combined Measurement of Brain Activation During a Motor Task Using fNIRS and OPM-MEG
- 1 June 2024
- Published by Springer Nature
- p. 170-177
Combined Measurement of Brain Activation During a Motor Task Using fNIRS and OPM-MEG
- Open Access
Simulation Study of Different OPM-MEG Measurement Components
Magnetoencephalography (MEG) is a neuroimaging technique that measures the magnetic fields of the brain outside of the head. In the past, the most suitable magnetometer for MEG was the superconducting quantum interference device (SQUID), but in recent years, a new type has also been used, the optically pumped magnetometer (OPM). OPMs can be configured to measure multiple directions of magnetic field simultaneously. This work explored whether combining multiple directions of the magnetic field lowers the source localization error of brain sources under various conditions of noise. We simulated dipolar-like sources for multiple configurations of both SQUID- and OPM-MEG systems. To test the performance of a given layout, we calculated the average signal-to-noise ratio and the root mean square of the simulated magnetic field; furthermore, we evaluated the performance of the dipole fit. The results showed that the field direction normal to the scalp yields a higher signal-to-noise ratio and that ambient noise has a much lower impact on its localization error; therefore, this is the optimal choice for source localization when only one direction of magnetic field can be measured. For a low number of OPMs, combining multiple field directions greatly improves the source localization results. Lastly, we showed that MEG sensors that can be placed closer to the brain are more suitable for localizing deeper sources.
Open AccessSimulation Study of Different OPM-MEG Measurement Components
Magnetoencephalography (MEG) is a neuroimaging technique that measures the magnetic fields of the brain outside of the head. In the past, the most suitable magnetometer for MEG was the superconducting quantum interference device (SQUID), but in recent years, a new type has also been used, the optically pumped magnetometer (OPM). OPMs can be configured to measure multiple directions of magnetic field simultaneously. This work explored whether combining multiple directions of the magnetic field lowers the source localization error of brain sources under various conditions of noise. We simulated dipolar-like sources for multiple configurations of both SQUID- and OPM-MEG systems. To test the performance of a given layout, we calculated the average signal-to-noise ratio and the root mean square of the simulated magnetic field; furthermore, we evaluated the performance of the dipole fit. The results showed that the field direction normal to the scalp yields a higher signal-to-noise ratio and that ambient noise has a much lower impact on its localization error; therefore, this is the optimal choice for source localization when only one direction of magnetic field can be measured. For a low number of OPMs, combining multiple field directions greatly improves the source localization results. Lastly, we showed that MEG sensors that can be placed closer to the brain are more suitable for localizing deeper sources.
- Open Access
- 19 January 2022
- Published by Public Library of Science (PLoS) in PLOS ONE
- Vol. 17 (1) , e0262669
Transforming and comparing data between standard SQUID and OPM-MEG systems
Optically pumped magnetometers (OPMs) have recently become so sensitive that they are suitable for use in magnetoencephalography (MEG). These sensors solve operational problems of the current standard MEG, where superconducting quantum interference device (SQUID) gradiometers and magnetometers are being used. The main advantage of OPMs is that they do not require cryogenics for cooling. Therefore, they can be placed closer to the scalp and are much easier to use. Here, we measured auditory evoked fields (AEFs) with both SQUID- and OPM-based MEG systems for a group of subjects to better understand the usage of a limited sensor count OPM-MEG. We present a theoretical framework that transforms the within subject data and equivalent simulation data from one MEG system to the other. This approach works on the principle of solving the inverse problem with one system, and then using the forward model to calculate the magnetic fields expected for the other system. For the source reconstruction, we used a minimum norm estimate (MNE) of the current distribution. Two different volume conductor models were compared: the homogeneous conducting sphere and the three-shell model of the head. The transformation results are characterized by a relative error and cross-correlation between the measured and the estimated magnetic field maps of the AEFs. The results for both models are encouraging. Since some commercial OPMs measure multiple components of the magnetic field simultaneously, we additionally analyzed the effect of tangential field components. Overall, our dual-axis OPM-MEG with 15 sensors yields similar information to a 62-channel SQUID-MEG with its field of view restricted to the right hemisphere.
- 19 January 2022
- Published by Public Library of Science (PLoS) in PLOS ONE
- Vol. 17 (1) , e0262669
Open AccessTransforming and comparing data between standard SQUID and OPM-MEG systems
Optically pumped magnetometers (OPMs) have recently become so sensitive that they are suitable for use in magnetoencephalography (MEG). These sensors solve operational problems of the current standard MEG, where superconducting quantum interference device (SQUID) gradiometers and magnetometers are being used. The main advantage of OPMs is that they do not require cryogenics for cooling. Therefore, they can be placed closer to the scalp and are much easier to use. Here, we measured auditory evoked fields (AEFs) with both SQUID- and OPM-based MEG systems for a group of subjects to better understand the usage of a limited sensor count OPM-MEG. We present a theoretical framework that transforms the within subject data and equivalent simulation data from one MEG system to the other. This approach works on the principle of solving the inverse problem with one system, and then using the forward model to calculate the magnetic fields expected for the other system. For the source reconstruction, we used a minimum norm estimate (MNE) of the current distribution. Two different volume conductor models were compared: the homogeneous conducting sphere and the three-shell model of the head. The transformation results are characterized by a relative error and cross-correlation between the measured and the estimated magnetic field maps of the AEFs. The results for both models are encouraging. Since some commercial OPMs measure multiple components of the magnetic field simultaneously, we additionally analyzed the effect of tangential field components. Overall, our dual-axis OPM-MEG with 15 sensors yields similar information to a 62-channel SQUID-MEG with its field of view restricted to the right hemisphere.
- Open Access
- 1 October 2021
- Published by Walter de Gruyter GmbH in Current Directions in Biomedical Engineering
- Vol. 7 (2) , 543-546
Avoiding non-linearity of optically pumped magnetometer MEG within an actively shielded two-layer mu-metal room
Some optically pumped magnetometer (OPM) sensors available for biomagnetic investigations have a linear range limited to +- 1 nT due to the specific properties of their open loop operation. In a two-layer magnetically shielded room of type Ak3b/Vacoshield Advanced with an external active compensation we studied how much sensor movement is allowed until amplitudes exceed the linearity range. Intentional movements were performed by a subject wearing an OPM-MEG sensor array. It was found that movements of 8 cm did yield non-linear amplitudes, but a reduction of the movement in half already preserves linearity. Despite movements, the heartbeat was found to generate a periodic signal, although the generating mechanism could not be identified so far.
- 1 October 2021
- Published by Walter de Gruyter GmbH in Current Directions in Biomedical Engineering
- Vol. 7 (2) , 543-546
Open AccessAvoiding non-linearity of optically pumped magnetometer MEG within an actively shielded two-layer mu-metal room
Some optically pumped magnetometer (OPM) sensors available for biomagnetic investigations have a linear range limited to +- 1 nT due to the specific properties of their open loop operation. In a two-layer magnetically shielded room of type Ak3b/Vacoshield Advanced with an external active compensation we studied how much sensor movement is allowed until amplitudes exceed the linearity range. Intentional movements were performed by a subject wearing an OPM-MEG sensor array. It was found that movements of 8 cm did yield non-linear amplitudes, but a reduction of the movement in half already preserves linearity. Despite movements, the heartbeat was found to generate a periodic signal, although the generating mechanism could not be identified so far.
- 30 November 2020
- Published by Springer Nature
- p. 254-261
Comparison Between Conventional SQUID Based and Novel OPM Based Measuring Systems in MEG
- 30 November 2020
- Published by Springer Nature
- p. 254-261
Comparison Between Conventional SQUID Based and Novel OPM Based Measuring Systems in MEG
- 12 December 2019
- Published by American Physical Society (APS) in Physical Review E
- Vol. 100 (6) , 062203
Proper spatial heterogeneities expand the regime of scale-free behavior in a lattice of excitable elements
Signatures of criticality, such as power law scaling of observables, have been empirically found in a plethora of real-life settings, including biological systems. The presence of critical states is believed to have many functional advantages and is associated with optimal operational abilities. Typically, critical dynamics arises in the proximity of phase transition points between absorbing disordered states (subcriticality) and ordered active regimes (supercriticality) and requires a high degree of fine tuning to emerge, which is unlikely to occur in real biological systems. In the present study we propose a rather simple, and biologically relevant mechanism that profoundly expands the critical-like region. In particular, by means of numerical simulation we show that incorporating spatial heterogeneities into the square lattice of map-based excitable oscillators broadens the parameter space in which the distribution of excitation wave sizes follows closely a power law. Most importantly, this behavior is only observed if the spatial profile exhibits intermediate-sized patches with similar excitability levels, whereas for large and small spatial clusters only marginal widening of the critical state is detected. Furthermore, it turned out that the presence of spatial disorder in general amplifies the size of excitation waves, whereby the relatively highest contributions are observed in the proximity of the critical point. We argue that the reported mechanism is of particular importance for excitable systems with local interactions between individual elements.
- 12 December 2019
- Published by American Physical Society (APS) in Physical Review E
- Vol. 100 (6) , 062203
Proper spatial heterogeneities expand the regime of scale-free behavior in a lattice of excitable elements
Signatures of criticality, such as power law scaling of observables, have been empirically found in a plethora of real-life settings, including biological systems. The presence of critical states is believed to have many functional advantages and is associated with optimal operational abilities. Typically, critical dynamics arises in the proximity of phase transition points between absorbing disordered states (subcriticality) and ordered active regimes (supercriticality) and requires a high degree of fine tuning to emerge, which is unlikely to occur in real biological systems. In the present study we propose a rather simple, and biologically relevant mechanism that profoundly expands the critical-like region. In particular, by means of numerical simulation we show that incorporating spatial heterogeneities into the square lattice of map-based excitable oscillators broadens the parameter space in which the distribution of excitation wave sizes follows closely a power law. Most importantly, this behavior is only observed if the spatial profile exhibits intermediate-sized patches with similar excitability levels, whereas for large and small spatial clusters only marginal widening of the critical state is detected. Furthermore, it turned out that the presence of spatial disorder in general amplifies the size of excitation waves, whereby the relatively highest contributions are observed in the proximity of the critical point. We argue that the reported mechanism is of particular importance for excitable systems with local interactions between individual elements.