In social networks, the balance theory has been studied by considering either the triple interactions between the links (structural balance) or the triple interaction of nodes and links (coevolutionary balance). In the structural balance theory, the links are not independent from each other, implying a global effect of this term and it leads to a discontinuous phase transition in the system's balanced states as a function of temperature. However, in the coevolutionary balance the links only connect two local nodes and a continuous phase transition emerges. In this paper, we consider a combination of both to understand which of these types of interactions will identify the stability of the network. We are interested to see how adjusting the robustness of each term versus the other might affect the system to reach a balanced state. We use statistical mechanics methods and the mean-field theory and also the Monte Carlo numerical simulations to investigate the behavior of the order parameters and the total energy of the system. We find the phase diagram of the system which demonstrates the competition of these two terms at different ratios against each other and different temperatures. The system shows a tricritical point above which the phase transition switches from continuous to discrete. Also the superiority of the local perspective is observed at low temperatures and the global view will be the dominant term in determining the stability of the system at higher temperatures.

Molecular interactions between liquid crystals (LCs) and reactive mesogens (RMs) at temperatures across the phase transition regions were comprehensively studied during photopolymerization-induced phase separation (PPIPS) beginning with raw mixtures until the formation of polymer network liquid crystals (PNLCs). Then, the molecules were found to be nonuniformly more and less mobile in response to temperature as PPIPS progressed. Optical birefringence and infrared absorption were carefully measured throughout PPIPS, using 4-cyano-4′-hexylbiphenyl (6CB) and 1,4bis-[4-(3-acryloyloxypropyloxy) benzoyloxy]-2-methylbenzene (RM257) as typical LCs and RMs. Microscopic views of thermoresponsive changes in the molecular orientation order of both LCs and RMs were obtained: LCs and RMs in raw mixtures interacted with one another but uniformly transformed their molecular orientation. Such interactions continuously change to become nonuniform with progress in PPIPS. At the incipient stages of PPIPS, RMs, which are polymerized but not completely networked, inhibit LCs from changing their molecular orientation and vice versa. As PPIPS progresses, some LCs become more mobile and some less mobile owing to RM constraints. The domain configuration of the submicrometer phase separation affects the thermoresponsive mobility of LCs and RMs, that is, LCs become more mobile in LC-richer areas. The quantitative knowledge here provides comprehensive insight that LCs and RMs are mutually constrained and that such interactive behavior varies nonuniformly as PPIPS progresses.

A wide-range (0 to 1044.0 g/${\mathrm{cm}}^3$ and 0 to ${10}^9$ K) equation-of-state (EOS) table for a ${\mathrm{CH}}_{1.72}{\mathrm{O}}_{0.37}{\mathrm{N}}_{0.086}$ quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. The simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

We present a scheme for amplifying an extreme-ultraviolet (XUV) seed isolated attosecond pulse via stimulated Raman scattering of a pulse-train pump. At sufficient seed and pump intensity, the amplification is nonlinear, and the amplitude of the seed pulse can reach that of the pump, one order of magnitude higher than the initial seed amplitude. In the linear amplification regime, we find that the spectral signature of the pump pulse train is imprinted on the spectrum of the amplified seed pulse. Since the spectral signature is imprinted with its frequency downshifted by the plasma frequency, it is possible to deduce the electron density in the region of interaction. This region can be of micrometer length scale longitudinally. By varying the delay between the seed and the pump, this scheme provides a local electron-density measurement inside solid-density plasmas that cannot be probed with optical frequencies, with micrometer resolution.

Large deviations in chaotic dynamics have potentially significant and dramatic consequences. We study large deviations of series of finite lengths $N$ generated by chaotic maps. The distributions generally display an exponential decay with $N$, associated with large-deviation (rate) functions. We obtain the exact rate functions analytically for the doubling, tent, and logistic maps. For the latter two, the solution is given as a power series whose coefficients can be systematically calculated to any order. We also obtain the rate function for the cat map numerically, uncovering strong evidence for the existence of a remarkable singularity of it that we interpret as a second-order dynamical phase transition. Furthermore, we develop a numerical tool for efficiently simulating atypical realizations of sequences if the chaotic map is not invertible, and we apply it to the tent and logistic maps.

We study the statistics of the first-passage time of a single run-and-tumble particle (RTP) in one spatial dimension, with or without resetting, to a fixed target located at $L>0$. First, we compute the first-passage time distribution of a free RTP, without resetting or in a confining potential, but averaged over the initial position drawn from an arbitrary distribution $p\left(x\right)$. Recent experiments used a noninstantaneous resetting protocol that motivated us to study in particular the case where $p\left(x\right)$ corresponds to the stationary non-Boltzmann distribution of an RTP in the presence of a harmonic trap. This distribution $p\left(x\right)$ is characterized by a parameter $\nu >0$, which depends on the microscopic parameters of the RTP dynamics. We show that the first-passage time distribution of the free RTP, drawn from this initial distribution, develops interesting singular behaviors, depending on the value of $\nu$. We then switch on resetting, mimicked by relaxation of the RTP in the presence of a harmonic trap. Resetting leads to a finite mean first-passage time and we study this as a function of the resetting rate for different values of the parameters $\nu$ and $b=L/c$, where $c$ is the position of the right edge of the initial distribution $p\left(x\right)$. In the diffusive limit of the RTP dynamics, we find a rich phase diagram in the $(b,\nu )$ plane, with an interesting reentrance phase transition. Away from the diffusive limit, qualitatively similar rich behaviors emerge for the full RTP dynamics.

The determination of the ionization of a system in the hot dense regime is a long standing issue. Recent studies have shown inconsistencies between standard predictions using average atom models and evaluations deduced from electronic transport properties computed with quantum molecular dynamics simulations [Bethkenhagen et al., Phys. Rev. Res. 2, 023260 (2020)]. Here, we propose a definition of the ionization based on its effect on the plasma structure as given by the pair distribution function (PDF), and on the concept of effective one-component plasma (eOCP). We also introduce a definition based on the total pressure and on a modelization of the electronic pressure. We show the equivalence of these definitions on two studies of carbon along the 100 eV isotherm and the 10 g/${\mathrm{cm}}^3$ isochor. Simulations along the 100 eV isotherm are obtained with the newly implemented Ext. First principles molecular dynamics (Fpmd) method in Abinit for densities ranging from 1 to 500 g/cm${}^3\mathrm{and}$ along the 10 g/${\mathrm{cm}}^3$ isochor with the recently published Spectral quadrature DFT (Sqdft) simulations, between 8 and 860 eV. The resulting ionizations are compared to the predictions of the average-atom code Qaam which is based on the muffin-tin approximation. A disagreement between the eOCP and the actual PDFs (non-OCP behavior) is interpreted as the onset of bonding in the system.

The work is devoted to experimental study of soft x-ray radiation spectra of laser-produced plasma in a wide spectral range of $5–100\AA$ at the Kamerton facility (GPI) with pulse duration 70 ps, pulse energy $1–5\mathrm{J}$, and wavelength 0.53 microns at which the laser intensity (power flux density) of $7\times {10}^{14}–3.5\times {10}^{15}\mathrm{W}/\mathrm{c}{\mathrm{m}}^2$ was achieved. A spectrograph was used, which had transmission diffraction gratings with a ratio of the elementary gap to the period of the structure of 0.25 and 0.41. Detection was performed on both UV-4 x-ray photo film and Fuji TR fluorescent imaging plate. Solid samples from Al, Si, Ti, Cu, Ta, and W were used as targets. The ionization states of the plasma corresponding to various electron temperatures were calculated, which made it possible to estimate the electron temperature by comparison of these calculation results with the experimentally obtained spectra. The estimated electron temperature, which depends on the laser pulse energy and the target material, varied within the range $100–450\mathrm{eV}$. To verify the correctness of the temperature estimations obtained by such comparison a numerical simulation of plasma radiation was carried out by the use of prismspect computer program. It was found that the results of this simulation are in good agreement with estimations on the base of experimentally obtained spectra. The analysis of these spectra showed that tantalum, tungsten, or titanium targets are the best candidates among the tested ones for the use of laser-produced plasma as a radiation source in the “water window” spectral range ($23–44\AA$) for applications in biology and medicine.

The Zeeman effect is used for measurement of magnetic fields in astrophysical and laboratory plasmas. Magnetic fields in atmospheres of magnetic white dwarf stars are in the range 40 kG–1 GG. The quadratic Zeeman effect results in the additional split and shift of lines for magnetic fields >2 MG. Hydrogen Balmer lines were studied in magnetic fields delivered by a 1 MA pulse power generator. The magnetic field was generated by rod loads 0.8–1 mm in diameter. A droplet of CH oil on the load center was the source of hydrogen. A low ionized oil layer was backlit by blackbody emission from the rod with a temperature of 0.5–0.6 eV. Zeeman splitting of $H$-alpha and $H$-beta absorption lines were with a grating spectrometer. A spectral shift of the central component of the triplet indicated the quadratic Zeeman effect in hydrogen lines.

Random perturbations applied in tandem to an ensemble of oscillating objects can synchronize their motion. We study multiple copies of an arbitrary dynamical system in a stable limit cycle, described via a standard phase reduction picture. The copies differ only in their arbitrary phases $\varphi$. Weak, randomly timed external impulses applied to all the copies can synchronize these phases over time. Beyond a threshold strength there is no such convergence to a common phase. Instead, the synchronization becomes erratic: successive impulses produce stochastic fluctuations in the phase distribution $q\left(\varphi \right)$, ranging from near-perfect to near-random synchronization. Here we show that the sampled entropies of these phase distributions themselves form a steady-state ensemble, whose average can be made arbitrarily negative by tuning the impulse strength. A random-walk description of the entropy's evolution accounts for the observed exponential distribution of entropies and for the stochastic synchronization phenomenon.

We consider coupled systems of advection-diffusion equations with initial and boundary conditions and determine conditions on the advection terms that allow us to obtain solutions that can be explicitly bounded above and below using the initial and boundary conditions. Given the advection terms, using our methodology one can easily check if such bounds can be obtained and then one can construct the necessary nonlinear transformation to allow the bounds to be determined. We apply this technique to determine bounding quantities for a number of examples. In particular, we show that the three-ion electroneutral Poisson-Nernst-Planck system of equations can be transformed into a system, which allows for the use of our techniques and we determine the bounding quantities. In addition, we determine the general form of advection terms that allow these techniques to be applied and show that our method can be applied to a very wide class of advection-diffusion equations.

We study a two-state model of attractive hard rods using the transfer matrix method, where the centers of the particles are confined to a straight line, but the orientations of the rods can be parallel or perpendicular to the confining line. The rods are modeled as hard rectangles with length $L$ and width $D$ and decorated with attractive sites at both ends of the rectangles. We find that the particles align parallel to the line and form long chains at low densities, while they turn out of the line and form a Tonks gas at high densities. With increasing the stickiness between the rods, the structural change between parallel and perpendicular states becomes stronger and the pressure vs density curve becomes almost a horizontal line at the transition pressure. We show that such a behavior is reminiscent of the first-order phase transition. This manifests in the validity of the lever rule of the phase transitions for very sticky cases.

We study a model of aligning self-propelled disks that nonreciprocally reorient the self-propulsion directions along the interparticle separation and towards the other disks. In the limit of small inertia and large softness, where conventional motility-induced phase separation is absent, we demonstrate that the homogeneous system at a small area fraction phase-separates into clusters and a low-density phase that, eventually, reenters the homogeneous phase with a monotonic increase in alignment strength. The disks inside the clusters move with a finite space-dependent speed, constantly shuttling between clusters through the surrounding low-density homogeneous phase while maintaining the hexatic structure properties within the clusters. The area fraction gradually increases from the periphery towards the center of the clusters with a negligible correlation of the velocity and propulsion direction inside the clusters. The novel collective behavior of reentrant phase separation is found to follow from both the limits of hard disks and extremely small inertia, tending towards the overdamped limit. However, important differences in the structural and dynamical properties are shown in the limit of hard disks and extremely small inertia, as compared to that for soft disks at finite inertia. We show that the cluster phase is associated with an effective temperature for a wide range of values of alignment strength, whereas an effective temperature is associated with the specific range of alignment in the low-density phase. We believe that the reentrant phase behavior in the limit of small area fraction and the remarkable properties of the clusters should be useful in understanding a wide range of physics issues, ranging from clogging and unclogging to information exchange and transport, in biological and synthetic self-propelled systems.

In computational models of particle packings with periodic boundary conditions, it is assumed that the packing is attached to exact copies of itself in all possible directions. The periodicity of the boundary then requires that all of the particles' images move together. An infinitely repeated structure, on the other hand, does not necessarily have this constraint. As a consequence, a jammed packing (or a rigid elastic network) under periodic boundary conditions may have a corresponding infinitely repeated lattice representation that is not rigid or indeed may not even be at a local energy minimum. In this manuscript, we prove this claim and discuss ways in which periodic boundary conditions succeed in capturing the physics of repeated structures and where they fall short.

Transducer arrays are a versatile tool for the contactless manipulation of spherical Rayleigh objects. Here we propose an analytical model for stable levitation of axisymmetric Mie objects through directly engineering the desired radiation force and torque. Acoustic contributions from multiple transducers are superimposed through the translation addition theorem, and the nonspherical objects are mapped into a sphere using the conformal transformation technique so that the scattered field can be asymptotically obtained. Then we give the acoustic radiation force and torque applied to a rigid nonspherical Mie object, which can be reconstructed as a series of quasiexplicit functions of the transducer (amplitude and phase) parameters. Through specifying the desired radiation force and torque exerted on the objects, a system of nonlinear equations is produced, which could be iteratively solved to retrieve appropriate transducer parameters that stabilize the object in an equilibrium position. Practically, we demonstrate several examples of stable levitation of a sphere, a spheroid, and a disk with an averaged radius of $a=7\mathrm{mm}$ (size parameter of $ka\approx 5.18$) above a transducer array. The absolute acoustic pressure field surrounding the objects simulated by the finite-element method is illustrated to verify the trapping results. The developed analytical model provides an alternative approach to retrieve the transducer parameters for levitating macroscopic nonspherical rigid objects, which may help design the systematic dynamical manipulation of Mie particles.

Perfectly matched layers (PMLs) are widely used in particle-in-cell simulations, in order to absorb electromagnetic waves that propagate out of the simulation domain. However, when charged particles cross the interface between the simulation domain and the PMLs, a number of numerical artifacts can arise. In order to mitigate these artifacts, we introduce a PML algorithm whereby the current deposited by the macroparticles in the PML is damped by an analytically derived optimal coefficient. The benefits of this algorithm are illustrated in practical simulations. In particular, it is shown that this algorithm is well suited for particles exiting the box in near-normal incidence, in the sense that the fields behave as if the exiting particle is propagating in an infinite vacuum.

The basic Lipkin-Meshkov-Glick model displays a second-order ground-state quantum phase transition and an excited-state quantum phase transition (ESQPT). The inclusion of an anharmonic term in the Hamiltonian implies a second ESQPT of a different nature. We characterize this ESQPT using the mean field limit of the model. The alternative ESQPT, associated with the changes in the boundary of the finite Hilbert space of the system, can be properly described using the order parameter of the ground-state quantum phase transition, the energy gap between adjacent states, the participation ratio, and the quantum fidelity susceptibility.

A particle diffusing in a two-dimensional channel of varying width $h\left(x\right)$ is considered. It is driven by a force of constant magnitude $f$, but random orientation across the channel. We suggest the projection technique to study the ratchet effect appearing in this system. Reducing the transverse coordinate, as well as the orientation of the force in the full-dimensional Fokker-Planck equation, we arrive at the generalized Fick-Jacobs equation, describing dynamics of the system in the longitudinal coordinate $x$ only. The additional effective potential $-\gamma \left(x\right)$, calculated within the mapping procedure, exhibits an increasing or decreasing part in the channel shaped by an asymmetric periodic $h\left(x\right)$, which determines the appearing ratchet current. As shown on a specific example, random driving in the transverse direction is much more effective than that in the longitudinal direction, at least for quickly flipping orientation of the force. Also, the transverse and the longitudinal driving push the rectified current in opposite directions along the same channel.

Neuromorphic computing exploits the dynamical analogy between many physical systems and neuron biophysics. Superconductor systems, in particular, are excellent candidates for neuromorphic devices due to their capacity to operate at great speeds and with low energy dissipation compared to their silicon counterparts. In this paper, we revisit a prior work on Josephson Junction-based neurons to identify the exact dynamical mechanisms underlying the system's neuronlike properties and reveal complex behaviors which are relevant for neurocomputation and the design of superconducting neuromorphic devices. Our paper lies at the intersection of superconducting physics and theoretical neuroscience, both viewed under a common framework—that of nonlinear dynamics theory.