Lipid-based dispersion of hydrophobic nanoparticles (NPs) not only gives fundamental insight into how nanomaterials distribute in live cells and organisms, but also provides a quite general route to designing nanocarrier agents in triggered drug delivery and medical imaging. It is not clearly understood how hydrophobic NPs arrange in the interior of a membrane. In this paper, with computer simulation techniques, we demonstrate that hydrophobic NPs having a diameter compared to the hydrophobic thickness of the membrane are capable of clustering in the hydrophobic interior of a cell membrane. Except from the isotropic aggregation, an unexpected linear arrangement of spherical NPs, which is still not found from experiments, is identified here. The free-energy costs associated with linear and isotropic aggregations are computed explicitly to interpret aggregation behavior and the obtained phase diagrams give us a comprehensive understanding of where linear aggregation is expected. In this work we also shows that NP size and membrane tension play key roles in determining the NP aggregate, while the effects of NP concentration and membrane curvature seem to be relatively weak.

A simple and reproducible analog experiment was used to simulate fracture formation in a low-permeability elastic solid during internal fluid/gas production, with the objective of developing a better understanding of the mechanisms that control the dynamics of fracturing, fracture opening and closing, and fluid transport. In the experiment, nucleation, propagation, and coalescence of fractures within an elastic gelatin matrix, confined in a Hele-Shaw cell, occurred due to ${\mathrm{CO}}_2$ production via fermentation of sugar, and it was monitored by optical means. We first quantified how a fracture network develops, and then how intermittent fluid transport is controlled by the dynamics of opening and closing of fractures. The gas escape dynamics exhibited three characteristic behaviors: (1) Quasiperiodic release of gas with a characteristic frequency that depends on the gas production rate but not on the system size. (2) A $1/f$ power spectrum for the fluctuations in the total open fracture area over an intermediate range of frequencies ($f$), which we attribute to collective effects caused by interaction between fractures in the drainage network. (3) A $1/f^2$ power spectrum was observed at high frequencies, which can be explained by the characteristic behavior of single fractures.

Synchronization is a ubiquitous phenomenon occurring in social, biological, and technological systems when the internal rythms of their constituents are adapted to be in unison as a result of their coupling. This natural tendency towards dynamical consensus has spurred a large body of theoretical and experimental research in recent decades. The Kuramoto model constitutes the most studied and paradigmatic framework in which to study synchronization. In particular, it shows how synchronization appears as a phase transition from a dynamically disordered state at some critical value for the coupling strength between the interacting units. The critical properties of the synchronization transition of this model have been widely studied and many variants of its formulations have been considered to address different physical realizations. However, the Kuramoto model has been studied only within the domain of classical dynamics, thus neglecting its applications for the study of quantum synchronization phenomena. Based on a system-bath approach and within the Feynman path-integral formalism, we derive equations for the Kuramoto model by taking into account the first quantum fluctuations. We also analyze its critical properties, the main result being the derivation of the value for the synchronization onset. This critical coupling increases its value as quantumness increases, as a consequence of the possibility of tunneling that quantum fluctuations provide.

Pyragas time-delayed feedback is a control scheme designed to stabilize unstable periodic orbits, which occur naturally in many nonlinear dynamical systems. It has been successfully implemented in a number of applications, including lasers and chemical systems. The control scheme targets a specific unstable periodic orbit by adding a feedback term with a delay chosen as the period of the unstable periodic orbit. However, in an experimental or industrial environment, obtaining the exact period or setting the delay equal to the exact period of the target periodic orbit may be difficult. This could be due to a number of factors, such as incomplete information on the system or the delay being set by inaccurate equipment. In this paper, we evaluate the effect of Pyragas control on the prototypical generic subcritical Hopf normal form when the delay is close to but not equal to the period of the target periodic orbit. Specifically, we consider two cases: first, a constant, and second, a linear approximation of the period. We compare these two cases to the case where the delay is set exactly to the target period, which serves as the benchmark case. For this comparison, we construct bifurcation diagrams and determine any regions where a stable periodic orbit close to the target is stabilized by the control scheme. In this way, we find that at least a linear approximation of the period is required for successful stabilization by Pyragas control.

Generic Hamiltonian systems have a mixed phase space where regions of regular and chaotic motion coexist. We present a method for constructing an integrable approximation to such regular phase-space regions including a nonlinear resonance chain. This approach generalizes the recently introduced iterative canonical transformation method. In the first step of the method a normal-form Hamiltonian with a resonance chain is adapted such that actions and frequencies match with those of the nonintegrable system. In the second step a sequence of canonical transformations is applied to the integrable approximation to match the shape of regular tori. We demonstrate the method for the generic standard map at various parameters.

We present a dynamically tunable mechanism of wave transmission in one-dimensional helicoidal phononic crystals in a shape similar to DNA structures. These helicoidal architectures allow slanted nonlinear contact among cylindrical constituents, and the relative torsional movements can dynamically tune the contact stiffness between neighboring cylinders. This results in cross-talking between in-plane torsional and out-of-plane longitudinal waves. We numerically demonstrate their versatile wave mixing and controllable dispersion behavior in both wavenumber and frequency domains. Based on this principle, a suggestion toward an acoustic configuration bearing parallels to a transistor is further proposed, in which longitudinal waves can be switched on and off through torsional waves.

We present a whole series of methods to alleviate the sign problem of the fermionic shadow wave function in the context of variational Monte Carlo. The effectiveness of our techniques is demonstrated on liquid ${}^3\mathrm{He}$. We found that although the variance is reduced, the gain in efficiency is restricted by the increased computational cost. Yet, this development not only extends the scope of the fermionic shadow wave function, but also facilitates highly accurate quantum Monte Carlo simulations previously thought not feasible.

A way to understand thermalization in an isolated system is to interpret it as an increase in entanglement between subsystems. Here we test this idea through a combination of analytical and Krylov-subspace-based numerical methods applied to a quantum gas of bosons. We find that the entanglement entropy of a subsystem is rapidly generated at the initial state of the evolution, to quickly approach the thermal value. Our results also provide an accurate numerical test of the eigenstate thermalization hypothesis (ETH), according to which a single energy eigenstate of an isolated system behaves in certain respects as a thermal state. In the context of quantum black holes, we propose that the ETH is a quantum version of the classical no-hair theorem.

Colloidal systems are considered good models of molecular glasses and we further explore the range of validity of this paradigm using a thermosensitive core-shell particle dispersion to study the aging response of a colloidal glass subsequent to both shear-melting and temperature (concentration)-jump perturbations in the vicinity of the glass transition concentration or temperature. Sequential creep experiments were used to probe the different aging responses of the system. The colloidal glass displays aging behavior after both types of perturbation and our results indicate that this colloidal glass is similar to a molecular glass, in that shift rates are found to be below unity and to decrease towards zero as the glass temperature (or concentration) is approached as temperature increases. However, the kinetics of the aging in the two cases are different indicating that the structural changes induced by the mechanical perturbation are different from those induced by the temperature or concentration jump—similar to findings on mechanical rejuvenation of molecular glasses. We also find differences between the colloidal glass and molecular glasses: In the case of the colloidal glass the structural recovery or equilibration times do not diverge, while the mechanical relaxation times do. On the other hand, for the molecular glass, both times change very rapidly with decreasing temperature, apparently towards a distant point of divergence.

Power-law distributions are typical macroscopic features occurring in almost all complex systems observable in nature. As a result, researchers in quantitative analyses must often generate random synthetic variates obeying power-law distributions. The task is usually performed through standard methods that map uniform random variates into the desired probability space. Whereas all these algorithms are theoretically solid, in this paper we show that they are subject to severe machine-dependent limitations. As a result, two dramatic consequences arise: (i) the sampling in the tail of the distribution is not random but deterministic; (ii) the moments of the sample distribution, which are theoretically expected to diverge as functions of the sample sizes, converge instead to finite values. We provide quantitative indications for the range of distribution parameters that can be safely handled by standard libraries used in computational analyses. Whereas our findings indicate possible reinterpretations of numerical results obtained through flawed sampling methodologies, they also pave the way for the search for a concrete solution to this central issue shared by all quantitative sciences dealing with complexity.

The mechanism of subnanosecond current front rise observed previously in the experiment in high-voltage pulse open discharge in helium is studied in kinetic particle-in-cell simulations. The Boltzmann equations for electrons, ions, and fast atoms are solved self-consistently with the Poisson equations for the electrical potential. The partial contributions to the secondary electron emission from the ions, fast atoms, photons, and electrons, bombarding the electrode, are calculated. In simulations, as in the experiment, the discharge glows between two symmetrical cathodes and the anode grid in the midplane at $P=6$ Torr and the applied voltage of 20 kV. The electron avalanche development is considered for two experimental situations during the last stage of breakdown: (i) with constant voltage and (ii) with decreasing voltage. For case (i), the subnanosecond current front rise is set by photons from the collisional excitation transfer reactions. For the case (ii), the energetic electrons swamp the cathode during voltage drop and provide the secondary electron emission for the subnanosecond current rise, observed in the experiment.

Stripe patterns, observed in a large variety of physical systems, often exhibit a slow nonequilibrium dynamics because ordering is impeded by the presence of topological defects. Using computer simulations based on a well-established model for stripe formation, we show that a slow dynamics and aging occur also in stripe patterns free of topological defects. For a wide range of noise strengths, the two-time orientation correlation function follows a scaling form that is typical for systems exhibiting a growing length scale. In our case, the underlying mechanism is the coarsening of orientation fluctuations, ultimately leading to power-law spatial correlations perpendicular to the stripes. Our results show that even for the smallest amount of noise, stripe phases without topological defects do not reach equilibrium. This constitutes an important aspect of the dynamics of modulated phases.

We address the problem of supratransmission of waves in a discrete nonlinear system, driven at one end by a periodic excitation at a frequency lying above the phonon band edge. In an experimental electrical transmission line made of 200 inductance-capacitance LC cells, we establish the existence of a voltage threshold for a supratransmission enabling the generation and propagation of cut-off solitons within the line. The decisive role of modulational instability in the onset and development of the process of generation of cut-off solitons is clearly highlighted. The phenomenon of dissipation is identified as being particularly harmful for the soliton generation, but we show that its impact can be managed by a proper choice of the amplitude of the voltage excitation of the system.

The scattering of kinks and low-frequency breathers of the nonlinear sine-Gordon (SG) equation on a spatially localized parity-time-symmetric perturbation (defect) with a balanced gain and loss is investigated numerically. It is demonstrated that if a kink passes the defect, it always restores its initial momentum and energy, and the only effect of the interaction with the defect is a phase shift of the kink. A kink approaching the defect from the gain side always passes, while in the opposite case it must have sufficiently large initial momentum to pass through the defect instead of being trapped in the loss region. The kink phase shift and critical velocity are calculated by means of the collective variable method. Kink-kink (kink-antikink) collisions at the defect are also briefly considered, showing how their pairwise repulsive (respectively, attractive) interaction can modify the collisional outcome of a single kink within the pair with the defect. For the breather, the result of its interaction with the defect depends strongly on the breather parameters (velocity, frequency, and initial phase) and on the defect parameters. The breather can gain some energy from the defect and as a result potentially even split into a kink-antikink pair, or it can lose a part of its energy. Interestingly, the breather translational mode is very weakly affected by the dissipative perturbation, so that a breather penetrates more easily through the defect when it comes from the lossy side, than a kink. In all studied soliton-defect interactions, the energy loss to radiation of small-amplitude extended waves is negligible.

The Kuramoto model with time-varying parameters has been extended to consider the effect of delay in couplings. A collective dynamics arises from the interplay between the time scales of the original system, the external forcing, and the delays. This complex low-dimensional dynamics is described, uncovering an echo effect near the synchronization threshold. Hence, the delayed couplings substantially alter the dynamics of what is an open system and should be taken into consideration, depending on the ensemble's evolution time scale. We also introduce a first-harmonic approximation for the evolution of the mean field under harmonic forcing, valid for any delays and forcing of a coherent population.

We study turbulent diffusion of chemically reacting gaseous admixtures in a developed turbulence. In our previous study [Phys. Rev. Lett. 80, 69 (1998)] using a path-integral approach for a delta-correlated in a time random velocity field, we demonstrated a strong modification of turbulent transport in fluid flows with chemical reactions or phase transitions. In the present study we use the spectral $\tau$ approximation that is valid for large Reynolds and Peclet numbers and show that turbulent diffusion of the reacting species can be strongly depleted by a large factor that is the ratio of turbulent and chemical times (turbulent Damköhler number). We have demonstrated that the derived theoretical dependence of a turbulent diffusion coefficient versus the turbulent Damköhler number is in good agreement with that obtained previously in the numerical modeling of a reactive front propagating in a turbulent flow and described by the Kolmogorov-Petrovskii-Piskunov-Fisher equation. We have found that turbulent cross-effects, e.g., turbulent mutual diffusion of gaseous admixtures and turbulent Dufour effect of the chemically reacting gaseous admixtures, are less sensitive to the values of stoichiometric coefficients. The mechanisms of the turbulent cross-effects differ from the molecular cross-effects known in irreversible thermodynamics. In a fully developed turbulence and at large Peclet numbers the turbulent cross-effects are much larger than the molecular ones. The obtained results are applicable also to heterogeneous phase transitions.

The structure and thermodynamics of strongly coupled dusty plasmas are investigated with the soft mean spherical approximation. This integral theory approach is analytically solvable for Yukawa pair interactions yielding a closed-form solution for the direct correlation function. The pair correlation function, the structure factor, and basic thermodynamic quantities are calculated for a wide range of parameters. Exact consistency between the “energy”-“virial” thermodynamic routes and approximate consistency between the “energy”-“compressibility” paths is demonstrated. Comparison with extensive molecular dynamics results is carried out and a remarkable agreement from the Coulomb limit to the strongly screened limit is revealed. The soft mean spherical approximation is concluded to be particularly well suited for the study of dusty plasma liquids, uniquely combining simplicity and accuracy.

In this paper we investigate the implementation of contact angles in the pseudopotential lattice Boltzmann modeling of wetting at a large density ratio ${\rho }_L/{\rho }_V=500$. The pseudopotential lattice Boltzmann model [X. Shan and H. Chen, Phys. Rev. E 49, 2941 (1994)] is a popular mesoscopic model for simulating multiphase flows and interfacial dynamics. In this model the contact angle is usually realized by a fluid-solid interaction. Two widely used fluid-solid interactions, the density-based interaction and the pseudopotential-based interaction, as well as a modified pseudopotential-based interaction formulated in the present paper are numerically investigated and compared in terms of the achievable contact angles, the maximum and the minimum densities, and the spurious currents. It is found that the pseudopotential-based interaction works well for simulating small static (liquid) contact angles $\theta <{90}^{\circ }$, however, it is unable to reproduce static contact angles close to ${180}^{\circ }$. Meanwhile, it is found that the proposed modified pseudopotential-based interaction performs better in light of the maximum and the minimum densities and is overall more suitable for simulating large contact angles $\theta >{90}^{\circ }$ as compared with the two other types of fluid-solid interactions. Furthermore, the spurious currents are found to be enlarged when the fluid-solid interaction force is introduced. Increasing the kinematic viscosity ratio between the vapor and liquid phases is shown to be capable of reducing the spurious currents caused by the fluid-solid interactions.

Hawking radiation has become experimentally testable thanks to the many analog systems which mimic the effects of the event horizon on wave propagation. These systems are typically dominated by dispersion and give rise to a numerically soluble and stable ordinary differential equation only if the rest-frame dispersion relation ${\Omega }^2\left(k\right)$ is a polynomial of relatively low degree. Here we present a new method for the calculation of wave scattering in a one-dimensional medium of arbitrary dispersion. It views the wave equation as an integral equation in Fourier space, which can be solved using standard and efficient numerical techniques.

Analog gravity experiments make feasible the realization of black hole space-times in a laboratory setting and the observational verification of Hawking radiation. Since such analog systems are typically dominated by dispersion, efficient techniques for calculating the predicted Hawking spectrum in the presence of strong dispersion are required. In the preceding paper, an integral method in Fourier space is proposed for stationary $1+1$-dimensional backgrounds which are asymptotically symmetric. Here, this method is generalized to backgrounds which are different in the asymptotic regions to the left and right of the scattering region.