2025-12-08 | | Total: 103
Magnetic shear across the polarity inversion line (PIL) plays an important role in the explosive nature of reconnection onset and in the equilibration of current sheets, acting as a source of free energy that can enhance or inhibit the onset process under certain conditions. In this study, we use a 2D PIC simulation to examine the local interaction between the reconnection guide field and thermodynamic variables during reconnection onset in a region of initially depleted thermal energy and enhanced magnetic energy in a large guide field background. We identify critical stages of the equilibration process, characterize intervals based on whether the pressure evolution is driven by changes in density or temperature, and discuss what these intervals imply about the evolution of local heat and work density. Finally, we examine power densities associated with electromagnetic field time evolution and electromagnetic energy transfer and compare to those related to thermodynamic changes.
We consider a pair of driver/witness electron bunches propagating in an ionized gas background a configuration similar to the one produced in a capillary discharge where a plasma oscillation has been excited by a driving pulse. We assume as in the plasma nonlinear regime that the plasma electrons behind the driver are completely expelled and an ellipsoidal cavity filled with ions only is formed. The fields are linear in both longitudinal and transverse directions, at least in the region of interest for particle acceleration, as the one produced by a uniform ion distribution within a uniformly charged ellipsoidal volume. The fields produced by the ions and experienced by a witness electron beam are purely electrostatic, being the ions at rest in the laboratory frame on the time scale of interest and it can be represented with the field distribution produced by a 3D charged ellipsoidal. The energy spread and emittance degradation has been studied by slicing the bunch in an array of cylinders and solving envelope equations for each bunch slice. The properties of transverse envelope and emittance oscillations and energy spread degradation have been analyzed together with the related matching conditions for optimal transport and acceleration.
There is an ongoing interest in using spectroscopy in inertial confinement fusion (ICF) experiments, where dopants such as krypton can provide vital information about the temperature and density of the imploding plasma. While the most advanced tools for calculating Stark profiles are computer simulation models (CSMs), their application to complex lineshapes under the extreme conditions of ICF experiments is computationally challenging. In this manuscript, we present results of several CSM realizations applied to the Stark shape of the krypton He-beta line and its satellites at ICF-relevant conditions (ne = 1e24 to 1e25 cm-3, Te = 3keV). We demonstrate that codes with the same underlying physics but different numerical approaches yield identical results and analyze the differences in the line profile caused by various physical effects.
The two-layer quasigeostrophic model (2LQG) and the Eady model are two idealized systems illustrating the baroclinic instability of atmospheric jets and ocean currents. The two setups share many ingredients -- background vertically sheared zonal flow of density-stratified fluid in a rapidly rotating frame -- while differing in complexity and dimensionality. The Eady model has a continuous vertical direction, with baroclinic turbulence induced by boundary potential vorticity (PV) gradients at top and bottom. By contrast, the 2LQG sytem typically models baroclinic instability induced by interior PV gradients. This distinction challenges our ability to clearly identify a couple of 'modes' through which the Eady dynamics could be inferred from a simpler 2LQG system. In the present study, we show that this difficulty can be circumvented in the turbulent regime arising for weak bottom drag. Namely, guided by the common organization of both systems into a gas of coherent vortices, we identify a quantitative mapping between the Eady and the 2LQG models. The mapping allows for parameter-free predictions of the eddy diffusivity of the Eady model based on the knowledge of the 2LQG diffusivity. We illustrate these results using numerical simulations of the Eady and 2LQG models with linear or quadratic bottom drag.
We introduce an alternative receiver architecture for deep-space optical communication, in which a single large aperture is replaced by an array of smaller ones with outputs combined coherently, employing phase stabilization based on photon counting events. We show that it allows to increase the signal to noise ratio, thus potentially attaining higher information transmission rates in the regime of large noise, typical for daytime communication. We analyze its practical performance by simulating pulse position modulation-based communication from the recently launched Psyche mission. Under nighttime conditions the achieved performance is comparable to that offered by a single large aperture, whereas in daytime conditions the single photon coherent beam combination architecture provides an advantage in the information transmission rate.
Functional near-infrared spectroscopy (fNIRS) is impacted by signal contamination from superficial hemodynamics. It is important to develop methods that account for such contamination and provide accurate measurements of cerebral hemodynamics. This work aims to investigate whether simulated data with two-layer or three-layer tissue models are able to reproduce in vivo data collected with dual-slope (DS) frequency-domain (FD) near-infrared spectroscopy (NIRS) on human subjects during brain activation. We performed Monte Carlo simulations to generate DS FD-NIRS data from two- and three-layer media with a range of layer thicknesses and optical properties. We collected in vivo data with DS FD-NIRS (source-detector distances: 25, 37 mm; wavelengths: 690, 830 nm; modulation frequency: 140 MHz) over the occipital lobe of human subjects during visual stimulation. Simulated and in vivo data were analyzed with diffusion theory for a homogeneous medium and results were compared for each DS FD-NIRS data type. We found that the main qualitative features of in vivo data could be reproduced by simulated data from a three-layer medium, with a second layer (representing the cerebrospinal fluid in the subarachnoid space) that is less absorbing and less scattering than the other two layers, and with a top layer thickness that represents the combined scalp and skull thickness. A three-layer model is a viable improvement over a homogeneous model to analyze DS FD-NIRS data (or any other fNIRS data) to generate more accurate measurements of cerebral hemodynamic changes without a need for large data sets for tomographic reconstructions.
We report the coherent excitation of the highly forbidden $^2S_{1/2} \rightarrow {^2}F_{7/2}$ clock transition in the odd isotope $^{173}\mathrm{Yb}^+$ with nuclear spin $I = 5/2$, and reveal the hyperfine-state-dependent, nuclear spin induced quenching of this transition. The inferred lifetime of the $F_e = 4$ hyperfine state is one order of magnitude shorter than the unperturbed ${^2}F_{7/2}$ clock state of $^{171}\mathrm{Yb}^+$. This reduced lifetime lowers the required optical power for coherent excitation of the clock transition, thereby reducing the AC Stark shift caused by the clock laser. Using a 3-ion Coulomb crystal, we experimentally demonstrate an approximately 20-fold suppression of the AC Stark shift, a critical improvement for the scalability of future multi-ion $\mathrm{Yb}^+$ clocks. Furthermore, we report the $|^2S_{1/2},F_g=3\rangle~\rightarrow~|^2F_{7/2},F_e=6\rangle$ unquenched reference transition frequency as $642.11917656354(43)$ THz, along with the measured hyperfine splitting and calculated quadratic Zeeman sensitivities of the ${^2}F_{7/2}$ clock state. Our results pave the way toward multi-ion optical clocks and quantum computers based on $^{173}\mathrm{Yb}^+$.
We study the dynamics of inertial particles in turbulence using datasets obtained from both direct numerical simulations and laboratory experiments of turbulent swirling flows. By analyzing time series of particle velocity increments at different scales, we show that their evolution is consistent with a Markov process across the inertial range. This Markovian character enables a coarse-grained description of particle dynamics through a Fokker-Planck equation, from which we can extract drift and diffusion coefficients directly from the data. The inferred coefficients reveal scale-dependent relaxation and noise amplitudes, indicative of inertial filtering and intermittency effects. Beyond the kinematic description, we analyze the thermodynamic properties of particle trajectories by computing the trajectory-dependent entropy production. We show that the statistics of entropy fluctuations satisfy both the Integral Fluctuation Theorem and, under certain conditions, the Detailed Fluctuation Theorem. These results establish a quantitative bridge between stochastic thermodynamics and particle-laden flows, and open the door to modeling turbulent transport using effective stochastic theories constrained by data and physical consistency.
We present a continuous nonlinear optimization model for the Spin Glass Problem (SGP), building on a classical result by Rosenberg (1972), which shows that for a class of multilinear polynomial problems the optimal values of the continuous relaxation and the corresponding discrete model coincide. Using the SGP as a case study, we provide a simple, problem-specific argument showing how any optimal solution returned by a continuous solver can be converted into an optimal discrete spin configuration, even when the solver outputs non-integer values. The relaxed model remains nonconvex and does not alter the inherent computational hardness of the problem, but it offers a direct and conceptually transparent continuous formulation that can be handled by modern global optimization software. Computational experiments on standard benchmark instances indicate that this approach can match, and in several cases surpass, recent integer programming linearization techniques, making it a practical and complementary tool for researchers working at the interface between statistical physics and combinatorial optimization.
In the scope of solar axion searches, detectors which are able to measure low energetic X-rays with high efficiency are required. For this purpose a detector based on the GridPix technology was built for the CAST experiment at CERN. The GridPix is a pixelised readout ASIC (Timepix) with a Micromegas-like gas amplification stage (grid) built photolithographically on top. In order to reduce the detector's background level, several hardware and software vetoes were implemented. Hardware-wise, these vetoes consist of a ring of six GridPixes around a central GridPix, a readout of the induced grid signal, and two scintillators. On the software side, multiple approaches to distinguish between background events and X-ray photons are implemented. Here, also the hardware features, like the six surrounding GridPixes, are used. The new detector was tested in a long ($3500\,\text{h}$) background data taking campaign. The performance of the new vetoes was evaluated. The detector performance itself, for low energetic X-rays, was also evaluated with a variable X-ray generator using eight different energies from 0 to $10\,\text{keV}$. The efficiency for very low energetic X-rays and the energy resolution was determined.
This letter proposes a high-power and high-efficiency GaN-HEMT voltage-controlled oscillator. The VCO consists of a coupled-line coupler, an inverse class-F amplifier, and a novel frequency-tunable stepped-impedance resonator. Using a harmonic control circuit and a parasitic parameter compensation circuit, the power amplifier operates in the inverse class-F state to achieve high efficiency. The feedback circuit uses a coupled-line coupler instead of the traditional coupling capacitor to control feedback power precisely. The measurement results show that the VCO with an oscillation frequency of 2.41-2.45 GHz achieves a maximum conversion efficiency of 74.5% at 2.44 GHz and an output power of 40.2 dBm. It is a candidate for the microwave source in a wireless power transmission system.
Understanding the Earth s internal structure remains a major challenge, as traditional geophysical methods face ambiguities in linking seismic observations to temperature, composition, or mass density variations. Atmospheric neutrinos offer a complementary probe: while traversing the Earth, they undergo flavor oscillations that depend on the local electron density, which reflects both mass density and composition.
Two-dimensional (2D) materials have been proposed, among many other applications, as a efficient tool for the separation of atomic and molecular species and their corresponding isotopes, given the confinement provided by their subnanometric dimensions. In this work we present three dimensional quantum wave packet calculations revealing an enhancement in the quantum transport in bilayer over monolayer graphdiyne membranes, one of the most popular 2D materials which is commonly employed for this purpose. Besides, resonances emerge superimposed over the typical monolayer profile for transmission probabilities, a feature that is general to other bilayer nanoporous 2D heterostructures and that shows a strong dependence on the interlayer separation.
We present exact analogies between the tautochrone problem of mechanics and the squeezed states of quantum optics, to optical lattices. Both phenomena emerge in the same physical system, that of waveguide arrays with non-uniform couplings. Extension to two dimensions yields Lissajous-type trajectories and multidirectional tautochrone focusing. Furthermore, we investigate the impact of Kerr nonlinearity and show that it determines the diffraction behavior, namely coherent-state-like or squeezed propagation. These quantum inspired classical lattices highlight the role of the coupling coefficients to beam engineering and light control in complex media.
Purpose: The Unadjusted Langevin Algorithm (ULA) in combination with diffusion models can generate high quality MRI reconstructions with uncertainty estimation from highly undersampled k-space data. However, sampling methods such as diffusion posterior sampling or likelihood annealing suffer from long reconstruction times and the need for parameter tuning. The purpose of this work is to develop a robust sampling algorithm with fast convergence. Theory and Methods: In the reverse diffusion process used for sampling the posterior, the exact likelihood is multiplied with the diffused prior at all noise scales. To overcome the issue of slow convergence, preconditioning is used. The method is trained on fastMRI data and tested on retrospectively undersampled brain data of a healthy volunteer. Results: For posterior sampling in Cartesian and non-Cartesian accelerated MRI the new approach outperforms annealed sampling in terms of reconstruction speed and sample quality. Conclusion: The proposed exact likelihood with preconditioning enables rapid and reliable posterior sampling across various MRI reconstruction tasks without the need for parameter tuning.
Owing to the spatially local nature of the constitutive equations typically used to model polymeric stresses, the differential operators governing the linearized dynamics of bounded viscoelastic shearing flows have singular points. As a result, the eigenspectra of such shearing flows contain, in addition to discrete eigenvalues, continuous spectra (CS) comprising singular eigenfunctions. A clear understanding of the theoretical CS loci is crucial in discriminating physically genuine (discrete) eigenvalues from the poorly approximated numerical CS. For rectilinear shear flows of Oldroyd-B fluids, the CS are a pair of line segments, with lengths equal to the base-state range of velocities. In this study, we provide the first comprehensive account of the nature of the CS for both rectilinear and curvilinear shearing flows of the FENE-P fluid. In stark contrast to the CS for the Oldroyd-B fluid mentioned above, we show analytically that there are up to six distinct continuous spectra for shearing flows of FENE-P fluids. When the finite extensibility parameter $L > 50$, as appropriate for large molecular weight polymers used in experiments, three of the CS are nearly identical, and independent of the solvent-to-solution viscosity ratio ($β$). The other three CS are $β$-dependent, with one of them being the analogue of the solvent (viscous) continuous spectrum in the Oldroyd-B fluid. The remaining two $β$-dependent CS are novel features of the FENE-P spectrum, and can have phase speeds outside the base range of velocities, including negative ones. The complexity of the CS predicted here for shearing flows of FENE-P fluids is expected to carry over to other nonlinear viscoelastic models that exhibit a shear-thinning rheology.
The Archimedes Tube Crystallizer (ATC) employs air-segmented flow in coiled tubes to achieve narrow residence time distributions for continuous crystallization. Taylor and Dean vortices drive particle suspension in this system. However, one-way coupled models fail to capture the fluid-particle feedback that becomes critical at higher loadings. We present a particle-resolved Direct Numerical Simulation (DNS) framework based on a Finite Element-Fictitious Boundary Method with hard-contact modeling of particle interactions. Simulations of L-alanine suspensions across varying particle sizes, solid contents, and rotational speeds are validated against experimental side-view imaging. Three quantitative metrics-axial distribution, radial index, and vertical asymmetry-are introduced to classify suspension regimes. The DNS results reproduce the experimentally observed flow map zones (green, yellow, red/yellow, red) and resolve subtle transitions such as rear loading and loss of vertical symmetry. This feasibility study demonstrates that DNS can reliably predict dense suspension behavior and provides a mechanistic foundation for crystallizer design.
Quantum algorithms have been identified as a potential means to accelerate computational fluid dynamics (CFD) simulations, with the lattice Boltzmann method (LBM) being a promising candidate for realizing quantum speedups. Here, we extend the recent quantum algorithm for the incompressible LBM to account for realistic fluid dynamics setups by incorporating walls, inlets, outlets, and external forcing. We analyze the associated complexity cost and show that these modifications preserve the asymptotic scaling, and potential quantum advantage, of the original algorithm. Moreover, to support our theoretical analysis, we provide a classical numerical study illustrating the accuracy, complexity, and convergence of the algorithm for representative incompressible-flow cases, including the driven Taylor-Green vortex, the lid-driven cavity flow, and the flow past a cylinder. Our results provide a pathway to accurate quantum simulation of nonlinear fluid dynamics, and a framework for extending quantum LBM to more challenging flow configurations.
We report measurements of zero- to ultra-low-field nuclear magnetic resonance (ZULF NMR) signals at frequencies of a few hertz with a diamond-based magnetic sensor. The sensing diamond is a truncated pyramid with 0.18 mm height and a 0.5 mm x 0.5mm base. The minimum stand-off distance is < 1 mm, and the sensor sensitivity is 13 pT/(Hz)^(1/2) at frequencies f above 5 Hz with 1/f -like behavior at lower frequencies. NMR signals were generated via signal amplification by reversible exchange (SABRE) parahydrogen-based hyperpolarization resulting in zero-field signals at 1.7 Hz and 3.4 Hz corresponding to the expected hetero-nuclear J-coupling pattern of acetonitrile. This work demonstrates a magnet-free platform for detecting chemically specific NMR signals at ultra-low frequencies paving the way for portable noninvasive diagnostics in microscopic sample volumes for biomedicine, industrial sensing through metal enclosures, and field-deployable quantum analytical devices.
We present a theoretical framework that unifies pump absorption, gain saturation, thermo-optic distortion, and cavity diffraction into a self-consistent model of multi-pass pumped solid-state lasers. By deriving a theoretical formulation of the nonlinear coupling of the superimposed pump energy and effective absorption, we prove an unique steady-state solution exists. Applied to the multi-pass Yb:YAG thin-disk module, the framework is quantitatively validated with experiments, reproducing a well matched absorption tendency, errors in output power, beam diameter and M^2 within 3.0%, 1.7%, and 0.05 respectively. This approach provides predictive guidelines for pump-power scaling and pass-number optimization in high-power lasers.
We consider a bimodal light field envelope propagating in a bulk medium characterized by competing cubic and quintic nonlinearities. The subfields are coupled by a cross-phase modulation term and experience effective attraction. We find dynamically stable stationary states which have two distinct flat-top regions with different intensities. These solutions represent half-vortices, where the first and second components are essentially different and, in particular, carry different topological charges: zero for one component and nonzero for the other. The typical propagation of an unstable half-vortex leads to the splitting of the central vortex core into several fragments which quasielastically interact with the boundary of the flat-top region. This behavior is interpreted as a self-bound solitary wave billiard, where the emerging fragments are the billiard balls and the flat-top region is the dynamically deforming table.
A medium, liquid in its positional degrees of freedom, exhibits in real space a homogeneous single-particle density and an isotropic linear response to perturbations of the density field. Under the hypothesis that, if the anisotropic dipolar interaction is a weak perturbation of the isotropic reference system potential, this condition persists in a dipolar liquid for any dipole configuration, a regularization of the dipolar potential can be introduced. This developments lead to the definition of a screened dipolar interaction that is isotropic, short-ranged, and ferroelectric-like. In a classical density functional theory framework the model is exactly solvable in the limit of infinite dimensions, whereas approximations can be obtained in the case of finite dimensions $d \geq 3$ in the optimized cluster expansion scheme. In both cases, it is shown that the system can support a ferroelectric phase transition. The study relates the emergence of ferroelectricity in dipolar liquids to a constitutive property of the liquid state, i.e. annealed positional disorder.
Nuclear quantum effects (NQEs) arising from the light mass of hydrogen can influence the structure and stability of hydrogen-bonded biomolecules, yet their role in determining peptide and protein folding remains unclear. Experiments show that substituting H$_2$O with D$_2$O often stabilizes folded states, but the microscopic mechanism associated with this phenomena remains unresolved. Through ab initio-level path-integral molecular dynamics simulations enabled by machine-learning interatomic potentials, we address the fundamental question of the role of NQEs in peptides by investigating both their overall impact and isotope substitution effects. Overall, NQEs systematically destabilize compact three-dimensional structures across peptide systems, independent of secondary structure type or side-chain interactions. Contrary to the conventional picture that places central importance on hydrogen bonds, we find that the dominant destabilization instead arises from the quantum C-H vibrations. In addition, we reveal microscopic insights into the stabilization of folded peptides upon H$_2$O to D$_2$O substitution, showing that the H/D isotope substitution of active peptide hydrogens, previously considered unimportant, produces free-energy changes within the range of experimentally observed shifts. These findings provide a new interpretation of isotope effects in biological systems, indicating that seemingly small H$\to$D substitutions within peptides can be as important as, or even outweigh, solvent contributions.
Accurate and efficient simulation of fluid-structure interaction (FSI) problems remains a central challenge in computational physics. High-order discontinuous Galerkin (DG) methods offer low numerical errors and excellent scalability on modern architectures, making them attractive for high-fidelity FSI simulations. This study presents a high-order immersed boundary method (IBM) for FSI problems which combines a volume-penalization approach with a high-order nodal DG solver. To improve near wall accuracy, an anisotropic p-adaptation strategy based on reinforcement learning is used to dynamically adjust the polynomial orders in the mesh elements located near the moving immersed boundaries. By doing so, we show enhanced accuracy with a limited increase in computational cost. Accurate evaluation of surface forces is achieved using symmetric high-order Gaussian quadrature on immersed boundaries. The proposed method is coupled with both rigid-body and elastic-structure solvers within a partitioned framework. Numerical validations using a pitching airfoil, stall flutter of an airfoil, and flow-induced vibration of an elastic beam behind a cylinder demonstrate high-order accuracy and robustness. These results indicate that the present approach provides an effective and scalable strategy for complex moving-boundary FSI simulations.
Developing high throughput applications of sonochemistry and sonoprocessing is an outstanding ultrasonic engineering challenge that continues to limit widespread industrial adoption. Conventional mass-produced Langevin-based technologies, such as the sonotrode or cleaning bath transducers, are not particularly well suited to treating large liquid volumes or flow-based systems, with a compromise between cavitation intensity and distribution through liquid bulk typically required. We report on the development of a tube transducer from a single element radially poled tubular piezoceramic, excited to generate an axially focused field. High-speed imaging and sonochemiluminescence are used to characterise the cavitation generated, which is also compared to the well-known activity at the tip of a sonotrode. Tube transducer and sonotrode sonications are then assessed for the material recycling application of graphite coating delamination from lithium-ion battery anode, both for intact and flaked anode sheets. The findings show that the tube transducer generates cavitation at sonotrode-like intensities or higher but distributed throughout the bore of the tube, with peak activity at the central axis. Prospects for developing tube transducer technology for high throughput flow-based applications are discussed.