Ion-beam cancer therapy is a promising technique to treat deep-seated tumors; however, for an accurate treatment planning, the energy deposition by the ions must be well known both in soft and hard human tissues. Although the energy loss of ions in water and other organic and biological materials is fairly well known, scarce information is available for the hard tissues (i.e., bone), for which the current stopping power information relies on the application of simple additivity rules to atomic data. Especially, more knowledge is needed for the main constituent of human bone, calcium hydroxyapatite (HAp), which constitutes 58% of its mass composition. In this work the energy loss of H and He ion beams in HAp films has been obtained experimentally. The experiments have been performed using the Rutherford backscattering technique in an energy range of 450–2000 keV for H and 400–5000 keV for He ions. These measurements are used as a benchmark for theoretical calculations (stopping power and mean excitation energy) based on the dielectric formalism together with the MELF-GOS (Mermin energy loss function-generalized oscillator strength) method to describe the electronic excitation spectrum of HAp. The stopping power calculations are in good agreement with the experiments. Even though these experimental data are obtained for low projectile energies compared with the ones used in hadron therapy, they validate the mean excitation energy obtained theoretically, which is the fundamental quantity to accurately assess energy deposition and depth-dose curves of ion beams at clinically relevant high energies. The effect of the mean excitation energy choice on the depth-dose profile is discussed on the basis of detailed simulations. Finally, implications of the present work on the energy loss of charged particles in human cortical bone are remarked.

The vibrational dynamics of the LiNC/LiCN molecular system is examined making use of the Riemannian geometry. Stability and chaoticity are analyzed, in this context, by means of the Jacobi-Levi-Civita equations, derived from the Jacobi metric, and its solutions. A dynamical indicator, called stability geometrical indicator, is introduced in order to ascertain the dynamical characteristics of stability and chaos in the molecule under study.

For generic 4D symplectic maps we propose the use of 3D phase-space slices, which allow for the global visualization of the geometrical organization and coexistence of regular and chaotic motion. As an example, we consider two coupled standard maps. The advantages of the 3D phase-space slices are presented in comparison to standard methods, such as 3D projections of orbits, the frequency analysis, and a chaos indicator. Quantum mechanically, the 3D phase-space slices allow for the comparison of Husimi functions of eigenstates of 4D maps with classical phase-space structures. This confirms the semiclassical eigenfunction hypothesis for 4D maps.

In this paper we demonstrate that two commonly used phenomenological post-Markovian quantum master equations can be derived without using any perturbative approximation. A system coupled to an environment characterized by self-classical configurational fluctuations, the latter obeying a Markovian dynamics, defines the underlying physical model. Both Shabani-Lidar equation [A. Shabani and D. A. Lidar, Phys. Rev. A 71, 020101(R) (2005)] and its associated approximated integrodifferential kernel master equation are obtained by tracing out two different bipartite Markovian Lindblad dynamics where the environment fluctuations are taken into account by an ancilla system. Furthermore, conditions under which the non-Markovian system dynamics can be unraveled in terms of an ensemble of measurement trajectories are found. In addition, a non-Markovian quantum jump approach is formulated. Contrary to recent analysis [L. Mazzola, E. M. Laine, H. P. Breuer, S. Maniscalco, and J. Piilo, Phys. Rev. A 81, 062120 (2010)], we also demonstrate that these master equations, even with exponential memory functions, may lead to non-Markovian effects such as an environment-to-system backflow of information if the Hamiltonian system does not commutate with the dissipative dynamics.

The percolation threshold for flow or conduction through voids surrounding randomly placed spheres is calculated. With large-scale Monte Carlo simulations, we give a rigorous continuum treatment to the geometry of the impenetrable spheres and the spaces between them. To properly exploit finite-size scaling, we examine multiple systems of differing sizes, with suitable averaging over disorder, and extrapolate to the thermodynamic limit. An order parameter based on the statistical sampling of stochastically driven dynamical excursions and amenable to finite-size scaling analysis is defined, calculated for various system sizes, and used to determine the critical volume fraction ${\varphi }_c=0.0317\pm 0.0004$ and the correlation length exponent $\nu =0.92\pm 0.05$.

We consider the mean territory covered by a particle that performs surface-mediated diffusion inside a spherical confining domain (in two and three dimensions) before exit through an opening on the surface. This quantity can be expressed in terms of the splitting probability between two targets on the surface. We derive a general formula that relates this splitting probability to the mean first passage time to a single target that has been recently calculated for such a surface-mediated diffusion process. This formula is exact for pointlike targets and is shown to be accurate for extended targets. The mean covered territory is then found and analyzed for an arbitrary extension of the exit region in both two- and three-dimensional spherical domains.

We extend a previously proposed field-theoretic self-consistent perturbation approach for the equilibrium dynamics of the Dean-Kawasaki equation presented in [Kim and Kawasaki, J. Stat. Mech. (2008) P02004]. By taking terms missing in the latter analysis into account we arrive at a set of three new equations for correlation functions of the system. These correlations involve the density and its logarithm as local observables. Our new one-loop equations, which must carefully deal with the noninteracting Brownian gas theory, are more general than the historic mode-coupling one in that a further approximation corresponding to Gaussian density fluctuations leads back to the original mode-coupling equation for the density correlations alone. However, without performing any further approximation step, our set of three equations does not feature any ergodic-nonergodic transition, as opposed to the historical mode-coupling approach.

We invoke mean-field density functional theory (DFT) to investigate the phase behavior of an amphiphilic fluid composed of a hard-sphere core plus a superimposed anisometric Lennard-Jones perturbation. The orientation dependence of the interactions consists of a contribution analogous to the interaction potential between a pair of “spins” in the classical, three-dimensional Heisenberg fluid and another one reminiscent of the interaction between (electric or magnetic) point dipoles. At fixed orientation both contributions are short-range in nature decaying as $r^{-6}$ ($r$ being the separation between the centers of mass of a pair of amphiphiles). Based upon two mean-field-like approximations for the pair correlation function that differ in the degree of sophistication we derive expressions for the phase boundaries between various isotropic and polar phases that we solve numerically by the Newton-Raphson method. For sufficiently strong coupling between the Heisenberg “spins” both mean-field approximations generate three topologically different and generic types of phase diagrams that are observed in agreement with earlier work [see, for example, Tavares et al., Phys. Rev. E 52, 1915 (1995)]. Whereas the dipolar contribution alone is incapable of stabilizing polar phases on account of its short-range nature it is nevertheless important for details of the phase diagram such as location of the gas-isotropic liquid critical point, triple, and tricritical points. By tuning the dipolar coupling constant suitably one may, in fact, switch between topologically different phase diagrams. Employing also Monte Carlo simulations in the isothermal-isobaric ensemble the general topology of the DFT phase diagrams is confirmed.

We develop a multicomponent phase-field model specially formulated to robustly simulate concentration variations from molar to atomic magnitudes across an interface, i.e., partition coefficients in excess of ${10}^{\pm 23}$ such as may be the case with species which are predominant in one phase and insoluble in the other. Substitutional interdiffusion on a normal lattice and concurrent interstitial diffusion are included. The composition in the interface follows the approach of Kim, Kim, and Suzuki [Phys. Rev. E 60, 7186 (1999)] and is compared to that of Wheeler, Boettinger, and McFadden [Phys. Rev. A 45, 7424 (1992)] in the context of large partitioning. The model successfully reproduces analytical solutions for binary diffusion couples and solute trapping for the demonstrated cases of extremely large partitioning.

Under the effect of strong genetic drift, it is highly probable to observe gene fixation or gene loss in a population, shown by singular peaks on a potential landscape. The genetic drift-induced noise gives rise to two-time-scale diffusion dynamics on the bipeaked landscape. We find that the logarithmically divergent (singular) peaks do not necessarily imply infinite escape times or biological fixations by iterating the Wright-Fisher model and approximating the average escape time. Our analytical results under weak mutation and weak selection extend Kramers's escape time formula to models with B (Beta) function-like equilibrium distributions and overcome constraints in previous methods. The constructed landscape provides a coherent description for the bistable system, supports the quantitative analysis of bipeaked dynamics, and generates mathematical insights for understanding the boundary behaviors of the diffusion model.

We numerically show that extreme events induced by parameter mismatches or noise in coupled oscillatory systems can be anticipated and suppressed before they actually occur. We show this in a main system unidirectionally coupled to an auxiliary system subject to a negative delayed feedback. Each system consists of two electronic oscillators coupled in a master-slave configuration. Extreme events are observed in this coupled system as large and sporadic desynchronization events. Under certain conditions, the auxiliary system can predict the dynamics of the main system. We use this to efficiently suppress the extreme events by applying a direct corrective reset to the main system.

Until recently, deterministic nonautonomous oscillatory systems with stable amplitudes and time-varying frequencies were not recognized as such and have often been mistreated as stochastic. These systems, named chronotaxic, were introduced in Phys. Rev. Lett. 111, 024101 (2013). In contrast to conventional limit cycle models of self-sustained oscillators, these systems posses a time-dependent point attractor or steady state. This allows oscillations with time-varying frequencies to resist perturbations, a phenomenon which is ubiquitous in living systems. In this work a detailed theory of chronotaxic systems is presented, specifically in the case of separable amplitude and phase dynamics. The theory is extended by the introduction of chronotaxic amplitude dynamics. The wide applicability of chronotaxic systems to a range of fields from biological and condensed matter systems to robotics and control theory is discussed.

We show that macroscopic nonintegrable lattices of spins 1/2, which are often considered to be chaotic, do not exhibit the basic property of classical chaotic systems, namely, exponential sensitivity to small perturbations. We compare chaotic lattices of classical spins and nonintegrable lattices of spins 1/2 in terms of their magnetization responses to an imperfect reversal of spin dynamics known as Loschmidt echo. In the classical case, magnetization is exponentially sensitive to small perturbations with a characteristic exponent equal to twice the value of the largest Lyapunov exponent of the system. In the case of spins 1/2, magnetization is only power-law sensitive to small perturbations.

The structure of a shock wave in a rarefied polyatomic gas is studied on the basis of the theory of extended thermodynamics. Three types of the shock wave structure observed in experiments, that is, the nearly symmetric shock wave structure (type A, small Mach number), the asymmetric structure (type B, moderate Mach number), and the structure composed of thin and thick layers (type C, large Mach number), are explained by the theory in a unified way. The theoretical prediction of the profile of the mass density agrees well with the experimental data. The well-known Bethe-Teller theory of the shock wave structure in a polyatomic gas is reexamined in the light of the present theory.

In a plane Couette cell a thin fluid layer consisting of water is sheared between the sides of a transparent band at Reynolds numbers ranging from 300 to 1400. The length of the cell's flow channel is large compared to the film separation. To extract the flow velocity in the experiments, a correlation image velocimetry method is used on pictures recorded with a high-speed camera. The flow is recorded at a resolution that allows us to analyze flow patterns similar in size to the film separation. The fluid flow is then studied by calculating flow velocity autocorrelation functions. The turbulent patterns that arise on this scale above a critical Reynolds number of $\mathrm{Re}=360$ display characteristic patterns that are proven by use of the calculated velocity autocorrelation functions. The patterns are metastable and reappear at different positions and times throughout the experiments. Typically these patterns are turbulent rolls which are elongated in the stream direction, which is the direction in which the band is moving. Although the flow states are metastable they possess similarities to the steady Taylor vortices known to appear in circular Taylor Couette cells.

The collision between two counterpropagating dust acoustic solitary waves in a strongly coupled dusty plasma has been observed. The measured velocity and width of the solitary wave agree with the solution of the Korteweg–de Vries equation derived by using the generalized hydrodynamic model. The two counterpropagating solitary waves of equal amplitude merge into a single pulse with twice the individual soliton amplitude and then pass through each other. The solitons suffer a small time delay in propagation after collision. The measured delay time obtained from their trajectories is also presented.

We develop a general-purpose formulation, based on two-dimensional spectral integrals, for computing electromagnetic fields produced by arbitrarily oriented dipoles in planar-stratified environments, where each layer may exhibit arbitrary and independent anisotropy in both its (complex) permittivity and permeability tensors. Among the salient features of our formulation are (i) computation of eigenmodes (characteristic plane waves) supported in arbitrarily anisotropic media in a numerically robust fashion, (ii) implementation of an $hp$-adaptive refinement for the numerical integration to evaluate the radiation and weakly evanescent spectra contributions, and (iii) development of an adaptive extension of an integral convergence acceleration technique to compute the strongly evanescent spectrum contribution. While other semianalytic techniques exist to solve this problem, none have full applicability to media exhibiting arbitrary double anisotropies in each layer, where one must account for the whole range of possible phenomena (e.g., mode coupling at interfaces and nonreciprocal mode propagation). Brute-force numerical methods can tackle this problem but only at a much higher computational cost. The present formulation provides an efficient and robust technique for field computation in arbitrary planar-stratified environments. We demonstrate the formulation for a number of problems related to geophysical exploration.

Structures in recurrence plots (RPs), preserving the rich information of nonlinear invariants and trajectory characteristics, have been increasingly analyzed in dynamic discrimination studies. The conventional analysis of RPs is mainly focused on quantifying the overall diagonal and vertical line structures through a method, called recurrence quantification analysis (RQA). This study extensively explores the information in RPs by quantifying local complex RP structures. To do this, an approach was developed to analyze the combination of three major RQA variables: determinism, laminarity, and recurrence rate (DLR) in a metawindow moving over a RP. It was then evaluated in two experiments discriminating (1) ideal nonlinear dynamic series emulated from the Lorenz system with different control parameters and (2) data sets of human heart rate regulations with normal sinus rhythms ($n$ = 18) and congestive heart failure ($n$ = 29). Finally, the DLR was compared with seven major RQA variables in terms of discriminatory power, measured by standardized mean difference ($D_{\mathrm{SMD}}$). In the two experiments, DLR resulted in the highest discriminatory power with $D_{\mathrm{SMD}}$ = 2.53 and 0.98, respectively, which were 7.41 and 2.09 times the best performance from RQA. The study also revealed that the optimal RP structures for the discriminations were neither typical diagonal structures nor vertical structures. These findings indicate that local complex RP structures contain some rich information unexploited by RQA. Therefore, future research to extensively analyze complex RP structures would potentially improve the effectiveness of the RP analysis in dynamic discrimination studies.

Using the ferroin- and the bathoferroin-catalyzed Belousov-Zhabotinsky (BZ) reaction embedded in the sodium-bis (2-ethylhexyl) sulfosuccinate (AOT) water-in-oil microemulsion, we observed different patterns occurring in two different solvents, hexane and octane. Turing patterns were found in both solvents with ferroin. They differ in their interaction with coexisting bulk oscillations, such that a new excitation front was formed around the evolving Turing patterns in hexane. However, in octane, the bulk oscillation merged with the evolving patterns, forming a new excitation front, which propagated into two directions: towards the existing patterns and away from them. For the bathoferroin-catalyzed BZ reaction, patterns like dash waves, jumping waves, and bubble waves were found in both solvents having different wavelengths. A curvature dependence of the splitting and merging of dashes was found.

Hysteresis cycles are very important features of energy conversion and harvesting devices, such as batteries. The efficiency of these may be strongly affected by the physical size of the system. Here, we show that in systems which are small enough, the existence of physical boundaries which produce nonhomogeneities of the interaction potential gives rise to inflections and barriers in the associated free energy. This in turn brings on irreversible processes which can be triggered under suitable external conditions imposed by a heat bath. As an example, by controlling the temperature, the state of a small system may be impelled to oscillate between two different structural configurations or aggregation states avoiding equilibrium coexistence and therefore dissipating energy. This cyclical behavior associated with a hysteresis cycle may be prototypical of energy conversion, storage, or generating nanodevices, as exemplified by Li-ion insertion batteries.