Physics

2025-02-14 | | Total: 73

#1 Shadow Molecular Dynamics for a Charge-Potential Equilibration Model [PDF] [Copy] [Kimi] [REL]

Authors: Robert Stanton, Cagri Mehmet Kaymak, Anders M. N. Niklasson

We introduce a shadow molecular dynamics (MD) approach based on the Atom-Condensed Kohn-Sham second-order (ACKS2) charge-potential equilibration model. In contrast to regular flexible charge models, the ACKS2 model includes both flexible atomic charges and potential fluctuation parameters that allow for physically correct charge fragmentation and scaling of the polarizability. Our shadow MD scheme is based on an approximation of the ACKS2's flexible charge-potential energy function, in combination with extended Lagrangian Born-Oppenheimer MD. Utilizing this shadow charge-potential equilibration approach mitigates the costly overhead and stability problems associated with finding well-converged iterative solutions to the charges and potential fluctuations of the ACKS2 model in an MD simulation. Our work provides a robust and versatile framework for efficient, high-fidelity MD simulations of diverse physical phenomena and applications.

Subjects: Chemical Physics , Materials Science

Publish: 2025-02-13 18:53:22 UTC


#2 Second-order theory for multi-hinged directional wavemakers [PDF] [Copy] [Kimi] [REL]

Author: Andreas H. Akselsen

The second-order directional wavemaker theory for regular and irregular waves is extended to multi-hinged wavemakers and combined piston--flap wavemaker systems. Derived expressions enable second-order signal correction, common in single-hinged wavemakers, to be applied to multi-hinged systems. Multi-hinged wavemakers offer additional degrees of freedom, with different combinations of paddle motion producing the same progressive wave. This is here exploited to better understand wavemaker behaviour. Single-harmonic signals are computed for double-hinged wavemakers that suppress spurious waves without introducing double-harmonic motions. Surprisingly, these flap motions are almost always in opposite phase, with the larger draft found underneath the water surface. %contradicting to assumptions from earlier studies. Due to the opposing paddle phase, the double-hinged wavemaker draft is usually smaller than the corresponding single-hinged draft. The ability of thsee systems to suppress spurious waves with single-harmonic motion is verified experimentally. The wavemaker theory further supports an arbitrary number of flap hinges, enabling the approximation of a fully flexible wavemaker through piecewise-linear segments. This is demonstrated with an exponential wavemaker profile that does not generate any evanescent waves at linear order. Such a wavemaker is likely to limit wave breaking and cross-modes, but is found to produce spurious second-order waves of a magnitude comparable to a single flap. The presented solution is complete, intrinsically including return flow through the second-order zero mode. This return flow is found to precisely match the Stokes drift.

Subject: Fluid Dynamics

Publish: 2025-02-13 18:44:49 UTC


#3 Wireless and passive pressure detection using magneto-mechanical resonances in process engineering [PDF] [Copy] [Kimi] [REL]

Authors: Timo Merbach, Felix Kexel, Jonas Faltinath, Martin Möddel, Michael Schlüter, Tobias Knopp, Fabian Mohn

A custom-developed magneto-mechanical resonator (MMR) for wireless pressure measurement is investigated for potential applications in process engineering. The MMR sensor utilises changes in the resonance frequency caused by pressure on a flexible 3D printed membrane. The thickness of the printed membrane plays a crucial role in determining the performance and sensitivity of MMRs, and can be tailored to meet the requirements of specific applications. The study includes static and dynamic measurements to determine the pressure sensitivity and temporal resolution of the sensor. The results show a minimum sensitivity of $0.06~\text{Hz mbar}^{-1}$ and are in agreement with theoretical calculations and measurements. The maximum sensor readout frequency is $2~\text{Hz}$ in this study. Additionally, the temperature dependence of the sensor is investigated, revealing a significant dependence of the resonance frequency on temperature. The developed MMR offers a promising and versatile method for precise pressure measurements in process engineering environments.

Subject: Applied Physics

Publish: 2025-02-13 18:32:46 UTC


#4 Nonmodal amplitude equations [PDF] [Copy] [Kimi] [REL]

Authors: Yves-Marie Ducimetière, François Gallaire

We consider fluid flows for which the linearized Navier-Stokes operator is strongly non-normal. The responses of such flows to external perturbations are spanned by a generically very large number of (non-orthogonal) eigenmodes. They are therefore qualified as ``nonmodal" responses, to insist on the inefficiency of the eigenbasis to describe them. For the purpose of this article to reduce the dimension of the system in a weakly nonlinear regime, (eigen)modal reduction techniques, such as the center manifold, are thus inappropriate precisely because the leading-order dynamics cannot be restricted to a low-dimensional eigensubspace. On the other hand, it is often true that only a small number (we assume one) of singular modes is sufficient to reconstruct the nonmodal responses at the leading order. By adopting the latter paradigm, we propose a method to analytically derive a weakly nonlinear amplitude equation for the nonmodal response of a fluid flow to a small harmonic forcing, stochastic forcing, and initial perturbation, respectively. The present approach is simpler than the one we previously proposed, for neither the operator perturbation nor the ensuing compatibility condition were formally necessary. When applied to the stochastic response, the present method makes it possible to derive an amplitude equation that is substantially easier to solve than the one we proposed previously. Despite being concerned with the response to external disturbances of three different natures, the derivations of the amplitude equations all proceed by the same general principle. Eventually, the three derived amplitude equations are tested in three distinct flows. At extremely low numerical cost as compared to fully nonlinear techniques, they can predict the weakly nonlinear modification of the gains as the flow departs from the linear regime by increasing the amplitude of the external excitation.

Subjects: Fluid Dynamics , Mathematical Physics

Publish: 2025-02-13 18:30:18 UTC


#5 Cryogenic SiPMs for the Optical Readout of DarkSide-20k [PDF] [Copy] [Kimi] [REL]

Author: Giuseppe Matteucci

Silicon photomultipliers (SiPM) have gained significant traction as an alternative technology to the well-established photomultiplier tube (PMT), with numerous high-sensitivity experiments adopting them either complementarily or as a replacement for PMTs. SiPMs are an ideal match for low-background cryogenic applications, such as massive noble liquid experiments for dark matter direct detection, due to (i) the significant reduction of dark noise in cold environments, (ii) relatively low radioactive content, and (iii) scalable industrial production. For these reasons, the Global Argon Dark Matter Collaboration has committed to this technology for DarkSide-20k, its next experiment for the direct search of WIMP Dark Matter, currently in construction at LNGS Hall C. The development of a large-area cryogenic SiPM-based photon counter has culminated in the Photon Detector Unit (PDU), a compact photosensor measuring $20\times20$ cm$^2$ with $100$ cm$^2$ active surface per channel, based on SiPM technology from Fondazione Bruno Kessler and incorporating custom front-end electronics suited for cryogenics. More than 600 PDUs are being produced and tested in various collaboration facilities to construct the two $~10.5$ m$^2$ optical planes of the massive two-phase argon time projection chamber of DarkSide-20k and the optical readout of its veto system.

Subjects: Instrumentation and Detectors , High Energy Physics - Experiment

Publish: 2025-02-13 18:10:40 UTC


#6 Correlation based modeling of the ionospheric magnetic field [PDF] [Copy] [Kimi] [REL]

Authors: K. Ferrat, M. Schanner, V. Lesur, M. Holschneider, C. Stolle, J. Baerenzung

Patterns of the magnetic signature of ionospheric currents, generated from an empirical model based on satellite observations, are used to build a statistical correlation based model for ionospheric fields. In order to stabilize the dynamics and to take into account the dominant role of the sun, the fields are represented in solar magnetic coordinates. The covariance structure is analyzed and a second order process that approximates the full dynamics is generated. We show that for synthetic data observations located at the Earth observatories, the full ionospheric field pattern as observed on the earth surface can be reconstructed with good precision. As a proof of principle we provide a first application to the inversion based on real observatory data.

Subject: Geophysics

Publish: 2025-02-13 16:59:16 UTC


#7 Structural phase transitions between layered Indium Selenide for inte-grated photonic memory [PDF] [Copy] [Kimi] [REL]

Authors: Tiantian Li, Yong Wang, Wei Li, Dun Mao, Chris J. Benmore, Igor Evangelista, Huadan Xing, Qiu Li, Feifan Wang, Ganesh Sivaraman, Anderson Janotti, Stephanie Law, Tingyi Gu

The primary mechanism of optical memristive devices relies on the phase transitions between amorphous-crystalline states. The slow or energy hungry amorphous-crystalline transitions in optical phase-change materials are detrimental to the devices scalability and performance. Leveraging the integrated photonic platform, we demonstrate a single nanosecond pulse triggered nonvolatile and reversible switching between two layered structures of indium selenide (In2Se3). High resolution pair distribution function reveals the detailed atomistic transition pathways between the layered structures. With inter-layer shear glide and isosymmetric phase transition, the switching between alpha and beta structural states contain low re-configurational entropy, allowing reversible switching between layered structures. Broadband refractive index contrast, optical transparency, and volumetric effect in the crystalline-crystalline phase transition are experimentally characterized in molecular beam epitaxy-grown thin films and compared to ab initials calculations. The nonlinear resonator transmission spectra measure an incremental linear loss rate of 3.3 GHz introduced by 1.5 micrometer long In2Se3 covered lay-er, resulting from the combinations of material absorption and scattering.

Subject: Optics

Publish: 2025-02-13 16:38:40 UTC


#8 Recoil nuclear size corrections in hydrogenic systems [PDF] [Copy] [Kimi] [REL]

Authors: Krzysztof Pachucki, Vojtěch Patkóš, Vladimir A. Yerokhin

Formulas for the combined nuclear-recoil and finite-nuclear-size effects of order $(Z\,\alpha)^5$ and $(Z\,\alpha)^6$ are derived without any expansion in the nuclear charge radius $r_C$, making them applicable to both electronic and muonic atoms. The obtained results are particularly relevant for high-precision determinations of root-mean-square charge radii from muonic atom spectroscopy. We demonstrate that calculations of the atomic isotope shift based on the widely used Breit approximation give rise to an unphysical nuclear-size contribution that is linear in the nuclear charge radius $r_C$ at order $(Z\,\alpha)^5$. This spurious term vanishes in a full QED treatment, leaving the correct contribution quadratic in $r_C$. For electronic atoms, this quadratic term is significantly smaller than the spurious linear contribution.

Subject: Atomic Physics

Publish: 2025-02-13 16:21:24 UTC


#9 Drivers of cooperation in social dilemmas on higher-order networks [PDF] [Copy] [Kimi] [REL]

Authors: Onkar Sadekar, Andrea Civilini, Vito Latora, Federico Battiston

Understanding cooperation in social dilemmas requires models that capture the complexity of real-world interactions. While network frameworks have provided valuable insights to model the evolution of cooperation, they are unable to encode group interactions properly. Here, we introduce a general higher-order network framework for multi-player games on structured populations. Our model considers multi-dimensional strategies, based on the observation that social behaviours are affected by the size of the group interaction. We investigate dynamical and structural coupling between different orders of interactions, revealing the crucial role of nested multilevel interactions, and showing how such features can enhance cooperation beyond the limit of traditional models with uni-dimensional strategies. Our work identifies the key drivers promoting cooperative behaviour commonly observed in real-world group social dilemmas.

Subjects: Physics and Society , Computer Science and Game Theory , Social and Information Networks , Populations and Evolution

Publish: 2025-02-13 16:15:45 UTC


#10 The Third Generation of Nanogenerators: The Irreplaceable Potential Source Enabled by the Flexoelectric Nanogenerator [PDF] [Copy] [Kimi] [REL]

Authors: Shang Ru Li, Qi Kang Zhang, Xiao Xiong Wang

The electroneutrality assumption has long been adopted by scholars; however, this assumption may lead to an oversight of certain physical effects. Using derivations from a discontinuous medium, we have obtained an expression for the potential and energy of a many-body unipolar charge system, which corresponds well to its counterpart in a continuous medium. The compressed form of this expression suggests that compressing a macroscale charged body to the nanoscale can yield an enormous electric potential and energy, thereby establishing a concrete research framework for third-generation nanogenerators. This effect may serve as a crucial reference for understanding anomalous spatial electromagnetic distributions and divergent energy fields.

Subject: Applied Physics

Publish: 2025-02-13 16:04:28 UTC


#11 Measuring and controlling the birth of quantum attosecond pulses [PDF] [Copy] [Kimi] [REL]

Authors: Matan Even Tzur, Chen Mor, Noa Yaffe, Michael Birk, Andrei Rasputnyi, Omer Kneller, Ido Nisim, Ido Kaminer, Michael Krüger, Nirit Dudovich, Oren Cohen

The generation and control of extreme ultraviolet (XUV) radiation by high harmonic generation (HHG) have advanced ultrafast science, providing direct insights into electron dynamics on their natural time scale. Attosecond science has established the capability to resolve ultrafast quantum phenomena in matter by characterizing and controlling the classical properties of the high harmonics. Recent theoretical proposals have introduced novel schemes for generating and manipulating XUV HHG with distinct quantum features, paving the way to attosecond quantum optics. In this work, we transfer fundamental concepts in quantum optics into attosecond science. By driving the HHG process with a combination of an infrared bright squeezed vacuum (BSV, a non-classical state of light), and a strong coherent field, we imprint the quantum correlations of the input BSV onto both the ultrafast electron wavefunction and the harmonics' field. Performing in-situ HHG interferometry provides an insight into the underlying sub-cycle dynamics, revealing squeezing in the statistical properties of one of the most fundamental strong-field phenomena -- field induced tunneling. Our measurement allows the reconstruction of the quantum state of the harmonics through homodyne-like tomography, resolving correlated fluctuations in the harmonic field that mirror those of the input BSV. By controlling the delay between the two driving fields, we manipulate the photon statistics of the emitted attosecond pulses with sub-cycle accuracy. The ability to measure and control quantum correlations in both electrons and XUV attosecond pulses establishes a foundation for attosecond electrodynamics, manipulating the quantum state of electrons and photons with sub-cycle precision.

Subject: Optics

Publish: 2025-02-13 15:48:43 UTC


#12 Convex envelope method for T, p flash calculations for chemical systems with an arbitrary number of components and arbitrary aggregate states [PDF] [Copy] [Kimi] [REL]

Authors: Quirin Göttl, Natalie Rosen, Jakob Burger

T, p flash calculations determine the correct number of phases at phase equilibrium and their compositions for fixed temperature and pressure. They are essential for chemical process simulation and optimization. The convex envelope method (CEM) is an existing approach that employs the tangent plane criterion to determine liquid phase equilibria for chemical systems with an arbitrary number of components without providing the number of phases beforehand. This work extends the CEM to include also vapor and solid phases. Thus, any phase equilibrium of a given system with an arbitrary number of components and phases can be calculated over the whole composition space. The CEM results are presented for various vapor-liquid and solid-liquid phase equilibria examples of up to four components. We show how the CEM can be used for parameter fitting of gE-models. As an outlook, we demonstrate how the CEM can be combined with a machine learning-based tool for property prediction to construct phase equilibria.

Subject: Chemical Physics

Publish: 2025-02-13 15:27:36 UTC


#13 Nondestructive quantitative estimation of cross-sectional corrosion degree of rebar using self-magnetic flux leakage field variation [PDF] [Copy] [Kimi] [REL]

Authors: Junli Qiu, Weiping Zhang, Chao Jiang, Hong Zhang, Jianting Zhou

To accurately assess the structural performance of corroded reinforced concrete structures, obtaining precise information on the corrosion range and corrosion degree of rebar is crucial. In this study, based on a comprehensive analysis of extensive high-precision magnetic field and three-dimensional structural light scanning data of 21 corroded rebars, it was found that the self-magnetic flux leakage can accurately identify the corrosion range, with an error not exceeding 3%. A proposed quantitative index NHx of the self-magnetic flux leakage amplitude exhibits a linear correlation with the cross-sectional corrosion degree of rebar, whose probability density distribution can be accurately described using the Weibull distribution function. Utilizing the Weibull distribution function of NHx and a Bayesian model, automatically rapid quantification of the rebar's cross-sectional corrosion degree based on the non-destructive testing-derived NHx values can be conveniently realized. This self-magnetic flux leakage-based novel method for quantifying rebar's cross-sectional corrosion degree is accurate, efficient, and well-suited for practical engineering applications, providing robust support for a precise assessment of the structural performance of corroded reinforced concrete structures.

Subject: Applied Physics

Publish: 2025-02-13 15:24:08 UTC


#14 Deterministic fabrication of GaAs-quantum-dot micropillar single-photon sources [PDF] [Copy] [Kimi] [REL]

Authors: Abdulmalik A. Madigawa, Martin Arentoft Jacobsen, Claudia Piccinini, Paweł Wyborski, Ailton Garcia Jr., Saimon F. Covre da Silva, Armando Rastelli, Battulga Munkhbat, Niels Gregersen

This study investigates the performance of droplet-etched GaAs quantum dots (QDs) integrated into micropillar structures using a deterministic fabrication technique. We demonstrate a unity QD positioning yield across 74 devices and consistent device performance. Under p-shell excitation, the QD decay dynamics within the micropillars exhibit biexponential behavior, accompanied by intensity fluctuations limiting the source efficiency to < 4.5%. Charge stabilization via low-power above-band LED excitation effectively reduces these fluctuations, doubling the source efficiency to $\sim$ 9%. Moreover, we introduce suppression of radiation modes by introducing cylindrical rings theoretically predicted to boost the collection efficiency by a factor of 4. Experimentally, only a modest improvement is obtained, underscoring the influence of even minor fabrication imperfections for this advanced design. Our findings demonstrate the reliability of our deterministic fabrication approach in producing high-yield, uniform devices, while offering detailed insights into the influence of charge noise and complex relaxation dynamics on the performance.

Subjects: Optics , Quantum Physics

Publish: 2025-02-13 14:33:40 UTC


#15 Modelling of shattered pellet injection experiments on the ASDEX Upgrade tokamak [PDF] [Copy] [Kimi] [REL]

Authors: Ansh Patel, Akinobu Matsuyama, Gergely Papp, Michael Lehnen, J Artola, Stefan Jachmich, Emiliano Fable, Alexander Bock, Bernd Kurzan, Matthias Hölzl, Weikang Tang, Michael Dunne, Rainer Fischer, Paul Heinrich, The ASDEX Upgrade Team, The EUROfusion Tokamak Exploitation Team

In a shattered pellet injection (SPI) system the penetration and assimilation of the injected material depends on the speed and size distribution of the SPI fragments. ASDEX Upgrade (AUG) was recently equipped with a flexible SPI to study the effect of these parameters on disruption mitigation efficiency. In this paper we study the impact of different parameters on SPI assimilation with the 1.5D INDEX code. Scans of fragment sizes, speeds and different pellet compositions are carried out for single SPI into AUG H-mode plasmas. We use a semi-empirical thermal quench (TQ) onset condition to study the material assimilation trends. For mixed deuterium-neon pellets, smaller/faster fragments start to assimilate quicker. However, at the expected onset of the global reconnection event (GRE),larger/faster fragments end up assimilating more material. Variations in the injected neon content lead to a large difference in the assimilated neon for neon content below $< 10^{21}$ atoms. For larger injected neon content, a self-regulating mechanism limits the variation in the amount of assimilated neon. We use a back-averaging model to simulate the plasmoid drift during pure deuterium injections with the back-averaging parameter determined by a interpretative simulation of an experimental pure deuterium injection discharge. Again, larger and faster fragments are found to lead to higher assimilation with the material assimilation limited to the plasma edge in general, due to the plasmoid drift. The trends of assimilation for varying fragment sizes, speeds and pellet composition qualitatively agree with the previously reported experimental observations.

Subject: Plasma Physics

Publish: 2025-02-13 14:22:20 UTC


#16 Optimizing omnigenity like quasisymmetry for stellarators [PDF] [Copy] [Kimi] [REL]

Authors: Hengqian Liu, Guodong Yu, Caoxiang Zhu, Ge Zhuang

Omnigenous magnetic fields, where the bounce-averaged radial drift vanishes, offer a promising solution to confine charged particles in fusion devices, particularly in stellarators. However, nonquasisymmetric omnigenity has remained underexplored due to the absence of a general optimization method. We introduce a novel approach for optimizing omnigenity. With simplicity comparable to quasisymmetry (QS) optimization, the new method unifies both QS and non-QS omnigenity optimization and can be further generalized to optimize configurations beyond omnigenity. Precisely omnigenous configurations with exceptional confinement and unprecedented compactness have been realized. Moreover, novel configurations like pseudosymmetry and piecewise omnigenity have been directly optimized for the first time. These advances enable efficient explorations of practical stellarator designs with enhanced confinement and engineering feasibility.

Subject: Plasma Physics

Publish: 2025-02-13 14:14:18 UTC


#17 The recurrence of groups inhibits the information spreading under higher-order interactions [PDF] [Copy] [Kimi] [REL]

Authors: Liang Yuan, Jiao Wu, Kesheng Xu, Muhua Zheng

Modeling social systems as networks based on pairwise interactions between individuals offers valuable insights into the mechanisms underlying their dynamics. However, the majority of social interactions occur within groups of individuals, characterized by higher-order structures. The mechanisms driving group formation and the impact of higher-order interactions, which arise from group dynamics, on information spreading in face-to-face interaction networks remain insufficiently understood. In this study, we examine some representative human face-to-face interaction data and find the recurrent patterns of groups. Moreover, we extend the force-directed motion (FDM) model with the forces derived from similarity distances within a hidden space to reproduce the recurrent group patterns and many key properties of face-to-face interaction networks. Furthermore, we demonstrate that the FDM model effectively predicts information-spreading behaviors under higher-order interactions. Finally, our results reveal that the recurrence of triangular groups inhibits the spread of information in face-to-face interaction networks, and the higher-order interactions will make this phenomenon more pronounced. These findings represent a significant advancement in the understanding of group formation and may open new avenues for research into the effects of group interactions on information propagation processes.

Subjects: Physics and Society , Chaotic Dynamics

Publish: 2025-02-13 14:11:41 UTC


#18 Disentangling the components of a multiconfigurational excited state in isolated chromophore [PDF] [Copy] [Kimi] [REL]

Authors: Rodrigo Cezar de Campos Ferreira, Amandeep Sagwal, Jiří Doležal, Tomáš Neuman, Martin Švec

Studying the excited states of doublets is challenging for their typically multiconfigurational character. We employ light-scanning-tunneling microscopy (light-STM) to investigate photon-induced currents on a single open-shell PTCDA anion molecule placed into a plasmonic nanocavity between a tip and a substrate, irradiated by laser. Submolecular mapping reveals a zero-bias bidirectional photocurrent strongly varying with the lateral position of the tip apex above the molecule. We elucidate the mechanism in terms of a theoretical model in which a multiconfigurational doublet state is excited and decays back to the anion ground state through sequential electron transfers with the tip and the substrate. The correspondence of the experimental and theoretical contrast proves the correlated character of the excited state which can be described as a superposition of two dominating electronic configurations. By applying bipolar voltage on the junction with the molecule, we switch the dominant recombination pathway from one of the configurations to the other, effectively disentangling the multiconfigurational state individual components through visualization of their Dyson orbitals, as corroborated by theoretical modelling.

Subject: Atomic and Molecular Clusters

Publish: 2025-02-13 14:11:37 UTC


#19 Leveraging the Bias-Variance Tradeoff in Quantum Chemistry for Accurate Negative Singlet-Triplet Gap Predictions: A Case for Double-Hybrid DFT [PDF] [Copy] [Kimi] [REL]

Authors: Atreyee Majumdar, Raghunathan Ramakrishnan

Molecules that violate Hund's rule -- having first excited singlet state (S$_1$) below the triplet state (T$_1$) -- are rare yet promising as efficient light emitters. Their high-throughput identification demands exceptionally accurate excited-state modeling to minimize false positives and negatives. Benchmarking twelve S$_1$-T$_1$ energy gaps, we find that local variants of ADC(2) and CC2 deliver excellent accuracy and speed for screening medium-sized molecules. Notably, while double-hybrid DFT approximations (e.g., B2GP-PLYP and PBE-QIDH) exhibit high mean errors ($>100$ meV) despite very low standard deviations ($\approx10$ meV), exploring their parameter space reveals that a configuration with 75% exchange and 55% correlation reduces the mean error to below $5$ meV -- albeit with increased variance. Using this low-bias parameterization as an internal reference, we correct the systematic error while maintaining low variance, effectively combining the strengths of both low-bias and low-variance DFT parameterizations to enhance overall accuracy. Our findings suggest that low-variance DFT methods -- often overlooked due to high bias -- can serve as reliable tools for predictive modeling in first-principles molecular design.

Subject: Chemical Physics

Publish: 2025-02-13 13:47:16 UTC


#20 NextGenPB: an analytically-enabled super resolution and local (de)refinement Poisson-Boltzmann Equation solver [PDF] [Copy] [Kimi] [REL]

Authors: Vincenzo Di Florio, Patrizio Ansalone, Sergii V. Siryk, Sergio Decherchi, Carlo de Falco, Walter Rocchia

The Poisson-Boltzmann equation (PBE) is a relevant partial differential equation commonly used in biophysical applications to estimate the electrostatic energy of biomolecular systems immersed in electrolytic solutions. A conventional mean to improve the accuracy of its solution, when grid-based numerical techniques are used, consists in increasing the resolution, locally or globally. This, however, usually entails higher complexity, memory demand and computational cost. Here, we introduce NextGenPB, a linear PBE, adaptive-grid, FEM solver that leverages analytical calculations to maximize the accuracy-to-computational-cost ratio. Indeed, in NextGenPB (aka NGPB), analytical corrections at the surface of the solute enhance the solution's accuracy without requiring grid adaptation. This leads to more precise estimates of the electrostatic potential, fields, and energy at no perceptible additional cost. Also, we apply computationally efficient yet accurate boundary conditions by taking advantage of local grid de-refinement. To assess the accuracy of our methods directly, we expand the traditionally available analytical case set to many non-overlapping dielectric spheres. Then, we use an existing benchmark set of real biomolecular systems to evaluate the energy convergence concerning grid resolution. Thanks to these advances, we have improved state-of-the-art results and shown that the approach is accurate and largely scalable for modern high-performance computing architectures. Lastly, we suggest that the presented core ideas could be instrumental in improving the solution of other partial differential equations with discontinuous coefficients.

Subjects: Computational Physics , Soft Condensed Matter , Biological Physics , Chemical Physics

Publish: 2025-02-13 13:37:00 UTC


#21 Experimental demonstration of a beam shaping non-imaging metasurface [PDF] [Copy] [Kimi] [REL]

Authors: Kirstine Engell Sandager Poulsen, Xavier Zambrana-Puyalto, Rafael de la Fuente Herrezuelo, Villads Egede Johansen, Søren Raza

Light emitting diodes have superior performance over most other light sources, but the need for secondary optics to shape their illumination for specific applications yield bulky lighting products. Here, we present an approach to shaping light from incoherent sources, such as light emitting diodes, using non-imaging metasurfaces. We present a theoretical framework and a numerical tool for designing the metasurface phase, and use it to construct a proof-of-principle beam shaping metasurface. We demonstrate the optical performance of the fabricated sample, and perform numerical experiments to investigate the quality of the metasurface design. Our approach bridges the fields of non-imaging optics and metaoptics, and may enable metasurface applications for shaping coherent and incoherent light.

Subject: Optics

Publish: 2025-02-13 13:34:50 UTC


#22 Glacier data assimilation on an Arctic glacier: Learning from large ensemble twin experiments [PDF] [Copy] [Kimi] [REL]

Authors: Wenxue Cao, Kristoffer Aalstad, Louise S. Schmidt, Sebastian Westermann, Thomas V. Schuler

Glacier modeling is crucial for quantifying the evolution of cryospheric processes. At the same time, uncertainties hamper process understanding and predictive accuracy. Here, we suggest improving glacier mass balance simulations for the Kongsvegen glacier in Svalbard through the application of Bayesian data assimilation techniques in a set of large ensemble twin experiments. Noisy synthetic observations of albedo and snow depth, generated using the multilayer CryoGrid community model with a full energy balance, are assimilated using two ensemble-based data assimilation schemes: the particle batch smoother and the ensemble smoother. A comprehensive evaluation exercise demonstrates that the joint assimilation of albedo and snow depth improves the simulation skill by up to 86% relative to the prior in specific glacier regions. The particle batch smoother excels in representing albedo dynamics, while the ensemble smoother is particularly effective for snow depth under low snowfall conditions. By combining the strengths of both observations, the joint assimilation achieves improved mass balance simulations across different glacier zones using either assimilation scheme. This work underscores the potential of ensemble-based data assimilation methods for refining glacier models by offering a robust framework to enhance predictive accuracy and reduce uncertainties in cryospheric simulations. Further advances in glacier data assimilation will be critical to better understanding the fate and role of Arctic glaciers in a changing climate.

Subjects: Geophysics , Data Analysis, Statistics and Probability

Publish: 2025-02-13 13:28:03 UTC


#23 Natural van der Waals canalization lens for non-destructive nanoelectronic circuit imaging and inspection [PDF] [Copy] [Kimi] [REL]

Authors: Qingdong Ou, Shuwen Xue, Weiliang Ma, Jiong Yang, Guangyuan Si, Lu Liu, Gang Zhong, Jingying Liu, Zongyuan Xie, Ying Xiao, Kourosh Kalantar-Zadeh, Xiang Qi, Peining Li, Zhigao Dai, Huanyang Chen, Qiaoliang Bao

Optical inspection has long served as a cornerstone non-destructive method in semiconductor wafer manufacturing, particularly for surface and defect analysis. However, conventional techniques such as bright-field and dark-field scattering optics face significant limitations, including insufficient resolution and the inability to penetrate and detect buried structures. Atomic force microscopy (AFM), while offering higher resolution and precise surface characterization, is constrained by slow speed, limited to surface-level imaging, and incapable of resolving subsurface features. Here, we propose an approach that integrates the strengths of dark-field scattering optics and AFM by leveraging a van der Waals (vdW) canalization lens based on natural biaxial {\alpha}-MoO3 crystals. This method enables ultrahigh-resolution subwavelength imaging with the ability to visualize both surface and buried structures, achieving a spatial resolution of 15 nm and grating pitch detection down to 100 nm. The underlying mechanism relies on the unique anisotropic properties of {\alpha}-MoO3, where its atomic-scale unit cells and biaxial symmetry facilitate the diffraction-free propagation of both evanescent and propagating waves via a flat-band canalization regime. Unlike metamaterial-based superlenses and hyperlenses, which suffer from high plasmonic losses, fabrication imperfections, and uniaxial constraints, {\alpha}-MoO3 provides robust and aberration-free imaging in multiple directions. We successfully applied this approach to high-resolution inspection of buried nanoscale electronic circuits, offering unprecedented capabilities essential for next-generation semiconductor manufacturing.

Subject: Optics

Publish: 2025-02-13 13:20:21 UTC


#24 Wetting simulation of the porous structure of a heat pipe using an eXtended Discontinuous Galerkin Method and a Parameterized Level-Set [PDF] [Copy] [Kimi] [REL]

Authors: Irina Shishkina, Matthias Rieckmann, Martin Oberlack, Florian Kummer

We perform high-order simulations of two-phase flows in capillaries, with and without evaporation. Since a sharp-interface model is used, singularities can arise at the three-phase contact line, where the fluid-fluid interface interacts with the capillary wall. These singularities are especially challenging when a highly accurate, high-order method with very little numerical diffusion is used for the flow solver. In this work, we employ the eXtended Discontinuous Galerkin (XDG) method, which has a very high accuracy but a severe limit regarding e.g., the time-step restriction. To address this challenge and enhance the stability of our numerical method we introduce a novel approach for representing a moving interface in the case of two-phase flows. We propose a global analytical representation of the interface-describing level-set field, defined by a small set of time-dependent parameters. Noteworthy for its simplicity and efficiency, this method effectively addresses the inherent complexity of two-phase flow problems. Furthermore, it significantly improves numerical stability and enables the use of larger time steps, ensuring both reliability and computational efficiency in our simulations. We compare different analytic expressions for level-set representation, including the elliptic function and the fourth-order polynomial, and validate the method against established literature data for capillary rise, both with and without evaporation. These results highlight the effectiveness of our approach in resolving complex interfacial dynamics.

Subject: Fluid Dynamics

Publish: 2025-02-13 13:14:54 UTC


#25 Internal Aerodynamics of Supersonic Crossflows with Transverse Liquid Injection [PDF] [Copy] [Kimi] [REL]

Authors: Srinivas M V V, Arun Kumar Rajagopal, Lebonah B, Jegesh David M

This study experimentally investigates the internal aerodynamics of transverse liquid injection in a supersonic crossflow (Mach (M) = 2.1) using two configurations: single and tandem (8 mm spacing) at three injection mass flow rates. Back-lit imaging revealed classical jet breakup phenomena, including surface wave instabilities with increasing amplitudes along the jet boundary, leading to protrusions, breakup into large liquid clumps, and their disintegration into finer droplets under aerodynamic forces. The single injection exhibited the formation of large liquid clumps further downstream compared to the tandem injection. Schlieren imaging showed that at a low momentum flux ratio (J = 0.94), both configurations produced regular reflection (RR) of the bow shock wave from the top wall. Increasing J to 1.90 resulted in RR for the single injection, while the tandem injection transitioned to Mach reflection (MR). At J = 2.67, both configurations exhibited MR. The earlier RR-to MR transition in tandem injection is attributed to its higher jet penetration and spanwise spread, which reduce the downstream crossflow passage area, acting as a supersonic diffuser and increasing downstream pressure which is favorable for MR transition. Separation zones were observed at the bottom wall due to bow shock wave-boundary layer interaction, and at the side walls due to the interaction of the Mach stem of the MR structure with the walls. These interactions create complex flow regions dominated by vortex structures, significantly influencing the overall flow dynamics.

Subject: Fluid Dynamics

Publish: 2025-02-13 13:09:46 UTC