Researcher profile

Ruslan L. Davidchack

Ruslan L. Davidchack contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
11works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2020arXiv

A unified description of Surface Free Energy and Surface Stress

Even though the study of interfacial phenomena dates back to Laplace and was formalised by Gibbs, it appears that some concepts and relations among them are still causing some confusion and debates in the literature, particularly for interfaces involving solids. Moreover, ever since the Molecular Dynamics (MD) simulations have started to be widely used in the study of surface properties, these debates only intensified. In this work, we present a systematic description of the interfacial properties from the thermodynamic and statistical mechanics points of view. In particular, we link our derivations to MD simulations, describing precisely what different quantities represent and how they can be calculated. We do not follow the usual way that consists of describing the thermodynamics of the surfaces in general and then considering specific cases (e.g. liquid-liquid interface, liquid-solid interface). Instead, we present our analysis of various properties of surfaces in a hierarchical way, starting with the simplest case that we have identified: a single component liquid-vacuum interface, and then adding more and more complications when we progress to more complex interfaces involving solids. We propose that the term "surface tension" should not be used in the description of surfaces and interfaces involving solids, since its meaning is ambiguous. Only "Surface Free Energy" and "Surface Stress" are well defined and represent distinct, but related, properties of the interfaces. We demonstrate that these quantities, as defined in thermodynamics and measured in MD simulations, satisfy the Shuttleworth equation.

preprint2020arXiv

Pulsed Waveforms and Intermittently Nonlinear Filtering in Synthesis of Low-SNR and Covert Communications

In traditional spread-spectrum techniques, a wideband transmit signal is obtained by modulating a wideband carrier by a narrowband signal containing a relatively low-rate message. In the receiver, the respective demodulation/despreading restores the information-carrying narrowband signal. In this paper, we introduce an alternative approach, where the low-rate information is encoded directly into a wideband waveform of a given bandwidth, without physical "spreading" of the carrier's frequency. The main advantages of this approach lie in extended options for encoding the information, and in retaining a reversible control over the temporal and amplitude structures of the modulating wideband waveforms. Significant "excess bandwidth" (over that needed to carry the information) enables us to use allpass filters to manage statistical properties and time-domain appearances of these waveforms without changing their spectral composition. For example, a mixture of transmitted waveforms can be shaped as a low-crest-factor signal (e.g. to reduce the burden on the power amplifier), and/or made statistically indistinguishable from Gaussian noise (e.g. for covert transmissions and physical layer steganography), while the selected components of the received waveform can be transformed into high-crest-factor pulse trains suitable for multiplexing and/or low-SNR communications. Further, control over the temporal and amplitude structures of wideband waveforms carrying low-rate information enables effective use of nonlinear filtering techniques. Such techniques can be employed for robust real-time asynchronous extraction of the information, as well as for separation of wideband signal components with identical spectral content from each other. This can facilitate development of a large variety of low-SNR and covert communication configurations.

preprint2020arXiv

Utilizing Pileup Effect and Intermittently Nonlinear Filtering in Synthesis of Covert and Hard-to-Intercept Communication Links

We outline an approach to physical-layer steganography where the transmitted low-power stego messages are statistically indistinguishable from the Gaussian component of the channel noise (e.g. the thermal noise) observed in the same spectral band, and thus the channel noise itself serves as an effective cover signal. We also demonstrate how the apparent spectral and temporal properties of transmitted additional, higher-power cover signals (including those using the existing communication protocols) can be made to match those of the low-power stego payload and the Gaussian noise, providing extra layers of obfuscation for both the cover and the stego messages. We further illustrate how a specific combination of linear and nonlinear filtering can be used for effective separation of the cover, payload, and/or "friendly jamming" signals even when all transmissions have essentially the same spectral characteristics as well as temporal and amplitude structures, and when there are no explicit differences in the spectral and/or temporal allocations for the cover and the stego messages.

preprint2020arXiv

Utilizing Pulse Pileup Effect in Development of Robust Low-SNR Covert Communication Links

In contrast to other spread-spectrum techniques, wideband pulse trains with relatively low pulse arrival rates may be considered unsuitable for covert communications. The high crest factor of such trains can be extremely burdensome for the transmitter hardware, and it makes the pulse trains easily detectable even at very low signal-to-noise ratios. In addition, it may appear that sharing the wideband channel by multiple users would require explicit allocation of the pulse arrival times for each sub-channel, which would be impractical in most cases. On the other hand, messaging by wideband pulse trains has many appealing features. Among those are the ease of synchronous as well as asynchronous pulse detection, and on-the-fly channel reconfigurability (e.g. changing the spreading factor). Favorably, the crest factor of a pulse train, as well as its apparent temporal and amplitude structure, can be easily, and reversibly, controlled by simple linear filtering. For example, a transmitted pulse train can be made statistically indistinguishable from the Gaussian component of the channel noise (e.g. the thermal noise) observed in the same spectral band, while the received signal will be the designed high-crest-factor wideband pulse train. In this paper, we utilize the so-called pulse pileup effect to perform such reversible control of the pulse train structure, enabling a wider use of this approach for synthesis of robust low-SNR covert communication links. We place a particular focus on the synchronous pulse detection in the receiver, that provides a better utilization of the channel spectrum.

preprint2015arXiv

Reduction of SO(2) symmetry for spatially extended dynamical systems

Spatially extended systems, such as channel or pipe flows, are often equivariant under continuous symmetry transformations, with each state of the flow having an infinite number of equivalent solutions obtained from it by a translation or a rotation. This multitude of equivalent solutions tends to obscure the dynamics of turbulence. Here we describe the `first Fourier mode slice', a very simple, easy to implement reduction of SO(2) symmetry. While the method exhibits rapid variations in phase velocity whenever the magnitude of the first Fourier mode is nearly vanishing, these near singularities can be regularized by a time-scaling transformation. We show that after application of the method, hitherto unseen global structures, for example Kuramoto-Sivashinsky relative periodic orbits and unstable manifolds of travelling waves, are uncovered.

preprint2014arXiv

Discretization errors in molecular dynamics simulations with deterministic and stochastic thermostats

We investigate the influence of numerical discretization errors on computed averages in a molecular dynamics simulation of TIP4P liquid water at 300 K coupled to different deterministic (Nosé-Hoover and Nosé-Poincaré) and stochastic (Langevin) thermostats. We propose a couple of simple practical approaches to estimating such errors and taking them into account when computing the averages. We show that it is possible to obtain accurate measurements of various system quantities using step sizes of up to 70% of the stability threshold of the integrator, which for the system of TIP4P liquid water at 300 K corresponds to the step size of about 7 fs.

preprint2014arXiv

Parameterizing the surface free energy and excess adsorption of a hard-sphere fluid at a planar hard wall

he inhomogeneous structure of a fluid at a wall can be characterized in several ways. Within a thermodynamic description the surface free energy $γ$ and the excess adsorption $Γ$ are of central importance. For theoretical studies closed expression of $γ$ and $Γ$ can be very valuable; however, even for a well-studied model system such as a hard-sphere fluid at a planar hard wall, the accuracy of existing expressions for $γ$ and $Γ$, compared to precise computer simulation data, can still be improved. Here, we compare several known expressions for $γ$ and $Γ$ to the most precise computer simulation data. While good agreement is generally found at low to intermediate fluid densities, the existing parameterizations show significant deviation at high density. In this work, we propose new parameterizations for $γ$ and $Γ$ that agree with the simulation data within statistical error over the entire fluid density range.

preprint2012arXiv

Interfacial free energy of a hard-sphere fluid in contact with curved hard surfaces

Using molecular-dynamics simulation, we have calculated the interfacial free energy, γ, between a hard-sphere fluid and hard spherical and cylindrical colloidal particles, as functions of the particle radius R and the fluid packing fraction η= ρσ^3/6, where ρand σare the number density and hard-sphere diameter, respectively. These results verify that Hadwiger's theorem from integral geometry, which predicts that γfor a fluid at a surface, with certain restrictions, should be a linear combination of the average mean and Gaussian surface curvatures, is valid within the precision of the calculation for spherical and cylindrical surfaces up to ηabout 0.42. In addition, earlier results for γfor this system [Bryk, et al., Phys. Rev. E, 68, 031602 (2003)] using a geometrically-based classical Density Functional Theory are in excellent agreement with the current simulation results for packing fractions in the range where Hadwiger's theorem is valid. However, above ηabout 0.42, γ(R) shows significant deviations from the Hadwiger form indicating limitations to its use for high-density hard-sphere fluids. Using the results of this study together with Hadwiger's theorem allows one, in principle, to determine $γ$ for any sufficiently smooth surface immersed in a hard-sphere fluid.

preprint2009arXiv

Langevin Thermostat for Rigid Body Dynamics

We present a new method for isothermal rigid body simulations using the quaternion representation and Langevin dynamics. It can be combined with the traditional Langevin or gradient (Brownian) dynamics for the translational degrees of freedom to correctly sample the NVT distribution in a simulation of rigid molecules. We propose simple, quasi-symplectic second-order numerical integrators and test their performance on the TIP4P model of water. We also investigate the optimal choice of thermostat parameters.

preprint2009arXiv

On the state space geometry of the Kuramoto-Sivashinsky flow in a periodic domain

The continuous and discrete symmetries of the Kuramoto-Sivashinsky system restricted to a spatially periodic domain play a prominent role in shaping the invariant sets of its chaotic dynamics. The continuous spatial translation symmetry leads to relative equilibrium (traveling wave) and relative periodic orbit (modulated traveling wave) solutions. The discrete symmetries lead to existence of equilibrium and periodic orbit solutions, induce decomposition of state space into invariant subspaces, and enforce certain structurally stable heteroclinic connections between equilibria. We show, on the example of a particular small-cell Kuramoto-Sivashinsky system, how the geometry of its dynamical state space is organized by a rigid `cage' built by heteroclinic connections between equilibria, and demonstrate the preponderance of unstable relative periodic orbits and their likely role as the skeleton underpinning spatiotemporal turbulence in systems with continuous symmetries. We also offer novel visualizations of the high-dimensional Kuramoto-Sivashinsky state space flow through projections onto low-dimensional, PDE representation independent, dynamically invariant intrinsic coordinate frames, as well as in terms of the physical, symmetry invariant energy transfer rates.

preprint2008arXiv

On the use of stabilising transformations for detecting unstable periodic orbits in the Kuramoto-Sivashinsky equation

In this paper we develop further a method for detecting unstable periodic orbits (UPOs) by stabilising transformations, where the strategy is to transform the system of interest in such a way that the orbits become stable. The main difficulty of using this method is that the number of transformations, which were used in the past, becomes overwhelming as we move to higher dimensions (Davidchack and Lai 1999; Schmelcher et al. 1997, 1998). We have recently proposed a set of stabilising transformations which is constructed from a small set of already found UPOs (Crofts and Davidchack 2006). The main benefit of using the proposed set is that its cardinality depends on the dimension of the unstable manifold at the UPO rather than the dimension of the system. In a typical situation the dimension of the unstable manifold is much smaller than the dimension of the system so the number of transformations is much smaller. Here we extend this approach to high-dimensional systems of ODEs and apply it to the model example of a chaotic spatially extended system -- the Kuramoto-Sivashinsky equation. A comparison is made between the performance of this new method against the competing methods of Newton-Armijo (NA) and Levernberg-Marquardt (LM). In the latter case, we take advantage of the fact that the LM algorithm is able to solve under-determined systems of equations, thus eliminating the need for any additional constraints.