Haokui Jiang, Jean-Lou Pfister, Daniel Zhengyu Huang, Shunxiang Cao
{"title":"Koopman reduced-order modeling and analysis of flag flapping in the wake of a cylinder.","authors":"Haokui Jiang, Jean-Lou Pfister, Daniel Zhengyu Huang, Shunxiang Cao","doi":"10.1103/PhysRevE.111.045101","DOIUrl":"https://doi.org/10.1103/PhysRevE.111.045101","url":null,"abstract":"<p><p>We develop a Koopman reduced-order model (ROM) to analyze the instability mechanism and predict the hydrodynamic behavior for the flag flapping in the wake of a cylinder. The Koopman ROM is constructed using a kernel dynamic mode decomposition method and enhanced through a residual dynamical mode decomposition algorithm, which improves accuracy by identifying and eliminating spurious modes. Our analysis reveals a flow transition from the \"2S\" mode in the periodic phase to the \"2P\" mode in the quasiperiodic phase, with the main Koopman mode M_{1} providing insights into the instability mechanism. In the case of chaotic flapping at a Reynolds number of Re=1200, the Koopman ROM demonstrates high accuracy in predicting the chaotic fluid-structure interaction flow when comparing the fractal dimension d_{c} and the maximum Lyapunov exponent λ_{max} of the true and reconstructed flow. Additionally, we observe similar flag flapping and vortex shedding characteristics in the near-structure region throughout the investigated Reynolds number range Re∈[500,1200], leading to similar vorticity patterns for M_{1}. The flag flapping has a local minimum displacement at position x_{0}≈3.0 in the flag's displacement envelope. Notably, this position closely matches the critical position x_{c} obtained from global linear instability analysis at low Reynolds numbers. This conclusion aligns with the performance of the Koopman ROM using sparse measurement probes attached to the flag. The position of the probes influences the accuracy of the ROM, where higher accuracy corresponds to a larger displacement of the eigenfunction.</p>","PeriodicalId":20085,"journal":{"name":"Physical review. E","volume":"111 4-2","pages":"045101"},"PeriodicalIF":2.4,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144136410","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Physical meaning of principal component analysis for classical lattice systems with translational invariance.","authors":"Su-Chan Park","doi":"10.1103/PhysRevE.111.045301","DOIUrl":"https://doi.org/10.1103/PhysRevE.111.045301","url":null,"abstract":"<p><p>We explore the physical implications of applying principal component analysis (PCA) to translationally invariant classical systems defined on a d-dimensional hypercubic lattice. Using Rayleigh-Schrödinger perturbation theory, we demonstrate that the principal components are related to the reciprocal lattice vectors of the hypercubic lattice, and the corresponding eigenvalues are connected to the discrete Fourier transform of the sampled configurations. From a different perspective, we show that the PCA in question can be viewed as a numerical method for computing the ensemble average of the squared moduli of the Fourier transform of physical quantities. Our results also provide a way to determine approximately the principal components of a classical system with translational invariance without the need for matrix diagonalization.</p>","PeriodicalId":20085,"journal":{"name":"Physical review. E","volume":"111 4-2","pages":"045301"},"PeriodicalIF":2.4,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144136470","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Quantum entanglement of linearly coupled quantum harmonic oscillators.","authors":"D N Makarov, K A Makarova","doi":"10.1103/PhysRevE.111.044140","DOIUrl":"https://doi.org/10.1103/PhysRevE.111.044140","url":null,"abstract":"<p><p>Quantum harmonic oscillators coupled through coordinates and momenta, represented by the Hamiltonian H[over ̂]=∑_{i=1}^{2}(p[over ̂]_{i}^{2}/2m_{i}+m_{i}ω_{i}^{2}/2x_{i}^{2})+H[over ̂]_{int}, where the interaction of two oscillators H[over ̂]_{int}=ik_{1}x_{1}p[over ̂]_{2}+ik_{2}x_{2}p[over ̂]_{1}+k_{3}x_{1}x_{2}-k_{4}p[over ̂]_{1}p[over ̂]_{2}, are found in many applications of quantum optics, nonlinear physics, molecular chemistry, and biophysics. Despite this, there is currently no general solution to the Schrödinger equation for such a system. This is especially relevant for quantum entanglement of such a system in quantum optics applications. Here this problem is solved and it is shown that quantum entanglement depends on only one coefficient, R∈(0,1), which includes all the parameters of the system under consideration. It has been shown that quantum entanglement can be very large at certain values of this coefficient. The results obtained have a fairly simple analytical form, which facilitates analysis.</p>","PeriodicalId":20085,"journal":{"name":"Physical review. E","volume":"111 4-1","pages":"044140"},"PeriodicalIF":2.4,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144136503","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
K S Krishnamurthy, S Krishna Prasad, D S Shankar Rao, R J Mandle, C J Gibb, J Hobbs, N V Madhusudana
{"title":"Static and hydrodynamic periodic structures induced by ac electric fields in the antiferroelectric SmZ_{A} phase.","authors":"K S Krishnamurthy, S Krishna Prasad, D S Shankar Rao, R J Mandle, C J Gibb, J Hobbs, N V Madhusudana","doi":"10.1103/PhysRevE.111.045420","DOIUrl":"https://doi.org/10.1103/PhysRevE.111.045420","url":null,"abstract":"<p><p>We report the effect of AC electric fields in the range of 0.1-300 kHz on antiferroelectric SmZ_{A} layers of DIO in the bookshelf geometry. Significant results are (a) primary bifurcation into a quasistationary periodic instability with its voltage threshold U_{c} and wave vector q_{c} along the initial director being, respectively, quadratic and linear functions of f over 10-150 kHz, and with an azimuthal distortion of the director n which changes sign between adjacent stripes, (b) transition from the modulated state to a homogeneous state at higher voltages, and (c) third bifurcation into traveling wave periodic state on further rise in U in the region 10-40 kHz. We interpret these findings as follows. The low voltage instability is very similar to that seen in the higher temperature apolar nematic phase and is the electrohydrodynamic (EHD) instability possibly belonging to the region of dielectric inversion frequency. The azimuthal distortions of n result from an undulatory distortion of the SmZ_{A} layers in the book-shelf geometry. The intermediate homogeneous state of SmZ_{A} in which the periodic structure is absent results from a linear coupling between the layer polarization P and applied field E, giving rise to a scissoring type mutual P reorientation in adjacent layers. Finally, at even higher voltages, the medium goes over to a field-induced transition to the ferroelectric nematic, with the polarization following the AC field, and the periodic EHD instability being similar to that of the dielectric regime. The polar vector symmetry of the medium leads in general to traveling waves.</p>","PeriodicalId":20085,"journal":{"name":"Physical review. E","volume":"111 4-2","pages":"045420"},"PeriodicalIF":2.4,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144136527","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Simulating squirmers with smoothed particle dynamics.","authors":"Xinwei Cai, Kuiliang Wang, Gaojin Li, Xin Bian","doi":"10.1103/PhysRevE.111.045401","DOIUrl":"https://doi.org/10.1103/PhysRevE.111.045401","url":null,"abstract":"<p><p>Microswimmers play an important role in shaping the world around us. The squirmer is a simple model of a microswimmer whose cilia oscillations on its spherical surface induce an effective slip velocity to propel itself. The rapid development of computational fluid dynamics methods has markedly enhanced our capacity to study the behavior of squirmers in aqueous environments. Nevertheless, a unified methodology that can fully address the complexity of fluid-solid coupling at multiple scales and interface tracking for multiphase flows remains elusive, posing an outstanding challenge to the field. To this end, we investigate the potential of the smoothed particle dynamics (SPD) method as an alternative approach to simulating squirmers. The Lagrangian nature of the method allows it to effectively address the aforementioned difficulty. By introducing a novel treatment of the boundary condition and assigning appropriate slip velocities to the boundary particles, the SPD squirmer model is able to accurately represent a range of microswimmer types, including pushers, neutral swimmers, and pullers. We systematically validate the steady-state velocity of the squirmer, the resulting flow field, its hydrodynamic interactions with the surrounding environment, and the mutual collision of two squirmers. In the presence of Brownian motion, the model is also able to correctly calculate the velocity and angular velocity autocorrelation functions at the mesoscale. Finally, we simulate a squirmer within a multiphase flow by considering a droplet that encloses a squirmer and imposing a surface tension between the two flow phases. We find that the squirmer within the droplet exhibits different motion types. Since the proposed method is applicable to a wide range of complex scenarios, it has implications for a number of areas, including the design and application of micro and/or nano artificial swimmers, flow manipulation in microfluidic chips, and drug delivery in the biomedical field.</p>","PeriodicalId":20085,"journal":{"name":"Physical review. E","volume":"111 4-2","pages":"045401"},"PeriodicalIF":2.4,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144136558","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Spatial organization of multiple species of active particles interacting with an interface.","authors":"Love Grover, Rajeev Kapri, Abhishek Chaudhuri","doi":"10.1103/PhysRevE.111.045412","DOIUrl":"https://doi.org/10.1103/PhysRevE.111.045412","url":null,"abstract":"<p><p>We investigate the steady-state organization of active particles residing on an interface. Particle activity induces interface deformations, while the local shape of the interface guides particle movement. We consider multiple species of particles which can locally pull on the interface or push it. This coupled system exhibits a wide variety of behaviors, including clustering, anticlustering, diffusion, mixing, demixing, and localization. Our findings suggest that one can control surface properties by strategically adding or removing specific particle types. Furthermore, by adjusting particle activity levels, we can selectively disperse particle types, enabling precise manipulation of surface movement and geometry.</p>","PeriodicalId":20085,"journal":{"name":"Physical review. E","volume":"111 4-2","pages":"045412"},"PeriodicalIF":2.4,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144136562","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Revisiting classical nucleation theory: Insights into heterogeneous ice nucleation on nanoscale substrates.","authors":"Yufeng Liu, Jincheng Zeng, Yu Zhang, Jianyang Wu, Zhisen Zhang","doi":"10.1103/PhysRevE.111.044107","DOIUrl":"https://doi.org/10.1103/PhysRevE.111.044107","url":null,"abstract":"<p><p>Heterogeneous nucleation plays a pivotal role in the ice nucleation process. Within the classical nucleation theory (CNT) framework, the heterogeneous nucleation rate is proportional to the substrate surface area, typically assuming infinite substrate surfaces. However, when the substrate size approaches the nanoscale, the nucleation rate deviates significantly from CNT predictions. This study presents a novel theoretical model that distinguishes the nanoscale substrate into central and edge regions, attributing different contributions to ice nucleation. We hypothesize that the edge width equals the critical size of the nucleus (r_{c}) and validate this hypothesis using molecular dynamics (MD) simulations with the coarse-grained water model (mW model) on circular and rectangular substrates of varying sizes. Our results demonstrate that the edge region impedes heterogeneous ice nucleation, with the MD calculated nucleation rates aligning well with our model. Furthermore, the statistical edge width matches the critical nucleus size r_{c}. By incorporating this refined model, our findings reconcile the nucleation rates with CNT predictions, offering new insights into heterogeneous ice nucleation at the nanoscale.</p>","PeriodicalId":20085,"journal":{"name":"Physical review. E","volume":"111 4-1","pages":"044107"},"PeriodicalIF":2.4,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144136572","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
I Ioannou Sougleridis, O Richoux, V Achilleos, G Theocharis, D J Frantzeskakis
{"title":"Ring-shaped linear waves and solitons in a square lattice of acoustic waveguides.","authors":"I Ioannou Sougleridis, O Richoux, V Achilleos, G Theocharis, D J Frantzeskakis","doi":"10.1103/PhysRevE.111.044201","DOIUrl":"https://doi.org/10.1103/PhysRevE.111.044201","url":null,"abstract":"<p><p>We study the propagation of both low- and high-amplitude ring-shaped sound waves in a two-dimensional square lattice of acoustic waveguides with Helmholtz resonators. We show that the inclusion of the Helmholtz resonators suppresses the inherent anisotropy of the system in the low-frequency regime allowing for radially symmetric solutions. By employing the electroacoustic analogue approach and asymptotic methods we derive an effective cylindrical Korteweg-de Vries equation. Low-amplitude waveforms are self-similar structures of the Airy function profile, while high-amplitude ones are of the form of cylindrical solitons. Our analytical predictions are corroborated by results of direct numerical simulations, with a very good agreement between the two.</p>","PeriodicalId":20085,"journal":{"name":"Physical review. E","volume":"111 4-1","pages":"044201"},"PeriodicalIF":2.4,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144136576","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Liming Zhang, Lan Zhang, Shun Gao, Changwei Huang, Qionglin Dai
{"title":"Structural balance and evolution of cooperation in a population with hybrid interactions.","authors":"Liming Zhang, Lan Zhang, Shun Gao, Changwei Huang, Qionglin Dai","doi":"10.1103/PhysRevE.111.044309","DOIUrl":"https://doi.org/10.1103/PhysRevE.111.044309","url":null,"abstract":"<p><p>This study explores the evolution of cooperation in populations with mixed pairwise and three-body interactions, investigating the impact of higher-order interaction density ρ and individual interaction preference α. Our results reveal that sparse higher-order interactions markedly boost cooperation, exhibiting two critical phase transitions as ρ changes. These transitions underscore the delicate equilibrium needed for optimal cooperation, as excessive higher-order interactions can diminish returns. The preference parameter α significantly influences cooperation sustainability, with intermediate values maximizing cooperative outcomes, particularly when the temptation to defect r is not strong. Crucially, our findings demonstrate that hybrid social dilemmas structurally encode emergent cooperation pathways that are unattainable within homogeneous interaction frameworks, emphasizing the importance of modeling mixed interactions to capture real-world complexity. These insights offer valuable guidance for designing systems aimed at promoting cooperative behavior across social, ecological, and artificial domains.</p>","PeriodicalId":20085,"journal":{"name":"Physical review. E","volume":"111 4-1","pages":"044309"},"PeriodicalIF":2.4,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144136605","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Yan Xia, Zhaosheng Yu, Minkang Zhang, Zhaowu Lin, Zhenyu Ouyang
{"title":"Swimming dynamics of a spheroidal microswimmer near a wall.","authors":"Yan Xia, Zhaosheng Yu, Minkang Zhang, Zhaowu Lin, Zhenyu Ouyang","doi":"10.1103/PhysRevE.111.045106","DOIUrl":"https://doi.org/10.1103/PhysRevE.111.045106","url":null,"abstract":"<p><p>In this work, we investigate the swimming dynamics of a spheroidal squirmer near a flat wall for various aspect ratios using the direct-forcing fictitious domain method. Our results show that the swimming mode of a strong pusher undergoes the transition from either oscillating or escaping to crawling as the aspect ratio increases. A strong puller exhibits an opposite transition: from crawling to escaping and then to oscillating as the aspect ratio increases. The mechanism for the near-wall swimming behavior of a strong puller and pusher is explored by analyzing the hydrodynamic force and torque on a swimmer with its height and orientation fixed. The results indicate that both collision and hydrodynamic toques are important to the near-wall swimming state of the squirmer. Additionally, we found that the initial orientation angle and the release distance do not influence the swimming mode when the squirmer initially swims toward the wall at an angle smaller than -π/8.</p>","PeriodicalId":20085,"journal":{"name":"Physical review. E","volume":"111 4-2","pages":"045106"},"PeriodicalIF":2.4,"publicationDate":"2025-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144136582","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}