This PDF file contains the front matter associated with SPIE Proceedings Volume 10717, including the Title Page, Copyright information, Table of Contents, and Conference Committee listing.

The paper justifies 𝐸²-similarity of kinetic equation solutions to describe vacuum emergence of electronpositron plasma under the effect of strong “laser” fields, where 𝐸(𝑡) is the intensity of the strong time-dependent “laser” field. The boundaries of existence of this similarity were studied.

The problem of determining the thickness of thin films by the method of laser triangulation is considered. An expression is derived for the film thickness and the distance between the focused beams on the photo detector. The possibility of applying the chosen method for measuring thickness is in the range [0.1; 1] mm. We could resolve 2 individual light marks for a minimum film thickness of 0.23 mm. We resolved with the help of computer processing of photos with a resolution of 0.10 mm. The obtained results can be used in ophthalmology for express diagnostics during surgical operations on the corneal layer.

There are several trace gases such as N₂O, CO, CO₂, NO, H₂O, NO₂, NH₃, CH₄ etc. which have their absorption peaks in Mid-IR spectrum These gases strongly absorb in the mid-IR > 2.5 μm spectral region due to their fundamental rotational and vibrational transitions. In this work, we modelled and optimized three different kinds of waveguides such as rib, strip and slot based on silicon platform to obtain maximum evanescent field ratio. These waveguides are designed at 3.39 μm and 4.67 μm which correspond to the absorption line of methane (CH₄) and carbon monoxide (CO) respectively.

Kinetics of vibrationally-excited singlet oxygen O₂(a¹Δ,ν) in gas mixture O₃/N₂/CO₂ was studied using a pulse laser technique. Molecules O₂(a¹Δ,ν) were produced by laser photolysis of ozone at 266 nm. The O₃ molecules number density was followed using time-resolved absorption spectroscopy. It was found that an upper bound for the rate constant of chemical reaction O₂(a¹Δ,ν)+ O₃ is about 10^-15 cm³/s. The rate constants of O²(a¹Δ,ν= 1, 2 and 3) quenching by CO₂ are presented.

Rate constants for the removal of O₂ b¹∑⁺_g by collisions with O₂, N₂, CO₂ and H₂O have been determined at temperature 297 K. O₂(b¹ ∑⁺_g) was excited by pulses from a tunable dye laser, and the deactivation kinetics were followed by observing the temporal behavior of the b¹∑⁺_g - X³∑^-_g fluorescence. The removal rate constants for CO₂, N₂ and H₂O were not strongly dependent on temperature, and could be represented by the expressions k_CO2=(1.8±0.05)×10^-16; k_N2=(2.2 ± 0.2)×10^-15, and k_H2O=(6.12±0.67)×10^-12 cm³ molecule^-1 s^-1. Rate constant for O₂(b¹∑⁺ ) removal by O₂(X), being orders of magnitude lower, represented by the fitted expression k_O2=(3.67 ± 0.06)×10^-17 cm³ molecule^-1 s^-1. All of the rate constants measured at room temperature were found to be in good agreement with previously reported values.

Currently, lasers are widely used for surgery, medical diagnostics and oncology research. Unfortunately, most of the used laser sources have a significant drawback - the lack of operating wavelength tuning possibility, which imposes significant limitations on the investigation of biological tissues spectral properties and searching for the optimal mode of their treatment. Comparison between different promising mid-IR sources was made. We report on development of mid-infrared (mid-IR) tunable lasers based on the Cr2+:CdSe single-crystals. These lasers operate in CW mode with the maximum output power of up to 2 W and possible tuning range from 2.2 to 3.6 μm.

The possibilities of self-mixing interferometry for measuring nanodisplacement of a probe included in a near-field scanning microwave microscope have been considered. The features of the formation of a laser interference signal at current modulation of the wavelength of laser radiation have been investigated. Experimental responses of a semiconductor laser system included in scanning microwave microscope to control nanodisplacement of the probe have been demonstrated.To register the nanodisplacement of the probe, it is proposed to use the method of determining the stationary phase of a laser interference signal by low-frequency spectrum of a semiconductor laser. The change of the amplitudes of the spectral components in the spectrum of the interference signal due to creation of the standing wave in the external resonator of the laser self-mixing system has been shown. The form of the interference signal at current modulation of the radiation wavelength was experimentally obtained when the probe moves with a step of 80 nm. The results of measuring nanodisplacements of an electromagnetic translator STANDA 8MVT40-13 have been demonstrated. Deviation of the nanodisplacement of the proposed method does not exceed 15%.

The results of studies on the effect of temperature on the output spectral characteristics of continuous semiconductor lasers of the visible range are presented. The paper presents the results of studying the spectral-optical radiation parameters of semiconductor lasers, their coherence lengths, and the dependence of the position of the spectral peak of the wavelength on temperature. This is necessary for the selection of the most optimal laser in order to use it for medical ophthalmologic diagnosis. The experiment was carried out using semiconductor laser modules based on a laser diode. The spectra were recorded by using a two-channel automated spectral complex based on the MDR-23 monochromator. Spectral dependences on the temperature of semiconductor lasers are obtained, in the range from 300 to 370 K. The possibility of determining the internal damage to the stabilization of laser modules without opening the case is shown, but only with the use of their spectral characteristics. The obtained data allow taking into account temperature characteristics and further optimization of parameters of such lasers when used in medical practice, in particular, in ophthalmologic diagnostics.

Medium-power (10-100 W) THz-band gyrotrons operating in a continuous-wave (CW) mode are of great importance for many applications such as NMR spectroscopy with dynamic nuclear polarization (DNP/NMR), plasma diagnostics, nondestructive inspection, stand-off detection of radioactive materials, biomedical applications, etc. For all these applications, high frequency stability and tunability within 1-2 GHz frequency range is typically required. Apart from different existing techniques for frequency stabilization, phase locking has recently attracted strong interest. In this paper, we present the results of theoretical analysis and numerical simulation for several phase locking techniques: (a) phase locking by injection of the external driving signal; (b) mutual phase locking of two coupled gyrotrons; and (c) selfinjection locking by a wave reflected from the remote load.

We present a study of the features of generalized synchronization in laser models involving coupled class B lasers with delay, as well as the singularities of synchronization of large ensembles, arising from the antiphase synchronization of individual pairs of lasers.

This work is devoted to a theoretical investigation of the Ramsey method of detection of the coherent population trapping resonance in cold atomic clouds taking into account collective effects caused by finite optical depth of the considered clouds. The interaction of atoms with pulsed laser radiation is described in the formalism of density matrix by means of Maxwell-Bloch set of equations. The Ramsey signal of coherent population trapping resonance was calculated for the radiation passed through the medium and analyzed for different length of the atomic cloud. Also the population of excited level was calculated in dependence on the two-photon detuning and coordinate along the main optical axis. The light shift of sidebands and appearance of additional harmonics were discovered.

The silicate materials are perspective for different areas of laser physics and photonics. In this paper the comparison of the nonlinear conversion with the generation of the light harmonic in Al- and Ge-containing silicate materials is presented. The peculiarities of the processes of the light harmonic generation in dependence on the concentrations of the chemical components are discussed and the influences of the additional small inclusion of the elements of fifth group and the rare-earth elements are estimated.

Investigations of the multi-peaks scattering of the laser light on the micro-scale susceptibility gratings with small periodicities photo-induced in the various glass materials are presented. The observed pictures of the multi-peaks scattering of light in oxide samples show that the efficiencies of the processes of scattering can vary for the different chemical compositions. Experimental results are in agreement with the proposed theory of light scattering.

A quantum model of spin dynamics of spin-orbit coupled two-dimensional electron gas in the presence of strong high- frequency electromagnetic field is suggested. Interaction of electrons with optical phonons is taken into account in the second order of perturbation theory.

The theoretical possibility of retrieving the additional information on the dielectric properties of the nanoparticles material by single scattering in suspensions was studied. We have demonstrated a method of recreating the dielectric function of the material in the fundamental absorption band using the closed aperture z-scanning with the simultaneous Rayleigh scattering intensity measurements and the polarization control of an input laser beam. A possibility to recreate the form factor of the non-spherical particles or anisotropic nonlinear sensitivity for the sphere-like particles was also observed.

The manuscript addresses the advantages and possible applications of acousto-optic image spectral filtration in lowcoherence interferometry. In particular, an effective operation of acousto-optical tunable filters in combination with Michelson-type interferometers is shown. The results of original experiments are presented. It is demonstrated that amplitude and phase spatial distributions of light waves reflected from or transmitted through the object can be fast determined in contactless manner for any spectral intervals with use of the presented techniques.

The mandatory condition for efficient operation of an optically-pumped all-rare-gas laser (OPRGL) is the presence of rare gas metastable atoms in the discharge plasma with number density of the order of 10¹²-10¹³ cm^-3. This requirement mainly depends on the choice of a discharge system. In this study the number density values of argon metastable atoms were obtained in the condition of the dielectric-barrier discharge (DBD) at an atmospheric pressure.

The theoretical investigation results of disintegration effect of elliptic polarized shot probe pulses of electromagnetically induced transparency in the counterintuitive superposed elliptic polarized control field and in weak probe field approximation are presented. It is shown that this disintegration occurs because the probe field in the medium is the sum of two normal modes, which correspond to elliptic polarized pulses with different speeds of propagation. The polarization ellipses of normal modes have equal eccentricities and mutually perpendicular major axes. Major axis of polarization ellipse of one normal mode is parallel to polarization ellipse major axis of control field, and electric vector of this mode rotates in the opposite direction, than electric vector of the control field. The electric vector other normal mode rotates in the same direction that the control field electric vector. The normal mode speed of the first type aforementioned is less than that of the second type. The polarization characteristics of the normal mode depend uniquely on the polarization characteristics of elliptic polarized control field and remain changeless in the propagation process. The theoretical investigation is performed for Λ-scheme of degenerated quantum transitions between ³P₀, ³P₁⁰ and ³P₂ energy levels of ²⁰⁸Pb isotope.

Using the cabaret method numerical realization, long-time spatio-temporal dynamics of the electromagnetic field in a nonlinear ring fibre cavity with dispersion is investigated during the hundreds of round trips. Formation of both the temporal cavity solitons and irregular pulse trains is demonstrated and discussed.

It is well known that polarisation transformations can occur not only under sharp focusing conditions, but also in the paraxial regime. We investigate such spatial polarisation transformations of three-dimensional distribution in the focal region in the paraxial regime due to the addition of a lens by a spiral binary diffraction axicon. Theoretical analysis demonstrates the possibility of such transformations. We demonstrate that, in these cases, a zero intensity is formed on the beam axis before the focal plane, and an intensity peak is formed on the beam axis after the passage of the focal plane. The experimental results obtained are in good agreement with the simulation results. We believe that results can be useful in the field of optical manipulation of nano- and microscale objects, as well as in the field of material processing.

In this paper, the focusing of laser light at 532 nm by dielectric cylinders with a metal shells is studied by using COMSOL Multiphysics. The analysis of cylinder design which proposed multilayered shell shows that a microcylinder with a gold-silver (or silver-gold) shell can improve the focusing process, especially in the case of TM polarization. The microcylinder with thin internal silver layer of 1 nm and outside gold layer of 9 nm focus TE-polarized light to nanojet with maximal intensity of 5.65 a.u., full width and full length at half maximum of intensity of of 0.39λ and 0.72λ, respectively.

In this work, we fabricated and studied the performance of a 100×100-μm four-sector binary subwavelength grating polarizer in a golden film. It was experimentally demonstrated that a linearly polarized 532-nm Gaussian beam reflected at the polarizer was converted to an azimuthally polarized beam. Putting a spiral phase plate (SPP) with the topological charge n = 1 into the azimuthally polarized beam from the micropolarizer was experimentally shown to enable the conversion of the annular intensity pattern into a central intensity peak.

Artificial periodic structures offers an additional degree of freedom in manipulation of nonlinear waves and solitons. It is shown that periodical variation of the coefficients of nonlinear Schrödinger equation allows to change soliton eigenvalues. Pairwise interaction of solitons is demonstrated. The interacting solitons can be selected by appropriate period of the variation of the coefficients of nonlinear Schrödinger equation. The proposed approach for controlling soliton eigenvalues allows to manipulate soliton fission process.

The propagation of nanosecond laser pulses in solutions of large-scale carbon nanostructures has been theoretically and numerically investigated in this paper. The three-level ladder-type scheme was used to study the medium with RSA. Numerical simulation of the propagation of nanosecond laser pulses in RSA media was carried out using the classical transport equation for incoherent nonlinear interaction of radiation with the medium described by the proposed three-level model. It is shown that for sufficiently small relaxation times of the medium and high intensity of, a splitting of the transmitted pulse by two is possible.

We demonstrate a new type of circular Airy beams, the so-called azimuthally modulated circular Airy beams, generated by utilizing a diffraction element, whose transmission function is the sum of the transmission function of the element generating a “petal” pattern and the transmission function of the element generating a circular Airy beam. We experimentally investigate the propagation dynamics of such beams and demonstrate that their autofocusing and selfhealing properties are strongly dependent on the number of generated petals. These beams are a combination of a conventional circular Airy beam and vortex laser beams (or their superpositions). Using a spatial light modulator, we demonstrate that these beams have unique properties such as autofocusing, "nondiffractive" propagation and self-healing after passing through an obstacle. The experimental results are in good agreement with the simulation. We believe that these results can be very useful for lensless laser fabrication and laser manipulation techniques, as well as for development of new filament plasma multi-channel formation methods.

The properties of quantum networks based on generalized Tavis-Cummings models are theoretically investigated. We have calculated the information transfer success rate from one node to another in a simple model of a quantum network realized with two-level atoms placed in the cavities and interacting with an external laser field and cavity photons. The method of dynamical group of the Hamiltonian and technique of corresponding coherent states were used for investigation of the temporal dynamics of the two nodes model.

The probability of quantum transitions of a molecule between its states under the action of an electromagnetic field is represented as an integral over trajectories from a real alternating functional. A method is proposed for computing the integral using recurrence relations. The method is attached to describe the two-photon Rabi oscillations.

An exact solution of the problem of two-atom one- and two-mode Jaynes-Cummings model with intensity- dependent coupling is presented. Asymptotic solutions for system state vectors are obtained in the approximation of large initial coherent fields. The atom-field entanglement is investigated on the basis of the reduced atomic entropy dynamics. The possibility of the system being initially in a pure disentangled state to revive into this state during the evolution process for both models is shown. Conditions and times of disentanglement are derived.

We investigated the entanglement between two qubits interacting with two independent modes of lossless resonators taking into account the dipole-dipole interaction. The model with different qubit-field couplings and detunings is under consideration. We derived the exact solution for considered model and investigated the dependence of the atom-atom entanglement on the strength of the dipole-dipole interaction and other parameters of the considered system. The results showed that the presence of a sufficiently large dipole-dipole interaction leads to stabilization of initial atom-atom entanglement.

Two algorithms are presented that use the properties of twisted light to construct symmetric encryption schemes. Their possible physical realization is discussed.

On the basis of quantum microscopic approach we study the excitation dynamics of two motionless point atoms located near the perfectly conducting mirror. We have analyzed the spontaneous decay rate of individual atoms near the mirror as well as the strength of dipole-dipole interaction between different atoms. It is shown that the spontaneous decay rate of an excited atom significantly depends on the distance from this atom to the mirror. In the case when the interatomic separation is less or comparable with the wavelength of resonant radiation, the spontaneous decay dynamics of an excited atom is described by multi-exponential law. It depends both the interatomic separation and the spatial orientation of diatomic quasimolecule.

We have investigated the dynamics of entanglement of two identical qubits non-resonantly interacting with one mode of a thermal electromagnetic field in a lossless cavity in the presence of a direct dipole-dipole interaction. On the basis of the exact solution the evolution equation the negativity for qubits has been found. Numerical simulation of the negativity for various model parameters has been carried out. It has been shown that the interaction of qubits with the cavity thermal field can lead to their entanglement. It has been established that the detuning and dipole-dipole interaction of qubits can be used to manipulate and control the degree of their entanglement.

Nowadays the geometric approach in optics is often used to find out media parameters based on propagation paths of the rays because in this case it is a direct problem. However inverse problem in the framework of geometrized optics is usually not given attention. The aim of this work is to demonstrate the work of the proposed the algorithm in the framework of geometrized approach to optics for solving the problem of finding the propagation path of the electromagnetic radiation depending on environmental parameters. The methods of differential geometry are used for effective metrics construction for isotropic and anisotropic media. For effective metric space ray trajectories are obtained in the form of geodesic curves. The introduced algorithm is applied to well-known objects, Maxwell and Luneburg lenses. The similarity of results obtained by classical and geometric approach is demonstrated.

The article is devoted to the calculation of normal modes of the closed waveguides with an arbitrary filling ε, μ in the system of computer algebra Sage. Maxwell equations in the cylinder are reduced to the system of two bounded Helmholtz equations, the notion of weak solution of this system is given and then this system is investigated as a system of ordinary differential equations. The normal modes of this system are an eigenvectors of a matrix pencil. We suggest to calculate the matrix elements approximately and to truncate the matrix by usual way but further to solve the truncated eigenvalue problem exactly in the field of algebraic numbers. This approach allows to keep the symmetry of the initial problem and in particular the multiplicity of the eigenvalues. In the work would be presented some results of calculations.

The exactly soluble model of a train of zero-duration electromagnetic pulses interacting with a 1D atom with short-range interaction potential modelled by a δ-function is considered. The model is related to the up-to-date laser techniques providing the duration of pulses as short as a few attoseconds and the intensities higher than 10¹⁴ W/cm².

In this paper we investigate the solution of the problem of modeling the propagation of electromagnetic radiation in three-dimensional integrated optical structures, such as waveguide lenses. When propagating through three-dimensional waveguide structures the waveguide modes can be hybridized, so the mathematical model of their propagation must take into account the connection of TE- and TM-mode components. Therefore, an adequate consideration of hybridization of the waveguide modes is possible only in vector formulation of the problem. An example of three-dimensional structure that hybridizes waveguide modes is the Luneburg waveguide lens, which also has focusing properties.

If the waveguide lens has a radius of the order of several tens of wavelengths, its variable thickness at distances of the order of several wavelengths is almost constant. Assuming in this case that the electromagnetic field also varies slowly in the direction perpendicular to the direction of propagation, one can introduce a small parameter characterizing this slow varying and decompose the solution in powers of the small parameter. In this approach, in the zeroth approximation, scalar diffraction problems are obtained, the solution of which is less resource-consuming than the solution of vector problems. The calculated first-order corrections of smallness describe the connection of TE- and TM-modes, so the solutions obtained are weakly-hybridized modes.

The formulation of problems and methods for their numerical solution in this paper are based on the authors' research on waveguide diffraction on a lens in a scalar formulation.

We present new calculation schemes using high-order finite element method implemented on unstructured grids with triangle elements for solving boundary-value problems that describe axially symmetric quantum dots. The efficiency of the algorithms and software is demonstrated by benchmark calculations of the energy spectrum, the envelope eigenfunctions of electron, hole and exciton states, and the direct interband light absorption in conical and spheroidal impenetrable quantum dots.

The transformation optics uses geometrized Maxwell’s constitutive equations to solve the inverse problem of optics, namely to solve the problem of finding the parameters of the medium along the paths of propagation of the electromagnetic field. For the geometrization of Maxwell’s constitutive equations, the quadratic Riemannian geometry is usually used. This is due to the use of the approaches of the general relativity. However, there arises the question of the insufficiency of the Riemannian structure for describing the constitutive tensor of the Maxwell’s equations. The authors analyze the structure of the constitutive tensor and correlate it with the structure of the metric tensor of Riemannian geometry. It is concluded that the use of the quadratic metric for the geometrization of Maxwell’s equations is insufficient, since the number of components of the metric tensor is less than the number of components of the constitutive tensor. A possible solution to this problem may be a transition to Finslerian geometry, in particular, the use of the Berwald-Moor metric to establish the structural correspondence between the field tensors of the electromagnetic field.

The electric and magnetic surface plasmons in symmetric structures of metallic and dielectric layers are considered. The existence of backward and forward waves and the slow and fast plasmon-polaritons are obtained. It is shown that the anomalous negative dispersion in the structures with dissipation does not necessarily indicate the backward surface plasmons.

We consider the structures type of two-dimensional electron gas in the form of a thin conductive, in particular, graphene films described by tensor conductivity, which are isolated or located on the dielectric layers. The dispersion equation for hybrid modes, as well as scattering parameters. We show that free wave (eigenwaves) problem follow from the problem of diffraction when linking the amplitude of the current of the linear equations are unsolvable, i.e., the determinant of this system is zero. As a particular case the dispersion equation follow from the conditions of matching (with zero reflection coefficient).

In the paper, an oblique incidence of a plane electromagnetic wave on a planar magnetically anisotropic inhomogeneous layer is considered. We consider the case when all the components of the magnetic permeability tensor are non zero and vary with distance from the interface of media. The WKB method gives a matrix 4 × 4 solution for the projections of the electromagnetic wave fields during its propagation. The dependence of the cross-polarized components on the orientation of the anisotropic medium relative to the plane of incidence of the medium is analyzed.

The reflection of linearly polarized light at the boundary of a plane anisotropic layer is considered, the optical axis of which is located arbitrarily in relation to the plane of incidence. Matrix amplitude coefficients of reflection and transmission of light by a plane anisotropic uniaxial layer are obtained, the energy coefficients of reflection, as well as ellipsometric parameters of the reflected light, and also their angular spectra.

Electromagnetic waves propagation in the complex cavity with anisotropic hyperbolic metamaterial are investigated using direct calculation of modal field and dispersion equation. The transfer matrix method was adopted for arbitrary orientation of optical axis according to slab boundary. Increasing of the density of states in the cavity have show.

The profound research and physical applications of interactions of different types of waves with medium are very important. Particularly the most interesting sphere is for complex environments, which may be characterized by the increasing number of methods. Their objective analysis increased because of great applied significance. For the optical range it comes to considering the structure, the dimensions of the spatial inhomogeneity of which are comparable to the wavelength of the radiation.

Microfabricated vacuum-tube millimeter- and THz-band sources are of great interest for numerous applications such as communications, radar, sensors, imaging, etc. Recently, miniaturized sheet-beam traveling-wave tubes for sub-THz and THz operation have attracted a considerable interest. In this paper, we present the results of modeling and development of slow-wave structures (SWS) for medium power (10-100 W) traveling-wave tube (TWT) amplifiers and backwardwave oscillators (BWO) in near-THz frequency band. Different types of SWSs are considered, such as double-vane SWS for TWT with a sheet electron beam, a folded-waveguide SWS, and novel planar SWSs on dielectric substrates.

This paper investigates nontrivial temporal dynamics of the optical field in the nonlinear ring resonator (containing the metamaterial with the cubic nonlinearity of the Kerr type) driven by an external optical injection and subjected to a delayed feedback. Unstable periodic orbits are obtained for definite sets of parameters including the delay time, feedback strength, feedback phase, amplitude of injection field and detuning. The bifurcation diagram of nonlinear behavior of this optical system is constructed.

Results of the experimental study of changes in the bubble size statistics during aging of wet foams are discussed. It is proposed that the evolution of the bubble radii distributions can be described in terms of the one dimensional Fokker- Planck equation. The empirical distributions of the bubble radii exhibit a self-similarity of their shapes and can be transformed to a time-independent form using the radius renormalization. Analysis of obtained data allows us to suggest that the drift term of the Fokker-Planck equation dominates in comparison with the diffusion term in the case of aging of isolated quasi-stable wet foams.

Empirical data on the diffuse and collimated transmittance of aged liquid foams are discussed in terms of influence of mutual correlations in the scatter positions. This influence can be described introducing the static structure factor of a scattering system and occurs remarkable in the case of wet foams with gas bubbles as the basic scattering units. On the contrary, mutual correlations of basic scattering units (Plateau-Gibbs channels and vertices) in dry foams are negligible due to low values of their volume fraction. This causes dramatic changes of the scattering anisotropy of foam layers in the vicinity of the wet-to-dry transition. Some analogies can be drawn between this effect and a previously reported "optical inversion" of densely packed random media.

Dielectric function of wide-zone semiconductor nanoparticles (titanium dioxide) was studied under the condition of laser pumping at various wavelengths. A closed-aperture z-scan method with simultaneous measurements of the right-anglescattered intensity was used to retrieve the real and imaginary parts of dielectric function in the dependence on the pump intensity. It was found that the efficiency of dielectric function modulation by pumping light strongly depends on detuning of the wavelength of pumping light with respect to the fundamental absorption band of nanoparticles. The ColeCole diagrammatic technique was applied for interpretation of the pump-induced changes of the dielectric function in the optical range. Applicability of the Kramers-Kronig relations for description of the observed behavior of the dielectric function is discussed.

Features of diffusing light probes of aged liquid foams are discussed. These probes were carried out using broadband and laser radiation with the wavelength of 532 nm. Experimental data were obtained for aged samples of Gillette shaving cream as a model foamed substance. The coefficients of diffusion and collimated transmittance were applied for characterization of the structural changes in the aged wet foams. Also, the changes in the liquid fraction due to gravitational drainage were monitored using volumometric measurements in the course of the foam aging. Obtained empirical data on the diffuse and collimated transmittance were used as the reference values for correction of the spectral measurements in the visible range with a broadband source of probe light. The problem of correction of the collimated transmittance partially corrupted by the diffusing component of multiply scattered light is discussed.

We consider the practical realization of a new optical probe method of the random media which is defined as the reference-free path length interferometry with the intensity moments analysis. A peculiarity in the statistics of the spectrally selected fluorescence radiation in laser-pumped dye-doped random medium is discussed. Previously established correlations between the second- and the third-order moments of the intensity fluctuations in the random interference patterns, the coherence function of the probe radiation, and the path difference probability density for the interfering partial waves in the medium are confirmed. The correlations were verified using the statistical analysis of the spectrally selected fluorescence radiation emitted by a laser-pumped dye-doped random medium. Water solution of Rhodamine 6G was applied as the doping fluorescent agent for the ensembles of the densely packed silica grains, which were pumped by the 532 nm radiation of a solid state laser. The spectrum of the mean path length for a random medium was reconstructed.

Sol-gel transition was studied using the speckle correlometry method with a localized light source and spatial filtering of backscattered radiation. Water solutions of technical or food gelatin with added TiO₂ nanoparticles were used as studied objects. Structural transformation of "sol-gel" system was studied at various temperatures from 25°C to 50°C using analysis of the correlation and structure functions of speckle intensity fluctuations. The characteristic temperatures of “sol – gel” transition were evaluated for studied systems. Obtained results can be used for various applications in biomedicine and food industry.

The temporal evolution of the metastable and unstable foams had been studied. Diffusion wave spectroscopy was chosen as the diagnostic method, with calculation of the correlation time of the fluctuations in the intensity of the probing radiation. It was established that the correlation time increases with the time according to the power law with different parameters, depending on the type of the evolution and was found to be equal to 0.5 for the case of the metastable and to 2,52 for the unstable foam. It was also determined that the behaviour of the correlation time agrees well with the evolution of the characteristic dimensions of the scatterers in the form of bubbles in the medium, which can be used for contactless monitoring of the foaming processes in the production of the foam-like materials for various applications, for example, in the synthesis of the biocompatible polymer matrices - scaffolds.

In this paper, based on the apparatus of artificial neural networks, a technique for recognizing and classifying patterns corresponding to imaginary movements on electroencephalograms (EEGs) obtained from a group of untrained subjects was developed. The works on the selection of the optimal type, topology, training algorithms and neural network parameters were carried out from the point of view of the most accurate and fast recognition and classification of patterns on multi-channel EEGs associated with the imagination of movements. The influence of the number and choice of the analyzed channels of a multichannel EEG on the quality of recognition of imaginary movements was also studied, and optimal configurations of electrode arrangements were obtained. The effect of pre-processing of EEG signals is analyzed from the point of view of improving the accuracy of recognition of imaginary movements.

In this paper we found a correlation between the characteristics of a person revealed in classical psychological testing on the basis of Schulte tables, and its neurophysiological features of the functioning of the brain obtained from the time-frequency analysis of EEG. The results obtained are interesting from the point of view of the choice of training strategies for a particular individual. We believe that the obtained results are of interest for fundamental science and applied works of psychological testing and diagnostics. The study of such forming strategies on EEG data can be automated and do not require the work of highly skilled psychologists.

We have considered time-frequency and spatio-temporal structure of electrical brain activity, associated with real and imaginary movements based on the multichannel EEG recordings. We have found that along with wellknown effects of event-related desynchronization (ERD) in α/μ – rhythms and β – rhythm, these types of activity are accompanied by the either ERS (for real movement) or ERD (for imaginary movement) in low-frequency δ – band, located mostly in frontal lobe. This may be caused by the associated processes of decision making, which take place when subject is deciding either perform the movement or imagine it. Obtained features have been found in untrained subject which it its turn gives the possibility to use our results in the development of brain-computer interfaces for controlling anthropomorphic robotic arm.

In the paper we propose an approach based on artificial neural networks for recognition of different human brain states associated with distinct visual stimulus. Based on the developed numerical technique and the analysis of obtained experimental multichannel EEG data, we optimize the spatiotemporal representation of multichannel EEG to provide close to 97% accuracy in recognition of the EEG brain states during visual perception. Different interpretations of an ambiguous image produce different oscillatory patterns in the human EEG with similar features for every interpretation. Since these features are inherent to all subjects, a single artificial network can classify with high quality the associated brain states of other subjects.

This paper considers the possibility of classification of electroencephalogram (EEG) and electromyogram (EMG) signals corresponding to different phases of sleep and wakefulness of mice by the means of artificial neural networks. A feed-forward artificial neural network based on multilayer perceptron was created and trained on the data of one of the rodents. The trained network was used to read and classify the EEG and EMG data corresponding to different phases of sleep and wakefulness of the same mouse and other mouse. The results show a good recognition quality of all phases for the rodent on which the training was conducted (80–99%) and acceptable recognition quality for the data collected from the same mouse after a stroke.

In the present paper, the possibility of classification by artificial neural networks of a certain architecture of ambiguous images is investigated using the example of the Necker cube from the experimentally obtained EEG recording data of several operators. The possibilities of artificial neural network classification of ambiguous images are investigated in the different frequency ranges of EEG recording signals.

We study abilities of the wavelet-based multifractal analysis in recognition specific dynamics of electrical brain activity associated with real and imaginary movements. Based on the singularity spectra we analyze electroencephalograms (EEGs) acquired in untrained humans (operators) during imagination of hands movements, and show a possibility to distinguish between the related EEG patterns and the recordings performed during real movements or the background electrical brain activity. We discuss how such recognition depends on the selected brain region.

Characterization of dynamical features of chaotic oscillations from point processes is based on embedding theorems for non-uniformly sampled signals such as the sequences of interspike intervals (ISIs). This theoretical background confirms the ability of attractor reconstruction from ISIs generated by chaotically driven neuron models. The quality of such reconstruction depends on the available length of the analyzed dataset. We discuss how data resampling improves the reconstruction for short amount of data and show that this effect is observed for different types of mechanisms for spike generation.

We have proposed brain-computer interface (BCI) for the estimation of the brain response on the presented visual tasks. Proposed BCI is based on the EEG recorder Encephalan-EEGR-19/26 (Medicom MTD, Russia) supplemented by a special home-made developed acquisition software. BCI is tested during experimental session while subject is perceiving the bistable visual stimuli and classifying them according to the interpretation. We have subjected the participant to the different external conditions and observed the significant decrease in the response, associated with the perceiving the bistable visual stimuli, during the presence of distraction. Based on the obtained results we have proposed possibility to use of BCI for estimation of the human alertness during solving the tasks required substantial visual attention.

The most abundant non-neuronal cells in the brain, astrocytes, populate all parts of the central nervous system (CNS). Astrocytic calcium activity ranging from subcellular sparkles to intercellular waves is believed to be the key to a plethora of regulatory pathways in the central nervous system from synaptic plasticity to blood flow regulation. Modeling of the calcium wave initiation and transmission and their spatiotemporal dynamics is therefore an important step stone in understanding the crucial cogs of cognition. Astrocytes are active sensors of ongoing neuronal and synaptic activity, and neurotransmitters diffusing from the synaptic cleft make a strong impact on the astrocytic activity. Here we propose a model describing the patterns of calcium wave formation at a single cell level and discuss the interplay between astrocyte shape the calcium waves dynamics driven by local stochastic surges of glutamate simulating synaptic activity.

In this paper we study specific oscillatory patterns on EEG signals of WAG/Rij rats. These patterns are known as proepileptic because they occur in time period during the development of absence-epilepsy before fully-formed epileptic seizures. In the paper we analyze EEG signals of WAG/Rij rats with continuous wavelet transform and empirical mode decomposition in order to find particular features of epileptic spike-wave discharges and nonepileptic sleep spindles. Then we introduce proepileptic activity as patterns that combine features of epileptic and non-epileptic activity. We analyze proepileptic activity in order to specify its features and time-frequency structure.

The purpose of the paper is to propose a fully automated method of volume assessment of structures within human brain. Our statistical approach uses maximum interdependency principle for decision making process of measurements consistency and unequal observations. Detecting outliers performed using maximum normalized residual test. We propose a statistical model which utilizes knowledge of tissues distribution in human brain and applies partial data restoration for precision improvement. The approach proposes completed computationally efficient and independent from segmentation algorithm used in the application.

Here we studied the role of cerebral lymphatic system in the brain clearing using intraparenchymal injection of Evans Blue and gold nanorods assessed by optical coherent tomography and fluorescence microscopy. Our data clearly show that the cerebral lymphatic system plays an important role in the brain cleaning via meningeal lymphatic vessels but not cerebral veins. Meningeal lymphatic vessels transport fluid from the brain into the deep cervical node, which is the first anatomical “station” for lymph outflow from the brain. The lymphatic processes underlying brain clearing are more slowly vs. peripheral lymphatics. These results shed light on the lymphatic mechanisms responsible for brain clearing as well as interaction between the intra- and extracranial lymphatic compartment.

Transparenchymal transport attracted the attention of many research groups after the discovery of glymphatic mechanism for the brain drainage in 2012. While the main facts of rapid transport of substances across the parenchyma are well established experimentally, specific mechanisms that drive this drainage are just hypothezised but not proved yed. Moreover, the number of modeling studies show that the pulse wave powered mechanism is unlikely able to perform pumping as suggested. Thus, the problem is still open. In addition, new data obtained under the conditions of intensionally opened blood brain barrier shows the presence of equally fast transport in opposite durection. In our study we investigate the possible physical mechanisms for rapid transport of substances after the opening of blood-brain barrier under the conditions of zero net flow.

We have proposed a wavelet-domain de-noising technique for imaging of human brain malignant glioma by optical coherence tomography (OCT). It implies OCT image decomposition using the direct fast wavelet transform, thresholding of the obtained wavelet spectrum and further inverse fast wavelet transform for image reconstruction. By selecting both wavelet basis and thresholding procedure, we have found an optimal wavelet filter, which application improves differentiation of the considered brain tissue classes – i.e. malignant glioma and normal/intact tissue. Namely, it allows reducing the scattering noise in the OCT images and retaining signal decrement for each tissue class. Therefore, the observed results reveals the wavelet-domain de-noising as a prospective tool for improved characterization of biological tissue using the OCT.

The development of novel technologies for the imaging of meningeal lymphatic vessels is one of the amazing trends of biophotonics thanks to discovery of brain lymphatics over several years ago. However, there is the limited technologies exist for the study of lymphatics in vivo because lymphatic vessels are transparent with a low speed flow of lymph. Here we demonstrate the successful application of fluorescent microscopy for the imaging of lymphatic system in the mouse brain in vivo.

The assessment of functioning microcirculatory network implies usage of adequate tools for testing the network responses on local changes of vessels state. While there are well-developed and widely used methods, such as focal application of vasoactive substances, or electric stimulation, there is a need for a non-destructive (and ideally – non-contact) and local method of impact a single vessel in order to trigger the network responce.

In this paper, we investigate the possibility of applying the effect of a reversible change in the diameter of a blood vessel caused by laser radiation as a functional test of a microcirculatory system. For this purpose, we combine this effect with the method of micro-PIV (particle image velocimetry), which provides information on both the dynamics of blood flow in neighboring segments and the changes in their diameters.

The estimation of cerebral hemodynamic signals has an important role for monitoring the stage of neurological diseases. Functional Near-Infrared Spectroscopy (fNIRS) can be used for monitoring of brain activities. fNIRS utilizes light in the near-infrared spectrum (650-1000 nm) to study the response of the brain vasculature to the changes in neural activities, called neurovascular coupling, within the cortex when cognitive activation occurs. The neurovascular coupling may be disrupted in the brain pathological condition. Therefore, we can also use fNIRS to diagnosis brain pathological conditions or to monitor the efficacy of related treatments. The Cold pressor test (CPT), followed by immersion of dominant hand or foot in the ice water, can induce cortical activities. The perception of pain induced by CPT can be related to cortical neurovascular coupling. Hence, the variation of cortical hemodynamic signals during CPT can be an indicator for studying neurovascular coupling. Here, we study the effect of pain induced by CPT on the temporal variation of concentration of oxyhemoglobin [HbO₂] and deoxyhemoglobin [Hb] in the healthy brains. We use fNIRS data collected on forehead during a CPT from 11 healthy subjects, and the average data are compared with post-stimulus pain rating scores. The results show that the variation of [Hb] and [HbO₂] are positively correlated with self-reported scores during the CPT. These results depict that fNIRS can be potentially applied to study the decoupling of neurovascular process in brain pathological conditions.

Endothelial cells are cells lining the inner surface of the blood and lymphatic vessels, they separate the bloodstream from the deeper layers of the vascular wall. Earlier endothelium was considered only as a passive barrier between blood and tissues. However, it has now become apparent that endothelial cells, specifically reacting to different molecular signals generated locally and remotely, perform a variety of functions. Simulation of large vascular networks requires the development of specialized models of autoregulation of vascular tone. On the one hand, such models should have a strong support for cellular dynamics, on the other - be as computationally efficient as possible. A model of a two-dimensional cylindrical array of endothelial cells is proposed on the basis of the integral description by means of the whole-cell CVC. The process of propagation of hyperpolarizing and depolarizing pulses is investigated depending on the statistics of cell distribution between the two main types. Endothelial cells are considered as a dynamic system possessing bistability. Based on the articles, the results of the distribution of the resting-potential values were repeated, the propagation of the hyperpolarizing pulse was observed, the endothelial cell chain supported the propagation of the wave switching to a hyperpolarized state, and then the return wave returned to its original state.

Blood platelets play a pivotal role in blood coagulation and in other normal and pathological processes. The understanding of fundamental mechanisms underlying their functions is very important for diagnostics and treatment. Single-cell experiments are needed for this purpose, which are complicated by insufficient spatiotemporal precision of conventional activation protocols. We present an approach to trigger single platelet activation optically, without the need of reagent mixing. This is achieved using photolabile compound, which rapidly delivers epinephrine upon UV irradiation. We demonstrated the applicability of the technique to rapidly induce platelet activation for studying dynamics of activation. The presented method may give novel fundamental knowledge about platelet functions and facilitate current research of their ability to deliver drugs to tumors or vascular injury sites.

The paper is devoted to the analysis of dynamic of interactions between signals of autonomic circulatory regulation. We investigated two-hour experimental records of 30 healthy people. Phase synchronization was studied using the signals of the electrocardiogram and the photoplethysmogram of vessels. We found the presence of long synchronous intervals in some subjects. For analysis of the dynamic we calculated autocorrelation functions. The analysis made it possible to reveal indirect signs of the influence of the humoral regulation system.

Despite the variety of existing methods of recording the human pulse and a solid history of their development, there is still considerable interest in this topic. The development of new non-contact methods, based on advanced image processing, caused a new wave of interest in this issue. We present a simple but quite effective method for analyzing the mechanical pulsations of blood vessels lying close to the surface of the skin. Our technique is a modification of imaging (or remote) photoplethysmography (i-PPG). We supplemented this method with the addition of a laser light source, which made it possible to use other methods of searching for the proposed pulsation zone. During the testing of the method, several series of experiments were carried out with both artificial oscillating objects as well as with the target signal source (human wrist). The obtained results show that our method allows correct interpretation of complex data. To summarize, we proposed and tested an alternative method for the search and analysis of pulsating vessels.

Physiological signals often contain various bad segments that occur due to artifacts, failures of the recording equipment or varying experimental conditions. The related experimental data need to be preprocessed to avoid such parts of recordings. In the case of few bad segments, they can simply be removed from the signal and its analysis is further performed. However, when there are many extracted segments, the internal structure of the analyzed physiological process may be destroyed, and it is unclear whether such signal can be used in diagnostic-related studies. In this paper we address this problem for the case of cerebral vessels dynamics. We perform analysis of simulated data in order to reveal general features of quantifying scaling features of complex signals with distinct correlation properties and show that the effects of data loss are significantly different for experimental data with long-range correlations and anti-correlations. We conclude that the cerebral vessels dynamics is significantly less sensitive to missed data fragments as compared with signals with anti-correlated statistics.

A model of human cardiovascular system is proposed to describe the main heart rhythm, influence of autonomous regulation on frequency and strength of heart contractions and resistance of arterial vessels; process of formation of arterial pressure during systolic and diastolic phases; influence of respiration; synchronization between loops of autonomous regulation. The proposed model is used to simulate the dynamics of heart rate and arterial pressure during passive transition from supine to upright position. Results of mathematical modeling are compared to original experimental data.

The most promising method for the quantitative determination of cardiovascular tone indicators and of cerebral hemodynamics indicators is the method of impedance plethysmography. The accurate determination of these indicators requires the correct identification of the characteristic points in the thoracic impedance plethysmogram and the cranial impedance plethysmogram respectively. An algorithm for automatic analysis of these plethysmogram is presented. The algorithm is based on the hard temporal relationships between the phases of the cardiac cycle and the characteristic points of the plethysmogram. The proposed algorithm does not require estimation of initial data and selection of processing parameters. Use of the method on healthy subjects showed a very low detection error of characteristic points.

Laser Speckle Contrast Analysis (LASCA) is one of the most powerful modern methods for revealing blood dynamics. The experimental design and theory for this method are well established, and the computational recipie is often regarded to be trivial. However, the achieved performance and spatial resolution may considerable differ for different implementations. We comprise a minireview of known approaches to the spatial laser speckle contrast data processing and their realization in MATLAB code providing an explicit correspondence to the mathematical representation, a discussion of available implementations. We also present the algorithm based on the 2D Haar wavelet transform, also supplied with the program code. This new method provides an opportunity to introduce horizontal, vertical and diagonal speckle contrasts; it may be used for processing highly anisotropic images of vascular trees. We provide the comparative analysis of the accuracy of vascular pattern detection and the processing times with a special attention to details of the used MATLAB procedures.

In this paper, we study the complex multi-scale network of nonlocally coupled oscillators for the appearance of chimera states. Chimera is a special state in which, in addition to the asynchronous cluster, there are also completely synchronous parts in the system. We show that the increase of nodes in subgroups leads to the destruction of the synchronous interaction within the common ring and to the narrowing of the chimera region.

We study the features of formation and evolution of spatiotemporal chaotic regime generated by autonomous pacemakers in excitable deterministic and stochastic bistable active media using the example of the FitzHugh – Nagumo biological neuron model under discrete medium conditions. The following possible mechanisms for the formation of autonomous pacemakers have been studied: 1) a temporal external force applied to a small region of the medium, 2) geometry of the solution region (the medium contains regions with Dirichlet or Neumann boundaries). In our work we explore the conditions for the emergence of pacemakers inducing target patterns in a stochastic bistable excitable system and propose the algorithm for their analysis.

The paper considers the phenomena of competition in multiplex network whose structure evolves corresponding to dynamics of it’s elements, forming closed loop of self-learning with the aim to reach the optimal topology. Numerical analysis of proposed model shows that it is possible to obtain scale-invariant structures for corresponding parameters as well as the structures with homogeneous distribution of connections in the layers. Revealed phenomena emerges as the consequence of the self-organization processes related to structure-dynamical selflearning based on homeostasis and homophily, as well as the result of the competition between the network’s layers for optimal topology. It was shown that in the mode of partial and cluster synchronization the network reaches scale-free topology of complex nature that is different from layer to layer. However, in the mode of global synchronization the homogeneous topologies on all layer of the network are observed. This phenomenon is tightly connected with the competitive processes that represent themselves as the natural mechanism of reaching the optimal topology of the links in variety of real-world systems.

Computed tomography is a great source of biomedical data because it allows a detailed exploration of complex anatomical structures. Some structures are not visible on CT scans, and some are hard to distinguish due to partial volume effect. CT datasets require preprocessing before using them as anatomical models in a simulation system. The work describes segmentation and data transformation methods for an anatomical model creation from the CT data. The result models may be used for visual and haptic rendering and drilling simulation in a virtual surgery system.

We discuss the problem of quantifying power-law statistics of complex processes from short signals. Based on the analysis of electroencephalograms (EEG) we compare three interrelated approaches which enable characterization of the power spectral density (PSD) and show that an application of the detrended fluctuation analysis (DFA) or the wavelet-transform modulus maxima (WTMM) method represents a useful way of indirect characterization of the PSD features from short data sets. We conclude that despite DFA- and WTMM-based measures can be obtained from the estimated PSD, these tools outperform the standard spectral analysis when characterization of the analyzed regime should be provided based on a very limited amount of data.

In this paper we study the spiking behaviour of a neuronal network consisting of Rulkov elements. We find that the regularity of this behaviour maximizes at a certain level of environment noise. This effect referred to as coherence resonance is demonstrated in a random complex network of Rulkov neurons. An external stimulus added to some of neurons excites them, and then activates other neurons in the network. The network coherence is also maximized at the certain stimulus amplitude.

Biomedical volumetric data visualization techniques for the exploration purposes are well developed. Most of the known methods are inappropriate for surgery simulation systems due to lack of realism. A segmented data visualization is a well-known approach for the visualization of the structured volumetric data. The research is focused on improvement of the segmented data visualization technique by the aliasing problems resolution and the use of material transparency modeling for better semitransparent structures rendering.

A starlike network of non-identical phase oscillators is considered that contains the hub and tree rays each having a single node. In such network effect of indirect synchronization control is reported: changing the natural frequency and the coupling strength of one of the peripheral oscillators one can switch on an off the synchronization of the others. The controlling oscillator at that is not synchronized with them and has a frequency that is approximately four time higher then the frequency of the synchronization. The parameter planes showing a corresponding synchronization tongue are represented and time dependencies of phase differences are plotted for points within and outside of the tongue.

We study influence of time delay in coupling on the dynamics of two coupled multimode optoelectronic oscillators. We reveal the structure of main synchronization region on the parameter plane and main bifurcations leading to synchronization and multistability formation. The dynamics of the system is studied in a wide range of values of control parameters.

A method for the reconstruction of time-delayed feedback system is investigated, which is based on the detection of synchronous response of a slave time-delay system with respect to the driving from the master system under study. The structure of the driven system is similar to the structure of the studied time-delay system, but the feedback circuit is broken in the driven system. The method efficiency is tested using short and noisy data gained from an electronic chaotic oscillator with time-delayed feedback.

The bone drilling is an essential part of surgeries in ENT and Dentistry. A proper training of drilling machine handling skills is impossible without proper modelling of the drilling process. Utilization of high precision methods like FEM is limited due to the requirement of 1000 Hz update rate for haptic feedback. The study presents a mathematical model of the drilling process that accounts the properties of materials, the geometry and the rotation rate of a burr to compute the removed material volume. The simplicity of the model allows for integrating it in the high-frequency haptic thread. The precision of the model is enough for a virtual surgery system targeted on the training of the basic surgery skills.