We explore the Casimir-Lifshitz free-energy theory for surface freezing of methane gas hydrates near the freezing point of water. The theory enables us to explore different pathways, resulting in anomalous (stabilizing) ice layers on methane hydrate surfaces via energy minimization. Notably, we will contrast the gas hydrate material properties, under which thin ice films can form in water vapor, with those previously predicted to be required in the presence of liquid water. It is predicted that methane hydrates in water vapor near the freezing point of water nucleate ice films, instead of water films.

The ion velocity distribution functions of thermonuclear plasmas generated by spherical laser direct drive implosions are studied using deuterium-tritium (DT) and deuterium-deuterium (DD) fusion neutron energy spectrum measurements. A hydrodynamic Maxwellian plasma model accurately describes measurements made from lower temperature $(<10 \mathrm{keV})$, hydrodynamiclike plasmas, but is insufficient to describe measurements made from higher temperature more kineticlike plasmas. The high temperature measurements are more consistent with Vlasov-Fokker-Planck (VFP) simulation results which predict the presence of a bimodal plasma ion velocity distribution near peak neutron production. These measurements provide direct experimental evidence of non-Maxwellian ion velocity distributions in spherical shock driven implosions and provide useful data for benchmarking kinetic VFP simulations.

Strong shocks are essential components in many high-energy-density environments such as inertial confinement fusion implosions. However, the experimental measurements of the spatial structures of such shocks are sparse. In this paper, the soft x-ray emission of a shock front in a helium gas mixture ($90\%$ helium, $10\%$ neon) and a pure neon gas was spatially resolved using an imaging spectrometer. We observe that the shock width in the helium mixture gas is about twice as large as in the pure neon gas. Moreover, they exhibit different precursor layers, where electron temperature greatly exceeds ion temperature, extending for more than $\sim 350\mathrm{\mu }\mathrm{m}$ with the helium gas mixture but less than $30\mathrm{\mu }\mathrm{m}$ in the pure neon. At the shock front, calculations show that the electrons are strongly collisional with mean-free path two orders of magnitude shorter than the characteristic length of the shock. However, the helium ions can reach a kinetic regime as a consequence of their mean-free path being comparable to the scale of the shock. A radiation-hydrodynamic simulation demonstrates the impact of thermal conduction on the formation of the precursors with charge state, Z, playing a major role in heat flow and the precursor formation in both the helium mixture and the pure neon gases. Particle-in-cell simulations are also performed to study the ion kinetic effects on the formation of the observed precursors. A group of fast-streaming ions is observed leading the shock only in the helium gas mixture. Both effects explain the longer precursor layer in the helium shock.

Magnetized target fusion approach to inertial confinement fusion involves the formation of strong shocks that travel along a magnetized plasma. Shocks, which play a dominant role in thermalizing the upstream kinetic energy generated in the implosion stage, are seldom free from perturbations, and they wrinkle in response to upstream or downstream disturbances. In Z-pinch experiments, significant plasma instability mitigation was observed with pre-embedded axial magnetic fields. To isolate effects, in this work we theoretically study the impact of perpendicular magnetic fields on the planar shock dynamics for different equations of state. For fast magnetosonic shocks in ideal gases, it was found that the magnetic field amplifies the intensity of the perturbations when $\gamma >2$ or it weakens them when $\gamma <2$. Weak shocks have been found to be stable regardless of the magnetic plasma intensity and gas compressibility; however, for sufficiently strong shocks the magnetic fields can promote a neutral stability/SAE at the shock if the adiabatic index is higher than $1+\sqrt{2}$. Results have been validated with numerical simulations performed with the FLASH code.

We study the buckling of pressurized spherical shells by Monte Carlo simulations in which the detailed balance is explicitly broken—thereby driving the shell to be active, out of thermal equilibrium. Such a shell typically has either higher (active) or lower (sedate) fluctuations compared to one in thermal equilibrium depending on how the detailed balance is broken. We show that, for the same set of elastic parameters, a shell that is not buckled in thermal equilibrium can be buckled if turned active. Similarly a shell that is buckled in thermal equilibrium can unbuckle if sedated. Based on this result, we suggest that it is possible to experimentally design microscopic elastic shells whose buckling can be optically controlled.

We consider the minimal thermodynamic cost of an individual computation, where a single input $x$ is mapped to a single output $y$. In prior work, Zurek proposed that this cost was given by $K\left(x\right|y)$, the conditional Kolmogorov complexity of $x$ given $y$ (up to an additive constant that does not depend on $x$ or $y$). However, this result was derived from an informal argument, applied only to deterministic computations, and had an arbitrary dependence on the choice of protocol (via the additive constant). Here we use stochastic thermodynamics to derive a generalized version of Zurek's bound from a rigorous Hamiltonian formulation. Our bound applies to all quantum and classical processes, whether noisy or deterministic, and it explicitly captures the dependence on the protocol. We show that $K\left(x\right|y)$ is a minimal cost of mapping $x$ to $y$ that must be paid using some combination of heat, noise, and protocol complexity, implying a trade-off between these three resources. Our result is a kind of “algorithmic fluctuation theorem” with implications for the relationship between the second law and the Physical Church-Turing thesis.

We examine the quench dynamics of an extended Su-Schrieffer-Heeger (SSH) model involving long-range hopping that can hold multiple topological phases. Using winding number diagrams to characterize the system's topological phases geometrically, it is shown that there can be multiple winding number transition paths for a quench between two topological phases. The dependence of the quench dynamics is studied in terms of the survival probability of the fermionic edge modes and postquench transport. For two quench paths between two topological regimes with the same initial and final topological phase, the survival probability of edge states is shown to be strongly dependent on the winding number transition path. This dependence is explained using energy band diagrams corresponding to the paths. Following this, the effect of the winding number transition path on transport is investigated. We find that the velocities of maximum transport channels varied along the winding number transition path. This variation depends on the path we choose, i.e., it increases or decreases depending upon the path. An analysis of the coefficient maps, energy spectrum, and spatial structure of the edge states of the final quench Hamiltonian provides an understanding of the path-dependent velocity variation phenomenon.

Reported in the paper are results of unsteady three-dimensional direct numerical simulations of laminar and turbulent, lean hydrogen-air, complex-chemistry flames propagating in forced turbulence in a box. To explore the eventual influence of thermodiffusive instability of laminar flames on turbulent burning velocity, (i) a critical length scale ${\mathrm{\Lambda }}_n$ that bounds regimes of unstable and stable laminar combustion is numerically determined by gradually decreasing the width $\mathrm{\Lambda }$ of computational domain until a stable laminar flame is obtained, and (ii) simulations of turbulent flames are performed by varying the width from $\mathrm{\Lambda }<{\mathrm{\Lambda }}_n$ (in this case, the instability is suppressed) to $\mathrm{\Lambda }>{\mathrm{\Lambda }}_n$ (in this case, the instability may grow). Moreover, simulations are performed either using mixture-averaged transport properties (low Lewis number flames) or setting diffusivities of all species equal to heat diffusivity of the mixture (equidiffusive flames), with all other things being equal. Obtained results show a significant increase in turbulent burning velocity $U_T$ when the boundary $\mathrm{\Lambda }={\mathrm{\Lambda }}_n$ is crossed in weak turbulence, but almost equal values of $U_T$ are computed at $\mathrm{\Lambda }<{\mathrm{\Lambda }}_n$ and $\mathrm{\Lambda }>{\mathrm{\Lambda }}_n$ in moderately turbulent flames characterized by a Karlovitz number equal to 3.4 or larger. These results imply that thermodiffusive instability of laminar premixed flames substantially affects burning velocity in weak turbulence only, in line with a simple criterion proposed by Chomiak and Lipatnikov (Phys. Rev. E 107, 015102, (2023)).

We analyze the statistical properties of radiative transitions for a molecular system possessing discrete, equally spaced, energy levels, interacting with thermal radiation at constant temperature. A radiative fluctuation-dissipation theorem is derived and the particle velocity distribution analyzed. It is shown analytically that, neglecting molecular collisions, the velocity distribution function cannot be Gaussian, as the equilibrium value for the kurtosis $\kappa$ is different from $\kappa =3$. A Maxwellian velocity distribution can be recovered in the limit of small radiative friction.

Energy transport is a fundamental physical process that plays a prominent role in the function and performance of myriad systems and technologies. Recent experimental measurements have shown that subjecting a macroscale system to a time-periodic temperature gradient can increase thermal conductivity in comparison to a static temperature gradient. Here, we theoretically examine this mechanism in a nanoscale model by applying a stochastic Langevin framework to describe the energy transport properties of a particle connecting two heat baths with different temperatures, where the temperature difference between baths is oscillating in time. Analytical expressions for the energy flux of each heat bath and for the system itself are derived for the case of a free particle and a particle in a harmonic potential. We find that dynamical effects in the energy flux induced by temperature oscillations give rise to complex energy transport hysteresis effects. The presented results suggest that applying time-periodic temperature modulations is a potential route to control energy storage and release in molecular devices and nanosystems.

We study pairwise interactions between localized topological structures in chiral magnetic and cholesteric liquid crystal (CLC) systems confined in the planar geometry. Our calculations for magnetics are based on the lattice model that takes into account the bulk and surface anisotropies along with the exchange and the Dzyaloshinskii-Moriya interactions. In CLC cells, these anisotropies describe the energy of interaction with an external magnetic or electric field and the anchoring energy assuming that the magnetic or electric anisotropy is negative and the boundary conditions are homeotropic. We have selected the region of the phase diagram, where various localized solitonlike structures, including skyrmion tubes, torons, and leeches, embedded in the ground state of the $z$-cone (conical phase) coexist, and carried out numerical analysis of the distance dependencies of the effective intersoliton interaction potentials. For skyrmions and torons, the potentials are found to be attractive in the large separation region. It turned out that for these potentials, the effects of axial asymmetry are negligible. By contrast, it turned out that for the intermediate structures between the skyrmions and torons known as the leeches, the leech-leech potentials generally depend on the orientation of the intersoliton separation vector and their large distance parts may become repulsive at certain directions of the vector. All the potentials have the short distance repulsive parts and the local minima located at the equilibrium separations. It is found that the skyrmion-skyrmion potential has an additional metastable configuration shifted towards the short-distance region.

We consider pairs of Gaussian orthogonal ensemble matrices which are correlated with each other and subject to additive and multiplicative rank-one perturbations. We focus on the regime of parameters in which the finite-rank perturbations generate outliers in the spectrum of the matrices. We investigate the statistical correlation (i.e., the typical overlap) between the eigenvectors associated to the outlier eigenvalues of each matrix in the pair, as well as the typical overlap between the outlier eigenvector of one matrix with the eigenvectors in the bulk of the spectrum of the other matrix. We discuss implications of these results for the signal recovery problem for spiked matrices, as well as for problems of high-dimensional random landscapes.

In recent work [Rodrigues et al., Phys. Rev. E 100, 022121 (2019)], evidence was found that the surface adsorption transition of interacting self-avoiding trails (ISATs) placed on the square lattice displays a nonuniversal behavior at the special adsorption point (SAP) where the collapsing polymers adsorb. In fact, different surface exponents ${\varphi }^{\left(s\right)}$ and $1/{\delta }^{\left(s\right)}$ were found at the SAP depending on whether the surface orientation is horizontal (HS) or diagonal (DS). Here, we revisit these systems and study other ones, through extensive Monte Carlo simulations, considering much longer trails than previous works. Importantly, we demonstrate that the different exponents observed in the reference above are due to the presence of a surface-attached-globule (SAG) phase in the DS system, which changes the multicritical nature of the SAP and is absent in the HS case. By considering a modified horizontal surface (mHS), on which the trails are forbidden from having two consecutive steps, resembling the DS situation, a stable SAG phase is found in the phase diagram, and both DS and mHS systems present similar $1/{\delta }^{\left(s\right)}$ exponents at the SAP, namely, $1/{\delta }^{\left(s\right)}\approx 0.44$, whereas $1/{\delta }^{\left(s\right)}\approx 0.34$ in the HS case. Intriguingly, while ${\varphi }^{\left(s\right)}\approx 1/{\delta }^{\left(s\right)}$ is found for the DS and HS scenarios, as expected, in the mHS case ${\varphi }^{\left(s\right)}$ is about 10% smaller than $1/{\delta }^{\left(s\right)}$. These results strongly indicate that at least two universality classes exist for the SAPs of adsorbing ISATs on the square lattice.

A one-dimensional nonlinear dynamical system of gauge-coupled intrasite excitations and lattice vibrations on an infinite one-dimensional regular lattice is studied. The system as a whole is shown to be integrable in the Lax sense and it admits the exact four-component analytical solutions. Two mutually $\mathcal{PT}$-conjugated symmetry broken solutions are explicitly isolated in the framework of the Darboux-Bäcklund dressing technique. Each of the obtained four-component solutions demonstrates the pronounced interplay between the interacting subsystems in the form of an essentially nonlinear superposition of two principally distinct types of traveling waves characterized by two physically distinct spatial scales and by two distinct running velocities. Depending on the relationships between the spatial scaling parameters the system can manifest itself in three qualitatively distinct dynamical regimes referred to as the monopole regime, dipole regime, and threshold regime. The threshold value of the localization parameter separating the monopole and dipole dynamical regimes is strictly established in terms of basic physical parameters. The phenomenon of dipole-monopole crossover for the spatial distribution of pseudoexcitons is shown to initiate the partial splitting between the pseudoexcitonic and vibrational subsystems in the threshold dynamical regime specified by the threshold value of the localization parameter. This partial splitting is manifested by the complete elimination of one pseudoexcitonic component accompanied by the actual conversion of another pseudoexcitonic component into a pseudoexcitonic chargeless half mode. The integrable nonlinear pseudoexciton-phonon system under study is expected to be applicable for modeling the nonlinear dynamical properties of properly designed $\mathcal{PT}$-symmetric metamaterials.

The stochastic motion of a particle with long-range correlated increments (the moving phase) which is intermittently interrupted by immobilizations (the trapping phase) in a disordered medium is considered in the presence of an external drift. In particular, we consider trapping events whose times follow a scale-free distribution with diverging mean trapping time. We construct this process in terms of fractional Brownian motion with constant forcing in which the trapping effect is introduced by the subordination technique, connecting “operational time” with observable “real time.” We derive the statistical properties of this process such as non-Gaussianity and nonergodicity, for both ensemble and single-trajectory (time) averages. We demonstrate nice agreement with extensive simulations for the probability density function, skewness, kurtosis, as well as ensemble and time-averaged mean-squared displacements. We place a specific emphasis on the comparisons between the cases with and without drift.

In neuronal systems, inhibition contributes to stabilizing dynamics and regulating pattern formation. Through developing mean-field theories of neuronal models, using complete graph networks, inhibition is commonly viewed as one “control parameter” of the system, promoting an absorbing phase transition. Here, we show that, for low connectivity sparse networks, inhibition weight is not a control parameter of the absorbing transition. We present analytical and simulation results using generic stochastic integrate-and-fire neurons that, under specific restrictions, become other simpler stochastic neuron models common in literature, which allows us to show that our results are valid for those models as well. We also give a simple explanation about why the inhibition role depends on topology, even when the topology has a dimensionality greater than the critical one. The absorbing transition independence of the inhibitory weight may be an important feature of a sparse network, as it will allow the network to maintain a near-critical regime, self-tuning average excitation, but at the same time have the freedom to adjust inhibitory weights for computation, learning, and memory, exploiting the benefits of criticality.

We numerically study the dynamic state of a low-Reynolds-number turbulent channel flow from the viewpoints of symbolic dynamics and nonlinear forecasting. A low-dimensionally (high-dimensionally) chaotic state of the streamwise velocity fluctuations emerges at a viscous sublayer (logarithmic layer). The possible presence of the chaotic states is clearly identified by orbital instability-based nonlinear forecasting and ordinal partition transition network entropy in combination with the surrogate data method.

We study transition paths in energy landscapes over multicategorical Potts configurations using the mean-field approach introduced by Mauri et al. [Phys. Rev. Lett. 130, 158402 (2023)]. Paths interpolate between two fixed configurations or are anchored at one extremity only. We characterize the properties of “good” transition paths realizing a trade-off between exploring low-energy regions in the landscape and being not too long, such as their entropy or the probability of escape from a region of the landscape. We unveil the existence of a phase transition separating a regime in which paths are stretched in between their anchors from another regime where paths can explore the energy landscape more globally to minimize the energy. This phase transition is first illustrated and studied in detail on a mathematically tractable Hopfield-Potts toy model, then studied in energy landscapes inferred from protein sequence data.