2025-12-30 | | Total: 156
This study examines the applicability of two leading-edge dynamic stall criteria, namely, the maximum magnitudes of the leading-edge suction parameter (LESP) and the boundary enstrophy flux (BEF), in a moderately compressible flow regime. While previously shown to predict stall onset ahead of dynamic stall vortex (DSV) formation in incompressible and mildly compressible regimes, these criteria are assessed here at a Reynolds number of $1 \times 10^6$ and freestream Mach numbers between 0.3 and 0.5. Unsteady RANS simulations indicate that DSV formation occurs in close temporal proximity to the attainment of the stall criteria. However, at the highest Mach number considered, stronger shock interaction effects with the shear layer leads to DSV formation prior to the criteria being reached, reducing their predictive accuracy. These findings suggest that while the criteria remain effective at lower Mach numbers, their definitions require modification in compressible regimes where strong shock interactions significantly influence the stall process.
A standard procedure to achieve accurate, precise, and fast polarization measurement is to choose analyzing and generating polarization states that yield an $\ell^2$-condition number optimized instrument matrix. This strategy works well for rotating-waveplate systems, where the accessible polarization states trace a curve on the Poincaré sphere and the corresponding optimization problem is generally well posed. However, it becomes degenerate for liquid-crystal-based systems, which can generate arbitrary polarization states, and whose additional degrees of freedom allow the optimization of metrics beyond the $\ell^2$-condition number. Leveraging this unique advantage of liquid-crystal polarimeters, we introduce additional performance measures derived from alternative norms and error distributions computed via Monte Carlo simulations to inform the design of measurement schemes. We then experimentally demonstrate their effectiveness in suppressing errors, paving the way for more robust and efficient polarization measurements.
Achieving self-consistent performance predictions for ITER requires integrated modeling of core transport and divertor power exhaust under realistic impurity conditions. We present results from the first systematic power-flow and impurity-content study for the ITER 15 MA baseline scenario constrained directly by existing SOLPS-ITER neon-seeded divertor solutions. Using the OMFIT STEP workflow, stationary temperature and density profiles are predicted with TGYRO for $1.5 \le Z_{\rm eff} \le 2.5$, and the corresponding power crossing the separatrix $P_{\rm sep}$ is evaluated. We find that $P_{\rm sep}$ varies by more than a factor of 1.7 across this scan and matches the $\sim 100$~MW SOLPS-ITER prediction when $Z_{\rm eff} \simeq 1.6$ or when auxiliary heating is reduced to $\sim 75\%$ of nominal. Rotation-sensitivity studies show that plausible variations in toroidal flow magnitude modify $P_{\rm sep}$ by $\lesssim 20\%$, while AURORA modeling confirms that charge-exchange radiation inside the separatrix is dynamically negligible under predicted ITER neutral densities. These results identify a restricted compatibility window, $Z_{\rm eff} \approx 1.6$--1.75 and $0.75 \lesssim f_{P_{\rm aux}} \le 1.0$, in which core transport predictions remain aligned with neon-seeded divertor protection targets. This self-consistent, model-constrained framework provides actionable guidance for impurity control and auxiliary-heating scheduling in early ITER operation and supports future whole-device scenario optimization.
We present a unified theoretical framework for the transfer of angular momentum from a Bessel wave of twisted light to a fully absorbing disk of finite radius. Exact expressions for the orbital angular momentum density and the total angular momentum transmitted to the disk are obtained for both paraxial and nonparaxial regimes. By varying the beam wavelength, polarization, and cone angle, several experimentally relevant regimes of angular momentum transfer are identified. In the extreme nonparaxial regime, the absorbed angular momentum displays a staircase-like de- pendence on the object size, which can be interpreted as a geometric Hall-type response of the twisted field. The results suggest potential applications for controlled angular momentum transfer and size-sensitive probing of absorbing objects.
Additive manufacturing offers efficient production of plastic scintillators with nontrivial geometries using vat polymerization, allowing fabrication of geometries which would be difficult or even impossible to produce using conventional subtractive manufacturing. This work presents a novel photocurable scintillator formula that includes coumarin 450 as a tertiary dye to enable high-resolution 3D printing via the manipulation of the 405 nm cure light. Bulk photocured and 3D printed (with and without tertiary dye) samples were compared through observational assessment and spectral response. All samples showed pulse shape discrimination between neutron and gamma events. Inclusion of the tertiary dye has minimal impact on emission spectrum and light output, but significant impact on print resolution as shown by comparison of printed high-complexity geometries and feature resolution test objects. With the use of a cure-limiting dye, unsupported features, such as freestanding pillars, were resolvable down to 0.7 mm. Even finer resolution at or below 0.1 mm was achieved in fully supported, integrated structures printed with off-the-shelf 405 nm desktop 3D printer. Scintillators demonstrated a light output up to 50% of EJ-200 with a PSD figure of merit up to 1.35 at 0.9-1.1 MeVee.
In this paper, a semi-analytical model based on the Wiener-Hopf technique is developed to predict the turbulent boundary layer trailing edge noise from serrated edges, aiming to account for the correct three-dimensional noise source and propagation effects. The scattered surface pressure over a semi-infinite flat plate is first obtained using the Green's function developed for the acoustic scattering by a serrated edge. A radiation integral over the flat plate of a finite size is subsequently performed to obtain the far-field noise using Amiet's approach, capturing the correct three-dimensional source and propagation effects. The model is subsequently validated by comparing it against the two-dimensional Wiener-Hopf-based model under various serration sizes and frequencies. Far-field spectral predictions show close agreement between the three- and two-dimensional models at moderate observer distances around $r/c=1$, where $r$ and $c$ represent the observer distance and the chord length, respectively. However, unlike the two-dimensional model, the present model successfully captures the far-field $1/r$ decay in noise amplitudes. In addition, the predicted directivity agrees well with the two-dimensional model at most observer angles, but also captures the correct dipolar behaviour at upstream angles and additional high-frequency lobes due to interference patterns induced by the finite flat plate. Compared to the previous three-dimensional serrated models, the present model is based on the Wiener-Hopf technique and achieves a speed-up ratio of two orders of magnitude. It is hoped that such a model may be used to enable an efficient numerical optimisation of the serration shape in realistic applications.
MulAtoLEG (Multi-Atom Liouville Equation Generator) is an open-source Mathematica package for generating Liouville superoperators and Liouville equations, specialized for multilevel atomic systems comprising an arbitrary number of atoms. This scheme is based on an extension to multilevel atomic systems, originally developed by Lehmberg [R. H. Lehmberg, Phys. Rev. A 2, 883 (1970)] as an adjoint master equation for ensembles of two-level emitters and later reformulated by Genes [M. Reitz, C. Sommer and C. Genes, PRX Quantum 3, 010201 (2022)] as a master equation. The package facilitates the generation of equations for complex transition configurations in alkali atoms. Although primarily designed for atomic systems, it can also generate the master and adjoint master equations for general Hamiltonians and Lindbladians. In addition, it includes functionalities to construct the differential equations in the dressed-state basis, where, in many cases, the non-unitary evolution operator can be determined explicitly. To maximize computational efficiency, the package leverages Mathematica's vectorization and sparse linear algebra capabilities. Since MulAtoLEG produces exact equations without approximations, the feasible system size is naturally limited by the available computational resources.
We introduce a new approach to deriving approximate analytical solutions of a harmonic oscillator damped by purely nonlinear, or combinations of linear and nonlinear damping forces. Our approach is based on choosing a suitable trial solution, i.e. an ansatz, which is the product of the time-dependent amplitude and the oscillatory (trigonometric) function that has the same frequency but different initial phase, compared to the undamped case. We derive the equation for the amplitude decay using the connection of the energy dissipation rate with the power of the total damping force and the approximation that the amplitude changes slowly over time compared to the oscillating part of the ansatz. By matching our ansatz to the initial conditions, we obtain the equations for the corresponding initial amplitude and initial phase. Here we demonstrate the validity of our approach in the case of damping quadratic in velocity, Coulomb damping, and a combination of the two, i.e. in this paper we consider purely nonlinear damping, while the dynamics with combinations of damping linear in velocity and nonlinear damping will be analyzed in a follow-up paper. In the case of damping quadratic in velocity, by comparing our approximate analytical solutions with the corresponding numerical solutions, we find that our solutions excellently describe the dynamics of the oscillator in the regime of weak to moderately strong quadratic damping. In the case of Coulomb damping, as well as in the case of a combined Coulomb and quadratic damping, our approximate analytical solutions agree well with the corresponding numerical solutions until the last few half-periods of the motion. Therefore, for these two cases, we introduce improved variants of our approximate solutions which describe the dynamics well until the very end.
Photonic topological insulators (PTIs) offer robust platforms for light manipulation, but reconfigurable control of their topological properties without degrading performance remains a major challenge. While phase-change materials (PCMs) provide large refractive index modulation, widely used materials such as Ge2Sb2Te5 (GST) have been successfully deployed in commercial applications including optical data storage. However, they exhibit significant optical absorption in their crystalline state, which poses a challenge for transmissive photonic devices such as PTIs where high transparency is essential. Here, we overcome this fundamental limitation by integrating the ultra-low-loss PCM antimony triselenide (Sb2Se3) onto a silicon-based 2D PTI. We achieve submicron-scale selective patterning of Sb2Se3 on a photonic crystal for the first time, and demonstrate a topological phase transition induced by the material phase change. Owing to the transparency of Sb2Se3 in both its amorphous and crystalline states, a high Q-factor on the order of 10^3 is preserved-representing nearly an order-of-magnitude improvement over previous GST-based devices. This work resolves the absorption-loss bottleneck in reconfigurable PTIs and paves the way for practical, low-loss, tunable topological photonic devices.
The design of chiral metasurfaces with tailored optical properties remains a central challenge in nanophotonics due to the highly nonlinear relationship between geometry and chiroptical response. Machine-learning-assisted optimization pipelines have recently emerged as efficient tools to accelerate this process, yet their performance strongly depends on the choice of neural-network (NN) architecture. In this work, we integrate the NeuroEvolution of Augmenting Topologies (NEAT) algorithm into an established deep-learning optimization framework for dielectric chiral metasurfaces. NEAT autonomously evolves both network topology and connection weights, enabling task-specific architectures without manual tuning, whereas the reinforcement-learning strategy in our framework evolves knowledge of the solution space and fine-tunes a model's weights in parallel. Using a pipeline-produced dataset of 9,600 simulated GaP metasurface geometries, we evaluate NEAT under varying input dimensionalities, feature-scaling methods, and data sizes. With standardized feature scaling yielding the most consistent performance for both examined output dimensionalities, the relatively compact NEAT-evolved NN models, when integrated into the full optimization pipeline, achieve similar or improved predictive accuracy and generalization compared to initially employed dense few-layer perceptrons. Accordingly, these resource-efficient models successfully perform inference of metasurfaces exhibiting strong circular dichroism in the visible spectrum, allowing for transfer learning between simulated and experimental data. This approach demonstrates a scalable path toward adaptive, self-configuring machine-learning frameworks for automated photonic design both standalone and as building block for agentic artificial intelligence (AI).
Early cancer detection relies on invasive tissue biopsies or liquid biopsies limited by biomarker dilution. In contrast, tumour-derived extracellular vesicles (EVs) carrying biomarkers like melanoma-associated antigen-A (MAGE-A) are highly concentrated in the peri-tumoral interstitial space, offering a promising near-field target. However, at micrometre scales, EV transport is governed by stochastic diffusion in a low copy number regime, increasing the risk of false negatives. We theoretically assess the feasibility of a smart-needle sensor detecting MAGE-A-positive microvesicles near a tumour. We use a hybrid framework combining particle-based Brownian dynamics (Smoldyn) to quantify stochastic arrival and false negative probabilities, and a reaction-diffusion PDE for mean concentration profiles. Formulating detection as a threshold-based binary hypothesis test, we find a maximum feasible detection radius of approximately 275 micrometers for a 6000 s sensing window. These results outline the physical limits of proximal EV-based detection and inform the design of minimally invasive peri-tumoral sensors.
Structured optical beams possess rich spatial features that are commonly characterized using entropic measures of field complexity. However, such measures do not directly quantify the operational usefulness of optical structure for parameter estimation and sensing. Here we introduce Fisher information as an operational metric to assess the metrological content of structured optical fields. By treating the measured intensity distribution as a statistical object, we define Fisher information with respect to physically relevant parameters, such as transverse displacement. We demonstrate that optical modes with comparable Shannon entropy can exhibit markedly different Fisher information, revealing sensitivity features associated with nodal structure and local curvature. Using Hermite--Gaussian modes as minimal test cases, we show that increasing modal order systematically enhances Fisher information. We then extend the analysis to two widely used families in structured light: Laguerre--Gaussian vortex beams and finite-energy Bessel--Gauss beams. Across these representative families, Fisher information provides a unified and experimentally accessible criterion for comparing structured optical fields in sensing applications.
In many biological networks the responses of individual elements are ambiguous. We consider a scenario in which many sensors respond to a shared signal, each with limited information capacity, and ask that the outputs together convey as much information as possible about an underlying relevant variable. In a low noise limit where we can make analytic progress, we show that individually ambiguous responses optimize overall information transmission.
The remarkable advancements in our capacity to synthesise structured light have facilitated the generation of any desired optical landscapes and even controlling the spatial distribution of light propagating through optically complex media such as multimode fibres. The availability of precisely defined structured light at the extremity of an exceedingly narrow and flexible cable holds the potential to stimulate a diverse range of highly sought-after applications, encompassing rapid communication, quantum computing, and, notably, imaging. What we lack in reaching these aspirations is the resilience of such light transport to deformations of the waveguide. Although recent theoretical investigations have delineated the attributes of ideal multimode fibres capable of deformation-enduring conveyance of structured light, tangible fibres possessing this indispensable trait to a practical extent remain elusive. Our study takes a deep dive into the precision of commercially available multimode fibres with the highest probability of demonstrating the phenomenon. We identified minuscule imperfections in their refractive index distribution, examined how these affect light transport when the fibre is deformed, and studied their implications for imaging applications. Our investigation has confirmed that these imperfections are indeed responsible for the undesirable alterations introduced into the output structured light fields during bending. Finally, as an alternative to standard graded-index fibres, manufactured by drawing silica-based preforms, we present narrow multimode waveguides in which the refractive-index profile has been established by ion exchange. These waveguides indeed exhibit previously unseen resilience of structured light transport even under severe deformation conditions and aptly fulfil the requirements of imaging applications.
Bone marrow (BM) metabolic quantification with 18F-fluorodeoxyglucose (FDG) positron emission tomography (PET) is of broad clinical significance for accurate assessment of BM at staging and follow-up, especially when immunotherapy is involved. However, current methods of quantifying BM may be inaccurate because the volume defined to measure bone marrow may also consist of a fraction of trabecular bone in which 18F-FDG activity is negligible, resulting in a potential underestimation of true BM uptake. In this study, we demonstrate this bone-led tissue composition effect and propose a bone fraction correction (BFC) method using X-ray dual-energy computed tomography (DECT) material decomposition. This study included ten scans from five cancer patients who underwent baseline and follow-up dynamic 18F-FDG PET and DECT scans using the uEXPLORER total-body PET/CT system. The voxel-wise bone volume fraction was estimated from DECT and then incorporated into the PET measurement formulas for BFC. The standardized uptake value (SUV), 18F-FDG delivery rate K1, and net influx rate Ki values in BM regions were estimated with and without BFC and compared using the statistical analysis. The results first demonstrated the feasibility of performing voxel-wise material decomposition using DECT for metabolic BM imaging. With BFC, the SUV, K1, and Ki values significantly increased by an average of 13.28% in BM regions compared to those without BFC (all P<0.0001), indicating the impact of BFC for BM quantification. Parametric imaging with BFC further confirmed regional analysis. Our study using DECT suggests current SUV and kinetic quantification of BM are likely underestimated in PET due to the presence of a significant bone volume fraction. Incorporating tissue composition information through BFC may improve BM metabolic quantification.
Adiabatic topological pumping promises robust transport of energy and information, but its speed is fundamentally limited by the instantaneous adiabatic condition, which demands prohibitively slow parameter variations. Here we develop a paradigm shift from instantaneous to global adiabaticity. We derive a global adiabatic criterion (GAC), which sets an absolute fidelity bound by controlling the root-mean-square value of nonadiabaticity factor. We further introduce a fluctuation-suppression acceleration criterion, which minimizes spatial inhomogeneity and allows us to safely increase the mean nonadiabaticity. Experimentally, we implement this principle in femtosecond-laser-written photonic Su-Schrieffer-Heeger waveguide arrays via scalable power-law coupling modulation. Our accelerated topological pumping achieves >0.95 fidelity over a fivefold reduced device length compared to the conventional scheme, exhibits the predicted linear scaling with the system size, and maintains robust performance across a >400 nm bandwidth. This principle of GAC provides a universal design rule for fast, compact, and robust adiabatic devices across quantum and classical topological platforms.
Local high strain in solid rocket motor grains is a primary cause of structural failure. However, traditional numerical simulations are computationally expensive, and existing surrogate models cannot explicitly establish geometric models and accurately capture high-strain regions. Therefore, this paper proposes an adaptive graph network, GrainGNet, which employs an adaptive pooling dynamic node selection mechanism to effectively preserve the key mechanical features of structurally critical regions, while concurrently utilising feature fusion to transmit deep features and enhance the model's representational capacity. In the joint prediction task involving four sequential conditions--curing and cooling, storage, overloading, and ignition--GrainGNet reduces the mean squared error by 62.8% compared to the baseline graph U-Net model, with only a 5.2% increase in parameter count and an approximately sevenfold improvement in training efficiency. Furthermore, in the high-strain regions of debonding seams, the prediction error is further reduced by 33% compared to the second-best method, offering a computationally efficient and high-fidelity approach to evaluate motor structural safety.
The efficiency of recombination is of crucial importance for the existence of ultracold plasmas, particularly, the ones formed in the magneto-optical traps. Unfortunately, a straightforward simulation of the recombination encounters the problem of huge difference in the spatial and temporal scales for free and bound motion of the electrons. As a result, only the "virtual" electron-ion pairs are usually reproduced in such simulations, and it is necessary to employ some additional criteria to identify them with the recombined atoms (this might be a minimal number of revolutions of the electron about the nearest ion or a maximal distance between them). It is the aim of this paper to present the first successful ab initio simulation of the recombination without any auxiliary assumptions. We employed a special algorithm, which was based on: (i) using the "scalable" reference frame, co-moving with the expanding plasma, (ii) dynamical choice of the number of "mirror" cells, taking into account in calculation of the Coulomb sums, and (iii) accurate treatment of the singular interparticle interactions, without any truncation or "softening" of the Coulomb forces. Then, the recombination events are identified by a series of sharp equidistant peaks in the kinetic and/or potential energies for a sample of particles, which are caused by the captured electrons passing near the pericenters of their orbits; and this is confirmed by a detailed inspection of the particle trajectories. Thereby, we were able to trace formation of the real - rather than "virtual" - electron-ion pairs. The total efficiency of recombination for the realistic experimental conditions was found to be about 20%, which is in perfect agreement both with the laboratory measurements and with the earlier semi-empirical simulations.
This letter presents the design and implementation of a compact high-efficiency octave microwave rectifier. A key highlight is the novel segmented impedance matching method, a unique approach that expands the rectifier bandwidth. The diode reactance is initially regulated by a series short-ended microstrip line. Impedance-compensated structures, characterized by varying admittance properties across an extensive frequency range, partition the operating frequency band into two segments based on the input impedance, thereby minimizing impedance variation. Ultimately, the input impedance is matched by a novel triple-band matching network. An octave rectifier was fabricated and measured. Results demonstrate that the rectifier achieves over 50% efficiency over 1.3-2.55 GHz fractional bandwidth 64.9% at 0-dBm RF input power. Even with a decrease in input power to -10 dBm, the rectifier maintains over 30% efficiency.
Physics-Informed Neural Networks (PINNs) are a methodology that aims to solve physical systems by directly embedding PDE constraints into the neural network training process. In electromagnetism, where well-established methodologies such as FDTD and FEM already exist, new methodologies are expected to provide clear advantages to be accepted. Despite their mesh-free nature and applicability to inverse problems, PINNs can exhibit deficiencies in terms of accuracy and energy metrics when compared to FDTD solutions. This study demonstrates hybrid training strategies can bring PINNs closer to FDTD-level accuracy and energy consistency. This study presents a hybrid methodology addressing common challenges in wave propagation scenarios. The causality collapse problem in time-dependent PINN training is addressed via time marching and causality-aware weighting. In order to mitigate the discontinuities that are introduced by time marching, a two-stage interface continuity loss is applied. In order to suppress loss accumulation, which is manifested as cumulative energy drift in electromagnetic waves, a local Poynting-based regularizer has been developed. In the developed PINN model, high field accuracy is achieved with an average 0.09\% $NRMSE$ and 1.01\% $L^2$ error over time. Energy conservation is achieved on the PINN side with only a 0.024\% relative energy mismatch in the 2D PEC cavity scenario. Training is performed without labeled field data, using only physics-based residual losses; FDTD is used solely for post-training evaluation. The results demonstrate that PINNs can achieve competitive results with FDTD in canonical electromagnetic examples and are a viable alternative.
We report the experimental realization of a LiNbO3 metasurface for electro-optic modulation of light polarization in the telecommunication band. High-Q quasi-bound states in the continuum are emploied to enhance the modulation of amplitude and phase of an impinging beam by a driving electric field, leading to efficient polarization rotation and conversion. We quantified modulation effects under a CMOS-compatible bias at 1 MHz frequency, achieving a variation of 5% in the Stokes parameters and a variation of the polarization ellipse angles of about 3° for the transmitted light. These results demonstrate that dynamic polarization and phase modulation can be attained in a compact platform, highlighting the potential of high-Q resonant LiNbO3 metasurfaces for enhanced light-matter interaction in subwavelength electro-optic devices.
In our previous study [Medipati \textit{et al}., (2025) \textit{J. Fluid Mech}. \textbf{1014}, A34] \cite{medipati2025elliptic}, a detailed experimental investigation is performed on the elliptical liquid jets in a supersonic cross-flow ($M_{\infty}$ = 2.5), focusing on the effect of orifice aspect ratio ($AR$ = spanwise dimension/streamwise dimension) on the atomization mechanism for a fixed momentum flux ratio ($J$). In this paper, we present experimental studies that show the influence of $J$ on the jet breakup mechanism, shock structures, and unsteady interactions for each $AR$. A wide range of $J$ values (1.5 to 9.7) and three $AR$ cases (0.3, 1, and 3.3) are chosen for the study. We find that in the case of lower $J$, the jet exhibits large unsteadiness, with larger wavelength Rayleigh-Taylor (RT) waves on the windward surface. In contrast, as the $J$ increases, the unsteadiness decreases, smaller and more regular RT wavelength is formed due to the enhanced drag resulting from the reduced jet deflection. However, irrespective of $J$, in the case of $AR$ = 0.3 and 1, the primary atomization mechanism is due to the formation of Kelvin-Helmholtz instabilities (KHI) on the lateral surfaces. Furthermore, in the case of lower $J$, the shock waves formed upstream of the jet are highly corrugated with significant variations in time. The intense interaction of the liquid jet with the oncoming boundary layer streaks, in the case of lower $J$, is the primary source of large-scale unsteadiness. These findings highlight the significance of $J$ on the atomization mechanism in supersonic cross-flow.
We present a space-time extension of a conservative Cartesian cut-cell finite-volume method for two-phase diffusion in prescribed-motion geometries. The formulation follows a two-fluid approach: one scalar field is solved in each phase with discontinuous material properties, coupled by sharp interface conditions enforcing flux continuity and jump laws. To handle moving boundaries on a fixed Cartesian grid, the discrete balance is written over phase-restricted space-time control volumes, whose geometric moments (swept volumes and apertures) are used as weights in the finite-volume operators. This construction naturally accounts for the creation and destruction of cut cells (fresh/dead-cell events) and yields strict discrete conservation. The resulting scheme retains the algebraic structure of the static cut-cell formulation while incorporating motion through local geometric weights and interface coupling operators. A series of verification and validation tests in two and three dimensions demonstrate super-linear accuracy in space, robust behavior under repeated topology changes and conservation across strong coefficient jumps and moving interfaces. The proposed space-time cut-cell framework provides a conservative building block for multiphase transport in evolving geometries and a foundation for future free-boundary extensions such as Stefan-type phase change.
The static and dynamic electric dipole polarizabilities of the $5s^2~^1\!S_0$ and $5s5p~^3\!P_{0,2}$ states of Sr atoms are calculated using the relativistic configuration interaction plus the many-body perturbation theory (RCI+MBPT) method. Magic wavelengths are determined for the transitions $5s^2~^1\!S_0\rightarrow 5s5p~^3\!P_{0}$, $5s^2~^1\!S_0\rightarrow 5s5p~^3\!P_{2}$, and $5s5p~^3\!P_0\rightarrow 5s5p~^3\!P_{2}$. A comprehensive study is conducted on the dependence of magic wavelengths on the angle between the laser polarization and the magnetic field. Furthermore, the conditions for realizing triple magic trapping at 813.4~nm for the $5s^2~^1\!S_0$, $5s5p~^3\!P_{0}$ and $5s5p~^3\!P_{2}$ states are investigated. In the case of linearly polarized light, when the angle ($θ_p$) between the laser polarization direction and the magnetic field is $79.1(0.7)^\circ$, triple magic trapping for the $5s^2~^1\!S_0$, $5s5p~^3\!P_{0}$, and $5s5p~^3\!P_{2}~M=0$ states can be achieved. This result agrees well with the recent experimental measurement (78.49(3)$^\circ$)[Phys. Rev. Lett. 135, 143401 (2025)]. Meanwhile, triple magic trapping involving the $5s5p~^3\!P_{2}~M=2$ state can be achieved when $θ_p= 37.4(0.3)^\circ$. The conditions for achieving triple magic trapping with circularly and arbitrarily elliptically polarized light are also presented.
Multicomponent vesicles suspended in viscoelastic fluids are crucial for understanding a variety of physiological processes. In this work, we develop a continuum surface force (CSF) phase-field model to investigate the hydrodynamics of inextensible multicomponent vesicles in viscoelastic fluid flows with inertial forces. Our model couples a fluid field comprising both Newtonian and Oldroyd-B fluids, a surface concentration field representing the multicomponent distribution on the vesicle membrane, and a phase-field variable governing the membrane evolution. The viscoelasticity effect of extra stress is well incorporated into the full Navier-Stokes equations in the fluid field. The surface concentration field is determined by Cahn-Hilliard equations, while the membrane evolution is governed by a nonlinear advection-diffusion equation. The membrane is coupled to the surrounding fluid through the continuum surface force (CSF) framework. To ensure stable numerical solutions of the highly nonlinear multi-field model, we employ a residual-based variational multiscale (RBVMS) method for the Navier-Stokes equations, a Streamline-Upwind Petrov-Galerkin (SUPG) method for the Oldroyd-B equations, and a standard Galerkin finite element framework for the remaining equations. The system of PDEs is solved using an implicit, monolithic scheme based on the generalized-$α$ time integration method. To enhance spatial accuracy, we employ isogeometric analysis (IGA). We present a series of two-dimensional numerical examples in shear and Poiseuille flows to elucidate the influence of membrane composition and fluid viscoelasticity on the hydrodynamics of multicomponent vesicles.