A survey of recent applications of the renormalization group approach to the study of thermodynamic Casimir forces near critical points is given. The successes and limitations of approaches based on the small expansion are explained. The large- limit of the theory on a -dimensional strip of thickness bounded by free surfaces is expressed exactly in terms of the eigenvalues and eigenfunctions of a self-consistent Schrödinger equation. Numerically exact results for the scaling functions of the -dependent part of the excess free energy and the Casimir force are presented for . Apart from brute-force Monte Carlo simulations of lattice model, the limit provides so far the only single framework that is (i) capable to deal simultaneously with confined critical fluctuations, confined Goldstone modes, and dimensional crossover and (ii) exhibits qualitatively all features of the temperature-dependent Casimir force known from wetting experiments on He and Monte Carlo simulations of the XY model on a strip.

Temperature Green function formalism is applied to Fermi gas with density-density interaction. New field variables corresponding to critical modes are introduced. Usual mean-field theory can be applied in the new variables. A phase transition is investigated in the framework of the quantum-field renormalization group and -- expansion. As a result it is obtained that the character of the phase transition depends on fermion spin.

Using the field theoretic renormalization group technique we have considered the influence of helicity (violation of spatial parity) on the diffusion processes of a passively advected vector field in the framework of the kinematic MHD turbulence. The influence of the helicity on the turbulent magnetic Prandtl number is analyzed. A nontrivial impact of the internal structure of the advected field on the diffusion rate is analytically described.

The models E and F of critical dynamics are studied in the presence of random velocity fluctuations. Stochastic Navier-Stokes equation is used for generation of velocity fluctuations. Using Martin-Siggia-Rose procedure field theoretic model is obtained. The power counting and analysis of ultraviolet divergences are performed near the critical dimension of the models. Both of them are renormalized in the one-loop approximation using dimensional regularization in minimal subtraction scheme. Influence of hydrodynamic fluctuations on the stability near phase transition point is discussed.

The method of searching the embeddings of Riemannian spaces with a high enough symmetry in a flat ambient space is suggested. It is based on a procedure of construction of surfaces with a given symmetry. The method is used for the construction of the embeddings for the Schwarzschild solution in a six-dimensional ambient space and Friedmann-Robertson-Walker solutions in five-dimensional ambient space. All known variants of embeddings for these spaces are constructed using this method. Two new embeddings for Schwarzschild solution are constructed also. One of these new embeddings is asymptotically flat, while the other embeddings in a six-dimensional ambient space do not have this property. The asymptotically flat embedding can be used in the analysis of the many-body problem, as well as for the development of gravity description as a theory of a surface in a flat ambient space.
/In collaboration with A.A. Sheykin/

It is shown that the energy of two colliding particles in the erhosphere of the rotating black hole close to the horizon can get values close to the Planckean scale.This occurs due to the possibility of very small radial velocity of one of the colliding particle called critical. The mechanism of getting critical particles by previous scattering and receiving critical momentum is analysed.The possibility to explain the observations of Ultra high energetic cosmic particles by the Auger group is discuused.

The Anti-de Sitter space in Beltrami coordinates is mostly suitable for modern cosmology because it has (quasi "Big Bang") singularity. It is possible to find the initial invariant distributions of the matter in the phase space, compute its dependance on time and compare the result with modern observations.

The exactly solvable totally asymmetric models of the low dimensional non-equilibrium physics on a ring, namely the totally asymmetric simple exclusion process and the totally asymmetric simple zero range process, are considered. The Quantum Inverse Method allows to calculate the scalar products of the state vectors of the models and to represent the answers in the determinantal form. It is shown that the eigenvectors of the models form a complete orthogonal basis. The projections of the state vectors on a stationary states are studied.

The structure of a cotangent bundle is investigated for quantum linear groups and . Using a q-version of the Cayley-Hamilton theorem we construct an extension of the algebra of differential operators on (otherwise called the Heisenberg double) by spectral values of the matrix of right invariant vector fields. We calculate an evolution operator for the model of q-deformed isotropic top introduced by A.Alekseev and L.Faddeev. The evolution operator is not uniquely defined and we present two possible expressions for it. The first one is a Riemann theta function in the spectral variables. The second one is an almost free motion evolution operator in terms of logarithms of the spectral variables. Connection between the two operators is given by a modular functional equation for Riemann theta function. Relations to the construction of the Faddeev's modular double will be discussed.

The fermionization method is applied to the vertex models with boundary conditions. The six vertex model with the domain wall boundary conditions has been studied in detail. The two-point correlation function and the generating functional are given in an explicit form in the free fermion point.

I would like to advocate the new approach to gauge theories put forward by Z.Bern, L.Dixon, D.Kosower et al which allows to perform calculations of gauge invariant quantities without fixing the gauge, without ghosts, etc and drastically simplifies the calculations. I will present some explicit examples and discuss the subtleties.

The problem of propagation of photons and massive vector mesons in the presence of Lorenz and CPT invariance violating medium (Chern-Simons action) will be dicussed when the parity-odd medium is separated from the vacuum along a hyperplane. The solutions in different half-spaces will be matched both in the case of space-like and time-like Chern-Simons action by means of the Bogolubov transformations. Accordingly the two different Fock vacua (zero particle states) happen to be mutually coherent states. The boundary mirror/transparency properties will be illustrated. They can help to register local parity violation in a finite volume of heavy ion fireball with nontrivial axial charge and/or of a star with cold axion condensate degrading to its boundary.

Annihilation processes, where the reacting particles are influenced by some external advective field, are one of the simplest examples of nonlinear statistical systems. This type of processes can be observed in miscellaneous chemical, biological or physical systems. In low space dimensions usual description by means of kinetic rate equation is not sufficient and the effect of density fluctuations must be taken into account. Using perturbative renormalization group we study the influence of random velocity field on the kinetics of single-species annihilation reaction at and below its critical dimension . The advecting velocity field is modelled by the self-similar in space Gaussian variable finite correlated in time (Obukhov-Kraichnan model). Effect of the compressibility of velocity field is taken into account and the model is analyzed near its critical dimension by means of three-parameter expansion in and . Here is the deviation from the Kolmogorov scaling, is the deviation from the (critical) space dimension and is the deviation from the parabolic dispersion law. Depending on the value of these exponents and the value of compressiblity parameter , the studied model can exhibit various asymptotic (long-time) regimes corresponding to the infrared fixed points of the renormalization group. The possible regimes are summarized and the decay rates for the mean particle number are calculated in the next-to-leading order of the pertubartion theory.

The use of field-theoretical renormalization group for studying the critical behavior of two-dimensional models is discussed. RG functions of 2D Euclidean n-vector theory are written down up to the five-loop terms. Numerical estimates obtained from these series by means of Pade-Borel-Leroy resummation are presented and compared with their exact counterparts known for . Then, pseudo-\epsilon expansions for the fixed-point location , critical exponents and effective sextic coupling are derived from the five-loop series. The expansions obtained are shown to be more "friendly" than original RG series: the higher-order coefficients of pseudo- expansions for , and turn out to be much smaller than their RG analogs. This enables one to resum the pseudo- expansions using simple Pade approximants, without addressing Borel-Leroy transformation. Moreover, numerical estimates extracted from the pseudo- expansions are found to be closer to the known exact values than those obtained from Pade-Borel-Leroy resummed five-loop RG series.

Using perturbative renormalization group we study the influence of random velocity field on the critical behaviour of directed bond percolation process near its second-order phase transition between absorbing and active phase. Kraichnan model with finite correlation time for modelling advecting velocity field is considered. The field-theoretic renormalization group is applied to obtain information about possible large scale behavior. The model is analyzed near its critical dimension by means of three-parameter expansion in , where is the deviation from the Kolmogorov scaling, the deviation from the critical space dimension and is the deviation from the parabolic dispersion law for the velocity correlator. Fixed points with corresponding regions of stability are evaluated to the second order in the perturbation scheme.

The solution of the famous Yang-Baxter equation is usually called R-matrix. We consider the general scheme of the construction of R-matrix as integral operator. Everything is illustrated on various examples starting from the simplest one related to conformal group in two and four dimensions and finishing by elliptic Sklyanin algebra.

Several modern methods for calculation of complicated master integrals are considered. The results for two-loop two-point and three-point diagrams are presented.

The approach to construction of generalized oscillators, conтected with polynomials orthogonal on the real line, proposed by one of the authors early, is extended to the case of bivariate polynomials. As an interesting example we consider bivariate Chebyshev polynomials connected with representations of sl(3) algebra. The polynomials under consideration form a complete system of eigenfunctions for a quadratic hamiltonian of the constructed two-dimensional Chebyshev oscillator. We show that the concidered hamiltonian splits into the sum of two commuting operators. So it is possible to consider this two-dimensional oscillator as a system of two noninteracting oscillators. The generators of the first oscillator algebra leave invariant all N-particles subspaces while the ladder operators of the second oscillator transform neighboring N-particles subspaces one onto another. The representation of the ladder operators for two-dimensional Chebyshev oscillator by differential operators are given. /In collaboration with E.V.Damaskinski/

We treat electron transport in graphene monolayer close to the charge neutrality point in the framework of the Dirac model with unstable vacuum. It is shown that the electron-hole creation by a slowly varying electric field is crucial for the understanding the conductivity of graphene at low carrier density beyond linear response in dc. Self-consistent solutions of the Dirac-Maxwell equations are found in the strong-field approximation and it is shown that such solutions can explain the experimental super linear I-V of graphene devices that were recently reported. /In collaboration with D.M. Gitman, and N. Yokomizo from University of Sao Paulo, Brazil/

The proton form factors at large momentum transfer are dominated by two contributions which are associated with the hard and soft rescattering respectively. Motivated by a very active experimental form factor program at intermediate values of momentum transfers, , where an understanding in terms of only a hard scattering mechanism cannot yet be expected, we investigate the soft spectator scattering contribution using soft collinear effective theory (SCET). Within such description, the form factor is characterized, besides the hard scale , by a hard-collinear scale , which arises due to presence of soft spectators, with virtuality (~GeV), such that . In case of nucleon FFs the soft spectator scattering mechanism contributes at leading power accuracy and therefore it must be considered in the systematic QCD factorization approach. We discuss how to generalize the well known Brodsky-Lepage collinear factorization in order to include the soft spectator contribution using the SCET technique. We carefully investigated the factorization of the soft and collinear modes in this case and demonstrate that even for the FF the pure collinear factorization could not be valid due to the mixing of the soft and collinear contributions. As a result this allows one to put specific constrains on the end-point behavior of the nucleon distribution amplitude.

We consider noncommutativity of space-time which is covariant under the action of quantum-deformed Poincare algebra. The consistency with quantum deformation of Poincare algebra leads to braided noncommutative field theory, with braid factors determined by the universal R-matrix. The technique will be applied to the kappa-deformation of free quantum fields.

The notion of supersymmetry is used with at least three meanings. They are explained and several applications in solid state physics are given.

We develop a general approach to the calculation of kinematic corrections , in hard processes which involve momentum transfer from the initial to the final hadron state. As the principal result, the complete expression is derived for the time-ordered product of two electromagnetic currents that includes all kinematic corrections to twist-four accuracy. The results are immediately applicable e.g. to the studies of deeply-virtual Compton scattering.

We compare the predictions of high density QCD with the LHC data on multiplicities,jet quenching,elliptic flow and long range rapidity correlations Please, write formulae using TeX notation, e.g:

Physical quantity known as the helicity is a signature of parity violation in developed (magneto)hydrodynamic turbulence. The helicity similarly to the energy is time invariant in inviscid fluid. It occurs naturally in a wide variety of geophysical flows. In three-dimensional turbulence, contrary to the 2D turbulence, there is a joint direct cascade of both energy and helicity simultaneously from large to small scales. It can affect celebrated Kolmogorov scaling and change non-universal amplitudes like Kolmogorov constant in kinetic energy specrum and effective viscosity coefficient. A set of self-consistent equations in a statistical model of fully developed homogeneous isotropic turbulence, which is based on the principle of maximal randomness of the velocity field with a given energy spectral flux, is derived. The formal solution of these equations yields Kolmogorov exponents, but this solution leads to negative Kolmogorov constant and negative viscosity, which thereafter induces the instability of turbulent flows. It has been demonstrated, that stable solution can be obtained by means of specific mechanism of spontaneous parity breaking. More precisely, the Kolmogorov solution becomes stable (both and effective viscosity become positive) if one takes into account all allowable tensor structures in the statistical averages of velocity fluctuations including pseudotensor contributions generated by helicity. This solution predicts a large, closed to limit value ( equal to one), helical coefficient in inertial range. The relationship obtained between and confirms this conclusion for the experimental value of .

Within the framework of the field theoretic renormalization group (RG), a stochastic equation, proposed earlier for describing a randomly growing interface, is studied. Correct analysis of the model requires introducing an infinite number of counterterms and corresponding coupling constants (charges). The one-loop expression for the countertems can be obtained in a closed form, which allows to derive and analyze the corresponding RG equation. The analysis shows that the infinite-dimensional space of charges involves a two-dimensional surface of fixed points of the RG equations. It involves a region of infrared stability, and the model can exhibit a scale invariance with non-universal critical dimensions of the interface height and the time, satisfying a certain exact relation.

We outline the spectacular developments in the field of multiloop renormalization group from yearly 70-th till present time, which are heavily based on a proper use of versatile possibilities of the Bogouibov R-operation within dimensional regularization framework. Then we discuss future perspectives and open problems of multiloop RG calculations. In particular, we discuss the possibility of extending the 30-years old analytical calculation of the five-loop beta function in the model on the six-loop level.

Sector decomposition method proved to be a great tool for the numerical computation of Feynman integrals. This approach significantly improves the convergence of the integrals at the cost of a substantial increase in the integrand (proportional to number of sectors). Therefore, one of the most important tasks is to find strategies, which give the minimal number of sectors. For diagrams at zero external momenta, the best of known strategies is a Speer sectors strategy. The proposed algorithm can reduce the number of sectors compared to Speer sectors for the most of Feynman diagrams of this type and hence reduce the time required for decomposition and integration.

Casimir energy of two gratings separated by a vacuum slit is expressed in terms of Rayleigh coefficients. Comparison of theory and experiments is performed.

As a basic conception for the modeling of interaction of a macroscopic material body with quantum fields is considered the proposed by K.Symanzik approach. Its application in quantum electrodynamics enables one to establish the most general form of the action functional describing the interaction of 2-dimensional material surface with photon field. The models making it possible to calculate for thin films from non-ideal conducting material the Casimir energy, Casimir-Polder potential, characteristics of scattering processes and investigation of magneto- and electrostatic phenomena are presented. The specific of regularization and renormalization procedures used by calculations and the physical meaning of obtained results are discussed.

The approach to construction of generalized oscillators, conтected with polynomials orthogonal on the real line, proposed by one of the authors early, is extended to the case of bivariate polynomials. As an interesting example we consider bivariate Chebyshev polynomials connected with representations of sl(3) algebra. The polynomials under consideration form a complete system of eigenfunctions for a quadratic hamiltonian of the constructed two-dimensional Chebyshev oscillator. We show that the concidered hamiltonian splits into the sum of two commuting operators. So it is possible to consider this two-dimensional oscillator as a system of two noninteracting oscillators. The generators of the first oscillator algebra leave invariant all N-particles subspaces while the ladder operators of the second oscillator transform neighboring N-particles subspaces one onto another. The representation of the ladder operators for two-dimensional Chebyshev oscillator by differential operators are given. /In collaboration with V.V.Borzov/

Composite superconformal string model is constructed to describe hadron amplitudes and hadron spectrum. It belongs to a new class of dual (string) models to be out of usual approaches. This model considers string scale alpha' as hadron scale is of order 1 GeV^(-2). There is supersymmetry on world two-dimensional surface only. Interacton amplitudes for pi-, K-mesons and nucleons are calculated within the framework of this model. The way of elimination of degeneracy in parity for baryons is found. (Joint report with A.N.Semenova.)

Recent years, the models of beyond standard model physics based on the hypothesis that our universe is a four-dimensional space-time hypersurface (3-brane) embedded in a fundamental multi-dimensional space have become quite popular. The influence of gravity is especially interesting, playing an important role in a (de) localization of matter fields on the brane. The question arises, under what circumstances the localization of matter fields with spin zero on a brane is still possible when the minimal interaction with gravity is present? We consider a model of the domain wall formation with finite thickness (”thick” branes) by two scalar fields with quartic interaction characterized by mass and gravity in five-dimensional noncompact space-time, but with modified behavior of a metric factor in the bulk on both sides of the brane, i.e. with different anti-de Sitter geometries. Brane formation is possible when scalar field vacuum configurations have nontrivial topology. There are two types of these configurations corresponding to massless phase and massive phase controlled by mass parameter . Due to mixing of scalar and gravitational degrees of freedom nonperturbative corrections of equations on second variation arise. This changes the spectrum in the scalar sector dramatically in comparison with model without gravity. Zero-mode exists in the case of asymmetric linear term with coefficient or in the case of symmetric linear term (defect of cosmological constant) with arbitrary positive coefficient or negative coefficient satisfying the same condition as in the case of asymmetric term. Also there is a light scalar state. Its leading order mass is the same as in the similar model without gravity.
/In collaboration with Alexander and Vladimir Andrianov/