In the context of an ultracold Fermi gas, we derive conservation laws for mass, energy and momentum based on a generalized nonlocal Boltzmann equation with gradient corrections in the collision integral. The corrections are expressed in terms of effective collision duration, particle displacement and changes of total momentum and energy. Their origin is in the in-medium $T$ matrix. Using variations of the optical theorem, we show that in the collision integral the particle-hole symmetry can be recast into a form of collision symmetry amenable to semiclassical simulation. Pauli-blocked collisions are distinguished from Bose-stimulated nondissipative ones; the latter are not present in the absence of gradient corrections. Consolidating with the microscopic theory, we extract local conservation laws for a general time-dependent interaction potential, and demonstrate how both types of collisions affect densities and flows of conserving quantities. Comparison is made with the approach of Nozières and Schmitt-Rink in the limit of thermal equilibrium. Under approximations used for normal-state ultracold Fermi gases interacting via Feshbach resonances we demonstrate the effect of the collision delay on the shear viscosity.

High-density percolation is the formation of a system spanning cluster of vertices with at least $m$ occupied neighbors. We discuss high-density percolation on the modified Bethe lattice in terms of the theory of large random graphs with arbitrary degree distributions. Using the formalism of generating functions, we derive expressions for the cluster size distribution, the percolation threshold, the percolation probability, and the mean size of finite clusters. We show that the critical exponents $\beta =\gamma =1$. Additionally, numerical solutions and simulation results for the percolation probability and mean size of finite clusters are compared for illustration.

We provide a concrete expression for the phase-space distribution function at nonequilibrium steady state under a constant thermal gradient, which is a typical system of the nonequilibrium molecular dynamics simulation of heat conduction. First, the phase-space distribution function of all particles in a local volume is formulated. Our formulation explicitly takes into account the entropy production due to the change in equilibrium thermodynamic variables in addition to the traditional entropy production described by the spatial gradients and fluxes of equilibrium thermodynamic variables. This treatment is necessary to explain the nonequilibrium response of a quantity that has no equilibrium correlation with mass and heat fluxes and is essential to correctly deduce one-particle distribution functions from the all-particle one. From the all-particle distribution function, we derive the Green–Kubo relations that express the one-particle distribution functions of density and velocity in terms of equilibrium correlation functions and verify these expressions using the molecular dynamics simulation of a Lennard-Jones liquid. These nonequilibrium one-particle distribution functions are sufficiently tractable for practical use, such as for the analytical evaluation of the nonequilibrium average of physical quantities.

Topological heterogeneities of social networks have a strong impact on the individuals embedded in those networks. One of the interesting phenomena driven by such heterogeneities is the friendship paradox (FP), stating that the mean degree of one's neighbors is larger than the degree of oneself. Alternatively, one can use the median degree of neighbors as well as the fraction of neighbors having a higher degree than oneself. Each of these reflects on how people perceive their neighborhoods, i.e., their perception models, hence how they feel peer pressure. In our paper, we study the impact of perception models on the FP by comparing three versions of the perception model in networks generated with a given degree distribution and a tunable degree-degree correlation or assortativity. The increasing assortativity is expected to decrease network-level peer pressure, while we find a nontrivial behavior only for the mean-based perception model. By simulating opinion formation, in which the opinion adoption probability of an individual is given as a function of individual peer pressure, we find that it takes the longest time to reach consensus when individuals adopt the median-based perception model compared to other versions. Our findings suggest that one needs to consider the proper perception model for better modeling human behaviors and social dynamics.

This study proposes an innovative approach to estimate relative biological effectiveness (RBE) of fast neutrons using the Geant4 toolkit. The Geant4-DNA version cannot track heavy ions below 0.5 MeV/nucleon. In order to explore the impact of this issue, secondary particles are simulated instead of the primary low-energy neutrons. The Evaluated Nuclear Data File library is used to determine the cross sections for the elastic and inelastic interactions of neutrons with water and to find the contribution of each secondary particle spectrum. Two strategies are investigated in order to find the best possible approach and results. The first one takes into account only light particles, protons produced from elastic scattering, and $\alpha$ particles from inelastic scattering. Geantino particles are shot instead of heavy ions; hence all heavy ions are considered in the simulations, though their physical effects on DNA not. The second strategy takes into account all the heavy and light ions, although heavy ions cannot be tracked down to very low energies ($E<0.5$ MeV/nucleon). Our model is based on the combination of an atomic resolution DNA geometrical model and a Monte Carlo simulation toolkit for tracking particles. The atomic coordinates of the DNA double helix are extracted from the Protein Data Bank. Since secondary particle spectra are used instead of simulating the interaction of neutrons explicitly, this method reduces the computation times dramatically. Double-strand break induction is used as the end point for the estimation of the RBE of fast neutrons. ${}^{60}\mathrm{Co} \gamma$ rays are used as the reference radiation quality. Both strategies succeed in reproducing the behavior of the ${\mathrm{RBE}}_{\text{max}}$ as a function of the incident neutron energy ranging from 0.1 to 14 MeV, including the position of its peak. A comparison of the behavior of the two strategies shows that for neutrons with energies less than 0.7 MeV, the effect of heavy ions would not be very significant, but above 0.7 MeV, heavy ions have an important role in neutron RBE.

The mechanical behavior of polymers such as their probability density function (PDF), strain energy, and entropic force, has long been described by the non-Gaussian statistical model. Non-Gaussian models are often approximated by the Kuhn-Grün (KG) distribution function, which is derived from the first-order approximation of the complex Rayleigh's exact Fourier integral distribution. The KG function is widely accepted in polymer physics, where the non-Gaussian theory is often used to describe energy of the chains with various flexibility ratios. However, the KG function is shown to be relevant only for long chains and becomes extremely inaccurate for chains with fewer than 40 segments. In comparison to the KG model, other approximations of the non-Gaussian statistical model are often less accurate, and those with higher accuracy are usually too complex to be implemented in large-scale simulations. Here a new accurate approximation of the non-Gaussian PDF, entropic force, and strain energy of a single chain subsequently is developed to describe the mechanics of a polymer chain. With a similar level of complexity, the presented approximations of non-Gaussian PDF, strain energy, and entropic force are at least 10 times more accurate than KG approximations and thus are an excellent alternative option to be used in micromechanical constitutive models.

Braids composed of two interwoven polymer chains exhibit a “buckling” transition whose origin has been explained through the onset of plectonemic structures. Here we study, by a combination of simulation and analytics, the dynamics of plectoneme formation and their statistics in steady state. The introduction of an order parameter—the plectonemic fraction—allows us to map out the phase boundary between the straight-braid phase and the plectonemic one. We then monitor the formation and the growth of plectonemes, observing events typical of phase separation kinetics for liquid-gas systems (fusion, fission, and one-dimensional Ostwald ripening) but also of DNA supercoiling dynamics (plectonemic hopping). Finally, we propose a stochastic field theory for the coupled dynamics of twist and local writhe which explains the phenomenology found with Brownian dynamics simulations as well as the power laws underlying the coarsening of plectonemes.

The dynamical structure factor (DSF) of strongly coupled ions in dense plasmas with partially and strongly degenerate electrons is investigated. The main focus is on the impact of electronic correlations (nonideality) on the ionic DSF. The latter is computed by carrying out molecular dynamics (MD) simulations with a screened ion-ion interaction potential. The electronic screening is taken into account by invoking the Singwi-Tosi-Land-Sjölander approximation, and it is compared to the MD simulation data obtained considering the electronic screening in the random phase approximation and using the Yukawa potential. We find that electronic correlations lead to lower values of the ion-acoustic mode frequencies and to an extension of the applicability limit with respect to the wave-number of a hydrodynamic description. Moreover, we show that even in the limit of weak electronic coupling, electronic correlations have a nonnegligible impact on the ionic longitudinal sound speed. Additionally, the applicability of the Yukawa potential with an adjustable screening parameter is discussed, which will be of interest, e.g., for the interpretation of experimental results for the ionic DSF of dense plasmas.

Amorphous solids have excess soft modes in addition to the phonon modes described by the Debye theory. Recent numerical results show that if the phonon modes are carefully removed, the density of state of the excess soft modes exhibit universal quartic scaling, independent of the interaction potential, preparation protocol, and spatial dimensions. We hereby provide a theoretical framework to describe this universal scaling behavior. For this purpose, we extend the mean-field theory to include the effects of finite-dimensional fluctuation. Based on a semiphenomenological argument, we show that mean-field quadratic scaling is replaced by the quartic scaling in finite dimensions. Furthermore, we apply our formalism to explain the pressure and protocol dependence of the excess soft modes.

We introduce a model based on the one-dimensional nonlinear Schrödinger equation with critical (quintic) or supercritical self-focusing nonlinearity. We demonstrate that a family of solitons, which are unstable in this setting against the critical or supercritical collapse, is stabilized by pinning to an attractive defect, that may also include a parity-time ($\mathcal{PT}$)-symmetric gain-loss component. The model can be realized as a planar waveguide in nonlinear optics, and in a super-Tonks-Girardeau bosonic gas. For the attractive defect with the delta-functional profile, a full family of the pinned solitons is found in an exact analytical form. In the absence of the gain-loss term, the solitons' stability is investigated in an analytical form too, by means of the Vakhitov-Kolokolov criterion; in the presence of the $\mathcal{PT}$-balanced gain and loss, the stability is explored by means of numerical methods. In particular, the entire family of pinned solitons is stable in the quintic (critical) medium if the gain-loss term is absent. A stability region for the pinned solitons persists in the model with an arbitrarily high power of the self-focusing nonlinearity. A weak gain-loss component gives rise to intricate alternations of stability and instability in the system's parameter plane. Those solitons which are unstable under the action of the supercritical self-attraction are destroyed by the collapse. On the other hand, if the self-attraction-driven instability is weak and the gain-loss term is present, unstable solitons spontaneously transform into localized breathers, while the collapse does not occur. The same outcome may be caused by a combination of the critical nonlinearity with the gain and loss. Instability of the solitons is also possible when the $\mathcal{PT}$-symmetric gain-loss term is added to the subcritical nonlinearity. The system with self-repulsive nonlinearity is briefly considered too, producing completely stable families of pinned localized states.

Cortical circuits operate in an inhibition-dominated regime of spiking activity. Recently, it was found that spiking circuit models in this regime can, despite disordered connectivity and asynchronous irregular activity, exhibit a locally stable dynamics that may be used for neural computation. The lack of existing mathematical tools has precluded analytical insight into this phase. Here we present analytical methods tailored to the granularity of spike-based interactions for analyzing attractor geometry in high-dimensional spiking dynamics. We apply them to reveal the properties of the complex geometry of trajectories of population spiking activity in a canonical model of locally stable spiking dynamics. We find that attractor basin boundaries are the preimages of spike-time collision events involving connected neurons. These spike-based instabilities control the divergence rate of neighboring basins and have no equivalent in rate-based models. They are located according to the disordered connectivity at a random subset of edges in a hypercube representation of the phase space. Iterating backward these edges using the stable dynamics induces a partition refinement on this space that converges to the attractor basins. We formulate a statistical theory of the locations of such events relative to attracting trajectories via a tractable representation of local trajectory ensembles. Averaging over the disorder, we derive the basin diameter distribution, whose characteristic scale emerges from the relative strengths of the stabilizing inhibitory coupling and destabilizing spike interactions. Our study provides an approach to analytically dissect how connectivity, coupling strength, and single-neuron dynamics shape the phase space geometry in the locally stable regime of spiking neural circuit dynamics.

A clinical study of tinnitus patients found promising results using a noninvasive therapy. We introduce a dynamical model to explore both the onset of tinnitus and the effects of coordinated reset therapy. Our model extends an existing theory of individual outer hair cell dynamics to include their mutual interaction, and considers how sustained activity can inhibit the natural recovery exhibited by normal (healthy) individuals. The model is investigated through numerical simulations and shows behavior broadly similar to that reported in the clinical study.

The turbulent boundary layer flow measured by two-dimensional time-resolved particle image velocimetry is analyzed using the discrete orthogonal wavelet method. The Reynolds number of the turbulent boundary layer based on the friction velocity is ${\mathrm{Re}}_{\tau }=235$. The flow field is decomposed into a number of wavelet levels which have different characteristic scales. The velocity statistics and coherent structures at different wavelet levels are investigated. It is found that the fluctuation intensities and their peak locations differ for varying scales. The proper orthogonal decomposition (POD) of different wavelet components reveals a cascade of scales of coherent structures, especially the small-scale ones that are usually difficult to be identified in POD modes of the undecomposed flow field. The interactions among the scales are investigated in terms of large-scale amplitude modulations of the small-scale structures. In previous studies the velocity fluctuations are separated into two parts, the large scale and the small scale, divided usually by the boundary layer thickness. In the present study, however, the scales smaller than the boundary layer thickness are further separated. Therefore, the modulation analysis is a refined investigation that differentiates the modulation effects on separated small scales. The results reveal that the modulation effects vary among the small scales.

A theory is developed to model the nonlinear dynamics of two coupled bubbles inside a spherical liquid-filled cavity surrounded by an elastic medium. The aim is to study how the conditions of full confinement affect the coupled oscillations of the bubbles. To make the problem amenable to analytical consideration, the bubbles are assumed to be located on a diameter of the cavity, which makes the problem axisymmetric. Equations for the pulsation and translation motion of the bubbles are derived by the Lagrangian formalism. The derived equations are used in numerical simulations. The behavior of two bubbles in a cavity is compared with the behavior of the same bubbles in an unbounded liquid. It is found that both forced and free oscillations of two bubbles in a cavity occur differently than those in an unbounded liquid. In particular, it is shown that the eigenfrequencies of a two-bubble system in a cavity are different from those in an unbounded liquid.

We present a local approach to the study of scale-to-scale energy transfers in magnetohydrodynamic (MHD) turbulence. This approach is based on performing local averages of the physical fields, which amounts to filtering scales smaller than some parameter $\ell$. A key step in this work is the derivation of a local Kármán-Howarth-Monin relation which provides a local form of Politano and Pouquet's 4/3 law, without any assumption of homogeneity or isotropy. Our approach is exact and nonrandom, and we show its connection to the usual statistical results of turbulence. Its implementation on data obtained via a three-dimensional direct numerical simulation of the forced incompressible MHD equations from the John Hopkins turbulence database constitutes the main part of our study. First, we show that the local Kármán-Howarth-Monin relation holds well. The space statistics of local cross-scale transfers are studied next, their means and standard deviations being maximum at inertial scales and their probability density functions (PDFs) displaying very wide tails. Events constituting the tails of the PDFs are shown to form structures of strong transfers, either positive or negative, which can be observed over the whole available range of scales. As $\ell$ is decreased, these structures become more and more localized in space while contributing to an increasing fraction of the mean energy cascade rate. Finally, we highlight their quasi-one-dimensional (filamentlike) or quasi-two-dimensional (sheetlike or ribbonlike) nature and show that they appear in areas of strong vorticity or electric current density.

We consider a harmonic charging field as an energy charger for the quantum battery, which consists of an ensemble of two-level atoms. The charging of noninteracting atoms is completely fulfilled, which exhibits a substantial improvement over previous static charging fields. Involving the repulsive interactions of atoms, the fully charging is achieved with shorter charged period over the noninteracting case, yielding an advantage for the charging. Excluding the charging field, a quantum phase transition is induced by the attractive atom-atom interactions, and the interacting atoms become degenerate in the ground state. We find that the degenerate states play a negative role in the charging due to the gapless energies. The atoms with strong attractive interactions can not be charged completely, which is accompanied by a drop of the maximum stored energy.

We consider a model of a particle trapped in a harmonic optical trap but with the addition of a nonconservative radiation induced force. This model is known to correctly describe experimentally observed trapped particle statistics for a wide range of physical parameters, such as temperature and pressure. We theoretically analyze the effect of nonconservative force on the underlying steady state distribution as well as the power spectrum for the particle position. We compute perturbatively the probability distribution of the resulting nonequilibrium steady states for all dynamical regimes underdamped through to overdamped and give expressions for the associated currents in phase space (position and velocity). We also give the spectral density of the trapped particle's position in all dynamical regimes and for any value of the nonconservative force. Signatures of the presence of nonconservative forces are shown to be particularly strong for the underdamped regime at low frequencies.

Transition from a split to a forward kinetic energy cascade system is explored in the context of rotating turbulence using direct numerical simulations with a three-dimensional isotropic random force uncorrelated with the velocity field. Our parametric study covers confinement effects in high-aspect-ratio domains and a broad range of rotation rates. The data presented here add substantially to previous works, which, in contrast, focused on smaller and shallower domains. Results indicate that for fixed geometrical dimensions the Rossby number acts as a control parameter, whereas for a fixed Rossby number the product of the domain size along the rotation axis and the forcing wave number governs the amount of energy that cascades inversely. The regime transition criterion hence depends on both control parameters.

A rotating wheel experimental facility was developed to investigate incompressible Rayleigh-Taylor instability in elastic-plastic materials. A soft solid (mayonnaise) was chosen as the elastic-plastic material for experiments; material properties that include shear modulus and yield strength were fully characterized using a vane spindle type rheometer. Initial perturbations of varying amplitudes and wavelengths were generated on the interface of the soft solid using sinusoidal cutting guides. A backlit imaging technique was used in conjunction with a high-speed camera to track the motion of the interface at various phases of the instability. Results for both two- and three-dimensional perturbations were compared to study the acceleration required for instability and the growth after the interface yielded. Exponential growth rates were observed after instability was reached with trends of increasing growth rates for lower initial amplitudes. It was found that the acceleration required for instability increased when initial amplitude and wavelength decreased. Three-dimensional interfaces were found to be more stable. For both cases, a decrease in initial amplitude produced a more stable interface that increased the threshold acceleration required for the instability. Critical amplitude conditions for instability were calculated and compared with various analytical models and other experimental results.

Early afterdepolarizations (EADs) are pathological voltage oscillations in cardiomyocytes that have been observed in response to a number of pharmacological agents and disease conditions. Phase-2 EADs consist of small voltage fluctuations during the plateau of an action potential, typically under conditions in which the action potential is elongated. Although a single-cell behavior, EADs can lead to tissue-level arrhythmias. Much is currently known about the biophysical mechanisms (i.e., the roles of ion channels and intracellular ${\mathrm{Ca}}^{2+}$ stores) for the various forms of EADs, due partially to the development and analysis of mathematical models. This includes the application of slow-fast analysis, which takes advantage of timescale separation inherent in the system to simplify its analysis. We take this further, using a minimal three-dimensional model to demonstrate that phase-2 EADs are canards formed in the neighborhood of a folded node singularity. This allows us to predict the number of EADs that can be produced for a given parameter set, and provides guidance on parameter changes that facilitate or inhibit EAD production. With this approach, we demonstrate why periodic stimulation, as occurs in intact heart, preferentially facilitates EAD production when applied at low frequencies. We also explain the origin of complex alternan dynamics that can occur with intermediate-frequency stimulation, in which varying numbers of EADs are produced with each pulse. These revelations fall out naturally from an understanding of folded node singularities, but are difficult to glean from knowledge of the biophysical mechanism for EADs alone. Therefore, understanding the canard mechanism is a useful complement to understanding of the biophysical mechanism that has been developed over years of experimental and computational investigations.