Physics
See recent articles
Showing new listings for Monday, 11 November 2024
- [1] arXiv:2411.05008 [pdf, html, other]
-
Title: Electromagnetic stress-energy tensor is a direct consequence of 6-dimensional left-and right-helical representations of the Lorentz groupSubjects: Classical Physics (physics.class-ph); Optics (physics.optics)
We introduce 6-dimensional irreducible vector and helical representations of the Lorentz group especially tailored for transforming 6-component electromagnetic and magnetic field vectors, and we show that mixture of the density matrices of helical states is equivalent to the electromagnetic stress-energy tensor.
- [2] arXiv:2411.05009 [pdf, html, other]
-
Title: A Study of Performance Portability in Plasma Physics SimulationsComments: 15 pages, 8 figures, this is a pre-print to be published in the Latin America High Performance Computing Conference (CARLA) 2024 proceedingsSubjects: Plasma Physics (physics.plasm-ph); Distributed, Parallel, and Cluster Computing (cs.DC); Performance (cs.PF); Computational Physics (physics.comp-ph)
The high-performance computing (HPC) community has recently seen a substantial diversification of hardware platforms and their associated programming models. From traditional multicore processors to highly specialized accelerators, vendors and tool developers back up the relentless progress of those architectures. In the context of scientific programming, it is fundamental to consider performance portability frameworks, i.e., software tools that allow programmers to write code once and run it on different computer architectures without sacrificing performance. We report here on the benefits and challenges of performance portability using a field-line tracing simulation and a particle-in-cell code, two relevant applications in computational plasma physics with applications to magnetically-confined nuclear-fusion energy research. For these applications we report performance results obtained on four HPC platforms with server-class CPUs from Intel (Xeon) and AMD (EPYC), and high-end GPUs from Nvidia and AMD, including the latest Nvidia H100 GPU and the novel AMD Instinct MI300A APU. Our results show that both Kokkos and OpenMP are powerful tools to achieve performance portability and decent "out-of-the-box" performance, even for the very latest hardware platforms. For our applications, Kokkos provided performance portability to the broadest range of hardware architectures from different vendors.
- [3] arXiv:2411.05015 [pdf, html, other]
-
Title: Phase-field models for ductile fatigue fractureSubjects: Computational Physics (physics.comp-ph)
Fatigue fracture is one of the main causes of failure in structures. However, the simulation of fatigue crack growth is computationally demanding due to the large number of load cycles involved. Metals in the low cycle fatigue range often show significant plastic zones at the crack tip, calling for elastic-plastic material models, which increase the computation time even further. In pursuit of a more efficient model, we propose a simplified phase-field model for ductile fatigue fracture, which indirectly accounts for plasticity within the fatigue damage accumulation. Additionally, a cycle-skipping approach is inherent to the concept, reducing computation time by up to several orders of magnitude.
Essentially, the proposed model is a simplification of a phase-field model with elastic-plastic material behavior. As a reference, we therefore implement a conventional elastic-plastic phase-field fatigue model with nonlinear hardening and a fatigue variable based on the strain energy density, and compare the simplified model to it. Its approximation of the stress-strain behavior, the neglect of the plastic crack driving force and consequential range of applicability are discussed.
Since in fact the novel efficient model is similar in its structure to a phase-field fatigue model we published in the past, we include this older version in the comparison, too. Compared to this model variant, the novel model improves the approximation of the plastic strains and corresponding stresses and refines the damage computation based on the Local Strain Approach. For all model variants, experimentally determined values for elastic, plastic, fracture and fatigue properties of AA2024 T351 aluminum sheet material are employed. - [4] arXiv:2411.05018 [pdf, html, other]
-
Title: Biquaternion Signal Processing for Nonlinear UltrasonicsComments: 22 pages, 12 figuresSubjects: Optics (physics.optics); High Energy Physics - Lattice (hep-lat)
Localization and classification of scattered nonlinear ultrasonic signatures in 2 dimensional complex damaged media using Time Reversal based Nonlinear Elastic Wave Spectroscopy (TR-NEWS) approach is extended to 3 dimensional complex damaged media. In (2+1)D, i.e. space 2 dimensional time 1 dimensional spacetime, we used quaternion bases for analyses, while in (3+1)D, we use biquaternion bases.
The optimal weight function of the path of ultrasonic wave in (3+1)D lattice is obtained by using the Echo State Network (ESN) which is a Machine Learning technique. The hysteresis effect is incorporated by using the Preisach-Mayergoyz model. - [5] arXiv:2411.05019 [pdf, html, other]
-
Title: Enhancing Accuracy and Feature Insights in Hydration Free Energy Predictions for Small Molecules with Machine LearningSubjects: Chemical Physics (physics.chem-ph); Computational Physics (physics.comp-ph)
The accurate prediction of solvation free energy is of significant importance as it governs the behavior of solutes in solution. In this work, we apply a variety of machine learning techniques to predict and analyze the alchemical free energy of small molecules. Our methodology incorporates an ensemble of machine learning models with feature processing using the K-nearest neighbors algorithm. Two training strategies are explored: one based on experimental data, and the other based on the offset between molecular dynamics (MD) simulations and experimental measurements. The latter approach yields a substantial improvement in predictive accuracy, achieving a mean unsigned error (MUE) of 0.64 kcal/mol. Feature analysis identifies molecular geometry and topology as the most critical factors in predicting alchemical free energy, supporting the established theory that surface tension is a key determinant. Furthermore, the feature analysis of offset results highlights the relevance of charge distribution within the system, which correlates with the inaccuracies in force fields employed in MD simulations and may provide guidance for improving force field designs. These results suggest that machine learning approaches can effectively capture the complex features governing solvation free energy, offering novel pathways for enhancing predictive accuracy.
- [6] arXiv:2411.05021 [pdf, other]
-
Title: Considerations and recommendations from the ISMRM Diffusion Study Group for preclinical diffusion MRI: Part 3 -- Ex vivo imaging: data processing, comparisons with microscopy, and tractographyKurt G Schilling, Amy FD Howard, Francesco Grussu, Andrada Ianus, Brian Hansen, Rachel L C Barrett, Manisha Aggarwal, Stijn Michielse, Fatima Nasrallah, Warda Syeda, Nian Wang, Jelle Veraart, Alard Roebroeck, Andrew F Bagdasarian, Cornelius Eichner, Farshid Sepehrband, Jan Zimmermann, Lucas Soustelle, Christien Bowman, Benjamin C Tendler, Andreea Hertanu, Ben Jeurissen, Marleen Verhoye, Lucio Frydman, Yohan van de Looij, David Hike, Jeff F Dunn, Karla Miller, Bennett A Landman, Noam Shemesh, Adam Anderson, Emilie McKinnon, Shawna Farquharson, Flavio Dell' Acqua, Carlo Pierpaoli, Ivana Drobnjak, Alexander Leemans, Kevin D Harkins, Maxime Descoteaux, Duan Xu, Hao Huang, Mathieu D Santin, Samuel C. Grant, Andre Obenaus, Gene S Kim, Dan Wu, Denis Le Bihan, Stephen J Blackband, Luisa Ciobanu, Els Fieremans, Ruiliang Bai, Trygve B Leergaard, Jiangyang Zhang, Tim B Dyrby, G Allan Johnson, Julien Cohen-Adad, Matthew D Budde, Ileana O JelescuComments: Part 3 of 3 in "Considerations and recommendations for preclinical diffusion MRI"Subjects: Medical Physics (physics.med-ph)
Preclinical diffusion MRI (dMRI) has proven value in methods development and validation, characterizing the biological basis of diffusion phenomena, and comparative anatomy. While dMRI enables in vivo non-invasive characterization of tissue, ex vivo dMRI is increasingly being used to probe tissue microstructure and brain connectivity. Ex vivo dMRI has several experimental advantages that facilitate high spatial resolution and high signal-to-noise ratio (SNR) images, cutting-edge diffusion contrasts, and direct comparison with histological data as a methodological validation. However, there are a number of considerations that must be made when performing ex vivo experiments. The steps from tissue preparation, image acquisition and processing, and interpretation of results are complex, with many decisions that not only differ dramatically from in vivo imaging of small animals, but ultimately affect what questions can be answered using the data. This work concludes a 3-part series of recommendations and considerations for preclinical dMRI. Herein, we describe best practices for dMRI of ex vivo tissue, with a focus on image pre-processing, data processing and model fitting, and tractography. In each section, we attempt to provide guidelines and recommendations, but also highlight areas for which no guidelines exist (and why), and where future work should lie. We end by providing guidelines on code sharing and data sharing, and point towards open-source software and databases specific to small animal and ex vivo imaging.
- [7] arXiv:2411.05033 [pdf, html, other]
-
Title: Diagonalization without Diagonalization: A Direct Optimization Approach for Solid-State Density Functional TheoryTianbo Li, Min Lin, Stephen Dale, Zekun Shi, A. H. Castro Neto, Kostya S. Novoselov, Giovanni VignaleSubjects: Chemical Physics (physics.chem-ph); Machine Learning (cs.LG); Computational Physics (physics.comp-ph)
We present a novel approach to address the challenges of variable occupation numbers in direct optimization of density functional theory (DFT). By parameterizing both the eigenfunctions and the occupation matrix, our method minimizes the free energy with respect to these parameters. As the stationary conditions require the occupation matrix and the Kohn-Sham Hamiltonian to be simultaneously diagonalizable, this leads to the concept of ``self-diagonalization,'' where, by assuming a diagonal occupation matrix without loss of generality, the Hamiltonian matrix naturally becomes diagonal at stationary points. Our method incorporates physical constraints on both the eigenfunctions and the occupations into the parameterization, transforming the constrained optimization into an fully differentiable unconstrained problem, which is solvable via gradient descent. Implemented in JAX, our method was tested on aluminum and silicon, confirming that it achieves efficient self-diagonalization, produces the correct Fermi-Dirac distribution of the occupation numbers and yields band structures consistent with those obtained with SCF methods in Quantum Espresso.
- [8] arXiv:2411.05035 [pdf, html, other]
-
Title: FLAIM: A reduced volume ignition model for the compression and thermonuclear burn of spherical fuel capsulesAbd Essamade Saufi, Hannah Bellenbaum, Martin Read, Nicolas Niasse, Sean Barrett, Nicholas Hawker, Nathan Joiner, David ChapmanComments: 60 pages, 24 figures, part of the special issue "12th International Conference on Inertial Fusion Science and Applications"Journal-ref: Volume 53, December 2024, 101159Subjects: Plasma Physics (physics.plasm-ph)
We present the "First Light Advanced Ignition Model" (FLAIM), a reduced model for the implosion, adiabatic compression, volume ignition and thermonuclear burn of a spherical DT fuel capsule utilising a high-Z metal pusher. FLAIM is characterised by a highly modular structure, which makes it an appropriate tool for optimisations, sensitivity analyses and parameter scans. One of the key features of the code is the 1D description of the hydrodynamic operator, which has a minor impact on the computational efficiency, but allows us to gain a major advantage in terms of physical accuracy. We demonstrate that a more accurate treatment of the hydrodynamics plays a primary role in closing most of the gap between a simple model and a general 1D rad-hydro code, and that only a residual part of the discrepancy is attributable to the heat losses. We present a detailed quantitative comparison between FLAIM and 1D rad-hydro simulations, showing good agreement over a large parameter space in terms of temporal profiles of key physical quantities, ignition maps and typical burn metrics.
- [9] arXiv:2411.05053 [pdf, html, other]
-
Title: Characterization of the LUNA neutron detector array for the measurement of the 13C(a,n)16O reactionL. Csedreki, G.F. Ciani, J. Balibrea-Correa, A. Best, M. Aliotta, F. Barile, D. Bemmerer, A. Boeltzig, C.Broggini, C.G. Bruno, A. Caciolli, F. Cavanna, T. Chillery, P. Colombetti, P. Corvisiero, T. Davinson, R. Depalo, A. Di Leva, Z. Elekes, F. Ferraro, E. M. Fiore, A. Formicola, Zs. Fulop, G. Gervino, A. Guglielmetti, C. Gustavino, Gy. Gyurky, G. Imbriani, Z. Janas, M. Junker, I. Kochanek, M. Lugaro, P. Marigo, E. Masha, C. Mazzocchi, R. Menegazzo, V. Mossa, F. R. Pantaleo, V.Paticchio, R. Perrino, D.Piatti, P. Prati, L. Schiavulli, K. Stockel, O. Straniero, T. Szucs, M. P. Takacs, F.Terrasi, S. ZavatarelliComments: 8 pages, 5 figures, to be published in NIMAJournal-ref: Nuclear Inst. and Methods in Physics Research, A 994 (2021) 165081Subjects: Instrumentation and Detectors (physics.ins-det); Nuclear Experiment (nucl-ex)
We introduce the LUNA neutron detector array developed for the investigation of the 13C(a,n)16O reaction towards its astrophysical s-process Gamow peak in the low-background environment of the Laboratori Nazionali del Gran Sasso (LNGS). Eighteen 3He counters are arranged in two different configurations (in a vertical and a horizontal orientation) to optimize neutron detection effciency, target handling and target cooling over the investigated energy range Ea;lab = 300 - 400 keV (En = 2.2 - 2.6 MeV in emitted neutron energy). As a result of the deep underground location, the passive shielding of the setup and active background suppression using pulse shape discrimination, we reached a total background rate of 1.23 +- 0.12 counts/hour. This resulted in an improvement of two orders of magnitude over the state of the art allowing a direct measurement of the 13C(a,n)16O cross-section down to Ea;lab = 300 keV. The absolute neutron detection efficiency of the setup was determined using the 51V(p,n)51Cr reaction and an AmBe radioactive source, and completed with a Geant4 simulation. We determined a (34+-3) % and (38+-3) % detection efficiency for the vertical and horizontal configurations, respectively, for En = 2.4 MeV neutrons.
- [10] arXiv:2411.05090 [pdf, html, other]
-
Title: Rarefied neutral flow in orificed hollow cathodesComments: 25 pages, 27 figures, 4 tables. Submitted to Journal of Applied PhysicsSubjects: Fluid Dynamics (physics.flu-dyn); Plasma Physics (physics.plasm-ph)
A parametric study is conducted to quantify the effect of the keeper electrode geometry on the neutral flow quantities within orificed hollow cathodes, prior to cathode ignition. The keeper impinges directly on the flow out of the cathode orifice and its geometry strongly influences the product between the pressure in the orifice-keeper region and the distance between cathode and keeper, $P_{ko}\cdot D_{ko}$, a key parameter for successful cathode ignition. A representative cathode equipped with a keeper is simulated using the Direct Simulation Monte Carlo method. The numerical model is first validated with computational results from the literature, and a parametric study is then conducted. Parameters include the cathode pressure-diameter in the range of 1-5 Torr-cm and the following geometric ratios: cathode orifice-to-cathode inner radii in the range of 0.1-0.7, keeper orifice-to-cathode orifice radii in the range of 1-5, and keeper distance-to-cathode-orifice diameter, in the range of 0.5-10. If both keeper and cathode have identical orifice radii, the flow remains subsonic in the cathode orifice-to-keeper region. In most cases, however, the flow becomes underexpanded and supersonic, and the static pressure within the cathode orifice-keeper region is, on average, 4% that of the upstream pressure value. The orifice-keeper region pressure increases with either a decrease in the keeper orifice diameter or an increase in the distance between the cathode orifice plate and the keeper plate, in agreement with literature data. Trends are explained through control-volume-based conservation laws. The ratio of ignition-to-nominal mass flow rates is found to be in the range of 9-120, with a most probable value of 40, in agreement with literature data. This suggests that heater-less cathode ignition at a minimum DC voltage may be achieved by increasing the input mass flow rate by a factor of 40.
- [11] arXiv:2411.05094 [pdf, other]
-
Title: Experimental Investigation of Variations in Polycrystalline Hf0.5Zr0.5O2 (HZO)-based MFIMComments: 6 pages, 8 figuresSubjects: Applied Physics (physics.app-ph); Materials Science (cond-mat.mtrl-sci)
Device-to-device variations in ferroelectric (FE) hafnium oxide-based devices pose a crucial challenge that limits the otherwise promising capabilities of this technology. Earlier simulation-based studies have identified polarization (P) domain nucleation and polycrystallinity as key contributors to these variations. In this work, we experimentally investigate the effect of these two factors on remanent polarization (PR) variation in Hf0.5Zr0.5O2 (HZO) based metal-ferroelectric-insulator-metal (MFIM) capacitors for different set voltages (VSET) and FE thicknesses (TFE). Our measurements reveal a non-monotonic behavior of PR variations with VSET, which is consistent with previous simulation-based predictions. For low and high-VSET regions, we find that PR variations are dictated primarily by saturation polarization (PS) variations and are associated with the polycrystallinity in HZO. Our measurements also reveal that PR variations peak near the coercive voltage (VC), defined as the mid-VSET region. We attribute the increase of PR variation around VC to the random nature and sharp P switching associated with domain nucleation, which is dominant near VC. Further, we observe a reduction in the peak PR variation as HZO thickness (TFE) is scaled. We validate our arguments by establishing the correlation between the measured values of PR with VC and PS. Our results display that a strong correlation exists between PR and VC in the mid-VSET region and between PR and PS in the low and high-VSET regions across various TFE.
- [12] arXiv:2411.05168 [pdf, html, other]
-
Title: Disentangling enhanced diffusion and ballistic motion of excitons coupled to Bloch surface waves with molecular dynamics simulationsComments: 23 pages, 3 figuresSubjects: Chemical Physics (physics.chem-ph)
Placing an organic material on top of a Bragg mirror can significantly enhance exciton transport. Such enhancement has been attributed to strong coupling between the evanescent Bloch surface waves (BSW) on the mirror, and the excitons in the material. In this regime, the BSW and excitons hybridize into Bloch surface wave polaritons (BSWP), new quasi-particles with both photonic and excitonic character. While recent experiments unveiled a mixed nature of the enhanced transport, the role of the material degrees of freedom in this process remains unclear. To clarify their role, we performed atomistic molecular dynamics simulations of an ensemble of Methylene blue molecules, a prototype organic emitter, strongly coupled to a BSW. In contrast to the established static models of polaritons, even with disorder included, our dynamic simulations reveal a correlation between the photonic content of the BSWP and the nature of the transport. In line with experiment, we find ballistic motion for polaritons with high photonic character, and enhanced diffusion if the photonic content is low. Our simulations furthermore suggest that the diffusion is due to thermally activated vibrations that drive population transfer between the stationary dark states and mobile bright polaritonic states.
- [13] arXiv:2411.05202 [pdf, html, other]
-
Title: Particle Levitation Velocimetry for boundary layer measurements in high Reynolds number liquid helium turbulenceComments: 12 pages, 6 figuresSubjects: Fluid Dynamics (physics.flu-dyn); Other Condensed Matter (cond-mat.other)
Understanding boundary layer flows in high Reynolds number (Re) turbulence is crucial for advancing fluid dynamics in a wide range of applications, from improving aerodynamic efficiency in aviation to optimizing energy systems in industrial processes. However, generating such flows requires complex, power-intensive large-scale facilities. Furthermore, the use of local probes, such as hot wires and pressure sensors, often introduces disturbances due to the necessary support structures, compromising measurement accuracy. In this paper, we present a solution that leverages the vanishingly small viscosity of liquid helium to produce high Re flows, combined with an innovative Particle Levitation Velocimetry (PLV) system for precise flow-field measurements. This PLV system uses magnetically levitated superconducting micro-particles to measure the near-wall velocity field in liquid helium. Through comprehensive theoretical analysis, we demonstrate that the PLV system enables quantitative measurements of the velocity boundary layer over a wall unit range of $44\le y^{+}\le 4400$, with a spatial resolution that, depending on the particle size, can reach down to about 10~$\mu$m. This development opens new avenues for exploring turbulence structures and correlations within the thin boundary layer that would be otherwise difficult to achieve.
- [14] arXiv:2411.05240 [pdf, html, other]
-
Title: Reducing data resolution for better super-resolution: Reconstructing turbulent flows from noisy observationSubjects: Fluid Dynamics (physics.flu-dyn); Dynamical Systems (math.DS)
A super-resolution (SR) method for the reconstruction of Navier-Stokes (NS) flows from noisy observations is presented. In the SR method, first the observation data is averaged over a coarse grid to reduce the noise at the expense of losing resolution and, then, a dynamic observer is employed to reconstruct the flow field by reversing back the lost information. We provide a theoretical analysis, which indicates a chaos synchronization of the SR observer with the reference NS flow. It is shown that, even with noisy observations, the SR observer converges toward the reference NS flow exponentially fast, and the deviation of the observer from the reference system is bounded. Counter-intuitively, our theoretical analysis shows that the deviation can be reduced by increasing the lengthscale of the spatial average, i.e., making the resolution coarser. The theoretical analysis is confirmed by numerical experiments of two-dimensional NS flows. The numerical experiments suggest that there is a critical lengthscale for the spatial average, below which making the resolution coarser improves the reconstruction.
- [15] arXiv:2411.05249 [pdf, html, other]
-
Title: A High-Order Analytical Extension of the Corrected Smagorinsky Model for Non-Equilibrium Turbulent FlowComments: 25 pagesSubjects: Fluid Dynamics (physics.flu-dyn); Mathematical Physics (math-ph)
This study presents an extension of the corrected Smagorinsky model, incorporating advanced techniques for error estimation and regularity analysis of far-from-equilibrium turbulent flows. A new formulation that increases the model's ability to explain complex dissipative processes in turbulence is presented, using higher-order Sobolev spaces to address incompressible and compressible Navier-Stokes equations. Specifically, a refined energy dissipation mechanism that provides a more accurate representation of turbulence is introduced, particularly in the context of multifractal flow regimes. Furthermore, we derive new theoretical results on energy regularization in multifractal turbulence, contributing to the understanding of anomalous dissipation and vortex stretching in turbulent flows. The work also explores the numerical implementation of the model in the presence of challenging boundary conditions, particularly in dynamically evolving domains, where traditional methods struggle to maintain accuracy and stability. Theoretical demonstrations and analytical results are provided to validate the proposed framework, with implications for theoretical advances and practical applications in computational fluid dynamics. This approach provides a basis for more accurate simulations of turbulence, with potential applications ranging from atmospheric modeling to industrial fluid dynamics.
- [16] arXiv:2411.05256 [pdf, html, other]
-
Title: Radiopurity measurements of liquid scintillator for the COSINE-100 UpgradeJ. Kim, C. Ha, S.H. Kim, W.K. Kim, Y.D. Kim, Y.J. Ko, E.K. Lee, H. Lee, H.S. Lee, I.S. Lee, J. Lee, S.H. Lee, S.M. Lee, Y.J. Lee, G.H. YuSubjects: Instrumentation and Detectors (physics.ins-det); High Energy Physics - Experiment (hep-ex)
A new 2,400 L liquid scintillator has been produced for the COSINE-100 Upgrade, which is under construction at Yemilab for the next COSINE dark matter experiment phase. The linear-alkyl-benzene-based scintillator is designed to serve as a veto for NaI(Tl) crystal targets and a separate platform for rare event searches. We measured using a sample consisting of a custom-made 445 mL cylindrical Teflon container equipped with two 3-inch photomultiplier tubes. Analyses show activity levels of $0.091 \pm 0.042$ mBq/kg for $^{238}$U and $0.012 \pm 0.007$ mBq/kg for $^{232}$Th.
- [17] arXiv:2411.05264 [pdf, html, other]
-
Title: A Stochastic Conceptual Model for the Coupled ENSO and MJOComments: 26 pages, 12 figuresSubjects: Atmospheric and Oceanic Physics (physics.ao-ph); Dynamical Systems (math.DS)
Understanding the interactions between the El Nino-Southern Oscillation (ENSO) and the Madden-Julian Oscillation (MJO) is essential to studying climate variabilities and predicting extreme weather events. Here, we develop a stochastic conceptual model for describing the coupled ENSO-MJO phenomenon. The model adopts a three-box representation of the interannual ocean component to characterize the ENSO diversity. For the intraseasonal atmospheric component, a low-order Fourier representation is used to describe the eastward propagation of the MJO. We incorporate decadal variability to account for modulations in the background state that influence the predominant types of El Nino events. In addition to dynamical coupling through wind forcing and latent heat, state-dependent noise is introduced to characterize the statistical interactions among these multiscale processes, improving the simulation of extreme events. The model successfully reproduces the observed non-Gaussian statistics of ENSO diversity and MJO spectra. It also captures the interactions between wind, MJO, and ENSO.
- [18] arXiv:2411.05284 [pdf, other]
-
Title: Intra-operative Optimal Flow Diverter Selection for Intracranial aneurysm treatment using Angiographic parametric imaging: Feasibility studyJournal-ref: Proc. SPIE 12468, Medical Imaging 2023: Biomedical Applications in Molecular, Structural, and Functional Imaging, 124680P (10 April 2023)Subjects: Medical Physics (physics.med-ph)
During intracranial aneurysm (IA) treatment with Diverters (FDs), the device/parent artery diameters ratio may influence the ability of the device to induce aneurysm healing response. Oversized FDs are safer to deploy but may not induce enough hemodynamic resistance to ensure aneurysm occlusion. Methods based on Computational Fluid Dynamics (CFD) could allow optimal device selection but are time-consuming and inadequate for intra-operative guidance. To address this limitation, we propose to investigate a method for optimal FD selection using Angiographic Parametric Imaging (API) and machine learning (ML). We selected 128 pre-treatment angiographic sequences of IAs which demonstrated full occlusion at six months follow-up. For each IA, we extracted five API parameters from the aneurysm dome and normalized them to the feeding artery corresponding parameters. We dichotomized the dataset based on the FD/ proximal artery diameter ratio as undersized, if the ratio<1 or if multiple FDs were used and oversized otherwise. Single API parameter and ML analysis were used to determine whether API parameters could be used to determine the need for FD under-sizing (i.e., increased flow resistance). Classification accuracy was assessed using area under the receiver operator characteristic (AUROC). In total we identified 51 and 77 cases for the undersized and oversized cohorts respectively. Single API parameter analysis yielded an inadequate AUROC ~0.5 while machine learning using all five API parameters yielded and AUROC of 0.72.
- [19] arXiv:2411.05287 [pdf, other]
-
Title: In-Silico Analysis of Curve Fitting in Angiographic Parametric Imaging in Intracranial Aneurysms to Reduce Patient and Interventionalist-Induced ErrorsJournal-ref: Medical Imaging 2024: Physics of Medical Imaging, 139, (2024)Subjects: Medical Physics (physics.med-ph)
In Angiographic Parametric Imaging (API), accurate estimation of parameters from Time Density Curves (TDC) is crucial. However, these estimations are often marred by errors arising from factors such as patient motion, procedural preferences, image noise, and injection variability. While fitting methods like gamma-variate fitting offer a solution to recover incomplete or corrupted TDC data, they might also introduce unforeseen biases. This study investigates the trade-offs and benefits of employing gamma-variate fitting on virtual angiograms to enhance the precision of API biomarkers. Utilizing Computational Fluid Dynamics (CFD) in patient specific 3D geometries, we generated a series of high-definition virtual angiograms at distinct inlet velocities: 0.25m/s, 0.35m/s, and 0.45m/s. These velocities were investigated across injection durations ranging from 0.5s to 2.0s. From these angiograms, TDCs for aneurysms and their corresponding inlets were constructed. To emulate typical clinical challenges, we introduced noise, simulated patient motion, and generated temporally incomplete data sets. These modified TDCs underwent gamma-variate fitting. We quantified both the original and fitted TDC curves using standard angiography metrics such as Cross-Correlation (Cor), Time to Peak (TTP), Mean Transit Time (MTT), Peak Height (PH), Area Under the Curve (AUC), and Maximum Gradient (Max-Gr) for a comprehensive comparison. TDCs enhanced by gamma-variate fitting exhibited a robust correlation with vascular flow dynamics. Our results affirm that gamma-variate fitting can adeptly restore TDCs from fragmentary sequences, elevating the precision of derived API parameters.
- [20] arXiv:2411.05321 [pdf, html, other]
-
Title: Trapping of Single Atoms in Metasurface Optical Tweezer ArraysComments: 14 pages, 5 main figures, 8 extended data figuresSubjects: Atomic Physics (physics.atom-ph); Quantum Gases (cond-mat.quant-gas); Applied Physics (physics.app-ph); Optics (physics.optics); Quantum Physics (quant-ph)
Optical tweezer arrays have emerged as a key experimental platform for quantum computation, quantum simulation, and quantum metrology, enabling unprecedented levels of control over single atoms and molecules. Existing methods to generate tweezer arrays mostly rely on active beam-shaping devices, such as acousto-optic deflectors or liquid-crystal spatial light modulators. However, these approaches have fundamental limitations in array geometry, size, and scalability. Here we demonstrate the trapping of single atoms in optical tweezer arrays generated via holographic metasurfaces. We realize two-dimensional arrays with more than 250 tweezer traps, arranged in arbitrary geometries with trap spacings as small as 1.5 um. The arrays have a high uniformity in terms of trap depth, trap frequency, and positional accuracy, rivaling or exceeding existing approaches. Owing to sub-micrometer pixel sizes and high pixel densities, holographic metasurfaces open a path towards optical tweezer arrays with >100,000 traps.
- [21] arXiv:2411.05326 [pdf, other]
-
Title: Implementation of spectral methods on Ising machines: toward flow simulations on quantum annealerComments: 26 pages, 8 figuresJournal-ref: Fluid Dynamics Research (2024)Subjects: Fluid Dynamics (physics.flu-dyn)
We investigate the possibility and current limitations of flow computations using quantum annealers by solving a fundamental flow problem on Ising machines. As a fundamental problem, we consider the one-dimensional advection-diffusion equation. We formulate it in a form suited to Ising machines (i.e., both classical and quantum annealers), perform extensive numerical tests on a classical annealer, and finally test it on an actual quantum annealer. To make it possible to process with an Ising machine, the problem is formulated as a minimization problem of the residual of the governing equation discretized using either the spectral method or the finite difference method. The resulting system equation is then converted to the Quadratic Unconstrained Binary Optimization (QUBO) form though quantization of variables. It is found in the numerical tests using a classical annealer that the spectral method requiring smaller number of variables has a particular merit over the finite difference method because the accuracy deteriorates with the increase of the number of variables. We also found that the computational error varies depending on the condition number of the coefficient matrix. In addition, we extended it to a two-dimensional problem and confirmed its fundamental applicability. From the numerical test using a quantum annealer, however, it turns out that the computation using a quantum annealer is still challenging due largely to the structural difference from the classical annealer, which leaves a number of issues toward its practical use.
- [22] arXiv:2411.05327 [pdf, html, other]
-
Title: Social balance in directed networksSubjects: Physics and Society (physics.soc-ph)
Social networks inherently exhibit complex relationships that can be positive or negative, as well as directional. Understanding balance in these networks is crucial for unraveling social dynamics, yet traditional theories struggle to incorporate directed interactions. This perspective presents a comprehensive roadmap for understanding balance in signed directed networks, extending traditional balance theory to account for directed interactions. Balance is indicated by the enrichment of higher-order patterns like triads compared to an adequate null model, where the network is randomized with some key aspects being preserved. Finding appropriate null models has been a challenging task even without considering directionality, which largely expands the space of potential null models. Recently, it has been shown that in the undirected case both the network topology and the signed degrees serve as key factors to preserve. Therefore, we introduce a maximally constrained null model that preserves the directed topology as well as node-level features given by the signed unidirectional, reciprocated, and conflicting node degrees. Our null model is based on the maximum-entropy principle and reveals consistent patterns across large-scale social networks. We also consider directed generalizations of balance theory and find that the observed patterns are well aligned with two proposed directed notions of strong balance. Our approach not only unveils balance in signed directed networks but can also serve as a starting point towards generative models of signed directed social networks, advancing our understanding of complex social systems and their dynamics.
- [23] arXiv:2411.05336 [pdf, html, other]
-
Title: Backscattering-Immune Floquet Conversion in Ring ModulatorsComments: 5 pages, 4 figures, LetterSubjects: Optics (physics.optics)
Backscattering in micro-ring cavities induces mode mixing and limits device performance. Existing methods to mitigate backscattering often involve complex fabrication processes or are insufficient for complete suppression. In this work, we introduce a novel method to eliminate backscattering by operating the cavity at an exceptional point (EP). By engineering non-conservative coupling between degenerate clockwise (CW) and counter-clockwise (CCW) modes, we achieve chiral transmission that prevents degeneracy lifting and suppresses unwanted mode coupling. Unlike previous approaches that rely on precise gain-loss balance or complex structures, our method utilizes non-conservative coupling between the counterpropgating cavity modes. Using this method, we further show significant enhancement in the cavity performance in Floquet mode conversion efficiency at the EP. Our highly adaptable approach enables seamless integration into various photonic platforms with electro-optic modulators. This advancement mitigates backscattering and improves the precision of light-matter interactions, offering promising applications in quantum communication and information processing.
- [24] arXiv:2411.05339 [pdf, html, other]
-
Title: Spectral Width of Maximum Deposition Eigenchannels in Diffusive MediaRohin E. McIntosh, Arthur Goetschy, Nicholas Bender, Alexey Yamilov, Chia Wei Hsu, Hasan Yilmaz, Hui CaoComments: 13 pages, 5 figuresSubjects: Optics (physics.optics)
The maximum deposition eigenchannel provides the largest possible power delivery to a target region inside a diffusive medium by optimizing the incident wavefront of a monochromatic beam. It originates from constructive interference of scattered waves, which is frequency sensitive. We investigate the spectral width of maximum deposition eigenchannels over a range of target depths using numerical simulations of a 2D diffusive system. Compared to tight focusing into the system, power deposition to an extended region is more sensitive to frequency detuning. The spectral width of enhanced delivery to a large target displays a rather weak, non-monotonic variation with target depth, in contrast to a sharp drop of focusing bandwidth with depth. While the maximum enhancement of power deposited within a diffusive system can exceed that of power transmitted through it, this comes at the cost of a narrower spectral width. We investigate the narrower deposition width in terms of the constructive interference of transmission eigenchannels within the target. We further observe that the spatial field distribution inside the target region decorrelates slower with spectral detuning than power decay of the maximum deposition eigenchannel. Additionally, absorption increases the spectral width of deposition eigenchannels, but the depth dependence remains qualitatively identical to that without absorption. These findings hold for any diffusive waves, including electromagnetic waves, acoustic waves, pressure waves, mesoscopic electrons, and cold atoms.
- [25] arXiv:2411.05351 [pdf, html, other]
-
Title: Extreme narrow band in Moir\'e Photonic time crystalSubjects: Optics (physics.optics)
The Moiré superlattice has attracted growing interest in the electromagnetic and optical communities. Here, we extend this concept to time-varying photonic systems by superposing two binary modulations on the refractive index with different modulation periods, i.e., the Moiré photonic time crystal (PTC). Such a Moiré PTC leads to extreme narrow bands in momentum space which supports temporal localized modes, exhibiting periodically self-reconstructing pulse in time domain. We investigate how the modulation parameters change the bandstructure of the Moiré PTC and the temporal localization behavior. Moreover, we explore mode-locking mechanism in frequency space in the Moiré PTC, which points towards potential applications in mode-locked lasers with tunable time width of the emitted pulses. Our work therefore extends the study of PTC to complex modulation patterns, and unveils new possibility in wave manipulation with time-varying systems.
- [26] arXiv:2411.05352 [pdf, html, other]
-
Title: An atomic Fabry-Perot interferometer-based acceleration sensor for microgravity environmentsComments: 17 pages, 9 figuresSubjects: Atomic Physics (physics.atom-ph); Quantum Physics (quant-ph)
We investigate the use of an atomic Fabry-Perot interferometer (FPI) with a pulsed non-interacting Bose-Einstein condensate (BEC) source as a space-based acceleration sensor. We derive an analytic approximation for the device's transmission under a uniform acceleration, which we use to compute the device's attainable acceleration sensitivity using the classical Fisher information. In the ideal case of a high-finesse FPI and an infinitely narrow momentum width atomic source, we find that when the total length of the device is constrained to small values, the atomic FPI can achieve greater acceleration sensitivity than a Mach-Zender (MZ) interferometer of equivalent total device length. Under the more realistic case of a finite momentum width atomic source, We identify the ideal cavity length that gives the best sensitivity. Although the MZ interferometer now offers enhanced sensitivity within currently-achievable experimental parameter regimes, our analysis demonstrates that the atomic FPI holds potential as a promising alternative in the future, provided that narrow momentum width atomic sources can be engineered.
- [27] arXiv:2411.05380 [pdf, html, other]
-
Title: Exploring Quantum Entanglement: A Tutorial on Spin Dynamics and Entanglement Transfer with Python ModulesComments: 23 pages, 7 figures; python modules available from the authors upon requestSubjects: Physics Education (physics.ed-ph)
Quantum entanglement is a captivating phenomenon in quantum physics, characterized by intricate and non-classical correlations between particles. This phenomenon plays a crucial role in quantum computing and measurement processes. In this tutorial we explore the dynamics of quantum systems with up to three spins, providing an introductory guide to understanding how entanglement evolves and transfers within such systems. Through detailed examples, simulations, and analyses, the tutorial offers insights into the fundamental principles of entanglement; . We also provide \code{python} modules for reproducing the presented results and as a basis for further projects. The target audience of this tutorial is physics enthusiasts among high school students and students in their first semesters.
- [28] arXiv:2411.05392 [pdf, html, other]
-
Title: Long-term stable laser injection locking for quasi-CW applicationsSubjects: Optics (physics.optics); Quantum Gases (cond-mat.quant-gas); Atomic Physics (physics.atom-ph)
Generating high output powers while achieving narrow line single mode lasing are often mutual exclusive properties of commercial laser diodes. For this reason, efficient and scalable amplification of narrow line laser light is still a major driving point in modern laser system designs. Commonly, injection locking of high-power semiconductor laser diodes are used for this purpose. However, for many laser diodes it is very challenging to achieve stable operation of the injection locked state due to a complex interplay of non-linearities and thermal effects. Different approaches of active or passive stabilization usually require a large overhead of optical and electrical equipment and are not generally applicable. In our work we present a passive stabilization scheme, that is generally applicable, technically easy to implement and extremely cost-effective. It is based on the externally synchronized automatic acquisition of the optimal injection state. Central to our simple but powerful scheme is the management of thermalization effects during lock acquisition. By periodical relocking, spectrally pure amplified light is maintained in a quasi-CW manner over long timescales. We characterize the performance of our method for laser diodes amplifying 671 nm light and demonstrate the general applicability by confirming the method to work also for laser diodes at 401 nm, 461 nm and 689 nm. Our scheme enables the scaled operation of injection locks, even in cascaded setups, for the distributed amplification of single frequency laser light.
- [29] arXiv:2411.05394 [pdf, html, other]
-
Title: Relativistic and electron-correlation effects in static dipole polarizabilities for group 12 elementsComments: arXiv admin note: substantial text overlap with arXiv:2410.01493Subjects: Atomic Physics (physics.atom-ph)
In this study, we report a detailed calculation of the static dipole polarizabilities for group 12 elements using the finite-field approach combined with the relativistic coupled-cluster method, including single, double, and perturbative triple excitations. We examine three types of relativistic effects on dipole polarizabilities: scalar-relativistic, spin-orbit coupling (SOC), and fully relativistic Dirac-Coulomb contributions. The final recommended polarizability values, with their associated uncertainties, are $37.95 \pm 0.72$ for Zn, $45.68 \pm 1.16$ for Cd, $34.04 \pm 0.67$ for Hg, and $27.92 \pm 0.24$ for Cn. Our results align closely with the recommended values in the 2018 Table of static dipole polarizabilities for neutral atoms [Mol. Phys. \textbf{117}, 1200 (2019)], while providing reduced uncertainties for Cd and Cn. The analysis indicates that scalar-relativistic effects are the dominant relativistic contributions to atomic dipole polarizabilities for these atoms, with SOC effects found to be negligible. Furthermore, we evaluate the influence of electron correlation across all relativistic regimes, underscoring its critical role in the precise determination of dipole polarizabilities.
- [30] arXiv:2411.05397 [pdf, other]
-
Title: Analytic model for sheath-plasma resonance in inverted fireballsComments: This original article contains 24 pages, 15 figures, etc. It has potential applications in the areas if accepted. It is for your kind consideration and needful actionSubjects: Plasma Physics (physics.plasm-ph)
The sheath plasma resonance (SPR) in an inverted fireball (IFB) system is semi-analytically investigated using a generalized hydrodynamic isothermal model formalism. It incorporates the constitutive ionic fluid viscosity, inter-species collisions, and geometric curvature effects. The SPR stability is studied for an anodic (hollow, meshed) IFB for the first time against the traditional cathode-plasma arrangements of regular electrode (solid, smooth) fireballs. The SPR develops near a spherical electrode enclosed by a plasma sheath amid a given electric potential. A generalized linear quartic dispersion relation (DR) with diverse plasma multi-parametric coefficients is methodically derived using a standard normal mode analysis. The mathematical construction of the obtained DR roots confirms that only one feasible nonzero frequency mode exists (emerging in the IFB). This root existence is confirmed both analytically and numerically. This consequent SPR creates trapped acoustic fluctuations in the IFB plasmas because of the internal reflections at the sheath plasma boundary. Also, sensible parametric changes in the SPR features, with both plasma density and viscosity, are seen. A local condition for the SPR excitation and its subsequent transition to collective standing wave-like patterns in the IFBs is illustrated. A fair corroboration of our results with the earlier SPR experimental observations of standing wave-like eigenmode patterns (evanescent) strengthens the reliability of our study alongside new applicability.
- [31] arXiv:2411.05414 [pdf, other]
-
Title: Antenne metasurface {\`a} polarisation circulaireAlejandro Arroyo (SU), Massimiliano Casaletti (SU), Romain Contreres, Alexandre Piche, Hélène RousselComments: in French language, JNM, Scienceconf, Jun 2022, Limoges (Universit{é}), FranceSubjects: Applied Physics (physics.app-ph)
A new approach using scalar metasurfaces for the design of linearly polarized antennas is presented. The proposed method is based on the construction of the surface impedance Zs using a technique called "phase-matching," which employs the sum of two circular polarizations in phase opposition. This process allows for the achievement of good performance of the synthesized antenna, such as the reduction of side lobe levels and the attainment of an almost symmetric main lobe regardless of the pointing direction. Numerical and measurement results are also presented.
- [32] arXiv:2411.05426 [pdf, other]
-
Title: Air entrainment by turbulent plunging jets: effect of jet roughness revisitedSubjects: Fluid Dynamics (physics.flu-dyn)
The amount of air entrained by vertical water jets impacting a large pool is revisited. To test available phenomenological models, new data on the jet deformation at impact and on the entrained air flow rate were collected both on a small-scale (height of fall H about 1 m, jet diameter D0 = 7.6 mm) and a large-scale (H up to 9 m, D0 up to 213 mm) facilities. Conditions for which jet break-up occurred were not considered. For short heights of fall (H less than a few D0), the jet deformation remains smaller than 0.1 jet diameter, and the entrained air flow rate happens to grow as Ui^3/2, where Ui is the jet velocity at impact. This scaling agrees with the air film model proposed by Sene, 1988. At larger fall heights, even though conditions leading to jet break-up were avoided, the jets exhibited complex topologies, including strong deformations and/or interface stripping and/or jet aeration. Further, the roughness model initiated by Henderson, McCarthy and Molloy, 1970 which stipulates that the entrained air flow rate corresponds to the air trapped within jet corrugations, was found valid for these conditions. More precisely, for corrugated jets, the effective roughness amounts to the maximum jet deformation (as measured from the 90% detection probability on the diameter pdf) or equivalently to about two times the total deformation of one side of the jet (where the total deformation of one side of the jet is experimentally evaluated as the standard deviation of the position of one jet edge). However, for jets experiencing strong stripping or aeration (the latter being identified by a threshold on the growth of the jet diameter with the falling distance), the effective roughness amounts to about 0.8 times the maximum jet deformation or equivalently to 1.1 times the total deformation of one side of the jet. Compared with corrugated jets, the effective roughness is thus diminished by half.
- [33] arXiv:2411.05427 [pdf, other]
-
Title: Interpreting Purcell Enhancement of Non-Hermitian Metasurfaces with Spectral ParametersComments: 21 pages, 7 figures, 1 tableSubjects: Optics (physics.optics)
The Purcell effect describes the enhancement of the spontaneous emission rate of an emitter near a resonant structure. However, evaluating the Purcell factor quantitatively and empirically is difficult due to the difficulties in measuring the electromagnetic nearfield of an optical resonance for calculation of the exact effective modal volume, especially with non-Hermitian resonators. Therefore, we propose a new analytical approach to circumvent the need to measure the nearfield and predict the Purcell enhancement with an analytical model and spectrally measurable parameters. Our proposed model predicts the averaged Purcell enhancement by metasurfaces on a photoluminescence medium, and is verified with experimental measurements and numerical simulations of nanoparticle arrays coupled to a fluorescent thin film. The model directly analyzes the photoluminescence enhancement and extraction efficiency of metasurface, and can be generalized to work with arbitrarily-shaped photoluminescent medium that is coupled to a resonator. This discovery will facilitate optimization of metasurfaces that support open cavity modes for efficient extraction of enhanced luminescence.
- [34] arXiv:2411.05428 [pdf, html, other]
-
Title: Exact Non-Local Hydrodynamics Predict Rarefaction EffectsSubjects: Fluid Dynamics (physics.flu-dyn); Mathematical Physics (math-ph)
We combine the theory of slow spectral closure for linearized Boltzmann equations with Maxwell's kinetic boundary conditions to derive non-local hydrodynamics with arbitrary accommodation. Focusing on shear-mode dynamics, we obtain explicit steady state solutions in terms of Fourier integrals and closed-form expressions for the mean flow and the stress. We demonstrate that the exact non-local fluid model correctly predicts several rarefaction effects with accommodation, including the Couette flow and thermal creep in a plane channel.
- [35] arXiv:2411.05455 [pdf, other]
-
Title: Cities Reconceptualized: Unveiling Hidden Uniform Urban Shape through Commute Flow Modeling in Major US CitiesSubjects: Physics and Society (physics.soc-ph); Applications (stat.AP)
Urban development is shaped by historical, geographical, and economic factors, presenting challenges for planners in understanding urban form. This study models commute flows across multiple U.S. cities, uncovering consistent patterns in urban population distributions and commuting behaviors. By embedding urban locations to reflect mobility networks, we observe that population distributions across redefined urban spaces tend to approximate log-normal distributions, in contrast to the often irregular distributions found in geographical space. This divergence suggests that natural and historical constraints shape spatial population patterns, while, under ideal conditions, urban organization may naturally align with log-normal distribution. A theoretical model using preferential attachment and random walks supports the emergence of this distribution in urban settings. These findings reveal a fundamental organizing principle in urban systems that, while not always visible geographically, consistently governs population flows and distributions. This insight into the underlying urban structure can inform planners seeking to design efficient, resilient cities.
- [36] arXiv:2411.05459 [pdf, html, other]
-
Title: Multidimensional quantum dynamics with explicitly correlated Gaussian wave packets using Rothe's methodComments: 15 pages, 11 figuresSubjects: Chemical Physics (physics.chem-ph); Computational Physics (physics.comp-ph); Quantum Physics (quant-ph)
In a previous publication [J. Chem. Phys., 161, 044105 (2024)], it has been shown that Rothe's method can be used to solve the time-dependent Schrödinger equation (TDSE) for the hydrogen atom in a strong laser field using time-dependent Gaussian wave packets. Here, we generalize these results, showing that Rothe's method can propagate arbitrary numbers of thawed, complex-valued, explicitly correlated Gaussian functions (ECGs) with dense correlation matrices for systems with varying dimensionality. We consider the multidimensional Henon-Heiles potential, and show that the dynamics can be quantitatively reproduced using only 30 Gaussians in 2D, and that accurate spectra can be obtained using 20 Gaussians in 2D and 30 to 40 Gaussians in 3D and 4D. Thus, the relevant multidimensional dynamics can be described at high quality using only a small number of ECGs that give a very compact representation of the wave function. This efficient representation, along with the demonstrated ability of Rothe's method to propagate Gaussian wave packets in strong fields and ECGs in complex potentials, paves the way for accurate molecular dynamics calculations beyond the Born-Oppenheimer approximation in strong fields.
- [37] arXiv:2411.05475 [pdf, html, other]
-
Title: 3D-Printed Dual-Polarized Magneto-Electric Dipole Antenna with Wideband High Isolation for Full-Duplex ApplicationsComments: This work has been submitted to the IEEE for possible publicationSubjects: Applied Physics (physics.app-ph); Systems and Control (eess.SY)
The paper introduces a novel dual-port dual-polarized magneto-electric dipole (MED) antenna with orthogonal Gamma and inverted-Gamma shape probes, which was fabricated by means of an additive 3D metal printing process. Electromagnetic wave simulation and RF measurement report a resonance bandwidth from 3 GHz to 4 GHz at both MED's ports with respect to a standing wave ratio of less than 2. The cross-polarization isolation (XPI) between the MED's ports was also measured to be greater than 50 dB across its entire resonance bandwidth. The paper also thoroughly examines the impact of misalignments in the polarization of the MED probes on the XPI level. The broadband resonance and excellent isolation between the MED ports make it a strong candidate for a full-duplex wireless transceiver in network infrastructure.
- [38] arXiv:2411.05477 [pdf, html, other]
-
Title: Generating surrogate temporal networks from mesoscale building blocksSubjects: Physics and Society (physics.soc-ph); Data Analysis, Statistics and Probability (physics.data-an)
Surrogate networks can constitute suitable replacements for real networks, in particular to study dynamical processes on networks, when only incomplete or limited datasets are available. As empirical datasets most often present complex features and interplays between structure and temporal evolution, creating surrogate data is however a challenging task, in particular for data describing time-resolved interactions between agents. Here we propose a method to generate surrogate temporal networks that mimic such observed datasets. The method is based on a decomposition of the original dataset into small temporal subnetworks encoding local structures on a short time scale. These are used as building blocks to generate a new synthetic temporal network that will hence inherit the shape of local interactions from the dataset. Moreover, we also take into account larger scale correlations on structural and temporal dimension, using them to inform the process of assembling the building blocks. We showcase the method by generating surrogate networks for several datasets of social interactions and comparing them to the original data on two complementary aspects. First, we show that the surrogate data possess complex structural and temporal features similar to the ones of the original data. Second, we simulate several dynamical processes, describing respectively epidemic spread, opinion formation and emergence of norms in a population, and compare the outcome of these processes on the generated and original datasets. We describe the method in detail and provide an implementation so that it can be easily used in future works based on temporally evolving networks.
- [39] arXiv:2411.05496 [pdf, html, other]
-
Title: Combining Bayesian Optimization, SVD and Machine Learning for Advanced Optical DesignSubjects: Optics (physics.optics); Applied Physics (physics.app-ph); Computational Physics (physics.comp-ph)
The design and optimization of optical components, such as Bragg gratings, are critical for applications in telecommunications, sensing, and photonic circuits. To overcome the limitations of traditional design methods that rely heavily on computationally intensive simulations and large datasets, we propose an integrated methodology that significantly reduces these burdens while maintaining high accuracy in predicting optical response. First, we employ a Bayesian optimization technique to strategically select a limited yet informative number of simulation points from the design space, ensuring that each contributes maximally to the model's performance. Second, we utilize singular value decomposition to effectively parametrize the entire reflectance spectra into a reduced set of coefficients, allowing us to capture all significant spectral features without losing crucial information. Finally, we apply XGBoost, a robust machine learning algorithm, to predict the entire reflectance spectra from the reduced dataset. The combination of Bayesian optimization for data selection, SVD for full-spectrum fitting, and XGBoost for predictive modeling provides a powerful, generalizable framework for the design of optical components.
- [40] arXiv:2411.05499 [pdf, html, other]
-
Title: On nonlinear transitions, minimal seeds and exact solutions for the geodynamoComments: 21 pages, 8 figuresSubjects: Fluid Dynamics (physics.flu-dyn); Earth and Planetary Astrophysics (astro-ph.EP); Geophysics (physics.geo-ph)
Nearly fifty years ago, Roberts (1978) postulated that Earth's magnetic field, which is generated by turbulent motions of liquid metal in its outer core, likely results from a subcritical dynamo instability characterised by a dominant balance between Coriolis, pressure and Lorentz forces. Here we numerically explore the generation of subcritical geomagnetic fields using techniques from optimal control and dynamical systems theory to uncover the nonlinear dynamical landscape underlying dynamo action. Through nonlinear optimisation, via direct-adjoint looping, we identify the minimal seed - the smallest magnetic field that attracts to a nonlinear dynamo solution. Additionally, using the Newton-hookstep algorithm, we converge stable and unstable travelling wave solutions to the governing equations. By combining these two techniques, complex nonlinear pathways between attracting states are revealed, providing insight into a potential subcritical origin of the geodynamo. This paper showcases these methods on the widely studied benchmark of Christensen et al. (2001), laying the foundations for future studies in more extreme and realistic parameter regimes. We show that the minimal seed reaches a nonlinear dynamo solution by first attracting to an unstable travelling wave solution, which acts as an edge state separating a hydrodynamic solution from a magnetohydrodynamic one. Furthermore, by carefully examining the choice of cost functional, we establish a robust optimisation procedure that can systematically locate dynamo solutions on short time horizons with no prior knowledge of its structure.
- [41] arXiv:2411.05519 [pdf, html, other]
-
Title: Surface-acoustic-wave driven silicon microfluidic chips for acoustic tweezing of motile cells and viscoelastic microbeadsSubjects: Fluid Dynamics (physics.flu-dyn)
Acoustic tweezers comprising a surface acoustic wave chip and a disposable silicon microfluidic chip are potentially advantageous to stable and cost-ffective acoustofluidic experiments while avoiding the cross-contamination by reusing the surface acoustic wave chip and disposing of the microfluidic chip. For such a device, it is important to optimize the chip-to-chip bonding and the size and shape of the microfluidic chip to enhance the available acoustic pressure. In this work, aiming at studying samples with the size of a few tens of microns, we explore the device structure and assembly method of acoustic tweezers. By using a polymer bonding layer and shaping the silicon microfluidic chip via deep reactive ion etching, we were able to attain the acoustic pressure up to 2 MPa with a corresponding acoustic radiation pressure of 0.2 kPa for 50 MHz ultrasound, comparable to reported values at lower ultrasound frequencies. We utilized the fabricated acoustic tweezers for non-contact viscoelastic deformation experiments of soft matter and trapping of highly motile cells. These results suggests that the feasibility of the hybrid chip approach to attaining the high acoustic force required to conduct acoustomechanical testing of small soft matters and cells.
- [42] arXiv:2411.05525 [pdf, other]
-
Title: The influence of geometry and specific electronic and nuclear energy deposition on ion-stimulated desorption from thin self-supporting membranesSubjects: Chemical Physics (physics.chem-ph); Applied Physics (physics.app-ph)
We investigate the dependence of the yield of positive secondary ions created upon impact of primary He, B and Ne ions on geometry and electronic and nuclear energy deposition by the projectiles. We employ pulsed beams in the medium energy regime and a large position-sensitive, time-of-flight detection system to ensure accurate quantification. As a target, we employ a single crystalline Si(100) self-supporting 50 nm thick membrane thus featuring two identical surfaces enabling simultaneous measurements in backscattering and transmission geometry. Electronic sputtering is identified as the governing mechanism for the desorption of hydrogen and molecular species found on the surfaces. Nevertheless, larger energy deposition to the nuclear subsystem by heavier projectiles as well as due to the directionality of the collision cascade appears to act in synergy with the electronic energy deposition leading to an overall increase in secondary ion yields. A higher yield of ions sputtered from the matrix is observed in transmission geometry only for B and Ne ions, consistent with the observed role of nuclear stopping.
- [43] arXiv:2411.05528 [pdf, html, other]
-
Title: Prediction of Mode Structure Using A Novel Physics-Embedded Neural ODE MethodComments: 21 pages, 3 figures, and 4 tables. Prepared for submission to Physics of PlasmasSubjects: Plasma Physics (physics.plasm-ph)
We designed a new artificial neural network by modifying the neural ordinary differential equation (NODE) framework to successfully predict the time evolution of the 2D mode profile in both the linear growth and nonlinear saturated stages. Starting from the magnetohydrodynamic (MHD) equations, simplifying assumptions were applied based on physical properties and symmetry considerations of the energetic-particle-driven geodesic acoustic mode (EGAM) to reduce complexity. Our approach embeds physical laws directly into the neural network architecture by exposing latent differential states, enabling the model to capture complex features in the nonlinear saturated stage that are difficult to describe analytically, and thus, the new artificial neural network is named as ExpNODE (Exposed latent state Neural ODE). ExpNODE was evaluated using a data set generated from first-principles simulations of the EGAM instability, focusing on the pre-saturated stage and the nonlinear saturated stage where the mode properties are most complex. Compared to state-of-the-art models such as ConvLSTM, ExpNODE with physical information not only achieved lower test loss but also converged faster during training. Specifically, it outperformed ConvLSTM method in both the 20-step and 40-step prediction horizons, demonstrating superior accuracy and efficiency. Additionally, the model exhibited strong generalization capabilities, accurately predicting mode profiles outside the training data set. Visual comparisons between model predictions and ground truth data showed that ExpNODE with physical information closely captured detailed features and asymmetries inherent in the EGAM dynamics that were not adequately captured by other models. These results suggest that integrating physical knowledge into neural ODE frameworks enhances their performance, and provides a powerful tool for modeling complex plasma phenomena.
- [44] arXiv:2411.05530 [pdf, other]
-
Title: Digitization of 93.3-keV recoilless $\gamma$-radiationComments: arXiv admin note: text overlap with arXiv:2410.16157Subjects: Optics (physics.optics)
We propose a technique for transforming the intensity of quasi-monochromatic recoilless radiation with a photon energy of 93.3 keV, emitted by radioactive Mössbauer sources 67Ga or 67Cu, into a sequence of short pulses with individually and independently controlled, on demand, the moments of appearance of pulses, as well as the peak intensity, duration and shape of each pulse. The technique is based on the transmission of recoilless (Mössbauer) 93.3-keV photons from the source through a medium containing resonantly absorbing 67Zn nuclei. The pulses are formed due to the rapid reciprocating movement of the source relative to the absorber at specified moments in time along the direction of photon propagation at a distance not exceeding half the radiation wavelength. The produced sequences of $\gamma$-ray pulses are similar to the digitization of information carried by electromagnetic waves. They can also be used to develop Mössbauer spectroscopy of the atomic and subatomic structures, as well as open novel opportunities for $\gamma$-ray quantum optics.
- [45] arXiv:2411.05535 [pdf, html, other]
-
Title: Dynamic equilibrium of electrochemical bubbles growing on micro-electrodesComments: 16 pages, 13 figuresSubjects: Fluid Dynamics (physics.flu-dyn)
In gas evolving electrolysis, bubbles grow at electrodes due to a diffusive influx from oversaturation generated locally in the electrolyte by the electrode reaction. When considering electrodes of micrometer-size resembling catalytic islands, bubbles are found to approach dynamic equilibrium states at which they neither grow nor shrink. Such equilibrium states are found at low oversaturation for both, pinning and expanding wetting regimes of the bubbles and are based on the balance of local influx near the bubble foot and global outflux. Unlike the stability of pinned nano-bubbles studied earlier, the Laplace pressure plays a minor role only. We extend the analytical solution of Zhang & Lohse (2023) by taking into account the non-uniform distribution of dissolved gas around the bubble obtained from direct numerical simulation. This allows us to identify the parameter regions of bubble growth, dissolution and dynamic equilibrium as well as to analyze the stability of the equilibrium states. Finally, we draw conclusions on how to possibly enhance the efficiency of electrolysis.
- [46] arXiv:2411.05546 [pdf, html, other]
-
Title: Isolated Attosecond $\gamma$-Ray Pulse Generation with Transverse Orbital Angular Momentum Using Intense Spatiotemporal Optical Vortex LasersSubjects: Plasma Physics (physics.plasm-ph)
An isolated attosecond vortex $\gamma$-ray pulse is generated by using a relativistic spatiotemporal optical vortex (STOV) laser in particle-in-cell simulations. A $\sim$ 300-attosecond electron slice with transverse orbital angular momentum (TOAM) is initially selected and accelerated by the central spatiotemporal singularity of the STOV laser. This slice then collides with the laser's reflected Gaussian-like front from a planar target, initiating nonlinear Compton scattering and resulting in an isolated, attosecond ($\sim$ 300 as), highly collimated ($\sim$ 4$\degree$), ultra-brilliant ($\sim 5\times 10^{24}$ photons/s/mm$^2$/mrad$^2$/0.1\%BW at 1 MeV) $\gamma$-ray pulse. This STOV-driven approach overcomes the significant beam divergence and complex two-laser requirements of prior Gaussian-based methods while introducting TOAM to the attosecond $\gamma$-ray pulse, which opens avenues for ultrafast imaging, nuclear excitation, and detection applications.
- [47] arXiv:2411.05576 [pdf, other]
-
Title: On-chip Moir\'e Optical Skyrmion Clusters with Nanoscale DynamicsSubjects: Optics (physics.optics); Applied Physics (physics.app-ph)
Skyrmions are topological defects belonging to nontrivial homotopy classes in particle theory, and are recognized as a suitable unit in the high-density, low-dissipation microelectronic devices in condensed matter physics. Their remarkably stable topology has been observed in electromagnetic waves recently. For the evanescent fields near a surface, this has been realized so far only for elementary optical skyrmions, with a fixed skyrmion number. Here we introduce the concept of moiré optical skyrmion clusters-multiskyrmions are nested to form a large optical skyrmion cluster-crystallized or quasi-crystallized as a consequence of the twisted nanostructures. The rapid inverting of optical skyrmion number is achieved in the imperfectly aligned composite nanostructures. This moiré optical skyrmion interaction mechanism is described by a lattice model. Further, the nucleation and collapse of optical skyrmion are studied, where their nanoscale dynamics are revealed with a tiny change of the twist angle. The sudden reversal of the on-chip skyrmion can serve as a precise beacon of the relative alignment deviation between twisted composite nanostructures.
- [48] arXiv:2411.05579 [pdf, html, other]
-
Title: Electron vortex beams for chirality probing at the nanoscaleJournal-ref: Phys. Rev. Applied 22 (2024) 054017Subjects: Optics (physics.optics)
In this work we propose a method for probing the chirality of nanoscale electromagnetic near fields utilizing the properties of a coherent superposition of free-electron vortex states in electron microscopes. Electron beams optically modulated into vortices carry orbital angular momentum, thanks to which they are sensitive to the spatial phase distribution and topology of the investigated field. The sense of chirality of the studied specimen can be extracted from the spectra of the electron beam with nanoscale precision owing to the short picometer de Broglie wavelength of the electron beam. We present a detailed case study of the interaction of a coherent superposition of electron vortex states and the optical near field of a golden nanosphere illuminated by circularly polarized light as an example, and we examine the chirality sensitivity of electron vortex beams on intrinsically chiral plasmonic nanoantennae.
- [49] arXiv:2411.05585 [pdf, html, other]
-
Title: Long-range and dead-zone free dual-comb ranging for the interferometric tracking of moving targetsSandro L. Camenzind, Lukas Lang, Benjamin Willenberg, Justinas Pupeikis, Hayk Soghomonyan, Robert Presl, Pabitro Ray, Andreas Wieser, Ursula Keller, Christopher R. PhillipsComments: Supplementary Information includedSubjects: Optics (physics.optics)
Dual-comb ranging has emerged as an effective technology for long-distance metrology, providing absolute distance measurements with high speed, precision, and accuracy. Here, we demonstrate a dual-comb ranging method that utilizes a free-space transceiver unit, enabling dead-zone-free measurements and simultaneous ranging with interchanged comb roles to allow for long-distance measurements even when the target is moving. It includes a GPU-accelerated algorithm for real-time signal processing and a free-running single-cavity solid-state dual-comb laser with a carrier wavelength $\lambda_c \approx$ 1055 nm, a pulse repetition rate of 1 GHz and a repetition rate difference of 5.06 kHz. This combination offers a fast update rate and sufficient signal strength to reach a single-shot time-of-flight precision of around 0.1 $\mu$m (i.e. $< \lambda_c/4$) on a cooperative target placed at a distance of more than 40 m. The free-running laser is sufficiently stable to use the phase information for interferometric distance measurements, which improves the single-shot precision to $<$20 nm. To assess the ranging accuracy, we track the motion of the cooperative target when moved over 40 m and compare it to a reference interferometer. The residuals between the two measurements are below 3 $\mu$m. These results highlight the potential of this approach for accurate and dead-zone-free long-distance ranging, supporting real-time tracking with nm-level precision.
- [50] arXiv:2411.05589 [pdf, other]
-
Title: Predicting Resistive Pulse Signatures in Nanopores by Accurately Modeling Access RegionsSubjects: Applied Physics (physics.app-ph); Biological Physics (physics.bio-ph)
Resistive pulse sensing has been widely used to characterize and count single particles in solution moving through channels under an electric bias, with nanoscale pores more recently providing enough spatial resolution for nucleic acid sequencing at the single-molecule level. At its core, this technique relies on measuring the drop in ionic current through the pore induced by the passage of a molecule and, through conductance models, translating the blockage signal to molecular dimensions. However, there exists no model considering the resistive contributions of the pore exterior, i.e. the access regions, when obstructed by a molecule. This is becoming increasingly important for low aspect ratio pores, with the advent of 2D materials and ultrathin membranes. In this work, a general method by which to model the resistance of the access regions of a pore in the presence of an insulating obstruction is presented. Thin oblate spheroidal slices are used to partition access regions and infer their conductance when blocked by differently shaped objects. We show that our model accurately estimates the blocked-state conductance of 2D and finite-length pores as a function of the distance from the pore in the presence of simple obstructions geometries (e.g. cylindrical and spherical objects) or complex structures (i.e. sequence of simple obstruction sub-units). The model is further shown to capture off-axis effects by predicting deeper blockages for obstructions offset from the pore's central axis. A web-based tool is created to predict the electrical signatures of a wide range of molecule geometries translocating through differently shaped pores. The introduced model will help guide experimental designs and thus presents a straightforward way to extend the quantification of the resistive pulse technique at the nanoscale.
- [51] arXiv:2411.05592 [pdf, html, other]
-
Title: Turbulent Boundary Layer in a 3-Element High-LiftWing: Coherent Structures IdentificationComments: DLES13 2022 conference proceeding paperSubjects: Fluid Dynamics (physics.flu-dyn)
A wall-resolved large-eddy simulation (LES) of the fluid flow around a 30P30N airfoil is conducted at a Reynolds number of Rec=750,000 and an angle of attack (AoA) of 9 degrees. The simulation results are validated against experimental data from previous studies and further analyzed, focusing on the suction side of the wing main element. The boundary layer development is investigated, showing characteristics typical of a zero-pressure-gradient turbulent boundary layer (ZPG TBL). In particular, the boundary layer exhibits limited growth, and the outer peak of the streamwise Reynolds stresses is virtually absent, distinguishing it from an adverse-pressure-gradient turbulent boundary layer (APG TBL). A proper orthogonal decomposition (POD) analysis is performed on a portion of the turbulent boundary layer, revealing a significant energy spread across higher-order modes. Despite this, TBL streaks are identified, and the locations of the most energetic structures correspond to the peaks in the Reynolds stresses.
- [52] arXiv:2411.05610 [pdf, html, other]
-
Title: Creating compact localized modes for robust sound transport via singular flat band engineeringSubjects: Applied Physics (physics.app-ph)
We experimentally demonstrate the emergence of flat-band-induced compact-localized modes in acoustic Kagome lattices. Compact localized states populate singular dispersion bands characterized by band crossing, where a quadratic and a flat-band dispersion coalesce into a singularity. These conditions enable intriguing wave phenomena when the Hilbert Schmidt quantum distance, measuring the strength of the singularity, is nonzero. We report numerically and experimentally the formation of compact localized states (CLS), extremely localized in space and protected by dispersion flatness. In our system of coupled acoustic waveguides, sound waves are confined to propagate within tightly localized sites positioned both at the boundaries and within the interior of the lattice, achieving broadband and sustained confinement over time. This framework opens new avenues for the manipulation and transport of information through sound waves, with potential application in mechanics and acoustics, including communication, signal processing, and sound isolation. This work also expands the exploration of flat-band lattice physics within the realm of acoustics.
- [53] arXiv:2411.05615 [pdf, html, other]
-
Title: New methods of neutrino and anti-neutrino detection from 0.115 to 105 MeVNickolas Solomey, Mark Christl, Brian Doty, Jonathan Folkerts, Brooks Hartsock, Evgen Kuznetsco, Robert McTaggart, Holger Meyer, Tyler Nolan, Greg Pawloski, Daniel Reichart, Miguel Rodriguez-Otero, Dan Smith, Lisa SolomeyComments: Contribution to the 25th International Workshop on Neutrinos from AcceleratorsSubjects: Instrumentation and Detectors (physics.ins-det); High Energy Physics - Experiment (hep-ex); Nuclear Experiment (nucl-ex)
We have developed a neutrino detector with threshold energies from ~0.115 to 105 MeV in a clean detection mode almost completely void of accidental backgrounds. It was initially developed for the NASA $\nu$SOL project to put a solar neutrino detector very close to the Sun with 1,000 to 10,000 times higher solar neutrino flux than on Earth. Similar interactions have been found for anti-neutrinos, which were initially intended for Beta decay neutrinos from reactors, geological sources, or for nuclear security applications. These techniques work at the 1 to 100 MeV region for neutrinos from the ORNL Spallation Neutron Source or low energy accelerator neutrino and anti-neutrino production targets less than $\sim$100 MeV. The identification process is clean, with a double pulse detection signature within a time window between the first interaction producing the conversion electron or positron and the secondary gamma emission 100 ns to ~1 $\mu$s, which removes most accidental backgrounds. These new modes for neutrino and anti-neutrino detection of low energy neutrinos and anti-neutrinos could allow improvements to neutrino interaction measurements from an accelerator beam on a target.
- [54] arXiv:2411.05617 [pdf, html, other]
-
Title: Identifying spatially-localized instability mechanisms using sparse optimizationComments: 9 pages, 2 figuresSubjects: Fluid Dynamics (physics.flu-dyn); Mathematical Physics (math-ph)
Recent investigations have established the physical relevance of spatially-localized instability mechanisms in fluid dynamics and their potential for technological innovations in flow control. In this letter, we show that the mathematical problem of identifying spatially-localized optimal perturbations that maximize perturbation-energy amplification can be cast as a sparse (cardinality-constrained) optimization problem. Unfortunately, cardinality constrained optimization problems are non-convex and combinatorially hard to solve in general. To make the analysis viable within the context of fluid dynamics problems, we propose an efficient iterative method for computing sub-optimal spatially-localized perturbations. Our approach is based on a generalized Rayleigh quotient iteration algorithm followed by a variational renormalization procedure that reduces the optimality gap in the resulting solution. The approach is demonstrated on a sub-critical plane Poiseuille flow at Re = 4000, which has been a benchmark problem studied in prior investigations on identifying spatially-localized flow structures. Remarkably, we find that a subset of the perturbations identified by our method yield a comparable degree of energy amplification as their global counterparts. We anticipate our proposed analysis tools will facilitate further investigations into spatially-localized flow instabilities, including within the resolvent and input-output analysis frameworks.
- [55] arXiv:2411.05621 [pdf, html, other]
-
Title: Second-order hyperfine correction to H, D, and $^3$He energy levelsComments: 11 pagesSubjects: Atomic Physics (physics.atom-ph)
The complete second-order hyperfine-interaction correction is calculated for centroid energy levels of H, D, and $^3$He atoms. For $^3$He, the corrections of $-2.075$ kHz and $-0.305$ kHz beyond the leading hyperfine-mixing contribution are obtained for the $2^1S$ and $2^3S$ states, respectively. These results shift the nuclear charge radii difference derived from the $^3$He - $^4$He isotope shift and largely resolve the previously reported disagreement between the muonic and electronic helium determinations [Y. van der Werf et al., arXiv:2306.02333 (2023); K. Schuhmann et al., arXiv:2305.11679 (2023)].
- [56] arXiv:2411.05626 [pdf, html, other]
-
Title: Can Efficient Fourier-Transform Techniques Favorably Impact on Broadband Computational Electromagnetism?Subjects: Computational Physics (physics.comp-ph); Numerical Analysis (math.NA)
In view of recently demonstrated joint use of novel Fourier-transform techniques and effective high-accuracy frequency domain solvers related to the Method of Moments, it is argued that a set of transformative innovations could be developed for the effective, accurate and efficient simulation of problems of wave propagation and scattering of broadband, time-dependent wavefields. This contribution aims to convey the character of these methods and to highlight their applicability in computational modeling of electromagnetic configurations across various fields of science and engineering.
- [57] arXiv:2411.05630 [pdf, html, other]
-
Title: Collisional charging of dust particles by suprathermal plasmas. II -- Regularized Kappa distributionsComments: Submitted for publication in Physics of PlasmasSubjects: Plasma Physics (physics.plasm-ph); Solar and Stellar Astrophysics (astro-ph.SR)
We study the effects of the velocity distribution functions of the plasma particles on the equilibrium charge of dust grains, acquired through inelastic collisions of the particles with the grains. This paper is the second in a series of two papers on the subject. Here, we consider the charging process when the plasma particles are statistically described by the recently proposed regularized Kappa distribution functions, which allow for extreme suprathermal states, characterized by extremely low values of the kappa index, previously forbidden to the standard Kappa distributions, whose effects on dust charging were studied in Paper I of this series. We analyse the effects that extreme suprathermal states of the plasma particles have on dust charging and verify conditions for the uncommon result of positive equilibrium charge, employing two different models for the regularized Kappa distributions.
- [58] arXiv:2411.05631 [pdf, html, other]
-
Title: Physics-constrained coupled neural differential equations for one dimensional blood flow modelingSubjects: Fluid Dynamics (physics.flu-dyn); Machine Learning (cs.LG)
Computational cardiovascular flow modeling plays a crucial role in understanding blood flow dynamics. While 3D models provide acute details, they are computationally expensive, especially with fluid-structure interaction (FSI) simulations. 1D models offer a computationally efficient alternative, by simplifying the 3D Navier-Stokes equations through axisymmetric flow assumption and cross-sectional averaging. However, traditional 1D models based on finite element methods (FEM) often lack accuracy compared to 3D averaged solutions. This study introduces a novel physics-constrained machine learning technique that enhances the accuracy of 1D blood flow models while maintaining computational efficiency. Our approach, utilizing a physics-constrained coupled neural differential equation (PCNDE) framework, demonstrates superior performance compared to conventional FEM-based 1D models across a wide range of inlet boundary condition waveforms and stenosis blockage ratios. A key innovation lies in the spatial formulation of the momentum conservation equation, departing from the traditional temporal approach and capitalizing on the inherent temporal periodicity of blood flow. This spatial neural differential equation formulation switches space and time and overcomes issues related to coupling stability and smoothness, while simplifying boundary condition implementation. The model accurately captures flow rate, area, and pressure variations for unseen waveforms and geometries. We evaluate the model's robustness to input noise and explore the loss landscapes associated with the inclusion of different physics terms. This advanced 1D modeling technique offers promising potential for rapid cardiovascular simulations, achieving computational efficiency and accuracy. By combining the strengths of physics-based and data-driven modeling, this approach enables fast and accurate cardiovascular simulations.
- [59] arXiv:2411.05640 [pdf, html, other]
-
Title: Adaptive Dissipation in the Smagorinsky Model for Turbulence in Boundary-Driven FlowsComments: 9 pagesSubjects: Fluid Dynamics (physics.flu-dyn)
This paper enhances the classic Smagorinsky model by introducing an innovative, adaptive dissipation term that adjusts dynamically with distance from boundary regions. This modification addresses a known limitation of the standard model over dissipation near boundaries thereby improving accuracy in turbulent flow simulations in confined or wall-adjacent areas. We present a rigorous theoretical framework for this adaptive model, including two foundational theorems. The first theorem guarantees existence and uniqueness of solutions, ensuring that the model is mathematically well-posed within the adaptive context. The second theorem provides a precise bound on the energy dissipation rate, demonstrating that dissipation remains controlled and realistic even as boundary effects vary spatially. By allowing the dissipation coefficient to decrease near boundary layers, this approach preserves the finer turbulent structures without excessive smoothing, yielding a more physically accurate representation of the flow. Future work will focus on implementing this adaptive model in computational simulations to empirically verify the theoretical predictions and assess performance in scenarios with complex boundary geometries.
- [60] arXiv:2411.05644 [pdf, other]
-
Title: Polarization-independent metasurfaces based on bound states in the continuum with high Q-factor and resonance modulationXingye Yang, Alexander Antonov, Andreas Aigner, Thomas Weber, Yohan Lee, Tao Jiang, Haiyang Hu, Andreas TittlSubjects: Optics (physics.optics); Applied Physics (physics.app-ph)
Metasurfaces offer a powerful platform for effective light manipulation, which is crucial for advanced optical technologies. While designs of polarization-independent structures have reduced the need for polarized illumination, they are often limited by either low Q factors or low resonance modulation. Here, we design and experimentally demonstrate a metasurface with polarization-independent quasi-bound state in the continuum (quasi-BIC), where the unit cell consists of four silicon squares arranged in a two-dimensional array and the resonance properties can be controlled by adjusting the edge length difference between different squares. Our metasurface experimentally achieves a Q factor of approximately 100 and a resonance modulation of around 50%. This work addresses a common limitation in previous designs, which either achieved high Q factors exceeding 200 with a resonance modulation of less than 10%, leading to challenging signal-to-noise ratio requirements, or achieved strong resonance modulation with Q factors of only around 10, limiting light confinement and fine-tuning capabilities. In contrast, our metasurface ensures that the polarization-independent signal is sharp and distinct within the system, reducing the demands on signal-to-noise ratio and improving robustness. Experiments show the consistent performance across different polarization angles. This work contributes to the development of versatile optical devices, enhancing the potential for the practical application of BIC-based designs in areas such as optical filtering and sensing.
- [61] arXiv:2411.05647 [pdf, html, other]
-
Title: Telecom wavelength quantum dots interfaced with silicon-nitride circuits via photonic wire bondingUlrich Pfister, Daniel Wendland, Florian Hornung, Lena Engel, Hendrik Hüging, Elias Herzog, Ponraj Vijayan, Raphael Joos, Erik Jung, Michael Jetter, Simone L. Portalupi, Wolfram H. P. Pernice, Peter MichlerSubjects: Optics (physics.optics); Quantum Physics (quant-ph)
Photonic integrated circuits find ubiquitous use in various technologies, from communication, to computing and sensing, and therefore play a crucial role in the quantum technology counterparts. Several systems are currently under investigation, each showing distinct advantages and drawbacks. For this reason, efforts are made to effectively combine different platforms in order to benefit from their respective strengths. In this work, 3D laser written photonic wire bonds are employed to interface triggered sources of quantum light, based on semiconductor quantum dots embedded into etched microlenses, with low-loss silicon-nitride photonics. Single photons at telecom wavelengths are generated by the In(Ga)As quantum dots which are then funneled into a silicon-nitride chip containing single-mode waveguides and beamsplitters. The second-order correlation function of g(2)(0) = 0.11+/-0.02, measured via the on-chip beamsplitter, clearly demonstrates the transfer of single photons into the silicon-nitride platform. The photonic wire bonds funnel on average 28.6+/-8.8% of the bare microlens emission (NA = 0.6) into the silicon-nitride-based photonic integrated circuit even at cryogenic temperatures. This opens the route for the effective future up-scaling of circuitry complexity based on the use of multiple different platforms.
- [62] arXiv:2411.05657 [pdf, html, other]
-
Title: Electromagnetic Multipole Theory for Two-dimensional PhotonicsSubjects: Optics (physics.optics)
We develop a full-wave electromagnetic (EM) theory for calculating the multipole decomposition in two-dimensional (2-D) structures consisting of isolated, arbitrarily shaped, inhomogeneous, anisotropic cylinders or a collection of such. To derive the multipole decomposition, we first solve the scattering problem by expanding the scattered electric field in divergenceless cylindrical vector wave functions (CVWF) with unknown expansion coefficients that characterize the multipole response. These expansion coefficients are then expressed via contour integrals of the vectorial components of the scattered electric field evaluated via an electric field volume integral equation (EFVIE). The kernels of the EFVIE are the products of the tensorial 2-D Green's function (GF) expansion and the equivalent 2-D volumetric electric and magnetic current densities. We validate the theory using the commercial finite element solver COMSOL Multiphysics. In the validation, we compute the multipole decomposition of the fields scattered from various 2-D structures and compare the results with alternative formulations. Finally, we demonstrate the applicability of the theory to study an emerging photonics application on oligomers-based highly directional switching using active media. This analysis addresses a critical gap in current literature, where multipole theories exist primarily for three-dimensional (3-D) particles of isotropic materials. Our work enhances the understanding and utilization of the optical properties of 2-D, inhomogeneous, and anisotropic cylindrical structures, contributing to advancements in photonic and meta-optics technologies.
- [63] arXiv:2411.05660 [pdf, html, other]
-
Title: Stratification governs the Existence of Surface-Intensified Eastward Jets in Turbulent Gyres without Bottom FrictionSubjects: Fluid Dynamics (physics.flu-dyn); Atmospheric and Oceanic Physics (physics.ao-ph)
This study examines the role of stratification in the formation and persistence of eastward jets (like the Gulf Stream and Kuroshio currents). Using a wind-driven, two-layer quasi-geostrophic model in a double-gyre configuration, we construct a phase diagram to classify flow regimes. The parameter space is defined by a criticality parameter $\xi$, which controls the emergence of baroclinic instability, and the ratio of layer depths $\delta$, which describes the surface intensification of stratification. Eastward jets detaching from the western boundary are observed when $\delta \ll 1$ and $\xi \sim 1$, representing a regime transition from a vortex-dominated western boundary current to a zonostrophic regime characterized by multiple eastward jets. Remarkably, these surface-intensified patterns emerge without considering bottom friction. The emergence of the coherent eastward jet is further addressed with complementary 1.5-layer simulations and explained through both linear stability analysis and turbulence phenomenology. In particular, we show that coherent eastward jets emerge when the western boundary layer is stable, and find that the asymmetry in the baroclinic instability of eastward and westward flows plays a central role in the persistence of eastward jets, while contributing to the disintegration of westward jets.
- [64] arXiv:2411.05685 [pdf, html, other]
-
Title: Beyond Pairwise Interactions: Unveiling the Role of Higher-Order Interactions via Stepwise ReductionComments: 14 pages, 6 figuresSubjects: Physics and Society (physics.soc-ph); Social and Information Networks (cs.SI)
Complex systems, such as economic, social, biological, and ecological systems, usually feature interactions not only between pairwise entities but also among three or more entities. These multi-entity interactions are known as higher-order interactions. Hypergraph, as a mathematical tool, can effectively characterize higher-order interactions, where nodes denote entities and hyperedges represent interactions among multiple entities. Meanwhile, all higher-order interactions can also be projected into a number of lower-order interactions or even some pairwise interactions. Whether it is necessary to consider all higher-order interactions, and whether it is with little loss to replace them by lower-order or even pairwise interactions, remain a controversial issue. If the role of higher-order interactions is insignificant, the complexity of computation and the difficulty of analysis can be drastically reduced by projecting higher-order interactions into lower-order or pairwise interactions. We use link prediction, a fundamental problem in network science, as the entry point. Specifically, we evaluate the impact of higher-order interactions on link predictive accuracy to explore the necessity of these structures. We propose a method to decompose the higher-order structures in a stepwise way, thereby allowing to systematically explore the impacts of structures at different orders on link prediction. The results indicate that in some networks, incorporating higher-order interactions significantly enhances the accuracy of link prediction, while in others, the effect is insignificant. Therefore, we think that the role of higher-order interactions varies in different types of networks. Overall, since the improvement in predictive accuracy provided by higher-order interactions is significant in some networks, we believe that the study of higher-order interactions is both necessary and valuable.
- [65] arXiv:2411.05686 [pdf, html, other]
-
Title: Comprehensive Study on the Slat Noise of 30P30N High-Lift Airfoil Basd on High-Order Wall-Resolved Large-Eddy SimulationSubjects: Fluid Dynamics (physics.flu-dyn)
This study presents wall-resolved large-eddy simulations (WRLES) of a high-lift airfoil, based on high-order flux reconstruction (FR) commercial software Dimaxer, which runs on consumer level GPUs. A series of independence tests are conducted, including various Ffowcs Williams-Hawkings sampling surfaces, different mesh densities, simulations at 4th and 5th order accuracies, and varying spanwise lengths, to establish best practice for predicting slat noise through high-order WRLES. The results show excellent agreement with experimental data while requiring significantly fewer computational resources than traditional second-order methods. An investigation on the effects of Reynolds number (Re) is performed by scaling the airfoil size, with Reynolds numbers ranging from 8.55e5 to a real aircraft level of 1.71e7. By applying simple scaling through Strouhal number (St), spanwise correction, and distance from the receiver, the far-field noise spectra for different Reynolds numbers can be coincided. Additionally, simulations are performed at four angles of attack: 3°, 5.5°, 9.5°, and 14°. The results indicate that higher angles of attack lead to a less intense feedback loop, resulting in lower tonal noise frequencies and reduced noise amplitude. The maximum noise reduction observed is over 14dB when comparing 14° to 3°. Furthermore, an improved formula is proposed to enhance the prediction of slat noise tonal frequencies and to better elucidate the mechanism behind tonal noise generation.
- [66] arXiv:2411.05687 [pdf, html, other]
-
Title: Quadratic Zeeman and Electric Quadrupole Shifts in Highly Charged IonsComments: 14 pages, 1 figureSubjects: Atomic Physics (physics.atom-ph)
Recent advances in high-precision spectroscopy of highly charged ions necessitate an understanding of energy shifts of ionic levels caused by external electric and magnetic fields. Beyond the well-known Stark and linear Zeeman shifts, trapped ions may also exhibit quadratic Zeeman and electric quadrupole shifts. In this contribution, we present a systematic approach for the theoretical analysis of these shifts for arbitrary many-electron ions. Based on the derived expressions and making use of the multiconfigurational Dirac-Fock approach, we performed calculations of quadratic Zeeman shift coefficients and quadrupole moments for various ionic states in Ca$^{14+}$, Ni$^{12+}$ and Xe$^{q+}$ ions. These ions attract particular interest for ongoing and future experiments in optical clocks and tests of fundamental physics.
- [67] arXiv:2411.05690 [pdf, other]
-
Title: Multiple-partition cross-modulation programmable metasurface empowering wireless communicationsJun Wei Zhang, Zhen Jie Qi, Li Jie Wu, Wan Wan Cao, Xinxin Gao, Zhi Hui Fu, Jing Yu Chen, Jie Ming Lv, Zheng Xing Wang, Si Ran Wang, Jun Wei Wu, Zhen Zhang, Jia Nan Zhang, Hui Dong Li, Jun Yan Dai, Qiang Cheng, Tie Jun CuiSubjects: Applied Physics (physics.app-ph)
With the versatile manipulation capability, programmable metasurfaces are rapidly advancing in their intelligence, integration, and commercialization levels. However, as the programmable metasurfaces scale up, their control configuration becomes increasingly complicated, posing significant challenges and limitations. Here, we propose a multiple-partition cross-modulation (MPCM) programmable metasurface to enhance the wireless communication coverage with low hardware complexity. We firstly propose an innovative encoding scheme to multiply the control voltage vectors of row-column crossing, achieving high beamforming precision in free space while maintaining low control hardware complexity and reducing memory requirements for coding sequences. We then design and fabricate an MPCM programmable metasurface to confirm the effectiveness of the proposed encoding scheme. The simulated and experimental results show good agreements with the theoretically calculated outcomes in beam scanning across the E and H planes and in free-space beam pointing. The MPCM programmable metasurface offers strong flexibility and low complexity by allowing various numbers and combinations of partition items in modulation methods, catering to diverse precision demands in various scenarios. We demonstrate the performance of MPCM programmable metasurface in a realistic indoor setting, where the transmissions of videos to specific receiver positions are successfully achieved, surpassing the capabilities of traditional programmable metasurfaces. We believe that the proposed programmable metasurface has great potentials in significantly empowering the wireless communications while addressing the challenges associated with the programmable metasurface's design and implementation.
- [68] arXiv:2411.05709 [pdf, html, other]
-
Title: Exact solutions of the harmonically confined Vicsek modelComments: 13 pages, 10 figures, RevTexSubjects: Biological Physics (physics.bio-ph); Statistical Mechanics (cond-mat.stat-mech); Mathematical Physics (math-ph); Adaptation and Self-Organizing Systems (nlin.AO)
The discrete time Vicsek model confined by a harmonic potential explains many aspects of swarm formation in insects. We have found exact solutions of this model without alignment noise in two or three dimensions. They are periodic or quasiperiodic (invariant circle) solutions with positions on a circular orbit or on several concentric orbits and exist for quantized values of the confinement. There are period 2 and period 4 solutions on a line for a range of confinement strengths and period 4 solutions on a rhombus. These solutions may have polarization one, although there are partially ordered period 4 solutions and totally disordered (zero polarization) period 2 solutions. We have explored the linear stability of the exact solutions in two dimensions using the Floquet theorem and verified the stability assignements by direct numerical simulations.
- [69] arXiv:2411.05723 [pdf, other]
-
Title: On-chip rewritable phase-change metasurface for programmable diffractive deep neural networksSubjects: Optics (physics.optics)
Photonic neural networks capable of rapid programming are indispensable to realize many functionalities. Phase change technology can provide nonvolatile programmability in photonic neural networks. Integrating direct laser writing technique with phase change material (PCM) can potentially enable programming and in-memory computing for on-chip photonic neural networks. Sb2Se3 is a newly introduced ultralow-loss phase change material with a large refractive index contrast over the telecommunication transmission band. Compact, low-loss, rewritable, and nonvolatile on-chip phase-change metasurfaces can be created by using direct laser writing on a Sb2Se3 thin film. Here, by cascading multiple layers of on-chip phase-change metasurfaces, an ultra-compact on-chip programmable diffractive deep neural network is demonstrated at the wavelength of 1.55um and benchmarked on two machine learning tasks of pattern recognition and MNIST (Modified National Institute of Standards and Technology) handwritten digits classification and accuracies comparable to the state of the art are achieved. The proposed on-chip programmable diffractive deep neural network is also advantageous in terms of power consumption because of the ultralow-loss of the Sb2Se3 and its nonvolatility which requires no constant power supply to maintain its programmed state.
- [70] arXiv:2411.05727 [pdf, other]
-
Title: Cascade hot carriers via broad-band resonant tunnelingSubjects: Applied Physics (physics.app-ph)
Extraction of hot carriers (HCs) over the band-edge is a key to harvest solar energy beyond Shockley-Queisser limit1. Graphene is known as a HC-layered material due to phonon bottleneck effect near Dirac point, but limited by low photocarrier density2. Graphene/transition metal dichalcogenide (TMD) heterostructures circumvent this issue by ultrafast carrier transfer from TMD to graphene2,3. Nevertheless, efficient extraction of photocurrent by means of HCs together with carrier multiplication (CM) is still missing. Here, we introduce an ultrathin broadband resonant tunneling (BRT) barrier, TiOX to efficiently extract photocurrent with simultaneous CM and HC measurements in MoS2/graphene/TiOX heterostructure. The BRT layer gives rise to boosting open circuit voltage which is linearly proportional to incident photon energy. Meanwhile, short circuit current rises rapidly over 2Eg with obvious CM feature. This was explained by defining the joint density of states between graphene and TiOX layer over positive and negative voltage. The broadband resonant tunneling states inherently constructed from oxidation states varying from Ti3+ to Ti4+ allow the ultrafast HCs to efficiently transfer from graphene to TiOX layer. We find that the number of available tunneling states is directly proportional to short circuit current, which is well corroborated with TiOX and MoS2 thickness variance. We obtained an optimum thickness of BRT layer of ~2.8 nm, yielding cascade open circuit voltage as high as ~0.7 V, two orders of magnitude higher than that without BRT layer to reach a record efficiency of 5.3% with improved fill factor owing to synergistic HC and CM conversion under 1-SUN with long-term stability.
- [71] arXiv:2411.05737 [pdf, html, other]
-
Title: Accurate Unsupervised Photon Counting from Transition Edge Sensor SignalsNicolas Dalbec-Constant, Guillaume Thekkadath, Duncan England, Benjamin Sussman, Thomas Gerrits, Nicolás QuesadaComments: Code is publicly available at this https URLSubjects: Instrumentation and Detectors (physics.ins-det); Optics (physics.optics); Quantum Physics (quant-ph)
We compare methods for signal classification applied to voltage traces from transition edge sensors (TES) which are photon-number resolving detectors fundamental for accessing quantum advantages in information processing, communication and metrology. We quantify the impact of numerical analysis on the distinction of such signals. Furthermore, we explore dimensionality reduction techniques to create interpretable and precise photon number embeddings. We demonstrate that the preservation of local data structures of some nonlinear methods is an accurate way to achieve unsupervised classification of TES traces. We do so by considering a confidence metric that quantifies the overlap of the photon number clusters inside a latent space. Furthermore, we demonstrate that for our dataset previous methods such as the signal's area and principal component analysis can resolve up to 16 photons with confidence above $90\%$ while nonlinear techniques can resolve up to 21 with the same confidence threshold. Also, we showcase implementations of neural networks to leverage information within local structures, aiming to increase confidence in assigning photon numbers. Finally, we demonstrate the advantage of some nonlinear methods to detect and remove outlier signals.
- [72] arXiv:2411.05741 [pdf, html, other]
-
Title: Electromagnetic symmetry dislocationsSubjects: Optics (physics.optics)
Singular optics aims to understand and manipulate light's topological defects, pioneered by the discovery that phase vortex lines, strands of destructive interference, naturally occur in scalar wave fields. Monochromatic electromagnetic fields, however, are described by complex three-dimensional vectors that make individual scalar phase vortices in their vector components, which depend on the choice of co-ordinate basis, less meaningful. Instead, polarisation singularities can capture the vector texture of complicated, even non-paraxial light, with separate spatial descriptions for the electric $\mathbf{E}$ and magnetic $\mathbf{H}$ fields. But polarisation textures, too, are basis-dependent, because the laws of electromagnetism can be expressed not only by separate $\mathbf{E}$ and $\mathbf{H}$ fields, but by linear combinations of the two. We instead propose fundamental, basis-independent topological features generic in monochromatic electromagnetic fields: one- and two-dimensional structures that relate to time-averaged symmetries, including parity, duality and time-reversal, held locally by the combined electric and magnetic field polarisation geometry.
- [73] arXiv:2411.05748 [pdf, html, other]
-
Title: Multi-Dimensional Reconfigurable, Physically Composable Hybrid Diffractive Optical Neural NetworkComments: 7 pagesSubjects: Optics (physics.optics); Artificial Intelligence (cs.AI); Hardware Architecture (cs.AR)
Diffractive optical neural networks (DONNs), leveraging free-space light wave propagation for ultra-parallel, high-efficiency computing, have emerged as promising artificial intelligence (AI) accelerators. However, their inherent lack of reconfigurability due to fixed optical structures post-fabrication hinders practical deployment in the face of dynamic AI workloads and evolving applications. To overcome this challenge, we introduce, for the first time, a multi-dimensional reconfigurable hybrid diffractive ONN system (MDR-HDONN), a physically composable architecture that unlocks a new degree of freedom and unprecedented versatility in DONNs. By leveraging full-system learnability, MDR-HDONN repurposes fixed fabricated optical hardware, achieving exponentially expanded functionality and superior task adaptability through the differentiable learning of system variables. Furthermore, MDR-HDONN adopts a hybrid optical/photonic design, combining the reconfigurability of integrated photonics with the ultra-parallelism of free-space diffractive systems. Extensive evaluations demonstrate that MDR-HDONN has digital-comparable accuracy on various task adaptations with 74x faster speed and 194x lower energy. Compared to prior DONNs, MDR-HDONN shows exponentially larger functional space with 5x faster training speed, paving the way for a new paradigm of versatile, composable, hybrid optical/photonic AI computing. We will open-source our codes.
- [74] arXiv:2411.05751 [pdf, other]
-
Title: Latest progress on the reduced-order particle-in-cell scheme: I. refining the underlying formulationComments: 29 pages, 24 figuresSubjects: Plasma Physics (physics.plasm-ph); Computational Physics (physics.comp-ph)
The particle-in-cell (PIC) method is a well-established and widely used kinetic plasma modelling approach that provides a hybrid Lagrangian-Eulerian approach to solve the plasma kinetic equation. Despite its power in capturing details of the underlying physics of plasmas, conventional PIC implementations are associated with a significant computational cost, rendering their applications for real-world plasma science and engineering challenges impractical. The acceleration of the PIC method has thus become a topic of high interest, with several approaches having been pursued to this end. Among these, the concept of reduced-order (RO) PIC simulations, first introduced in 2023, provides a uniquely flexible and computationally efficient framework for kinetic plasma modelling - characteristics verified extensively in various plasma configurations. In this two-part article, we report the latest progress achieved on RO-PIC. Part I article revisits the original RO-PIC formulation and introduces refinements that substantially enhance the cost-efficiency and accuracy of the method. We discuss these refinements in comparison against the original formulation, illustrating the progression to a "first-order" implementation from the baseline "zeroth-order" one. In a detailed step-by-step verification, we first test the newly updated reduced-dimension Poisson solver (RDPS) in the first-order RO-PIC against its zeroth-order counterpart using test-case Poisson problems. Next, comparing against the zeroth-order version, we examine the performance of the complete first-order RO-PIC code in two-dimensional plasma problems. The detailed verifications demonstrate that the improvements in the RO-PIC formulation enable the approach to provide full-2D-equivalent results at a substantially lower (up to an order of magnitude) computational cost compared to the zeroth-order RO-PIC.
- [75] arXiv:2411.05756 [pdf, html, other]
-
Title: Advantages of the adoption of a generalized flame displacement velocity as a central element of flamelet theorySubjects: Fluid Dynamics (physics.flu-dyn)
In combustion theory, flames are usually described in terms of the dynamics of iso-surfaces of a specific scalar. The flame displacement speed is then introduced as a local variable quantifying the progression of these iso-surfaces relative to the flow field. While formally defined as a scalar, the physical meaning of this quantity allows relating it with a vector pointing along the normal direction of the scalar iso-surface. In this work, this one-dimensional concept is extended by the introduction of a generalized flame displacement velocity vector, which is associated with the dynamics of iso-surfaces of two generic scalars, $\alpha$ and $\beta$. It is then shown how a new flamelet paradigm can be built around this velocity vector, which leads to a very compact and generic set of two-dimensional flamelet equations for thermochemical quantities and the conditioning scalar gradients, $g_{\alpha} = \lvert \nabla \alpha \rvert$ and $g_{\beta} = \lvert \nabla \beta \rvert$. The most important features of the developed framework are discussed in the context of partially-premixed flames, which provides significant insights into several aspects of the theory, including the nature of the different contributions to the flamelet equations for the conditioning scalar gradients and the fact that different flamelet coordinate systems (orthogonal and non-orthogonal) can be characterized by the same flame displacement velocity vector. This approach opens an entire spectrum of possibilities for the definition of new two-dimensional composition spaces, which represents a very promising basis for the development of new variants of flamelet theory.
- [76] arXiv:2411.05759 [pdf, other]
-
Title: Latest progress on the reduced-order particle-in-cell scheme: II. Quasi-3D implementation and verificationComments: 24 pages, 21 figuresSubjects: Plasma Physics (physics.plasm-ph); Computational Physics (physics.comp-ph)
Across many plasma applications, the underlying phenomena and interactions among the involved processes are known to exhibit three-dimensional characteristics. Furthermore, the global properties and evolution of plasma systems are often determined by a process called inverse energy cascade, where kinetic plasma processes at the microscopic scale interact and lead to macroscopic coherent structures. These structures can have a major impact on the stability of plasma discharges, with detrimental effects on the operation and performance of plasma technologies. Kinetic particle-in-cell (PIC) methods offer a sufficient level of fidelity to capture these processes and behaviors. However, three-dimensional PIC simulations that can cost-effectively overcome the curse of dimensionality and enable full-scale simulations of real-world time significance have remained elusive. Tackling the enormous computational cost issue associated with conventional PIC schemes, the computationally efficient reduced-order (RO) PIC approach provides a viable path to 3D simulations of real-size plasma systems. This part II paper builds upon the improvements to the RO-PIC's underpinning formulation discussed in part I and extends the novel "first-order" RO-PIC formulation to 3D. The resulting Quasi-3D (Q3D) implementation is rigorously verified in this paper, both at the module level of the Q3D reduced-dimension Poisson solver (RDPS) and at the global PIC code level. The plasma test cases employed correspond to 3D versions of the 2D configurations studied in Part I, including a 3D extension to the Diocotron instability problem. The detailed verifications of the Q3D RO-PIC confirm that it maintains the expected levels of cost-efficiency and accuracy, demonstrating the ability of the approach to indistinguishably reproduce full-3D simulation results at a fraction of the computational cost.
- [77] arXiv:2411.05772 [pdf, html, other]
-
Title: Plasma Seismology: Fully Exploiting the Information Contained in Velocity Space of Kinetic Plasmas using the Morrison G TransformComments: 9 pages, 2 figures, 34 referencesSubjects: Plasma Physics (physics.plasm-ph); Solar and Stellar Astrophysics (astro-ph.SR); Space Physics (physics.space-ph)
Weakly collisional plasmas contain a wealth of information about the dynamics of the plasma in the particle velocity distribution functions, yet our ability to exploit fully that information remains relatively primitive. Here we aim to present the fundamentals of a new technique denoted Plasma Seismology that aims to invert the information from measurements of the particle velocity distribution functions at a single point in space over time to enable the determination of the electric field variation over an extended spatial region. The fundamental mathematical tool at the heart of this technique is the Morrison $G$ Transform. Using kinetic numerical simulations of Langmuir waves in a Vlasov-Poisson plasma, we demonstrate the application of the standard Morrison $G$ Transform, which uses measurements of the particle velocity distribution function over all space at one time to predict the evolution of the electric field in time. Next, we introduce a modified Morrison $G$ Transform which uses measurements of the particle velocity distribution function at one point in space over time to determine the spatial variation of the electric field over an extended spatial region. We discuss the limitations of this approach, particularly for the numerically challenging case of Langmuir waves. The application of this technique to Alfven waves in a magnetized plasma holds the promise to apply the technique to existing spacecraft particle measurement instrumentation to determine the electric fields over an extended spatial region away from the spacecraft.
New submissions (showing 77 of 77 entries)
- [78] arXiv:2411.05017 (cross-list from q-bio.NC) [pdf, html, other]
-
Title: Time in a bottle. A psychophysics study of human time perception through agingComments: 6 pages, 1 figureSubjects: Neurons and Cognition (q-bio.NC); Biological Physics (physics.bio-ph)
Time perception is crucial for a coherent human experience. As life progresses, our perception of the passage of time becomes increasingly non-uniform, often feeling as though it accelerates with age. While various causes for this phenomenon have been theorized, a comprehensive mathematical and theoretical framework remains underexplored. This study aims to elucidate the mechanisms behind perceived time dilation by integrating classical and revised psychophysical theorems with a novel mathematical approach. Utilizing Weber-Fechner laws as foundational elements, we develop a model that transitions from exponential to logarithmic functions to represent changes in time perception across the human lifespan. Our results indicate that the perception of time shifts significantly around the age of mental maturity, aligning with a proposed inversion point where sensitivity to temporal stimuli decreases, eventually plateauing out at a constant rate. This model not only explains the underlying causes of time perception changes but also provides analytical values to quantify this acceleration. These findings offer valuable insights into the cognitive and neurological processes influencing how we experience time as we go through life.
- [79] arXiv:2411.05302 (cross-list from eess.IV) [pdf, html, other]
-
Title: Adaptive Whole-Body PET Image Denoising Using 3D Diffusion Models with ControlNetSubjects: Image and Video Processing (eess.IV); Computer Vision and Pattern Recognition (cs.CV); Medical Physics (physics.med-ph)
Positron Emission Tomography (PET) is a vital imaging modality widely used in clinical diagnosis and preclinical research but faces limitations in image resolution and signal-to-noise ratio due to inherent physical degradation factors. Current deep learning-based denoising methods face challenges in adapting to the variability of clinical settings, influenced by factors such as scanner types, tracer choices, dose levels, and acquisition times. In this work, we proposed a novel 3D ControlNet-based denoising method for whole-body PET imaging. We first pre-trained a 3D Denoising Diffusion Probabilistic Model (DDPM) using a large dataset of high-quality normal-dose PET images. Following this, we fine-tuned the model on a smaller set of paired low- and normal-dose PET images, integrating low-dose inputs through a 3D ControlNet architecture, thereby making the model adaptable to denoising tasks in diverse clinical settings. Experimental results based on clinical PET datasets show that the proposed framework outperformed other state-of-the-art PET image denoising methods both in visual quality and quantitative metrics. This plug-and-play approach allows large diffusion models to be fine-tuned and adapted to PET images from diverse acquisition protocols.
- [80] arXiv:2411.05310 (cross-list from cond-mat.mtrl-sci) [pdf, other]
-
Title: Unconventional and Powerful Ion Sources for Solid-State Ion Exchange, Cu2SO4 and Cu3PO4: Exemplified by Synthesis of Metastable \beta-CuGaO2 from Stable \beta-LiGaO2Comments: 17 pages, 8 figuresSubjects: Materials Science (cond-mat.mtrl-sci); Chemical Physics (physics.chem-ph)
This study introduces a new method for synthesizing Cu+-containing metastable phases through ion exchange. Traditionally, CuCl has been used as a Cu+ ion source for solid-state ion exchanges; however, its thermodynamic driving force is often insufficient for complete ion exchange with Li+-containing precursors. First-principles calculations have identified Cu2SO4 and Cu3PO4 as more powerful alternatives, providing a higher driving force than CuCl. It has been experimentally demon-strated that these ion sources can open up new reaction pathways through experimental ion exchanges, such as from \beta-LiGaO2 to \beta-CuGaO2, which were previously unattainable. An important perspective provided by this study is that the poten-tial of such basic compounds to act as powerful ion sources has been overlooked, and that they were identified through straightforward first-principles calculations. This work presents the initial strategic design of an ion exchange reaction by exploring suitable ion sources, thereby expanding the potential for synthesizing metastable materials.
- [81] arXiv:2411.05319 (cross-list from quant-ph) [pdf, html, other]
-
Title: Dual axis atomic magnetometer and gyroscope enabled by nuclear spin perturbationSubjects: Quantum Physics (quant-ph); Applied Physics (physics.app-ph); Atomic Physics (physics.atom-ph)
Alkali-noble-gas comagnetometers have become an essential tool for tests of fundamental physics and offer a compact platform for precision gyroscopy. They are, however, limited by technical noise at low frequencies, commonly due to their limited suppression of magnetic noise. Here we investigate a new method for co-magnetometry between a single noble gas and alkali species. While similar to well-known devices using self-compensation, our scheme introduces magnetic pulses that controllably perturb the noble gas and pulsed optical pumping to polarise the alkali atoms. These applied pulses allow our scheme to measure, rather than just suppress, the effect of magnetic noise thereby offering reduced cross-talk. We show numerically that our scheme retrieves four signals (rotations and magnetic fields on two transverse axes) with similar sensitivity to a single axis device. We also present a proof-of-principle experiment based on a 87Rb-129Xe cell. Our data shows a low magnetic-rotation cross-talk of $0.2 \pm 0.1\mu$Hz$/$pT, which is already on par with the most sensitive devices relying on self-compensation.
- [82] arXiv:2411.05347 (cross-list from cond-mat.soft) [pdf, html, other]
-
Title: Bending Elasticity of the reversible Freely Jointed ChainComments: Submitted to The Journal of Chemical Physics. 14 pages, 11 figuresSubjects: Soft Condensed Matter (cond-mat.soft); Biological Physics (physics.bio-ph)
The freely jointed chain model with reversible hinges (rFJC) is the simplest theoretical model that captures reversible transitions of the local bending stiffness along the polymer chain backbone, e.g. helix-coil-type of local conformational changes or changes due to the binding/unbinding of ligands). In this work, we analyze the bending fluctuations and the bending response of a grafted rFJC in the Gibbs (fixed-force) ensemble. We obtain a recursion relation for the partition function of the grafted rFJC under bending force, which allows, in principle, exact-numerical calculation of the behavior of a rFJC of arbitrary size. In contrast to stretching, we show that under sufficiently stiff conditions, the differential bending compliance and the mean fraction of closed hinges are non-monotonic functions of the force. We also obtain the persistence length $L_p$ of the rFJC, the moments $\langle R^2 \rangle$ (mean-square end-to-end distance), and $\langle z^2 \rangle$ (mean-square transverse deflection) for the discrete chain and take the continuum limit. The tangent vector auto-correlation decays exponentially, as in the wormlike chain model (WLC). Remarkably, the expression of $\langle R^2 \rangle$ as a function of the contour length $L$ becomes the same as that in the WLC. In the thermodynamic limit, we have calculated the exact bending response analytically. As expected, for $L\gg L_p$, the boundary conditions do not matter, and the bending becomes equivalent to stretching. In contrast, for $L_p\gg L$, we have shown the non-monotonicity of the bending response (the compliance and mean fraction of closed hinges).
- [83] arXiv:2411.05420 (cross-list from cs.LG) [pdf, html, other]
-
Title: WeatherGFM: Learning A Weather Generalist Foundation Model via In-context LearningXiangyu Zhao, Zhiwang Zhou, Wenlong Zhang, Yihao Liu, Xiangyu Chen, Junchao Gong, Hao Chen, Ben Fei, Shiqi Chen, Wanli Ouyang, Xiao-Ming Wu, Lei BaiSubjects: Machine Learning (cs.LG); Artificial Intelligence (cs.AI); Computer Vision and Pattern Recognition (cs.CV); Atmospheric and Oceanic Physics (physics.ao-ph)
The Earth's weather system encompasses intricate weather data modalities and diverse weather understanding tasks, which hold significant value to human life. Existing data-driven models focus on single weather understanding tasks (e.g., weather forecasting). Although these models have achieved promising results, they fail to tackle various complex tasks within a single and unified model. Moreover, the paradigm that relies on limited real observations for a single scenario hinders the model's performance upper bound. In response to these limitations, we draw inspiration from the in-context learning paradigm employed in state-of-the-art visual foundation models and large language models. In this paper, we introduce the first generalist weather foundation model (WeatherGFM), designed to address a wide spectrum of weather understanding tasks in a unified manner. More specifically, we initially unify the representation and definition of the diverse weather understanding tasks. Subsequently, we devised weather prompt formats to manage different weather data modalities, namely single, multiple, and temporal modalities. Finally, we adopt a visual prompting question-answering paradigm for the training of unified weather understanding tasks. Extensive experiments indicate that our WeatherGFM can effectively handle up to ten weather understanding tasks, including weather forecasting, super-resolution, weather image translation, and post-processing. Our method also showcases generalization ability on unseen tasks.
- [84] arXiv:2411.05430 (cross-list from cond-mat.soft) [pdf, html, other]
-
Title: Manifestations of the possible thermodynamic origin of water's anomalies in non-classical vapor nucleation at negative pressuresSubjects: Soft Condensed Matter (cond-mat.soft); Chemical Physics (physics.chem-ph)
Over the years, various scenarios -- such as the stability-limit conjecture (SLC), two critical point (TCP), critical point-free (CPF), and singularity-free (SF) -- have been proposed to explain the thermodynamic origin of supercooled waters anomalies. However, direct experimental validation is challenging due to the rapid phase transition from metastable water. In this study, we explored whether the phase transition pathways from metastable water provide insight into the thermodynamic origin of these anomalies. Using a classical density functional theory approach with realistic theoretical water models, we examined how different thermodynamic scenarios influence vapor nucleation kinetics at negative pressures. Our findings show significant variations in nucleation kinetics and mechanism during both isobaric and isochoric cooling. In the TCP scenario, the nucleation barrier increases steadily during isobaric cooling, with a slight decrease near the Widom line at lower temperatures (Ts). In contrast, the SF scenario shows a monotonic increase in the nucleation barrier. For the CPF scenario, we observed a non-classical mechanism, such as wetting-mediated nucleation (where the growing vapor nucleus is wetted by the intermediate low-density liquid phase) and the Ostwald step rule at low temperatures. Isochoric cooling pathways also revealed notable differences in T-dependent nucleation barrier trends between the TCP and CPF scenarios. Overall, this study underscores the importance of analyzing phase transition kinetics and mechanism to understand the precise thermodynamic origin of supercooled waters anomalies.
- [85] arXiv:2411.05462 (cross-list from math.OC) [pdf, html, other]
-
Title: Tracking disturbances in transmission networksSubjects: Optimization and Control (math.OC); Mathematical Physics (math-ph); Applied Physics (physics.app-ph)
We study the nonlinear inverse source problem of detecting, localizing and identifying unknown accidental disturbances on forced and damped transmission networks. A first result is that strategic observation sets are enough to guarantee detection of disturbances. To localize and identify them, we additionally need the observation set to be absorbent. If this set is dominantly absorbent, then detection, localization and identification can be done in "quasi real-time". We illustrate these results with numerical experiments.
- [86] arXiv:2411.05468 (cross-list from quant-ph) [pdf, html, other]
-
Title: Noncommuting common causes revisitedComments: 34+7 pages, 2 figuresSubjects: Quantum Physics (quant-ph); Mathematical Physics (math-ph); History and Philosophy of Physics (physics.hist-ph)
In this paper, we revisit the concept of noncommuting common causes; refute two objections raised against them, the triviality objection and the lack of causal explanatory force; and explore how their existence modifies the EPR argument. More specifically, we show that 1) product states screening off all quantum correlations do not compromise noncommuting common causal explanations; 2) noncommuting common causes can satisfy the law of total probability; 3) perfect correlations can have indeterministic noncommuting common causes; and, as a combination of the above claims, 4) perfect correlations can have noncommuting common causes which are both nontrivial and satisfy the law of total probability.
- [87] arXiv:2411.05505 (cross-list from astro-ph.EP) [pdf, html, other]
-
Title: Planetesimal gravitational collapse in a gaseous environment: Thermal and dynamic evolutionComments: 12 pages, 12 figures, accepted for publication in A&ASubjects: Earth and Planetary Astrophysics (astro-ph.EP); Geophysics (physics.geo-ph)
Planetesimal formation models often invoke the gravitational collapse of pebble clouds to overcome various barriers to grain growth and propose processes to concentrate particles sufficiently to trigger this collapse. On the other hand, the geochemical approach for planet formation constrains the conditions for planetesimal formation and evolution by providing temperatures that should be reached to explain the final composition of planetesimals, the building blocks of planets. To elucidate the thermal evolution during gravitational collapse, we used numerical simulations of a self-gravitating cloud of particles and gas coupled with gas drag. Our goal is to determine how the gravitational energy relaxed during the contraction is distributed among the different energy components of the system, and how this constrains a thermal and dynamical planetesimal's history. We identify the conditions necessary to achieve a temperature increase of several hundred kelvins, and as much as 1600 K. Our results emphasise the key role of the gas during the collapse.
- [88] arXiv:2411.05518 (cross-list from cond-mat.mtrl-sci) [pdf, other]
-
Title: Piezovalley effect and magnetovalley coupling in altermagnetic semiconductorsComments: 14 pages, 6 figuresSubjects: Materials Science (cond-mat.mtrl-sci); Computational Physics (physics.comp-ph)
Clarifying the physical origin of valley polarization and exploring new ferrovalley materials are conducive to the application of valley degrees of freedom in the field of information storage. Here, we explore two new-type altermagnetic semiconductors (monolayers Nb2Se2O and Nb2SeTeO) with above-room-temperature Néel temperature based on first-principles calculations. It reveals that uniaxial strain induces valley polarization without spin-orbital coupling (SOC) in altermagnets owing to the piezovalley effect, while uniaxial compressive strain transforms from ferrovalley semiconductor to semimetal, half metal and metal. Moreover, moderate biaxial strain renders Janus monolayer Nb2SeTeO to robust Dirac-like band dispersion. The SOC and intrinsic in-plane magnetocrystalline anisotropy energy induces Dirac-like altermagnets to generate apparent valley polarization through magnetovalley coupling. In terms of SOC perturbation, we elucidate the physical mechanism behind in-plane-magnetization induced valley polarization and point out the magnitude of valley polarization is positively correlated with the square of SOC strength and negatively correlated with the bandgap. The present work reveals the physical origin of valley polarization in altermagnets and expands the application of ferrovalley at room temperature in valleytronics.
- [89] arXiv:2411.05536 (cross-list from cs.LG) [pdf, html, other]
-
Title: Towards Active Flow Control Strategies Through Deep Reinforcement LearningComments: ECOMMAS 2024 conference proceeding paperSubjects: Machine Learning (cs.LG); Fluid Dynamics (physics.flu-dyn)
This paper presents a deep reinforcement learning (DRL) framework for active flow control (AFC) to reduce drag in aerodynamic bodies. Tested on a 3D cylinder at Re = 100, the DRL approach achieved a 9.32% drag reduction and a 78.4% decrease in lift oscillations by learning advanced actuation strategies. The methodology integrates a CFD solver with a DRL model using an in-memory database for efficient communication between
- [90] arXiv:2411.05542 (cross-list from quant-ph) [pdf, html, other]
-
Title: Quantum magnetometry of transient signals with a time resolution of 1.1 nanosecondsKonstantin Herb, Laura A. Völker, John M. Abendroth, Nicholas Meinhardt, Laura van Schie, Pietro Gambardella, Christian L. DegenComments: 11 pages including methodsSubjects: Quantum Physics (quant-ph); Mesoscale and Nanoscale Physics (cond-mat.mes-hall); Materials Science (cond-mat.mtrl-sci); Instrumentation and Detectors (physics.ins-det)
Quantum magnetometers based on spin defects in solids enable sensitive imaging of various magnetic phenomena, such as ferro- and antiferromagnetism, superconductivity, and current-induced fields. Existing protocols primarily focus on static fields or narrow-band dynamical signals, and are optimized for high sensitivity rather than fast time resolution. Here, we report detection of fast signal transients, providing a perspective for investigating the rich dynamics of magnetic systems. We experimentally demonstrate our technique using a single nitrogen-vacancy (NV) center magnetometer at room temperature, reaching a best-effort time resolution of 1.1 ns, an instantaneous bandwidth of 0.9 GHz, and a time-of-flight precision better than 20 ps. The time resolution can be extended to the picosecond range by use of on-chip waveguides. At these speeds, NV quantum magnetometers will become competitive with time-resolved synchrotron X-ray techniques. Looking forward, adding fast temporal resolution to the spatial imaging capability further promotes single-spin probes as powerful research tools in spintronics, mesoscopic physics, and nanoscale device metrology.
- [91] arXiv:2411.05543 (cross-list from q-bio.PE) [pdf, html, other]
-
Title: Evolution of cooperation in a three-strategy game combining snowdrift and stag hunt gamesComments: 10 pages, 6 figuresSubjects: Populations and Evolution (q-bio.PE); Physics and Society (physics.soc-ph)
This study aimed to investigate the evolutionary dynamics of a three-strategy game that combines snowdrift and stag hunt games. This game is motivated by an experimental study, which found that individual solution lowers cooperation levels. Agents adopting this option aim to address a problem to the extent necessary to remove negative impact on themselves, although they do not free ride on cooperation effort provided by others. This property of the individual solution is similar to that of option defection in the stag hunt. Thus, the role of the interplay of defection in the snowdrift game and individual solution was examined in this study. The well-mixed population has two asymptotically stable rest points, one wherein the individual solution occupies the population, and the other wherein cooperation and defection coexist. The interactions on a square lattice enlarge the parameter region wherein cooperation survives, and the three strategies often coexist. The scrutinization of the evolutionary process shows that multiple mechanisms lead to the coexistence of the three strategies depending on parameter values. Our analysis suggests that considering the individual solution adds complexity to the evolutionary process, which might contribute to our understanding on the evolution of cooperation.
- [92] arXiv:2411.05587 (cross-list from cond-mat.soft) [pdf, html, other]
-
Title: Smart navigation through a rotating barrier: Deep reinforcement learning with application to size-based separation of active microagentsComments: 12 pages, 5 figuresSubjects: Soft Condensed Matter (cond-mat.soft); Statistical Mechanics (cond-mat.stat-mech); Biological Physics (physics.bio-ph); Computational Physics (physics.comp-ph)
We employ deep reinforcement learning methods to investigate shortest-time navigation strategies for smart active Brownian particles (agents) that self-propel through a rotating, localized potential barrier in an otherwise static and viscous fluid background. The particle motion is prescribed to begin from a specified start point and terminate at a specified destination point. The potential barrier is modeled as a repulsive Gaussian potential with a finite range that rotates with given angular velocities around a fixed center within the plane of motion. We use an advantage actor-critic approach to train the agents and demonstrate that, by leveraging this approach, the rotating potential can be utilized for size-based sorting and separation of the agents. In other words, agents of different (hydrodynamic) radii reach the end-point at sufficiently well-separated times, enabling their separation. The efficiency of particle separation in this context is discussed using specific criteria. We also study the effect of rotational Brownian noise on the quality of the proposed size-based sorting mechanism. Additionally, we demonstrate how the use of noise-induced training can enhance this mechanism in a noisy background within the range of parameters explored in our work. Reinforcement learning in the context of active particles thus offers promising avenues to unravel optimal navigation strategies within complex environments, with potential applications that can be utilized in microfluidic settings.
- [93] arXiv:2411.05704 (cross-list from astro-ph.HE) [pdf, html, other]
-
Title: Understanding Streaming Instabilities in the Limit of High Cosmic Ray Current DensitySubjects: High Energy Astrophysical Phenomena (astro-ph.HE); Plasma Physics (physics.plasm-ph)
A critical component of particle acceleration in astrophysical shocks is the non-resonant (Bell) instability, where the streaming of cosmic rays (CRs) leads to the amplification of magnetic fields necessary to scatter particles. In this work we use kinetic particle-in-cells simulations to investigate the high-CR current regime, where the typical assumptions underlying the Bell instability break down. Despite being more strongly driven, significantly less magnetic field amplification is observed compared to low-current cases, an effect due to the anisotropic heating that occurs in this regime. We also find that electron-scale modes, despite being fastest growing, mostly lead to moderate electron heating and do not affect the late evolution or saturation of the instability.
- [94] arXiv:2411.05711 (cross-list from cond-mat.mtrl-sci) [pdf, html, other]
-
Title: How grain structure evolution affects kinetics of a solid-state reaction: a case of interaction between iridium and zirconium carbideComments: 27 pages, 7 figuresSubjects: Materials Science (cond-mat.mtrl-sci); Chemical Physics (physics.chem-ph)
This work investigates the solid-state reaction between iridium and zirconium carbide, resulting in the formation of carbon and $\mathrm{ZrIr}_{3}$ -- an intermetallic compound of great interest for modern high-temperature materials science. We have found a transition of kinetic regimes in this reaction: from linear kinetics (when the chemical reaction is a limiting stage) at 1500 and 1550°C to `non-parabolic kinetics' at 1600°C. Non-parabolic kinetics is characterized by thickness of a product layer being proportional to a power of time less than 1/2. The nature of non-parabolic kinetics was still an open question, which motivated us to develop a model of this kinetic regime. The proposed model accounts for the grain growth in the product phase and how it leads to the time dependence of the interdiffusion coefficient. We have obtained a complete analytic solution for this model and an equation that connects the grain-growth exponent and the power-law exponent of non-parabolic kinetics. The measurements of the thickness of the product layer and the average grain size of the intermetallic phase confirm the results of the theoretical solution.
- [95] arXiv:2411.05717 (cross-list from cond-mat.soft) [pdf, html, other]
-
Title: Using optical tweezer electrophoresis to investigate clay nanoplatelet adsorption on Latex microspheres in aqueous mediaComments: 32 pages, 14 figures and supporting information includedSubjects: Soft Condensed Matter (cond-mat.soft); Applied Physics (physics.app-ph)
The adsorption of charged clay nanoplatelets plays an important role in stabilizing emulsions by forming a barrier around the emulsion droplets and preventing coalescence. In this work, the adsorption of charged clay nanoplatelets on a preformed Latex microsphere in an aqueous medium is investigated at high temporal resolution using optical tweezer-based single-colloid electrophoresis. Above a critical clay concentration, charged clay nanoplatelets in an aqueous medium self-assemble gradually to form gel-like networks that become denser with increasing medium salinity. In a previous publication [R. Biswas et. al., Soft Matter, 2023, 19, 24007-2416], some of us had demonstrated that a Latex microsphere, optically trapped in a clay gel medium, is expected to attach to the network strands of the gel. In the present contribution, we show that for different ionic conditions of the suspending medium, the adsorption of clay nanoplatelets increases the effective surface charge on an optically trapped Latex microsphere while also enhancing the drag experienced by the latter. Besides the ubiquitous contribution of non-electrostatic dispersion forces in driving the adsorption process, we demonstrate the presence of an electrostatically-driven adsorption mechanism when the microsphere was trapped in a clay gel. These observations are qualitatively verified via cryogenic field emission scanning electron microscopy and are useful in achieving colloidal stabilisation, for example, during the preparation of clay-armoured Latex particles in Pickering emulsion polymerisation.
- [96] arXiv:2411.05728 (cross-list from quant-ph) [pdf, html, other]
-
Title: Quantum hyperspins: Highly nonclassical collective behavior in quantum optical parametric oscillatorsComments: 6 pages, 3 figuresSubjects: Quantum Physics (quant-ph); Computational Physics (physics.comp-ph); Optics (physics.optics)
We report on the emergence of a highly non-classical collective behavior in quantum parametric oscillators, which we name quantum hyperspin, induced by a tailored nonlinear interaction. This is the second quantized version of classical multidimensional spherical spins, as XY spins in two dimensions, and Heisenberg spins in three dimensions. In the phase space, the quantum hyperspins are represented as spherical shells whose radius scales with the number of particles in a way such that it cannot be factorized even in the limit of large particle number. We show that the nonlinearly coupled quantum oscillators form a high-dimensional entangled state that is surprisingly robust with respect to the coupling with the environment. Such a behavior results from a properly engineered reservoir. Networks of entangled quantum hyperspins are a new approach to quantum simulations for applications in computing, Ising machines, and high-energy physics models. We analyze from first principles through ab initio numerical simulations the properties of quantum hyperspins, including the interplay of entanglement and coupling frustration.
- [97] arXiv:2411.05761 (cross-list from math.NA) [pdf, html, other]
-
Title: The Helmholtz Dirichlet and Neumann problems on piecewise smooth open curvesComments: 27 pages, 39 figuresSubjects: Numerical Analysis (math.NA); Computational Physics (physics.comp-ph)
A numerical scheme is presented for solving the Helmholtz equation with Dirichlet or Neumann boundary conditions on piecewise smooth open curves, where the curves may have corners and multiple junctions. Existing integral equation methods for smooth open curves rely on analyzing the exact singularities of the density at endpoints for associated integral operators, explicitly extracting these singularities from the densities in the formulation, and using global quadrature to discretize the boundary integral equation. Extending these methods to handle curves with corners and multiple junctions is challenging because the singularity analysis becomes much more complex, and constructing high-order quadrature for discretizing layer potentials with singular and hypersingular kernels and singular densities is nontrivial. The proposed scheme is built upon the following two observations. First, the single-layer potential operator and the normal derivative of the double-layer potential operator serve as effective preconditioners for each other locally. Second, the recursively compressed inverse preconditioning (RCIP) method can be extended to address "implicit" second-kind integral equations. The scheme is high-order, adaptive, and capable of handling corners and multiple junctions without prior knowledge of the density singularity. It is also compatible with fast algorithms, such as the fast multipole method. The performance of the scheme is illustrated with several numerical examples.
Cross submissions (showing 20 of 20 entries)
- [98] arXiv:2211.12017 (replaced) [pdf, html, other]
-
Title: On the Threshold of Drop Fragmentation under Impulsive AccelerationComments: 40 pages, 22 figuresSubjects: Fluid Dynamics (physics.flu-dyn)
Secondary fragmentation of an impulsively accelerated drop depends on fluid properties and velocity of the ambient. The critical Weber number $(\mathit{We}_\mathit{cr})$, the minimum Weber number at which a drop undergoes non-vibrational breakup, depends on fluid density ratio $(\rho)$, the drop $(\mathit{Oh}_d)$, and the ambient $(\mathit{Oh}_o)$ Ohnesorge numbers. The current study uses VoF based interface-tracking multiphase flow simulations to quantify the effect of different non-dimensional groups on the threshold at which secondary fragmentation occur. For $\mathit{Oh}_d \leq 0.1$, a decrease in $\mathit{Oh}_d$ significantly influences the breakup morphology, plume formation, and the resulting $\mathit{We}_\mathit{cr}$. $\rho$ and $\mathit{Oh}_o$, were found to influence the balance between the pressure differences between the poles and the periphery, and the shear stresses on the upstream surface. These external forces induce flow inside the initially spherical drop, resulting in deformation into pancakes and eventually the breakup morphology of forward/backward bag. The evolution pathways of the drop morphology based on their non-dimensional groups have been charted. With the inclusion of the new observations, the traditional $\mathit{We}_\mathit{cr}-\mathit{Oh}_d$ plot, used for illustrating the dependence of critical Weber number on $\mathit{Oh}_d$, was found to be inadequate in predicting the minimum initial $\mathit{We}$ required to undergo fragmentation. A new non-dimensional parameter $C_{breakup}$ is derived based on the competition between forces driving drop deformation and the forces resisting drop deformation. Tested on the available experimental data and current simulations, $C_{breakup}$ is found to be a robust predictor for the threshold of drop fragmentation.
- [99] arXiv:2309.14027 (replaced) [pdf, html, other]
-
Title: TomOpt: Differential optimisation for task- and constraint-aware design of particle detectors in the context of muon tomographyGiles C. Strong, Maxime Lagrange, Aitor Orio, Anna Bordignon, Florian Bury, Tommaso Dorigo, Andrea Giammanco, Mariam Heikal, Jan Kieseler, Max Lamparth, Pablo Martínez Ruíz del Árbol, Federico Nardi, Pietro Vischia, Haitham ZaraketComments: V3: Updated to published version; 29 pages content + appendix + refsJournal-ref: Mach. Learn.: Sci. Technol. 5 035002 (2024)Subjects: Instrumentation and Detectors (physics.ins-det); High Energy Physics - Experiment (hep-ex); Machine Learning (stat.ML)
We describe a software package, TomOpt, developed to optimise the geometrical layout and specifications of detectors designed for tomography by scattering of cosmic-ray muons. The software exploits differentiable programming for the modeling of muon interactions with detectors and scanned volumes, the inference of volume properties, and the optimisation cycle performing the loss minimisation. In doing so, we provide the first demonstration of end-to-end-differentiable and inference-aware optimisation of particle physics instruments. We study the performance of the software on a relevant benchmark scenario and discuss its potential applications. Our code is available on Github.
- [100] arXiv:2311.14174 (replaced) [pdf, html, other]
-
Title: Rigorous Hydrodynamics from Linear Boltzmann Equations and Viscosity-Capillarity BalanceSubjects: Fluid Dynamics (physics.flu-dyn); Mathematical Physics (math-ph)
An exact closure for hydrodynamic variables is rigorously derived from the linear Boltzmann kinetic equation. Our approach, based on spectral theory, structural properties of eigenvectors and the theory of slow manifolds, allows us to define a unique, optimal reduction in phase space close to equilibrium. The hydrodynamically constrained system induces a modification of entropy that ensures pure viscous dissipation on the hydrodynamic manifold, which is interpreted as a non-local variant of Korteweg's theory of viscosity-capillarity balance. The rigorous hydrodynamic equations are exemplified on the Knudsen minimum paradox in a channel flow.
- [101] arXiv:2402.01029 (replaced) [pdf, html, other]
-
Title: Response Theory via Generative Score ModelingComments: In press. Includes supplementary material in the file this http URLJournal-ref: Phys. Rev. Lett. 2024Subjects: Data Analysis, Statistics and Probability (physics.data-an); Machine Learning (cs.LG)
We introduce an approach for analyzing the responses of dynamical systems to external perturbations that combines score-based generative modeling with the Generalized Fluctuation-Dissipation Theorem (GFDT). The methodology enables accurate estimation of system responses, including those with non-Gaussian statistics. We numerically validate our approach using time-series data from three different stochastic partial differential equations of increasing complexity: an Ornstein-Uhlenbeck process with spatially correlated noise, a modified stochastic Allen-Cahn equation, and the 2D Navier-Stokes equations. We demonstrate the improved accuracy of the methodology over conventional methods and discuss its potential as a versatile tool for predicting the statistical behavior of complex dynamical systems.
- [102] arXiv:2403.12084 (replaced) [pdf, html, other]
-
Title: Singular parametric oscillators from the one-parameter Darboux transformation of the classical harmonic oscillatorComments: 16 pages, 10 figures, 35 references, published version, a few remaining misprints correctedJournal-ref: Annals of Physics 470 (2024) 169830Subjects: Classical Physics (physics.class-ph); Exactly Solvable and Integrable Systems (nlin.SI)
The singular parametric oscillators obtained from the one-parameter Darboux deformation/transformation effected upon the classical harmonic oscillator are introduced and discussed in some detail using sin(omega_0 t) and cos(omega_0 t) as seed solutions. The corresponding Ermakov-Lewis integrability problem of these parametric oscillators is also studied. It is shown that the Ermakov-Lewis invariants do not depend on the deformation parameter and are singularity-free.
- [103] arXiv:2404.07416 (replaced) [pdf, html, other]
-
Title: Highly spin-polarized molecules via collisional microwave pumpingSubjects: Chemical Physics (physics.chem-ph); Atomic Physics (physics.atom-ph); Quantum Physics (quant-ph)
We propose a general technique to produce cold spin-polarized molecules, in which rotationally excited states are first populated by coherent microwave excitation, and then allowed to spin-flip and relax via collisional quenching, which populates a single final spin state. We illustrate the high selectivity of the technique for $^{13}$C$^{16}$O molecules immersed in a cold buffer gas of helium atoms, achieving a high degree ($\geq$95\%) of nuclear spin polarization at 1 K.
- [104] arXiv:2404.10945 (replaced) [pdf, html, other]
-
Title: Dynamics of an oscillatory boundary layer over a sediment bedSubjects: Fluid Dynamics (physics.flu-dyn)
We investigate the dynamics of an oscillatory boundary layer developing over a bed of collisional and freely evolving sediment grains. We perform Euler-Lagrange simulations at Reynolds numbers $\mathrm{Re}_\delta= 200$, 400, and 800, density ratio $\rho_p/\rho_f = 2.65$, Galileo number $\mathrm{Ga} = 51.9$, maximum Shields numbers from $2.79 \times 10^{-2}$ to $4.40 \times 10^{-2}$, and Keulegan-Carpenter number from $134.5$ to $538.0$. We show that the dynamics of the oscillatory boundary layer and sediment bed are strongly coupled due to two mechanisms: (I) bed permeability, which leads to flow penetration deep inside the sediment layer, a slip velocity at the bed-fluid interface, and the expansion of the boundary layer, and (II) particle motion, which leads to rolling-grain ripples at $\mathrm{Re}_\delta = 400$ and sheet flow at $\mathrm{Re}_\delta = 800$. While at $\mathrm{Re}_\delta = 200$ the sediment bed remains static during the entire cycle, the permeability of the bed-fluid interface causes a thickening of the boundary layer and reduction of the bed shear stress by 38%. With increasing $\mathrm{Re}_\delta$, the particles become mobile, which leads to rolling-grain ripples at $\mathrm{Re}_\delta = 400$ and suspended sediment at $\mathrm{Re}_\delta = 800$. Due to their feedback force on the fluid, the mobile sediment particles cause greater velocity fluctuations in the fluid. In addition, greater flow intrusion causes the bed shear stress to decrease by 54% at $\mathrm{Re}_\delta = 400$and 79% at $\mathrm{Re}_\delta = 800$, compared to the flow over a smooth and impermeable wall.
- [105] arXiv:2404.11832 (replaced) [pdf, html, other]
-
Title: Prolate spheroids settling in a quiescent fluid: clustering, microstructures and collisionsSubjects: Fluid Dynamics (physics.flu-dyn)
In this study, we investigate the sedimentation of prolate spheroids in a quiescent fluid by means of the particle-resolved direct numerical simulation. With the increase of the particle volume fraction $\phi$ from $0.1\%$ to $10\%$, we observe a non-monotonic variation of the mean settling velocity of particles, $\langle V_s \rangle$. By virtue of the Voronoi analysis, we find that the degree of particle clustering is highest when $\langle V_s \rangle$ reaches the local maximum at $\phi=1\%$. Under the swarm effect, clustered particles are found to preferentially sample downward fluid flows in the wake regions, leading to the enhancement of the settling speed. As for lower or higher volume fraction, the tendency of particle clustering and the preferential sampling of downward flows are attenuated. Hindrance effect becomes predominant when the volume fraction exceeds 5\% and reduces $\langle V_s \rangle$ to less than the isolated settling velocity. Particle orientation plays a minor role in the mean settling velocity, although individual prolate particles still tend to settle faster in suspensions when they deviate more from the broad-side-on alignment. Moreover, we also demonstrate that particles are prone to form column-like microstructures in dilute suspensions under the effect of wake-induced hydrodynamic attractions. The radial distribution function is higher at a lower volume fraction. As a result, the collision rate scaled by the particle number density decreases with the increasing volume fraction. By contrast, as another contribution to the particle collision rate, the relative radial velocity for nearby particles shows a minor degree of variation due to the lubrication effect.
- [106] arXiv:2405.01377 (replaced) [pdf, html, other]
-
Title: Nonlinearity-induced symmetry breaking in a system of two parametrically driven Kerr-Duffing oscillatorsJournal-ref: New Journal of Physics, 26, 103020 (2024)Subjects: Classical Physics (physics.class-ph); Mesoscale and Nanoscale Physics (cond-mat.mes-hall); Optics (physics.optics)
We study the classical dynamics of a system comprising a pair of Kerr-Duffing nonlinear oscillators, which are coupled through a nonlinear interaction and subjected to a parametric drive. Using the rotating wave approximation (RWA), we analyze the steady-state solutions for the amplitudes of the two oscillators. For the case of almost identical oscillators, we investigate separately the cases in which only one oscillator is parametrically driven and in which both oscillators are simultaneously driven. In the latter regime, we demonstrate that even when the parametric drives acting on the two oscillators are identical, the system can transition from a stable Nesymmetric solution to a broken-symmetry solution as the detuning is varied.
- [107] arXiv:2407.09309 (replaced) [pdf, html, other]
-
Title: Dynamic-Mode Decomposition of Geostrophically Balanced Motions from SWOT AltimetryTakaya Uchida, Yadidya Badarvada, Karl E. Lapo, Xiaobiao Xu, Jeffrey J. Early, Brian K. Arbic, Dimitris Menemenlis, Luna Hiron, Eric P. Chassignet, Jay F. Shriver, Maarten C. BuijsmanSubjects: Atmospheric and Oceanic Physics (physics.ao-ph); Fluid Dynamics (physics.flu-dyn)
The decomposition of oceanic flow into its balanced and unbalanced motions carries theoretical and practical significance for the oceanographic community. These two motions have distinct dynamical characteristics and affect the transport of tracers differently from one another. The launch of Surface Water and Ocean Topography (SWOT) satellite provides a prime opportunity to diagnose the surface balanced and unbalanced motions on a global scale at an unprecedented spatial resolution. Here, we apply dynamic-mode decomposition (DMD), a linear-algebraic data-driven method, to a tidally-forced numerical simulation and one-day-repeat SWOT observations of sea-surface height (SSH) in the Gulf Stream extension. DMD is able to separate out the spatial modes associated with sub-inertial periods from super-inertial periods. The sub-inertial modes of DMD can be used to extract geostrophically balanced motions from SSH fields, which have an imprint of internal tides and gravity waves. We utilize the statistical relation between relative vorticity and strain rate as the metric to gauge the extraction of geostrophy.
- [108] arXiv:2409.10483 (replaced) [pdf, other]
-
Title: Efficient polarizable QM/MM using the direct reaction field Hamiltonian with electrostatic potential fitted multipole operatorsSubjects: Chemical Physics (physics.chem-ph)
Electronic polarization and dispersion are decisive actors in determining interaction energies between molecules. These interactions have a particularly profound effect on excitation energies of molecules in complex environments, especially when the excitation involves a significant degree of charge reorganisation. The direct reaction field (DRF) approach, which has seen a recent revival of interest, provides a powerful framework for describing these interactions in quantum mechanics/molecular mechanics (QM/MM) models of systems, where a small subsystem of interest is described using quantum chemical methods and the remainder is treated with a simple MM force field. In this paper we show how the DRF approach can be combined with the electrostatic potential fitted (ESPF) multipole operator description of the QM region charge density, which significantly improves the efficiency of the method, particularly for large MM systems, and for typical calculations effectively eliminates the dependence on MM system size. We also show how the DRF approach can be combined with fluctuating charge descriptions of the polarizable environment, as well as previously used atom-centred dipole-polarizability based models. We further show that the ESPF-DRF method provides an accurate description of molecular interactions in both ground and excited electronic states of the QM system and apply it to predict the gas to aqueous solution solvatochromic shifts in the UV/visible absorption spectrum of acrolein.
- [109] arXiv:2410.04362 (replaced) [pdf, html, other]
-
Title: What is a physical theory? Philosophers do have an answer by distinguishing two forms of empirical adequacyComments: 20 pagesSubjects: History and Philosophy of Physics (physics.hist-ph)
What is a physical theory? Although this intriguing question has been addressed from many different perspectives, some physicists and philosophers of physics have implicitly or explicitly embraced a philosophically-neutral definition of a physical theory, independently of the philosophical position endorsed with respect to it. Considering some objections against this view, I shall argue that the most appropriate definition of a physical theory already presupposes some commitments shared by a philosophical position associated with scientific realism. As we shall see, what physical theories and scientific realist positions have in common is the commitment of satisfying a non-factive notion of empirical adequacy, whilst a factive notion of empirical adequacy shall be solely associated with scientific realism. Based on this factive/non-factive distinction, we shall finally present a case study in physics in order to show that a well-known foundational problem associated with the concept of manifest time can be dissolved once the non-factive notion of empirical adequacy (as opposed to factive one) is endorsed.
- [110] arXiv:2410.05862 (replaced) [pdf, html, other]
-
Title: Deliberation Among Informed Citizens: The Value of Exploring Alternative Thinking FramesComments: 28+2 pages, 3 figuresSubjects: Physics and Society (physics.soc-ph)
We investigate the potential of deliberation to create consensus among fully-informed citizens. Our approach relies on two cognitive assumptions: i. citizens need a thinking frame (or perspective) to consider an issue; and ii. citizens cannot consider all relevant perspectives simultaneously: they are incompatible in the mind. These assumptions imply that opinions are intrinsically contextual. Formally, we capture contextuality in a simple quantum-like cognitive model. We consider a binary voting problem, in which two citizens with incompatible thinking frames and initially opposite voting intentions deliberate under the guidance of a benevolent facilitator. We find that when citizens consider alternative perspectives, their opinion may change. When the citizens' perspectives are two-dimensional and maximally uncorrelated, the probability for consensus after two rounds of deliberation reaches 75%; and this probability increases proportionally with the dimensionality (namely, the richness) of the perspectives. When dealing with a population of citizens, we also elaborate a novel rationale for working in subgroups. The contextuality approach delivers a number of insights. First, the diversity of perspectives is beneficial and even necessary for deliberations to overcome initial disagreement. Second, successful deliberation demand the active participation of citizens in terms of "putting themselves in the other's shoes". Third, well-designed procedures managed by a facilitator are necessary to secure increased probability for consensus. A last insight is that the richness of citizens' thinking frames is beneficial, while the optimal strategy entails focusing deliberation on a properly reduced problem.
- [111] arXiv:2410.14370 (replaced) [pdf, html, other]
-
Title: Design rules for low-insertion-loss magnonic transducersRóbert Erdélyi, György Csaba, Levente Maucha, Felix Kohl, Björn Heinz, Johannes Greil, Markus Becherer, Philipp Pirro, Ádám PappSubjects: Applied Physics (physics.app-ph)
We present a computational framework for the design of magnonic transducers, where waveguide antennas generate and pick up spin-wave signals. Our method relies on the combination of circuit-level models with micromagnetic simulations and allows simulation of complex geometries in the magnonic domain. We validated our model with experimental measurements, which showed good agreement witch the predicted scattering parameters of the system. Using our model we identified scaling rules of the antenna radiation resistance and we show strategies to maximize transduction efficiency between the electric and magnetic domains. We designed a transducer pair on YIG with 5dB insertion loss in a 100 MHz band, an unusually low value for micron-scale spin-wave devices. This demonstrates that magnonic devices can be very efficient and competitive in RF applications.
- [112] arXiv:2410.17367 (replaced) [pdf, html, other]
-
Title: Generalizing Geometric Partition Entropy for the Estimation of Mutual Information in the Presence of Informative OutliersSubjects: Physics and Society (physics.soc-ph); Information Theory (cs.IT); Chaotic Dynamics (nlin.CD)
The recent introduction of geometric partition entropy brought a new viewpoint to non-parametric entropy quantification that incorporated the impacts of informative outliers, but its original formulation was limited to the context of a one-dimensional state space. A generalized definition of geometric partition entropy is now provided for samples within a bounded (finite measure) region of a d-dimensional vector space. The basic definition invokes the concept of a Voronoi diagram, but the computational complexity and reliability of Voronoi diagrams in high dimension make estimation by direct theoretical computation unreasonable. This leads to the development of approximation schemes that enable estimation that is faster than current methods by orders of magnitude. The partition intersection ($\pi$) approximation, in particular, enables direct estimates of marginal entropy in any context resulting in an efficient and versatile mutual information estimator. This new measure-based paradigm for data driven information theory allows flexibility in the incorporation of geometry to vary the representation of outlier impact, which leads to a significant broadening in the applicability of established entropy-based concepts. The incorporation of informative outliers is illustrated through analysis of transient dynamics in the synchronization of coupled chaotic dynamical systems.
- [113] arXiv:2411.02411 (replaced) [pdf, other]
-
Title: SK-PINN: Accelerated physics-informed deep learning by smoothing kernel gradientsSubjects: Computational Physics (physics.comp-ph)
The automatic differentiation (AD) in the vanilla physics-informed neural networks (PINNs) is the computational bottleneck for the high-efficiency analysis. The concept of derivative discretization in smoothed particle hydrodynamics (SPH) can provide an accelerated training method for PINNs. In this paper, smoothing kernel physics-informed neural networks (SK-PINNs) are established, which solve differential equations using smoothing kernel discretization. It is a robust framework capable of solving problems in the computational mechanics of complex domains. When the number of collocation points gradually increases, the training speed of SK-PINNs significantly surpasses that of vanilla PINNs. In cases involving large collocation point sets or higher-order problems, SK-PINN training can be up to tens of times faster than vanilla PINN. Additionally, analysis using neural tangent kernel (NTK) theory shows that the convergence rates of SK-PINNs are consistent with those of vanilla PINNs. The superior performance of SK-PINNs is demonstrated through various examples, including regular and complex domains, as well as forward and inverse problems in fluid dynamics and solid mechanics.
- [114] arXiv:2411.03389 (replaced) [pdf, other]
-
Title: Neurons for Neutrons: A Transformer Model for Computation Load Estimation on Domain-Decomposed Neutron Transport ProblemsComments: 28 pages, 14 figuresSubjects: Computational Physics (physics.comp-ph); Artificial Intelligence (cs.AI)
Domain decomposition is a technique used to reduce memory overhead on large neutron transport problems. Currently, the optimal load-balanced processor allocation for these domains is typically determined through small-scale simulations of the problem, which can be time-consuming for researchers and must be repeated anytime a problem input is changed. We propose a Transformer model with a unique 3D input embedding, and input representations designed for domain-decomposed neutron transport problems, which can predict the subdomain computation loads generated by small-scale simulations. We demonstrate that such a model trained on domain-decomposed Small Modular Reactor (SMR) simulations achieves 98.2% accuracy while being able to skip the small-scale simulation step entirely. Tests of the model's robustness on variant fuel assemblies, other problem geometries, and changes in simulation parameters are also discussed.
- [115] arXiv:2411.03721 (replaced) [pdf, other]
-
Title: Self-Healing in Dielectric Capacitors: a Universal Method to Computationally Rate Newly Introduced Energy Storage DesignsSubjects: Applied Physics (physics.app-ph); Materials Science (cond-mat.mtrl-sci)
Metal-film dielectric capacitors provide lump portions of energy on demand. While the capacities of various capacitor designs are comparable in magnitude, their stabilities make a difference. Dielectric breakdowns - micro-discharges - routinely occur in capacitors due to the inevitable presence of localized structure defects. The application of polymeric dielectric materials featuring flexible structures helps obtain more uniform insulating layers. At the modern technological level, it is impossible to completely avoid micro-discharges upon device exploitation. Every micro-discharge results in the formation of a soot channel, which is empirically known to exhibit semiconductor behavior. Because of its capability to conduct electricity, the emerged soot channels harm the subsequent capacitor performance and decrease the amount of stored energy. The accumulation of the soot throughout a dielectric capacitor ultimately results in irreversible overall failure. In the context of the dielectric breakdown, self-healing designates a range of chemical processes, which spontaneously rearrange the atoms in the soot channels to partially return their insulative function. We developed a universal method capable of rating new capacitor designs including electrode and polymer material and their proportions. We found the best-performing designs produce abundant volatile by-products after micro-discharge, whereas the soot samples exhibit lower electronic conductivities. We proved the capability of the theoretical method to rate the empirical performance of the known capacitors. The method relies on various electronic-structure simulations and potential landscape explorations. The reported advance opens an impressive avenue to computationally probe thousands of hypothetical capacitor designs.
- [116] arXiv:2411.04350 (replaced) [pdf, html, other]
-
Title: Lagrangian Drifter Path Identification and Prediction: SINDy vs Neural ODESubjects: Atmospheric and Oceanic Physics (physics.ao-ph)
In this study, we investigate the performance of the sparse identification of nonlinear dynamics (SINDy) algorithm and the neural ordinary differential equations (ODEs) in identification of the underlying mechanisms of open ocean Lagrangian drifter hydrodynamics with possible applications in coastal and port hydrodynamic processes. With this motivation we employ two different Lagrangian drifter datasets acquired by National Oceanic and Atmospheric Administration (NOAA)'s surface buoys with proper World Meteorological Organization (WMO) numbers. In the SINDy approach, the primary goal is to identify the drifter paths of buoys using ordinary differential equation sets with a minimal number of sparse coefficients. In the neural ODE approach, the goal is to identify the derivative of the hidden state of a neural network (NN). Using the acquired data, we examine the applicability of the SINDy and the neural ODE algorithms in identification of the drifter trajectories comparatively. We propose that while both of the algorithms may give acceptable results for open ocean, the SINDy-based algorithmic approach can predict the Lagrangian drifter paths more accurately and consistently at least for the datasets investigated and parameters selected. A discussion of our findings with potential applications in search and rescue missions in the open ocean, their limitations and applicability are also presented.
- [117] arXiv:1207.0258 (replaced) [pdf, html, other]
-
Title: Control of probability flow in Markov chain Monte Carlo -- Nonreversibility and liftingComments: 10 pages, 9 figuresJournal-ref: J. Chem. Phys. 161, 174107 (2024)Subjects: Statistical Mechanics (cond-mat.stat-mech); Mathematical Physics (math-ph); Numerical Analysis (math.NA); Data Analysis, Statistics and Probability (physics.data-an); Computation (stat.CO); Methodology (stat.ME)
The Markov chain Monte Carlo (MCMC) method is widely used in various fields as a powerful numerical integration technique for systems with many degrees of freedom. In MCMC methods, probabilistic state transitions can be considered as a random walk in state space, and random walks allow for sampling from complex distributions. However, paradoxically, it is necessary to carefully suppress the randomness of the random walk to improve computational efficiency. By breaking detailed balance, we can create a probability flow in the state space and perform more efficient sampling along this flow. Motivated by this idea, practical and efficient nonreversible MCMC methods have been developed over the past ten years. In particular, the lifting technique, which introduces probability flows in an extended state space, has been applied to various systems and has proven more efficient than conventional reversible updates. We review and discuss several practical approaches to implementing nonreversible MCMC methods, including the shift method in the cumulative distribution and the directed-worm algorithm.
- [118] arXiv:2307.04948 (replaced) [pdf, html, other]
-
Title: Viscous tweezers: controlling particle orientation with viscosityComments: 5 pages, 4 figures + appendicesJournal-ref: Physical Review Research 6, L042039 (2024)Subjects: Soft Condensed Matter (cond-mat.soft); Fluid Dynamics (physics.flu-dyn)
Control of particle motion is generally achieved by applying an external field that acts directly on each particle. Here, we propose a global way to manipulate the motion of a particle by dynamically changing the properties of the fluid in which it is immersed. We exemplify this principle by considering a small particle sinking in an anisotropic fluid whose viscosity depends on the shear axis. In the Stokes regime, the motion of an immersed object is fully determined by the viscosity of the fluid through the mobility matrix, which we explicitly compute for a pushpin-shaped particle. Rather than falling upright under the force of gravity, as in an isotropic fluid, the pushpin tilts to the side, sedimenting at an angle determined by the viscosity anisotropy axis. By changing this axis, we demonstrate control over the pushpin orientation as it sinks, even in the presence of noise, using a closed feedback loop. This strategy to control particle motion, that we dub viscous tweezers, could be experimentally realized in systems ranging from polyatomic fluids under external fields to chiral active fluids of spinning particles by suitably changing their direction of global alignment or anisotropy.
- [119] arXiv:2312.16182 (replaced) [pdf, html, other]
-
Title: Tunneling-like wave transmission in non-Hermitian lattices with mirrored nonreciprocitySubjects: Mesoscale and Nanoscale Physics (cond-mat.mes-hall); Optics (physics.optics)
We report a peculiar tunneling phenomenon that occurs in lattices with nonreciprocal couplings. The nonreciprocity holds for an inner portion of the lattice, constituting a non-Hermitian interface between outer Hermitian sections. The couplings are mirrored about the interface center. As a standalone system that was widely studied in recent years, each section of the interface supports the non-Hermitian skin effect, in which modes are accumulated at one boundary. Here, we investigate what happens to a wave that propagates along the lattice and hits the interface. The skin mode accumulation, which effectively constitutes a barrier, forbids wave penetration into the interface, but surprisingly, under certain conditions the wave is transmitted to the other side, keeping the interface dark, as if the wave invisibly tunneled through it. Remarkably, the tunneling is independent of the interface length, and a perfect transmission can be achieved independently of frequency and nonreciprocity strength. We derive the phenomenon both for quantum and classical systems, and realize it experimentally in an active topoelectric metamaterial. Our study fosters the research of wave tunneling through other types of non-Hermitian interfaces, which may also include nonlinearities, time-dependence and more.
- [120] arXiv:2403.12018 (replaced) [pdf, html, other]
-
Title: Exceptional points of any order in a generalized Hatano-Nelson modelSubjects: Mesoscale and Nanoscale Physics (cond-mat.mes-hall); Optics (physics.optics); Quantum Physics (quant-ph)
Exceptional points (EPs) are truly non-Hermitian (NH) degeneracies where matrices become defective. The order of such an EP is given by the number of coalescing eigenvectors. On the one hand, most work focuses on studying $N$th-order EPs in $(N\leq4)$-dimensional NH Bloch Hamiltonians. On the other hand, some works have remarked on the existence of EPs of orders scaling with systems size in models exhibiting the NH skin effect. In this work, we introduce a new type of EP and provide a recipe on how to realize EPs of arbitrary order not scaling with system size. We introduce a generalized version of the paradigmatic Hatano-Nelson model with longer-range hoppings. The EPs existing in this system show remarkable physical features: Their associated eigenstates have support on a subset of sites and exhibit the NH skin effect, which can be tuned to localize on the opposite end of the chain compared to all remaining states. Furthermore, the EPs are robust against generic perturbations in the hopping strengths as well as against a specific form of on-site disorder.
- [121] arXiv:2403.18088 (replaced) [pdf, html, other]
-
Title: Discretize first, filter next: learning divergence-consistent closure models for large-eddy simulationComments: 53 pagesSubjects: Numerical Analysis (math.NA); Fluid Dynamics (physics.flu-dyn)
We propose a new neural network based large eddy simulation framework for the incompressible Navier-Stokes equations based on the paradigm "discretize first, filter and close next". This leads to full model-data consistency and allows for employing neural closure models in the same environment as where they have been trained. Since the LES discretization error is included in the learning process, the closure models can learn to account for the discretization. Furthermore, we employ a divergence-consistent discrete filter defined through face-averaging and provide novel theoretical and numerical filter analysis. This filter preserves the discrete divergence-free constraint by construction, unlike general discrete filters such as volume-averaging filters. We show that using a divergence-consistent LES formulation coupled with a convolutional neural closure model produces stable and accurate results for both a-priori and a-posteriori training, while a general (divergence-inconsistent) LES model requires a-posteriori training or other stability-enforcing measures.
- [122] arXiv:2405.11414 (replaced) [pdf, html, other]
-
Title: High-Resolution Agent-Based Modeling of Campus Population Behaviors for Pandemic Response PlanningComments: 15 pages, 6 figures; to be published in the Proceedings of PPAM 2024Subjects: Computers and Society (cs.CY); Computational Engineering, Finance, and Science (cs.CE); Multiagent Systems (cs.MA); Physics and Society (physics.soc-ph)
This paper reports a case study of an application of high-resolution agent-based modeling and simulation to pandemic response planning on a university campus. In the summer of 2020, we were tasked with a COVID-19 pandemic response project to create a detailed behavioral simulation model of the entire campus population at Binghamton University. We conceptualized this problem as an agent migration process on a multilayer transportation network, in which each layer represented a different transportation mode. As no direct data were available about people's behaviors on campus, we collected as much indirect information as possible to inform the agents' behavioral rules. Each agent was assumed to move along the shortest path between two locations within each transportation layer and switch layers at a parking lot or a bus stop, along with several other behavioral assumptions. Using this model, we conducted simulations of the whole campus population behaviors on a typical weekday, involving more than 25,000 agents. We measured the frequency of close social contacts at each spatial location and identified several busy locations and corridors on campus that needed substantial behavioral intervention. Moreover, systematic simulations with varying population density revealed that the effect of population density reduction was nonlinear, and that reducing the population density to 40-45% would be optimal and sufficient to suppress disease spreading on campus. These results were reported to the university administration and utilized in the pandemic response planning, which led to successful outcomes.
- [123] arXiv:2405.12022 (replaced) [pdf, html, other]
-
Title: Lattice physics approaches for neural networksComments: 16 pages, 2 figuresSubjects: Neurons and Cognition (q-bio.NC); Applied Physics (physics.app-ph)
Modern neuroscience has evolved into a frontier field that draws on numerous disciplines, resulting in the flourishing of novel conceptual frames primarily inspired by physics and complex systems science. Contributing in this direction, we recently introduced a mathematical framework to describe the spatiotemporal interactions of systems of neurons using lattice field theory, the reference paradigm for theoretical particle physics. In this note, we provide a concise summary of the basics of the theory, aiming to be intuitive to the interdisciplinary neuroscience community. We contextualize our methods, illustrating how to readily connect the parameters of our formulation to experimental variables using well-known renormalization procedures. This synopsis yields the key concepts needed to describe neural networks using lattice physics. Such classes of methods are attention-worthy in an era of blistering improvements in numerical computations, as they can facilitate relating the observation of neural activity to generative models underpinned by physical principles.
- [124] arXiv:2406.01005 (replaced) [pdf, html, other]
-
Title: Transfer Matrix Approach for Topological Edge StatesComments: 12 pages, 9 figuresSubjects: Mesoscale and Nanoscale Physics (cond-mat.mes-hall); Optics (physics.optics)
We suggest and develop a novel approach for describing topological properties of a periodic system purely from the transfer matrix associated to a unit cell. Our approach uses the Iwasawa decomposition to parametrise the transfer matrix uniquely in terms of three real numbers. This allows us to obtain simple conditions for the existence of topologically protected edge states and to provide a visual illustration of all possible solutions. In order to demonstrate our method in action, we apply it to study some generalisations of the Su-Schrieffer-Heeger (SSH) model, such as the tetramer SSH4 model and a dimerised one-dimensional photonic crystal. Finally, we also obtained a simple pictorial proof of the Zak phase bulk-edge correspondence for any one dimensional system using this approach.
- [125] arXiv:2407.20572 (replaced) [pdf, html, other]
-
Title: Scale locality of information flow in turbulenceComments: 16 pages, 5 figures. Minor revisions are made in ver. 2Subjects: Statistical Mechanics (cond-mat.stat-mech); Fluid Dynamics (physics.flu-dyn)
Turbulent fluctuations exhibit universal scaling laws that are independent of large-scale statistics. It is often explained that such universality is caused by the loss of information about large-scale statistics during the cascade process. In our previous study [T. Tanogami and R. Araki, Phys. Rev. Research 6, 013090 (2024)], we applied information thermodynamics to turbulence and proved that information of large-scale turbulent fluctuations is propagated to small scales. As a first step toward understanding how universality emerges at small scales under the influence of the information flow from large scales, here we investigate the scale locality of the information flow for the Sabra shell model. First, we analytically show that the information flow can be decomposed into scale-local and scale-nonlocal parts. Then, we prove that the scale-nonlocal part can be ignored compared to the scale-local part. This result implies that the information transfer from large to small scales occurs mainly through scale-local interactions, which is consistent with the scale locality of the energy cascade.
- [126] arXiv:2408.01823 (replaced) [pdf, html, other]
-
Title: Taming Uncertainty in a Complex World: The Rise of Uncertainty Quantification -- A Tutorial for BeginnersComments: 33 pages, 15 figures, 1 table. To appear in the Notices of the American Mathematical Society. The code to the examples in the supplementary document can be found in this https URLSubjects: Dynamical Systems (math.DS); Atmospheric and Oceanic Physics (physics.ao-ph); Data Analysis, Statistics and Probability (physics.data-an)
This paper provides a tutorial about uncertainty quantification (UQ) for those who have no background but are interested in learning more in this area. It exploits many very simple examples, which are understandable to undergraduates, to present the ideas of UQ. Topics include characterizing uncertainties using information theory, UQ in linear and nonlinear dynamical systems, UQ via data assimilation, the role of uncertainty in diagnostics, and UQ in advancing efficient modeling. The surprisingly simple examples in each topic explain why and how UQ is essential. Both MATLAB and Python codes are made available for these simple examples.
- [127] arXiv:2408.03635 (replaced) [pdf, html, other]
-
Title: Scale-to-Scale Information Flow Amplifies Turbulent FluctuationsComments: 7 pages, 2 figures (+ supplemental material 8 pages, 2 figures). arXiv admin note: text overlap with arXiv:2407.20572. Minor revisions are made in ver. 2Subjects: Statistical Mechanics (cond-mat.stat-mech); Fluid Dynamics (physics.flu-dyn)
In three-dimensional turbulence, information of turbulent fluctuations at large scales is propagated to small scales. Here, we investigate the relation between the information flow and turbulent fluctuations described by a shell model. We first establish a connection between the information flow and phase-space contraction rate. From this relation, we then prove an inequality between the information flow and turbulent fluctuations, which suggests that the information flow from large to small scales amplifies turbulent fluctuations at small scales. This inequality can also be interpreted as a quantification of Landau's objection to the universality of turbulent fluctuations. We also discuss differences between the information flow and the Kolmogorov-Sinai entropy.
- [128] arXiv:2408.06093 (replaced) [pdf, html, other]
-
Title: Effective and efficient modeling of the hydrodynamics for bacterial flagellaSubjects: Soft Condensed Matter (cond-mat.soft); Biological Physics (physics.bio-ph)
The hydrodynamic interactions among bacterial cell bodies, flagella, and surrounding boundaries are essential for understanding bacterial motility in complex environments. In this study, we demonstrate that each slender flagellum can be modeled as a series of spheres, and that the interactions between these spheres can be accurately characterized using a resistance matrix. This approach allows us to effectively and efficiently evaluate the propulsive effects of the flagella. Notably, our investigation into bacterial motility near a colloidal sphere reveals significant discrepancies between results derived from the twin multipole moment and those obtained through resistive force theory. Consequently, neglecting the hydrodynamic interactions among cell bodies, flagella, and colloidal spheres may lead to substantial inaccuracies. Our model simplifies bacteria into a series of spheres, making it well-suited for examining bacterial motility near spherical boundaries, as well as the nonlinear deformation dynamics of elastic flagella.
- [129] arXiv:2408.11427 (replaced) [pdf, html, other]
-
Title: Role of momentum in the generator-coordinate method applied to barrier penetrationComments: 7 pages, 4 figures. To appear in Phys. Rev. CSubjects: Nuclear Theory (nucl-th); Atomic Physics (physics.atom-ph); Chemical Physics (physics.chem-ph); Quantum Physics (quant-ph)
Nuclear fission at barrier-top energies is conventionally modeled by a one-dimensional Schrödinger equation applied to internal fission channels, but that treatment is hard to justify in the configuration-interaction approach to nuclear Hamiltonians. Here we show that inclusion of states of finite momentum by the Generator Coordinate Method (GCM) considerably extends the range of energies at which GCM-based Hamiltonians could reproduce the Schrödinger treatment. The transmission probabilities for crossing the barrier are calculated by a discrete version of Kohn's variational method, which may also be useful for other systems of interacting fermions.
- [130] arXiv:2410.12873 (replaced) [pdf, html, other]
-
Title: The exact travelling wave solutions of a KPP equationComments: pdfLaTeX, 5 pages. The paper was additionally editedSubjects: Populations and Evolution (q-bio.PE); Mathematical Physics (math-ph); Biological Physics (physics.bio-ph)
We obtain the exact analytical traveling wave solutions of the Kolmogorov-Petrovskii-Piskunov equation with the reaction term belonging to the class of functions, which includes that of the (generalized) Fisher equation, for the particular values of the waves speed. Additionally we obtain the the exact analytical traveling wave solutions of the generalized Burgers-Huxley equation.
- [131] arXiv:2410.21046 (replaced) [pdf, html, other]
-
Title: Deep Learning the Forecast of Galactic Cosmic-Rays SpectraComments: 8 pages, 7 figures, 1 tableSubjects: High Energy Astrophysical Phenomena (astro-ph.HE); Solar and Stellar Astrophysics (astro-ph.SR); High Energy Physics - Phenomenology (hep-ph); Space Physics (physics.space-ph)
We introduce a novel deep learning framework based on Long Short-Term Memory (LSTM) networks to predict galactic cosmic-ray spectra on a one-day-ahead basis by leveraging historical solar activity data, overcoming limitations inherent in traditional transport models. By flexibly incorporating multiple solar parameters such as the heliospheric magnetic field, solar wind speed, and sunspot numbers, our model achieves accurate short-term and long-term predictions of cosmic-ray flux. The addition of historical cosmic-ray flux data significantly enhances prediction accuracy, allowing the model to capture complex dependencies between past and future flux variations. Additionally, the model reliably predicts full cosmic-ray spectra for different particle species, enhancing its utility for comprehensive space weather forecasting. Our approach offers a scalable, data-driven alternative to traditional physics-based methods, ensuring robust daily and long-term forecasts. This work opens avenues for advanced models that can integrate broader observational data, with significant implications for space weather monitoring and mission planning.
- [132] arXiv:2410.23893 (replaced) [pdf, html, other]
-
Title: DiffBatt: A Diffusion Model for Battery Degradation Prediction and SynthesisHamidreza Eivazi, André Hebenbrock, Raphael Ginster, Steffen Blömeke, Stefan Wittek, Christoph Herrmann, Thomas S. Spengler, Thomas Turek, Andreas RauschComments: 15 pages, 6 figuresJournal-ref: Foundation Models for Science Workshop, 38th Conference on Neural Information Processing Systems (NeurIPS 2024)Subjects: Machine Learning (cs.LG); Chemical Physics (physics.chem-ph)
Battery degradation remains a critical challenge in the pursuit of green technologies and sustainable energy solutions. Despite significant research efforts, predicting battery capacity loss accurately remains a formidable task due to its complex nature, influenced by both aging and cycling behaviors. To address this challenge, we introduce a novel general-purpose model for battery degradation prediction and synthesis, DiffBatt. Leveraging an innovative combination of conditional and unconditional diffusion models with classifier-free guidance and transformer architecture, DiffBatt achieves high expressivity and scalability. DiffBatt operates as a probabilistic model to capture uncertainty in aging behaviors and a generative model to simulate battery degradation. The performance of the model excels in prediction tasks while also enabling the generation of synthetic degradation curves, facilitating enhanced model training by data augmentation. In the remaining useful life prediction task, DiffBatt provides accurate results with a mean RMSE of 196 cycles across all datasets, outperforming all other models and demonstrating superior generalizability. This work represents an important step towards developing foundational models for battery degradation.
- [133] arXiv:2411.01071 (replaced) [pdf, other]
-
Title: Towards Nanoscale and Element-Specific Lattice Temperature Measurements using Core-Loss Electron Energy-Loss SpectroscopyComments: 33 pages, 5 main text figures, 1 tableSubjects: Materials Science (cond-mat.mtrl-sci); Chemical Physics (physics.chem-ph)
Measuring nanoscale local temperatures, particularly in vertically integrated and multi-component systems, remains challenging. Spectroscopic techniques like X-ray absorption and core-loss electron energy-loss spectroscopy (EELS) are sensitive to lattice temperature, but understanding thermal effects is nontrivial. This work explores the potential for nanoscale and element-specific core-loss thermometry by comparing the Si L2,3 edge's temperature-dependent redshift against plasmon energy expansion thermometry (PEET) in a scanning TEM. Using density functional theory (DFT), time-dependent DFT, and the Bethe-Salpeter equation, we ab initio model both the Si L2,3 and plasmon redshift. We find that the core-loss redshift is due to bandgap reduction from electron-phonon renormalization. Our results indicate that despite lower core-loss signal intensity and thus accuracy compared to PEET, core-loss thermometry still has important advantages. Specifically, we show that the Varshni equation easily interprets the core-loss redshift, which avoids plasmon spectral convolution for PEET in complex materials. We also find that core-loss thermometry is more accurate than PEET at modeling thermal lattice expansion unless the temperature-dependent effective mass and dielectric constant are known. Furthermore, core-loss thermometry has the potential to measure nanoscale heating in multi-component materials and stacked interfaces with elemental specificity at length scales smaller than the plasmon's wavefunction.
- [134] arXiv:2411.03467 (replaced) [pdf, html, other]
-
Title: A high resolution simulation of protoplanetary disk turbulence driven by the vertical shear instabilitySubjects: Solar and Stellar Astrophysics (astro-ph.SR); Earth and Planetary Astrophysics (astro-ph.EP); Fluid Dynamics (physics.flu-dyn)
A high resolution fourth-order Padé scheme is used to simulate locally isothermal 3D disk turbulence driven by the vertical shear instability (VSI) using 268.4 M points. In the early non-linear period of axisymmetric VSI, angular momentum transport by vertical jets creates correlated N-shaped radial profiles of perturbation vertical and azimuthal velocity. This implies dominance of positive perturbation vertical vorticity layers and a recently discovered angular momentum staircase with respect to radius ($r$). These features are present in 3D in a weaker form. The 3D flow consists of vertically and azimuthally coherent turbulent shear layers containing small vortices with all three vorticity components active. Previously observed large persistent vortices in the interior of the domain driven by the Rossby wave instability are absent. We speculate that this is due to a weaker angular momentum staircase in 3D in the present simulations compared to a previous simulation. The turbulent viscosity parameter $\alpha(r)$ increases linearly with $r$. At intermediate resolution, the value of $\alpha(r)$ at midradius is close to that of a previous simulation. The specific kinetic energy spectrum with respect to radial wavenumber has a power law region with exponent $-1.84$, close to the value $-2$ expected for shear layers. The spectrum with respect to azimuthal wavenumber has a $-5/3$ region and lacks a $-5$ region reported in an earlier study. Finally, it is found that axisymmetric VSI has artifacts at late times, including a very strong angular momentum staircase, which in 3D is present weakly in the disk's upper layers.