2026-02-10 | | Total: 152
Neural ordinary differential equations (neural ODEs) can effectively learn dynamical systems from time series data, but their behavior on graph-structured data remains poorly understood, especially when applied to graphs with different size or structure than encountered during training. We study neural ODEs ($\mathtt{nODE}$s) with vector fields following the Barabási-Barzel form, trained on synthetic data from five common dynamical systems on graphs. Using the $\mathbb{S}^1$-model to generate graphs with realistic and tunable structure, we find that degree heterogeneity and the type of dynamical system are the primary factors in determining $\mathtt{nODE}$s' ability to generalize across graph sizes and properties. This extends to $\mathtt{nODE}$s' ability to capture fixed points and maintain performance amid missing data. Average clustering plays a secondary role in determining $\mathtt{nODE}$ performance. Our findings highlight $\mathtt{nODE}$s as a powerful approach to understanding complex systems but underscore challenges emerging from degree heterogeneity and clustering in realistic graphs.
We study the dimensionality of two-dimensional Kolmogorov flows over a wide range of Reynolds numbers and forcing wavenumbers $k_f=\{2,4,8\}$ using two complementary approaches: convolutional autoencoders and a Kaplan-Yorke estimation based on Lyapunov analysis. As the Reynolds number increases, two distinct transitions are observed: the first corresponds to the destabilization of a periodic orbit, while the second marks the saturation of the large-scale motions. When expressed in terms of the forcing Reynolds number, these transitions occur at nearly the same value for all forcing wavenumbers, suggesting a universal scaling with respect to the forcing scale. By filtering the data to retain only the large-scale range ($k < k_f$), we show that the dimensionality estimated by the autoencoders also saturates at the second transition, implying that once the large scales are fully developed, the subsequent increase in dynamical activity occurs predominantly at smaller scales. At higher Reynolds numbers, the Kaplan-Yorke dimension ceases to grow, revealing its limited sensitivity to the nonlinear interactions that dominate in this regime. Both the Kaplan-Yorke saturation dimension and the filtered large-scale dimensionalities exhibit a linear dependence on $k_f$, indicating that the number of active degrees of freedom scales with the forcing scale rather than with the total number of available Fourier modes.
Cooperative collective dynamics is a principal determinant of the ability of synthetic micromotors to perform specific functions. However, realizing controllable and predictable collective behavior in complex physiological environments remains a significant challenge. Here, we show that collections of enzyme-coated colloids can be designed as various chemical logic gates, which subsequently can be organized into functional logic circuits. These circuits take environmental information as input signals and process it to produce output chemical species needed to achieve specific goals. The chemical computation performed by the circuit endows the active colloidal system with the ability to sense its surroundings and autonomously coordinate its collective motion. The results of simulations of several examples are presented, where self-assembled colloidal circuits can identify invasive threats by their signals, produce and deliver chemicals to the targets to suppress their activity. The results of this work can aid in the design of experimental chemical logic circuits through micromotor self-assembly that autonomously respond to environmental cues to execute specific tasks.
I argue that spacetime singularities entail no ontological commitment to material entities. First, I show that Penrose's singularity theorem is best understood as a theorem of incompleteness, it demonstrates the failure of specific spacetime models within General Relativity (or any theory incorporating the Raychaudhuri equation) under certain general conditions. Although this has been done before, I adopt a novel approach based on differentiating between physical and purely formal assumptions in the axiomatic foundation of general relativity. Next, I compare Penrose's result with Gödel's incompleteness theorem, highlighting key similarities and differences. Finally, I draw philosophical conclusions regarding the limits and prospects of our epistemic reconstructions of the physical world.
Dynamic volumetric MRI provides valuable information on in vivo motion and biomechanics, with applications spanning cardiac, musculoskeletal, or pulmonary imaging, amongst others. Developing reconstruction methods for time-resolved volumetric MRI is challenging due to the inherently slow acquisition process of MRI, which makes it an active area of research. However, in vivo validation of these methods remains challenging due to the lack of publicly available datasets with fully sampled ground-truth images. Here, we present a publicly available in vivo dataset designed to facilitate the development and validation of dynamic volumetric MRI reconstruction algorithms. Controlled and repeatable deformations of the muscles in the thigh were induced using a pneumatic pressure cuff, enabling the acquisition of both undersampled dynamic data and fully sampled validation images. The dataset comprises multichannel undersampled k-space data from nine healthy volunteers across four different dynamic deformations, with fully sampled validation data for one deformation. Additionally, an anatomical reference scan and muscle segmentation masks are provided for each subject. To illustrate a possible image reconstruction and validation approach, a binning-based reconstruction was performed on the undersampled data from six dynamic repetitions. The resulting images were consistent with the corresponding fully sampled validation images. This dataset offers possibilities for validating and advancing time-resolved volumetric MRI reconstruction methods.
Due to their interactive nature, serious games offer valuable opportunities for supporting learning in educational contexts. Recent advances in large language models (LLMs) have further opened the door to new forms of personalized scaffolding in education. In this study, we combine both worlds and study three types of AI-based scaffolding designs in a serious game: (i) no scaffolding, (ii) chat-based (verbal) scaffolding provided by an AI-based non-player character (NPC), and (iii) combined chat-(verbal) and action-based (visual) scaffolding in which the AI may both try to explain or demonstrate the next step towards a solution. The scaffolding conditions are embedded in Qookies, a serious game designed to introduce fundamental concepts of quantum technologies. A total of 152 school students, university students, and members of the general public were randomly assigned to one of the three conditions. The results show that all groups experience significant learning gains, confirming the overall effectiveness of the serious game itself. No significant differences in learning outcomes emerged between scaffolding conditions. However, intrinsic cognitive load was lower in the combined chat-and-action (verbal+visual) scaffolding condition compared to the chat (verbal)-only condition, suggesting that visual demonstrations may offer more accessible support. Interaction analyses further revealed that players engaged with the AI character primarily for level-related questions and action recommendations, while deeper interactions were relatively rare.
Advanced tokamak scenarios often feature weak or reversed magnetic shear configurations. In this study, the hybrid kinetic-MHD model implemented in the NIMROD code is used to investigate the effects of reversed magnetic shear on internal kink and fishbone mode in a circular shaped limiter tokamak. In the absence of energetic particles (EPs), the mode growth rate initially increases and then decreases as the magnetic shear changes from positive to negative, indicating stabilizing effects of the reversed magnetic shear on the internal kink mode. In the presence of EPs, when the reversed magnetic shear region is sufficiently narrow, the transition from internal kink/fishbone modes to double kink/fishbone modes takes place, and the stabilizing effects of the reversed magnetic shear can significantly dominate the destabilization of EPs. For non-resonant modes, the EP beta fraction $β_f$ for excitation increases with $q_{min}$, concurrent with progressively lower growth rates in non-resonant fishbone modes. When the equilibrium profile has an internal transport barrier (ITB), broader ITB widths suppress internal kink modes more effectively, whereas steeper temperature gradients strengthen EP stabilization.
As humans make ambitious efforts toward long-duration activities beyond Earth, new challenges will continue to emerge that highlight the need for governance frameworks capable of managing shared resources and technical standards in order to sustain human life in these hostile environments. Earth-based governance models of cooperative sovereignty can inform governance mechanisms for future Mars settlements, particularly regarding inter-settlement relations and the technical coordination required for multiple independent settlements to coexist. This study analyzes the International Telecommunication Union (ITU) and the Universal Postal Union (UPU), two of the oldest international organizations, which have successfully established evolving standards across sovereign nations. This analysis of the development and governance structures of these two organizations, and how they resolved key sovereignty issues, reveals principles that could be applicable to future settlements beyond Earth, particularly on Mars. Key insights include the strategic necessity of institutional neutrality, the management of asymmetric power relations, and the governance of shared resources under conditions of mutual vulnerability. The study distinguishes between a "Survival Layer" of technical standards essential for immediate safety and an "Operational Layer" governing economic and political activities, suggesting different governance approaches for each. Although some of these examples of cooperative sovereignty on Earth might not be sufficient for Mars due to its unique environment, lessons from the ITU and UPU case studies offer valuable strategies for designing flexible and sustainable governance models that can function from inception through explicit Earth-based coordination.
Structured light (SL)-induced polarization perception presents a promising method for the early detection of macular diseases such as age-related macular degeneration. We investigated the test-retest reliability of a psychophysical task using SL-based stimuli to induce entoptic patterns in individuals with healthy vision. Twenty-eight participants underwent thorough eye examinations to confirm they had healthy eyes and good vision (logMAR BCVA ~0.00). Of these, 25 participants (n=50 eyes) aged 21 to 75 completed two identical tasks separated by 1 to 14 days. Using SL-based stimuli that produced a rotating entoptic pattern containing 22 azimuthal brushes, we measured the retinal eccentricity threshold $R_T$ at which participants could reliably identify the direction of rotation by varying the size of a central obstruction. This threshold reflects the visual angle of the pattern for each participant. We calculated the reliability coefficient (intraclass correlation coefficient, ICC) using a two-way mixed-effects model and conducted a Bland-Altman analysis to assess test-retest reliability. The ICC was 0.83 [95% CI: 0.62 - 0.93] for right eyes (RE) and 0.93 [95% CI: 0.84 - 0.97] for left eyes (LE), indicating good reliability. The Bland-Altman analysis showed a mean difference of -0.32° (SD: 1.51) and -0.12° (SD: 1.03) for RE and LE respectively between the first and second sessions, with limits of agreement ranging from -3.28° to 2.64° and -2.14° to 1.89° for RE and LE respectively, confirming strong agreement and no significant bias. These results demonstrate that SL-based psychophysical tasks are a reliable method for assessing polarization perception, potentially improving screening for diseases affecting macular health.
We present a systematic methodology for realizing accurate asymmetric all-angle transmission in nonlocal metasurfaces. As a representative example, we derive closed-form susceptibility conditions for exact first-order spatial differentiation of unity numerical aperture, clarifying the role of each underlying balance. We provide rigorous and detailed designs of physically meaningful structures that directly feature such susceptibilities: a conceptual local meta-atom and a realistic nonlocal multilayered printed circuit board (PCB). Importantly, the latter leverages an intricate system of nearfield coupling beyond standard homogenization. Validated in simulations, our results provide a general and modular route to high-resolution asymmetric nonlocal metasurfaces for optical analog processing.
Water-repellent properties of superhydrophobic surfaces make them promising for anti-icing and deicing applications. Through experimental visualization of frozen sessile droplets undergoing melting on superhydrophobic surfaces, we identify a melting mode with the unmelted ice layer deposited at the bottom of the melting droplet, even though the density of ice is lower than that of water. In the deposited mode of the melting process, the time required for the frozen droplet to melt completely is much shorter than that in the floating mode. Force analysis shows that the melted fluid flows along the gas-liquid interface toward the top of the melting droplet, thereby exerting force and then suppressing the upward movement of the unmelted ice layer. Moreover, the flow within the liquid film formed between the unmelted ice layer and the heating wall is dominated by the viscous force, which has a lubrication effect and maintains the deposition of the unmelted ice layer. High heating temperature, large contact angle, and low particle concentration are helpful for the occurrence of the deposited mode.
We conduct numerical simulations of rocket plume impingement on a lunar landing surface using two complementary frameworks: the Boltzmann equation, which naturally captures rarefied gas dynamics, and the Navier Stokes (NS) equations, the conventional workhorse for turbulent flow simulations. We show that subtle rarefaction effects, long considered negligible in turbulent regimes, can become locally dominant within shear layers where viscous stresses predicted by the NS constitutive relation undergo sign reversals. This phenomenon, which we term constitutive degeneracy, produces order-one relative errors in predicted surface shear stress and heat flux. Our results demonstrate that turbulence can expose hidden limits of NS equations with broad implications for high-speed aerodynamics and planetary exploration.
We report the first empirical efficiency measurement of germanium-based thermophotovoltaic devices under high-temperature, high-irradiance conditions using a high view-factor calorimetric setup. Two TPV cell architectures were fabricated on p-type, highly doped (10^17 cm-3) Ge substrates, differing only in rear contact configuration. A standard device with a gold rear contact achieves a peak efficiency of 7.3 % and a power density of 1.77 W/cm2 at an emitter temperature of 1480 C, while a PERC-type device reaches 6.3 % efficiency and 1.22 W/cm2 at 1426 C. The superior performance of the standard device is attributed to lower series resistance, whereas the PERC design exhibits slightly higher efficiency at lower emitter temperatures (4.0 % vs. 3.8 % at 1150 C) due to enhanced rear-surface reflectivity. A detailed TPV model has been developed and validated across both device architectures. The model identifies out-of-band optical losses as the dominant efficiency-limiting mechanism, primarily caused by strong free-carrier absorption in the highly doped Ge substrate. Using this model, we predict device performance under idealized spectral conditions commonly assumed in prior literature. For a simulated AlN/W spectrally selective emitter, efficiencies as high as 22.3 % at 1800 C are obtained, consistent with previous semi-empirical predictions. In contrast, when previously reported Ge devices are modeled under the realistic graphite emitter spectrum used here, projected efficiencies decrease to as low as 8.1 % at 1480 C. These results show that earlier projections remain valid but idealized and underscore the importance of emitter spectral engineering and substrate optimization. Finally, we present the first direct comparison of Ge and InGaAs TPV devices under identical conditions, demonstrating the superior performance of InGaAs while confirming the cost-driven competitiveness of Ge.
Plasma sheaths play a central role in plasma-surface interactions, yet their thickness remains challenging to measure experimentally. Although classical analytical models such as the Child-Langmuir (CL) sheath model provide clear predictions for sheath thickness, experimental validation has been limited because most diagnostics either rely on indirect, multi-step inference (e.g., Langmuir probes) or require invasive and technically demanding techniques. In this work, we demonstrate that the plasma impedance probe (PIP), when operated with a controlled DC bias, enables relatively direct, model-informed measurements of sheath thickness that are reasonably straightforward to implement experimentally. Across a range of discharge conditions, biased-PIP sheath thickness measurements are found to follow CL scaling closely, requiring a single, consistent empirical correction factor of $α\approx 0.74$ to reconcile the measured thickness with CL predictions. Concurrent measurements of plasma density and electron damping show that probe biasing does not significantly perturb the bulk plasma density, supporting the validity of the biased-PIP approach. Building on this validation, we leverage the empirically determined $α$ factor to extend the floating (unbiased) PIP analysis to obtain model-dependent estimates of electron temperature and plasma potential without electrical biasing. A side-by-side comparison demonstrates close agreement between floating-PIP results and those obtained from a biased Langmuir probe. Taken together, these results establish the PIP as a complementary diagnostic to the Langmuir probe, expanding the range of accessible plasma measurements while providing experimental support for classical sheath models.
Understanding how cells escape from embedded spheroids requires a mechanical framework linking stress generation within cells, across cells, and between cells and the surrounding extracellular matrix (ECM). We develop such a framework by coupling a 3D vertex model of a spheroid to a fibrous ECM network and deriving a 3D Cauchy stress tensor for deformable polyhedral cells, enabling direct cell-level stress quantification in three dimensions. We analyze maximum shear stress in solid-like and fluid-like spheroids: solid-like spheroids exhibit broader stress distributions and radial stress gradients, while fluid-like spheroids show lower stresses with weak spatial organization. Cell shape anisotropy is not generically aligned with principal stress directions, indicating that morphology alone is an unreliable proxy for mechanical state. We further demonstrate strain stiffening at the single-cell level, where elongation produces nonlinear increases in maximum shear stress, allowing boundary cells in otherwise low-stress, fluid-like spheroids to transiently generate forces sufficient to remodel the matrix. To connect strain-induced stress amplification to invasion modes, we introduce an extended 3D vertex model with explicit, tunable cell-cell adhesion springs. In this minimal mechanical framework, single-cell breakout results from strain stiffening combined with reduced adhesion, whereas multi-cell streaming additionally requires anisotropic adhesion strengthened along the elongation axis and weakened orthogonally. Together, these results identify distinct mechanical pathways coupling cell strain, stress amplification, and adhesion organization to spheroid invasion.
Radial asymmetry in a quadrupole mass filter (QMF), introduced by symmetric displacement of a diagonally opposite rod pair, modifies the confinement potential and alters the ion stability characteristics. In this work, the influence of such radial asymmetry on QMF operation in the second stability zone is investigated through simulations. Using an existing potential formulation for a radially asymmetric QMF, the stability diagram in the second stability zone is extracted for the first time, revealing a systematic shift of the stability apex with the asymmetry parameter. Radial asymmetry introduces additional multipole components, notably an octupole term whose magnitude scales with the asymmetry parameter and whose sign depends on the DC polarity applied to the displaced rods. Transmission simulations show that the transmission peak shifts in accordance with the displaced stability apex, while the resolution exhibits a strong dependence on both the asymmetry parameter and the DC polarity. Comparable maximum resolution is obtained for inward and outward displacements under suitable polarity configurations, with outward displacement providing higher transmission efficiency due to an increased effective aperture. These results demonstrate that controlled radial asymmetry provides an effective means to tailor the resolution and transmission characteristics of QMFs operating in the second stability zone.
Pumped storage type hydropower plants play an important role in mitigating real-time energy flexibility. Reversible pump-turbines undergo extreme operating conditions such as runaway and speed-no-load. Very limited studies are undertaken to understand the stochastic flow under these conditions in the reversible pump-turbine. The present study investigates the unsteady vortical flow, its transportation, and interaction with the blades at speed-no-load. Large eddy simulations are conducted in both turbine and pump modes. The computational domain contains 120 million nodes. Numerical results provided evidence of a large longitudinal vortex that develops on the high-pressure side of the blade, and transports into the blade passage and develops the unsteady "string of swirls". The results also showed another "string of swirls" in the draft tube, where flow in the center is reversible (pumping). The resulting flow instability is very high, and it has the potential to induce fatigue damage to the blades.
Chiroptical signals, optical responses sensitive to molecular handedness, are rapidly suppressed by multiple scattering, fundamentally limiting their use in turbid media. Here we show that coupling molecular chirality to the topological structure of twisted light generates a protected phase observable that survives strong scattering. When Laguerre-Gaussian beams carrying orbital angular momentum propagate through chiral media, spin-orbit interaction converts circular birefringence into an azimuthal rotation of the helical wavefront. Remarkably, this chiral phase memory persists at scattering strengths that fully depolarize conventional beams, with the rotation magnitude preserved quantitatively between transparent solutions and strongly scattering tissue. The sign of the azimuthal rotation encodes molecular handedness: opposite enantiomers produce mirror-symmetric phase maps even after multiple scattering. Differential measurements between conjugate topological charges isolate the chiral contribution while cancelling achiral background, enabling the resolution of refractive-index changes of order 10-6. These results establish topological phase observables as robust carriers of weak chiral light-matter interactions in complex media, opening new routes for chiroptical spectroscopy and sensing beyond the ballistic-photon regime.
We investigate the third harmonic response of platinum ditelluride (PtTe$_2$), a Dirac semimetal belonging to the transition metal dichalcogenides class. Due to its topological properties, this material has drawn a lot of attention, particularly because it hosts type-II (super-critically tilted) Dirac fermions in the $\rm A-Γ-\rm A$ high symmetry direction. Adopting a low-energy model fitted directly from density functional theory band structure simulations, we calculate analytically the nonlinear conductivity. We observe that third-order optical nonlinearities are efficiently modulated by the cones tilting, which produces a significant enhancement of the nonlinear susceptibility. Our results, besides shedding light on topological transitions of platinum ditelluride, are relevant for future nanophotonic devices exploiting the tunable nonlinear properties of type-II Dirac fermions.
We report an integrated platform for efficient second-order nonlinear interactions based on linearly uncoupled resonators and angular phase matching. The proposed architecture confines phase control to a limited section of the device, maximizing field enhancement and effective nonlinear interaction length while simultaneously reducing the overall footprint. As an example we show the results for an AlGaAs-on-insulator structure demonstrating a photon-pair generation rate of 3.16 GHz/mW in the continuous-wave regime and 5.89 MHz under pulsed pumping. The generated biphoton state exhibits a Schmidt number K=1.02, indicating nearly uncorrelated photon pairs. The compact and reconfigurable nature of this approach, together with its independence from material-specific poling techniques, makes it applicable to a broad class of integrated $χ_2$ nonlinear platforms.
Understanding the thermal behaviour of superconducting accelerator magnets is essential to ensure stable and reliable operation. This work presents an extension of the Finite Element Quench Simulator (FiQuS) Multipole module to include collar and pole regions of accelerator magnets, which influences the overall thermal response. A thermal thin shell approximation (TSA), which is shown to be effective in previous works, is employed to model thermal insulation layers efficiently, replacing an insulation surface mesh. The main novelty of this work lies in the development of a method to model the thermal connection between the magnet winding and the collar and pole regions via the TSA. To assess the accuracy and computational efficiency of the novel method, temperature and field variations are computed for a current ramp-up scenario. The thermal solution is coupled to a fully resolved magnetodynamic solution to capture the interaction between thermal and electromagnetic behaviour. The results obtained with the TSA are then compared to classical finite element (FE) solutions with explicitly meshed insulation domains. The TSA predicts the maximum temperature within 2-4% of the reference solution while substantially reducing mesh complexity and achieving up to a 5 times speed-up in computation time. While the TSA has traditionally been employed for short-duration quench simulations with high heat fluxes between magnet turns, these results demonstrate its reliability and efficiency for current ramp scenarios with low heat fluxes, significantly expanding its application range beyond what has been previously reported in the literature. To illustrate potential applications of this new functionality, conduction cooling through the collar region is studied, comparing different cooling configurations and collar materials.
A fully coupled implicit finite-volume algorithm for incompressible viscoelastic interfacial flows is proposed, whereby the viscoelasticity of the flow is described by an upper-convected Maxwell constitutive model, including limited extensibility and shear-thinning behaviour. The governing equations describing the conservation of continuity and momentum, as well as the constitutive model are discretized using standard finite-volume methods and are solved for pressure, velocity and the polymer stress tensor in a single linear system of equations. Treating all terms of the linearized and discretized governing equations implicit in velocity, pressure and/or the components of the polymer stress tensor, a tightly coupled system of equations is obtained. The interface separating the interacting bulk phases and the surface tension acting at the fluid interface are modelled using a state-of-the-art front-tracking method. We demonstrate the capabilities of the proposed numerical framework with four representative test cases, including the deformation of a viscoelastic droplet in shear flow at large Weissenberg numbers of up to Wi=10^4, and the jump discontinuity of the rise velocity of a bubble rising in a viscoelastic liquid as a result of a "negative wake". Contrary to previous studies using segregated algorithms, the proposed fully coupled implicit algorithm does not apply or require a log-conformation approach to predict these flows. Overall, the fully implicit coupled front-tracking formulation provides a robust framework to reliable numerical predictions of strongly elastic interfacial flows at large Weissenberg numbers.
Multicore fibers (MCF) are perspective media for telecommunications, sensing, imaging and laser technologies. Here, the effect of beam stirring between weakly coupled cores is observed for sub-nanosecond transform-limited pulses of several kW peak power propagating in ~10 m long 7-core fiber, for the first time to our knowledge. In contrast to low-power domain where the output power distribution in the cores is random with large fluctuations sensitive to fiber disturbances, at high power of the input pulse injected in the central core the output power becomes equalized between the cores with fluctuations reduced to <5% being insensitive to disturbances. Similar behavior is observed in cut-back experiments showing that equi-partition is approached at a distance of ~5 m. The performed modeling describes well the experimental results and clarifies mechanisms of the new effect reasoned by a large nonlinear phase shift changing along the pulse and thereby resulting in statistical averaging over the pulse length of multidirectional power transfer processes between cores, thus leading to the robust equilibrium (equi-partitition for hexagonal MCF topology). At the same time, the combined output beam measured in a far field takes a stable bell-shaped profile instead of speckled beam at low powers, similar to the beam self-cleaning effect in multimode fibers.
Future large-scale magnetic confinement fusion reactors require operational regimes that can avoid extreme heat fluxes onto the plasma-facing components. One promising regime is the X-point radiator (XPR), which relies on a highly radiative, cold and dense plasma volume forming above the X-point, and which can be accessed via impurity seeding. Experimentally, the height of the XPR can be controlled by adjusting the seeding rate and heating power. This contribution presents axisymmetric (2D) simulations of the XPR regime in ASDEX Upgrade using the nonlinear MHD code JOREK extended with a kinetic particle framework for the main species neutrals and nitrogen impurities. With the time-dependent simulations, the progression from attached divertors to a complete detachment with the XPR formation is shown, highlighting the effects of the neutrals and impurities separately. Amidst this progression, the formation and the loss of the high-field-side high-density are observed. After the XPR is well-formed at the height of 6.8 cm, the fuelling and seeding rates are adjusted so that the XPR remains stationary. From the stationary case, the seeding rate is then changed to see how the XPR location reacts. By increasing and decreasing the seeding rate, the XPR responds by moving upwards and downwards, respectively. These simulations show JOREK's capability of simulating time-varying XPR, which will provide a baseline for the transition to 3D simulations, so the MHD activities and their interaction with the XPR can be studied.
Contrast variation small-angle neutron scattering (CV-SANS) has been widely employed for nano structural analysis of multicomponent systems. In CV-SANS experiments, scattering intensities of samples with different scattering co\ ntrasts are decomposed into partial scattering functions, corresponding to structure of each component and cross-correlation between different components, by singular value decomposition (SVD). However, the estimation of partial scattering functions with small absolute values often suffers from instability due to the significant differences in the singular values. In this paper, we propose a remedy for this instability by introducing the Tikhonov regularization, which ensures more stable reconstruction of the partial scattering functions.