Modelling of the Electrical Membrane Potential for Concentration Polarization Conditions

Based on Kedem–Katchalsky formalism, the model equation of the membrane potential (Δψs) generated in a membrane system was derived for the conditions of concentration polarization. In this system, a horizontally oriented electro-neutral biomembrane separates solutions of the same electrolytes at different concentrations. The consequence of concentration polarization is the creation, on both sides of the membrane, of concentration boundary layers. The basic equation of this model includes the unknown ratio of solution concentrations (Ci/Ce) at the membrane/concentration boundary layers. We present the calculation procedure (Ci/Ce) based on novel equations derived in the paper containing the transport parameters of the membrane (Lp, σ, and ω), solutions (ρ, ν), concentration boundary layer thicknesses (δl, δh), concentration Raileigh number (RC), concentration polarization factor (ζs), volume flux (Jv), mechanical pressure difference (ΔP), and ratio of known solution concentrations (Ch/Cl). From the resulting equation, Δψs was calculated for various combinations of the solution concentration ratio (Ch/Cl), the Rayleigh concentration number (RC), the concentration polarization coefficient (ζs), and the hydrostatic pressure difference (ΔP). Calculations were performed for a case where an aqueous NaCl solution with a fixed concentration of 1 mol m⁻³ (Cl) was on one side of the membrane and on the other side an aqueous NaCl solution with a concentration between 1 and 15 mol m⁻³ (Ch). It is shown that (Δψs) depends on the value of one of the factors (i.e., ΔP, Ch/Cl, RC and ζs) at a fixed value of the other three.     

The Double-Sided Information Bottleneck Function

A double-sided variant of the information bottleneck method is considered. Let (X,Y) be a bivariate source characterized by a joint pmf PXY. The problem is to find two independent channels PU|X and PV|Y (setting the Markovian structure U→X→Y→V), that maximize I(U;V) subject to constraints on the relevant mutual information expressions: I(U;X) and I(V;Y). For jointly Gaussian X and Y, we show that Gaussian channels are optimal in the low-SNR regime but not for general SNR. Similarly, it is shown that for a doubly symmetric binary source, binary symmetric channels are optimal when the correlation is low and are suboptimal for high correlations. We conjecture that Z and S channels are optimal when the correlation is 1 (i.e., X=Y) and provide supporting numerical evidence. Furthermore, we present a Blahut–Arimoto type alternating maximization algorithm and demonstrate its performance for a representative setting. This problem is closely related to the domain of biclustering.     

Hypercontractive Inequalities for the Second Norm of Highly Concentrated Functions,and Mrs. Gerber’s-Type Inequalities for the Second Rényi Entropy

Let Tϵ, 0≤ϵ≤1/2, be the noise operator acting on functions on the boolean cube {0,1}n. Let f be a distribution on {0,1}n and let q>1. We prove tight Mrs. Gerber-type results for the second Rényi entropy of Tϵf which take into account the value of the qth Rényi entropy of f. For a general function f on {0,1}n we prove tight hypercontractive inequalities for the ℓ2 norm of Tϵf which take into account the ratio between ℓq and ℓ1 norms of f.     

Estimating Drift Parameters in a Sub-Fractional Vasicek-Type Process

This study deals with drift parameters estimation problems in the sub-fractional Vasicek process given by dxt=θ(μ−xt)dt+dStH, with θ>0, μ∈R being unknown and t≥0; here, SH represents a sub-fractional Brownian motion (sfBm). We introduce new estimators θ^ for θ and μ^ for μ based on discrete time observations and use techniques from Nordin–Peccati analysis. For the proposed estimators θ^ and μ^, strong consistency and the asymptotic normality were established by employing the properties of SH. Moreover, we provide numerical simulations for sfBm and related Vasicek-type process with different values of the Hurst index H.     

Four-Objective Optimizations of a Single Resonance Energy Selective Electron Refrigerator

According to the established model of a single resonance energy selective electron refrigerator with heat leakage in the previous literature, this paper performs multi-objective optimization with finite-time thermodynamic theory and NSGA-II algorithm. Cooling load (R¯), coefficient of performance (ε), ecological function (ECO¯), and figure of merit (χ¯) of the ESER are taken as objective functions. Energy boundary (E′/kB) and resonance width (ΔE/kB) are regarded as optimization variables and their optimal intervals are obtained. The optimal solutions of quadru-, tri-, bi-, and single-objective optimizations are obtained by selecting the minimum deviation indices with three approaches of TOPSIS, LINMAP, and Shannon Entropy; the smaller the value of deviation index, the better the result. The results show that values of E′/kB and ΔE/kB are closely related to the values of the four optimization objectives; selecting the appropriate values of the system can design the system for optimal performance. The deviation indices are 0.0812 with LINMAP and TOPSIS approaches for four-objective optimization (ECO¯−R¯−ε−χ¯), while the deviation indices are 0.1085, 0.8455, 0.1865, and 0.1780 for four single-objective optimizations of maximum ECO¯, R¯, ε, and χ¯, respectively. Compared with single-objective optimization, four-objective optimization can better take different optimization objectives into account by choosing appropriate decision-making approaches. The optimal values of E′/kB and ΔE/kB range mainly from 12 to 13, and 1.5 to 2.5, respectively, for the four-objective optimization.     

Contrarian Voter Model under the Influence of an Oscillating Propaganda: Consensus,Bimodal Behavior and Stochastic Resonance

We study the contrarian voter model for opinion formation in a society under the influence of an external oscillating propaganda and stochastic noise. Each agent of the population can hold one of two possible opinions on a given issue—against or in favor—and interacts with its neighbors following either an imitation dynamics (voter behavior) or an anti-alignment dynamics (contrarian behavior): each agent adopts the opinion of a random neighbor with a time-dependent probability p(t), or takes the opposite opinion with probability 1−p(t). The imitation probability p(t) is controlled by the social temperature T, and varies in time according to a periodic field that mimics the influence of an external propaganda, so that a voter is more prone to adopt an opinion aligned with the field. We simulate the model in complete graph and in lattices, and find that the system exhibits a rich variety of behaviors as T is varied: opinion consensus for T=0, a bimodal behavior for T<Tc, an oscillatory behavior where the mean opinion oscillates in time with the field for T>Tc, and full disorder for T≫1. The transition temperature Tc vanishes with the population size N as Tc≃2/lnN in complete graph. In addition, the distribution of residence times tr in the bimodal phase decays approximately as tr−3/2. Within the oscillatory regime, we find a stochastic resonance-like phenomenon at a given temperature T*. Furthermore, mean-field analytical results show that the opinion oscillations reach a maximum amplitude at an intermediate temperature, and that exhibit a lag with respect to the field that decreases with T.     

Four-Objective Optimization of an Irreversible Stirling Heat Engine with Linear Phenomenological Heat-Transfer Law

This paper combines the mechanical efficiency theory and finite time thermodynamic theory to perform optimization on an irreversible Stirling heat-engine cycle, in which heat transfer between working fluid and heat reservoir obeys linear phenomenological heat-transfer law. There are mechanical losses, as well as heat leakage, thermal resistance, and regeneration loss. We treated temperature ratio x of working fluid and volume compression ratio λ as optimization variables, and used the NSGA-II algorithm to carry out multi-objective optimization on four optimization objectives, namely, dimensionless shaft power output P¯s, braking thermal efficiency ηs, dimensionless efficient power E¯p and dimensionless power density P¯d. The optimal solutions of four-, three-, two-, and single-objective optimizations are reached by selecting the minimum deviation indexes D with the three decision-making strategies, namely, TOPSIS, LINMAP, and Shannon Entropy. The optimization results show that the D reached by TOPSIS and LINMAP strategies are both 0.1683 and better than the Shannon Entropy strategy for four-objective optimization, while the Ds reached for single-objective optimizations at maximum P¯s, ηs, E¯p, and P¯d conditions are 0.1978, 0.8624, 0.3319, and 0.3032, which are all bigger than 0.1683. This indicates that multi-objective optimization results are better when choosing appropriate decision-making strategies.     

On the Validity of Detrended Fluctuation Analysis at Short Scales

Detrended Fluctuation Analysis (DFA) has become a standard method to quantify the correlations and scaling properties of real-world complex time series. For a given scale ℓ of observation, DFA provides the function F(ℓ), which quantifies the fluctuations of the time series around the local trend, which is substracted (detrended). If the time series exhibits scaling properties, then F(ℓ)∼ℓα asymptotically, and the scaling exponent α is typically estimated as the slope of a linear fitting in the logF(ℓ) vs. log(ℓ) plot. In this way, α measures the strength of the correlations and characterizes the underlying dynamical system. However, in many cases, and especially in a physiological time series, the scaling behavior is different at short and long scales, resulting in logF(ℓ) vs. log(ℓ) plots with two different slopes, α1 at short scales and α2 at large scales of observation. These two exponents are usually associated with the existence of different mechanisms that work at distinct time scales acting on the underlying dynamical system. Here, however, and since the power-law behavior of F(ℓ) is asymptotic, we question the use of α1 to characterize the correlations at short scales. To this end, we show first that, even for artificial time series with perfect scaling, i.e., with a single exponent α valid for all scales, DFA provides an α1 value that systematically overestimates the true exponent α. In addition, second, when artificial time series with two different scaling exponents at short and large scales are considered, the α1 value provided by DFA not only can severely underestimate or overestimate the true short-scale exponent, but also depends on the value of the large scale exponent. This behavior should prevent the use of α1 to describe the scaling properties at short scales: if DFA is used in two time series with the same scaling behavior at short scales but very different scaling properties at large scales, very different values of α1 will be obtained, although the short scale properties are identical. These artifacts may lead to wrong interpretations when analyzing real-world time series: on the one hand, for time series with truly perfect scaling, the spurious value of α1 could lead to wrongly thinking that there exists some specific mechanism acting only at short time scales in the dynamical system. On the other hand, for time series with true different scaling at short and large scales, the incorrect α1 value would not characterize properly the short scale behavior of the dynamical system.     

Electron Kinetic Entropy across Quasi-Perpendicular Shocks

We use Magnetospheric Multiscale (MMS) data to study electron kinetic entropy per particle Se across Earth’s quasi-perpendicular bow shock. We have selected 22 shock crossings covering a wide range of shock conditions. Measured distribution functions are calibrated and corrected for spacecraft potential, secondary electron contamination, lack of measurements at the lowest energies and electron density measurements based on plasma frequency measurements. All crossings display an increase in electron kinetic entropy across the shock ΔSe being positive or zero within their error margin. There is a strong dependence of ΔSe on the change in electron temperature, ΔTe, and the upstream electron plasma beta, βe. Shocks with large ΔTe have large ΔSe. Shocks with smaller βe are associated with larger ΔSe. We use the values of ΔSe, ΔTe and density change Δne to determine the effective adiabatic index of electrons for each shock crossing. The average effective adiabatic index is 〈γe〉=1.64±0.07.     

Analyses of pp,Cu–Cu,Au–Au and Pb–Pb Collisions by Tsallis-Pareto Type Function at RHIC and LHC Energies

The parameters revealing the collective behavior of hadronic matter extracted from the transverse momentum spectra of π+, π−, K+, K−, p, p¯, Ks0, Λ, Λ¯, Ξ or Ξ−, Ξ¯+ and Ω or Ω¯+ or Ω+Ω¯ produced in the most central and most peripheral gold–gold (Au–Au), copper–copper (Cu–Cu) and lead–lead (Pb–Pb) collisions at 62.4 GeV, 200 GeV and 2760 GeV, respectively, are reported. In addition to studying the nucleus–nucleus (AA) collisions, we analyzed the particles mentioned above produced in pp collisions at the same center of mass energies (62.4 GeV, 200 GeV and 2760 GeV) to compare with the most peripheral AA collisions. We used the Tsallis–Pareto type function to extract the effective temperature from the transverse momentum spectra of the particles. The effective temperature is slightly larger in a central collision than in a peripheral collision and is mass-dependent. The mean transverse momentum and the multiplicity parameter (N0) are extracted and have the same result as the effective temperature. All three extracted parameters in pp collisions are closer to the peripheral AA collisions at the same center of mass energy, revealing that the extracted parameters have the same thermodynamic nature. Furthermore, we report that the mean transverse momentum in the Pb–Pb collision is larger than that of the Au–Au and Cu–Cu collisions. At the same time, the latter two are nearly equal, which shows their comparatively strong dependence on energy and weak dependence on the size of the system. The multiplicity parameter, N0 in central AA, depends on the interacting system’s size and is larger for the bigger system.     

Quantum Linear System Algorithm for General Matrices in System Identification

Solving linear systems of equations is one of the most common and basic problems in classical identification systems. Given a coefficient matrix A and a vector b, the ultimate task is to find the solution x such that Ax=b. Based on the technique of the singular value estimation, the paper proposes a modified quantum scheme to obtain the quantum state |x〉 corresponding to the solution of the linear system of equations in O(κ2rpolylog(mn)/ϵ) time for a general m×n dimensional A, which is superior to existing quantum algorithms, where κ is the condition number, r is the rank of matrix A and ϵ is the precision parameter. Meanwhile, we also design a quantum circuit for the homogeneous linear equations and achieve an exponential improvement. The coefficient matrix A in our scheme is a sparsity-independent and non-square matrix, which can be applied in more general situations. Our research provides a universal quantum linear system solver and can enrich the research scope of quantum computation.     

New Bounds and a Generalization for Share Conversion for 3-Server PIR

Private Information Retrieval (PIR) protocols, which allow the client to obtain data from servers without revealing its request, have many applications such as anonymous communication, media streaming, blockchain security, advertisement, etc. Multi-server PIR protocols, where the database is replicated among the non-colluding servers, provide high efficiency in the information-theoretic setting. Beimel et al. in CCC 12’ (further referred to as BIKO) put forward a paradigm for constructing multi-server PIR, capturing several previous constructions for k≥3 servers, as well as improving the best-known share complexity for 3-server PIR. A key component there is a share conversion scheme from corresponding linear three-party secret sharing schemes with respect to a certain type of “modified universal” relation. In a useful particular instantiation of the paradigm, they used a share conversion from (2,3)-CNF over Zm to three-additive sharing over Zpβ for primes p1,p2,p where p1≠p2 and m=p1·p2, and the relation is modified universal relation CSm. They reduced the question of the existence of the share conversion for a triple (p1,p2,p) to the (in)solvability of a certain linear system over Zp, and provided an efficient (in m,logp) construction of such a sharing scheme. Unfortunately, the size of the system is Θ(m2) which entails the infeasibility of a direct solution for big m’s in practice. Paskin-Cherniavsky and Schmerler in 2019 proved the existence of the conversion for the case of odd p1, p2 when p=p1, obtaining in this way infinitely many parameters for which the conversion exists, but also for infinitely many of them it remained open. In this work, using some algebraic techniques from the work of Paskin-Cherniavsky and Schmerler, we prove the existence of the conversion for even m’s in case p=2 (we computed β in this case) and the absence of the conversion for even m’s in case p>2. This does not improve the concrete efficiency of 3-server PIR; however, our result is promising in a broader context of constructing PIR through composition techniques with k≥3 servers, using the relation CSm where m has more than two prime divisors. Another our suggestion about 3-server PIR is that it’s possible to achieve a shorter server’s response using the relation CSm′ for extended Sm′⊃Sm. By computer search, in BIKO framework we found several such sets for small m’s which result in share conversion from (2,3)-CNF over Zm to 3-additive secret sharing over Zpβ′, where β′>0 is several times less than β, which implies several times shorter server’s response. We also suggest that such extended sets Sm′ can result in better PIR due to the potential existence of matching vector families with the higher Vapnik-Chervonenkis dimension.     

Search for Sub-Solar Mass Binaries with Einstein Telescope and Cosmic Explorer

A possible detection of sub-solar mass ultra-compact objects would lead to new perspectives on the existence of black holes that are not of astrophysical origin and/or pertain to formation scenarios of exotic ultra-compact objects. Both possibilities open new perspectives for better understanding of our universe. In this work, we investigate the significance of detection of sub-solar mass binaries with components mass in the range: 10−2M⊙ up to 1M⊙, within the expected sensitivity of the ground-based gravitational waves detectors of third generation, viz., the Einstein Telescope (ET) and the Cosmic Explorer (CE). Assuming a minimum of amplitude signal-to-noise ratio for detection, viz., ρ=8, we find that the maximum horizon distances for an ultra-compact binary system with components mass 10−2M⊙ and 1M⊙ are 40 Mpc and 1.89 Gpc, respectively, for ET, and 125 Mpc and 5.8 Gpc, respectively, for CE. Other cases are also presented in the text. We derive the merger rate and discuss consequences on the abundances of primordial black hole (PBH), fPBH. Considering the entire mass range [10−2–1]M⊙, we find fPBH<0.70 (<0.06) for ET (CE), respectively.     

A Compressed Sensing Measurement Matrix Construction Method Based on TDMA for Wireless Sensor Networks

Compressed sensing theory has been widely used for data aggregation in WSNs due to its capability of containing much information but with light load of transmission. However, there still exist some issues yet to be solved. For instance, the measurement matrix is complex to construct, and it is difficult to implement in hardware and not suitable for WSNs with limited node energy. To solve this problem, a random measurement matrix construction method based on Time Division Multiple Access (TDMA) is proposed based on the sparse random measurement matrix combined with the data transmission method of the TDMA of nodes in the cluster. The reconstruction performance of the number of non-zero elements per column in this matrix construction method for different signals was compared and analyzed through extensive experiments. It is demonstrated that the proposed matrix can not only accurately reconstruct the original signal, but also reduce the construction complexity from O(MN) to O(d2N) (d≪M), on the premise of achieving the same reconstruction effect as that of the sparse random measurement matrix. Moreover, the matrix construction method is further optimized by utilizing the correlation theory of nested matrices. A TDMA-based semi-random and semi-deterministic measurement matrix construction method is also proposed, which significantly reduces the construction complexity of the measurement matrix from O(d2N) to O(dN), and improves the construction efficiency of the measurement matrix. The findings in this work allow more flexible and efficient compressed sensing for data aggregation in WSNs.     

Constraints on Non-Flat Starobinsky f(R) Dark Energy Model

We study the viable Starobinsky f(R) dark energy model in spatially non-flat FLRW backgrounds, where f(R)=R−λRch[1−(1+R2/Rch2)−1] with Rch and λ representing the characteristic curvature scale and model parameter, respectively. We modify CAMB and CosmoMC packages with the recent observational data to constrain Starobinsky f(R) gravity and the density parameter of curvature ΩK. In particular, we find the model and density parameters to be λ−1<0.283 at 68% C.L. and ΩK=−0.00099−0.0042+0.0044 at 95% C.L., respectively. The best χ2 fitting result shows that χf(R)2≲χΛCDM2, indicating that the viable f(R) gravity model is consistent with ΛCDM when ΩK is set as a free parameter. We also evaluate the values of AIC, BIC and DIC for the best fitting results of f(R) and ΛCDM models in the non-flat universe.     

Penalty Virtual Element Method for the 3D Incompressible Flow on Polyhedron Mesh

In this paper, a penalty virtual element method (VEM) on polyhedral mesh for solving the 3D incompressible flow is proposed and analyzed. The remarkable feature of VEM is that it does not require an explicit computation of the trial and test space, thereby bypassing the obstacle of standard finite element discretizations on arbitrary mesh. The velocity and pressure are approximated by the practical significative lowest equal-order virtual element space pair (Xh,Qh) which does not satisfy the discrete inf-sup condition. Combined with the penalty method, the error estimation is proved rigorously. Numerical results on the 3D polygonal mesh illustrate the theoretical results and effectiveness of the proposed method.     

λ-Deformation: A Canonical Framework for Statistical Manifolds of Constant Curvature

This paper systematically presents the λ-deformation as the canonical framework of deformation to the dually flat (Hessian) geometry, which has been well established in information geometry. We show that, based on deforming the Legendre duality, all objects in the Hessian case have their correspondence in the λ-deformed case: λ-convexity, λ-conjugation, λ-biorthogonality, λ-logarithmic divergence, λ-exponential and λ-mixture families, etc. In particular, λ-deformation unifies Tsallis and Rényi deformations by relating them to two manifestations of an identical λ-exponential family, under subtractive or divisive probability normalization, respectively. Unlike the different Hessian geometries of the exponential and mixture families, the λ-exponential family, in turn, coincides with the λ-mixture family after a change of random variables. The resulting statistical manifolds, while still carrying a dualistic structure, replace the Hessian metric and a pair of dually flat conjugate affine connections with a conformal Hessian metric and a pair of projectively flat connections carrying constant (nonzero) curvature. Thus, λ-deformation is a canonical framework in generalizing the well-known dually flat Hessian structure of information geometry.     

The Longitudinal Plasma Modes of κ-Deformed Kaniadakis Distributed Plasmas Carrying Orbital Angular Momentum

Based on plasma kinetic theory, the dispersion and Landau damping of Langmuir and ion-acoustic waves carrying finite orbital angular momentum (OAM) were investigated in the κ-deformed Kaniadakis distributed plasma system. The results showed that the peculiarities of the investigated subjects relied on the deformation parameter κ and OAM parameter η. For both Langmuir and ion-acoustic waves, dispersion was enhanced with increased κ, while the Landau damping was suppressed. Conversely, both the dispersion and Landau damping were depressed by OAM. Moreover, the results coincided with the straight propagating plane waves in a Maxwellian plasma system when κ=0 and η→∞. It was expected that the present results would give more insight into the trapping and transportation of plasma particles and energy.