Exergy is a very important thermodynamic quantity in several fields such as economy, engineering, and ecology and yet it has attracted little attention in pure physics. One of the main problems of the currently used definition of exergy is its dependence on an arbitrarily chosen reference state, which is the thermodynamic state of a reservoir the system is supposedly in contact with. In this paper, starting from a very general definition of exergy, a formula is derived for the exergy balance of a general open continuous medium without any reference to an external environment. A formula is also derived for the most suitable thermodynamic parameters of the Earth atmosphere when seen as an external environment in the usual exergy applications.

We analyzed the translational and rotational Brownian motion of aggregates of micrometer-sized silica spheres under microgravity conditions and in rarefied gas. The experimental data was collected in the form of high-speed recordings using a long-distance microscope as part of the ICAPS (Interactions in Cosmic and Atmospheric Particle Systems) experiment on board of the sounding rocket flight Texus-56. Our data analysis shows that the translational Brownian motion can be used to determine the mass and translational response time of each individual dust aggregate. The rotational Brownian motion additionally provides the moment of inertia and the rotational response time. A shallow positive correlation between mass and response time was found as predicted for aggregate structures with low fractal dimensions. Translational and rotational response times are roughly in agreement. Using the mass and the moment of inertia of each aggregate, we determined the fractal dimension of the aggregate ensemble. Slight deviations from the pure Gaussian one-dimensional displacement statistics were found in the ballistic limit for both the translational and rotational Brownian motion.

We consider the instability and stability of periodic stationary solutions to the classical ${\varphi }^4$ equation numerically. In the superluminal regime, the model possesses dnoidal and cnoidal waves. The former are modulationally unstable and the spectrum forms a figure eight intersecting at the origin of the spectral plane. The latter can be modulationally stable, and the spectrum near the origin in that case is represented by vertical bands along the purely imaginary axis. The instability of the cnoidal states in that case stems from elliptical bands of complex eigenvalues far from the spectral plane origin. In the subluminal regime, there exist only snoidal waves which are modulationally unstable. Considering the subharmonic perturbations, we show that the snoidal waves in the subluminal regime are spectrally unstable with respect to all subharmonic perturbations, while for the dnoidal and cnoidal waves in the superluminal regime, the transition between the spectrally stable state and the spectrally unstable state occurs through a Hamiltonian Hopf bifurcation. The dynamical evolution of the unstable states is also considered, leading to some interesting localization events on the spatio-temporal backgrounds.

The propagation of acoustic and elastic waves in time-varying, spatially homogeneous media can exhibit different phenomena when compared to traditional spatially varying, temporally homogeneous media. In the present work, the response of a one-dimensional phononic lattice with time-periodic elastic properties is studied with experimental, numerical and theoretical approaches in both linear and nonlinear regimes. The system consists of repelling magnetic masses with grounding stiffness controlled by electrical coils driven with electrical signals that vary periodically in time. For small-amplitude excitation, in agreement with linear theoretical predictions, wave-number band gaps emerge. The underlying instabilities associated to the wave-number band gaps are investigated with Floquet theory and the resulting parametric amplification is observed in both theory and experiments. In contrast to genuinely linear systems, large-amplitude responses are stabilized via the nonlinear nature of the magnetic interactions of the system, and results in a family of nonlinear time-periodic states. The bifurcation structure of the periodic states is studied. It is found the linear theory correctly predicts parameter values from which the time-periodic states bifurcate from the zero state. In the presence of an external drive, the parametric amplification induced by the wave-number band gap can lead to bounded and stable responses that are temporally quasiperiodic. Controlling the propagation of acoustic and elastic waves by balancing nonlinearity and external modulation offers a new dimension in the realization of advanced signal processing and telecommunication devices. For example, it could enable time-varying, cross-frequency operation, mode- and frequency-conversion, and signal-to-noise ratio enhancements.

We consider growing directed network models that aim at minimizing the weighted connection expenses while at the same time favoring other important network properties such as weighted local node degrees. We employed statistical mechanics methods to study the growth of directed networks under the principle of optimizing some objective function. By mapping the system to an Ising spin model, analytic results are derived for two such models, exhibiting diverse and interesting phase transition behaviors for general edge weight, inward and outward node weight distributions. In addition, the unexplored cases of negative node weights are also investigated. Analytic results for the phase diagrams are derived showing even richer phase transition behavior, such as first-order transition due to symmetry, second-order transitions with possible reentrance, and hybrid phase transitions. We further extend previously developed zero-temperature simulation algorithm for undirected networks to the present directed case and for negative node weights, and we can obtain the minimal cost connection configuration efficiently. All the theoretical results are explicitly verified by simulations. Possible applications and implications are also discussed.

A diffusive trajectory drawn by the generalized Langevin equation (GLE) for a colloidal particle evokes a random fractal of a static polymer configuration. This article proposes a static GLE-like description that enables the generation of a single configuration of a polymer chain with the noise formulated to satisfy the static fluctuation-response relation (FRR) along a one-dimensional chain structure but not along a temporal coordinate. A remarkable point is qualitative differences and similarities in the FRR formulation between the static and the dynamic GLEs. Guided by the static FRR, we further make analogous arguments in light of stochastic energetics and the steady-state fluctuation theorem.

Electrophoretic (EP) mobility reversal is commonly observed for strongly charged macromolecules in multivalent salt solutions. This curious effect takes place, e.g., when a charged polymer, such as DNA, adsorbs excess counterions so that the counterion-dressed surface charge reverses its sign, leading to the inversion of the polymer drift driven by an external electric field. In order to characterize this seemingly counterintuitive phenomenon that cannot be captured by electrostatic mean-field theories, we adapt here a previously developed strong-coupling-dressed Poisson-Boltzmann approach to the cylindrical geometry of the polyelectrolyte-salt system. Within the framework of this formalism, we derive an analytical polymer mobility formula dressed by charge correlations. In qualitative agreement with polymer transport experiments, this mobility formula predicts that the increment of the monovalent salt, the decrease of the multivalent counterion valency, and the increase of the dielectric permittivity of the background solvent suppress charge correlations and increase the multivalent bulk counterion concentration required for EP mobility reversal. These results are corroborated by coarse-grained molecular dynamics simulations showing how multivalent counterions induce mobility inversion at dilute concentrations and suppress the inversion effect at large concentrations. This re-entrant behavior, previously observed in the aggregation of like-charged polymer solutions, calls for verification by polymer transport experiments.

The dynamics of a harmonically trapped three-dimensional Yukawa ball of charged dust particles immersed in plasma is investigated as function of external magnetic field and Coulomb coupling parameter using molecular dynamics simulation. It is shown that the harmonically trapped dust particles organize themselves into nested spherical shells. The particles start rotating in a coherent order as the magnetic field reaches a critical value corresponding to the coupling parameter of the system of dust particles. The magnetically controlled charged dust cluster of finite size undergoes a first-order phase transition from disordered to ordered phase. At sufficiently high coupling and strong magnetic field, the vibrational mode of this finite-sized charged dust cluster freezes, and the system retains only rotational motion.

The Beltrami state in a single-species (electron or ion) ideal plasma in the context of massive electromagnetism has been explored. The inclusion of photon mass, essentially treating the massive photon field as a mobile fluid in the ideal plasma vortical dynamics, has resulted in a triple curl Beltrami state of the magnetic vector potential $\vec{A}$. A variational principle is constructed which shows this state can also be obtained by constrained minimization of the system's energy with appropriate helicity invariants. Such a state is endowed with three different length scales; one of which is the system length, and the other two are species' skin depth and photon Compton wavelength, respectively. An analytical solution of this state in the cylindrical geometry is presented, which is the linear combination of three single Beltrami states. Possible observational signatures of this state in astrophysical and laboratory settings are also discussed.

The emergence of a power-law distribution for the energy released during an earthquake is investigated in several models. Generic features are identified which are based on the self-affine behavior of the stress field prior to an event. This field behaves at large scale as a random trajectory in one dimension of space and a random surface in two dimensions. Using concepts of statistical mechanics and results on the properties of these random objects, several predictions are obtained and verified, in particular the value of the power-law exponent of the earthquake energy distribution (the Gutenberg-Richter law) as well as a mechanism for the existence of aftershocks after a large earthquake (the Omori law).

Localization behaviors of Laplacian eigenvectors of complex networks furnish an explanation to various dynamical phenomena of the corresponding complex systems. We numerically examine roles of higher-order and pairwise links in driving eigenvector localization of hypergraphs Laplacians. We find that pairwise interactions can engender localization of eigenvectors corresponding to small eigenvalues for some cases, whereas higher-order interactions, even being much much less than the pairwise links, keep steering localization of the eigenvectors corresponding to larger eigenvalues for all the cases considered here. These results will be advantageous to comprehend dynamical phenomena, such as diffusion, and random walks on a range of real-world complex systems having higher-order interactions in better manner.

When a ferrofluid is magnetized in a strong magnetic field, and then the field is switched off, the magnetization decays from its saturation value to zero. The dynamics of this process are controlled by the rotations of the constituent magnetic nanoparticles, and for the Brownian mechanism, the respective rotation times are strongly influenced by the particle size and the magnetic dipole-dipole interactions between the particles. In this work, the effects of polydispersity and interactions on the magnetic relaxation are studied using a combination of analytical theory and Brownian dynamics simulations. The theory is based on the Fokker-Planck-Brown equation for Brownian rotation and includes a self-consistent, mean-field treatment of the dipole-dipole interactions. The most interesting predictions from the theory are that, at short times, the relaxation of each particle type is equal to its intrinsic Brownian rotation time, while at long times, each particle type has the same effective relaxation time, which is longer than any of the individual Brownian rotation times. Noninteracting particles, though, always relax at a rate controlled only by the Brownian rotation times. This illustrates the importance of including the effects of polydispersity and interactions when analyzing the results from magnetic relaxometry experiments on real ferrofluids, which are rarely monodisperse.

Fractal properties on self-affine surfaces of films growing under nonequilibrium conditions are important in understanding the corresponding universality class. However, measurement of the surface fractal dimension has been intensively investigated and is still very problematic. In this work, we report the behavior of the effective fractal dimension in the context of film growth involving lattice models believed to belong to the Kardar-Parisi-Zhang (KPZ) universality class. Our results, which are presented for growth in a $d$-dimensional substrate ($d=1,2$) and use the three-point sinuosity (TPS) method, show universal scaling of the measure $M$, which is defined in terms of discretization of the Laplacian operator applied to the height of the film surface, $M=t^{\delta }g\left[\mathrm{\Theta }\right]$, where $t$ is the time, $g\left[\mathrm{\Theta }\right]$ is a scale function, $\delta =2\beta , \mathrm{\Theta }\equiv \tau t^{-1/z}, \beta$, and $z$ are the KPZ growth and dynamical exponents, respectively, and $\tau$ is a spatial scale length used to compute $M$. Importantly, we show that the effective fractal dimensions are consistent with the expected KPZ dimensions for $d=1,2$, if $\mathrm{\Theta }\lesssim 0.3$, which include a thin film regime for the extraction of the fractal dimension. This establishes the scale limits in which the TPS method can be used to accurately extract effective fractal dimensions that are consistent with those expected for the corresponding universality class. As a consequence, for the steady state, which is inaccessible to experimentalists studying film growth, the TPS method provided effective fractal dimension consistent with the KPZ ones for almost all possible $\tau$, i.e., $1\lesssim \tau <L/2$, where $L$ is the lateral size of the substrate on which the deposit is grown. In the growth of thin films, the true fractal dimension can be observed in a narrow range of $\tau$, the upper limit of which is of the same order of magnitude as the correlation length of the surface, indicating the limits of self-affinity of a surface in an experimentally accessible regime. This upper limit was comparatively lower for the Higuchi method or the height-difference correlation function. Scaling corrections for the measure $M$ and the height-difference correlation function are studied analytically and compared for the Edwards-Wilkinson class at $d=1$, yielding similar accuracy for both methods. Importantly, we extend our discussion to a model representing diffusion-dominated growth of films and find that the TPS method achieves the corresponding fractal dimension only at steady state and in a narrow range of the scale length, compared to that found for the KPZ class.

Despite recent advances in pore-scale modeling of two-phase flow through porous media, the relative strengths and limitations of various modeling approaches have been largely unexplored. In this work, two-phase flow simulations from the generalized network model (GNM) [Phys. Rev. E 96, 013312 (2017); Phys. Rev. E 97, 023308 (2018)] are compared with a recently developed lattice-Boltzmann model (LBM) [Adv. Water Resour. 116, 56 (2018); J. Colloid Interface Sci. 576, 486 (2020)] for drainage and waterflooding in two samples—a synthetic beadpack and a micro-CT imaged Bentheimer sandstone—under water-wet, mixed-wet, and oil-wet conditions. Macroscopic capillary pressure analysis reveals good agreement between the two models, and with experiments, at intermediate saturations but shows large discrepancy at the end-points. At a resolution of 10 grid blocks per average throat, the LBM is unable to capture the effect of layer flow which manifests as abnormally large initial water and residual oil saturations. Critically, pore-by-pore analysis shows that the absence of layer flow limits displacement to invasion-percolation in mixed-wet systems. The GNM is able to capture the effect of layers, and exhibits predictions closer to experimental observations in water and mixed-wet Bentheimer sandstones. Overall, a workflow for the comparison of pore-network models with direct numerical simulation of multiphase flow is presented. The GNM is shown to be an attractive option for cost and time-effective predictions of two-phase flow, and the importance of small-scale flow features in the accurate representation of pore-scale physics is highlighted.

Almost every quantum circuit is built with two-qubit gates in the current stage, which are crucial to the quantum computing in any platform. The entangling gates based on Mølmer-Sørensen schemes are widely exploited in the trapped-ion system, with the utilization of the collective motional modes of ions and two laser-controlled internal states, which are served as qubits. The key to realize high-fidelity and robust gates is the minimization of the entanglement between the qubits and the motional modes under various sources of errors after the gate operation. In this work, we propose an efficient numerical method to search high-quality solutions for phase-modulated pulses. Instead of directly optimizing a cost function, which contains the fidelity and the robustness of the gates, we convert the problem to the combination of linear algebra and the solution to quadratic equations. Once a solution with the gate fidelity of 1 is found, the laser power can be further reduced while searching on the manifold where the fidelity remains 1. Our method largely overcomes the problem of the convergence and is shown to be effective up to 60 ions, which suffices the need of the gate design in current trapped-ion experiments.

This paper shows that, in one dimension, due to its topology, a closed-loop product of short-time propagators is always positive, despite the fact that each antisymmetric free fermion propagator can be of either sign.

Network correlation dimension governs the distribution of network distance in terms of a power-law model and profoundly impacts both structural properties and dynamical processes. We develop new maximum likelihood methods which allow us robustly and objectively to identify network correlation dimension and a bounded interval of distances over which the model faithfully represents structure. We also compare the traditional practice of estimating correlation dimension by modeling as a power law the fraction of nodes within a distance to a proposed alternative of modeling as a power law the fraction of nodes $\mathit{at}$ a distance. In addition, we illustrate a likelihood ratio technique for comparing the correlation dimension and small-world descriptions of network structure. Improvements from our innovations are demonstrated on a diverse selection of synthetic and empirical networks. We show that the network correlation dimension model accurately captures empirical network structure over neighborhoods of substantial size and span and outperforms the alternative small-world network scaling model. Our improved methods tend to lead to higher estimates of network correlation dimension, implying that prior studies could have produced or utilized systematic underestimates of dimension.

The stress-strain curve of two-dimensional frictional dispersed grains interacting with a harmonic potential without considering the dynamical slip under a finite strain is determined by using eigenvalue analysis of the Hessian matrix. After the configuration of grains is obtained, the stress-strain curve based on the eigenvalue analysis is in almost perfect agreement with that obtained by the simulation, even if there are plastic deformations caused by stress avalanches. Unlike the naive expectation, the eigenvalues in our model do not indicate any precursors to the stress-drop events.

Droplets suspended in fluids flowing through microchannels are often encountered in different contexts and scales, from oil extraction down to microfluidics. They are usually flexible and deform as a product of the interplay between flexibility, hydrodynamics, and interaction with confining walls. Deformability adds distinct characteristics to the nature of the flow of these droplets. We simulate deformable droplets suspended in a fluid at a high volume fraction flowing through a cylindrical wetting channel. We find a discontinuous shear thinning transition, which depends on the droplet deformability. The capillary number is the main dimensionless parameter that controls the transition. Previous results have focused on two-dimensional configurations. Here we show that, in three dimensions, even the velocity profile is different. To perform this study, we improve and extend to three dimensions a multicomponent lattice Boltzmann method which prevents the coalescence between the droplets.

Strong field physics close to or above the Schwinger limit are typically studied with vacuum as initial condition or by considering test particle dynamics. However, with a plasma present initially, quantum relativistic mechanisms such as Schwinger pair creation are complemented by classical plasma nonlinearities. In this work we use the Dirac-Heisenberg-Wigner formalism to study the interplay between classical and quantum mechanical mechanisms in the regime of ultrastrong electric fields. In particular, the effects of initial density and temperature on the plasma oscillation dynamics are determined. Finally, comparisons with competing mechanisms such as radiation reaction and Breit-Wheeler pair production are made.