**Maxim V. Shamolin**, Lomonosov Moscow State University, Institute of Mechanics, Russian Federation

**Abstract:**In this activity, we systematize some results on the study of the equations of spatial motion of dynamically symmetric fixed rigid bodies-pendulums located in a nonconservative force fields. The form of these equations is taken from the dynamics of real fixed rigid bodies placed in a homogeneous flow of a medium. In parallel, we study the problem of a spatial motion of a free rigid body also located in a similar force fields. Herewith, this free rigid body is influenced by a nonconservative tracing force; under action of this force, either the magnitude of the velocity of some characteristic point of the body remains constant, which means that the system possesses a nonintegrable servo constraint, or the center of mass of the body moves rectilinearly and uniformly; this means that there exists a nonconservative couple of forces in the system. Earlier, the author already proved the complete integrability of the equations of a plane-parallel motion of a fixed rigid body-pendulum in a homogeneous flow of a medium under the jet flow conditions when the system of dynamical equations possesses a first integral, which is a transcendental (in the sense of the theory of functions of a complex variable, i.e., it has essential singularities) function of quasi-velocities. It was assumed that the interaction of the medium with the body is concentrated on a part of the surface of the body that has the form of a (one-dimensional) plate. In sequel, the planar problem was generalized to the spatial (three-dimensional) case, where the system of dynamical equations has a complete set of transcendental first integrals. It was assumed that the interaction of the homogeneous medium flow with the fixed body (the spherical pendulum) is concentrated on a part of the body surface that has the form of a planar (two-dimensional) disk. In this activuty, the results relate to the case where all interaction of the homogeneous flow of a medium with the fixed body is concentrated on that part of the surface of the body, which has the form of a two-dimensional disk, and the action of the force is concentrated in a direction perpendicular to this disk. These results are systematized and are presented in invariant form.

**Acknowledgements:**Lomonosov Moscow State University

1

**Marcin Lawnik**, Silesian University of Technology, Faculty of Applied Mathematics, Poland

**Abstract:**The logistic map is commonly used in, for example, chaos based cryptography. However, its properties do not render a safe construction of encryption algorithms. Thus, the scope of the paper is a proposal of the generalization of a logistic map by means of a well-recognized family of chaotic transformations. In the next step, an analysis of Lyapunov exponent and the distribution of a iterative variable are studied. The obtained results confirm that the analyzed model can safely and effectively replace a classic logistic map for applications involving chaotic cryptography.

2

**Bartek Wierzba**, Rzeszow University of Technology, Materials Science, Poland

**Abstract:**The complete understanding of the phenomenological process related to the evolution of the ternary multiphase systems is still lacking. In this paper the ternary interdiffusion process in multiphase system is discussed by means of numerical modeling. The entropy production principle is proposed to chose the proper diffusion path during the process (e.g. internal oxidation). The simulations present the local entropy production curve that determines the diffusion path in Ni-Al-Ti ternary diffusion couple.

3

**Shabnam Jamshidzadeh**, Young Researchers and Elite Club, Ardabil Branch, Islamic Azad University, Ardabil, Iran, Mathematics, Iran (Islamic Republic of)

**Abstract:**In this paper, we applied the (G'/G)-expansion method to construct more general exact solutions of three special types of Boussinesq equation: namely Boussinesq equation, improved Boussinesq equation and variant Boussinesq equation, where the French scientist \textit{Joseph Valentin Boussinesq} (1842–-1929) described in the 1870’s model equations for the propagation of long waves on the surface of water with a small amplitude. Our work is motivated by the fact that the (G'/G)-expansion method provides not only more general forms of solutions but also periodic, solitary waves and rational solutions. The method appears to be easier and faster by means of a symbolic computation.

4

**Spyridon Vossos**, NKUA, Chemistry, Greece

**Elias Vossos**, NKUA, Physics, Greece

**Abstract:**Relativity Theory (RT) and the corresponding Relativistic Quantum Mechanics (RQMs) are the fundamental theories of physics. Special Relativity (SR) relates the frames of Relativistic Inertial observers (RIOs), through Linear Spacetime Transformation (LSTT) of linear spacetime. Classic Relativity (CR) uses real spacetime endowed with Lorentz metric and two frames of RIOs with parallel spatial axes are related through Lorentz Boost (LB). This cancels the transitive attribute in parallelism, when three RIOs are related, because LB is not closed transformation, causing Thomas Rotation. In this presentation, we consider a closed LSTT of Complex Spacetimes (CSTs), so there is no necessity for spatial axes rotation and all the frames are chosen having parallel spatial axes. The solution is expressed by a 4x4 matrix with elements containing components of the complex velocity of one RIO (βx, βy, βz) wrt another RIO and two functions of its norm (h, λ). Moreover, it emerges an equation (condition) relating the above parameters. So, it emerges a class of metrics which are in accordance with the condition. Choosing one of the acceptable metrics or an invariant quantity (isometry), the problem is totally solved. Sometimes the elements of the 4x4 matrix are complex numbers and the corresponding spacetime is necessarily complex and sometimes not. Thus, we have infinite number of closed LSTTs, each one with the corresponding special relativity theory. The first solution is Galileo Transformation (GT) endowed with Euclidean metric of real space and the invariant time, in which any other closed LSTT is reduced, if one RIO has small velocity wrt another RIO. A second solution is Vossos transformation (VT) endowed with Lorentz metric of complex spacetime and invariant spacetime interval (or equivalently invariant speed of light in vacuum – c), which produce the theory of Euclidean Complex Relativistic Mechanics (ECRMs). Some other closed LSTTs endowed with their acceptable metric of spacetime and their invariant quantity, are presented. In case that we relate observers with different metric of spacetime, we have the case of General Relativity (GR). The application of this kind of transformations to the SR and GR is obvious. But, the results may be applied to any complex or real spaces of dimension four endowed with the same or different metric, whose elements (four-vectors) have spatial part (vector) with Euclidean metric.

5

**Ioannis ZOIS**, PPC, TRSC, Greece

**Abstract:**We present some ideas for a possible Noncommutative Topological Quantum Field Theory (NCTQFT for short) and Noncommutative Floer Homology (NCFH for short). Our motivation is two-fold and it comes both from physics and mathematics: On the one hand we argue that NCTQFT is the correct mathematical framework for a quantum field theory of all known interactions in nature (including gravity). On the other hand we hope that a possible NCFH will apply to practically every 3-manifold (and not only to homology 3-spheres as ordinary Floer Homology currently does). The two motivations are closely related since, at least in the commutative case, Floer Homology Groups constitute the space of quantum observables of (3+1)-dim Topological Quantum Field Theory. Towards this goal we define some new invariants for 3-manifolds using the space of taut codim-1 foliations modulo coarse isotopy along with various techniques from noncommutative geometry.

6

**Reza Abazari**, Young Researchers and Elite Club, Ardabil Branch, Islamic Azad University, Ardabil, Iran, Department of Mathematics, Iran (Islamic Republic of)

**Abstract:**In this paper, the (G'/G)-expansion method is proposed to construct the exact traveling solutions of a second order wave equation of KdV type: u_t-u_{xxx}+\frac{1}{8}u_x^3-u_x(p\,e^u+q\,e^{-u})=0, where pq=/= 0. Our work is motivated by the fact that the (G'/G)-expansion method provides not only more general forms of solutions but also periodic and solitary waves. By using this method, new exact solutions involving parameters, expressed by three types of functions which are hyperbolic, trigonometric and rational function solutions, are obtained. The method appears to be easier and faster by means of a symbolic computation system.

7

**Jaykov Foukzon**, Israel Institute of Technology, Haifa, Israel., Center for Mathematical Sciences, Israel Institute of Technology, Haifa, Israel., Israel

**Abstract:**In his famous thought experiment,Schrôdinger(1935) imagined a cat that measures the value of an quantum mechanical observable with its life. Since Schrödinger's time, no any interpretations or extensions of quantum mechanics have been proposed which gives clear unambiguous answers to the questions posed by Schrödinger's cat of how long superpositions last and when (or whether) they collapse? In this paper appropriate modifications of quantum mechanics are proposed. We claim that canonical interpretation of the wave function ψ=c₁ψ_1+ c₂ψ_ 2 is correct only when the supports the wave functions ψ_1 and ψ_2 essentially overlap. When the wave functions ψ_1 and ψ_2 have separated supports (as in the case of the experiment that we are considering in this paper) we claim that canonical interpretation of the wave function ψ=c₁ψ_1+ c₂ψ_2 is no longer valid for a such cat state. Possible solution of the Schrödinger's cat paradox are considered.We pointed out that the collapsed state of the cat always shows definite and predictable outcomes even if cat also consists of a superposition [16]-[17] : |cat>=c₁|live cat>+c₂|death cat>.

**Acknowledgements:**-

8

**BOUSSAHEL MOUNIR**, M'sila University, Physics, Algeria

**Abstract:**The construction of the (331) model was successfully performed within a framework of the graded Lie algebra NCG. In input, the Dirac representation can be the most general, by adopting the condition of Wulkenhaar, the (331) gauge models are obtained with new physical properties in addition to classical models, if not mathematical model building in the case of NCG itself is the source of a new physics especially concerning the Higgs.

9

**Hossein Ghaffarnejad**, Semnan University, Physics, Iran (Islamic Republic of)

**Abstract:**Two dimensional analogue of vacuum sector of the Brans Dicke gravity [1] is studied to obtain dynamics of anisotropic spherically symmetric perfect fluid. Solutions of dynamical field equations are obtained in terms of time and radial coordinates. In static regime the obtained solutions describe dark matter fluid regime with state equation $\gamma=\frac{p(\rho)}{\varrho}=-0.25.$ For non-static regime the fluid can be treat as regular matter with positive barotropic index $\gamma>0.$ Evaluation of total mass of the fluid leads to choose particular values on the Brans Dick parameter as $\omega>\frac{2}{3};\omega<-1$. Positions of the horizons (event and apparent) of the collapsing fluid are obtained in its static regime. In case $\omega>0$ the apparent horizon is covered by event horizon where the cosmic censorship hypothesis is still valid. \\ In second part of the paper we obtain de Broglie pilot wave of our fluid model. It can be describe particles ensemble which are distinguished from each other via different values $\omega.$ Incident current density of particles ensemble is evaluated on the event and apparent horizons describing which can be called as the `Hawking radiation`. The quantum potential is calculated on the event (apparent) horizon which is independent (dependent) to values of the Brans Dicke parameter $\omega.$

10

**Won-Kwang Park**, Kookmin University, Mathematics, Korea, Republic of

**Abstract:**Throughout various works, it has been confirmed that only a small number of directions of incident fields is needed for applying topological derivative. In this contribution, we consider topological derivative based technique with a small number of such directions for imaging of thin, curve-like dielectric inhomogeneities embedded in a homogeneous domain and explore a suitable condition via various results of numerical simulations.

11

**Liang Zhang**, Jiangsu Normal University, School of Mechanical and Electrical Engineering, China

**Zhi-quan Liu**, Chinese Academy of Sciences, , China

**liang zhang**, Jiangsu Normal University, School of Mechanical and Electrical Engineering, China

**Yu-tong Ji**, Jiangsu Vocational Institute of Architectural Technology, , China

**Abstract:**Anand constitutive relation of SnAgCu and SnAgCu-nano Al solders were studied under uniaxial tension, and the constitutive model was used in the finite element simulation to analyze the stress-strain response of lead-free solder joints in 3D IC devices. The results showed that the nine parameters of the Anand model can be determined from separated constitutive relations and experimental results. Based on Anand model, the finite element method was selected to calculate the stress-strain response of lead-free solder joints, it was found that in the 3D IC device the maximum stress-strain concentrated in the concern solder joints, the stress-strain of SnAgCu-nano Al solder joints was lower than that of SnAgCu solder joints, which represented that the addition of nano Al particles can enhance the reliability of lead-free solder joints in 3D IC devices.

12

**Jorge Linares**, Université de Versailles-St. Quentin en Yvelines, GEMaC, France

**Abstract:**C. Jureschi, J. Linares, P. R. Dahoo, Yasser Alayli ABSTRACT Monte Carlo methods are one of the most utilized methods by researches in order to model probabilistic or stochastic systems whose analytical solutions are too complicated or impossible to obtain. In this paper we present the Monte Carlo entropic sampling [1] applied to an Ising-like model for 2D and 3D system in order to show the interaction influence of the edge molecules of system with their local environment. To take into account the interaction between the edge molecules with the environment, we use the Hamiltonian system proposed in a previous work [2]: We show that, as for the 1D [2], 2D [3] spin crossover (SCO) systems, the origin of multi steps transition in 3D SCO is the effect of the edge interaction molecules with local environment. Another important result worth noting is the co-existence of step transition with hysteresis and without hysteresis and that by increasing the value of the edge interaction, L, the transition is shifted to the lower temperatures: it means that the role of edge interaction is equivalent to a applied negative pressure. We analyze also, in this contribution, the role of the short- and long-range interaction, J respectively G, with respect to the environment interaction, L in equation (1). 1. Linares, J.; Enachescu, C.; Boukheddaden, K.; Varret, F., Monte Carlo entropic sampling applied to spin crossover solids: the squareness of the thermal hysteresis loop. Polyhedron 2003, 22, (14-17), 2453-2456. 2. Chiruta, D.; Jureschi, C.-M.; Linares, J.; Dahoo, P.; Garcia, Y.; Rotaru, A., On the origin of multi-step spin transition behaviour in 1D nanoparticles. The European Physical Journal B 2015, 88, (9), 1-5. 3. Jureschi, C.-M.; Pottier, B.-L.; Linares, J.; Dahoo, P.R.; Alayli, Y.; Rotaru, A., Simulation of multi-steps thermal transition in 2D spin-crossover nanoparticles. Physica B: Condensed Matter,

**Acknowledgements:**CHAIR Materials Simulation and Engineering, UVSQ, Université Paris Saclay, 78035 Versailles Cedex, France.

13

**Jorge Linares**, Université de Versailles-St. Quentin en Yvelines, GEMaC, France

**Catalin Jureschi**, Université de Versailles-St. Quentin en Yvelines, LISV, France

**Daniel Chiruta**, Université de Versailles-St. Quentin en Yvelines, GEMaC, France

**Yasser Alayli**, Université de Versailles-St. Quentin en Yvelines, LISV, France

**Corneliu Turcu**, University of Suceava, Electrical Engineering and Computer Science, Romania

**Pierre Dahoo**, Université de Versailles-St. Quentin en Yvelines, LATMOS CNRS-UMR 8190, France

**Abstract:**Spin crossover (SCO) [1] is a symbolic example of molecular bistability, in which a diamagnetic low-spin (LS) and a paramagnetic high spin (HS) states are switchable by different physical causes such as magnetic field, pressure, light, and temperature or electrical filed [1]. The SCO materials are formed by 3d transition metal ions (3d4 - 3d7) in an octahedral coordination environment, which display a spin transition between the two states. In recent years, the SCO phenomenon has drawn an increasing attention from physics and chemist points of view, because of potential applications in future devices as sensor, display or data storage devices. In this paper we use a recent model called atom-phonon coupling model [2,3] to explain and illustrate the behavior of a linear chain of molecules. It is well known that besides the thermodynamic parameters which influence the hysteretic behavior of SCO compounds, the system’s architecture state configuration also plays a determinant role. The analysis of the system’s behavior has been studied using Free Energy method, and Monte Carlo Metropolis (MC) method which take into account the phonon contribution. In particular we compare the MC algorithm and the dynamic-matrix method and we determine how the thermal behavior of a 1D SCO system varies as a function of different factors. Furthermore we have blocked the edge atoms of the chain in high spin state to study the system’s behavior. The variation of entropy and heat capacity are analyzed in the ferro-like case for both architectures.

**Acknowledgements:**CHAIR Materials Simulation and Engineering, UVSQ, Université Paris Saclay is gratefully acknowledged

14

**Jaykov Foukzon**, Israel Institute of Technology, Haifa, Israel., Center for Mathematical Sciences, Israel Institute of Technology, Haifa, Israel., Israel

**Abstract:**Exact quasi-classical asymptotic of solutions to the n-dimensional Schrödinger equation with a non-Hermitian Hamiltonian [1] beyond WKB-theory and beyond Maslov's canonical operator theory is presented. We obtain an important generalization of exact quasi-classical asymptotic given in paper [2]. We apply our results to continuous measurements of forced and damped quantum oscillators. [1] M. B.Mensky, Continuous Quantum Measurements and Path Integrals, CRC Press, 1 Jan 1993 - Science - 188 pp. [2] J.Foukzon, E.Men'kova, A.Potapov and S.Podosenov, New exact quasi-classical asymptotic beyond WKB approximation and beyond Maslov formal expansion.

15

**Ju Chunhua**, Harbin Institute of Technology, Department of Materials Science, China

**Abstract:**GRP pipes for its excellent performance has been widely used for oil and gas gathering pipe. Because CO2 is used as displacing agent, causing corrosion problems CO2 FRP pipe, so the research GFRP pipe corrosion mechanism and performance in CO2 media evolution is important. Based on material corrosion theory summary FRP pipe used in CO2 environment, we analyze the main form and the specific process material corrosion, the use of aromatic amine curing and acid anhydride curing two fiberglass pipe as a test tube, analog GFRP pipe containing crude oil, the next CO2 and water medium corrosion test, by surface morphology material, the glass transition temperature changes, and performance to evaluate physical corrosion mechanism of FRP in CO2 in; the test results show that CO2 environment, glass and steel addition to the visible Exterior changes, its glass transition temperature, Pap hardness, ring stiffness, density, and curing of the resin content showed a downward trend.

16

**Marco Zoli**, University of Camerino, School of Science and Technology, Italy

**Abstract:**Many experiments have demonstrated that DNA exhibits sequence-dependent bending and torsional flexibility which influence the DNA-protein interplay and are relevant for biological processes like transcription and replication in living cells. DNA flexibility is tested by a path integral method [1-4] which analyses the propensity of short DNA sequences to convert to the circular form. The partition function for the sub-ensemble of closed molecules is computed by imposing chain ends boundary conditions both on the radial fluctuations and on the angular degrees of freedom. The cyclization probability, the J-factor, proves to be highly sensitive to the model potential, mostly to the non-linear stacking parameters. We find that the J -factor generally decreases by reducing the sequence length ( N ) and, more significantly, below N = 100 base pairs. However, even for very small molecules, the J-factors remain sizeable in line with recent data provided by fluorescence resonance energy transfer assays. Large bending angles between adjacent base pairs and anharmonic stacking appear as the causes of the helix flexibility at short length scales. I also discuss some preliminary results regarding the single molecule stretching under the effect of an external force which samples the DNA mechanical properties at the nano-metric scale. [1] M. Zoli, The Journal of Chemical Physics 141, 174112 (2014). [2] M. Zoli, Soft Matter 10, 4304-4311 (2014). [3] M. Zoli, Journal of Physics: Conference Series 574, 012094 (2015). [4] M. Zoli, Europhysics Letters 110, 18001 (2015).

17

**Leonid Son**, Ural Federal University, Physical - technical institute, Russian Federation

**Dmitry Shulgin**, Ural Federal University, Physical - technical institute, Russian Federation

**Olga Ogluzdina**, Ural Federal University, Physical - technical institute, Russian Federation

**Abstract:**A wide variety of systems may be described by specific dependence, which is known as logistic curve, or S-curve, between the internal characteristic and the external parameter. Linear feedback between these two values may be suggested for a wide set of systems also. In present paper, we suggest a bifurcation behavior for systems with both features, and discuss it for two cases, which are the Ising magnet in external field, and the development of innovative enterprise.

**Acknowledgements:**The work was supported by Russian Ministry of Science and Education (Government job 2014/392 projects 2391 and 1177), and by Goverment of the Russian Federation, act 211 contract № 02.A03.21.0006.

18

**Rubén Danilo BOURDON GARCÍA**, Saint Matthew Fundation for the High Education, Technology Development Center - CDTEC, Colombia

**Javier BURGOS SALCEDO**, Ciinas Corporation, Research and Development, Colombia

**Abstract:**Actually, there is a huge necessity to generate novel strategies to social-ecological systems analyses for resolving global sustainability problems. This dissertation has as main purpose the application of the formal concept analysis to formalize the theory of Augusto Angel Maya, who without doubt, was one of the most important environmental philosophers in South America; Angel Maya proposed and established that Ecosystem - Culture relations, instead Human - Nature ones, is determinant in our understanding and management of natural resources. Based on this, a concept lattice, formal concepts, subconcept - super- concept relations, partially ordered sets, supremum and infimum of the lattice and implications between attributes (Duquenne - Guigues base), were determined for the ecosystem - culture relations.

19

**Vladimir Zeitlin**, Ecole Normale Supérieure, Laboratory of Dynamical Meteorology, France

**Abstract:**Thermodynamics of the moist air is very complicated so, for the purposes of weather and climate modeling, various simplifications and parameterzations are being used in general circulation models of the atmosphere. Dynamics of atmospheric jets and vortices which, together with inertia-gravity waves, constitute the principal dynamical entities of large-scale atmospheric motions, is well-represented in the framework of one- or multi-layer rotating shallow water models, which are obtained by vertically averaging of full "primitive" equations. We show that by adding a simple relaxational parameterization of precipitations and coupling the precipitation with convective fluxes with the help of moist enthalpy conservation results in moist-convective rotating shallow water model ( Bouchut, Lambaerts, Lapeyre and Zeitlin, 2009, Phys. Fluids, v. 21, 116604; Lambaerts, Lapeyre, Zeitlin and Bouchut, 2011, Phys Fluids, v. 23, 046603) which provides an optimal compromise between simplicity and fidelity of reproducing the moist-precipitating dynamics. Numerical simulations with a well-balanced shock-resolving finite-volume numerical scheme (Bouchut and Zeitlin, 2010, Disc. Cont. Dyn. Sys. B, v.13, 739) capture surprisingly well the life-cycle of the moist baroclinic instability (Lambaerts, Lapeyre and Zeitlin, J. Atmos. Sci. 2012, v. 69, 1405). The life-cycle of tropical cyclones is also well-reproduced within the model (N. Lahaye and V. Zeitlin, 2016, J. Atmos. Sci. , to appear). We will give further examples of the applications of the model.

20

**Jaqueline Silva**, Universidade Federal dos Vales do Jequitinhonha e Mucuri, Instituto de Ciência, Engenharia e Tecnologia, Brazil

**Abstract:**We present a computational model for the tree growth process in flooded areas of the Amazon Forest that contemplates seed dispersion processes and influence of light in the model. These processes, strongly influenced by the annual food, occur early in the plant's life cycle and affect the distribution, structure and dynamics of tree populations. The study of these processes is very important to understand the dynamics of the flooded ecosystems and their sustainable management. To study the spatial distribution of seeds in a region and the growth of germinated seeds, we define an influence region for a mature tree in the growth of the seeds, the growth neighborhood of young trees and the probability of a seed to fall at any point of this neighborhood. We make the hypothesis that every year each mature tree produces affixed number of seeds. The seeds are distributed in accordance with a probability that varies with the distance from a mature tree. Seed-germination and the growth process depend on access to the solar light. In particular, since the canopies of mature-trees in a neighborhood of a young-tree contain sufficient dense leaves and twigs structure, the light shall vary and influence the growth of surrounding trees. Indeed, for most species, the nearer a young-tree is from a mature-tree, smaller will be its possibility to survive. In some cases it may occur the death of these individuals by shading.

21

**Tomasz Blachowicz**, Silesian University of Technology, Institute of Physics, Poland

**Andrea Ehrmann**, Bielefeld University of Applied Sciences, Department of Engineering Sciences and Mathematics, Germany

**Abstract:**While the properties of conductive fibers and coatings on textiles can easily be measured and calculated, magnetic coatings of fibers, yarns and fabrics still lack descriptions of their physical properties. Since magnetic textiles can be used for a variety of applications, from magnetic filters to invisible water-marks to magnetic coils and sensors, simulations would be supportive to understand and utilize their properties. The articles gives an overview of different coatings on textile fibers, varying the magnetic materials as well as the fiber composition, giving rise to the interactions between neighboring coated fibers. In this way, it is possible to understand the strong shape anisotropy which must be taken into account when the magnetic properties of textiles are to be tailored. Additionally, the differences between several possible magnetic coating materials become visible. This study can help adjusting the magnetic properties of textile fabrics to a desired application.

22

**Tomasz Blachowicz**, Silesian University of Technology, Institute of Physics, Poland

**Andrea Ehrmann**, Bielefeld University of Applied Sciences, Department of Engineering Sciences and Mathematics, Germany

**Abstract:**Magnetic nano-structures can be used in various applications. Due to their possible utilization in data storage media, examinations of nano-structured systems often aim at decreasing the pattern size, in order to enhance the possible information density in a given area. Since this scaling process is limited by the resolution of the lithography process which is used to produce the nano-particles, it is important to determine the influence of erroneous shape modifications on the magnetic properties, such as magnetization reversal processes and coercive fields. For this, a square nano-wire system from permalloy has been simulated using Magpar. In a former work, changes of the wire diameter have been shown to result in different magnetization reversal mechanisms and significantly altered coercive fields. In a new project, the intersections of the wires – which are most susceptible to undesired shape modifications – have been changed by adding or subtracting parts. Additionally, the wire intersections have been separated step by step, resulting in a qualitatively changed angular dependence of the coercive fields. Similar experiments have been performed for nano-squares with walls of rectangular cross-section. This study allows for estimation of the reliability of magnetic properties of nano-structures with respect to undesired shape modifications in the lithography process.

23

**Alexey Kondikov**, Peter the Great St. Petersburg Polytechnic University, Institute of physics, nanotechnology and telecommunications, Russian Federation

**Vladimir Chaldyshev**, Peter the Great St. Petersburg Polytechnic University, , Russian Federation

**Pavel Tonkaev**, Peter the Great St.Petersburg Polytechnic University, Institute of physics, nanotechnology and telecommunications, Russian Federation

**Abstract:**Software was developed for quick numerical calculations and graphic display of the absorption, reflection and transmittance spectra of two-dimensional systems of small conductive particles, which allowed us to make instant comparison with experimental data. A lattice-gas model was used to simulate nearly distributed particles, and the coherent-potential approximation is applied to obtain a solution to the problem of interacting particles. The Delphi programming environment was used.

24

**Pavel Tonkaev**, Peter the Great St.Petersburg Polytechnic University, Institute of physics, nanotechnology and telecommunications, Russian Federation

**Alexey Kondikov**, Peter the Great St. Petersburg Polytechnic University, Institute of physics, nanotechnology and telecommunications, Russian Federation

**Vladimir Chaldyshev**, Peter the Great St. Petersburg Polytechnic University, , Russian Federation

**Abstract:**Measurements of the optical reflection spectra from periodic structures with two quantum wells in the elementary cell have been done. The dependencies of the light reflection on the angle of the light incidence, polarization and temperature were studied. An analysis of the experimental data showed that the pattern with 60 cells is a good Bragg reflector with reflectivity more than 90% in the maximum of the spectral band.

25

**Alessandra Magno**, Universidade Federal de Juiz de Fora, Departamento de Ciência da Computação, Brazil

**João Vitor De Sá Hauck**, Universidade Federal de Juiz de Fora, Departamento de Ciência da Computação, Brazil

**Itamar De Oliveira**, Universidade Federal de Juiz de Fora, Departamento de Ciência da Computação, Brazil

**Abstract:**The main goal of this work is to add and analyse an equation that represents the volume in a dynamical model of the mammalian cell cycle proposed by Gérard and Goldbeter (2011). The cell division occurs when the cyclinB/Cdk1 complex is totally degraded and it reaches a minimum value. At this point, the cell is divided into two newborn daughter cells and each one will contain the half of the cytoplasmic content of the mother cell. The equations of our base model are only valid if the cell volume, where the reactions occur, is constant. Whether the cell volume is not constant, that is, the rate of change of its volume with respect to time is explicitly taken into account in the mathematical model, then the equations of the original model are no longer valid. Therefore, every equations were modified from the mass conservation principle for considering a volume that changes with time. Through this approach, the cell volume affects all model variables. Two different dynamic simulation methods were accomplished: deterministic and stochastic. In the stochastic simulation, the volume affects every model's parameters which have molar unit, whereas in the deterministic one, it is incorporated into the differential equations. In deterministic simulation, the biochemical species may be in concentration units, while in stochastic simulation such species must be converted to number of molecules which are directly proportional to the cell volume. In an effort to understand the influence of the new equation an stability analysis was performed. This elucidates how the growth factor impacts the stability of the model's limit cycles. In order to find the approximated solution of the deterministic model, the fourth order Runge Kutta method was implemented. As for the stochastic model, the Gillespie's Direct Method was used. In conclusion, a more precise model, in comparison to the base model, was created for the cell cycle as it now takes into consideration the cell volume variation.

26

**Julio Espinoza-Ortiz**, Federal University of Goias, Physics Department , Brazil

**F. C. Bauke**, Universidade Estadual de Campinas, Instituto de Física Gleb Wataghin, Brazil

**R. E. Lagos**, Universidade Estadual Paulista, Departamento de Física, IGCE, Brazil

**Abstract:**We consider a charged Brownian particle bounded by an harmonic potential, embedded in a Markovian heat bath and driven from equilibrium by external electric and magnetic fields. We develop a quaternionic-like (or Pauli spinor-like) representation, hitherto exploited in classical Lorentz related dynamics. Within this formalism, in a very straight forward and elegant fashion, we compute the exact solution for the resulting generalized Langevin equation, for the case of a constant magnetic field. For the case the source electromagnetic fields satisfy Maxwell's equations, yielding spinor-like Mathieu equations, we compute the solutions within the JWKB approximation. With the solutions at hand we further compute spatial, velocities and crossed time correlations. In particular we study the (kinetically defined) nonequilbrium temperature. Therefore, we can display the system's time evolution towards equilibrium or towards non equilibrium (steady or not) states.

**Acknowledgements:**JSEO would like to thank the support of the Goiás Research Foundation - FAPEG

27

**Pedro Alberto**, University of Coimbra, Department of Physics, Portugal

**Manuel Malheiro**, Instituto Tecnológico de Aeronáutica, DCTA, 12228-900 São José dos Campos, , Brazil

**Tobias Frederico**, Instituto Tecnológico de Aeronáutica, DCTA, 12228-900 São José dos Campos, , Brazil

**Antonio Castro**, Departamento de Fisica e Quimica, Universidade Estadual Paulista, 12516-410 Guaratinguetá, , Brazil

**Abstract:**A generalization of pseudospin and spin symmetries, the SU(2) symmetries of Dirac equation with scalar and vector mean-field potentials originally found independently in the 70's by Smith and Tassie, and Bell and Ruegg, is proposed. As relativistic symmetries, they have been extensively researched and applied to several physical systems for the last 18 years. The main feature of these symmetries is the suppression of the spin-orbit coupling either in the upper or lower components of the Dirac spinor, thereby turning the respective second-order equations into Schrödinger-like equations, i.e, without a matrix structure. In this paper we use the original formalism of Bell and Ruegg to derive general requirements for the Lorentz structures of potentials in order to have these SU(2) symmetries in the Dirac equation, again allowing for the suppression of the matrix structure of the second-order equation of either the upper or lower components of the Dirac spinor. Furthermore, we derive equivalent conditions for spin and pseudospin symmetries with 2- and 1-dimensional potentials and list some possible candidates for 3, 2, and 1 dimensions. We suggest applications for physical systems in three and two dimensions, namely electrons in graphene.

**Acknowledgements:**Universidade Estadual Paulista, Guaratinguetá Campus, FAPESP thematic project 2015/26258-4 and CNPq.

28

**Ayse Humeyra Bilge**, Kadir Has University, Faculty of Engineering and Natural Sciences, Turkey

**Abstract:**Let y(t) be a sigmoidal function, i.e, a monotone increasing curve with horizontal asymptotes and such that all derivatives vanish at infinity. Let t_n be the point where the n'th derivative of y(t) reach its global extremum. Under fairly general assumptions on the Fourier and Hilbert transforms of y(t), we prove that the sequence {t_n} is convergent and call it "the critical point of the sigmoidal curve". In the context of phase transitions, the limit point is interpreted as a junction point of two different regimes where all derivatives undergo their higest rate of change. The existence of a critical point implies in particular a preferred choice of the origin of the time axis and an intrinsic definition of the even and odd components of a sigmoidal function. The location of the critical point agrees qualitatively with the location of the gel point in the sol-gel transition as discussed in previous work.

29

**Ju Chunhua**, Harbin Institute of Technology, Department of Materials Science, China

**Abstract:**GRP pipes for its excellent performance has been widely used for oil and gas gathering pipe. Because CO2 is used as displacing agent, causing corrosion problems CO2 FRP pipe, so the research GFRP pipe corrosion mechanism and performance in CO2 media evolution is important. Based on material corrosion theory summary FRP pipe used in CO2 environment, we analyze the main form and the specific process material corrosion, the use of aromatic amine curing and acid anhydride curing two fiberglass pipe as a test tube, analog GFRP pipe containing crude oil, the next CO2 and water medium corrosion test, by surface morphology material, the glass transition temperature changes, and performance to evaluate physical corrosion mechanism of FRP in CO2 in; the test results show that CO2 environment, glass and steel addition to the visible Exterior changes, its glass transition temperature, Pap hardness, ring stiffness, density, and curing of the resin content showed a downward trend.

30

**Ju Chunhua**, Harbin Institute of Technology, Department of Materials Science, China

**Abstract:**GRP pipes for its excellent performance has been widely used for oil and gas gathering pipe. Because CO2 is used as displacing agent, causing corrosion problems CO2 FRP pipe, so the research GFRP pipe corrosion mechanism and performance in CO2 media evolution is important. Based on material corrosion theory summary FRP pipe used in CO2 environment, we analyze the main form and the specific process material corrosion, the use of aromatic amine curing and acid anhydride curing two fiberglass pipe as a test tube, analog GFRP pipe containing crude oil, the next CO2 and water medium corrosion test, by surface morphology material, the glass transition temperature changes, and performance to evaluate physical corrosion mechanism of FRP in CO2 in; the test results show that CO2 environment, glass and steel addition to the visible Exterior changes, its glass transition temperature, Pap hardness, ring stiffness, density, and curing of the resin content showed a downward trend.

31

**Ju Chunhua**, Harbin Institute of Technology, Department of Materials Science, China

**Abstract:**GRP pipes for its excellent performance has been widely used for oil and gas gathering pipe. Because CO2 is used as displacing agent, causing corrosion problems CO2 FRP pipe, so the research GFRP pipe corrosion mechanism and performance in CO2 media evolution is important. Based on material corrosion theory summary FRP pipe used in CO2 environment, we analyze the main form and the specific process material corrosion, the use of aromatic amine curing and acid anhydride curing two fiberglass pipe as a test tube, analog GFRP pipe containing crude oil, the next CO2 and water medium corrosion test, by surface morphology material, the glass transition temperature changes, and performance to evaluate physical corrosion mechanism of FRP in CO2 in; the test results show that CO2 environment, glass and steel addition to the visible Exterior changes, its glass transition temperature, Pap hardness, ring stiffness, density, and curing of the resin content showed a downward trend.

32

**Grigorios Megariotis**, National Technical University of Athens, Chemical Engineering, Greece

**Georgios G. Vogiatzis**, National Technical University of Athens, School of Chemical Engineering, Greece

**Ludwig Schneider**, Georg-August-Universität, , Germany

**Marcus Müller**, Georg-August-Universität, , Germany

**Doros N. Theodorou**, National Technical University of Athens, , Greece

**Abstract:**A new methodology and the corresponding C++ code for mesoscopic simulations of elastomers are presented. The test system, crosslinked cis-1,4-polyisoprene, is simulated with a Brownian Dynamics/kinetic Monte Carlo algorithm as an assemblage of beads, each consisting of 5-10 Kuhn segments. Three kinds of beads (nodal points) are considered: end-points, internal beads, and crosslinks. Polymer networks are created from linear melts by either random crosslinking or end-grafting. From a thermodynamic point of view, the system is described by a Helmholtz free-energy containing contributions from entropic springs between successive beads along a chain, slip-springs representing entanglements between beads on different chains, and non-bonded interactions. All beads execute Brownian motion in the high friction limit. The free-energy density of non-bonded interactions is derived from the Sanchez-Lacombe equation of state and the isothermal compressibility of the polymer network is predicted correctly. The ends of each slip-spring hop between adjacent beads along the chains it connects. These thermally activated hops are tracked by a microscopically reversible Kinetic Monte Carlo scheme. Moreover, slip-spring creation/destruction takes place at dangling chain ends in the polymer network. The methodology is employed for the calculation of the stress relaxation function from simulations of several microseconds at equilibrium, as well as for the prediction of stress-strain curves of crosslinked polymer networks under deformation.

**Acknowledgements:**The major part of this work was funded by Volkswagen Foundation in the context of the project Mesoscopic Simulations of Viscoelastic Properties of Networks. In addition, it was funded by the European Union through the projects COMPNANOCOMP and COMPNANOCOMP-DPI under grant number 295355 and Limmat Foundation. Grigorios Megariotis and Doros N. Theodorou also thank the Limmat Foundation for giving them the opportunity to continue this research on polymer networks. Georgios G. Vogiatzis thanks the Alexander S. Onassis Public Benefit Foundation for a doctoral scholarship.

33

**William Yamashita**, Universidade Federal de Juiz de Fora, , Brazil

**Lucy Takahashi**, Universidade Federal de Juiz de Fora, Mathematics, Brazil

**Grigori Chapiro**, Universidade Federal de Juiz de Fora, Mathematics, Brazil

**Abstract:**The global incidences of dengue and, more recently, zica virus have increased the interest in studying and understanding the mosquito population dynamics. Understanding this dynamics is important for public health in countries where climatic and environmental conditions are favorable for the propagation of these diseases. This work is based on the study of nonlinear mathematical models dealing with the life cycle of the dengue mosquito using partial differential equations [1,2]. We investigate the existence of traveling wave solutions using semi-analytical method combining dynamical systems techniques and numerical integration [3]. Obtained solutions are validated through direct numerical simulations using finite difference schemes. [1] I. Freire, M. Torrisi, Symmetry methods in mathematical modeling of Aedes aegypti dispersal dynamics, Nonlinear Anal.-Real 14 (3) (2013) 1300 – 1307. [2] L. Takahashi, N. Maidana, W. Ferreira Jr., P. Pulino, H. Yang, Mathematical models for the aedes aegypti dispersal dynamics: travelling waves by wing and wind, B. Math. Biol. 67 (3) (2005) 509–528. [3] W. M. S. Yamashita, L. T. Takahashi, G. Chapiro, Traveling wave solutions for the dispersive models describing population dynamics of Aedes Aegypti. Submitted, 2016

**Acknowledgements:**This work was supported, in part, by FAPEMIG

34

**Pablo Martin**, Universidad de Antofagasta, Physics, Chile

**Jorge Olivares-Funes**, Universidad de Antofagasta, Mathematics, Chile

**Luis Cortes-Vega**, Antofagasta University, Mathematics Department, Chile

**Abstract:**The computation of Bessel functions is very important in several areas of Physics [3], Mathematics [5] and Chemistry [7], though the radius of convergence of the power series is infinite, however the number of terms required to calculate the functions in big values are many. There are several approximations to this functions [2-5], however the accuracy is not enough for some applications or the range of approximation is narrow. Here via a Multipoints Quasirational Approximations (MPQA Technique) [6] which have better accuracy than previous ones and the degrees of the polynomials in the approximation are only second order. Each approximation can be used with good accuracy for any positive value of the variable. The simplest approximation is obtained with one rational function, quotient of two polynomials combined with some auxiliar functions, and the most accurate approximations are obtained using two rational functions.

35

**Ioannis ZOIS**, PPC, TRSC, Greece

**Abstract:**The motivation for this work came from an attempt to give an alternative definition for the meter, the SI unit for measuring length. As a starting point towards this goal we present the underlying theory behind our approach which uses ideas from quantum field theory and non-commutative geometry, in particular the notion of an odd K-cycle which is based on the Dirac operator (and its inverse, the Dirac propagator). Using physics terminology, the key point in our strategy is this: instead of measuring ordinary length in space-time we measure the “algebraic (or spectral) length” in the space of quantum states of some fermion acted upon by the Dirac propagator.

36

**Jorge Linares**, Université de Versailles-St. Quentin en Yvelines, GEMaC, France

**Salim Guerroudj**, Université de Versailles-St. Quentin en Yvelines, Département de Sciences Physique, France

**Kamel Boukheddaden**, Université de Versailles-St. Quentin en Yvelines, GEMaC, France

**Abstract:**There has been a growing in a recent years, of the interest to design spin crossover (SCO) nanoparticles [1], with controlled shape and environment, due to their potential applications as sensors, bistable and reversible memories or nano-actuators. The experimental studies confirmed the role of the surface effects leading to unusual and non-trivial size dependence of their bistable character. In particular in some cases, the SCO nanoparticles may show a cooperative switching accompanied with a thermal hysteresis at small sizes while they exhibit gradual conversions (no phase transitions) at bigger sizes. This goes against common sense of usual thermodynamics. Even thought, in some special cases, re-entrant phase transitions are obtained with size. The mechanisms of these intriguing behaviors remain however still misunderstood and so their modeling becomes necessary. In this contribution we investigate, the above cited effects on SCO materials, by means of an Ising-like Hamiltonian [2, 3] including short and long-range interactions as well as the interaction of the edge molecules of the system with their local environment. The Hamiltonian is solved in the frame of Monte Carlo Metropolis procedure for 2D (N_x×N_y) rectangular-shaped and 3D parallelepiped rectangle (N_x×N_y×N_z) SCO nanoparticles. The thermal dependence of the high-spin fraction is then derived for different particle sizes and the results are discussed with the help of analytical predictions allowing to evaluate the size dependence of the transition temperature and the conditions of occurrence of the re-entrant phase transition. In practice, the latter emerged from the MC simulations below some critical nanoparticle size and for well identified conditions on the model parameter values. References. [1] Thibault C., Molnar G., Salmon L., Bousseksou A., Vieu C., Langmuir 2010, 26, 1557-1560 [2] Bousseksou A., Nasser J., Linares J., Boukheddaden K., Varret F., Journal de Physique I, 1992, 2, 1381-1403 [3] Linares J., Jureschi C., Boukheddaden K., Magnetochemistry, in press

**Acknowledgements:**CHAIR Materials Simulation and Engineering, UVSQ, Université Paris Saclay is gratefully acknowledged.

37

**Sabit Bekov**, Eurasian National University, General and Theoretical Physics, Kazakhstan

**Abstract:**In this paper, we consider a new model with a fermionic field that is non-minimally coupled to gravity in the framework of F(T) gravity for the Friedman–Robertson–Walker metric. In order to determine the forms of the coupling and potential function of fermionic field for the considered model, we use the Hojman symmetry approach. Furthermore, we search for the cosmological solution of the model. It is shown that the fermionic field can explain the modern evolution of the universe.

38

**Orestis George Ziogos**, National Technical University of Athens, School of Chemical Engineering, Greece

**Grigorios Megariotis**, National Technical University of Athens, Chemical Engineering, Greece

**Doros N. Theodorou**, National Technical University of Athens, , Greece

**Abstract:**This study concerns atomistic and coarse-grained molecular dynamics simulations of pristine Hexabenzocoronene (HBC) molecular crystals. HBC falls in the category of polyaromatic hydrocarbons that find numerous applications in the field of organic electronics. The HBC molecule is initially simulated in an all-atom representation by employing the well-known AMBER force field. The atomistic model reproduces fairly well the structural experimental properties of this system and thus can be used as a reliable starting point for the development of a coarse-grained model in a bottom-up coarse-graining approach. All simulations were conducted at a temperature of 300K, where the system is found in the crystalline phase. In particular, the system consists of molecules that self-organize into well aligned molecular stacks which in turn create a perfect monoclinic molecular crystal. The coarse-grained model is developed by applying Iterative Boltzmann Inversion, a systematic coarse-graining method which reproduces a set of target (atomistic) radial distribution functions and intramolecular distributions at the coarser level of description. This model allows the simulation of HBC crystals over longer time and length scales. In particular, upon coarse-graining, the integration time step of the equations of motion increases from 1fs to 12fs, while the number of interaction sites per molecule is reduced from 60 to 6. The polyaromatic core of HBC at the coarse-grained level is replaced by six superatoms of the same type lying in the same plane. The crystalline phase is analyzed in terms of the Saupe tensor and its eigenvalues.

**Acknowledgements:**This research has been co-financed by the European Union (European Social Fund – ESF) and Greek national funds through the Operational Programme ‘Education and Lifelong Learning’ of the National Strategic Reference Framework (NSRF) – Research Funding Programme: THALIS. Investing in knowledge society through the European social fund [grant number MIS 379436] and by the IKY Fellowships of Excellence for Postgraduate Studies in Greece – Siemens Programme [award number 11126/13a]. Grigorios Megariotis and Doros N. Theodorou also thank the Limmat Foundation.

39

**Petr Bour**, Faculty of Nuclear Sciences and Physical Engineering of the Czech Technical University in Prague, Department of Mathematics , Czech Republic

**Vaclav Kus**, Faculty of Nuclear Sciences and Physical Engineering of the Czech Technical University in Prague, Department of Mathematics , Czech Republic

**Abstract:**Binary decision trees are widely used tool for unsupervised classification of high-dimensional data, for example among particle physicists. We present our proposal of the supervised binary divergence decision tree with nested separation method based on the generalized linear models or on the kernel density estimation procedure. A key insight we provide is the clustering driven only by a few selected physical variables. The proper selection consists of the variables achieving the maximal divergence measure between two different classes. Further we apply our method to Monte Carlo data set simulating measured observations from the ttbar lepton+jets decay channel at DZERO experiment in Fermilab. The efficiency of our algorithm achieves 90 % AUC. We also introduce the modification of statistical tests applicable to weighted data sets in order to test homogeneity of the Monte Carlo simulation and real data. We provide a justification of powers of these modified tests by presenting their p-values for weighted Monte Carlo samples .

**Acknowledgements:**This work was supported by the grant SGS15/214/OHK4/3T/14 and INGO II - LG15047.

40

**Pantea Davoudifar**, Research Institute for Astronomy and Astrophysics of Maragha , Astroparticle Physics , Iran (Islamic Republic of)

**Keihanak Rowshan Tabari**, Research Institute for Astronomy and Astrophysics of Maragha (RIAAM), Astroparticle physics Group, Iran (Islamic Republic of)

**Abstract:**A model of turbulent galactic magnetic fields was developed, in which the type of magnetic fields were considered to be Kolmogorov. We tested the effect of this model on an isotropically distributed flux of ultra high energy cosmic ray in the extragalactic space. To do this, a giant halo of Galaxy was considered. Regular and random components of the Galactic Magnetic Fields were considered to have the mean observed relevant values and also satisfy a Kolmogorov field type. The deviation from isotropy then were calculated using this model and the results were discussed to show how isotropic is the flux of ultra high energy cosmic rays in the extragalactic space.

**Acknowledgements:**This work has been supported financially by research Institute for Astronomy and Astrophysics of Maragha (RIAAM)-Maragha, under research project No. 1/4165-23

41

**Pantea Davoudifar**, Research Institute for Astronomy and Astrophysics of Maragha , Astroparticle Physics , Iran (Islamic Republic of)

**Keihanak Rowshan Tabari**, Research Institute for Astronomy and Astrophysics of Maragha (RIAAM), Astroparticle physics Group, Iran (Islamic Republic of)

**Abstract:**Using a Galactic magnetic field model in combination with an extragalactic model of magnetic fields and considering the observed flux of ultra high energy cosmic rays; an anisotropic flux of these particles were estimated in the intergalactic space. To do this, we have considered the ultra high energy events were all originated from extragalactic sources in the distances not more than 75 mega parsecs. Then using our results, the origin of ultra high energy cosmic rays in the energies ranged from 10^18 up to the highest energies were discussed for the possible extragalactic point sources.

**Acknowledgements:**This work has been supported financially by research Institute for Astronomy and Astrophysics of Maragha (RIAAM)-Maragha, under research project No. 1/4165-27

42

**Ayse Humeyra Bilge**, Kadir Has University, Faculty of Engineering and Natural Sciences, Turkey

**Abstract:**The Istanbul Straight (Bosphorus) is a 31 km long narrow channel in the north-south direction joining the Black Sea and the Sea of Marmara basins. The high salinity Mediterranean water flows at the bottom layer from south to north while the low salinity Black Sea water flows at the top layer in the reverse direction. Increased river fluxes to the Black Sea during the spring and strong southward/northward wind systems influence and disturb the stratified flow in the Istanbul straight, causing blockages and current reversals under extreme conditions. The south and north exits of the straight are marked by sills at depths 33 and 61 meter respectively, forming hydrological control points. There is a third control point 10 km north of the southern entrance, caused by a contraction. Salinity profiles show that the flow along this 10 km part of the straight is a nearly horizontal stratified flow [Ozsoy, Iorio, Gregg, Backhaus, J. Marine Systems 31 (2001) 99-135]. There is about 30 cm difference between the sea levels at the north and south entrance, causing a net southward flux throughout the year. Excessive river runoffs to the Black Sea, that are more effective in summer, cause variations in the fluxes of both layers but the effects of atmospheric forcing on flux variations is less accentuated [Jarozs, Teage, Book, Besiktepe, Geophys. Res. Lett. 38 (2011) L21608]. In this work we study the variations in the upper layer due to wind effects under the assumption of uniform flow with constant flux, as a toy model at least for the region between the two south control points of the Bosphorus. Such effects have been analysed in an experimental study [Plate and Goodwin, Proceedings Santa Barbara Speciality Conference (1965)] where it is shown that sufficiently prolonged strong winds lead to increased velocities and decreased height in open channel flows. We use an analytical model to describe the steady-state water velocity profile induced by wind effects and use the constant flux assumption to compute the steady state layer height by determining the boundaries of the flux integral. The results agree with the available data for the water and wind speeds at observed blocakge and current reversal events.

43

**Andrei Lebed**, University of Arizona, Department of Physics, United States

**Abstract:**We determine active gravitational mass operator of the simplest composite quantum body - a hydrogen atom - within the semiclassical approach to the Einstein equation for a gravitational field. We show that the expectation value of the mass is equivalent to energy for stationary quantum states, as expected. On the other hand, it occurs that, for superpositions of stationary quantum states with constant expectation values of energy, the expectation values of the gravitational mass exhibit time-dependent oscillations. This breaks the equivalence between active gravitational mass and energy and can be observed as a macroscopic effect for a macroscopic ensemble of coherent quantum states of the atoms. We discuss experimental aspects of the above-mentioned breakdown and stress that the corresponding experiment could be the first direct observation of quantum effects in General Relativity.

**Acknowledgements:**This work was partially supported by the NSF under Grant DMR-1104512.

44

**Hiroaki Nakamura**, National Institute for Fusion Science, Department of Helical Plasma Reseach, Japan

**Abstract:**It is one of unsolved issues to find the plasma facing material (PFM) which endure heat flow from core plasma in order to realize nuclear fusion. Tungsten (W) is a good candidate material for PFM. Therefore, a lot of experimental researchers in plasma-wall-interaction research field are gathering reams of physical information of tungsten under irradiation of noble gas or hydrogen. In addition to the experimental measurement, numerical simulation on an atomic scale is also performed to elucidate a behavior of W. In this paper, we pick up nanostructured W, which is generated around the W surface when helium (He) is irradiated on an appropriate condition. For example, nano-filament structure, called fuzz, and bubble structure were discovered by the experiment. From a view point of the plasma confinement, it is necessary to evaluate the endurance of the nanostructure W against the plasma irradiation. To achieve the request, we gather the sputtering yield of W using binary-collision-approximation (BCA) simulation. Then we estimate the endurance of W which depends on the thickness of nanostructure layer.

**Acknowledgements:**This work was supported by JSPS KAKENHI Grant Numbers (15K06650, and 26330237), MEXT KAKENHI Grant Number (25249132 and 15H04142), and the NIFS Collaboration Research programs (NIFS14KNTS028 and NIFS14KBAS016).

45

**Vladimir Kashurnikov**, National Research Nuclear University MEPhI, , Russian Federation

**Andrey Krasavin**, National Research Nuclear Univewrsity MEPhI, Physical and Technical Problems of Metrology, Russian Federation

**Yaroslav Zhumagulov**, National Research Nuclear University MEPhI, , Russian Federation

**Abstract:**The momentum distribution of charge carriers in two-dimentional FeAs-clusters modeling iron-based superconductors was calculated within the limits of the two-orbital model. The calculation was implemented by use of quantum world-line Monte Carlo algorithm. Matsubara Green's function was calculated for clusters with sizes up to 10x10 FeAs-cells. Profiles of the distributions are presented for the entire Brillouin zone. The data indicate the presence of a jump near the Fermi level, which is the evidence of the Fermi-type distribution. The influence on the momentum distribution of the value of the interaction parameter, cluster size and temperature is discussed.

46

**Chehung Wei**, Tatung University, Mechanical Engineering, Taiwan

**Li-Chung Chang**, Tatung Univerisity, Mechanical Engineering, Taiwan, Province Of China

**Abstract:**The effect of substrate roughness on the film hardness is often modeled by a periodic roughness distribution to simplify the analysis. However, in real world application the roughness is random rather than periodic. Therefore, it is essential to assess the discrepancy between simplified and real world modeling. To address this problem, the hardness is evaluated by a nanoindentation finite element analysis using diamond-like carbon (DLC) film on a rough silicon substrate. Arithmetic average roughness Ra is chosen to represent the substrate roughness. Two types of roughness distributions, periodic and random distributions are used to represent the simplified and real world problem, respectively. The comparison is based on fact that the average roughness magnitude in random distribution is the same as that in periodic distribution. Six roughness magnitudes 5 nm, 20 nm, 35 nm, 50 nm, 200 nm, 350 nm and 500 nm are utilized to investigate the effect of roughness distribution. The results show the hardness is affected by the substrate roughness variation. For periodic roughness distribution, the hardness exhibits a synchronized pattern irrespective of the roughness magnitude. For random roughness, the randomness in hardness variation is not significant when Ra is less than 50 nm. When Ra is above 200 nm, the hardness variation becomes larger and exhibits no specific pattern. The local hardness in low Ra might be larger than that in high Ra. This indicates for random roughness, Ra alone can not predicate the overall hardness behavior. Another finding is, the average hardness in periodic roughness is larger than that in random roughness for same average Ra. This implies roughness modeling by periodic roughness tends to overestimate the hardness. In conclusion, the type of roughness distribution has an effect in local and overall hardness. The local roughness might affect the hardness measurement.

47

**Irina Knyazeva**, Pulkovo observatory, sector of mathematical modeling in geophysics, Russian Federation

**Irina Makarenko**, Newcastle University, UK, , United Kingdom

**Alina Vdovina**, Newcastle University, UK, , United Kingdom

**Nikolay Makarenko**, Pulkovo observatory, , Russian Federation

**Fedor Urtiev**, Pulkovo observatory, , Russian Federation

**Abstract:**Analysing dynamics of the photospheric magnetic field of the Sun is one of the most important problems in Solar Physics. Different estimates of the complexity of magnetograms of the Sun Active Regions (AR) are used to predict the time and the strength of the solar flares. Complexity of random fields (including magnetic fields) is described by the structure and behavior of their critical points. A magnetogram is a highly variable discrete image with very large number of local extrema. We use an idea of extraction of stable critical points within a framework of the scale-space proposed in [Steve Gu et al, 2010]. The convolution of the image with a Gaussian kernel and the difference between the received images allow to get a stable estimation of the Laplacian. A critical graph is constructed using maxima and minima of the Laplacian. This graph is not a Morse graph since it does not have saddle points. In the process of evolution of an AR the number and the values of local extrema are changing. Therefore the structure of a critical network is also changing. So the dynamics of critical graphs can be used for diagnosis of dynamical regimes of ARs. As a numerical descriptor we propose to use the so-called spectral gap. It is the difference between the two largest eigenvalues of the discrete Laplacian of the graph constructed for critical networks. Here we provide results for several ARs.

48

**Ioannis Mathioudakis**, National Technical University of Athens, School of Chemical Engineering, Greece

**Georgios G. Vogiatzis**, National Technical University of Athens, School of Chemical Engineering, Greece

**Christos Tzoumanekas**, National Technical University of Athens, School of Chemical Engineering, Greece

**Doros N. Theodorou**, National Technical University of Athens, , Greece

**Abstract:**The interaction energetics, molecular packing, entanglement network properties, segmental dynamics, and elastic constants of atactic polystyrene – amorphous silica nanocomposites, in molten and glassy state are studied via molecular simulations using two interconnected levels of representation: (a) A coarse – grained level, wherein each polystyrene repeat unit is mapped onto a single “superatom” and the silica nanoparticle is viewed as a single sphere. Equilibration at all length scales at this level is achieved via connectivity – altering Monte Carlo simulations. (b) An atomistic level, where the polymer chains are represented in terms of a united – atom forcefield and the silica nanoparticle is represented in terms of a simplified, fully atomistic model. Initial configurations for atomistic Molecular Dynamics (MD) simulations are obtained by reverse mapping well-equilibrated coarse – grained configurations. By analyzing microcanonical atomistic MD trajectories, the polymer density profile is found to exhibit layering in the vicinity of the nanoparticle surface. An estimate of the enthalpy of mixing between polymer and nanoparticles, derived from the atomistic simulations, compares favorably with available experimental values. The dynamics of polystyrene (in neat and filled melt systems) is characterized in terms of bond orientation and dihedral angle time autocorrelation functions. At low concentration in the molten polymer matrix, silica nanoparticles are found to bring about a slight deceleration of the segmental dynamics close to their surface relative to the bulk polymer. Well – equilibrated coarse – grained long – chain configurations are reduced to entanglement networks via topological analysis with the CreTA algorithm, yielding a slightly lower density of entanglements in the filled than in the neat systems. Atomistic melt configurations are glassified by MD cooling. The Young’ and bulk moduli of the resulting glassy nanocomposites are computed through an analysis of stress fluctuations in the undeformed state and through explicit mechanical deformation by MD, showing a stiffening of the polymer in the presence of nanoparticles. Atomistic simulation results for the elastic constants are compared to the predictions of continuum micromechanical calculations invoked in homogenization models of the overall mechanical behavior of heterogeneous materials.

**Acknowledgements:**Part of this work was funded by the European Union through the project NanoModel under the grant number 211778.

49

**Natalya Shipulya**, Peter the Great St.Petersburg Polytechnic University, , Russian Federation

**Stepan Konakov**, Peter the Great St.Petersburg Polytechnic University, , Russian Federation

**Valeria Krzhizhanovskaya**, University of Amsterdam, The Netherlands; ITMO University and St. Petersburg Polytechnic University, , Russian Federation

**Abstract:**In this work we present the results of analytical modeling and 3D computer simulation of microfluidic Wheatstone bridge, which is used for high accuracy measurements or precision systems. A geometrical configuration of microchannels was selected based on analytical estimations. The detailed 3D model is based on Navier-Stokes equations for a laminar flow of liquid in microchannels with low Reynolds number. We proposed and simulated a new method of a bridge balancing process by changing the microchannel geometry. This process is based on the “etching in microchannel” technology we developed earlier. Our method implements a precise control of flow rate and flow direction in the ‘‘bridge’’ microchannel. An advantage of our approach is ability to work without any control valves and other active electronic systems, which are usually used for bridge balancing. Using an analytical model and results of 3D computer simulation we investigated the behavior of the Wheatstone bridge with different model parameters. We found a relation between the channel resistance and flow rate through the bridge, and calculated a pressure drop across the system under different total flow rates and viscosities. Some ideas about application of described Wheatstone bridge in accuracy measurements and complex precision microfluidic systems were suggested for future work.

50

**Pierre Dahoo**, Université de Versailles-St. Quentin en Yvelines, LATMOS, France

**Raphael Puig**, CHAIR Materials Simulation and Engineering UVSQ, Université Paris Saclay, 78035 Versailles Cedex, France, , France

**Azzedine Azzedine Lakhlifi**, Institut UTINAM-UMR 6213 CNRS- Université de Franche-Comté, Observatoire de Besançon, 25010 Besançon Cedex, France, , France

**Constantin Meis**, INSTN- CEA-Saclay – 91191 Gif-sur-Yvette Cedex, France, , France

**Abstract:**: Clathrates are characterized by a compact assembly of low temperature nano-cages that naturally contain an atom or molecule that stabilizes their formation as cubic sI and sII or hexagonal sH structures. Infrared vibrational frequencies of triatomic carbon dioxide CO2 trapped in the nano cages are shifted (less than 0.5% of the observed gas phase frequencies as a result of a slight perturbation of the intramolecular electronic potential of the CO2 trapped molecule) from those of gas phase and the spectra are simpler because the rotational and translational degrees of freedom are cage-hindered. The energy levels (eigen states of the Hamiltonian operator) of the different degrees of freedom of CO2 are calculated within the frame of the Born Oppenheimer approximation [1][2]. The potential energy hypersurface experienced by the molecule in its nanocage of type n can then be written in terms of the minimum of the potential energy (the equilibrium configuration); the vibrational dependence in the normal vibrational coordinates {Q}; the orientational dependence in Euler angles(Omega=(phi, theta, khi),khi=0 for a linear molecule);the vibration-orientation coupling term, which induce the relaxation channels and n characterize the site (n=1 for small and n=2 for or large cage). The the displacement vectors (centre of mass) of of water molecules of the clathrate and the trapped molecule are calculated from the Green functions of the perfect crystal. In an attempt to simulate the relaxation of absorbed infrared energy in a clathrate, the different channels through which energy can spread inside, after CO2 absorbs infrared radiation are depicted. The results of these calculations are to be connected to the effects of global warming on the stability of clathrates. Depending on the competition between radiative and non-radiative relaxation, a limited temperature rise warming effect of the clathrate may occur as in low temperature matrices because of stimulated emission. This depends on how isolated CO2 is inside the nanocages of a clathrate. Results are compared from the point of view of the trapping site effect, the shifts observed for the different fundamental vibrational modes and the effect of coupling on the relaxation constants of intramolecular energy levels. [1] A. Lakhlifi, P.R. Dahoo, E. Dartois, E. Chassefière et O. Mousis, Modeling IR spectra of CO2 isotopologues and trapped In type I clathrate EPOV 2012: From Planets to Life - Colloquium of the CNRS Interdisciplinary Initiative "Planetary Environments and Origins of Life", Paris, France, November 29-30, 2012, février, vol 2, 03008, EDP Sciences, BIO Web of Conferences, M. Ollivier et M.-C. Maurel, eds., ISBN: 978-2-7598-1180-9, [2]P.R Dahoo., A. Lakhlifi, S. Picaud, C. Thomas, E. Dartois and E. Chassefière, Spectroscopy of CO2 modelled in Clathrates, 56th ICASS, U. Alberta, Edmonton, Canada, 15-18 August 2010.

**Acknowledgements:**CHAIR MSE UVSQ

51

**Pierre Dahoo**, Université de Versailles-St. Quentin en Yvelines, LATMOS, France

**Jorge Linares**, Université de Versailles-St. Quentin en Yvelines, GEMaC, France

**Christian Chong**, LISV, Université de Versailles St. Quentin en Yvelines, 78035 Versailles, France, , France

**Philippe Pougnet**, VALEO, 14 avenue des Béguines - BP 68532 Cergy, 95892 Cergy Pontoise Cedex, France, , France

**Constantin Meis**, INSTN- CEA-Saclay – 91191 Gif-sur-Yvette Cedex, France, , France

**Abdelkhalak El Hami**, INSA ROUEN Laboratoire de Mécanique, Pôle Technologique du Madrillet Avenue de l'Université, BP 8-76801 Saint Etienne du Rouvray, , France

**Abstract:**A Monte Carlo (MC) simulation of a 2D microscopic Ising model is developed to study the effect of thermal cycling on the lead free interconnect alloy of a mechatronics device. In order to understand the physical degradation processes of this alloy leading to failure, the impact of of temperature cycling on the growth and crack propagation [1] is studied. This work focuses on specific mechatronics systems in which die attach reliability needs to be improved. To determine the physical processes leading to failure, an ABV model [2] (Metals A and B and void V) is developed. A Monte Carlo (MC) simulation of this 2D microscopic Ising model is used to study the evolution of a homogeneous structure of a binary alloy containing voids at a particular temperature. In this approach a three states Ising Hamiltonian is used which depends on 3 parameters K,J and U and a fictitious spin σ value which can be set to -1, 0, 1. K, J and U are parameters function of the interaction energy parameters between the metal atoms A, B and void V. In order to model the observed diffusion of intermetallics in interconnexion materials, the interaction parameters EAA, EBB and EAB are adjusted and the time evolution of this material pseudo-spin configuration is calculated. The results are discussed in terms of intensity of these parameters and void percentage. 1: Pierre Richard Dahoo, Nadim Alayli, Armelle Girard, Philippe Pougnet, Ky-Lim Tan and Jean-Michel Morelle, Reliability in Mechatronic Systems from TEM, SEM and SE Material Analysis, in Reliability and Materials Issues of Semiconductor Optical and Electrical Devices and Materials, edited by O. Ueda, M. Fukuda, S. Pearton, E. Piner, P. Montanegro, Mater. Res. Soc. Symp. Proc. Warrendale, PA, 1195, B07-04.R1 (2010). 2: Yaldram, K. and Binder, K., Monte-Carlo simulation of phase-separation and clustering in the ABVmodel. Journal of Statistical Physics, 62, 161, 1991

**Acknowledgements:**CHAIR MSE UVSQ

52

**Humberto Híjar**, La Salle University Mexico, Engineering School, Mexico

**Abstract:**We propose a model for an active Brownian system that exhibits one-dimensional directed motion. This system consists of two Brownian spherical particles that interact through an elastic potential and have time-dependent radii. We suggest an algorithm by which the sizes of the particles can be varied, such that the center of mass of the system is able to move at an average constant speed in one direction. The dynamics of the system is studied theoretically using a Langevin model, as well as from Brownian Dynamics simulations. We also present an irreversible thermodynamics approach for the study of this system.

**Acknowledgements:**Authors acknowledges La Salle University Mexico for support under grant NEC-04/15

53

**Alexander Zuevsky**, Academy of Science of the Czech Republic, Mathematics, Czech Republic

**Abstract:**We give a conformal block description for the Zhu reduction formulae for vertex operator algebras with formal parameters associated to local coordinate on a genus two Riemann surface formed by sewing two initial tori in a geometric procedure. An example of Heisenberg vertex operator algebra is given.

54

**Wutthinan Boonsen**, Prince of Songkla University, Physics, Thailand

**Chutintorn Punwong**, Prince of Songkla University, Physics, Thailand

**Abstract:**Photo-induced electron transfer (PET) reaction has been observed and widely studied in many biological systems such as in proteins and oligonucleotides. PET between fluorescent dyes and amino acids or DNA bases can lead to fluorescence quenching of the dyes. Such quenching reaction can affect the fluorescence properties especially when the dye is employed to probe the conformation of the biomolecular systems. The most effective quencher among nucleotides is guanine, the easiest oxidizable base. 7-methoxycoumarin (C-3H) is one of the coumarin derivatives, widely used as a fluorescent dye. It can be reduced by guanine in the PET process. In this work, the molecular simulation is employed to elucidate the underlying mechanism of the C-3H fluorescence quenching induced by PET. Quantum mechanical method is used because PET involves electron transfer from a DNA quencher to a fluorescent dye in its excited state. The modified semiempirical method is selected for faster calculations especially in the dynamics simulations. The standard semiempirical method is incorporated with the floating occupation molecular orbital-configuration interactions approach to describe the multi-reference nature of an electron in the PET reaction. The dynamics simulations on the first and second singlet excited states (S1 and S2) of the dye-guanine cluster have been performed. The nuclear quantum dynamical characteristics of the cluster is described by the full multiple spawning (FMS) method, which solves the electronic and nuclear Schrödinger equations simultaneously. The obtained results indicate that the PET reaction pathway likely starts from S2 of the cluster and pass through its transient S1 state to the ground state. The similar pathway is also found in other coumarin dye (Barmen, et al. J. Phys. Chem. A 2013, 117, 3945-3953), which also emphasizes the role of hydrogen bonding. From our simulations, we likewise found that the formation of hydrogen bond between C-3H and guanine is necessary for the PET to occur.

**Acknowledgements:**Prince of Songkla University research grant number SCI560506S

55

**Ji Hwan Cha**, Ewha Womans University, Department of Statistics, Korea, Republic Of

**Abstract:**A cascading failure is a failure in a system of interconnected parts, in which the breakdown of one element can lead to the subsequent collapse of the others. Cascading failures usually begin when one part of the system fails. When this happens, remaining components must then take up the slack for the failed component. This in turn overloads these components, causing them to fail as well, prompting additional components to fail one after another. In this paper, we discuss the lifetimes of systems composed of components having interdependence and cascading effect. For this, a class of bivariate distributions is constructed based on the notion of conditional failure rate. We show that the members in a suitably defined subclass of the developed bivariate distributions satisfy the bivariate lack of memory property. Based on the developed class of distributions, we study the lifetimes of systems having interdependence and cascading effect. We compare the lifetimes of systems having interdependence and cascading effect with those of systems having independent components.

56

**Vladimir Kahsurnikov**, National Research Nuclear University (Moscow Engineering Physics Institute), Physics of the Solid State and Nanosystems, Russian Federation

**Anastasiya Maksimova**, National Research Nuclear University (Moscow Engineering Physics Institute), Physics of the Solid State and Nanosystems, Russian Federation

**Igor Rudnev**, National Research Nuclear University (Moscow Engineering Physics Institute), Physics of the Solid State and Nanosystems, Russian Federation

**Dmitry Odintsov**, National Research Nuclear University (Moscow Engineering Physics Institute), Physics of the Solid State and Nanosystems, Russian Federation

**Abstract:**The Monte Carlo simulations have been used to study the dynamics of three-dimensional vortex system in layered high-temperature superconductors with internal nanosized ferromagnetic defects. The magnetization processes in the array of ferromagnetic defects under the field of Abrikosov vortices have been taken into account. The current-voltage characteristics under applying of external dc magnetic field have been obtained. It has been shown that S-type nonlinearity of E(J) curve which had previously been obtained for two-dimensional vortex system also appears in more realistic anisotropic three-dimensional case. The effect of vortex-vortex inter-layer coupling on the S-nonlinearity has been analyzed. The wave of defects magnetization in three-dimensional case arising for the vortex flow has been shown.

**Acknowledgements:**Research was done with the financial support of RSF under grant № 14-22-00098

57

**Oleg Shkolik**, Ural Federal University, Graduate School of Economics and Management, Russian Federation

**Larisa Chirkova**, Ural Federal University, Institute of Physics and Technology, Russian Federation

**Polina Chirkova**, Ural Federal University, Graduate School of Economics and Management, Russian Federation

**Abstract:**Developing (underdeveloped) countries are territories of slow economic growth (catch-up growth). Perspectives of their economic growth largely depend on developing and introducing financial and technological innovations in the sphere of the financial markets. The level and quality of those innovations should enable provision of faster growth of the financial sector of the national economy by rising stability and effectiveness of the financial institutions. Powerful and stable financial sector is the basic element for attracting investments and upsurge of liquidity in the economic system of a developing country that aims to have developed economy. Intellectual capital is the most important of the fundamental factors of production in the financial sphere. It is a catalytic element of the process of the economic development. From this position, the researchers’ collective develops and presents a mathematical model which characterizes the connection between the intellectual capital and financial results of the commercial activity of financial institutions. The model is applied in the analysis of the activity of financial institutions that are part of the EEU.

58

**Eugen Anitas**, Joint Institute for Nuclear Research, Bogoliubov Laboratory of Theoretical Physics, Russian Federation

**Abstract:**The small-angle scattering (SAS) structure factor from a new model of a 3D deterministic fractal in which the relative positions and the number of structural units vary with iteration number is calculated. It is shown that, depending on the relative positions of scattering units inside the fractal, we can obtain various types of power-law successions, such as: “convex/concave” - when the absolute value of the scattering exponent of the first power-law decay is higher/smaller than that of the subsequent power-law decay, or any combination of them (i.e. convex-concave or concave-convex). The obtained results can explain experimental SAS (neutron or X-rays) data which are characterized by a succession of power-law decays of arbitrary length.

59

**Arman Kussainov**, al-Farabi Kazakh National University, Physics and Technology, Kazakhstan

**Natalya Pya**, Nazarbayev University, Department of Mathematics, School of Science and Technology , Kazakhstan

**Abstract:**We have assessed the potential application of the neutron monitor hardware as a random number generator for the normal distribution. The data from acquisition channels with no apparent changes in the signal level, besides the diurnal and other variations typical for this type of data, were chosen as the retrospective model. Variational Monte Carlo method for the quantum harmonic oscillator was used to test the targeted quality of our data. The data delivery rate is of importance and if the conventional minute based neutron counts' averages are insufficient we could always settle for an efficient seed generator to feed into the faster algorithmic random number generator. Scaling our data to the zero mean and variance of one was sufficient to obtain a stable standard normal random variate. Distributions under consideration are passable under all available normality tests. Inverse transform sampling is suggested to use as a source of the uniform random numbers.

**Acknowledgements:**This research was supported by grant №3824/ГФ4 provided by the Science Committee at the Ministry of Science and Education of Republic of Kazakhstan to the principal investigator at the National Nanotechnology Laboratory of Open Type, Physics and Technology Department, al-Farabi Kazakh National University.

60

**Georgiy Nujdov**, Institute of ionic, plazma and laser technology, laboratory of physics of multiphase systems, Uzbekistan

**Valeriy Krivorotov**, Institute of ionic, plazma and laser technology, laboratory of physics of multiphase systems, Uzbekistan

**Sirojiddin Mirzaev**, Institute of ionic, plazma and laser technology, laboratory of physics of multiphase systems, Uzbekistan

**Abstract:**The article presents the results of quantum-chemical calculations of the potential energy of the lattice during a motion of fluoride ions in different areas of the lattice of the rare earth trifluoride LnF3 both dielectric and superionic phases. The profiles of potential energy of the lattice during such a movement are presented and compared. Also shown, that the efficiency of ion transport in the surface layer of the crystal lattice is significantly higher than in the bulk.

61

**Marcin Schroeder**, Akita International University, , Japan

**Abstract:**Neither Ralph Landauer’s declaration that “Information Is Physical”, nor John A. Wheeler’s slogan “It from Bit” or extensive research that followed made information authentically physical. Similarly, attempts to naturalize computation did not make it a natural process. Information and computation became associated with physical phenomena and processes, but nothing would be lost from our understanding of physics, if we drop the terms “information” or "computation" and stay within standard conceptual framework of the discipline using orthodox terms, such as “entropy”. In the case of information, the reason was the lack of a clearly defined concept of information with ontological status of entity independent from underlying physical entities. In the case of computation, it was the lack of autonomy of computation from some forms of human intervention and the lack of recognition that natural processes occur through interactions, not goal oriented actions. In the earlier publications the author provided a conceptual framework for information and computation free from these deficiencies together with an appropriate mathematical formalism and a slightly generalized version of a Turing machine based on interaction (s-machine, which becomes Turing’s a-machine when we consider a special case of one-way actions). This generalization makes it possible to consider evolution of the machine in terms of dynamical changes of both the tape and the head through mutual interaction. The present paper presents more detailed study of the dynamics of information involved in interactions of information systems together with issues related to the invariance with respect to transformations (symmetry). The latter is of special importance, as it is only invariance with respect to reference frame that can give the status of an (objective) physical entity.

62

**Halina Grushevskaya**, Belarusian State University, Physics Department, Belarus

**Nina Krylova**, Belarusian State University, Physics Department, Belarus

**Ihor Lipnevich**, Belarusian State University, Physics Department, Belarus

**Abstract:**In this paper electrocapillary effects at an interphase boundary are subject to modeling by means of Finsler-Lagrange structures which have been proposed in [1-4]. We have shown that such effects are the dependency of system behavior on velocity and the non-conservation of particle mass. A developed simulation technique has been applied to describe the formation of the Langmuir--Blodgett monomolecular layer (monolayer) of amphiphilic molecules on the air / water subphase interface. Experimental data for compressibility of the monolayer in first-order phase transition state are well-approximated by obtained results of the mathematical modelling. [1] V. Balan, H.V. Grushevskaya, N.G. Krylova, M. Neagu, A. Oana. On the Berwald-Lagrange scalar curvature in the structuring process of the LB-monolayer. Applied Sciences. Vol. 15, 30-42 (2013). [2] H.V. Grushevskaya, N.G. Krylova. Effects of Finsler geometry in physics surface phenomena: case of monolayer systems. Hypercomplex Numbers in Geometry and Physics. Vol. 8, 128-146 (2011). [3] M. Neagu, N.G. Krylova, H.V. Grushevskaya. Jet theoretical Yang-Mills energy in the geometric dynamics of 2D-monolayer. J. Math. Phys. Vol. 54, 031508 (2013). [4] V. Balan, H. Grushevskaya, N. Krylova. Finsler geometry approach to thermodynamics of first order phase transitions in monolayers. Differential Geometry - Dynamical Systems. Vol. 17, 24-31 (2015).

63

**Tetsuya Takaishi**, Hiroshima University of Economics, Faculty of Economics, Japan

**Abstract:**We study the stock market instability by cross correlations constructed from the return time series of 366 stocks traded on the Tokyo Stock Exchange from January 5, 1998 to December 30, 2013. To investigate the dynamical evolution of the cross correlations the correlation matrices are calculated with a rolling window of 400 days. To quantify the volatile market stages where the potential risk is high we apply the principal components analysis and measure the cumulative risk fraction (CRF) which is the system variance associated with the first few principal components. From the CRF we detected three volatile market stages corresponding to the bankruptcy of Lehman Brothers, Tohoku Region Pacific Coast Earthquake and FRB QE3 reduction observation in the period we studied. We further apply the random matrix theory for the risk analysis and found that the first eigenvector is more equally delocalized when the market is volatile.

64

**Nikolaos Dimakis**, Universidad Austral de Chile, Instituto de Ciencias Fisicas y Matematicas, Chile

**Abstract:**In reparametrization invariant systems, such as mini-superspace Lagrangians, the existence of constraints can lead to the emergence of additional non-local integrals of motion in phase space. In the case of an FLRW flat/non-flat space-time minimally coupled to an arbitrary scalar field, we manage to use such conserved quantities to completely integrate the system of equations of motion. This is achieved without constraining the potential in any way. Thus, obtaining the most general solution that encompasses all possible cosmological scenarios which can be based on the existence of a scalar field.

**Acknowledgements:**The author acknowledges financial support by FONDECYT postdoctoral grant no. 3150016.

65

**Ana Lúcia Baccon**, Universidade Estadual de Ponta Grossa, Departamento de Matemática de Estatística, Brazil

**Jose Lunardi**, University of Glasgow, School of Physics and Astronomy, United Kingdom

**Abstract:**The characterization of the statistical distributions of observed weather data is of crucial importance for both constructing and validating weather models, such as the so-called Weather Generators. In such models the distribution for a given variable (such as daily precipitation or daily temperature, for instance) is often assumed to have a given shape, whose parameters for a given locality are then estimated from the aggregate historical data. For example, some models in the literature assume that data in the time series of daily temperatures in a given location are normally distributed within a month, whose mean and variance are estimated by pooling together several years of daily observations corresponding to that month. However, if there were changes in these parameters along the years then such a procedure would result in an inaccurate estimate for that distribution; what is worse, even if the variable were normally distributed in each month of the entire period such a change in the parameters along the years would cause a deviation from the normal shape in the aggregate distribution, and thus the shapes of the observed and the predicted distribution would not match. So, it is of crucial importance to test whether a given shape --and eventually the same parameters-- is in fact shared by all distributions of a weather variable within a given month and in several levels of aggregation (as, for instance, in the aggregate of all years of observation, or all months of a given year, or all weather stations). This is not a trivial task due to the fact that the time series for a single month is very short, and the usual statistical tools to test the distribution shapes for data in a single time series would not be suitable in this case. In this work we will address this issue by using a test procedure former introduced to investigate the shape of distributions of growth rates data (with a few time records) of business firms in large aggregates, and will apply such a procedure to test some alternative distribution shapes for balanced panels of weather data sets available in the software Wolfram Mathematica®.

**Acknowledgements:**JTL thanks CNPq/Brazil for partial support.

66

**Eugen Anitas**, Joint Institute for Nuclear Research, Bogoliubov Laboratory of Theoretical Physics, Russian Federation

**Abstract:**The fractal region of small-angle scattering intensity from monodisperse deterministic mass fractals with a single scale shows a complex superposition of maxima and minima on a simple power-law decay in a double logarithmic scale. The number of minima equals the number of fractal iterations, and the distance between their projections on 'q'-axis (where q is the magnitude of the scattering vector) is related to the fractal scaling factor. We show here that for certain classes of fractals the single minima are replaced by groups of minima, and the number of minima in each group is related to the number of centrosymmetric fractal substructures situated at various distances from the fractal's centre. These findings provide new insights about the structure of mass fractals and we illustrate them on a newly developed theoretical model based on 2D Cantor-like sets.

67

**Muddassir Ali**, Laurentian University, Bharti School of Engineering, Canada

**Redhouane Henda**, Laurentian University, Bharti School of Engineering, Canada

**Abstract:**Pulsed electron beam ablation (PEBA) has proven to be a promising and powerful technique for the growth of high quality thin films. Pulsed electron beam film deposition consists of many physical processes including target material heating, target ablation, plasma plume expansion, and film growth on a substrate. Plasma plume expansion into a vacuum or an ambient gas is a fundamental issue in PEBA as the quality of thin films deposited onto the substrate depends on the composition, energy and density of particles ejected from the target. In the present study, gas-dynamics equations are solved to investigate plasma expansion induced by interaction of a nanosecond electron beam pulse (~100 ns) with a graphite target in an argon atmosphere at reduced pressure. The effect of the electron beam efficiency and power density on plume expansion is assessed. The temperature, pressure, velocity and density profiles of the plasma plume are numerically simulated as a function of spatial dimension. The model is validated by comparing some of the obtained simulation results with experimental data available in the literature.

68

**Gerard Berginc**, THALES, R&T, France

**Abstract:**In this paper, we consider the coherent component of the electromagnetic wave field inside random media composed of dielectric or metallic nanoscale particles. The subject of our interest concerns a random medium with rough boundaries, consisting of a statistical ensemble of different scattering species and artificial material structures developed on base of dielectric or metallic resonant particles. The starting point of our theory is the multiple scattering theory, the averaged electric field satisfies a Dyson equation with a mass operator related to the effective dielectric permittivity of the homogenized structure. Quantum multiple scattering theory has been transposed in this electromagnetic case. We give a formal solution for the mass operator by introducing the T-matrix formalism. We show that the T-matrix satisfies a Lippman-Schwinger equation. Then, we introduce the Quasi-Crystalline Coherent Potential Approximation (QC-CPA), which takes into account the correlation between the particles with a pair-distribution function.The mass operator includes geometric effects, caused by resonant behavior due to the shape and size of particles, cluster effects because of correlations between particles. Significant modifications of particle scattering properties can be observed. The first intent of this paper is to establish new formulas for the effective dielectric constant which characterizes the coherent part of an electromagnetic wave propagating inside a random medium with randomly rough surfaces. In the high frequency limit the expression includes the vectorial case generalization of the result obtained by Keller, which has been derived, using a scalar theory. The equations are identical to equations obtained previously by Keller under conditions that the dyadic Green’s function is replaced by the scalar Green’s function. We obtained an important tensorial generalization of the notion of effective permittivity. Another aim of this paper is to discuss the contribution of light interaction with high-index dielectric nanoparticles supporting optically induced electric and magnetic Mie resonance. These nanoparticles are randomly distributed inside the film with rough boundaries. Near plasmon resonance frequencies, non-Rayleigh light scattering with unusual scattering diagrams can be observed. Attention is paid also to cooperative phenomena in this ensemble of randomly dispersed nanoparticles. We can demonstrate that dielectric structures allow to control both magnetic and electric components of light in a desirable way.

69

**Eva Oliveira**, Unicamp, , Brazil

**Mayara Maia**, Unicamp, , Brazil

**José Nunhez**, Unicamp, , Brazil

**Lubienska Cristina Lucas Jaquiê Ribeiro**, Unicamp, Saneamento, Brazil

**Abstract:**The Computational Fluid Dynamics (CFD) is used to numerically simulate the behavior of fluid flow, and all the laws that rule the study of fluids, both the energy and mass transfer, chemical reactions, hydraulic behavior, and other applications. Among these innumerable applications, the local behavior of a mixing tank can be estimated using the CFD techniques. In order to follow-up studies in the literature giving them increased credibility, this project aimed to determine whether both the scheduling Potency/Volume and Torque/Volume maintain the necessary level of agitation in the mix tank to extrapolate the volume 10 times and 100 times, using an axial impeller blades inclined at 45 degrees type (PBT 45), and a radial impeller (Rushton turbine) with the use of the CFD. With the results, was analyzed if a criterion is superior to the other in terms of maintaining the level of agitation in scheduling. The case study made it possible to achieve the main objective, deepen studies in Computational Fluid Dynamics. Therefore, managed to the largest domain in tools offered by CFX-ANSYS 14.0 software, allowing future work to be conducted.

70

**Xuesong Yuan**, University of Electronic Science and Technology of China, School of Physical Electronics, China

**Abstract:**In order to develop microwave and millimetre wave vacuum electron radiation source devices based on carbon nanotube field emission cold cathode, simulation model of carbon nanotube field emission cold cathode is investigated in this paper. Field emission current curve of a carbon nanotube cold cathode is firstly obtained according to experimental test. Using a numerical fitting method based on the curve, approximate Fowler-Nordheim constants A and B are calculated. Then, A and B are used in a 3D particle in cell simulation software and simulation model of carbon nanotube will be established. Carbon nanotube vacuum electron devices can be designed by using this model.

71

**Leonid Afremov**, Far Eastern Federal University (FEFU), Department of Theoretical and Experimental Physics, Russian Federation

**Aleksandr Petrov**, Far Eastern Federal University, School of Natural Sciences, Russian Federation

**Abstract:**The aim of this work is to modelling the dependence of the magnetic phase transition’s critical temperature on the NPs’ size and concentration of the “magnetic atoms”. The study was based on “average spin” method developed earlier [1] and it main principles are given below: 1) N “magnetic atoms” are distributed uniformly over M sites of cubic-shaped NP with L size with concentration p=N⁄M ; 2) fields of interaction h between atom’s spin magnetic moments are distributed randomly and there is direct exchange interaction between nearest neighbors only. Atom’s spin magnetic moments are oriented along axis oz (in the approximation of the Ising model). The calculation results provide a completely natural results: a) reducing of the NP’s size leads to a drop phase transition’s temperature as well as reducing of the ultrathin films’ thickness [2]; b) the critical temperature decreases with concentration of the “magnetic atoms” due to decreasing of the number of the interacting atoms. Moreover, at a certain critical concentration Curie temperature drops to zero - the system moves from an ordered to a disordered (paramagnetic) state. [1] Y. Kirienko, L. Afremov. Advanced Materials Research 472-475, 1827-1830 (2012). [2] T. Ambrose, C.L. Chien. Physical review letters 76(10), 1743-1746 (1996).

72

**Tatyana Gnitetskaya**, Far Eastern Federal University, School of Natural Sciences , Russian Federation

**Abstract:**In this paper the information model of intradisciplinary connections and semantic structures method are described. The information parameters, which we use in information model, are introduced. The question we would like to answer in this paper is – how to optimize the Physics Course’ content. As an example, the differences between entropy values in the contents of physics lecture with one topic but different logics of explanation are showed.

73

**Elias Vagenas**, Kuwait University, Department of Physics, Kuwait

**Abstract:**It had been shown earlier that the Generalized Uncertainty Principle (GUP), motivated from theories of quantum gravity, suggests quantization of measured lengths, areas and volumes at the fundamental level in flat spacetimes. In this work, we generalize this to curved spacetimes, by adding a weak gravitational field to the GUP modified Schrodinger, Klein-Gordon, and Dirac equations, and show that the discreteness results continue to hold. However, the nature of this new quantization is quite complex, although under proper limits, it reduces to cases without gravity. This supports the universal nature of quantum gravity effects.

74

**Halina Grushevskaya**, Belarusian State University, Physics Department, Belarus

**George Krylov**, Belarusian State University, Physics Department, Belarus

**Abstract:**Based on earlier developed quasi-relativistic self-consistent Dirac - Hartree - Fock field approximation for graphene [1,2], we perform simulations of charge transport properties of monolayer graphene both in optical case and near zero-frequency limit with accounting of spatial dispersion. Contributions of different nature to conductivity have been considered. Problem of monolayer graphene minimal conductivity is discussed in detail. [1] H.V. Grushevskaya and G.G. Krylov. Graphene: Beyond the Massless Dirac's Fermion Approach // In: Nanotechnology in the Security Systems, NATO Science for Peace and Security Series C: Environmental Security, J. Bonča and S. Kruchinin (eds.). (Springer, 2015). P. 21-31 [2] H.V. Grushevskaya, G. Krylov. Quantum Field Theory of Graphene with Dynamical Partial Symmetry Breaking // J. of Modern Phys. Vol. 5, Р. 984-994 (2014) [3] H.V. Grushevskaya, G. Krylov. Partially Breaking Pseudo-Dirac Band Symmetry in Graphene // J. Nonlin. Phen. in Complex Sys. Vol. 17, Р.86-96 ( 2014)

75

**Adrian Sotomayor**, Antofagasta University, Department of Mathematics, Chile

**Alvaro Restuccia**, Antofagasta University, Physics Department, Chile

**Abstract:**We present a local Bäcklund Wahlquist-Estabrook type supersymmetric transformation for the Korteweg-de Vries equation. As in the scalar case, such type of transformation generates infinite hierarchies of solutions by himself and also implicitly gives the associated (local) conserved quantities. A nice property is that every of such hierarchies admits a nonlinear superposition principle, starting for a initial solution, including as a particular case the multisolitonic solutions of the system. We comment about on a particular non-local of such transformations and also about on the related algebra of the non-local conserved quantities.

76

**Gillenhaal Beck**, Washington and Lee University, Physics and Engineering, United States

**Mohammad Abudayyeh**, Washington and Lee University, Physics and Engineering, United States

**Rebecca Melkerson**, Washington and Lee University, Physics and Engineering, United States

**Dan Mazilu**, Washington and Lee University, Physics and Engineering, United States

**Irina Mazilu**, Washington and Lee University, Physics and Engineering, United States

**Abstract:**We present experimental data and theoretical analysis which sheds new light on the rate and behavior of polymer/nanoparticle self-assembly. We study the adsorption of anionic silica nanoparticles in a colloidal suspension to cationic poly(diallyldimethylammonium chloride) (PDDA) polymer on glass and report on the time scale at which the surface coverage of the substrate reaches a steady-state when no further adsorption occurs. Our novel experimental method allows us to study particle self-assembly at time scales down to a few hundredths of a second, three orders of magnitude smaller than previously reported in literature. We compare our data with existing stochastic models in an analysis which both confirms the models’ prediction of pre-steady-state adsorption rates and allows us to incorporate real time into what was previously a theoretical calculation with arbitrary time units.

77

**Rebecca Melkerson**, Washington and Lee University, Physics and Engineering, United States

**Gillenhaal Beck**, Washington and Lee University, Physics and Engineering, United States

**Carlos M. da Fonseca**, Kuwait University, Mathematics, Kuwait

**Irina Mazilu**, Washington and Lee University, Physics and Engineering, United States

**Dan Mazilu**, Washington and Lee University, Physics and Engineering, United States

**Abstract:**We discuss the general methodology of finding exact solutions for two state systems using matrix theory. We analyze specific tridiagonal matrices with known spectra and propose new ones with eigenvalues that can be calculated exactly. We analyze the associated intractable eigenvectors numerically and present exact solutions for a specific experimental application, the ionic self assembly of nanoparticles. We also explore other applications of these tridiagonal matrices for a variety of other physical systems, such as epidemics, voter models, and traffic patterns.

78

**A. Zaidi**, Universiti Teknologi MARA, Mathematics, Malaysia

**Khairul Rasmani**, Universiti Teknologi MARA, Mathematics, Malaysia

**Abstract:**Demand for clean and treated water is increasing all over the world. Therefore it is crucial to conserve water for better use and to avoid unnecessary, excessive consumption or wastage of this natural resource. Classification of excessive domestic water consumption is a difficult task due to the complexity in determining the amount water usage per activity, especially as the data is known to vary between individuals. In this study, classification of excessive domestic water consumption is carried out using a well-known Fuzzy C-Means (FCM) clustering algorithm. Consumer data containing information on daily, weekly and monthly domestic water usage was employed for the purpose of classification. The result produced by the FCM clustering algorithm is compared with the result obtained from a statistical control chart which demonstrates the potential use of FCM clustering algorithm for the classification of domestic consumer water consumption data.

**Acknowledgements:**This research work is funded by the Ministry of Higher Education, Malaysia under the Fundamental Research Grant Scheme (FRGS) with reference number 600-RMI/FRGS TD 5/3 (1/2015). The authors also would like to thank the Institute of Research Management and Innovation (IRMI), Universiti Teknologi MARA, Malaysia.

79

**Daniele Pires Magalhães**, UNIVERSIDADE FEDERAL DE JUIZ DE FORA, PÓS GRADUAÇÃO EM MODELAGEM COMPUTACIONAL, Brazil

**Abstract:**The Calcium induced Calcium release (CICR) is mediated by the ryanodine receptor (RyR) on the sarcoplasmic reticulum (SR) and contributes to the global intracellular Calcium transient that drives cardiomyocyte contraction . In various disease settings, spontaneous Calcium release (SCR) from the SR gives rise to arrhythmic ectopic activity. SCR occurs in the form of waves of self propagating CICR. These waves result in Calcium dependent depolarizing currents generated by the electrogenic Na/Ca exchanger activity on the cell membrane giving rise to delayed afterdepolarizations (DADs) and consequently to triggered arrhythmias and Spontaneous Action Potential(SAPs). Our study is to improve the model Bondarenko to play DADs, SCR and SAPs. The main goal now is to replace or improve chain Markov Ryanodina to deliver results in accordance with existing experimental data.

**Acknowledgements:**UFJF, Capes e Fapemig.

80

**Inna Belashova**, Kazan Federal University , radiophysics, Russian Federation

**Vladimir Bochkarev**, Kazan federal university, Radiophysics, Russian Federation

**Abstract:**In this paper, nonlinear filtering algorithms for time series are tested using statistical modeling. The goal is to find a representation of the time sequence in a form of a wavelet series with a small number of non-zero coefficients (in other words to substantially reduce the dimensionality of the data). Thus, only statistically significant data remain in the filtered series. There are well-known efficient algorithms for nonlinear wavelet filtering for the case when the values of the time series have a normal distribution. In the case when distribution is not normal, good results can be obtained using the maximum likelihood estimation. We describe filtering using the criterion of maximum likelihood on the example of Poisson time series. This case is of great interest because such series are common for various problems of science and technology (for example, ranks of word usage frequency, number of sunspots, network traffic). We use different stochastic (genetic algorithms, simulated annealing) and deterministic optimization algorithms for direct optimization of the likelihood function. The algorithm was tested on model series, as well as on empirical data (ranks of rare words frequencies and number of sunspots were used). The conducted statistical modeling has shown that filtering using the maximum likelihood method can be significantly more efficient compared to the known algorithms for the case of Poisson series. The most promising methods of optimization were also selected for a given task.

81

**Pierre MILLET**, Université Paris-Sud, Chemistry, France

**Abstract:**The purpose of this communication is to present recent developments in gas-phase impedance spectroscopies and to outline their ability to provide accurate and detailed information on solid-gas reaction mechanisms. The discussion will focus on solid-hydrogen interactions, with a special emphasis on hydriding and permeation reactions. Basic principles of pneumato-chemical impedance spectroscopy (PIS) were introduced in 2005 [1]. PIS can be viewed as the direct transposition of electrochemical impedance spectroscopy (EIS) to solid-gas reactions. In PIS analysis, an analogy is made between pressure and electrical potential on one hand, gas flow and electric current on the other hand, and pneumato-chemical transfer functions are derived from gas phase measurements. Potentially, this spectroscopy can be used to analyze the dynamics of any solid-gas system including adsorption (surface) and absorption (bulk) phenomena, gas (H2) permeation across metallic membranes and electrocatalysis in gaseous fuel cells. Experiments can be performed using simple and conventional Sievert’s-type gas distribution apparatus. Measurements are well adapted to systems showing significant hysteresis. When applied to hydrogen absorption by intermetallic compounds [2] (a process of great practical interest for hydrogen storage applications) or to hydrogen permeation across metallic membranes [3], surface and bulk rate constants can be measured separately, opening the way to detailed kinetics analysis at microscopic levels and to system optimization. State-of-art as well as recent advances in the field will be presented and some perspectives will be discussed. references [1] P. Millet, Pneumato-chemical impedance spectroscopy : (I) principles, J. Phys. Chem., 109 (2005) 24016 – 24024 [2] P. Millet, Pneumato-chemical impedance spectroscopy : (II) dynamics of hydrogen sorption by metals, J. Phys. Chem., 109 (2005) 24025-24030 [3] C. Decaux, R. Ngameni, A. Ranjbari, S. Grigoriev, P. Millet, Dynamics of hydrogen permeation across metallic membranes, Int. J. Hydrogen Energy, 38 (2013) 8584–8589

82

**Evangelos Bakalis**, University of Bologna, Dipartimento di Chimica "G. Ciamician", Italy

**Abstract:**The physicochemical properties of a metallic carbon nanotube (SWCNT) are altered when either a molecule or a neutral distribution of point charges approaches in the vicinity of its external surface. The dipole moment term is responsible for the induced energy at the surface of the tube, which is not just the summation of individual contributions because of interference and screening effects. By using the Green’s expansion method the excess of dipole energy in the vicinity of the tube is analytically investigated. We calculate the total potential energy of the system. The most stable orientations of the dipole along the SWCNT are obtained, and for each one of them an analytical simple function is given. This function can be used either directly by experimentalists or as starting point in more sophisticated molecular modelling. The orientation of the dipole along to the z-axis of the tube is the most stable.

83

**Dérek Bomfim Prates**, UFVJM, ICET, Brazil

**Caio Jardim**, Federal University of the Valleys of Jequitinhonha and Mucuri, Institute of Science and Technology, Brazil

**Letícia Figueiredo**, UFVJM, , Brazil

**Jaqueline Silva**, UFVJM, ICET, Brazil

**Maurício Kritz**, LNCC, , Brazil

**Abstract:**Several mathematical models can be found in the literature for describing and analyzing epidemics. Some of these models allow the vaccination study of a determined population in an epidemic scenario. This work will address and compare two important vaccination strategies in a modified SIR model: a constant vaccination, where a certain parcel of the population is immediately and always vaccinated after its birth, and a vaccination in pulses, where in each time interval the vaccine is applied on a parcel of the population. Furthermore, an analysis of the epidemic behavior of different populations submitted to the same vaccination strategy, but whose vaccinated parcels are different, is performed. The analysis and results are performed through numerical solutions of the model and a conditional function that models the pulse vaccination periodicity is proposed.

**Acknowledgements:**We thank CNPQ for the scientiﬁc scholarship given to the undergraduate student who participates in this work.

84

**Seiki Saito**, National Institute of Technology, Kushiro College, Electrical Engineering, Japan

**Abstract:**Tungsten material attracts many researchers because it is planned to employ for divertor plates of nuclear fusion device. The material is also expected for the engineering use, because it is discovered that fibrous structure is generated when tungsten irradiated by helium plasma under specific condition. Moreover, it is known that bubble structure is formed as a precursor of the fibrous structure under helium plasma irradiation. The existence of the bubbles may affect the plasma-material interaction. In this study, therefore, binary-collision-approximation-based (BCA) simulation is performed for the investigation of the effects of the existence of helium bubbles on the retention process of incident hydrogen atoms.

85

**Wesley Oliveira**, Unicamp, , Brazil

**Marta Pires**, Unicamp, , Brazil

**Laura Canno**, Unicamp, , Brazil

**Lubienska Cristina Lucas Jaquiê Ribeiro**, Unicamp, Saneamento, Brazil

**Abstract:**The Computational Fluid Dynamics (CFD) is a tool used to numerically simulate fluid flow behavior, and all the laws that govern the study of fluids is the mass transfer and energy, chemical reactions, hydraulic behaviors, among others applications. This tool mathematical equations solves the problem in a specific manner over a region of interest, with predetermined boundary conditions on this region. This work is to study the flow channel through the CFD technique.

86

**Karen Akiyama**, Unicamp, , Brazil

**Caroline Bressan**, Unicamp, , Brazil

**Marta Pires**, Unicamp, , Brazil

**Laura Canno**, Unicamp, , Brazil

**Lubienska Cristina Lucas Jaquiê Ribeiro**, Unicamp, Saneamento, Brazil

**Abstract:**The issue of water pollution has worsened in recent times due to releases, intentional or not, of pollutants in natural water bodies. This causes several studies about the distribution of pollutants are carried out. The water quality models have been developed and widely used today as a preventative tool, ie to try to predict what will be the concentration distribution of constituent along a body of water in spatial and temporal scale. To understand and use such models, it is necessary to know some concepts of hydraulic high on their application, including the longitudinal dispersion coefficient. This study aims to conduct a theoretical and experimental study of the channel dispersion coefficient, yielding more information about their direct determination in the literature.

87

**Rodrigo Signoreti**, Unicamp, , Brazil

**Rafael Camargo**, Unicamp, , Brazil

**Laura Canno**, Unicamp, , Brazil

**Marta Pires**, Unicamp, , Brazil

**Lubienska Cristina Lucas Jaquiê Ribeiro**, Unicamp, Saneamento, Brazil

**Abstract:**Challenged with the high rate of leakage from water supply systems, these managers are committed to identify control mechanisms. In order to standardize and control the pressure Pressure Reducing Valves (VRP) are installed in the supply network, shown to be more effective and provide a faster return for the actual loss control measures. It is known that the control pressure is while controlling the occurrence of leakage. Usually the network is sectored in areas defined by pressure levels according to its topography, once inserted the VRP in the same system will limit the downstream pressure. This work aims to show the importance of VRP as loss reduction for tool.

88

**roger nakad**, Notre Dame University Louaizé, Department of Mathematics and Statistics, Lebanon

**Abstract:**We extend the Friedrich inequality for the eigenvalues of the Dirac operator on Spinc manifolds with boundary under four different boundary conditions. Limiting cases are then studied and many examples and applications in physics are given.

89

**Sumbul Sehar**, Government college university Lahore, Physics Department, Pakistan

**Abstract:**In many physical situations such as space or laboratory plasmas a hot low-density electron populations can be generated superimposed on the bulk cold population, resulting in a particle distribution function consisting of a dense cold part and hot superthermal tail. Space observations show that electron distributions are often observed with flat top at low energies and high energy tails. The appropriate distribution to model such non-Maxwellian features is the generalized (r,q) distribution function which in limiting forms can be reduced to kappa and Maxwellian distribution functions. In this study, Kinetic model is employed to study the electron-acoustic and ion-ion acoustic instabilities in four component plasma with generalized (r,q) distribution function for both magnetized and unmagnetized plasmas. Departure of plasma from Maxwellian distributions significantly alters the growth rates as compared to the Maxwellian plasma. Significant growth observed for highly non-Maxwellian distributions as well as plasmas with higher dense and hot electron population. Existence of weak damping is also established when the distribution contains broadened flat tops at the low energies or tends to be Maxwellian. These results may be applied in both experimental and space physics regimes.

90

**Abdeen Omer**, Energy Research Institute (ERI), Department of Renewable Energy, United Kingdom

**Abstract:**Globally buildings are responsible for approximately 40% of the total world annual energy consumption. Most of this energy is for the provision of lighting, heating, cooling and air conditioning. An increase in awareness of the environmental impact of CO2, NOx and CFCs emissions triggered a renewed interest in environmentally friendly cooling and heating technologies. Under the 1997 Montreal Protocol, governments agreed to phase out chemicals used as refrigerants that have the potential to destroy stratospheric ozone. It was therefore considered desirable to reduce energy consumption in order to decrease the rate of depletion of world energy reserves as well as the pollution to the environment. One way of reducing building energy consumption is to design buildings, which are more efficient in their use of energy for heating, lighting, cooling and ventilation. Passive measures, particularly natural or hybrid ventilation rather than air-conditioning, can dramatically reduce primary energy consumption. Therefore, promoting innovative renewable energy applications including the ground source energy may contribute to preservation of the ecosystem by reducing emissions at local and global levels. This will also contribute to the amelioration of environmental conditions by replacing conventional fuels with renewable energies that produce no air pollution or the greenhouse gases (GHGs). An approach is needed to integrate renewable energies in a way to achieve high building performance standards. However, because renewable energy sources are stochastic and geographically diffuse, their ability to match demand is determined by the adoption of one of the following two approaches: the utilisation of a capture area greater than that occupied by the community to be supplied, or the reduction of the community’s energy demands to a level commensurate with the locally available renewable resources. Ground source heat pump (GSHP) systems (also referred to as geothermal heat pump systems, earth-energy systems and GeoExchange systems) have received considerable attention in recent decades as an alternative energy source for residential and commercial space heating and cooling applications. The GSHP applications are one of three categories of geothermal energy resources as defined by ASHRAE and include high-temperature (>150°C) for electric power production, intermediate temperature (<150°C) for direct-use applications and GSHP applications (generally (<32°C). The GSHP applications are distinguished from the others by the fact that they operate at relatively low temperatures.

91

**Hafida CHEMOURI**, preparing school in sciences and technics tlemcen, chemistry, Algeria

**Abstract:**The hetero Diels Alder reaction (HAD) of 1-aza-1,3-butadiene with dimethylvinylamine is an efficient method for the synthesis of nitrogen heterocycles that constitute good intermediates for the preparation of complicated heterocycles and pharmacological systems. Experimentally, it has been found that the cycloaddition of 1-aza-1,3-butadiene derivatives (R = H, (CH3)CO, (CH3)CO–BH3) with dimethylvinylamine gives preferentially the ortho regioisomers17. Our aim in this work is to put in evidence the regioselectivity of these HDA reactions using four theoretical models based on the use of reactivity indices. The quantum chemical calculations were performed using the B3LYP/6-31G(d) implemented in Gaussian 09 suite of programs.

92

**Meriem MERAD - BENAROU**, university of Tlemcen , chemistry, Algeria

**Abstract:**The knowledge of the enzymes is capital, because it is catalyze the majority of the chemical reactions of the living organisms. The deterioration of the enzymes by specific inhibiters makes it possible to block the biochemical ways, only a particular site of enzymatic protein comes into contact with its substrate: it is the active site, At the active site, the behavior of atoms is quite different in the approach of the enzyme and substrate. In this approach, there are enzyme substrate interactions that govern the stability of an ES complex to produce the commodity according to the specificity of the enzyme. It is here that lies the difficulty of understanding the interactions between enzyme and substrate, starting from this idea we propose to elucidate the mechanism by molecular modeling. Our research objective is to study the inhibition of the phosphathase Cdc25 through the molecular modeling methods.

93

Use of two families species of Brassica and Solanaceae source in the fight against Alzheimer disease

**ABDELLI Imane**, a: Preparatory School in Science and Technology- Tlemecen-Algeria. b: University Abou-BakrBelkaid-Faculty of Science- Department of Chemistry- Laboratory of Natural Substances and Bioactive(LASNABIO)-Tlemcen-Algeria., Chemistry, Algeria

**Abstract:**Since always, the medicinal herbs were used to prevent or to look after various diseases and the current drugs have, for the majority of them, a natural origin. Research on natural substances has been a carrying topic for a few years and the pharmaceutical laboratories, always in search new compounds active, turn more and more to the identification and the characterization of molecules resulting from natural matrices, and, take as a starting point their molecular structure to imagine new drugs. The flavonoids are natural compounds, which can be met in a broad fruit variety and vegetables consumed daily by the human being. Besides their role in the plants pigmentation, some of these compounds present interests biological activities, such as anti-radicalizing and antioxidant actions. Their large structure varieties (9000 identified compounds) biological activities (antioxydant, anti-allergic, antiviral, anti-inflammatory and drug inhibiting) flavonoids make a particularly complex study object. Studies suggested that the oxidative stress has a critical role in the Alzheimer's disease (MA) is responsible with the β-amyloid neurotoxicity. Attenuation oxidative stress by the anti oxidant molecules is proposed as a potential treatment therapeutic in MA. This disease is characterized by the presence of β-amyloid deposits in the cerebral blood-vessels. The flavonoids can act in various ways in the processes of the stress oxidizing regulation, are known primarily for the protection of the cellular walls and of the blood-vessels they are regarded as good β-secretase inhibitors. Our item consists to study the inhibition of the β-secretase, enzyme implied in the Alzheimer's disease with a Brassica and Solanaceae family by molecular modeling methods. The purpose of this study is to decrease the formation of β-Amyloïd and consequently to delay its progression.

94

**Constantinos Vayenas**, University of Patras, LCEP , School of Engineering, Greece

**Abstract:**Abstract. We compute analytically the masses, binding energies and free energies of gravitationally bound Bohr-type states via the rotating relativistic lepton model which utilizes the de Broglie wavelength equation in conjunction with special relativity and Newton’s relativistic gravitational law. The latter uses the inertial-gravitational masses, rather than the rest masses, of the rotating particles [1, 2]. The model also accounts for the electrostatic charge-induced dipole interactions between a central charged lepton, which is usually a positron, with the rotating relativistic lepton ring. We use three rotating relativistic neutrinos to model baryons, two rotating relativistic neutrinos to model mesons and a rotating relativistic electron neutrino - positron (or electron) pair to model the W± bosons [3]. It is found that gravitationally bound states comprising three relativistic neutrinos have masses in the baryon mass range (∼ 0.9 to 2 GeV/c2), whereas states comprising two neutrinos have masses in the meson mass range (∼ 0.2 to 0.8 GeV/c2) [1, 2]. It is also found that the rest mass values of quarks are in good agreement with the heaviest neutrino mass value of 0.05 eV/c2 and that the mass of W± bosons (∼ 81 GeV/c2) corresponds

95